哈工大微波成像技术

1、HarbinInstituteofTechnology微波成像技术实验报告课程名称：微波成像技术院系：电子与信息工程学院姓名：学号：授课教师：哈尔滨工业大学年月曰实验一多散射点一维距离像产生实验要求：根据ISAR转台成像原理，选取四至五个（或更多）散射点产生不同时刻的一维距离像,使用X波段的FLM脉冲信号，信号波长取3cm，带宽为100MHz实验内容：1分析距离像随姿态角的变化2取出1点描述成像系统性能实验原理：假设雷达发射信号为：广、ts(t,t)=rectTp丿ej2兀（ft+1y池），其中rect(u)=UT2整数，T脉冲重复周期，t=mT为慢时间。mf为中心频率，T为脉宽，y为调频率，

2、t=t-mT为快时间，m为cp假设某点目标到雷达的距离为R，则该点目标接收到的雷达信号为:s(t,t)=Arectrm设参考距离为R，则参考信号为:refft-2R/c、j2兀rectrefeTref丿s(t,t)=refm2令R=R-R，则解线频调的差频输出为:Arefs(t,t)=s(t,t)-s*(t,t)ifmrmrefmf-AIt-2R；cTIp丿=Arect对快时间作傅里叶变换便可得到点目标的一维距离像:S(f,t)=ATsincT(f+2-R)eJ一-JcAifimppi实验过程及结果:1离像随姿态角的变化20001015001000-10500-20-1010200050010

3、00u/m雷达与目标的位置关系如右上方图所示，雷达位于原点，目标几何中心相对雷达的坐标为1500(1200,1600),单位米，可知目标距雷达初始距离为2km。假设RVP项可以忽略或已经被补偿掉，在特定t解线频调的差频输出可表示为:ms(t;t)=Arectfmii=1t2R/c,i根据上式设计目标回波矩阵，再对快时间进行傅里叶变换便可得到目标的一维距离像如下图：仿真目标由7个点目标构成，如左下方图所示:0.60.40.20-40-2002040距离单元转角为0deg)度幅化归0.60.40.20-40-2002040距离单元转角为0.01deg)0.60.40.20-40-2002040距离

4、单元转角为0.05deg)度幅化归0.60.40.20-40-2002040距离单元转角为0.1deg)0.60.40.20-40-2002040距离单元转角为1deg)度幅化归0.60.40.20-40-2002040距离单元转角为3deg)度幅化归度幅化归度幅化归由不同转角下多散射点目标的以为距离像可以看出，转角在O.OIrad内多散射点目标的距离像包络基本没有变化，相关性极强。转角在0.01之后相关性逐渐减弱，一维距离像包络幅度存明显变化，但散射点的位置基本不变，没有发生越距离单元走动。可以利用平均距离像避免多散射点一维距离像幅度变化对后续成像的影响。多散射点目标最大横向距离差为20m，

5、由AO单个散射点描述成像系统性能将多点目标改为单点目标，单点目标相对雷达的坐标为(1209,1612)。点目标距雷达的初始距离为2015m，相对参考点距离15m，其余参数不变。下列图为不同信噪比情况下点目标的一维距离像加窗和不加窗时的对比，图中实线为加矩形窗的距离像，虚线为加汉明窗时的距离像。幅度归日IIii432亠OOOOO甘花X亠co兽第幅度Jd日幅度闭日O4OO芯对S5S舄禹届股豆SNRAS距离/m(SNR=3)i/mSNR=-10皿ESI在没有噪声的情况加，加汉明窗成功地将第一旁瓣降低到-45dB左右，与此同时主瓣展宽约1.3倍。当SNR较大为20dB时汉明窗依然能够将旁瓣抑制到-40

6、dB左右，然而随着SNR减小汉明窗的旁瓣随之升高，SNR为3时汉明窗的旁瓣高度约为-25dB。在SNR信噪比逐渐减小的过程中加窗后的主瓣略有展宽但变化不大。实验二多散射点二维SAR条带式成像实验要求:使用RD方法进行机载SAR多散射点二维条带式成像，可以不考虑距离徒动，天线方位向尺寸为D=2m，仰角尺寸为W=3m，载机速度V=150m/s。场景中心距离雷达6km，斜视角自定。实验内容：1用RD方法仿真多点目标二维SAR图像2取出1点描述成像系统性能实验原理：对雷达接收的任意一点目标Q,设此点目标到飞行航线的垂直距离(或称最近距离)为R,到B雷达相位中心的瞬时斜距为R(t;R)，则雷达接收的基频

