MATLAB的滑动空间距平

kongfeng0824

function StaAno = StaFiltAno( location, stavalue, winsize )
    StaAno = zeros(length(stavalue),1);
    %% the distance of each station
    s = length(location);
    DisMat = zeros(s,s);
    for m = 1:s
        for n = m:s
            tempdis = sqrt(( location(1,m) - location(1,n) )^2 + ( location(2,m) - location(2,n) )^2 );
            DisMat(m,n) = tempdis;
            DisMat(n,m) = tempdis;
            clear temp*;
        end
    end
   
    %%
    for m = 1:s
        tempdis = DisMat(m,:);
        tempindex = find(tempdis < winsize);
        tempvalue = stavalue(tempindex);
        StaAno(m) = stavalue(m) - mean(tempvalue);
    end
end