声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1090|回复: 1

[HHT] 怎么不能运行啊?

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

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

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

x
前几天在这上下载了个名为几个函数的文件,发现里面的程序运行不了啊.原程序如下:
% dt=0.001;
% T=1;
% n=T/dt;
% t=0:dt:T;
% xt=2*sin(2*pi*15*t)+4*sin(2*pi*10*t).*sin(2*pi*t/10)+sin(2*pi*5*t);
%
原始信号
function imf=TryEmdt(h1,t)

flag=0;

ll=0;

xtT=h1;

while max(abs(h1))>0.5

uu=0;

while(1)

uu=uu+1

[bg,bt,cg,ct,N2]=MaxMin(h1,t);

l1=length(bt);

l2=length(ct);
   

if (l1==1)||(l2==1)

flag=1;
%
只有一个极值点的情况

break

end
   

N1=ZeroNum(h1);

[By,Cy]= Enovelope(bt,bg,ct,cg,t);

Avg=(By+Cy)/2;
%
极大极小平均

h=h1-Avg;

if ((abs(N2-N1)<=1)&&(max(abs(Avg))<0.001))||(uu>800)

ll=ll+1

imf(ll,:)=h1;

break

else

h1=h;

end

end

if (flag==1)

flag=0;


break

else

xtt=xtT-imf(ll,:);

h1=xtt;

xtT=h1;

end

end

imf(ll+1,:)=h1;
%
残余模态
显示h1没有定义是怎么回事啊?
回复
分享到:

使用道具 举报

 楼主| 发表于 2011-10-20 16:23 | 显示全部楼层
回复 1 # 李清志 的帖子

还有里面的N2是什么?Enovelope函数是什么作用啊?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-25 18:17 , Processed in 0.061461 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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