求教Matlab的 gabor变换

kristycreasy

N1=512;
fs=2000;
t=(0:N1-1)/fs;
Fm=zeros(3,length(t));
F=0;
f=[300,500,700];
A=[1,1,1];
d=[15,15,15];

for m=1:3;
    Fm(m,: )=A(m)*exp(-3.14*d(m)*t+i*(2*pi*f(m)*t));
end

F=sum(Fm);
F=F+awgn(F,inf);
FR=real(F);

plot(FR);


Q=1;
N=32;
Nh=39;
%tt=(0:N-1)/fs;
%h1=exp(-3.14*15*tt);
g=window(@gausswin,Nh); %如何用以指数衰减的函数做窗函数啊？[attachimg]50486
g=g/norm(g);
[tfr,dgr,gam] =tfrgabor(FR,N,Q,g);