大家来帮我看看我这个程序的hilbert谱和瞬时频率

longyi2003

这几个图的效果都好差啊，，各位大牛们，，能告诉我这是什么原因造成的吗？
程序如下：

1. clear;
   
2. N=3200;
   
3. t=linspace(0,1,N);
   
4. delta=t(2)-t(1); % 采样周期
   
5. fs=1/delta; % 采样频率
   
6. x=50*sin(2*pi*50*t);
   
7. y=15*sin(2*pi*150*t);
   
8. y1=8*sin(2*pi*250*t);
   
9. y2=5*sin(2*pi*350*t);
   
10. y3=2*sin(2*pi*450*t);
    
11. y4=1.8*sin(2*pi*550*t);
    
12. z=x+y+y1+y2+y3+y4;
    
13. % plot(t,z);
    
14. imf=emd1(z);
    
15. cemd_visu(z,1:length(z),imf);
    
16. 
    
17. [B,f,tt]=hhspectrum(imf);
    
18. [im,tt,Cenf]=toimage(B,f);
    
19. disp_hhs(im,[],fs);
    
20. colormap(flipud(gray)) %绘出Hilbert谱
    
21. 
    
22. %画出边际谱
    
23. %N=length(Cenf);%设置频率点数 %完全从理论公式出发。网格化后中心频率很重要，大家从连续数据变为离散的角度去思考，相信应该很容易理解
    
24. for k=1:size(im,1)
    
25. bjp(k)=sum(im(k,:))*1/fs;
    
26. end
    
27. figure;
    
28. plot(Cenf(1,:)*fs,bjp); % 作边际谱图 进行求取Hilbert谱时频率已经被抽样成具有一定窗长的离散频率，所以此时的频率轴已经是中心频率
    
29. xlabel('频率 / Hz');
    
30. ylabel('幅值');
    
31. 
    
32. s = size(imf);
    
33. k = s(1);
    
34. % figure
    
35. %imf分量幅度
    
36. for j=1:k
    
37. %imf1的Hilbert变换
    
38. xn(j,:)=hilbert(imf(j,:));
    
39. xr(j,:)=real(xn(j,:));
    
40. xi(j,:)=imag(xn(j,:));
    
41. A(j,:)=sqrt(xr(j,:).^2+xi(j,:).^2);
    
42. % subplot(k,1,j);plot(t,A(j,:));
    
43. % plot(t,xn(j,:));
    
44. % title(['imf',int2str(j)]);
    
45. % hold on;
    
46. end
    
47. % hold off;
    
48. %imf1的瞬时相位
    
49. % figure,
    
50. for j=1:k
    
51. P(j,:)=atan2(xi(j,:),xr(j,:));
    
52. % subplot(k,1,j);plot(t,P(j,:));
    
53. % title('imf瞬时相位')
    
54. end
    
55. 
    
56. 
    
57. %imf1瞬时频率
    
58. figure,
    
59. for j=1:k;
    
60. xh(j,:)=unwrap(P(j,:));
    
61. % fs=6400;
    
62. xhd(j,:)=fs*diff(xh(j,:))/(2*pi);
    
63. % subplot(k,1,j);
    
64. plot(t(1:N-1),xhd(j,:));title('imf瞬时频率')
    
65. 
    
66. hold on;
    
67. end
    
68. hold off

复制代码

shuihai707

模态混叠，产生的原因：1、信号中混有间断信号会导致EMD模态混叠，间断（间歇）信号是指高频小幅值信号；2、噪声（1和2可以看成一种情况）；3、几个频率共存的混合信号中,不能将某几个频率信号正确分离。显然你的问题是由第3种情况造成的。第3种情况要满足下述条件：

星星心

shuihai707 发表于 2013-5-8 17:08

                               
登录/注册后可看大图

模态混叠，产生的原因：1、信号中混有间断信号会导致EMD模态混叠，间断（间歇）信号是指高频小幅值信号；2、 ...

那楼主的这个问题需要怎么解决呀？

longyi2003

照这么说，HHT变换不能用在谐波检测上了？
还有个问题，下图我对信号直接做FFT变换后的频谱图，效果比HHT的频谱好多了，，基本上没误差。

shuihai707

longyi2003 发表于 2013-5-8 21:12

                               
登录/注册后可看大图

照这么说，HHT变换不能用在谐波检测上了？
还有个问题，下图我对信号直接做FFT变换后的频谱图，效果比HHT的 ...

说明你的信号是平稳信号，HHT主要是针对非平稳信号，如果是平稳信号就用FFT，效果还好，何必用HHT呢？不同方法具有不同的对象。

shuihai707

星星心 发表于 2013-5-8 20:36

                               
登录/注册后可看大图

那楼主的这个问题需要怎么解决呀？

暂时没有很好的解决办法，Gring提出用掩膜信号办法，EMD程序里有提到这个问题，论坛也有这方面的帖子，你可以找找，具体效果怎么样，自己试试吧，对于这类频率比较接近的信号，EMD不能很好的解决。

shuihai707

longyi2003 发表于 2013-5-8 21:12

                               
登录/注册后可看大图

照这么说，HHT变换不能用在谐波检测上了？
还有个问题，下图我对信号直接做FFT变换后的频谱图，效果比HHT的 ...

网上有些论文是用HHT用于谐波检测，你可以看看，他们是怎么解决的，我不做这方面，不是很了解。

longyi2003

shuihai707 发表于 2013-5-8 23:47

                               
登录/注册后可看大图

暂时没有很好的解决办法，Gring提出用掩膜信号办法，EMD程序里有提到这个问题，论坛也有这方面的帖子，你 ...

能不能用频带滤波把信号的频率分成几段，然后对各段进行emd分解。

shuihai707

longyi2003 发表于 2013-5-10 14:21

                               
登录/注册后可看大图

能不能用频带滤波把信号的频率分成几段，然后对各段进行emd分解。

试试吧，可能不行，分段越多，端点效应就越明显了，个人理解，可以尝试。

yangfengdehaizi

学习学习

jinhuaqiang

学习学习