求（短时傅立叶变换的程序）

xiaozuov

多谢！！！！！！！！！！！！！！！

特殊符号

　　我也需要，那位大哥有，小弟多谢了。

　　我的邮箱是：zhouql_0991@163.com

zhao197642

论坛里有

zhao197642

　　x=exp(j*pi*k*t.^2)

　　clear,clc,close all

　　figure(1)

　　k=4;T=5;

　　fc=k*T;

　　fs=3*fc;

　　Ts=1/fs;

　　N=T/Ts;

　　x=zeros(1,N);

　　t=0:N-1;

　　x=exp(j*k*pi*(t*Ts).^2);

　　% x=awgn(x,-3,'measured');

　　subplot(221)

　　plot(t*Ts,real(x))

　　X=fft(x);

　　X=fftshift(X);

　　subplot(222)

　　plot((t-N/2)*fs/N,abs(X))

　　Nw=20;

　　L=Nw/2;

　　Tn=(N-Nw)/L+1;

　　nfft=32;

　　TF=zeros(Tn,nfft);

　　for i=1:Tn

　　xw=x((i-1)*10+1:i*10+10);

　　temp=fft(xw,nfft);

　　temp=fftshift(temp);

　　TF(i,:)=temp;

　　end

　　subplot(223)

　　fnew=((1:nfft)-nfft/2)*fs/nfft;

　　tnew=(1:Tn)*L*Ts;

　　[F,T]=meshgrid(fnew,tnew);

　　mesh(F,T,abs(TF))

　　subplot(224)

　　contour(F,T,abs(TF))

xiaozuov

　　多谢zhao19764

　　我的毕业设计还差这个程序了

simon21

时频分析工具箱中就有

xiaozuov

　　多谢哈

huanying8426

多谢，找了好长时间了。

mlh_2002

我就是需要这个来研究