ISAR成像原理附matlabsimulation文档格式.docx
《ISAR成像原理附matlabsimulation文档格式.docx》由会员分享,可在线阅读,更多相关《ISAR成像原理附matlabsimulation文档格式.docx(5页珍藏版)》请在冰豆网上搜索。
雷达位于直角坐标系的原点,目标位于(50000,25000)处,且目标沿x轴负方向以50m/s的速度飞行。
成像角度为2.865度,方位向有128次回波采样。
`eMrP`
在经过距离向的脉冲压缩之后,进行包络对齐处理,采用的是相邻相关法,为了减小积累误差的影响,在第六次回波之前的对齐处理是与此前所有已对齐好的回波做相关,在取其偏移量的平均值,在第六次以后的是与其前面的五次对齐。
使用这种对齐方法得到包络对齐后的一维距离向如图1所示;
由于有积累误差的存在,使得对齐结果产生偏移,即单点目标的一维距离向]689Q%D
之后是初相校正:
因只有一个点目标,采用的单特显点方法。
可是现在用判断单特显点的方法得到的居然根本不是包络所在处的特显点,所以就采用了有包络存在距离单元的平均距离单元作为初相的参考点。
;
I!
+lx3[
最后是方位向FFT得到二维ISAR图像。
-A~;
MGY
图片:
1.JPG
4DQ07w
Re:
ISAR距离-多普勒单点目标成像算法仿真对应程序
clearall;
B/kn&
^z$|~
closeall;
Jd"
s~n<
>
K
~U6YN_W
C=3e8;
%光速?
$_4R<
Fc=1e10;
%载频ptni'
W3
lamda=C/Fc;
%工作波长810<
1NP
Tr=51.2e-6;
%发射脉宽zTT
Br=1e7;
%发射信号频宽T&
"
dBoUq>
G
Kr=Br/Tr;
pQ(eF0KG
theta=2.865/180*pi;
%成像所需转角m&
0BbyE.z
x0=50000;
%雷达位于坐标系原点,初始时刻目标位置C!
k9JAa$Z
y0=25000;
EZumJ."
r0=(x0^2+y0^2)^0.5;
%目标与雷达的距离vcJb\LW
v=50;
%目标速度D6)Cjc>
a
JBEgiQ/
pr=1.2*Br;
s?
x>
Yl%
sampr=1/pr;
skDk/-*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;
%方位向采样间隔tGO[A#9a
tm=linspace(0,time_acct,128);
eMl]tdrI
M=length(tm);
=Q8H]F
N=length(sn);
0_j{ig
q]i(CaKh
x=x0-v*tm;
Ns!
3-Y
R=sqrt(x.^2+y0^2);
%回波信号模型0fPHh>
u
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]aK'
Srf=exp(j*pi*Kr*sn.^2);
%参考信号4Ay`rG
Sr=ifty(fty(Srnm).*(ones(M,1)*conj(fty(Srf))));
:
e&
n.i^
%fori=1:
1281]%W\RHxo
%Sr1(i,:
)=fft(Srnm(i,:
)).*conj(fft(Srf));
n-hvh-ZO
%endzD)IU_GWa
128eg2U+g4
%Sr(i,:
)=ifft(Sr1(i,:
));
fKTDt%
%endb\?
7?
g
%对距离压缩后的信号进行包络对齐YA"
Ti9-EV
Sr2=zeros(M,N);
~&
DB!
6*
Sr2com=zeros(M,N);
SLdN.4idK
Sr3=zeros(M,N*8);
=sE2}/g
Sr3=(ifft(fft(abs(Sr).'
N*8))).'
;
%插值(一维距离像)Xq,UV
Sr4=(ifft(fft(Sr.'
%回波复包络g5tjj.
Sr2(1,:
)=Sr(1,:
);
t|UnG
temp1=zeros(1,8*N);
%中间变量IpWy)B>
Fl3
temp2=zeros(1,8*N);
Lj&
1K~U
%index_max=zeros(1,8*N);
d8x<
|D
fori=1:
127%包络对齐~K~b`|1
diff=0;
c=A(o
b=(Sr3(i+1,:
.D8~)ZWN
if(i<
6)3h9_Sz8
forj=1:
iYgE]d?
_h
a=(Sr3(j,:
`Yx-~y5X
c=xcorr(a,b);
Sl3Kp_Z
maxv=max(c);
.FIt.XPzv
index_max=find(c==maxv);
_PXGAS
diff=index_max-8*N+diff;
(jU/Wj!
q
end7I_1#O
diff=floor(diff/i);
#yW.o'
S+
elseforj=i-4:
i9z5\*bs
-1qZqU$h
ecr886
&
|v{#,ymeb
3gb|x?
M]}l^m>
L
endyGU.AM
diff=floor(diff/5);
^v_J"
-{
endNF6X-,cd
1hl]W+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=%H'
p
temp2(1:
diff)=Sr4(i+1,8*N-diff+1:
/QQ8.8=5
end7_`L$<
-n
Sr3(i+1,:
)=temp1;
!
)N|J$FU
Sr4(i+1,:
)=temp2;
CdY8#+"
%indexr=zeros(1,M);
#kyl?
E
Mr#A_RZ2_~@
%temp=abs(Sr2(i,:
$[g_=Z
%maxr=max(temp);
r!
7Y'
|
%indexr(i)=find(maxr==temp);
oK2pM18
%endb/t
end|g%mP1O
Sr2=(ifft(fft(Sr3.'
N))).'
%从插值变换回原序列XAPYpBgm
Sr2com=(ifft(fft(Sr4.'
3\AM=`
indexr=zeros(1,M);
\"
.L
indexr_old=zeros(1,M);
{XgnZ`*
M%统计各次回波的峰值_d@=nK)
temp=abs(Sr2(i,:
r[s!
F=^
temp_old=abs(Sr(i,:
G/vC~6x
maxr=max(temp);
`2`fiKm
max_old=max(temp_old);
Q/ZkW
indexr(i)=find(maxr==temp);
Tr&
}$kird
indexr_old(i)=find(max_old==temp_old);
XX6T$pA6
end0)|Q6*E>
F2~%zNe
%tempa=zeros(1,N);
*|?
x86~
%tempb=zeros(1,N);
R^*K6Ad
N$Z?
\>
K0i
%tempa=(sum(abs(Sr2com))./M).^2;
]3B%8
%tempb=sum(abs(Sr2com).^2)./M;
5\Fz!
%endQvLZg
%index_specl=find((1-tempa./tempb)<
0.12);
%特显点的判别准则在这里不适用index_spec=sum(indexr);
%初相矫正(只有一个点目标,采用的单特显点方法,但是现在的用%判断单特显点的方法居然得到的根本不是包络所在处的特显点,"
&
~0T#
%所以就采用了包络所在距离单元的平均距离单元作为初相的参考点)e-5?
p~>
index_spec=floor(index_spec/M);
ir16
match_fctr=Sr2com(:
index_spec)./abs(Sr2com(:
index_spec));
g~S)aU\:
match_fctr=match_fctr*ones(1,N);
GiFXX
Sra=Sr2com.*conj(match_fctr);
cK`"
lxO
Srapd=ftx(Sra);
x;
\wY'
figure,imagesc(abs(Srapd));
7=XL!
:
P
t=abs(Srapd);
4~z?
t=write(t,'
1.bmp'
($nrqAv4
k=1;
WelcomeTo
Download!
!
欢迎您的下载,资料仅供参考!