fs=1000; %采样频率 N=1024; %采样点数 n=0:N-1; t=n/fs; f0=100; %信号频率 x=sin(2*pi*f0*t); y=abs(fft(x,N)); %傅里叶变换后画出幅度谱 plot(y);
