<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta charset="utf-8">
<meta name="generator" content="pdf2htmlEX">
<meta http-equiv="X-UA-Compatible" content="IE=edge,chrome=1">
<link rel="stylesheet" href="https://static.pudn.com/base/css/base.min.css">
<link rel="stylesheet" href="https://static.pudn.com/base/css/fancy.min.css">
<link rel="stylesheet" href="https://static.pudn.com/prod/directory_preview_static/624f338f74bc5c010528612d/raw.css">
<script src="https://static.pudn.com/base/js/compatibility.min.js"></script>
<script src="https://static.pudn.com/base/js/pdf2htmlEX.min.js"></script>
<script>
try{
pdf2htmlEX.defaultViewer = new pdf2htmlEX.Viewer({});
}catch(e){}
</script>
<title></title>
</head>
<body>
<div id="sidebar" style="display: none">
<div id="outline">
</div>
</div>
<div id="pf1" class="pf w0 h0" data-page-no="1"><div class="pc pc1 w0 h0"><img class="bi x0 y0 w1 h1" alt="" src="https://static.pudn.com/prod/directory_preview_static/624f338f74bc5c010528612d/bg1.jpg"><div class="t m0 x1 h2 y1 ff1 fs0 fc0 sc0 ls0 ws0">%%%%%%%%%%%<span class="_ _0"></span><span class="ff2">载入信号</span></div><div class="t m0 x1 h2 y2 ff1 fs0 fc0 sc0 ls0 ws0">x=load('1.txt');%<span class="_ _1"> </span><span class="ff2">产生信号</span></div><div class="t m0 x1 h2 y3 ff1 fs0 fc0 sc0 ls0 ws0">N=length(x); <span class="_ _2"> </span>%<span class="_ _3"></span><span class="ff2">采样点数</span></div><div class="t m0 x1 h2 y4 ff1 fs0 fc0 sc0 ls0 ws0">fs=2000; <span class="_ _4"> </span>%<span class="_ _5"></span><span class="ff2">采样频率</span></div><div class="t m0 x1 h2 y5 ff1 fs0 fc0 sc0 ls0 ws0">dt=1<span class="_ _6"></span>/fs; <span class="_ _7"> </span>%<span class="_ _3"></span><span class="ff2">采样时间间隔</span></div><div class="t m0 x1 h2 y6 ff1 fs0 fc0 sc0 ls0 ws0">t=(0:N-1)*dt; <span class="_ _8"> </span>%<span class="_ _3"></span><span class="ff2">产生时间序列</span></div><div class="t m0 x1 h3 y7 ff1 fs1 fc0 sc0 ls0 ws0">%%%%%%%%EMD </div><div class="t m0 x1 h3 y8 ff1 fs1 fc0 sc0 ls0 ws0">imf=emd(x); </div><div class="t m0 x1 h4 y9 ff1 fs2 fc0 sc0 ls0 ws0">% EMD.M (EMD<span class="_ _9"> </span><span class="ff2">程序<span class="_ _a"></span></span>) </div><div class="t m0 x1 h3 ya ff1 fs1 fc0 sc0 ls0 ws0">% G. Rilling, July 2002</div><div class="t m0 x1 h3 yb ff1 fs1 fc0 sc0 ls0 ws0">%</div><div class="t m0 x1 h3 yc ff1 fs1 fc0 sc0 ls0 ws0">% computes EMD (Empirical Mode Decomposition) according to:</div><div class="t m0 x1 h3 yd ff1 fs1 fc0 sc0 ls0 ws0">%</div><div class="t m0 x1 h3 ye ff1 fs1 fc0 sc0 ls0 ws0">% N. E. Huang et al., "The empirical mode decomposition and the </div><div class="t m0 x1 h3 yf ff1 fs1 fc0 sc0 ls0 ws0">% <span class="_ _b"></span>Hilbert <span class="_ _c"> </span>spectrum <span class="_ _d"> </span>for <span class="_ _e"> </span>non-linear <span class="_ _f"> </span>and <span class="_ _10"> </span>non <span class="_ _10"> </span>stationary <span class="_ _11"> </span>time <span class="_ _12"> </span>series <span class="_ _d"> </span>analysis," </div><div class="t m0 x1 h3 y10 ff1 fs1 fc0 sc0 ls0 ws0">% Proc. Royal Soc. London A, Vol. 454, pp. 903-995, 1998</div><div class="t m0 x1 h3 y11 ff1 fs1 fc0 sc0 ls0 ws0">%</div><div class="t m0 x1 h3 y12 ff1 fs1 fc0 sc0 ls0 ws0">% with variations reported in:</div><div class="t m0 x1 h3 y13 ff1 fs1 fc0 sc0 ls0 ws0">%</div><div class="t m0 x1 h5 y14 ff1 fs1 fc0 sc0 ls0 ws0">% G. Rilling, P. Flandrin and P. Gon?alv<span class="_ _13"> </span><span class="ff2">èｓ</span></div><div class="t m0 x1 h3 y15 ff1 fs1 fc0 sc0 ls0 ws0">% "On Empirical Mode Decomposition and its algorithms"</div><div class="t m0 x1 h3 y16 ff1 fs1 fc0 sc0 ls0 ws0">% IEEE-EURASIP Workshop on Nonlinear Signal and Image Processing</div><div class="t m0 x1 h3 y17 ff1 fs1 fc0 sc0 ls0 ws0">% NSIP-03, Grado (I), June 2003</div><div class="t m0 x1 h3 y18 ff1 fs1 fc0 sc0 ls0 ws0">%</div><div class="t m0 x1 h3 y19 ff1 fs1 fc0 sc0 ls0 ws0">% stopping criterion for sifting : </div><div class="t m0 x1 h3 y1a ff1 fs1 fc0 sc0 ls0 ws0">% at each point : mean amplitude < threshold2*envelope amplitude</div><div class="t m0 x1 h3 y1b ff1 fs1 fc0 sc0 ls0 ws0">% &</div><div class="t m0 x1 h3 y1c ff1 fs1 fc0 sc0 ls0 ws0">% mean of boolean array ((mean amplitude)/(envelope amplitude) > </div><div class="t m0 x1 h3 y1d ff1 fs1 fc0 sc0 ls0 ws0">threshold) < tolerance</div><div class="t m0 x1 h3 y1e ff1 fs1 fc0 sc0 ls0 ws0">% &</div><div class="t m0 x1 h3 y1f ff1 fs1 fc0 sc0 ls0 ws0">% |#zeros-#extrema|<=1</div><div class="t m0 x1 h3 y20 ff1 fs1 fc0 sc0 ls0 ws0">%</div><div class="t m0 x1 h3 y21 ff1 fs1 fc0 sc0 ls0 ws0">% inputs: - x : analysed signal (line vector)</div><div class="t m0 x1 h3 y22 ff1 fs1 fc0 sc0 ls0 ws0">% - t (optional) : sampling times (line vector) (default : </div><div class="t m0 x1 h3 y23 ff1 fs1 fc0 sc0 ls0 ws0">1:length(x))</div><div class="t m0 x1 h3 y24 ff1 fs1 fc0 sc0 ls0 ws0">% - stop (optional) : threshold, threshold2 and tolerance </div><div class="t m0 x1 h3 y25 ff1 fs1 fc0 sc0 ls0 ws0">(optional)</div><div class="t m0 x1 h3 y26 ff1 fs1 fc0 sc0 ls0 ws0">% for sifting stopping criterion </div><div class="t m0 x1 h3 y27 ff1 fs1 fc0 sc0 ls0 ws0">% default : [0.05,0.5,0.05]</div><div class="t m0 x1 h3 y28 ff1 fs1 fc0 sc0 ls0 ws0">% - <span class="_ _14"> </span>tst <span class="_ _15"> </span>(optional) <span class="_ _16"> </span>: <span class="_ _17"> </span>if <span class="_ _15"> </span>equals <span class="_ _18"> </span>to <span class="_ _19"> </span>1 <span class="_ _1a"> </span>shows <span class="_ _12"> </span>sifting <span class="_ _1b"> </span>steps <span class="_ _1c"> </span>with <span class="_ _1d"> </span>pause</div><div class="t m0 x1 h3 y29 ff1 fs1 fc0 sc0 ls0 ws0">% if equals to 2 no pause</div><div class="t m0 x1 h3 y2a ff1 fs1 fc0 sc0 ls0 ws0">%</div><div class="t m0 x1 h3 y2b ff1 fs1 fc0 sc0 ls0 ws0">% outputs: - imf : intrinsic mode functions (last line = residual)</div></div><div class="pi" data-data='{"ctm":[1.002956,0.000000,0.000000,1.002956,0.000000,0.000000]}'></div></div>
</body>
</html>
<div id="pf2" class="pf w0 h0" data-page-no="2"><div class="pc pc2 w0 h0"><img class="bi x0 y0 w1 h1" alt="" src="https://static.pudn.com/prod/directory_preview_static/624f338f74bc5c010528612d/bg2.jpg"><div class="t m0 x1 h3 y2c ff1 fs1 fc0 sc0 ls0 ws0">% - ort : index of orthogonality</div><div class="t m0 x1 h3 y2d ff1 fs1 fc0 sc0 ls0 ws0">% - nbits : number of iterations for each mode</div><div class="t m0 x1 h3 y2e ff1 fs1 fc0 sc0 ls0 ws0">%</div><div class="t m0 x1 h3 y2f ff1 fs1 fc0 sc0 ls0 ws0">% calls: - extr (finds extrema and zero-crossings)</div><div class="t m0 x1 h3 y30 ff1 fs1 fc0 sc0 ls0 ws0">% - io : computes the index of orthogonality</div><div class="t m0 x1 h3 y7 ff1 fs1 fc0 sc0 ls0 ws0">function<span class="_ _1e"> </span> [imf,ort,nbits] = emd(x,t,stop,tst);</div><div class="t m0 x1 h3 y31 ff1 fs1 fc0 sc0 ls0 ws0">% default for stopping</div><div class="t m0 x1 h3 y32 ff1 fs1 fc0 sc0 ls0 ws0">defstop = [0.05,0.5,0.05];</div><div class="t m0 x1 h3 yb ff1 fs1 fc0 sc0 ls0 ws0">if<span class="_ _e"> </span>(nargin==1)</div><div class="t m0 x1 h3 yc ff1 fs1 fc0 sc0 ls0 ws0"> t = 1:length(x);</div><div class="t m0 x1 h3 yd ff1 fs1 fc0 sc0 ls0 ws0"> stop = defstop;</div><div class="t m0 x1 h3 ye ff1 fs1 fc0 sc0 ls0 ws0"> tst = 0;</div><div class="t m0 x1 h3 yf ff1 fs1 fc0 sc0 ls0 ws0">end</div><div class="t m0 x1 h3 y11 ff1 fs1 fc0 sc0 ls0 ws0">if<span class="_ _e"> </span>(nargin==2)</div><div class="t m0 x1 h3 y12 ff1 fs1 fc0 sc0 ls0 ws0"> stop = defstop;</div><div class="t m0 x1 h3 y13 ff1 fs1 fc0 sc0 ls0 ws0"> tst = 0;</div><div class="t m0 x1 h3 y14 ff1 fs1 fc0 sc0 ls0 ws0">end</div><div class="t m0 x1 h3 y16 ff1 fs1 fc0 sc0 ls0 ws0">if<span class="_ _e"> </span> (nargin==3)</div><div class="t m0 x1 h3 y17 ff1 fs1 fc0 sc0 ls0 ws0"> tst=0;</div><div class="t m0 x1 h3 y18 ff1 fs1 fc0 sc0 ls0 ws0">end</div><div class="t m0 x1 h3 y1a ff1 fs1 fc0 sc0 ls0 ws0">S = size(x);</div><div class="t m0 x1 h3 y1b ff1 fs1 fc0 sc0 ls0 ws0">if<span class="_ _e"> </span> ((S(1) > 1) & (S(2) > 1)) | (length(S) > 2)</div><div class="t m0 x1 h3 y1c ff1 fs1 fc0 sc0 ls0 ws0"> error(<span class="_ _1f"> </span>'x must have only one row or one column'<span class="_ _20"> </span>)</div><div class="t m0 x1 h3 y1d ff1 fs1 fc0 sc0 ls0 ws0">end</div><div class="t m0 x1 h3 y1f ff1 fs1 fc0 sc0 ls0 ws0">if<span class="_ _e"> </span> S(1) > 1</div><div class="t m0 x1 h3 y20 ff1 fs1 fc0 sc0 ls0 ws0"> x = x';</div><div class="t m0 x1 h3 y21 ff1 fs1 fc0 sc0 ls0 ws0">end</div><div class="t m0 x1 h3 y23 ff1 fs1 fc0 sc0 ls0 ws0">S = size(t);</div><div class="t m0 x1 h3 y24 ff1 fs1 fc0 sc0 ls0 ws0">if<span class="_ _e"> </span> ((S(1) > 1) & (S(2) > 1)) | (length(S) > 2)</div><div class="t m0 x1 h3 y25 ff1 fs1 fc0 sc0 ls0 ws0"> error(<span class="_ _1f"> </span>'t must have only one row or one column'<span class="_ _21"> </span>)</div><div class="t m0 x1 h3 y26 ff1 fs1 fc0 sc0 ls0 ws0">end</div><div class="t m0 x1 h3 y28 ff1 fs1 fc0 sc0 ls0 ws0">if<span class="_ _e"> </span> S(1) > 1</div><div class="t m0 x1 h3 y29 ff1 fs1 fc0 sc0 ls0 ws0"> t = t';</div><div class="t m0 x1 h3 y2a ff1 fs1 fc0 sc0 ls0 ws0">end</div></div><div class="pi" data-data='{"ctm":[1.002956,0.000000,0.000000,1.002956,0.000000,0.000000]}'></div></div>
<div id="pf3" class="pf w0 h0" data-page-no="3"><div class="pc pc3 w0 h0"><img class="bi x0 y0 w1 h1" alt="" src="https://static.pudn.com/prod/directory_preview_static/624f338f74bc5c010528612d/bg3.jpg"><div class="t m0 x1 h3 y2c ff1 fs1 fc0 sc0 ls0 ws0">if<span class="_ _e"> </span> (length(t)~=length(x))</div><div class="t m0 x1 h3 y2d ff1 fs1 fc0 sc0 ls0 ws0"> error(<span class="_ _1f"> </span>'x and t must have the same length'<span class="_ _22"> </span>)</div><div class="t m0 x1 h3 y2e ff1 fs1 fc0 sc0 ls0 ws0">end</div><div class="t m0 x1 h3 y30 ff1 fs1 fc0 sc0 ls0 ws0">S = size(stop);</div><div class="t m0 x1 h3 y33 ff1 fs1 fc0 sc0 ls0 ws0">if<span class="_ _e"> </span> ((S(1) > 1) & (S(2) > 1)) | (S(1) > 3) | (S(2) > 3) | (length(S) > 2)</div><div class="t m0 x1 h3 y7 ff1 fs1 fc0 sc0 ls0 ws0"> error(<span class="_ _1f"> </span>'stop <span class="_ _23"> </span>must <span class="_ _24"> </span>have <span class="_ _25"> </span>only <span class="_ _15"> </span>one <span class="_ _26"> </span>row <span class="_ _25"> </span>or <span class="_ _27"> </span>one <span class="_ _26"> </span>column <span class="_ _28"> </span>of <span class="_ _19"> </span>max <span class="_"> </span>three <span class="_ _29"> </span>elements'<span class="_ _c"> </span>)</div><div class="t m0 x1 h3 y34 ff1 fs1 fc0 sc0 ls0 ws0">end</div><div class="t m0 x1 h3 y32 ff1 fs1 fc0 sc0 ls0 ws0">if<span class="_ _e"> </span> S(1) > 1</div><div class="t m0 x1 h3 ya ff1 fs1 fc0 sc0 ls0 ws0"> stop = stop';</div><div class="t m0 x1 h3 yb ff1 fs1 fc0 sc0 ls0 ws0"> S = size(stop);</div><div class="t m0 x1 h3 yc ff1 fs1 fc0 sc0 ls0 ws0">end</div><div class="t m0 x1 h3 ye ff1 fs1 fc0 sc0 ls0 ws0">if<span class="_ _e"> </span> S(2) < 3</div><div class="t m0 x1 h3 yf ff1 fs1 fc0 sc0 ls0 ws0"> stop(3)=defstop(3);</div><div class="t m0 x1 h3 y10 ff1 fs1 fc0 sc0 ls0 ws0">end</div><div class="t m0 x1 h3 y12 ff1 fs1 fc0 sc0 ls0 ws0">if<span class="_ _e"> </span> S(2) < 2</div><div class="t m0 x1 h3 y13 ff1 fs1 fc0 sc0 ls0 ws0"> stop(2)=defstop(2);</div><div class="t m0 x1 h3 y14 ff1 fs1 fc0 sc0 ls0 ws0">end</div><div class="t m0 x1 h3 y16 ff1 fs1 fc0 sc0 ls0 ws0">sd = stop(1);</div><div class="t m0 x1 h3 y17 ff1 fs1 fc0 sc0 ls0 ws0">sd2 = stop(2);</div><div class="t m0 x1 h3 y18 ff1 fs1 fc0 sc0 ls0 ws0">tol = stop(3);</div><div class="t m0 x1 h3 y1a ff1 fs1 fc0 sc0 ls0 ws0">if<span class="_ _e"> </span> tst</div><div class="t m0 x1 h3 y1b ff1 fs1 fc0 sc0 ls0 ws0"> figure</div><div class="t m0 x1 h3 y1c ff1 fs1 fc0 sc0 ls0 ws0">end</div><div class="t m0 x1 h3 y1e ff1 fs1 fc0 sc0 ls0 ws0">% maximum number of iterations</div><div class="t m0 x1 h3 y1f ff1 fs1 fc0 sc0 ls0 ws0">MAXITERATIONS=2000;</div><div class="t m0 x1 h3 y21 ff1 fs1 fc0 sc0 ls0 ws0">% maximum number of symmetrized points for interpolations</div><div class="t m0 x1 h3 y22 ff1 fs1 fc0 sc0 ls0 ws0">NBSYM = 2;</div><div class="t m0 x1 h3 y24 ff1 fs1 fc0 sc0 ls0 ws0">lx = length(x);</div><div class="t m0 x1 h3 y26 ff1 fs1 fc0 sc0 ls0 ws0">sdt(lx) = 0;</div><div class="t m0 x1 h3 y27 ff1 fs1 fc0 sc0 ls0 ws0">sdt = sdt+sd;</div><div class="t m0 x1 h3 y28 ff1 fs1 fc0 sc0 ls0 ws0">sd2t(lx) = 0;</div><div class="t m0 x1 h3 y29 ff1 fs1 fc0 sc0 ls0 ws0">sd2t = sd2t+sd2;</div><div class="t m0 x1 h3 y2b ff1 fs1 fc0 sc0 ls0 ws0">% maximum number of extrema and zero-crossings in residual</div></div><div class="pi" data-data='{"ctm":[1.002956,0.000000,0.000000,1.002956,0.000000,0.000000]}'></div></div>
<div id="pf4" class="pf w0 h0" data-page-no="4"><div class="pc pc4 w0 h0"><img class="bi x0 y0 w1 h1" alt="" src="https://static.pudn.com/prod/directory_preview_static/624f338f74bc5c010528612d/bg4.jpg"><div class="t m0 x1 h3 y2c ff1 fs1 fc0 sc0 ls0 ws0">ner = lx;</div><div class="t m0 x1 h3 y2d ff1 fs1 fc0 sc0 ls0 ws0">nzr = lx;</div><div class="t m0 x1 h3 y2f ff1 fs1 fc0 sc0 ls0 ws0">r = x;</div><div class="t m0 x1 h3 y30 ff1 fs1 fc0 sc0 ls0 ws0">imf = [];</div><div class="t m0 x1 h3 y33 ff1 fs1 fc0 sc0 ls0 ws0">k = 1;</div><div class="t m0 x1 h3 y34 ff1 fs1 fc0 sc0 ls0 ws0">% iterations counter for extraction of 1 mode</div><div class="t m0 x1 h3 y31 ff1 fs1 fc0 sc0 ls0 ws0">nbit=0;</div><div class="t m0 x1 h3 ya ff1 fs1 fc0 sc0 ls0 ws0">% total iterations counter</div><div class="t m0 x1 h3 yb ff1 fs1 fc0 sc0 ls0 ws0">NbIt=0;</div><div class="t m0 x1 h3 yd ff1 fs1 fc0 sc0 ls0 ws0">while<span class="_ _2a"> </span> ner > 2</div><div class="t m0 x2 h3 yf ff1 fs1 fc0 sc0 ls0 ws0">% current mode</div><div class="t m0 x1 h3 y10 ff1 fs1 fc0 sc0 ls0 ws0"> m = r;</div><div class="t m0 x2 h3 y12 ff1 fs1 fc0 sc0 ls0 ws0">% mode at previous iteration</div><div class="t m0 x1 h3 y13 ff1 fs1 fc0 sc0 ls0 ws0"> mp = m;</div><div class="t m0 x1 h3 y15 ff1 fs1 fc0 sc0 ls0 ws0"> sx = sd+1;</div><div class="t m0 x2 h3 y17 ff1 fs1 fc0 sc0 ls0 ws0">% tests if enough extrema to proceed</div><div class="t m0 x1 h3 y18 ff1 fs1 fc0 sc0 ls0 ws0"> test = 0;</div><div class="t m0 x1 h3 y1a ff1 fs1 fc0 sc0 ls0 ws0"> [indmin,indmax,indzer] = extr(m);</div><div class="t m0 x1 h3 y1b ff1 fs1 fc0 sc0 ls0 ws0"> lm=length(indmin);</div><div class="t m0 x1 h3 y1c ff1 fs1 fc0 sc0 ls0 ws0"> lM=length(indmax);</div><div class="t m0 x1 h3 y1d ff1 fs1 fc0 sc0 ls0 ws0"> nem=lm + lM;</div><div class="t m0 x1 h3 y1e ff1 fs1 fc0 sc0 ls0 ws0"> nzm=length(indzer);</div><div class="t m0 x1 h3 y20 ff1 fs1 fc0 sc0 ls0 ws0"> j=1;</div><div class="t m0 x2 h3 y22 ff1 fs1 fc0 sc0 ls0 ws0">% sifting loop</div><div class="t m0 x2 h3 y23 ff1 fs1 fc0 sc0 ls0 ws0">while<span class="_ _2b"> </span>( <span class="_ _17"> </span>mean(sx <span class="_ _2c"> </span>> <span class="_ _10"> </span>sd) <span class="_ _12"> </span>> <span class="_ _2d"> </span>tol <span class="_ _2e"> </span>| <span class="_ _27"> </span>any(sx <span class="_ _2f"> </span>> <span class="_"> </span>sd2) <span class="_ _28"> </span>| <span class="_ _19"> </span>(abs(nzm-nem)>1)) <span class="_ _14"> </span>& <span class="_"> </span>(test </div><div class="t m0 x1 h3 y24 ff1 fs1 fc0 sc0 ls0 ws0">== 0) & nbit<MAXITERATIONS</div><div class="t m0 x3 h3 y26 ff1 fs1 fc0 sc0 ls0 ws0">if<span class="_ _30"> </span>(nbit>MAXITERATIONS/5 & mod(nbit,floor(MAXITERATIONS/10))==0)</div><div class="t m0 x1 h3 y27 ff1 fs1 fc0 sc0 ls0 ws0"> disp([<span class="_ _31"> </span>'mode <span class="_ _32"> </span>'<span class="_ _33"> </span>,int2str(k),<span class="_ _34"> </span>' <span class="_ _27"> </span>nombre <span class="_ _32"> </span>d <span class="_ _10"> </span>iterations <span class="_ _35"> </span>: <span class="_ _17"> </span>'<span class="_ _33"> </span>,int2str(nbit)])</div><div class="t m0 x1 h3 y28 ff1 fs1 fc0 sc0 ls0 ws0"> disp([<span class="_ _31"> </span>'stop parameter mean value : '<span class="_ _36"> </span>,num2str(s)])</div><div class="t m0 x3 h3 y29 ff1 fs1 fc0 sc0 ls0 ws0">end</div><div class="t m0 x3 h3 y2b ff1 fs1 fc0 sc0 ls0 ws0">% boundary conditions for interpolations :</div></div><div class="pi" data-data='{"ctm":[1.002956,0.000000,0.000000,1.002956,0.000000,0.000000]}'></div></div>
<div id="pf5" class="pf w0 h0" data-page-no="5"><div class="pc pc5 w0 h0"><img class="bi x0 y0 w1 h1" alt="" src="https://static.pudn.com/prod/directory_preview_static/624f338f74bc5c010528612d/bg5.jpg"><div class="t m0 x4 h3 y2d ff1 fs1 fc0 sc0 ls0 ws0">if<span class="_ _30"> </span> indmax(1) < indmin(1)</div><div class="t m0 x5 h3 y2e ff1 fs1 fc0 sc0 ls0 ws0">if<span class="_ _30"> </span> m(1) > m(indmin(1))</div><div class="t m0 x1 h3 y2f ff1 fs1 fc0 sc0 ls0 ws0"> lmax = fliplr(indmax(2:min(end,NBSYM+1)));</div><div class="t m0 x1 h3 y30 ff1 fs1 fc0 sc0 ls0 ws0"> lmin = fliplr(indmin(1:min(end,NBSYM)));</div><div class="t m0 x1 h3 y33 ff1 fs1 fc0 sc0 ls0 ws0"> lsym = indmax(1);</div><div class="t m0 x5 h3 y7 ff1 fs1 fc0 sc0 ls0 ws0">else</div><div class="t m0 x1 h3 y34 ff1 fs1 fc0 sc0 ls0 ws0"> lmax = fliplr(indmax(1:min(end,NBSYM)));</div><div class="t m0 x1 h3 y31 ff1 fs1 fc0 sc0 ls0 ws0"> lmin = [fliplr(indmin(1:min(end,NBSYM-1))),1];</div><div class="t m0 x1 h3 y32 ff1 fs1 fc0 sc0 ls0 ws0"> lsym = 1;</div><div class="t m0 x5 h3 ya ff1 fs1 fc0 sc0 ls0 ws0">end</div><div class="t m0 x3 h3 yb ff1 fs1 fc0 sc0 ls0 ws0">else</div><div class="t m0 x5 h3 yd ff1 fs1 fc0 sc0 ls0 ws0">if<span class="_ _30"> </span> m(1) < m(indmax(1))</div><div class="t m0 x1 h3 ye ff1 fs1 fc0 sc0 ls0 ws0"> lmax = fliplr(indmax(1:min(end,NBSYM)));</div><div class="t m0 x1 h3 yf ff1 fs1 fc0 sc0 ls0 ws0"> lmin = fliplr(indmin(2:min(end,NBSYM+1)));</div><div class="t m0 x1 h3 y10 ff1 fs1 fc0 sc0 ls0 ws0"> lsym = indmin(1);</div><div class="t m0 x5 h3 y11 ff1 fs1 fc0 sc0 ls0 ws0">else</div><div class="t m0 x1 h3 y12 ff1 fs1 fc0 sc0 ls0 ws0"> lmax = [fliplr(indmax(1:min(end,NBSYM-1))),1];</div><div class="t m0 x1 h3 y13 ff1 fs1 fc0 sc0 ls0 ws0"> lmin = fliplr(indmin(1:min(end,NBSYM)));</div><div class="t m0 x1 h3 y14 ff1 fs1 fc0 sc0 ls0 ws0"> lsym = 1;</div><div class="t m0 x5 h3 y15 ff1 fs1 fc0 sc0 ls0 ws0">end</div><div class="t m0 x3 h3 y16 ff1 fs1 fc0 sc0 ls0 ws0">end</div><div class="t m0 x3 h3 y18 ff1 fs1 fc0 sc0 ls0 ws0">if<span class="_ _30"> </span> indmax(end) < indmin(end)</div><div class="t m0 x5 h3 y19 ff1 fs1 fc0 sc0 ls0 ws0">if<span class="_ _30"> </span> m(end) < m(indmax(end))</div><div class="t m0 x1 h3 y1a ff1 fs1 fc0 sc0 ls0 ws0"> rmax = fliplr(indmax(max(end-NBSYM+1,1):end));</div><div class="t m0 x1 h3 y1b ff1 fs1 fc0 sc0 ls0 ws0"> rmin = fliplr(indmin(max(end-NBSYM,1):end-1));</div><div class="t m0 x1 h3 y1c ff1 fs1 fc0 sc0 ls0 ws0"> rsym = indmin(end);</div><div class="t m0 x5 h3 y1d ff1 fs1 fc0 sc0 ls0 ws0">else</div><div class="t m0 x1 h3 y1e ff1 fs1 fc0 sc0 ls0 ws0"> rmax = [lx,fliplr(indmax(max(end-NBSYM+2,1):end))];</div><div class="t m0 x1 h3 y1f ff1 fs1 fc0 sc0 ls0 ws0"> rmin = fliplr(indmin(max(end-NBSYM+1,1):end));</div><div class="t m0 x1 h3 y20 ff1 fs1 fc0 sc0 ls0 ws0"> rsym = lx;</div><div class="t m0 x5 h3 y21 ff1 fs1 fc0 sc0 ls0 ws0">end</div><div class="t m0 x3 h3 y22 ff1 fs1 fc0 sc0 ls0 ws0">else</div><div class="t m0 x5 h3 y23 ff1 fs1 fc0 sc0 ls0 ws0">if<span class="_ _30"> </span> m(end) > m(indmin(end))</div><div class="t m0 x1 h3 y24 ff1 fs1 fc0 sc0 ls0 ws0"> rmax = fliplr(indmax(max(end-NBSYM,1):end-1));</div><div class="t m0 x1 h3 y25 ff1 fs1 fc0 sc0 ls0 ws0"> rmin = fliplr(indmin(max(end-NBSYM+1,1):end));</div><div class="t m0 x1 h3 y26 ff1 fs1 fc0 sc0 ls0 ws0"> rsym = indmax(end);</div><div class="t m0 x5 h3 y27 ff1 fs1 fc0 sc0 ls0 ws0">else</div><div class="t m0 x1 h3 y28 ff1 fs1 fc0 sc0 ls0 ws0"> rmax = fliplr(indmax(max(end-NBSYM+1,1):end));</div><div class="t m0 x1 h3 y29 ff1 fs1 fc0 sc0 ls0 ws0"> rmin = [lx,fliplr(indmin(max(end-NBSYM+2,1):end))];</div><div class="t m0 x1 h3 y2a ff1 fs1 fc0 sc0 ls0 ws0"> rsym = lx;</div><div class="t m0 x5 h3 y2b ff1 fs1 fc0 sc0 ls0 ws0">end</div></div><div class="pi" data-data='{"ctm":[1.002956,0.000000,0.000000,1.002956,0.000000,0.000000]}'></div></div>
