https://wiki-pnb.eri.ucsb.edu/index.php?action=history&feed=atom&title=Moavg.m 2026-04-08T03:03:48Z Revision history for this page on the wiki MediaWiki 1.34.0 https://wiki-pnb.eri.ucsb.edu/index.php?title=Moavg.m&diff=209&oldid=prev 2010-08-05T23:46:51Z

<p>Created page with &#039;&lt;pre&gt; function mavg = moavg(flag,diff,r,c,matrix,trend,pts) %created by Erik Stassinos %3/18/2010 %function called by dspike to calculate moving average of points %passes in row …&#039;</p> <p><b>New page</b></p><div>&lt;pre&gt;<br /> function mavg = moavg(flag,diff,r,c,matrix,trend,pts)<br /> %created by Erik Stassinos<br /> %3/18/2010<br /> %function called by dspike to calculate moving average of points<br /> %passes in row and column position to do average around<br /> <br /> cnt = 0;<br /> sum = 0;<br /> <br /> width_2 = (flag)/2;<br /> width_2 = floor(width_2);<br /> %width_2 = 2<br /> %width_2 = floor(width_2);<br /> start = r;<br /> %r number of rows (row location)<br /> %c number of columns (column location)<br /> %matrix is up or downcast matrix<br /> %diff is array of differences calculated<br /> %----------------looping to find sum,cnt and trapflag-------------<br /> r;<br /> c;<br /> <br /> %--------moving avg for start of matrix<br /> <br /> if start &lt;= width_2<br /> for i = start:(start+width_2) % rows in column<br /> if ((i==start) || (diff(i,c) &gt; trend) || (matrix(i,c) &lt; -1000));<br /> else<br /> sum = sum + matrix(i,c);<br /> cnt = cnt +1;<br /> end<br /> if cnt ==0<br /> cnt;<br /> mavg = -9.9E35;<br /> %sprintf('inconclusive mavg calc, trapflag set')<br /> %break<br /> else<br /> sum;<br /> mavg = sum/cnt;<br /> end<br /> end<br /> end<br /> clear i<br /> %-------moving avg for end of matrix<br /> if start &gt;= (pts - width_2)<br /> for i = (start-width_2):start <br /> if ((i==start) || (diff(i,c) &gt; trend) || (matrix(i,c) &lt; -1000));<br /> else<br /> sum = sum + matrix(i,c);<br /> cnt = cnt+1;<br /> end<br /> if cnt ==0<br /> cnt;<br /> mavg = -9.9E35;<br /> %sprintf('inconclusive mavg calc, trapflag set')<br /> <br /> else<br /> mavg = sum/cnt;<br /> end<br /> end <br /> end<br /> clear i<br /> <br /> %-------moving avg for everything in between<br /> if start &gt; width_2 &amp;&amp; start &lt;= (pts - width_2)<br /> for i= (start - width_2):(start+width_2)<br /> if ((i==start) || (diff(i,c) &gt; trend) || (matrix(i,c) &lt; -1000));<br /> else<br /> sum = sum + matrix(i,c);<br /> matrixic = matrix(i,c);<br /> cnt = cnt+1;<br /> end<br /> if cnt ==0<br /> cnt;<br /> mavg = -9.9E35;<br /> %sprintf('inconclusive mavg calc, trapflag set')<br /> else<br /> mavg = sum/cnt;<br /> end<br /> end<br /> end<br /> clear i<br /> %----------------------end looping--------------------------<br /> <br /> &lt;/pre&gt;<br /> <br /> &lt;big&gt;[[moavg]]&lt;big&gt;</div>

128.111.101.185