neighbor.m

基于Matlab

function s2 = neighbor(s1,M,real_source)

if real_source
 C=2*[0:M-1]-M+1;  
 C=C/sqrt(mean(abs(C).^2)); 
else
 C=2*[0:sqrt(M)-1]-sqrt(M)+1;  
 C=ones(sqrt(M),1)*C;
 C=C+i*C';
 C=C(:);
 C=C/sqrt(mean(abs(C).^2)); 
end;

C=C(find(abs(s1-C)>100*eps));
d=round(100*abs(s1-C));
I=find(abs(d-min(d))==0);
[y,J]=min(d(I)+randn(size(d(I))));
s2=C(I(J));