7、信号在距离快时间-方位慢时间域mB(t-1域)可写为：mI，2R(t;R)a(t,Q)expC2R(t;R)atmBjntmBr(c丿am(c丿s(t,t;R)=mB4nexp-jjR(t;R)mB式中a()和a()分别为雷达线性调频(LFM)信号的窗函数和方位窗函数，前者在未加权时为ra矩形窗，后者除滤波加权外，还与天线波束形状有关，Y是发射的LFM信号的调频率，c为光速。对距离作匹配滤波压缩的参考函数为：srf(F)=a(F)exp(-jny九)rr快时间域的匹配滤波可在频率域采用FFT进行：s(t,t;R)=IFFTFFTs(f,t;R)FFTsref*(t,t;R)mBmBrmB若距离

8、向为矩形窗，则的接收信号通过上述处理后得：八s(t,t;R)=AsincmBr2R(t;R)Ba()expa-jTR(tm；RB)A为距离压缩后Q点的信号幅度，若线性调频信号的频带为B，有sinc(F)=sine伍BF)。r距离压缩完成后，下一步要进行方位处理，首先要检验距离徙动的影响，如为正侧视工作只要检验距离弯曲。在合成孔径期间，r=込匕(其中M通常取4或8)时距离弯曲可忽略。q32p2Ma假设上述条件满足，对最近距离为R的点目标Q，在t时刻的斜距为：BmR(t;R)=.；R2+(Vt)2=R+mBBmB2RB式中V为载机速度。距离快时间-方位慢时间域信号可写成：As(t,t;R)=Csi

9、ncmBr2RBca()exp.4n(jaL九I(Vt)2mRcBRB丄(Vt)2)+m2R丿B方位匹配滤波的参考函数为：sref(t;R)=a(t)exp(jnk12)amBrmam其中调频率为k=-空1，方位脉压在频率域进行比较方便，脉压后的输出为：a九RBs(t,t;R)=IFFTFFTs(2,t;R)-FFTsref*(t;R)mBmBamB若方位窗函数也是矩形，则上式可写成：sincGBt)dm八s(t,t;R)=Csinc用RD方法仿真多点目标二维SAR图像假设侧视角为30，地面目标由7个点目标构成，如左下方图所示:10iiIIiiII、气-5-rII1D-；-15一I600050

10、00400030002000-10DO-oO雷达与目标的位置关系如右上方图所示，图示坐标单位为米。雷达至场景中心线距离6km,雷达侧视角为30。确定雷达脉冲重复频率：为了不发生方位模糊，最低脉冲重复频率为九=2V/D=150Hz，为了不发生距离模糊，最高脉冲重复频率为RL/Wctan(90卩).OOAyrTTf=4.33MHz。选取PRF为600Hz。rh2R入0成像结果如下图所示:100150回波数据支撑区2002503003504004502004006008001000120014001600距离向距离压缩后的回波数据50100150200250300350400450200400600

11、80010001200距离向多散射点二维SAR条带式成像结果3单个散射点描述成像系统性能将多点目标改为单点目标，其余参数不变。距离向与方位向均不加窗的结果为：-100-100-；riL：L1瓯?OJQ2:1-.1-4!/II川密tl丁|fl1-nT7=距离问?041万伫向距离向与方位向均加汉明窗成像结果如图所示：单个散射点维朋尺条带式成像结果12jii匚_”hi-:加匡-60-万仿刁距离向-u-n-n-DJP融專u科-u-h.-n8uon-.-n4-u-6-u8-oon-2厂-,fi丫门茁ri4|Y|Ffl方位向距离向加汉明窗方位向不加窗:Q3p和方位向距离向U-6.O-S-nO-n-2-Os

