(完整word版)SAR合成孔径雷达图像点目标仿真报告(附matlab代码)

1、SAR图像点目标仿真报告徐一凡1SAR原理简介合成孔径雷达(SyntheticApertureRadar,简称SAR)是一种高分辨率成像雷达技术。它利用脉冲压缩技术获得高的距离向分辨率，利用合成孔径原理获得高的方位向分辨率，从而获得大面积高分辨率雷达图像。SAR回波信号经距离向脉冲压缩后，雷达的距离分辨率由雷达发射信号带宽决定:0=£，式中Q,表示雷达的距离分辨率，d表示雷达发射信号带宽，C表示光速。同2巧.样，SAR回波信号经方位向合成孔径后，雷达的方位分辨率由雷达方位向的多谱勒带宽决定：Pa=h,式中0“表示雷达的方位分辨率，3“表示雷达方位向多谱勒带宽，匕表示方位向SAR平台速

2、度。在小斜视角的情况下，方位分辨率近似表示为几=,其中£为方位向合成孔径的长度。2SAR的几何关系雷达位置和波束在地面覆盖区域的简单几何模型如图1所示。此次仿真考虑的是正侧视的条带式仿真，也就是说倾斜角为零，SAR波束中心和SAR平台运动方向垂直的情况。图1雷达数据获取的几何关系建立坐标系XYZ如图2所示，其中XOY平面为地平面；SAR平台距地平面高H,以速度V沿X轴正向匀速飞行；P点为SAR平台的位置矢量，设其坐标为(x,y,z)；T点为目标的位置矢量，设其坐标为(冯,)卜0);由几何关系，目标与SAR平台的斜距为：LUUR=PT=+(y-y7.)"+(z-zr)&quo

3、t;由图可知Hy=O,z=H、S=O;令x=v5,其中v为平台速度，s为慢时间变量(slowtmie),假设xT=vs,其中s表示SAR平台的x坐标为与的时刻；再令r=yH2+yT2,.LUUPT厂表示目标与SAR的垂直斜距，重写(1)式为:=R(s;r)=JF+S(s-So)Z?(5；r)就表示任意时刻$时，目标与雷达的斜距。一般情况卞，v|5-50|«r,于是通过傅里叶技术展开，可将式可近似写为：2门=閉G)可见，斜距是S和厂的函数，不同的目标，厂也不一样，但当目标距SAR较远时，在观测带内，可近似认为厂不变，即r=RO.图2：空间几何关系(a)正视图(b)侧视图图2(a)中，心

4、如表示合成孔径长度，它和合成孔径时间兀0厂的关系是=厂。(b)中，&为雷达天线半功率点波束角，0为波束轴线与Z轴的夹角，即波束视角，/?niui为近距点距离，Rmax为远距点距离，w为测绘带宽度，它们的关系为：Rnun=H-tg(O-s%)Rmax=H/g(0+)(4)W=7?max-7?mm3SAR的回波信号模型SAR在运动中以一定的周期(1/PRF)发射和接收信号，具体过程如图3所示。发射机以耳的时间发射蜩啾脉冲，然后切换天线开关接收回波信号。I脉神更复何隔为1/PRF咂«W-I发射h接收)*图3雷达发射脉冲串的时序当雷达不处于发射状态时，它接收3反射回波。发射和接收回波

5、的时间序列如图4所示。在机载情况下，每个回波可以在脉冲发射间隔内直接接收到。但是在星载情况下，由于距离过人，某个脉冲的回波要经过610个脉冲间隔才能接收到。这里仿真为了方便，默认为机载情况。脉沖凹彼时间图4脉冲雷达的发射与接收周期假设7；为chiip信号持续时间，下标厂表示距离向：PRF为重复频率，PRT为重复周期，等于"貯。接收序列中，航警表示发射第,个脉冲时，目标回波相对于发射序列的延时。雷达的发射序列数学表达式为式(5):X$(/)=工p(t-nPRT)(5)H=-0Cpit)=rect()ep7Krr严九式中，rect(g)表示矩形信号，K厂为距离向的chnp信号调频率，&#

