简谐波的边际谱

linqus · 发表于 2010-10-24 01:09

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

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

x

做了个简单的边际谱算例，
参考了论坛的相关帖子，
但是边际谱的结果不是太好。
请指点
谢谢。
ps：用的是新版emd包

linqus · 发表于 2010-10-24 14:55

直接贴出m文件吧。
上面的帖子没法修改。
clc;
clear all;
fa=10;
fb=20;
a=1;
b=1;
dt=0.002;
fs=1/dt;
t=0:dt:5;
tspan=dt*(length(t)-1);
pi=3.1415926;
x=a*sin(2*pi*fa*t)+b*sin(2*pi*fb*t);
imf=emd(x);
idx=size(imf);
han=idx(1)+1;
subplot(han,1,1)
plot(t,x)
ylabel('original')
%++++++++++++++++++++++++++++++++
for j=1:han-1
subplot(han,1,j+1)
plot(t,imf(j,:))
ylabel('imf')
end

[A,f,tt]=hhspectrum(imf); %HHT 时频谱计算
[im,tt]=toimage(A,f,tt,length(tt));
disp_hhs(im,[],fs);
ylim([0,30])

NN=size(im,1);
for k=1:NN
bjp(k)=sum(im(k,:))*1/fs;
end
fbjp=(0:NN-1)/NN*(fs/2);
%----------------------边际谱--------------------------
figure(3)
%bjp=bjp/max(bjp);
plot(fbjp,bjp)
xlim([0,30])
xlabel('频率/Hz')
ylabel('幅值')
title('边际谱')

Cena马 · 发表于 2010-10-24 15:27

新版EMD包？07年的还是？

linqus · 发表于 2010-10-27 22:35

用的是07年的emd包，
不知哪位可以看看这个小问题出在哪里，谢谢了

hub115 · 发表于 2011-1-7 16:15

问题不是由边际谱的程序造成的，而是来自EMD过程本身有一些问题

3060104124 · 发表于 2011-12-4 11:42

新手学习ING

zhangnan3509 · 发表于 2011-12-4 11:48

本帖最后由 zhangnan3509 于 2011-12-4 11:53 编辑

其实效果还算是不错的了，我没看你设定的两个叠加量，从IMF的情况来说，尤其是从HHT谱上看在20附近谱线都散了，10附近略好一些，所以边际谱上不可能得到特别好的效果，能成这算是不错了，这不是边际谱的问题

孙雅青 · 发表于 2015-3-11 17:33

好多内容都没有办法看，是为什么呀？

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

[HHT] 简谐波的边际谱

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