正(余)弦函数频谱图.docx

clear；cf;Fs=1000;f0=50;Ts=lFs1.=100O;A=20;n=O:l-l;t=n*Ts;%时间序列x=A*sin(2*pl*f*t);%1Ssubplot(2,l,l),p!ot(tp<);%绘出时域波形XIabe1('时间its);yabel('振幅)；title('时域波形,)；NFFT=2Mextpow2(1.);%求出的住FFT变换长度yy=fft(×,NFFT);%对信号进行伸电叶变换f=Fs2inspace(0zl,NFFT2+l);%即11(FsNFFT)其中(Fs/NFFT)是频率分辨率y=abs(yy)21.;%幅值处理subplot(2,l,2),plot(f,y(l:NFFT241);%绘出的均率变化的振幅XIabeI(濒率i<Hz;ylabel('振幅)titled嫉谱图)gridon;