tf2.m

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

clear;
t=20;
fs=8000;
f1=1000;
f2=2000;
le=fs*t;  
a=1/16;b=16;
N=a*fs;M=fs/b;
s=fix(le/fs);       %帧间重叠1/2,所以算出所需循环次数
hn=boxcar(M)';   %得到STFT分析窗,汉宁窗,帧长256点
T=1:fs*t;
d=sin(2*pi*f1*T/fs)+sin(2*pi*f2*T/fs);
for n=1:1:s       
    d1(1:M)=d((n-1)*N+1:(n-1)*N+M).*hn; %时域加窗 
    Xd(n,(1:M))=fft(d1,M);  %FFT  
end
[n,m]=size(Xd);
x=1:n;y=1:m;
mesh(y/m,x,abs(Xd));axis([0,0.5,0,20,0,100])