Z?=?peaks（512）;?
Z=Z*10;
SNR=0;
B=awgn（ZSNR）;%?pluse?gaussian?white?noise
U=cos（B）+1*i*sin（B）;
[widthheight]=size（U）;%錕斤拷取圖錕斤拷錕餃和高訛拷
UN?=?U;
real1=zeros（widthheight）;
imag1=zeros（widthheight）;

for?x=1:width
????for?y=1:height
???????
?????????????real1（xy）=?real（U（xy））;
?????????????imag1（xy）?=?imag（U（xy））;
????????????%?R=?sqrt（real1（xy）^2+imag1（xy）^2）;
?????????????%UN（xy）=U（xy）/R;%%錕斤拷一錕斤拷
????end
end
???????
piex=zeros（55）;
real2=zeros（55）;
imag2=zeros（55）;
r1?=?zeros（widthheight）;
for?x=3:width-2
????for?y=3:height-2
????????piex=UN（x-2:x+2y-2:y+2）;
????????real2（::）=real（piex）;
????????imag2（::）=imag（piex）;
????????real3=sum（sum（real2））/25;
????????imag3=sum（sum（imag2））/25;
????????r1=sqrt（real3^2+imag3^2）;%%錕斤拷模
????????r=r1;
????????%r3=abs（r3）;
????end
????
end

%?錕借定錕斤拷錕斤拷錕絞鹼拷錕斤拷錕斤拷錕?%?FFT錕斤拷錕津窗匡拷錕借定

m=15;
n=15;

armw=（m-1）/2;%錕斤拷錕節匡拷錕?
armh=（n-1）/2;%錕斤拷錕節高訛拷


%?錕借定錕斤拷通錕剿詫拷錕斤拷錕斤拷
w=zeros（mn）;
w（armw:armw+2armh:armh+2）=1;
%?錕斤拷始錕斤拷錕斤拷錕斤拷
result=ones（width-m+1height-n+1）;
i1=ones（mn）;
i2=i1;
s=i1;
windows=i1;
%?錕借定錕斤拷權系錕斤拷
%r=0.5;

%-------------------------------------------------------------------------
%