声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 3404|回复: 12

[HHT] 边际谱

[复制链接]
发表于 2009-10-28 10:59 | 显示全部楼层 |阅读模式

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

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

x
初涉EMD内容。在网上找到一些关于hht边际谱的内容,不大理解。用matlab如何求解边际谱?在论坛以前的帖子里找到一段小程序,有些地方不明白,想请教下大家:
clear;
fs=1000;    %fs为采样频率;
N=1000;     %采样点数(如何定?)
t=1/fs:1/fs:1;%
y1=2*sin(60*pi*t);
y2=5*sin(90*pi*t);
y=[y1;y2;zeros(size(y1))]; %IMF集(集合中为什么还要追加0行?)
%%%%%%%%%%%%%求边际谱
[A,fa,tt]=hhspectrum(y);%
[E,tt1]=toimage(A,fa,tt,length(tt));%
E=flipud(E);%(为什么要水平翻转?)
plot(tt1,E);
hold on;
for k=1:size(E,1)%
    bjp(k)=sum(E(k,:))*1/fs;%(bjp是什么量?)
end
f=(0:N-3)/N*(fs/2);
plot(f,bjp);
xlabel('频率 / Hz');
ylabel('幅值');
hold off;
回复
分享到:

使用道具 举报

 楼主| 发表于 2009-10-28 11:01 | 显示全部楼层
谢谢大家的答复
发表于 2009-11-2 12:07 | 显示全部楼层
 楼主| 发表于 2009-11-3 16:03 | 显示全部楼层

回复 板凳 bcyangbc 的帖子

这个帖子我也看过了,但是 还是没找到答案?谁能帮我一下,非常感激...
 楼主| 发表于 2009-11-4 10:07 | 显示全部楼层

【求助】hht边际谱问题

我用以下程序对一段语音进行emd分解,遂求边际谱,图如下,不知道对不对,哪有问题?这个是imf的边际谱吗?请高手指教
clear;
fs=8000;    %fs为采样频率;
t=1/fs:1/fs:1;
[s,fs]=wavread('0.5s/s1.wav');
imf=emd(s);
[A,fa,tt]=hhspectrum(imf);
[E,tt1]=toimage(A,fa,tt,length(tt));
E=flipud(E);
for k=1:size(E,1)
    bjp(k)=sum(E(k,:))*1/fs;
end
plot(bjp);
bjp.fig (32.91 KB, 下载次数: 14)
 楼主| 发表于 2009-11-4 10:10 | 显示全部楼层

hht边际谱问题

我用以下程序对一段语音进行emd分解,遂求边际谱,图如下,不知道对不对,哪有问题?这个是imf的边际谱吗?请高手指教
clear;
fs=8000;    %fs为采样频率;
t=1/fs:1/fs:1;
[s,fs]=wavread('0.5s/s1.wav');
imf=emd(s);
[A,fa,tt]=hhspectrum(imf);
[E,tt1]=toimage(A,fa,tt,length(tt));
E=flipud(E);
for k=1:size(E,1)
    bjp(k)=sum(E(k,:))*1/fs;
end
plot(bjp);
bjp.fig (32.91 KB, 下载次数: 54)

[ 本帖最后由 zhangnan3509 于 2009-12-15 09:23 编辑 ]
发表于 2009-12-11 00:48 | 显示全部楼层
本帖最后由 wdhd 于 2016-9-20 12:59 编辑

  你好:
  执行这一行都会出错,[E,tt1]=toimage(A,fa,tt,length(tt));
  请问是这个程序吗?谢
发表于 2009-12-15 15:49 | 显示全部楼层

回复 楼主 Rufter2007 的帖子

你这个是时频谱啊
 楼主| 发表于 2009-12-20 10:58 | 显示全部楼层

回复 板凳 leeting918 的帖子

?啊,...那该怎么求啊,谢谢指教,万分感激
发表于 2009-12-21 20:44 | 显示全部楼层

求助toimage.m文件

Rufter2007,你好,我刚刚接触hht,
我在执行这一行时出错,[E,tt1]=toimage(A,fa,tt,length(tt));
说是我的toimage.m文件中有错,你能给我传一份你的toimage.m文件吗?谢谢
 楼主| 发表于 2009-12-22 10:06 | 显示全部楼层

回复 5楼 chunhuajia 的帖子

提示什么错误呢?我这个文件也是下程序包是自带的喔。呵呵,大家共同学习下啦

[ 本帖最后由 Rufter2007 于 2009-12-22 10:08 编辑 ]

toimage.m

2.1 KB, 下载次数: 19

发表于 2009-12-22 10:34 | 显示全部楼层
运行这一行时:[E, tt1] = toimage(A,fa,tt,length(tt));
??? Error: File: toimage.m Line: 59 Column: 1
This statement is not inside any function.
(It follows the END that terminates the definition of
the function "toimage".)
      toimage.m中对应得行代码为:1t=length(t);
请教你一下,我刚刚学习这方面的知识,不是很懂。谢谢了
发表于 2009-12-22 10:39 | 显示全部楼层
我刚看了你的toimage.m文件了,跟我的一样呀,我也不知道为什么我的程序运行到这一行会出现错误,呵呵,帮我提点建议
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-6 07:51 , Processed in 0.065946 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表