cvmn.m

来自「C 程序 球面方位角算法： azimuth.cpp CRC校验快速算法： C」· M 代码 · 共 20 行

function Mm=cvmn(M,wh,ty)
%均值滤波：wh—窗宽;    ty=1  取邻域均值；	ty=0: 去除邻域均值（去除趋势线）
%格式: Mm=cvmn(M,wh,ty)
size=length(M);
L=floor(wh/2); R=L; A=L+1+R;
for(i=1:L)        Mm(i)=mean(M(1:i+R));  end
Mm(L+1)=sum(M(1:A));
for(i=L+2:size-R-1) 
    Mm(i)=Mm(i-1)-M(i-L-1)+M(i+R); 
end
Mm(L+1:size-R-1)=Mm(L+1:size-R-1)/(A);
for(i=size-R:size)  Mm(i)=mean(M(i-L:size));  end  

if(~ty) 
    [m,n]=size(M);
    if(m>1) M=M'; end
    Mm=M-Mm; 
end
return;