ISAR成像原理附matlabsimulation文档格式.docx

1、雷达位于直角坐标系的原点，目标位于（50000，25000）处，且目标沿x轴负方向以50m/s的速度飞行。成像角度为2.865度，方位向有128次回波采样。 eMrP 在经过距离向的脉冲压缩之后，进行包络对齐处理，采用的是相邻相关法，为了减小积累误差的影响，在第六次回波之前的对齐处理是与此前所有已对齐好的回波做相关，在取其偏移量的平均值，在第六次以后的是与其前面的五次对齐。使用这种对齐方法得到包络对齐后的一维距离向如图1所示；由于有积累误差的存在，使得对齐结果产生偏移，即单点目标的一维距离向 689Q%D 之后是初相校正：因只有一个点目标，采用的单特显点方法。可是现在用判断单特显点的方法得到的

2、居然根本不是包络所在处的特显点，所以就采用了有包络存在距离单元的平均距离单元作为初相的参考点。 ;I!+ lx3 最后是方位向FFT得到二维ISAR图像。 -A;MGY 图片:1.JPG4DQ07w Re:ISAR距离多普勒单点目标成像算法仿真 对应程序clear all; B/kn&z$| close all; JdsnK U6YN_ W C=3e8;%光速 ?$_4R Fc=1e10;%载频 ptniW3 lamda=C/Fc;%工作波长 810G Kr=Br/Tr; pQ（eF0KG theta=2.865/180*pi;%成像所需转角 m&0BbyE.z x0=50000;%雷达位于坐

3、标系原点，初始时刻目标位置 C!k9JAa$Z y0=25000; EZumJ.r0=（x02+y02）0.5;%目标与雷达的距离 vcJbLW v=50;%目标速度 D6）Cjca JBEgiQ/ pr=1.2*Br; s?xYl % sampr=1/pr; skD k/-*R sn=0:sampr:Tr; qK xj time_acct=r0*theta/v;%成像时间对应的目标运动时间 A）3_6 Ba=1/time_acct;%方位向的多普勒展宽 （GcKaUg8* %pa=2*Ba;%方位向的采样率 ?3PJVy?%sampa=1/pa;%方位向采样间隔 tGOA#9a tm=lin

4、space（0,time_acct,128）; eMltd rI M=length（tm）; =Q8HF N=length（sn）;0_jig qi（CaKh x=x0-v*tm; Ns!3- Y R=sqrt（x.2+y02）;%回波信号模型 0fPHhu tau=2*R/C; *）kx Dfast=ones（M,1）*sn-tau*ones（1,N）; _jq（D, phase=pi*Kr*Dfast.2-（4*pi/lamda）*（R*ones（1,N）; Y7kb1UG Srnm=exp（j*phase）; Z aKSrf=exp（j*pi*Kr*sn.2）;%参考信号 4AyrG S

5、r=ifty（fty（Srnm）.*（ones（M,1）*conj（fty（Srf）; :e&n.i %for i=1:128 1 %WRHxo % Sr1（i,:）=fft（Srnm（i,:）.*conj（fft（Srf）; n-hvh-ZO %end zD）IU_GWa 128 eg2U+g4 % Sr（i,:）=ifft（Sr1（i,:）; fKT Dt% %end b?7?g %对距离压缩后的信号进行包络对齐 YATi9-EV Sr2=zeros（M,N）; &DB!6* Sr2com=zeros（M,N）; SLdN.4idK Sr3=zeros（M,N*8）; =sE2/g Sr3

6、=（ifft（fft（abs（Sr）.,N*8）.;%插值（一维距离像） Xq,UV Sr4=（ifft（fft（Sr.%回波复包络 g5tjj. Sr2（1,:）=Sr（1,:）;t|Un G temp1=zeros（1,8*N）;%中间变量 IpWy）BFl3 temp2=zeros（1,8*N）; Lj&1KU %index_max=zeros（1,8*N）;d8x|D for i=1:127 %包络对齐 Kb|1 diff=0; c=A（o b=（Sr3（i+1,: .D8）ZWN if（iL end yGU .AM diff=floor（diff/5）; v_J- end NF6X-

7、 ,c d 1hlW+9 temp1（diff+1:8*N）=Sr3（i+1,1:8*N-diff）; z z*Wz temp2（diff+1:8*N）=Sr4（i+1,1: e # t3u_ if（diff=0） b9!.-8y temp1（1:diff）=Sr3（i+1,8*N-diff+1:8*N）; XC=%Hp temp2（1:diff）=Sr4（i+1,8*N-diff+1: /QQ8.8=5 end 7 _L$F2 %zNe %tempa=zeros（1,N）;*|?x86 %tempb=zeros（1,N）; R*K6Ad N $Z?K0i %tempa=（sum（abs（Sr

8、2com）./M）.2; 3B%8 %tempb=sum（abs（Sr2com）.2）./M; 5Fz!%end QvLZg %index_specl=find（1-tempa./tempb）index_spec=floor（index_spec/M）; ir16 match_fctr=Sr2com（:,index_spec）./abs（Sr2com（:,index_spec）; gS）aU:, match_fctr=match_fctr*ones（1,N）; Gi FXX Sra=Sr2com.*conj（match_fctr）; cKlxO Srapd=ftx（Sra）; x;wYfigure,imagesc（abs（Srapd）; 7=XL !:P t=abs（Srapd）; 4 z?t=write（t,1.bmp （$nrqAv4 k=1;Welcome ToDownload !欢迎您的下载，资料仅供参考！