ISAR成像原理附matlabsimulation

1、ISAR 成像原理 附 matlab simulation1 . 成像原理逆合成孔径雷达成像是指目标运动而雷达不动， 利用二者之间的相对运动形成的 弧形合成孔径来得到方位向的高分辨率，距离向的高分辨率依然是经过大带宽信 号的脉冲压缩得到的。同样是相对运动，但由于 ISAR成像中，对于目标的运动 情况不能确知，所以合成阵列的分布是不可能准确测量的， 好在为得到亚米级的 横向分辨率，雷达对目标视线的变化只要很小几度， 在这期间，由于目标的惰性， 其姿态变化不可能十分复杂。逆合成孔径雷达在另一些方面要比合成孔径雷达简单，主要是目标的尺寸比合成孔径雷达所要观测的场景小的多，一般目标不超过十几米，大的也

2、只有百余米， 当目标位于几十千米以外时，电波的平面波假设总是成立的。现在来分析一下目标相对于雷达的运动。 可将目标的运动分解为转动和平动，目 标平动是指该目标上的参考点沿运动轨迹移动， 而目标相对于雷达射线的姿态保 持不变；而转动分量是指目标围绕该参考点转动。要获得高的距离向分辨率，则发射的脉冲较窄，通常为纳秒级的，而回波序列的 时延变化常比脉冲宽度大得多，所以这时的时序脉冲包络在时间上是错开的 (由 于通常要求在成像时间内转动分量引起的散射点在距离向的走动不超过一个距 离分辨单元)，所以认为一维距离像的错开主要是平动分量的影响。在这种情况 下，目标上同一散射点的各个回波经过距离分辨以后，将处

3、于不同的距离单元， 因此必须对平动分量引起的包络时延进行补偿，否则无法对方位向进行分辨。通 常的做法是以某一次回波为基准，而将各次的包络对齐(保持原包络的振幅和相 位不变，只是位置搬移)，再比较各次回波的相位变化，从而得到点目标的多普 勒。经过包络对齐处理，各次回波的距离单元已基本对齐，各距离单元回波包络序列 的幅度和相位的横向变化基本正常。而各次回波中还包含平动分量表现出来的初 相，为了进行方位向分辨，ISAR需要目标相对于雷达有等效的姿态转动。在理 想的转台成像模式下，目标的运动只有相对于雷达的转动，易于得到目标的二维 高分辨率图像，而实际中，平动分量的存在(尤其在载频很高的情况下)，平移

4、运动在方位向引入了二次的相位分量， 必须要消除其影响，否则在进行方位向分 辨时(FFT),会造成图像方位向散焦。所以，成像的具体步骤有三：(1) 距离向脉压；(2) 运动补偿，包括包络对齐和初相校正；(3) 方位向分辨的获得。2 . 目标与雷达位置参数：雷达位于直角坐标系的原点，目标位于(50000, 25000)处，且目标沿x轴负方 向以50m/s的速度飞行。成像角度为2.865度，方位向有128次回波采样。在经过距离向的脉冲压缩之后，进行包络对齐处理，采用的是相邻相关法，为了 减小积累误差的影响，在第六次回波之前的对齐处理是与此前所有已对齐好的回 波做相关，在取其偏移量的平均值，在第六次以

5、后的是与其前面的五次对齐。 使 用这种对齐方法得到包络对齐后的一维距离向如图 1所示；由于有积累误差的存 在，使得对齐结果产生偏移，即单点目标的一维距离向之后是初相校正：因只有一个点目标，采用的单特显点方法。可是现在用判断单 特显点的方法得到的居然根本不是包络所在处的特显点， 所以就采用了有包络存 在距离单元的平均距离单元作为初相的参考点。最后是方位向FFT得到二维ISAR图像。图片：1.JPGRe:ISAR距离多普勒单点目标成像算法仿真对应程序clear all;close all;C=3e8;% 光速Fc=1e10;% 载频lamda=C/Fc;%工作波长Tr=51.2e-6;%发射脉宽B

