gabort.m

一个很不错的matlab 程序

clf;
clear all ;
close all; 
fs=100;   %采样率。
Ts=1/fs;   
t=-2:Ts:18;  
gass=2^(1/4)*exp(-pi*t.^2); %高斯函数
subplot(3,2,1),plot(gass);
title('高斯函数');
xlabel('t');
ylabel('幅度');
T=0:Ts:20;
ft=cos(T.^2+2*T)+cos(T.^2);  %信号函数
subplot(3,2,2),plot(ft);  
title('信号函数');
xlabel('time');
ylabel('幅度');
y=fft(ft);       %信号的FFT变换
amp=abs(y);
subplot(3,2,3);plot(amp);
title('信号的FFT变换');
xlabel('F(Hz)');
ylabel('幅度');
shl=100;           %窗平移每次点数
shn=(length(t)-1)/shl;   %求窗平移次数
for k=0:shn-1;
    gassc=2^(1/4)*exp(-pi*(t-k*shl*Ts).^2); %平移后的高斯函数
    yl=conv(ft,gassc);                %短时傅立叶变换
    y2=zeros(shn,4001);
    y2(k+1,:)=yl;
end
[F,T]=size(y2);
[F,T]=meshgrid(1:T,1:F);
subplot(3,2,4),mesh(F,T,abs(y2))
title('信号 T-F图');
xlabel('F(Hz)');
ylabel('time')
zlabel('幅度');
subplot(3,2,5),contour(F,T,abs(y2))