matlab生成正弦波，fft观测频谱

L=128; f1=100; fs=400; t=0:1/fs:(1/fs)*(L-1); si=sin(2*pi*f1*t); f=fft(si,L); f=20*log(abs(f))/log(10); ft=[0:(fs/L):fs/2]; f=f(1:length(ft)); subplot(3,2,1);plot(t,si); xlabel('时间（s）'); ylabel('幅度（v）

>>阅读原文<<