一个有关边际谱的问题

wang-wei

前几天再论坛上看到一个有关边际谱的程序，有的语句不太懂，请高手指点一下

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=[y1;y2;zeros(size(y1))]; %IMF集
%%%%%%%%%%%%%求边际谱
[A,fa,tt]=hhspectrum(y);
[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
f=(0:N-3)/N*(fs/2);
plot(f,bjp);
xlabel('频率 / Hz');
ylabel('幅值');

[E,tt1]=toimage(A,fa,tt,length(tt));这句是什么意思，其中的E是不是指toimage.m中的im, E(k,:)又是什么，为什么

边际谱就是bjp(k)=sum(E(k,:))*1/fs

civil1314

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

civil1314

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

wang-wei

我觉得边际谱就是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

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

civil1314

原帖由 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是用了一个很好的技术得出来的， ...

可以把你说的法国的那个im传上来么？
另外而不是bjp(k)=sum(E(k,:))*1/fs 是要积分用的啊。

wang-wei

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

cyp1967

谢谢指点！