请教matlab高手，帮我看看这个STFT程序问题

我爱计算机 · 发表于 2009-9-29 16:50

怎么会是频率呢？X是原始信号，只要取的点数和t一致，就没有问题了啊！

songzy41 · 发表于 2009-9-30 08:54

原帖由 cuixuepeng 于 2009-9-11 21:51 发表

 登录/注册后可看大图

你好，打挠你一下：
我用这个方法画出的图如下，t-f图中的横坐标不是时间吧，应该是[tfr,t,f]=tfrstft(x,1:length(x),1024)里面的频率分析点数吧？那怎么把它改成实际的时间呢？比如说我采样频率是100Hz，那么1024个 ...

把LZ的数据做了STFT分析，程序为
xx=load('11-4.txt');
fs=100;
N=length(xx);
t=(0:N-1)/fs;
subplot 211; plot(t,xx);
xlabel('时间'); ylabel('幅值')
title('信号波形');
x=xx(1:1024);
M=1024;
h=hanning(1023);
[tfr,t,f]=tfrstft(xx,1:M,M,h);
subplot 212;
imagesc(t/fs,f(1:512)*fs,abs(tfr(1:512,:)));
axis('xy');
xlabel('时间'); ylabel('频率')
title('信号的STFT频谱');
得图在下面，从图中看出基频约在8Hz左右，而且随时间的推移，基频也略有变化。

[ 本帖最后由 songzy41 于 2009-9-30 08:56 编辑 ]

bcyangbc · 发表于 2009-9-30 12:35

这个是傅里叶变换吧？而且还是显示的双边谱！横轴对应的是频率，竖轴是功率吧？

michaelbeck · 发表于 2009-12-16 19:56

我是刚开始用matlab做时频仿真，能否请教高人几个问题啊？你们说的这个时频工具箱我不会用，高人指导一下可以吗？

michaelbeck · 发表于 2009-12-16 20:00

求教一下高手，我刚下载了您提供的这个时频分析工具箱，但是不知道，怎么用，能不能赐教一下，我是刚开始用matlab做仿真的小菜鸟，我的Q是305432522，麻烦您了

我爱计算机 · 发表于 2010-1-9 21:05

在论坛中搜一下，关于MATLAB工具箱添加方法有很好的总结

hnyanhua · 发表于 2010-1-27 19:58

我的机器上面也没有tfrstft这个函数啊，麻烦前辈发给我好吗？谢谢哈
我的邮箱是：yanhua627693@126.com
非常感谢啦

hnyanhua · 发表于 2010-1-27 20:36

sorry
后来看了楼上发的帖子，在http://zhlong.ys168.com/下载到了(*^__^*) 嘻嘻……，遇到问题同大家一块探讨哈

hnyanhua · 发表于 2010-1-27 21:04

点击File，选set path，在弹出的对话框中点add with subfolders按钮，选中你的工具箱tftb-0.2，保存既可。任选工具箱中的一个命令，用help检验一下。

xiangyu537 · 发表于 2010-3-29 13:09

能否解释32楼 imagesc(t/fs,f(1:512)*fs,abs(tfr(1:512,:))); 这句话是画时频图的但是里面的三个参数的意义是什么呢？

xiangyu537 · 发表于 2010-3-29 15:53

clc
clear all;
SampFreq = 80;
% t= 0 : 1/SampFreq : 5;
t1= 1/SampFreq: 1/SampFreq : 2.5;
t2= 2.5 : 1/SampFreq : 5;
N=400;
t=(0:N-1)/SampFreq ;
c1 = 2 * pi * 10;
c2 = 2 * pi * 20;
Sig(1:200)=sin(c1*t1);
Sig(201:401)=sin(c2*t2);
%Sig = sin(c1*t1)+sin(c2*t2)
%figure(1)
%plot(t,Sig,'Line Width',2);
%N=length(S);
%subplot 211; plot(t,xx);
%xlabel('时间'); ylabel('幅值')
%title('信号波形');
%x=xx(1:1024);
nfft=400;
M=1024;
h=hanning(1023);
[tfr,t,f]=tfrstft(Sig,(1:400)/SampFreq,M,h);
%tfr=tfrstft(Sig)
figure(2);
imagesc(t/SampFreq ,f(1:512)*SampFreq ,abs(tfr(1:512,:)));
%axis('xy');
xlabel('时间'); ylabel('频率')
title('信号的STFT频谱');

??? Error using ==> tfrstft
X must have one column

Error in ==> text at 24
[tfr,t,f]=tfrstft(Sig,(1:400)/SampFreq,M,h);
运行时出错

xiangyu537 · 发表于 2010-3-29 15:54

请高手解答

华电机械 · 发表于 2013-5-20 15:56

zhlong 发表于 2007-5-13 22:35

登录/注册后可看大图

下载tftb时频分析工具箱
或是留下邮箱，我发给你

您有STFT和WVD的工具箱吗？本人现在急用，，在此感激您能发一份给我。邮箱：354773461@qq.com

yghit08 · 发表于 2013-5-20 17:05

华电机械发表于 2013-5-20 15:56

登录/注册后可看大图

您有STFT和WVD的工具箱吗？本人现在急用，，在此感激您能发一份给我。邮箱：

time_frequency.rar (509.94 KB, 下载次数: 1)
见附件。

华电机械 · 发表于 2013-5-20 17:14

yghit08 发表于 2013-5-20 17:05

登录/注册后可看大图

见附件。

十分感激您啊，，有空向您请教！

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

[综合] 请教matlab高手，帮我看看这个STFT程序问题

回复 30楼 cuixuepeng 的帖子

回复 27楼 cuixuepeng 的帖子

回复地板 zhlong 的帖子

回复 34楼 michaelbeck 的帖子

回复 37楼 hnyanhua 的帖子