wang-wei 发表于 2008-7-26 16:15

一个有关边际谱的问题

前几天再论坛上看到一个有关边际谱的程序，有的语句不太懂，请高手指点一下
clear;
fs=1000;    %fs为采样频率;tspan=2;
t=1/fs:1/fs:tspan;N=length(t);   %采样点数y1=5*sin(2*pi*241*t);
y2=3*sin(2*pi*73*t);
y=; %IMF集
%%%%%%%%%%%%%求边际谱
=hhspectrum(y);
=toimage(A,fa,tt,length(tt));
E=flipud(E)for k=1:size(E,1)
    bjp(k)=sum(E(k,:))*1/fs;
end
f=(0:N-3)/N*(fs/2);
plot(f,bjp);
xlabel('频率 / Hz');
ylabel('幅值'); =toimage(A,fa,tt,length(tt));这句是什么意思，其中的E是不是指toimage.m中的im, E(k,:)又是什么，为什么边际谱就是bjp(k)=sum(E(k,:))*1/fs

civil1314 发表于 2008-10-15 13:41

同样的问题。E指的是什么？toimage。m的输出没有E啊？
E是指im还是ff呢？不是很明白。
边际谱是对hilbert谱的时间积分，程序怎么体现出来的？
不明白。
望高手指点一下

civil1314 发表于 2008-10-19 22:42

自己给自己解答！呵呵。。

仔细看看这个帖子把！
http://forum.vibunion.com/forum/viewthread.php?tid=50010

wang-wei 发表于 2008-10-23 11:09

我觉得边际谱就是bjp(k)=sum(E(k,:)) 指某一频率对应的振幅总和，而不是bjp(k)=sum(E(k,:))*1/fs ; E就是im
另外频率是f0=(0:size(im,1)-1)*fs/size(im,1)+1/16;
法国人的程序中im是用了一个很好的技术得出来的，可以看子程序 toimage和DISP_HHS:@D

[ 本帖最后由 wang-wei 于 2008-10-23 11:12 编辑 ]

civil1314 发表于 2008-10-23 16:39

时频工具箱

振动论坛里面有的。你再找找。

civil1314 发表于 2008-10-23 16:43

原帖由 wang-wei 于 2008-10-23 11:09 发表 http://www.chinavib.com/forum/images/common/back.gif
我觉得边际谱就是bjp(k)=sum(E(k,:)) 指某一频率对应的振幅总和，而不是bjp(k)=sum(E(k,:))*1/fs ; E就是im
另外频率是f0=(0:size(im,1)-1)*fs/size(im,1)+1/16;
法国人的程序中im是用了一个很好的技术得出来的， ...
可以把你说的法国的那个im传上来么？
另外而不是bjp(k)=sum(E(k,:))*1/fs 是要积分用的啊。

wang-wei 发表于 2008-10-25 14:20

im和这的E是一样的，不过表示不同
=toimage(A,fa,tt,length(tt));即=toimage(A,fa,tt,length(tt));

cyp1967 发表于 2014-6-27 06:43

谢谢指点！{:{39}:}
页: [1]
查看完整版本: 一个有关边际谱的问题