6、163;为载频。雷达回波信号由发射信号波形，天线方向图，斜距，目标RCS,环境等因素共同决定，若不考虑坏境因素，则单点目标雷达回波信号可写成式(6)所示：»(0=E6叩(f-PRT-tJ其中，b表示点目标的雷达散射截面，w表示点目标天线方向图双向幅度加权，4表示载机发射第n个脉冲时，电磁波再次回到载机时的延时rn=$*R；°,带入式(6)中得：z't-nPKT-2Rs,r)IC»=工bwrect()H=YC丄rexp"K(/-;?PRT-2R(s;r)/C)2-4/rexpJR(s;r)expJ2X(/-/?PRTrjA式就是单点目标回波信号模型

7、，其中，亡xp归K(fnPRT2R(s;r)/C)订是chirp分量，它决定距离向分辨率；expJZ?(5；r)为多普勒分量，它决定方位向分辨A率。对于任意一个脉冲，回波信号可表示为式(8)所示：S(人S)=wr(r一2R(s;r)/C)wa($-咚)xcxp-)4疋/(s;r)/C-(&)expj7tKr(r-2/?(j;r)/C)2我们知道，由于Z?(5；r)随慢时间s的变化而变化，所以计算机记录到的回波数据存储形式如图5所示:转目趣理現财电裁页I(裔达沐冲长度,-_丄(_,05W152025»»«!>图5目标照射时间内，单个点目标回波能量在信号

8、处理器的二维存储器中的轨迹4距离徙动及校正根据图2可知，在倾斜角为零或很小的时候，目标与雷达的瞬时距离为R(s/),根据几何关系可知，/?(5;r)=y/r2+v2(s-s0)29根据泰勒级数展开可得:RZ)=yp由式(9)可知，不同慢时间对应着不同的/?($/),并且是一个双曲线形式或者近似为一个二次形式。如图5所示，同一目标的回波存储在计算机里不在同一直线上，存在距离徙动。从而定义距离徙动量：V“A/?(5,r)=(5-50)-(10)2r为了进行方位向的压缩，方位向的回波数据必须在同一条直线上，也就是说必须校正距离徙动W,r)o由式(10)可知，不同的最近距离r对应着不同的因此在时域处理

9、距离徙动会非常麻烦。因此，对方位向进行傅里叶变换，对距离向不进行变换，得到新的域。由于方位向的频率即为多普勒频率，所以这个新的域也称为距离多普勒域。将斜距R写成多普勒fa的函数，即尺(£,厂)。众所周知，对最近距离为r的点目标2VP，回波多普勒九是倾斜角&的函数，即/,=sm6>,斜距/?(»)=广/cos。，于是兄R(fa,r)=r/cos0=r/>/1-sm20”卜炼)'(U)r+-()2rfa28V12所以距离多普勒域中的我距离徙动为A7?(£9r)=A(-)2,y；i2,可发现它不随慢时间变换，8V同一最短距离r对应着相同人小的

10、距离徙动。因此在距离多普勒域对一个距离徙动校正就是对一组具有相同最短距离的点目标的距离徙动校正，这样可以节省运算量。为了对距离徙动进行校正，需要得到距离徙动单元，即距离徙动体现在存储单元中的移动数值，距离徙动单元可以表示为尺(£")/心，这个值通常为一个分数，由于存储单元都是离散的，所以不同通过在存储单元简单的移动得到准确的值。为了得到准确的徙动校正值，通常需要进行插值运算。本仿真采用了两种插值方法最近邻点插值和sine插值，卜面分别进行介绍。最近邻点插值法的优点是简单而快速，缺点是不够精确。R(£")/p=N+，其中N为整数部分，n为小数部分，整数部分

11、徙动可以直接通过平移消除，对于小数部分则通过四舍五入的方法变为0或者1,这样就可以得到较为精确的插值。Sine插值原理如下：在基带信号下，卷积核是sine函数A(x)=smcW=(12)插值信号为g=工gd(0sillC(x-i)(13)i即为所有输入样本的加权平均。可通过频域来理解，如图6所示，采样信号gd(i)的频谱Gd(/)等于以采样率重复的信号频谱。为了重建信号g(x),只需要一个周期频谱(如基带周期)，因此需要理想矩形低通滤波器在频域中提取基带频谱(如图6)所示。已知该理想滤波器在时域中是SH1C函数。由于频域相乘相当于时域卷积，故插值可以通过与sine核的卷枳来实现。图6理想低通滤

12、波器怎样对采样信号进行插值5点目标成像matIab仿真5.1距离多普勒算法距离多普勒算法(RDA)是在1976年至1978年为民用星载SAR提出的，它兼顾了成熟、简单、高效和精确等因素，至今仍是使用最广泛的成像算法。它通过距离和方位上的频域操作，到达了高效的模块化处理要求，同时又具有了一维操作的简便性。图7示意了RDA的处理流程。这里主要讨论小倾斜角及短孔径下的基本RDA处理框图。1当数据处在方位时域时，可通过快速卷枳进行距离压缩。也就是说，距离FFT后随即进行距离向匹配滤波，再利用距离IFFT完成距离压缩。回波信号为：%(人耳二A)-2尺($)/c(S-)(14)xcxp-J4;r厶R(s)

13、/cexpy>ATr(r-2R(s)/cy距离向压缩后的信号为:(15)Snt,s)=lFFTlSQUs)HUl)=pr!-2R(s)/cwa(5-5.)exp)4好0尺($)/cH(f)=牝cf£exp-归£exp-山龙几(16)2通过方位FFT将数据变换至距离多普勒域，多普勒中心频率估计以及人部分后续操作都在该域进行。方位向傅里叶变换后信号为：s«)="Wfe)=人卩-兰也巴(£-£)(17)cxexp-jexp7>CKa3在距离多普勒域进行随距离时间及方位频率变化的RCMC,该域中同一距离上的一组目标轨迹相互重合。R

