补充：互质阵列平滑算法的子函数

抱歉，刚才遗漏了，这是代码里ssp的函数定义

function crs = ssp(cr, K)
% spatial smoothing 
%Rzz=cr*cr'/n;
[M,MM]=size(cr);
N=M-K+1;
crs = zeros(K,K);
for  in =1:N
        Rzz=cr(in:in+K-1,:)*cr(in:in+K-1,:)';
        crs = crs + Rzz;
end
crs = crs / N;