还是边际谱问题 请高手帮助！

toniling310

我是名正在做毕业设计的大四学生 ，题目是hht在功率密度谱估计的应用，以前完全没有接触过hht，所以很多很多问题都不懂，现在做到边际谱了 程序很简单 就是把破凰的程序改了一下clear;
fs=1000;   
N=1000;     
t=1/fs:1/fs:1;
xn = sin(2*pi*80*t)+3*sin(2*pi*140*t) ;
imf=emd(xn);
[A,fa,tt]=hhspectrum(imf(1:end-1,:));
[E,tt1]=toimage(A,fa,tt,length(tt));
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('幅值');
但是出来的图却是这样:

这是怎么回事，我觉得谱峰应该很窄而且比较平滑才对，
还有边际谱最后能不能得到功率谱呢 ？

Rovis

主要的原因不是因为边际谱本身出的问题，而是EMD分解出来的IMF的质量有问题，现在正在尝试解决中

toniling310

我的直接就用了emd（xn）这个语句 是不是这样不妥呢？