14、CMC将距离徙动曲线拉直到与方位频率轴平行的方向。这里可以采用最近邻点插值法或者sine插值法，具体插值方法见前面。假设RCMC插值是精确的，信号变为：S«)=A)以一牛)w“T,)4fR厂(xeXp-jexpcK“4通过每一距离门上的频域匹配滤波实现方位压缩。为进行方位压缩，将RCMC后的S(,£)乘以频域匹配滤波器Ha：(fs)。亿；(/)=旳一丿力织(19)KaSg)=Sg)HM)4兀fR(20)=九几U-2凡/c比(丄-几)exp-丿c5最后通过方位IFFT将数据变换回时域，得到压缩后的复图像。复原后的图像为：sa(.(t,S)=IFFTsStJs)=A)几(12&

15、amp;/0)化($)xexpJ匹丛expJ2m,$宙达瓯始数据图8距离多普勒算法流程图5.2ChirpSeaIing算法距离多普勒算法具有诸多优点，但是距离多普勒算法有两点不足：首先，当用较长的核函数提高距离徙动校正（RCMC）精度时，运算量较人；其次，二次距离压缩（SRC）对方位频率的依赖性问题较难解决，从而限制了其对某些大斜视角和长孔径SAR的处理精度。ChirpScaling算法避免了RCMC中的插值操作，通过对Chirp信号进行频率调制，实现了对该信号的尺度变换或平移。图8显示了ChnpScalmg算法处理流程。这里主要讨论小倾斜角及短孔径下的基本CSA处理框图。主要步骤包括四次FF

16、T和三次相位相乘。1通过方位向FFT将数据变换到距离多普勒域。2通过相位相乘实现ChnpScaling操作，使所有目标的距离徙动轨迹一致化。这是第一步相位相乘。用以改变线调频率尺度的ClmpScaling二次相位函数为：0（f,fa;R）=expU©（£Rjag-2"仁;RJy（22）c3通过距离向FFT将数据变到二维频域。4通过与参考函数进行相位相乘，同时完成距离压缩、SRC和一致RCMC。这是第二步相位相乘。用于距离压缩，距离徙动校正的相位函数写为：=eXp）龙f；1-r°5心九）1+心）xexpU坯血c5通过距离向IFFT将数据变回到距离多普勒域。

17、6通过与随距离变化的匹配滤波器进行相位相乘，实现方位压缩。此外，由于步骤2中的ChirpScalmg操作，相位相乘中还需要附加一项相位校正。这是第三步相位相乘。补偿由ChirpScaling引起的剩余相位函数是：比4，£；心）=exp耳Rrj£,J-£2expJ0A（£九）<24）7.最后通过方位向IFFT将数据变回到二维时域，即SAR图像域。第一步相位相乘第匸步（最厉一步）相位相来笛达原始数撫斥熔敎抿SARfdV%SAR图似域趾离多許勒域频域距离多普勒域图SChirpScaling算法流程图简而言之，R-D算法是将徙动曲线逐一校正，CS算法是以某

18、一徙动曲线为参考，在Doppler域内消除不同距离门的徙动曲线的差异，令这些曲线成为一组相互“平行”的曲线,然后在二维频率域内统一的去掉距离徙动。通俗一点就是，RD算法是将弯曲的信号一根根掰直，而CS算法是先把所有信号都掰得一样弯，然后再统一掰直。6仿真结果6.1使用最近邻点插值的距离多普勒算法仿真结果本文首先对5个点目标的回波信号进行了仿真，5个点目标构成了矩形的4个顶点和中心，其坐标分别如下，格式为（方位向,距离向，后向反射系数）：09750110097501501000010102501100102501图9的上图是距离向压缩后的图像，从图中可以看到5条回波信号（其中有几条部分重合，但仍

