matlab sar 图像检测的程序

%SAR.m
%产生SAR的同波
clear all
thetaT=0; %波束宽度角
thetaT=thetaT*pi/180;%转换成弧度
thetaR=0; %波束高度角
thetaR=thetaR*pi/180;
c=3e8; %光速
fc=10e9; %载频
lambda=c/fc; %波长
%测绘带区域
X0=200;%方位向【-X0，X0】
Rtc=3000;
Rrc=3000;
Rc=（Rtc+Rrc）/2;
R0=150; %距离向[Rc-R0，Rc+R0】
%距离向（Range），快时间域
Tr=1.5e-6;  %LFM信号脉宽1.5us（200m）
Br=1.5e8; %LFM信号带宽150MHz
Kr=Br/Tr;   %调频斜率
Nr=512;     %快时间采样点数
r=Rc+linspace（-R0R0Nr）; %瞬时斜距
t=2*r/c; %时域序列
dt=R0*4/c/Nr; %采样周期
f=linspace（-1/2/dt1/2/dtNr）;  %频域序列
%方位向（Azimuth，Cross—Range），慢时间域
v=100; %SAR平台速度
D=1;%天线方位尺寸
Lsar=lambda*Rc/D; %合成孔径长度
Na=1024; %快时间采样点数
x=linspace（-X0X0Na）; %慢时域序列
u=x/v; %慢时间
du=2*X0/v/Na;%采样周期
fu=linspace（-1/2/du1/2/duNa）; %慢频域序列
ftdc=v*sin（thetaT）;
ftdr=-（v*cos（thetaT））^2/lambda/Rtc;
frdc=v*sin（thetaR）;
frdr=-（v*cos（thetaR））^2/lambda/Rrc;
fdc=ftdc+frdc; %距离走动增量
fdr=ftdr+frdr; %多普勒调频率
%目标位置
Ntar=3; %目标个数
Ptar=[Rrc01
      Rrc+50-501
      Rrc+50501]; %距离向坐标，方位向坐标，目标散射系数（sigma）
%产生同波
s_ut=zeros（NrNa）;
U=ones（Nr1）*u; %扩充为矩阵
T=t‘*ones（1Na）;
for i=1:1:Ntar
      rn=Ptar（i1）;xn=Ptar（i2）;sigma=Ptar（i3）;%获取点目标散射系数
      rtn=rn+Rtc-Rrc;
      RT=sqrt（rtn^2+（rtn*tan（thetaT）+xn-v*U）.^2）;
      RR=sqrt（rn^2+（rn*tan（thetaT）+xn-v*U）.^2）;
      R=RT+RR;
      DT=T-R/c; %快时间延迟
      phase=pi*Kr*DT.^2-2*pi/lambda*R;%同波相位
      s_ut=s_ut+sigma*exp（1i*phase）.*（abs（DT） end;
%距离压缩
p0_t=exp（1i*pi*Kr*（t-2*Rc/c）.^2）.*（abs（t-2*Rc/c）p0_f=fftshift（fft（fftshift（p0_t）））;
s_uf=fftshift（fft（fftshift（s_ut）））; %距离向FFT
src_uf=s_u