matlab SAR图像滤波

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;

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