Fs = 200;
t = 0:1/Fs:1;
fa = 5;
fb =37;
x = sin(2* pi*t*fa+pi/3)+ 0.6*sin(2* pi*t*fb);
x =awgn(x,20);
nfft = 1024; % comprimento FFT
X = fft(x,nfft); %Se comprimento de X for igual à nfft, a FFT é preenchida com zeros
Xshift=fftshift(X); % FFT centralizada
mx = abs(Xshift); % Magnitude da FFT de X centralizada
f = (-nfft/2:nfft/2-1)*Fs/nfft; % Frequencia igualmente espaçada
X1 = X(1:nfft/2); % Como a FFT simétrica, pegou-se somente a metade
mx1 = abs(X1); % Magnitude da FFT de X
f1 = (0:nfft/2-1)*Fs/nfft; % Frequencia igualmente espaçada
% Gráficos
figure(1);
plot(x);
title('Sinal Estudado');
xlabel('Tempo (s)');
ylabel('Amplitude');
figure(2);
plot(f,mx);
title('FFT do sinal - centralizada');
xlabel('Frequência (Hz)');
ylabel('Amplitude');
figure(3);
plot(f1,mx1);
title('FFT do sinal');
xlabel('Frequência (Hz)');
ylabel('Amplitude');