哈工大微波成像技术..

HarbinInstituteofTechnology微波成像技术实验报告课程名称：微波成像技术院系：电子与信息工程学院姓名：学号：授课教师：哈尔滨工业大学年月日实验一多散射点一维距离像产生实验要求：根据ISAR转台成像原理，选取四至五个（或更多）散射点产生不同时刻的一维距离像，使用X波段的FLM脉冲信号，信号波长取3cm，带宽为100MHz实验内容：1分析距离像随姿态角的变化2取出1点描述成像系统性能实验原理：假设雷达发射信号为：212ˆ2()ˆˆ(,)rectcjfttmptstteT，其中12121rect()0uuu，cf为中心频率，pT为脉宽，为调频率，ˆ=ttmT为快时间，m为整数，T脉冲重复周期，mTtm为慢时间。假设某点目标到雷达的距离为R，则该点目标接收到的雷达信号为：21222ˆ2ˆ2ˆ(,)rectcRRjfttccrmptRcsttAeT设参考距离为refR，则参考信号为：21222ˆ2ˆ2ˆ(,)rectrefrefcRRjfttccrefrefmreftRcstteT令refRRR，则解线频调的差频输出为：22*2444ˆ()ˆˆˆ(,)(,)(,)ˆ2rectrefcifmrmrefmRjRjtRjfRccccpsttsttstttRcAeeeT对快时间作傅里叶变换便可得到点目标的一维距离像：2244(,)sinc(2)cjRjfRccifimppiSftATTfReec实验过程及结果：1离像随姿态角的变化仿真目标由7个点目标构成，如左下方图所示：0500100015000500100015002000u/mv/m雷达与目标的位置关系如右上方图所示，雷达位于原点，目标几何中心相对雷达的坐标为（1200,1600），单位米，可知目标距雷达初始距离为2km。假设RVP项可以忽略或已经被补偿掉，在特定mt解线频调的差频输出可表示为：,24ˆ()+1ˆ2ˆ(;)rectrefciRNjtfRccifmiiptRcsttAeT根据上式设计目标回波矩阵，再对快时间进行傅里叶变换便可得到目标的一维距离像如下图：-40-200204000.20.40.6距离单元(转角为0deg)归一化幅度-40-200204000.20.40.6距离单元(转角为0.01deg)归一化幅度-40-200204000.20.40.6距离单元(转角为0.05deg)归一化幅度-40-200204000.20.40.6距离单元(转角为0.1deg)归一化幅度-40-200204000.20.40.6距离单元(转角为1deg)归一化幅度-40-200204000.20.40.6距离单元(转角为3deg)归一化幅度由不同转角下多散射点目标的以为距离像可以看出，转角在0.01rad内多散射点目标的距离像包络基本没有变化，相关性极强。转角在0.01之后相关性逐渐减弱，一维距离像包络幅度存明显变化，但散射点的位置基本不变，没有发生越距离单元走动。可以利用平均距离像避免多散射点一维距离像幅度变化对后续成像的影响。多散射点目标最大横向距离差为20m，由8L知最大相关积累角应约小于0.01度，这与实验结果是一致的。2单个散射点描述成像系统性能将多点目标改为单点目标，单点目标相对雷达的坐标为（1209,1612）。点目标距雷达的初始距离为2015m，相对参考点距离15m，其余参数不变。下列图为不同信噪比情况下点目标的一维距离像加窗和不加窗时的对比，图中实线为加矩形窗的距离像，虚线为加汉明窗时的距离像。在没有噪声的情况加，加汉明窗成功地将第一旁瓣降低到-45dB左右，与此同时主瓣展宽约1.3倍。当SNR较大为20dB时汉明窗依然能够将旁瓣抑制到-40dB左右，然而随着SNR减小汉明窗的旁瓣随之升高，SNR为3时汉明窗的旁瓣高度约为-25dB。在SNR信噪比逐渐减小的过程中加窗后的主瓣略有展宽但变化不大。实验二多散射点二维SAR条带式成像实验要求：使用RD方法进行机载SAR多散射点二维条带式成像，可以不考虑距离徒动，天线方位向尺寸为D=2m，仰角尺寸为W=3m，载机速度V=150m/s。场景中心距离雷达6km，斜视角自定。实验内容：1用RD方法仿真多点目标二维SAR图像2取出1点描述成像系统性能实验原理：对雷达接收的任意一点目标Q，设此点目标到飞行航线的垂直距离(或称最近距离)为BR，到雷达相位中心的瞬时斜距为(;)mBRtR，则雷达接收的基频信号在距离快时间-方位慢时间域(ˆmtt域)可写为：22(;)2(;)4πˆˆˆ(,;)(,)expπexp(;)ccmBmBmBrammBRtRRtRsttRatatQjtjRtR式中()ra和()aa分别为雷达线性调频(LFM)信号的窗函数和方位窗函数，前者在未加权时为矩形窗，后者除滤波加权外，还与天线波束形状有关，是发射的LFM信号的调频率，c为光速。对距离作匹配滤波压缩的参考函数为：2ˆˆˆ()()exp(π)rrsreftatjt快时间域的匹配滤波可在频率域采用FFT进行：*ˆˆˆ(,;)IFFT{FFT[(,;)]FFT[(,;)]}mBmBrmBsttRsttRsrefttR若距离向为矩形窗，则的接收信号通过上述处理后得：2(;)4πˆˆ(,;)sinc()exp(;)mBmBramBRtRsttRAtajRtRcA为距离压缩后Q点的信号幅度，若线性调频信号的频带为B，有ˆˆsinc()sinc()rtBt。距离压缩完成后，下一步要进行方位处理，首先要检验距离徙动的影响，如为正侧视工作，只要检验距离弯曲。在合成孔径期间，2232BrqaRrM（其中M通常取4或8）时距离弯曲可忽略。假设上述条件满足，对最近距离为BR的点目标Q，在mt时刻的斜距为：222()(;)()2mmBBmBBvtRtRRVtRR式中V为载机速度。距离快时间-方位慢时间域信号可写成：BmBaBmBrBmRVtRjacRVtcRtRtts2)(π4exp)()(2ˆCsinc);,ˆ(22方位匹配滤波的参考函数为：2(;)()exp(π)amBrmamsreftRatjkt其中调频率为22aBVkR，方位脉压在频率域进行比较方便，脉压后的输出为：*ˆˆ(,;)IFFT{FFT[(,;)]FFT[(;)]}mBmBamBsttRsttRsreftR若方位窗函数也是矩形，则上式可写成：2ˆˆ(,;)CsincsincBmBdmRsttRBtBtc对距离徙动不考虑的情况，通过对接收的二维信号，通过简单的在距离和方位分别进行线性调频信号的匹配滤波，就可实现对场景的二维成像。实验过程及结果：1用RD方法仿真多点目标二维SAR图像假设侧视角为30°，地面目标由7个点目标构成，如左下方图所示：雷达与目标的位置关系如右上方图所示，图示坐标单位为米。雷达至场景中心线距离6km，雷达侧视角为30°确定雷达脉冲重复频率：为了不发生方位模糊，最低脉冲重复频率为2/150RLfVDHz，为了不发生距离模糊，最高脉冲重复频率为0tan(90)=4.33MHz2RHWcfR。选取PRF为600Hz。成像结果如下图所示：回波数据支撑区距离向方位向200400600800100012001400160050100150200250300350400450距离向方位向距离压缩后的回波数据2004006008001000120050100150200250300350400450距离/m方位/m多散射点二维SAR条带式成像结果50005100520053005400550056005700580059006000-60-40-2002040603单个散射点描述成像系统性能将多点目标改为单点目标，其余参数不变。距离向与方位向均不加窗的结果为：距离向与方位向均加汉明窗成像结果如图所示：距离向加汉明窗方位向不加窗：距离向不加窗方位向加汉明窗：由实验结果可以看出，方位向加汉明窗后单散射点的二维像在方位向上旁瓣抵制-40dB左右同时主瓣有所展宽，距离向加汉明窗后单散射点的二维像在距离向上旁瓣抵制-40dB左右同时主瓣有所展宽，两个维度上的旁瓣存在交叉影响，但单个散射点时交叉影响程度不大。实验程序实验一：clearall;closeall%%散射点模型%posxy=[-10-5000510;00-50500]'*2;posxy=[912];scatter(posxy(:,1),posxy(:,2),'.')axisequalgridon%%参数设定M=512;%%%%距离向采样点数N=512;%%%%积累脉冲数PRF=600;%%%%脉冲重复频率D=size(posxy,1);%%%%散射点个数u0=1200;%%%%初始位置v0=1600;%%%%初始位置figure,scatter([posxy(:,1)+u0;0],[posxy(:,2)+v0;0])xlabel('u/m'),ylabel('v/m')R0=sqrt(u0.^2+v0.^2);%%%%目标初始位置至雷达距离C=3e8;%%%%电波传播速度lamda=0.03;%%%%X波段Fc=C/lamda;%%%%载频B=100e6;%%%%发射信号带宽Tp=20e-6;%%%%发射脉冲宽度gama=B/Tp;%%%%调频斜率fs=M/Tp;%%%%距离向采样率t_fast=(0:1/fs:(Tp-1/fs))-Tp/2;%%%%快时间t_slow=(0:1/PRF:(N-1)/PRF);%%%%慢时间mtheta=0.512/180*pi;s=zeros(N,M);thetai=(0:N-1)*mtheta/N;fori=1:N%%%%创建回波矩阵st=zeros(1,M);posxyi=[posxy(:,1)*cos(thetai(i))+posxy(:,2)*sin(thetai(i))...-posxy(:,1)*sin(thetai(i))+posxy(:,2)*cos(thetai(i))];posxyi(:,1)=posxyi(:,1)+u0;posxyi(:,2)=posxyi(:,2)+v0;Rri=sqrt(posxyi(:,1).^2+posxyi(:,2).^2);forii=1:Dst=st+exp(-sqrt(-1)*4*pi/C*(Fc+gama*(t_fast-R0/C))*(Rri(ii)-R0));ends(i,:)=st;s(i,:)=awgn(s(i,:),-10);end%%rangeprofiledecho1_ham=s*diag(hamming(M));%echo1_ham=s.*hamming(M)';echo1_rec=s;M2=16*M;echo2_ham=abs(fftshift(ifft(echo1_ham,M2,2),2));echo2_rec=abs(fftshift(ifft(echo1_rec,M2,2),2));%%加窗不加窗成像系统性能比较ff=((0:M2-1)-M2/2)/M2*fs/2/gama*C;figure,holdonplot(ff,20*log10(echo2_ham(1,:)/max(echo2_ham(1,:))),'--');plot(ff,20*log10(echo2_rec(1,:)/max(echo2_rec(