u6.m

来自「EMD方法的一個例子」· M 代码 · 共 38 行

a=0:pi:2*pi; b=-1.5:1.5:1.5; p=[0.05,0.8,0.9,0.18];
plct(x,y,a,b,p,1,1,'ro',0); hold on;
plot(xx,yy,'b')
hold off;
p(2)=0.6;
local=[0,2*pi];
plct(xx,C(1,:),a,b,p,1,1,'b-',0);hold on;
[maxx,maxy,minx,miny]=find_max_min(xx,C(1,:));
[maxxx,maxyy]=generalize(maxx,maxy,local);
[minxx,minyy]=generalize(minx,miny,local);
plot(maxxx,maxyy,'r--');
plot(minxx,minyy,'r--');
hold off;
p(2)=0.4;
plct(xx,C(2,:),a,b,p,1,1,'b-',0);hold on;
[maxx,maxy,minx,miny]=find_max_min(xx,C(2,:));
[maxxx,maxyy]=generalize(maxx,maxy,local);
[minxx,minyy]=generalize(minx,miny,local);
plot(maxxx,maxyy,'r--');
plot(minxx,minyy,'r--');
hold off;
p(2)=0.2;
plct(xx,C(3,:),a,b,p,1,1,'b-',0);hold on;
[maxx,maxy,minx,miny]=find_max_min(xx,C(3,:));
[maxxx,maxyy]=generalize(maxx,maxy,local);
[minxx,minyy]=generalize(minx,miny,local);
plot(maxxx,maxyy,'r--');
plot(minxx,minyy,'r--');
hold off;
p(2)=0;
plct(xx,C(4,:),a,b,p,1,1,'b-',0);hold on;
[maxx,maxy,minx,miny]=find_max_min(xx,C(4,:));
[maxxx,maxyy]=generalize(maxx,maxy,local);
[minxx,minyy]=generalize(minx,miny,local);
plot(maxxx,maxyy,'r--');
plot(minxx,minyy,'r--');
hold off;