%fs=1000;     %采樣率為1000
%t=0:1/fs:1;  %時(shí)長1s
%f1=50;         %信號(hào)頻率50Hz
fs=100;%設(shè)定采樣頻率
N=512;
n=0:N-1;
t=n/fs;
f1=20;
x=sin(2*pi*f1*t); %生成隨機(jī)噪聲的正弦波函數(shù)
x=x+randn(1,length(x));
subplot(211)
plot(x);
title('原正弦信號(hào)(含隨機(jī)噪聲)')
xlabel('序列(n)')
grid on

y=dfft(x);
n=0:length(y)-1;
f=fs*n/length(y);
subplot(212)
plot(f,abs(y));
title('原正弦信號(hào)(含隨機(jī)噪聲)的FFT')
xlabel('頻率Hz')
grid on