设为首页收藏本站联系我们
切换到宽版

声振论坛

 找回密码
 我要加入

QQ登录

只需一步，快速开始

声振论坛»论坛 › 工程技术 › 信号处理方法 › 大家有没有能用的EMD程序
返回列表 发新帖
查看: 2759|回复: 11

[HHT] 大家有没有能用的EMD程序

[复制链接]
maxingking
发表于 2005-12-30 10:03 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 wdhd 于 2016-9-21 11:11 编辑

　　各位老大帮帮忙,网上的EMD程序不全,大家有没有能用的EMD程序发给我一份,maxingcellar@etang.com 谢谢大家了
　　[ 本帖最后由 zhangnan3509 于 2007-6-22 15:01 编辑 ]
回复
分享到：

使用道具 举报

tlmttt
发表于 2006-1-14 19:31 | 显示全部楼层
关注中．
回复 支持 反对

使用道具 举报

电光石火
发表于 2006-2-26 16:38 | 显示全部楼层
本帖最后由 wdhd 于 2016-9-21 11:11 编辑

　　同样关注中...
回复 支持 反对

使用道具 举报

dnj0306
发表于 2006-2-27 15:29 | 显示全部楼层

我也要！！

本帖最后由 wdhd 于 2016-9-21 11:12 编辑

　　我也要EMD程序!要有的话也给我发一个!dnj0306@163.com
回复 支持 反对

使用道具 举报

simon21
发表于 2006-2-27 20:48 | 显示全部楼层

回复：（maxingking）[求助]大家有没有能用的EMD程序

本帖最后由 wdhd 于 2016-9-21 11:12 编辑

　　http://vib.hit.edu.cn/vibbbs/dispbbs.asp?boardid=47&id=1355

　　四楼
回复 支持 反对

使用道具 举报

simon21
发表于 2006-2-27 21:14 | 显示全部楼层

回复：（maxingking）[求助]大家有没有能用的EMD程序

本帖最后由 wdhd 于 2016-9-21 11:12 编辑

　　其中缺少的extr试试下面的代码

　　%extr

　　function [indmin,indmax,indzer]=extr(x)

　　[DATAfile DATApath]=uigetfile('*.txt','输入信号');

　　watchon;

　　FILENAME=[DATApath,DATAfile];

　　x=load(FILENAME);

　　numzer=1;

　　nummin=1;

　　nummax=1;

　　for g=2:(length(x)-1)

　　if (x(g-1)*x(g)<0)|(x(g-1)*x(g)==0)&(x(g-1)==0)

　　indzer(numzer)=g-1;

　　numzer=numzer+1;

　　elseif (x(g-1)*x(g)<0)|(x(g-1)*x(g)==0)&(x(g)==0)

　　indzer(numzer)=g;

　　numzer=numzer+1;

　　elseif (x(g-1)*x(g)<0)|(x(g-1)*x(g)==0)&(x(g-1)~=0)&(x(g)~=0)

　　indzer(numzer)=(2*g-1)/2;

　　numzer=numzer+1;

　　end

　　if (x(g-1)>x(g))&(x(g)>x(g+1))

　　continue;