6、r=1e7;%发射信号频宽Kr=Br/Tr;theta=2.865/180*pi;% 成像所需转角x0=50000;%雷达位于坐标系原点，初始时刻目标位置y0=25000;rO=(x0A2+y0A2)A0.5;%目标与雷达的距离v=50;%目标速度pr=1.2*Br;sampr=1/pr;sn=0:sampr:Tr;time_acct=r0*theta/v;% 成像时间对应的目标运动时间Ba=1/time_acct;% 方位向的多普勒展宽%pa=2*Ba;%方位向的采样率%sampa=1/pa;% 方位向采样间隔tm=linspace(0,time_acct,128);M=length(tm)

7、;N=length(sn);x=x0-v*tm;R=sqrt(x.A2+y0A2);% 回波信号模型tau=2*R/C;Dfast=ones(M,1)*sn-tau*ones(1,N);phase=pi*Kr*Dfast.A2-(4*pi/lamda)*(R*ones(1,N);Srnm=exp(j*phase);Srf=exp(j*pi*Kr*sn.A2);% 参考信号Sr=ifty(fty(Srnm).*(ones(M,1)*conj(fty(Srf);%for i=1:128%Sr1(i,:尸fft(Srnm(i,:).*conj(fft(Srf);%end%for i=1:128% S

8、r(i,:-ifft(Sr1(i,:);%end%对距离压缩后的信号进行包络对齐Sr2=zeros(M,N);Sr2 com=zeros(M,N);Sr3=zeros(M,N*8);Sr3=(ifft(fft(abs(Sr).,N*8).;%插值(一维距离像)Sr4=(ifft(fft(Sr.,N*8).;%回波复包络Sr2(1,:)=Sr(1,:);temp1=zeros(1,8*N);% 中间变量temp2=zeros(1,8*N);%index_max=zeros(1,8*N);for i=1:127 %包络对齐diff=0;b=(Sr3(i+1,:);if(i6)for j=1:ia=

9、(Sr3(j,:);c=xcorr(a,b);maxv=max(c);index_max=find(c=maxv);diff=index_max-8*N+diff;enddiff=floor(diff/i);else for j=i-4:ia=(Sr3(j,:);c=xcorr(a,b);maxv=max(c);index_max=find(c=maxv);diff=index_max-8*N+diff;enddiff=floor(diff/5);endtemp1(diff+1:8*N)=Sr3(i+1,1:8*N-diff);temp2(diff+1:8*N)=Sr4(i+1,1:8*N-d

10、iff);if(diff=0)temp1(1:diff)=Sr3(i+1,8*N-diff+1:8*N);temp2(1:diff)=Sr4(i+1,8*N-diff+1:8*N);endSr3(i+1,:)=temp1;Sr4(i+1,:)=temp2;%indexr=zeros(1,M);%for i=1:M%temp=abs(Sr2(i,:);%maxr=max(temp);%indexr(i尸find(maxr=temp);ii%endendSr2=(ifft(fft(Sr3.,N).;% 从插值变换回原序列Sr2com=(ifft(fft(Sr4.,N).;indexr=zeros(

11、1,M);indexr_old=zeros(1,M);for i=1:M%统计各次回波的峰值temp=abs(Sr2(i,:);temp old=abs(Sr(i,:);maxr=max(temp);max_old=max(temp_old);indexr(i)=find(maxr=temp);indexr_old(i尸find(max_old=temp_old);iiend%tempa=zeros(1,N);%tempb=zeros(1,N); %for i=1:N%tempa=(sum(abs(Sr2com)./M).A2;%tempb=sum(abs(Sr2com2)./M; %end%index_specl=find(1-tempa./tempb)0.12);%特显点的判别准则在这里不适用index_spec=sum(indexr);% 初相矫正(只有一个点目标，采用的单特显点方法，但是现在的用判断单特显点的方法居然得到的根本不是包络所在处的特显点，%所以就采用了包络所在距离单元的平均距离单元作为初相的参考点)index_spec=floor(index_spec/M);match_fctr