仿真信号短时傅里叶变换

2、仿真信号短时傅里叶变换通过函数算法产生一个非平稳信号sig=amgauss(128).*fmlin(128);并作出信号的“时间t—幅值A”时域波形图，波形图如figure（1）所示；为原信号的基础上设置h=hamming(127)的窗函数，窗函数的图像如figure（2）所示；然后对信号进行短时傅里叶变换，并将变换结果时频表示，作出“时间—频率”图像，图像如figure所示。（1）仿真信号处理程序：sig=amgauss(128).*fmlin(128);%%时域波形figure(1);plot(real(sig),'LineWidth',2);xlabel('时间t');ylabel('幅值A');%设置窗函数h=hamming(127);figure(2);plot(h)[tfr,t,f]=tfrstft(sig,1:128,128,h);%时频表示figure(3);contour(t,f(1:length(f)/2),abs(tfr(1:length(f)/2,:)));xlabel('时间t');ylabel('频率f')（2）信号图像：020406080100120140-1-0.8-0.6-0.4-0.200.20.40.60.81时间t幅值AFigure（1）时变信号时域波形图02040608010012014000.10.20.30.40.50.60.70.80.91Figure（2）窗函数的图像时间t频率f2040608010012000.050.10.150.20.250.30.350.40.45Figure（3）短时傅里叶变化后的图像