　　elseif (x(g-1)>x(g))&(x(g) indmin(nummin)=g;

　　nummin=nummin+1;

　　elseif (x(g-1) continue;

　　elseif (x(g-1)x(g+1))

　　indmax(nummax)=g;

　　nummax=nummax+1;

　　end

　　end

　　%subplot(2,3,1)

　　plot(indzer,x(indzer),'r-+')

　　title('indzer')

　　pause

　　% subplot(2,3,2)

　　plot(indmin,x(indmin),'b-o')

　　title('indmin')

　　pause

　　%subplot(2,3,3)

　　plot(indmax,x(indmax),'g-*')

　　title('indmax')

　　pause

　　%subplot(2,3,4)

　　plot(1:length(x),x,'k-s')

　　title('x')

　　pause

　　envmax = interp1(indmax,x(indmax),1:length(x),'spline');

　　envmin = interp1(indmin,x(indmin),1:length(x),'spline');

　　%subplot(2,3,5)

　　plot(1:length(x),envmin,'y-d')

　　title('envmin')

　　pause

　　%subplot(2,3,6)

　　plot(1:length(x),envmax,'c-h')

　　title('envmax')

　　pause

　　plot(1:length(x),mean((envmax+envmin)/2),'m-p')

　　title('mean((envmax+envmin)/2)')%extr

　　function [indmin,indmax,indzer]=extr(x)

　　[DATAfile DATApath]=uigetfile('*.txt','输入信号');

　　watchon;

　　FILENAME=[DATApath,DATAfile];

　　x=load(FILENAME);

　　numzer=1;

　　nummin=1;

　　nummax=1;

　　for g=2:(length(x)-1)

　　if (x(g-1)*x(g)<0)|(x(g-1)*x(g)==0)&(x(g-1)==0)

　　indzer(numzer)=g-1;

　　numzer=numzer+1;

　　elseif (x(g-1)*x(g)<0)|(x(g-1)*x(g)==0)&(x(g)==0)

　　indzer(numzer)=g;

　　numzer=numzer+1;

　　elseif (x(g-1)*x(g)<0)|(x(g-1)*x(g)==0)&(x(g-1)~=0)&(x(g)~=0)

　　indzer(numzer)=(2*g-1)/2;

　　numzer=numzer+1;

　　end

　　if (x(g-1)>x(g))&(x(g)>x(g+1))

　　continue;

　　elseif (x(g-1)>x(g))&(x(g) indmin(nummin)=g;

　　nummin=nummin+1;

　　elseif (x(g-1) continue;

　　elseif (x(g-1)x(g+1))

　　indmax(nummax)=g;

　　nummax=nummax+1;

　　end

　　end

　　%subplot(2,3,1)

　　plot(indzer,x(indzer),'r-+')

　　title('indzer')

　　pause

　　% subplot(2,3,2)

　　plot(indmin,x(indmin),'b-o')

　　title('indmin')

　　pause

　　%subplot(2,3,3)

　　plot(indmax,x(indmax),'g-*')

　　title('indmax')

　　pause

　　%subplot(2,3,4)

　　plot(1:length(x),x,'k-s')

　　title('x')

　　pause

　　envmax = interp1(indmax,x(indmax),1:length(x),'spline');

　　envmin = interp1(indmin,x(indmin),1:length(x),'spline');

　　%subplot(2,3,5)

　　plot(1:length(x),envmin,'y-d')

　　title('envmin')

　　pause

　　%subplot(2,3,6)

　　plot(1:length(x),envmax,'c-h')

　　title('envmax')

　　pause

　　plot(1:length(x),mean((envmax+envmin)/2),'m-p')

　　title('mean((envmax+envmin)/2)')%extr

　　function [indmin,indmax,indzer]=extr(x)

　　[DATAfile DATApath]=uigetfile('*.txt','输入信号');

　　watchon;

　　FILENAME=[DATApath,DATAfile];

　　x=load(FILENAME);

　　numzer=1;

　　nummin=1;

　　nummax=1;

　　for g=2:(length(x)-1)

　　if (x(g-1)*x(g)<0)|(x(g-1)*x(g)==0)&(x(g-1)==0)

　　indzer(numzer)=g-1;

　　numzer=numzer+1;

　　elseif (x(g-1)*x(g)<0)|(x(g-1)*x(g)==0)&(x(g)==0)

　　indzer(numzer)=g;

　　numzer=numzer+1;

　　elseif (x(g-1)*x(g)<0)|(x(g-1)*x(g)==0)&(x(g-1)~=0)&(x(g)~=0)

　　indzer(numzer)=(2*g-1)/2;

　　numzer=numzer+1;

　　end

　　if (x(g-1)>x(g))&(x(g)>x(g+1))

　　continue;

　　elseif (x(g-1)>x(g))&(x(g) indmin(nummin)=g;

　　nummin=nummin+1;

　　elseif (x(g-1) continue;

　　elseif (x(g-1)x(g+1))

　　indmax(nummax)=g;

　　nummax=nummax+1;

　　end

　　end

　　%subplot(2,3,1)

　　plot(indzer,x(indzer),'r-+')

　　title('indzer')

　　pause

　　% subplot(2,3,2)

　　plot(indmin,x(indmin),'b-o')

　　title('indmin')

　　pause

　　%subplot(2,3,3)

　　plot(indmax,x(indmax),'g-*')

　　title('indmax')

　　pause

　　%subplot(2,3,4)

　　plot(1:length(x),x,'k-s')

　　title('x')

　　pause

　　envmax = interp1(indmax,x(indmax),1:length(x),'spline');

　　envmin = interp1(indmin,x(indmin),1:length(x),'spline');

　　%subplot(2,3,5)

　　plot(1:length(x),envmin,'y-d')

　　title('envmin')

　　pause

　　%subplot(2,3,6)

　　plot(1:length(x),envmax,'c-h')

　　title('envmax')

　　pause

　　plot(1:length(x),mean((envmax+envmin)/2),'m-p')

　　title('mean((envmax+envmin)/2)')%extr

　　function [indmin,indmax,indzer]=extr(x)

　　[DATAfile DATApath]=uigetfile('*.txt','输入信号');

　　watchon;

　　FILENAME=[DATApath,DATAfile];

　　x=load(FILENAME);

　　numzer=1;

　　nummin=1;

　　nummax=1;

　　for g=2:(length(x)-1)

　　if (x(g-1)*x(g)<0)|(x(g-1)*x(g)==0)&(x(g-1)==0)

　　indzer(numzer)=g-1;

　　numzer=numzer+1;

　　elseif (x(g-1)*x(g)<0)|(x(g-1)*x(g)==0)&(x(g)==0)

　　indzer(numzer)=g;

　　numzer=numzer+1;

　　elseif (x(g-1)*x(g)<0)|(x(g-1)*x(g)==0)&(x(g-1)~=0)&(x(g)~=0)

　　indzer(numzer)=(2*g-1)/2;

　　numzer=numzer+1;

　　end

　　if (x(g-1)>x(g))&(x(g)>x(g+1))

　　continue;

　　elseif (x(g-1)>x(g))&(x(g) indmin(nummin)=g;

　　nummin=nummin+1;

　　elseif (x(g-1) continue;

　　elseif (x(g-1)x(g+1))

　　indmax(nummax)=g;

　　nummax=nummax+1;

　　end

　　end

　　%subplot(2,3,1)

　　plot(indzer,x(indzer),'r-+')

　　title('indzer')

　　pause

　　% subplot(2,3,2)

　　plot(indmin,x(indmin),'b-o')

　　title('indmin')

　　pause

　　%subplot(2,3,3)

　　plot(indmax,x(indmax),'g-*')

　　title('indmax')

　　pause

　　%subplot(2,3,4)

　　plot(1:length(x),x,'k-s')

　　title('x')

　　pause

　　envmax = interp1(indmax,x(indmax),1:length(x),'spline');

　　envmin = interp1(indmin,x(indmin),1:length(x),'spline');

　　%subplot(2,3,5)

　　plot(1:length(x),envmin,'y-d')

　　title('envmin')

　　pause

　　%subplot(2,3,6)

　　plot(1:length(x),envmax,'c-h')

　　title('envmax')

　　pause

　　plot(1:length(x),mean((envmax+envmin)/2),'m-p')

　　title('mean((envmax+envmin)/2)')%extr

　　function [indmin,indmax,indzer]=extr(x)

　　[DATAfile DATApath]=uigetfile('*.txt','输入信号');

　　watchon;

　　FILENAME=[DATApath,DATAfile];

　　x=load(FILENAME);

　　numzer=1;

　　nummin=1;

　　nummax=1;

　　for g=2:(length(x)-1)

　　if (x(g-1)*x(g)<0)|(x(g-1)*x(g)==0)&(x(g-1)==0)

　　indzer(numzer)=g-1;

　　numzer=numzer+1;

　　elseif (x(g-1)*x(g)<0)|(x(g-1)*x(g)==0)&(x(g)==0)

　　indzer(numzer)=g;

　　numzer=numzer+1;

　　elseif (x(g-1)*x(g)<0)|(x(g-1)*x(g)==0)&(x(g-1)~=0)&(x(g)~=0)

　　indzer(numzer)=(2*g-1)/2;

　　numzer=numzer+1;

　　end

　　if (x(g-1)>x(g))&(x(g)>x(g+1))

　　continue;

　　elseif (x(g-1)>x(g))&(x(g) indmin(nummin)=g;

　　nummin=nummin+1;

　　elseif (x(g-1) continue;

　　elseif (x(g-1)x(g+1))

　　indmax(nummax)=g;

　　nummax=nummax+1;

　　end

　　end

　　%subplot(2,3,1)

　　plot(indzer,x(indzer),'r-+')

　　title('indzer')

　　pause

　　% subplot(2,3,2)

　　plot(indmin,x(indmin),'b-o')

　　title('indmin')

　　pause

　　%subplot(2,3,3)

　　plot(indmax,x(indmax),'g-*')

　　title('indmax')

　　pause

　　%subplot(2,3,4)

　　plot(1:length(x),x,'k-s')

　　title('x')

　　pause

　　envmax = interp1(indmax,x(indmax),1:length(x),'spline');

　　envmin = interp1(indmin,x(indmin),1:length(x),'spline');

　　%subplot(2,3,5)

　　plot(1:length(x),envmin,'y-d')

　　title('envmin')

　　pause

　　%subplot(2,3,6)

　　plot(1:length(x),envmax,'c-h')

　　title('envmax')

　　pause

　　plot(1:length(x),mean((envmax+envmin)/2),'m-p')

　　title('mean((envmax+envmin)/2)')%extr

　　function [indmin,indmax,indzer]=extr(x)

　　[DATAfile DATApath]=uigetfile('*.txt','输入信号');

　　watchon;

　　FILENAME=[DATApath,DATAfile];

　　x=load(FILENAME);

　　numzer=1;

　　nummin=1;

　　nummax=1;

　　for g=2:(length(x)-1)

　　if (x(g-1)*x(g)<0)|(x(g-1)*x(g)==0)&(x(g-1)==0)

　　indzer(numzer)=g-1;

　　numzer=numzer+1;

　　elseif (x(g-1)*x(g)<0)|(x(g-1)*x(g)==0)&(x(g)==0)

　　indzer(numzer)=g;

　　numzer=numzer+1;

　　elseif (x(g-1)*x(g)<0)|(x(g-1)*x(g)==0)&(x(g-1)~=0)&(x(g)~=0)

　　indzer(numzer)=(2*g-1)/2;

　　numzer=numzer+1;

　　end

　　if (x(g-1)>x(g))&(x(g)>x(g+1))

　　continue;

　　elseif (x(g-1)>x(g))&(x(g) indmin(nummin)=g;

　　nummin=nummin+1;

　　elseif (x(g-1) continue;

　　elseif (x(g-1)x(g+1))

　　indmax(nummax)=g;

　　nummax=nummax+1;

　　end

　　end

　　%subplot(2,3,1)

　　plot(indzer,x(indzer),'r-+')

　　title('indzer')

　　pause

　　% subplot(2,3,2)

　　plot(indmin,x(indmin),'b-o')

　　title('indmin')

　　pause

　　%subplot(2,3,3)

　　plot(indmax,x(indmax),'g-*')

　　title('indmax')

　　pause

　　%subplot(2,3,4)

　　plot(1:length(x),x,'k-s')

　　title('x')

　　pause

　　envmax = interp1(indmax,x(indmax),1:length(x),'spline');

　　envmin = interp1(indmin,x(indmin),1:length(x),'spline');

　　%subplot(2,3,5)

　　plot(1:length(x),envmin,'y-d')

　　title('envmin')

　　pause

　　%subplot(2,3,6)

　　plot(1:length(x),envmax,'c-h')

　　title('envmax')

　　pause

　　plot(1:length(x),mean((envmax+envmin)/2),'m-p')

　　title('mean((envmax+envmin)/2)')
回复 支持 反对

使用道具 举报

zfw
发表于 2006-3-23 20:33 | 显示全部楼层
本帖最后由 wdhd 于 2016-9-21 11:12 编辑

　　各位兄弟好，请问谁有EMD的matlab程序，本人也在搞这个

　　想要个程序参考参考，多谢!my123_@sohu.com
回复 支持 反对

使用道具 举报

sealily
发表于 2006-3-24 14:06 | 显示全部楼层
我是个新手，请问高手，是不是要把那个缺少的extr代码复制到原来的EMD程序中才能使用？<BR>多谢了！<BR>还有，如果哪位高手有了完整的程序了，请发给我一个<BR>多谢！rzhang_838@163.com
回复 支持 反对

使用道具 举报

liugen_666
发表于 2006-3-24 19:57 | 显示全部楼层

仁兄!!!

本帖最后由 wdhd 于 2016-9-21 11:12 编辑

　　各位兄弟!!!

　　我也想要一个EMD程序嗄!!!

　　谢谢
回复 支持 反对

使用道具 举报

simon21
发表于 2006-3-25 21:19 | 显示全部楼层

回复：（sealily）我是个新手，请问高手，是不是要把那...

本帖最后由 wdhd 于 2016-9-21 11:13 编辑

　　以下是引用sealily在2006-3-24 14:06:35的发言：

　　我是个新手，请问高手，是不是要把那个缺少的extr代码复制到原来的EMD程序中才能使用?

　　多谢了!

　　还有，如果哪位高手有了完整的程序了，请发给我一个

　　多谢!rzhang_838@163.com

　　是的，建立extr.m文件，文件内容上面已经贴了
回复 支持 反对

使用道具 举报

snowwa
发表于 2006-4-3 09:28 | 显示全部楼层
本帖最后由 wdhd 于 2016-9-21 11:13 编辑

　　elseif (x(g-1)>x(g))&(x(g) indmin(nummin)=g;

　　nummin=nummin+1;

　　elseif (x(g-1) continue;

　　elseif (x(g-1)x(g+1))

　　indmax(nummax)=g;

　　nummax=nummax+1;

　　end

　　end

　　%subplot(2,3,1)

　　这段代码有问题啊
回复 支持 反对

使用道具 举报

wht0211
发表于 2010-7-30 15:29 | 显示全部楼层
emd.m急需中
回复 支持 反对

使用道具 举报

返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-26 02:35 , Processed in 0.075901 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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