一个EMD的程序，大家帮我看看结果对不对？

汀小汀

本人初学EMD
程序是根据论坛里面的一个程序改动的
仿真信号为x(t)=2sin(30πt)+4sin(20πt)sin(2πt/10)+sin(10πt)
请大家指教下，先谢了！  

1. 
   
2. fs=1000;
   
3. l=1;
   
4. t=1/fs:1/fs:l;
   
5. N=length(t);
   
6. x=2*sin(2*pi*15*t)+4*sin(2*pi*10*t).*sin(2*pi*t/10)+sin(2*pi*5*t);
   
7. plot(t,x);
   
8. xlabel('t/s');ylabel('x(t)');
   
9. imf=emd(x);
   
10. cemd_visu(x,1:length(x),imf);
    
11. xn1=hilbert(imf(1,:));
    
12. xr1=real(xn1);
    
13. xi1=imag(xn1);
    
14. A1=sqrt(xr1.^2+xi1.^2);
    
15. figure,subplot(3,1,1);plot(t,A1);
    
16. xlabel('t/s');ylabel('瞬时振幅');title('imf1');
    
17. xn2=hilbert(imf(2,:));
    
18. xr2=real(xn2);
    
19. xi2=imag(xn2);
    
20. A2=sqrt(xr2.^2+xi2.^2);
    
21. subplot(3,1,2);plot(t,A2);
    
22. xlabel('t/s');ylabel('瞬时振幅');title('imf2');
    
23. xn3=hilbert(imf(3,:));
    
24. xr3=real(xn3);
    
25. xi3=imag(xn3);
    
26. A3=sqrt(xr3.^2+xi3.^2);
    
27. subplot(3,1,3);plot(t,A3);
    
28. xlabel('t/s');ylabel('瞬时振幅');title('imf3');
    
29. P1=atan2(xi1,xr1);
    
30. xh1=unwrap(P1);
    
31. fs=1000;
    
32. xhd1=fs*diff(xh1)/(2*pi);
    
33. figure,subplot(3,1,1);plot(t(1:999),xhd1);
    
34. xlabel('t/s');ylabel('瞬时频率');title('imf1');
    
35. P2=atan2(xi2,xr2);
    
36. xh2=unwrap(P2);
    
37. fs=1000;
    
38. xhd2=fs*diff(xh2)/(2*pi);
    
39. subplot(3,1,2);plot(t(1:999),xhd2);
    
40. xlabel('t/s');ylabel('瞬时频率');title('imf2');
    
41. P3=atan2(xi3,xr3);
    
42. xh3=unwrap(P3);
    
43. fs=1000;
    
44. xhd3=fs*diff(xh3)/(2*pi);
    
45. subplot(3,1,3);plot(t(1:999),xhd3);
    
46. xlabel('t/s');ylabel('瞬时频率');title('imf3');
    
47. [A,f,t]=hhspectrum(imf);
    
48. [E,t,cenf]=toimage(A,f);
    
49. disp_hhs(E,[],1000);
    
50. for k=1:size(E,1)
    
51.     bjp(k)=sum(E(k,:))*1/l*1/fs;
    
52. end
    
53. H=size(E,1);
    
54. f=(0:H-1)/H*(fs/2);
    
55. figure();
    
56. plot(f,bjp);
    
57. xlabel('频率/HZ');
    
58. ylabel('幅值');title('边际谱')
    

复制代码

杨德昌

哥们，你首先把刻度值标上啊！

汀小汀

回复 2 # 杨德昌 的帖子

不好意思啊，马上改下图

汀小汀

回复 2 # 杨德昌 的帖子

汀小汀

x(t)

imf

瞬时幅值

瞬时频率

Hilbert谱

边际谱