clcclearallfs=input('采样频率为:fs=');N=1024;t=(0:N-1)/fs;y=0.5*sin(20*pi*t)+0.6*sin(180*pi*t)+0.8*cos(300*pi*t);subplot(2,1,1)plot(t,y);gridontitle(['抽样信号时域波形图,抽样频率fs=',num2str(fs)]);xlabel('t/ms');ylabel('y')f=(0:N-1)*fs/N;y1=fft(y,N);%对信号进行快速Fourier变换y2=abs(y1)/N*2subplot(2,1,2);plot(f(1:N/2),y2(1:N/2));gridontitle('信号的频谱图');xlabel('f/Hz');ylabel('Rex(f)');