12、o-404-UA-.o25n1向离5r距离向不加窗方位向加汉明窗:羽泮谢点二：tew拎港弍扃饥吉圭ytxjr喩客rllDIISI/III；-lzU4.1方也向U.；ill由实验结果可以看出，方位向加汉明窗后单散射点的二维像在方位向上旁瓣抵制-40dB左右同时主瓣有所展宽，距离向加汉明窗后单散射点的二维像在距离向上旁瓣抵制40dB左右同时主瓣有所展宽，两个维度上的旁瓣存在交叉影响，但单个散射点时交叉影响程度不大。实验程序实验一：clearall;closeall%散射点模型%posxy=-10-5000510;00-50500*2;posxy=912;scatter(posxy(:,1),pos

13、xy(:,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)RO=sqrt(uO42+vO42);%目标初始位置至雷达距离C=3e8;%电波传播速度lamda=0.03;%X波段Fc=C/lamda;%载频B=100e6;%发射信号带宽Tp=20e-6;%发射脉冲宽度gama=B/

14、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

15、(:,1)+u0;posxyi(:,2)=posxyi(:,2)+v0;Rri=sqrt(posxyi(:,1).A2+posxyi(:,2)42);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(iff

16、t(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(1,:);gridon;holdoff;axis(822-511)xlabel(距离/m(SNR=-10);ylabel(幅度/dB)%一维距离

17、像随角度的变化%figure,%subplot(321)%plot(echo2(2,200:280);%xlabel(距离单元(转角为Orad),ylabel(归一化幅度),axis(-infinf00.75)%subplot(322)%plot(echo2(12,200:280);%xlabel(距离单元(转角为0.01rad),ylabel(归一化幅度),axis(-infinf00.75)%subplot(323)%plot(echo2(103,200:280);%xlabel(距离单元(转角为0.1rad),ylabel(归一化幅度),axis(-infinf00.75)%subplo

18、t(324)%plot(echo2(203,200:280);%xlabel(距离单元(转角为0.2rad),ylabel(归一化幅度),axis(-infinf00.75)%subplot(325)%plot(echo2(303,200:280);%xlabel(距离单元(转角为0.3rad),ylabel(归一化幅度),axis(-infinf00.75)%subplot(326)%plot(echo2(403,200:280);%xlabel(距离单元(转角为0.4rad),ylabel(归一化幅度),axis(-infinf00.75)实验二：%雷达参数clearall;closeal

19、l;c=3e8;lamda=0.03;%X波段fc=c/lamda;%载频D=2;%雷达天线横向孔径L=3;%雷达天线纵向孔径Rb=6e3;%场景中心距phi=30/180*pi;%斜视角V=150;%载机飞行速度theta=lamda/L;%波束宽度beta=30/180*pi;%侧视角H=Rb*cos(beta);%载机高度Ls=theta*Rb;%合成孔径长度%距离向rhor=1;%方位向分辨率Tp=2.5e-6;%脉冲宽度B=c/2/rhor;%带宽Kr=B/Tp;%调频斜率Fr=1.2*B;%距离向采样频率Ymin=H*tan(beta-theta/2);Ymax=H*tan(bet

20、a+theta/2);gate=c/Fr/2;%距离分辨单元长度Rmin=floor(sqrt(YminA2+HA2)*gate;%最近距离Rmax=sqrt(YmaxA2+HA2);Nwg=ceil(Rmax-Rmin)/gate);%条带场景所需距离单元个数Np=floor(Tp*Fr);N_fast=Np+Nwg;N_fft=2Anextpow2(N_fast);t_fast=(-Np/2+1:Np/2+Nwg)/Fr;%快时间%方位向rhoa=D/2;%方位向分辨率Ka=2*V2/lamda/Rb;%多普勒调频率S=Ls/V;%有效观测标称时间Bd=abs(Ka)*S;%多普勒频宽PR

21、F=600;%方位向采样频率Tr=1/PRF;%脉冲重复间隔N_slow=ceil(2*Ls)/Tr/V);%方位向采样点数t_slow=(-N_slow/2:N_slow/2-1)*(1/PRF);%方位向时间域%回波数据仿真%P=100;00;-100;0-10;010;0-20;020;P=00;Ntarget=size(P,1);%figure,scatter(P(:,1),P(:,2),ok);%axisequal;axistight;gridon%pos=00H;Rb*sin(beta)+P(:,1)P(:,2)zeros(Ntarget,1);%figure,scatter3(p

22、os(:,1),pos(:,2),pos(:,3),ok);sr1=zeros(N_slow,N_fast);sr=zeros(N_slow,N_fast);R=zeros(1,N_slow);forn=1:Ntargetfork=1:N_slowR(1,k)=sqrt(Rb+P(n,2)八2+(t_slow(1,k)*V-P(n,1)八2);temp_win=find(abs(t_fast-2*(R(1,k)-Rmin)/c)Tp/2);win_r1=zeros(size(t_fast);win_r1(temp_win)=hamming(length(temp_win);%对回波数据加窗%w

23、in_r1(temp_win)=ones(size(temp_win);sr(k,:)=sr(k,:)+exp(-j*4*pi*R(1,k)/lamda+j*pi*Kr*(t_fast-2*(R(1,k)-Rb)/c).八2).*.win_r1.*(abs(t_slow(1,k)*V-P(n,1)(Ls/2);endend%figure,imagesc(real(sr);%title(回波数据支撑区)%xlabel(距离向)；ylabel(方位向);t=(-Np/2:Np/2-1)*(1/Fr);win_r2=hamming(length(t);%win_r2=ones(size(t);h=w

24、in_r2.*exp(j*pi*Kr*t.入2).*(abs(t)Tp/2);%距离向参考函数hf=fft(h,N_fft);fork=1:N_slowsrf(k,:)=fft(sr(k,:),N_fft);sof(k,:)=srf(k,:).*conj(hf);so(k,:)=ifft(sof(k,:);%距离向压缩后输出endfigure;plot(db(so(160,:)/max(so(200,:);%imagesc(abs(so(:,1:Nwg);%xlabel(距离向);ylabel(方位向);%title(距离压缩后的回波数据);figure;SSo=20*log10(abs(so)/max(max(a