声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1362|回复: 4

[综合] 一个EMD的程序,大家帮我看看结果对不对?

[复制链接]
发表于 2011-5-10 12:05 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 汀小汀 于 2011-5-10 12:19 编辑

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

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

x(t)

x(t)

imf

imf

瞬时幅值

瞬时幅值

瞬时频率

瞬时频率

Hilbert谱

Hilbert谱

边际谱

边际谱

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 2011-5-10 17:30 | 显示全部楼层
哥们,你首先把刻度值标上啊!
 楼主| 发表于 2011-5-10 22:27 | 显示全部楼层
回复 2 # 杨德昌 的帖子

不好意思啊,马上改下图
 楼主| 发表于 2011-5-10 22:37 | 显示全部楼层
回复 2 # 杨德昌 的帖子

 楼主| 发表于 2011-5-10 22:43 | 显示全部楼层

x(t)

x(t)

imf

imf

瞬时幅值

瞬时幅值

瞬时频率

瞬时频率

Hilbert谱

Hilbert谱

边际谱

边际谱

_(1RK87$ZKBT8CWLLMOP74W.jpg
}(X~V5X~]Q8D`BG)A3I)WFC.jpg
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-14 15:23 , Processed in 0.081305 second(s), 23 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表