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

回复 30楼 cuixuepeng 的帖子

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

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

原帖由 cuixuepeng 于 2009-9-11 21:51 发表 http://www.chinavib.com/forum/images/common/back.gif
你好,打挠你一下:
我用这个方法画出的图如下,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);
=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

回复 27楼 cuixuepeng 的帖子

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

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

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

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

回复 地板 zhlong 的帖子

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

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

回复 34楼 michaelbeck 的帖子

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

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

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

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

回复 37楼 hnyanhua 的帖子

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);
=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
=tfrstft(Sig,(1:400)/SampFreq,M,h);
运行时出错

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

请高手解答

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

zhlong 发表于 2007-5-13 22:35 static/image/common/back.gif
下载tftb时频分析工具箱
或是留下邮箱,我发给你

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

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

华电机械 发表于 2013-5-20 15:56 static/image/common/back.gif
您有STFT和WVD的工具箱吗? 本人现在急用,,在此感激您能发一份给我。 邮箱:


见附件。

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

yghit08 发表于 2013-5-20 17:05 static/image/common/back.gif
见附件。

十分感激您啊,,有空向您请教!
页: 1 2 [3] 4
查看完整版本: 请教matlab高手,帮我看看这个STFT程序问题