次端点镜像延拓问题

zt861217

这里有个次端点镜像延拓的程序，不知道怎么用

1. function imf = TryEmd(ip)
   
2. %--------------------------------------------------------------------------
   
3. %  
   
4. N=length(ip);
   
5. x=1:N;
   
6. [max,max_No,min,min_No]=MaxMin(ip);
   
7. 
   
8. %--------------------------------------------------------------------------
   
9. %----函数说明(152)：
   
10. %    输入ip--待分解数据序列，一维列向量；
    
11. %    输出imf--分解后得到的IMF分量，二维数组表示；
    
12. %--------------------------------------------------------------------------
    
13. %----调用函数：
    
14. %    Enovelope--次端点镜像法求解上、下包络线；(39)
    
15. %    MaxMin--求局部极值点及其位置；(24)
    
16. %    MaxMax--求数组极大值;(6)
    
17. %    ZeroNum--求过零点数目；(17)
    
18. %--------------------------------------------------------------------------
    
19. %----参数设置：
    
20.      PRE_yuzhi=0.005; %--偏差阈值；
    
21.      NUM1=5; %--分解得到IMF分量的最多数目；
    
22.      NUM2=800; %--提取IMF分量的最多筛分次数；
    
23. %--------------------------------------------------------------------------
    
24. flag=0;
    
25. I=0;
    
26. num1=0;
    
27. num2=0;
    
28. while(max_No(1)~=0 && min_No(1)~=0)
    
29.     if(num1>=NUM1) break; end
    
30.     while(1)
    
31.         num2=num2+1;
    
32.         [yy1,yy2]=enovelope(x,max,max_No,min,min_No);
    
33.         mean=(yy1+yy2)/2;
    
34.         h1=ip-mean;
    
35.         [max,max_No,min,min_No]=MaxMin(h1);
    
36.         if(max_No(1)==0 || min_No(1)==0)
    
37.             flag=1;
    
38.             break;
    
39.         end
    
40.         [yy1,yy2]=enovelope(x,max,max_No,min,min_No);
    
41.         mean=(yy1+yy2)/2;
    
42.         hmax=MaxMax(mean);
    
43.         N1=ZeroNum(h1);
    
44.         tiaojian1=hmax-PRE_yuzhi;
    
45.         tiaojian2=abs(length(max_No)+length(min_No)-N1);
    
46.         if((tiaojian1<0 && tiaojian2<=1) || num2<=NUM2)
    
47.             I=I+1;
    
48.             imf(I,:)=h1;
    
49.             break;
    
50.         else
    
51.             ip1=h1;
    
52.             [max,max_No,min,min_No]=MaxMin(ip1);
    
53.         end
    
54.     end
    
55.    
    
56.     if(flag==1)
    
57.         flag=0;
    
58.         break;
    
59.     else
    
60.         ip=ip-imf(I,:);
    
61.         num1=I;
    
62.         [max,max_No,min,min_No]=MaxMin(ip);
    
63.     end
    
64. end
    
65. imf(I+1,:)=ip;
    
66. end
    
67. %%%%%%%%%%%%%%
    
68. 
    
69. function [y_up,y_low]= Enovelope(x,max,max_No,min,min_No)
    
70. %--次端点镜像延拓法获得的上下包络线；
    
71. 
    
72. l=length(x);
    
73. l1=length(max);
    
74. N1=max_No(1)-1;
    
75. N2=length(x)-max_No(end);
    
76. 
    
77. l2=length(min);
    
78. N3=min_No(1)-1;
    
79. N4=length(x)-min_No(end);
    
80. max1_No(1)=1;
    
81. for i=2:(l1+1)
    
82.     max1_No(i)=max_No(i-1)+N1;
    
83. end
    
84. max1_No(l1+2)=l+N1+N2;
    
85. max1(1)=max(1);
    
86. for i=2:(l1+1)
    
87.     max1(i)=max(i-1);
    
88. end
    
89. max1(l1+2)=max(end);
    
90. min1_No(1)=1;
    
91. for i=2:(l2+1)
    
92.     min1_No(i)=min_No(i-1)+N3;
    
93. end
    
94. min1_No(l2+2)=l+N3+N4;
    
95. min1(1)=min(1);
    
96. 
    
97. for i=2:(l2+1)
    
98.     min1(i)=min(i-1);
    
99. end
    
100. min1(l2+2)=min(end);
     
101. 
     
102. x1=1:(l+N1+N2);
     
103. cs1=spline(max1_No,max1);%三次样条差值
     
104. y1=ppval(cs1,x1);%求出插值曲线
     
105. for i=1:l
     
106.     y_up(i)=y1(i+N1);
     
107. end
     
108. 
     
109. x2=1:(l+N3+N4);
     
110. cs2=spline(min1_No,min1);%三次样条差值
     
111. y2=ppval(cs2,x2);%求出插值曲线
     
112. for
     
113.     i=1:l
     
114.     y_low(i)=y2(i+N3);
     
115. end
     
116. end

复制代码