有谁做过stft的matlab程序啊，借我看看啊

yushouzhuo · 发表于 2006-9-18 16:35

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

有谁做过stft的matlab程序啊，能发给我看看么？发邮箱也行。cyc_19810215@yahoo.com.cn.谢谢了啊:handshake:handshake

jimin · 发表于 2006-9-18 16:46

发到matlab　讨论区吧
这里估计来回答问题的人不多的

FtpAdmin · 发表于 2006-9-19 10:10

看帖子http://forum.vibunion.com/forum/thread-17216-1-2.html

eight · 发表于 2007-2-18 22:42

贴一下代码如下：

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))

花瓣雨 · 发表于 2008-4-11 19:07

请问对新的时间坐标为什么要这样计算”Tn=(N-Nw)/L+1“？谢谢！

guoguo1983 · 发表于 2009-6-4 09:35

程序运行出错啊说k没定义:@(

friendchj · 发表于 2009-6-4 11:11

把第一行x=exp(j*pi*k*t.^2)
改成
% x=exp(j*pi*k*t.^2)
那是分析的对象

guoguo1983 · 发表于 2009-6-10 10:28

哦运行出结果了谢谢楼上了:@)

hotman007 · 发表于 2011-4-27 21:23

支持eight大大的无私奉献，谢谢了，收藏

song_bjtu · 发表于 2012-2-5 15:41

多谢分享，收藏

账号		自动登录	找回密码
密码			我要加入

[综合讨论] 有谁做过stft的matlab程序啊，借我看看啊

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复 4楼的帖子

回复 6楼 guoguo1983 的帖子

评分

[综合讨论] 有谁做过stft的matlab程序啊，借我看看啊

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复 4楼 的帖子

回复 6楼 guoguo1983 的帖子

评分

回复 4楼的帖子