tf1.m

来自「这是一个时频分析中关于gabor变换的程序」· M 代码 · 共 20 行

clear   
fs=8000; 
t1=10; 
t2=20; 
f1=1000;
f2=2000;
le=fs*t2;                                    
s=fix(le/128);       %帧间重叠1/2,所以算出所需循环次数
hn=hanning(256)';   %得到STFT分析窗,汉宁窗,帧长256点
T1=[1:fs*t1,zeros(1,fs*t1)];
T2=[zeros(1,fs*t1),fs*t1+1:fs*t2];
d=sin(2*pi*f1*T1/fs)+sin(2*pi*f2*T2/fs);
for n=1:1:s-3         
    d1(1:256)=d((n-1)*128+1:(n-1)*128+256).*hn; %时域加窗 
    Xd(n,(1:256))=fft(d1,256);  %FFT  
end
[n,m]=size(Xd);
x=1:n;y=1:m;
mesh(y/m,x/n*t2,abs(Xd)),axis([0,0.4,0,20,0,100])