声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1464|回复: 2

[HHT] 关于LMD的问题:极值点不足

[复制链接]
发表于 2013-7-22 17:43 | 显示全部楼层 |阅读模式

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

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

x
从论坛下载LMD代码,输入例子验证,发现错误,请高手指教。
fs=5000;
% x=@(t) (2+cos(90*t).*cos(500*t+1800.*t.*t));
t=0:1/fs:0.341;
% y=x(t);
y=cos(90*t)+cos(500*t);
subplot(5,1,1);plot(t,y);xlabel('原始信号');
[pf,a,si]=lmd1(y);
subplot(5,1,2);plot(t,pf(1,:));xlabel('PF1');
subplot(5,1,3);plot(t,pf(2,:));xlabel('PF2');
subplot(5,1,4);plot(t,pf(3,:));xlabel('PF3');
subplot(5,1,5);plot(t,pf(4,:));xlabel('残量信号')
输出结果显示:
??? Error using ==> boundary_conditions at 12
not enough extrema
Error in ==> envelope at 48
[tmin,tmax,xmin,xmax] = boundary_conditions(indmin,indmax,t,x,NBSYM);
Error in ==> zhaochun1 at 10
[envmin,envmax,envmoy,indmin,indmax,indzer] = envelope(t,h,'spline');
Error in ==> lmd1 at 18
    [pf,a,si]=zhaochun1(a,h,wucha1);
回复
分享到:

使用道具 举报

发表于 2013-7-26 10:27 | 显示全部楼层
仔细看看boundary_conditions函数,法国人的镜像延拓函数要求极值点数不能小于3.
if (length(indmin) + length(indmax) < 3)
error('not enough extrema')
所以分解的终止条件设置为极值点个数<3的话就会出现上述问题。

评分

1

查看全部评分

发表于 2013-8-17 16:41 | 显示全部楼层
什么是lmd啊?有相关的资料介绍吗?谢谢楼上

点评

搜索论坛, 局部均值分解(local mean decomposition)  发表于 2013-8-17 21:10
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-11 00:58 , Processed in 0.074137 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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