19、能看出来）目标回波信号存在明显的距离徙动，需要进行校正。图9的下图是通过最近邻点插值法校正后的图像，可以看出图像基本被校正为直线。图9距离向压缩后最近邻点插值的结果图10为进行方位向压缩后形成的图像，可以明显看出5个点目标，并且5个点目标构成了矩形的四个顶点及其中心。-KO300200100K7TTTY:-5.M7Mugau<cutD££CX;WTTY:EMu3WOQau»a&l>039XWMY2kRGftQVXQKiaKlXltC3»4Y:-aLMahin:C.XWO*M；fcaMT.QM7.au?X)B3»«O

20、4Y:VT衣hdn：-T.yJeXMcftQsr?,aM?.aii?a97Q9Ba991101LO2LOJ104105图10通过最近邻点插值生成的点目标图像6.2使用最近邻点插值的距离多普勒算法仿真结果图11上图为通过距离压缩后的图像，图11的卞图为通过smc插值法校正后的图像。400-200§04Dz200400距醫向压缩未校正距器徙动的图像0.950.960.970.980.991L011.021.031.041.05距离向X104距蛊向压缩,校正距醫徙动后的圉像400-200§04D小2004000.950.960.970.980.991L011.021.031.04

21、1.05距离向X104图11距离向压缩后sine插值的结果图12为进行方位向压缩后形成的图像，可以明显看出5个点目标，并且5个点目标构成了矩形的四个顶点及其中心。jcvtuy：-xvnwax402L44M0a4Q4Ci0*X1»4YgCUI33*1onX：£OT)mO4Y:RateX1£OT)mO4Y:VT34RGftMUQ411Q41Ja95Q%a97Q9Ba991L01104105图12通过sine插值生成的点目标图像6.3ChirpSealing算法仿真结果同样，在ChupScaling中，对5个点目标的回波信号进行了仿真，5个点目标构成了矩形的4个顶点和中

22、心，其坐标分别如下，格式为（方位向距离向后向反射系数）:1200011250-50112505011150-5011150501图13是仿真的雷达回波信号图仿真出来的信号500452535500050111002003004005006007008009001000方位向图13仿真出来的SAR回波信号图14是经过第一次相位校正之后，通过距离向压缩后的距离时域-方位时域信号图（ChirpScaling算法的七个步骤中并不包含该信号，该信号是将步骤2之后的信号通过方位向傅里叶逆变换，再进行距离向压缩得到的，只为了验证原理）。按照理论，该图中所有点的距离徙动都应该一样。从图中人致可看出，五个点的距离

23、徙动是差不多的。ChirpScaling后、经过距离向压缩，距离徙动一致501001520025300354004505001002003004005006007008009001000方位向图14ClurpScahng第2步之后、经过距离向压缩得到的图图15为步骤5之后，信号距离压缩，距离徙动校正之后的距离多普勒域中的信号图。消除距离徙动后的信号图15距离徙动校正之后的图图16为步骤6之后，消除相位偏移的图。相位枝正后的信号05050505112233445001002003004005006007008009001000方位向图16消除相位偏移的图图17为通过ChirpScaling算法生

24、成的点目标图像生U的須clear;clc;closeall;煤=勰常数定义C二3e8;勰雷达参数Fc=le9;lambda=C/Fc;勰日标区域参数Xmin二0;Xmax=50;Yc二10000;Y0二500;勰轨道参数V二100;图17通过ChirpScaling算法生成的点目标图像6.4几种算法比较本文讨论了距离多普勒算法和ChirpScaling算法，其中距离多普勒算法考虑了最近邻点插值和smc插值两种插值方法。距离多普勒算法兼顾了成熟、简单、高效和精确等因素，至今仍被广泛使用，但是距离多普勒算法有两点不足：首先，当用较长的核函数提高距离徙动校正(RCMC)精度时，运算量较大；其次，二次

25、距离压缩(SRC)对方位频率的依赖性问题较难解决，从而限制了其对某些大斜视角和长孔径SAR的处理精度。最近邻点插值的优点是速度快，该插值的运行时间为2.267137秒，缺点是不够精确；sine插值的优点是精确，该方法的运行时间为29.148728秒，缺点是速度慢：ChupScaling算法避免了插值运算，提高了速度，运行时间为0.323327秒，但是其算法较为复杂。勰二二二二二=勰文件名：NearSAR.m紈作者：徐一凡勰功能：合成孔径雷达距离多普勒算法点目标成像勰二二二二二=%光速%载频1GHz%波长%目标区域方位向范圉Xmin,Xmax%成像区域中线%目标区域距离向范Yc-Y0,Yc+YO

