% ===================================================================== % MATLAB code to do FFT of a signal with noise % By Mahesha MG, MIT Manipal % Date: 11/10/2012 % ===================================================================== Fs = 1000; % Sampling frequency T = 1/Fs; % Sample time L = 1000; % Length of signal t = (0:L-1)*T; % Time vector % Sum of a 50 Hz sinusoid and a 120 Hz sinusoid x = 0.7*sin(2*pi*50*t) + sin(2*pi*120*t); y = x + 2*randn(size(t)); % Sinusoids plus noise figure(1) plot(t,y) title('Signal with noise') figure(2) NFFT = 2^nextpow2(L); % Next power of 2 from length of y Y = fft(y,NFFT)/L; f = Fs/2*linspace(0,1,NFFT/2+1); % Plot single-sided amplitude spectrum. stem(f,2*abs(Y(1:NFFT/2+1))) title('Single-Sided Amplitude Spectrum of y(t)') xlabel('Frequency (Hz)') ylabel('|Y(f)|')

## Thursday, October 11, 2012

### MATLAB - FFT of a signal with noise

Subscribe to:
Post Comments (Atom)

## No comments:

## Post a Comment