合成孔径雷达点目标仿真MATLAB程序

合成孔径雷达成像系统点目标仿真源程序：clcclose allC=3e8； %光速Fc=1e9； ％载波频率lambda=C/Fc； ％波长％成像区域Xmin=0；Xmax=50;Yc=10000;Y0=500；％SAR 基本参数V=100; ％雷达平台速度H=0； ％雷达平台高度R0=sqrt（Yc^2+H^2)；D=4； ％天线孔径长度Lsar=lambda*R0/D; ％合成孔径长度Tsar=Lsar/V; %合成孔径时间Ka=-2＊V^2/lambda/R0；%线性调频率Ba=abs(Ka*Tsar)；PRF=2＊Ba； %脉冲重复频率PRT=1/PRF；ds=PRT； ％脉冲重复周期 Nslow=ceil((Xmax—Xmin+Lsar）/V/ds);％脉冲数Nslow=2^nextpow2(Nslow）； ％量化为 2 的指数sn=linspace((Xmin—Lsar/2)/V，(Xmax+Lsar/2）/V，Nslow); ％创建时间向量PRT=（Xmax-Xmin+Lsar）/V/Nslow; ％更新 PRF=1/PRT； ％ 更新脉冲重复频率fa=linspace（—0.5＊PRF，0。5＊PRF,Nslow);Tr=5e-6； ％脉冲宽度Br=30e6; %调频信号带宽Kr=Br/Tr; %调频率Fsr=2*Br； %快时间域取样频率dt=1/Fsr； ％快时间域取样间隔Rmin=sqrt（(Yc—Y0）^2+H^2）；Rmax=sqrt(（Yc+Y0）^2+H^2+(Lsar/2)^2）；Nfast=ceil(2＊（Rmax—Rmin）/C/dt+Tr/dt）;Nfast=2^nextpow2(Nfast）; tm=linspace(2*Rmin/C,2＊Rmax/C+Tr,Nfast）;dt=（2*Rmax/C+Tr—2＊Rmin/C）/Nfast; ％更新Fsr=1/dt；fr=linspace(-0.5＊Fsr，0.5*Fsr，Nfast);DY=C/2/Br； %距离分辨率DX=D/2; ％方位分辨率Ntarget=3; ％目标数目Ptarget=[Xmin，Yc,1 %目标位置 Xmin,Yc+10＊DY,1 Xmin+20＊DX，Yc+50＊DY,1］; K=Ntarget； %目标数目N=Nslow； ％慢时间采样数M=Nfast; %快时间采样数T=Ptarget; %目标位置%合成孔径回波仿真Srnm=zeros（N，M）；for k=1:1:K sigma=T(k，3）; Dslow=sn＊V-T（k,1）； R=sqrt(Dslow.^2+T（k，2）^2+H^2）; tau=2＊R/C； Dfast=ones（N，1)＊tm—tau'＊ones(1,M); phase=pi＊Kr＊Dfast。^2—（4＊pi/lambda）*（R’*ones(1,M）)； Srnm=Srnm+sigma＊exp(j＊phase).＊（0〈Dfast&Dfast〈Tr).＊(（abs(Dslow）