|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
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);
|
|