fft(matlab)1

N=512;    %长度
f1=100;   %频率
f2=105;   
Fs=400;         %采样频率
t=0:1/Fs:1/Fs*(N-1);        %产生时间序列
s=sin(2*pi*f1*t)+sin(2*pi*f2*t);   %叠加信号

f=fft(s,N);  %求频谱
f=20*log(abs(f))/log(10);       %换成dBW单位
ft=0:(Fs/N):Fs/2;           %转换以hZ单位为横坐标
f=f(1:length(ft));

subplot(211);plot(t,s);
xlabel('时间(s)');ylabel('幅度(v)');title('时域波形');
subplot(212);plot(ft,f);
xlabel('频率(hz)');ylabel('功率(dBW)');title('信号频谱');