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

简谐波的边际谱

做了个简单的边际谱算例，
参考了论坛的相关帖子，
但是边际谱的结果不是太好。
请指点
谢谢。
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

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

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()
xlabel('频率/Hz')
ylabel('幅值')
title('边际谱')

   

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

新版EMD包？07年的还是？

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

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

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

好多内容都没有办法看，是为什么呀？
页: [1]
查看完整版本: 简谐波的边际谱