求(短时傅立叶变换的程序)
多谢!!!!!!!!!!!!!!! 本帖最后由 wdhd 于 2016-9-14 10:22 编辑我也需要,那位大哥有,小弟多谢了。
我的邮箱是:zhouql_0991@163.com
论坛里有 本帖最后由 wdhd 于 2016-9-14 10:22 编辑
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;
=meshgrid(fnew,tnew);
mesh(F,T,abs(TF))
subplot(224)
contour(F,T,abs(TF))
本帖最后由 wdhd 于 2016-9-14 10:23 编辑
多谢zhao19764
我的毕业设计还差这个程序了
回复:(xiaozuov)求(短时傅立叶变换的程序)
时频分析工具箱中就有回复:(simon21)回复:(xiaozuov)求(短时傅立叶变...
本帖最后由 wdhd 于 2016-9-14 10:23 编辑多谢哈
多谢,找了好长时间了。 我就是需要这个来研究
页:
[1]