26、%成像宽度为2*Y0%SAR的运动速度100m/sH二5000;%高度5000mR0二sqrt(Yc"2+H"2);%最短距离號天线参数D=4;%方位向天线长度Lsar=lambda*RO/D;%SAR合成孔径长度，合成孔径雷达成像算法与实现P.100Tsar=Lsar/V;勰慢时间域参数Ka=-2*V2/lambda/R0;Ba=abs(Ka*Tsar);PRF=Ba;%SAR照射时间%多普勒频域调频率P.93%多普勒频率调制带宽%脉冲重复频率，PRF其实为多普勒频率的采样率，又为复频率，所以等于Ba.P.93PRT二1/PRF;ds二PRT;%脉冲重复时间%慢时域的时间

27、步氏Nslow=ceil(Xmax-Xmin+Lsar)/V7ds);%慢时域的采样数，ceil为取整函数，结合P.76的图理解Nslow=2"nextpow2(Nslow);%nextpow2为最壽近2的幕次函数，这里为fft变换做准备sn=linspace(Xmin-Lsar/2)/V,(Xmax+Lsar/2)/V,Nslow);%慢时间域的时间矩阵PRT=(Xmax-Xmin+Lsar)/V/Nslow;更新，周期减小了PRF二1/PRT;ds二PRT;勰快时间域参数设置制II于Nslow改变了，所以相应的一些参数也需要Tr=5e-6;Br=30e6;Kr=Br/Tr;Fsr

28、=2*Br;dt=l/Fsr;%脉冲持续时间5us%chirp频率调制带宽为30MHz%chirp调频率%快时域采样频率，为3倍的带宽%快时域采样间隔Rmin=sqrt(Yc-YO厂2+H"2);Rmax=sqrt(Yc+YO)2+H2+(Lsar/2)2):Nfast=ceil(2*(Rmax-Rmin)/C/dt+Tr/dt);%快时域的采样数量Nfast=2"nextpow2(Nfast);%更新为2的幕次，方便进行fft变换dt=(2*Rmax/C+Tr-2*Rmin/C)/Nfast;Fsr=l/dt;抵分辨率参数设置DY二C/2/Br;DX二D/2;勰点目标参数

29、设置Ntarget=5;%点目标格式x,y,反射系数sigmaZPtarget=Xmin,Yc-50*DY,1成一个矩形以及矩形的中心tm=linspace(2*Rmin/C,2*Rmax/C+Tr,Nfast);%快时域的离散时间矩阵%更新间隔%距离向分辨率%方位向分辨率%点目标的数量%点目标位置，这里设置了5个点目标，构Xmin+50*DX,Yc-50*DY,1Xmin+25*DX,Yc,1Xmin,Yc+5SDY,1Xmin+50*DX,Yc+50*DY,1;dispCParameters:J)%参数显示disp(,SamplingRateinfast-timedomain，);disp

30、(Fsr/Br)disp(,SamplingNumberinfast-timedomain，):disp(Nfast)disp('SamplingRateinslow-timedomain');disp(PRF/Ba)disp(,SamplingNumberinslow-timedomain，):disp(Nslow)disp('RangeResolution,);disp(DY)disp(，Cross-rangeResolution*);disp(DX)disp('SARintegrationlength*);disp(Lsar)dispC'Posit

31、ionoftargets):disp(Ptarget)%=勰生成回波信号K=Ntarget;%目标数目N=Nslow;%慢时域的采样数M=Nfast;%快时域的采样数T=Ptarget;%目标矩阵Srnm=zeros(N,M);%生成零矩阵存储回波信号fork=l:l:K%总共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,l)*tm-tau,*ones(1,

32、M);%(t-tau),其实就是时间矩阵，ones(N,1)和ones(1,M)都是为了将其扩展为矩阵phase=pi*Kr*Dfast."2-(4*pi/lambda)*(R5*ones(1,M)；%相位，公式参见P.96Srnm=Srnm+sigma*exp(j*phase)*(0Dfast&DfastTr)*(abs(Dslow)<Lsar/2)'*ones(l,M)；%由于是多个目标反射的回波,所以此处进行叠加end呱:號距离-多普勒算法开始眺距离向压缩tic;tr二tm-2*Rmin/C;Refr=exp(j*pi*Kr*tr.2).*(0<tr

33、&tr<Tr);Sr=ifty(fty(Srnm).*(ones(N,1)*conj(fty(Refr);Gr=abs(Sr);抵开始进行距离弯曲补偿正侧视没有距离走动项上要是因为斜距的变化引起回波包络的徙动抵补偿方法：最近邻域插值法，具体为：先变换到距离多普勒域，分别对单个像素点计算出距离徙动量，得到距离徙动量与距离分辨率的比值，號该比值可能为小数，按照四舍五入的方法近似为整数，而后在该像素点上减去徙动量勰方位向做fft处理再在频域做距离弯曲补偿Sa.RD二ftx(Sr);%方位向FFT变为距离多普域进行距离弯曲校正%距离徙动运算，由于是正侧视，fdc二0,只需要进行距离弯曲补

34、偿。Kp二1;%计算或者预设预滤波比h=waitbar(0,'最近邻域插值');抵首先计算距离迁移量矩阵forn=l：N%总共有N个方位采样form=l：M%每个方位采样上有M个距离采样delta_R=(1/8)*(lambda/V)(R0+(mM/2)*C/2/Fsr)*(n-N/2)*PRF/N)"2;E离迁移量P.160：(R0+(m-M/2)*C/2/Fsr):每个距离向点m的RO更新；(n-N/2)*PRF/N：不同方位向的多普勒频率不一样RN!C=2*delta_R*Fsr/C;%此处为delta_R/DY,距离徒动了几个距离单元delta_RMC二RMC

35、-round(RMC):%距离徒动量的小数部分辻m+round(RMC)>M%判断是否超出边界Sa_RD(n,m)=Sa_RD(n,M/2);elseifdeltjRMC二0.5%五入Sd_RD(n,m)=Sa_RD(n,m+round(RMC)+1);else%四舍Sd_RD(n,m)=Sa_RD(n,m+round(RMC);endendendwaitbar(n/N)endclose(h)%=二二二二=Sr_rmc=iftx(Sa_RD);%距离徙动校正后还原到时域Ga=abs(Sr_rmc);勰方位向压缩ta=sn-Xmin/V;Refa=exp(j*pi*Ka*ta.2).*(a

36、bs(ta)<Tsar/2);Sa=iftx(ftx(Sr_rmc).*(conj(ftx(Refa).'*ones(l,M);Gar=abs(Sa);toe;眺二二二二二二二二二二二二二二二二二=勰绘图colormap(gray);figure(l)subplot(211);row=tm*C/2-2008:col=sn*V-26;imagesc(row,col,255-Gr);%距离向压缩，未校正距离徙动的图像axis(Yc-YO,Yc+YO,Xmin-Lsar/2,Xmax+Lsar/2);xlabel('距离向'),ylabel('方位向')

37、,titleC距离向压缩，未校正距离徙动的图像')，subplot(212);imagesc(row,col,255Ga);%距离向压缩，校正距离徙动后的图像axis(Yc-YO,Yc+YO,Xmin-Lsar/2,Xmax+Lsar/2);xlabel('距离向'),ylabel('方位向'),titleC距离向压缩,校正距离徙动后的图像'),figure(2)colormap(gray);imagesc(row,col,255-Gar);%方位向压缩后的图像axis(Yc-YO,Yc+YO,Xmin-Lsar/2,Xmax+Lsar/2);x

38、label('距离向'),ylabel('方位向'),titleC方位向压缩后的图像')，勰二二二二二=二二二二=二二二=二二二二二二=二=二二=二二二二=二二二=二=二二=二=二二二=勰文件名：SincSAR.m紈作者：徐一凡勰功能：合成孔径雷达距离多普勒算法点目标成像勰二二二二二=二二二二=二二二=二二二二二二=二=二二=二二二二=二二二=二=二二=二=二二二=clear;clc;closeall;勰二二二二二=二二二二=二二二=二二二二二二=二=二二=二二二二=二二二=二=二二=二=二二二=勰常数定义C=3e8;勰雷达参数%光速Fc=le9;lam

39、bda=C/Fc;勰日标区域参数%载频1GHz%波长Xmin二0;%目标区域方位向范圉Xmin,XmaxXmax=50;Yc=10000;Y0二500;%成像区域中线%目标区域距离向范mYc-YO,Yc+YO%成像宽度为2*Y0勰轨道参数V二100;H二5000;R0=sqrt(Yc"2+H"2);紈天线参数D二4;%SAR的运动速度100m/s%高度5000m%最短距离%方位向天线长度Lsar=lambda*R0/D;%SAR合成孔径长度，合成孔径雷达成像算法与实现P.100Tsar=Lsar/V;%慢时间域参数%SAR照射时间Ka=-2*V2/lambda/R0;Ba=

40、abs(Ka*Tsar);PRF=Ba;%多普勒频域调频率P93%多普勒频率调制带宽%脉冲重复频率，PRF其实为多普勒频率的采样率，又为复频率，所以等于Ba.P.93PRT二1/PRF;ds二PRT;%脉冲重复时间%慢时域的时间步长Nslow=ceil(Xmax-Xmin+Lsar)/V/ds);%慢时域的采样数，ceil为取整函数，结合P.76的图理解Nslow=2"nextpow2(Nslow);%nextpow2为最常近2的幕次函数，这里为fft变换做准备sn=linspace(Xmin-Lsar/2)/V,(Xmax+Lsar/2)/V,Nslow);%慢时间域的时间矩阵PR

41、T=(Xmax-Xmin+Lsar)/V/Nslow;更新，周期减小了PRF二1/PRT;ds二PRT;勰快时间域参数设置%由于血100改变了，所以相应的一些参数也需要Tr=5e-6;Br=30e6;Kr=Br/Tr;Fsr=2*Br;dt=l/Fsr;%脉冲持续时间5us%chirp频率调制带宽为30MHz%chirp调频率%快时域采样频率，为3倍的带宽%快时域采样间隔Rmin=sqrt(Yc-YO)2+H2);Rmax=sqrt(Yc+YO)*2+H*2+(Lsar/2)*2);Nfast=ceil(2*(Rmax-Rmin)/C/dt+Tr/dt);%快时域的采样数量Nfast=2&qu

42、ot;nextpow2(Nfast);%更新为2的幕次，方便进行fft变换tm=linspace(2*Rmin/C,2*Rmax/C+Tr,Nfast);%快时域的离散时间矩阵dt=(2*Rmax/C+Tr-2*Rmin/C)/Nfast;%更新间隔Fsr=l/dt;號分辨率参数设置DY=C/2/Br;DX二D/2;勰点目标参数设置%距离向分辨率%方位向分辨率Ntarget=5;%点目标的数量%点目标格式x,y,反射系数sigmaZPtarget=Xmin,Yc-50*DY,1成一个矩形以及矩形的中心%点目标位置，这里设置了5个点目标，构Xmin+50*DX,Yc-50*DY,1Xmin+25

43、*DX,Yc,1Xmin,Yc+5SDY,1Xmin+50*DX,Yc+50*DY,1;dispCParameters:J)%参数显示disp(,SamplingRateinfast-timedomain，);disp(Fsr/Br)disp(,SamplingNumberinfast-timedomain，):disp(Nfast)disp(,SamplingRateinslow-timedomain');disp(PRF/Ba)dispC*SamplingNumberinslow-timedomain，):disp(Nslow)disp('RangeResolution,)

44、;disp(DY)disp(*Cross-rangeResolution'):disp(DX)disp('SARintegrationlength');disp(Lsar)disp('Positionoftargets*):disp(Ptarget)%=勰生成回波信号K=Ntarget;%目标数目N=Nslow;%慢时域的采样数M=Nfast;%快时域的采样数T=Ptarget;%目标矩阵Srnm=zeros(N,M);%生成零矩阵存储回波信号fork=l:l:K%总共K个目标sigma=T(k,3);%得到目标的反射系数Dslow=sn*V-T(k,1);%方

45、位向距离，投影到方位向的距离R=sqrt(Dslow."2+T(k,2)"2+H"2);%实际距离矩阵tau=2*R/C;%回波相对于发射波的延时Dfast=ones(N,l)*tm-tau,*ones(1,M);%(t-tau),其实就是时间矩阵，ones(N,1)和ones(1,M)都是为了将其扩展为矩阵phase=pi*Kr*Dfast."2-(4*pi/lambda)*(R5*ones(1,M)相位，公式参见P.96Srnm=Srnm+sigma*exp(j*phase)*(0<Dfast&Dfast<Tr)*(abs(Dsl

46、ow)<Lsar/2)'*ones(l,M);%由于是多个目标反射的回波,所以此处进行叠加end呱:勰距离-多普勒算法开始號距离向压缩tic;tr=tm-2*Rmin/C;Refr=exp(j*pi*Kr*tr."2).*(0<tr&tr<Tr);Sr=ifty(fty(Srnm).*(ones(N,l)*conj(fty(Refr):Gr=abs(Sr);抵开始进行距离弯曲补偿正侧视没有距离走动项上要是因为斜距的变化引起回波包络的徙动抵补偿方法：最近邻域插值法，具体为：先变换到距离多普勒域，分别对单个像素点计算出距离徙动量，得到距离徙动量与距离分辨

47、率的比值，抵该比值可能为小数，按照四舍五入的方法近似为整数，而后在该像素点上减去徙动量紈方位向做fft处理再在频域做距离弯曲补偿Sa_RD=ftx(Sr);%方位向FFT变为距离多普域进行距离弯曲校正%距离徙动运算，由于是正侧视，fdc二0,只需要进行距离弯曲补偿。Kp二1;%计算或者预设预滤波比紈第二种方法进行截断sine插值进行距离徒动校正h=waitbar(0,'Sine插值');P=4;%4点sine插值RNICmaxtix=zeros(N,M);forn=l:Nform=P：Mdelta_R=(1/8)*(lambda/V)2*(R0+(mM/2)*C/2/Fsr)*

48、(n-N/2)*PRF/N)“2;%首先计算距离迁移量计算方法就是把斜距变换到距离多普勒域就知道了RMC=2*de11a_R*Fsr/C:%距离徒动了几个距离单元delta_RIC=RMC-round(RMC);%距离徒动量的小数部分fori=-P/2:P/2-lifm+RMCi>M%判断是否超出边界RMCmaxtix(n,m)=RMCmaxtix(n,m)+Sa_RD(n,M)*sinc(pi*(-i-RMC);elseRMCmaxtix(n,m)=RMCmaxtix(n,m)+Sd_RD(n,m+round(RMC)+i)*sinc(pi*(-i+delta_RMC);endende

49、ndwaitbar(n/N)endclose(h)%=Sr_rme=iftx(RMCmaxtix);%距离徙动校正后还原到时域Ga=abs(Sr_rmc);勰方位向压缩ta=sn-Xmin/V;Refa=exp(j*pi*Ka*ta.2).(abs(ta)<Tsar/2);Sa=iftx(ftx(Sr_rmc).*(conj(ftx(Refa).'*ones(1,M);Gar=abs(Sa);toe;%=勰绘图colormap(gray);figure(l)subplot(211);row=tm*C/2-2008:col=sn*V-26;imagesc(row,col,255Gr

50、):%距离向压缩，未校正距离徙动的图像axis(Yc-YO,Yc+YO,Xmin-Lsar/2,Xmax+Lsar/2);xlabel('距离向'),ylabel('方位向')，titleC距离向压缩,未校正距离徙动的图像'),subplot(212);imagesc(row,col,255Ga):%距离向压缩，校正距离徙动后的图像&xis(Yc-YO,Yc+YO,Xmin-Lsar/2,Xmax+Lsar/21);xlabel('距离向'),ylabel('方位向')，titleC距离向压缩,校正距离徙动后的图像

51、'),figure(2)colormap(gray);imagesc(row,col,255-Gar);%方位向压缩后的图像axis(Yc-YO,Yc+YO,Xmin-Lsar/2,Xmax+Lsar/2);xlabel('距离向'),ylabel('方位向')，titleC方位向压缩后的图像')，%ChirpScaling算法%徐一凡clearall;clc;%距离向参数range：xdomainTr=200;%l时宽200mBr=l;%带宽1Kr=Br/Tr;%调频斜率Fc=4;%载频4Nfast=512;%为了快速运算Xc=1200;X0=

52、150;%定义距离向范围x二Xc+linspace(-X0,X0,Nfast);%x域序列：Xc-XO'Xc+XOdx=2*X0/Nfast;%定义步长kx=linspace(-l/dx/2,l/dx/2,Nfast);%kx域序列%方位向参数cross-range：ydomainTa二300;%时宽300m,合成孔径长度Ba=l;%带宽l(l/m)Ka二Fc/Xc;%调频斜率Ka二Ba/Ta二Fc/XcNslow二1024;%为了快速运算Y0二200;y=linspace(-Y0,Y0,Nslow);%y域序歹U:-YOYOdy=2*Y0/Nslow;ky=linspace(-l/d

53、y/2,l/dy/2,Nslow);%ky域序列勰日标几何关系targetgeometry紙坐标,y坐标，复后向散射系数Ptar=Xc,0,1+0jXc+50,"50,l+0jXc+50,50,HOjXc-50,-50,l+0jXc-50,50,HOj;disp(Positionoftargets*):disp(Ptar)勰生成SAR正交解调后的回波数据Srnm二zeros(Nfast,Nslow);N二size(Ptar,1);%目标个数h=waitbar(0,*SAR回波生成');fori=l:l:Nxn=Ptar(i,1);yn=Ptar(i,2);sigma二Ptar(i,3);%提取每个目标的信息X=x.'*ones(1,Nslow);%扩充为矩阵Y=ones(Nfast,l)*y；%扩充为矩阵DX=sqrt(xn2+(Y-yn).2);%中间变量phase二pi*Kr*(X-DX)."2-2*pi*Fc*DX;%回波相位Srnm=Srnm+sigma*exp(j*phase).*(abs(XDX)<Tr/2).*(