**
合成孔径雷达
**
1雷达的基本原理及应用
1.1基本原理
雷达,英文Radar(Radio detection and ranging ),是利用电磁波进行物体探测及测距的一种方法。通过发射电磁波对目标进行照射并接收回波,从而得到目标到电磁波发射点的距离,方位,海拔等信息。雷达的应用范围很广泛,例如:天气雷达,防空雷达,以及阅兵时出现的火控雷达等。
1.2雷达的天线
雷达可以“嗅到”目标的气息,是因为他有一双明亮的大眼睛——天线。
如图1所示,从天线一段到另一端点间线段的长度表示增益的大小,最长的那条线段叫做最大增益,而处增益与最大增益的夹角叫做天线波束角。超出波束角的范围,天线增益越来越小。天线增益和天线波束角之间的关系为:
1.3雷达的分辨率
描述雷达性能好坏的因素之一叫做分辨率,见图2。它表现了雷达能否精确的识别一个目标。这里,我们主要讨论距离分辨率和方位分辨率。
方位分辨率的好坏取决于天线波束角的大小和目标与发射点的距离。天线波束角越小,代表雷达扫描一次该区域时,扫过的区域越小,则更能精确地区分目标。而且,目标与电磁波发射点距离越近,雷达探测一次扫过的区域越小,区分目标越精确。雷达的方位分辨率为:
距离分辨率:在雷达图像中,当两个目标位于同一方位角时,但与雷达的距离不同时,二者被雷达区分出来的最小距离则是距离分辨率。它通常定义为:当较近目标回波脉冲的后沿(下降沿)与较远目标回波的前沿(上升沿)刚好重合时,作为可分辨的极限。如图3,此时两目标间的距离就是距离分辨率。
2 合成孔径雷达系统
我们设计了一个合成孔径雷达系统,它采用了基于条带式方式的距离多普勒算法。合成孔径雷达是利用与目标作用相对运动的小孔径天线,把在不同位置接收的回波进行相干处理,从而获得较高分辨力的成像雷达。我们设计了直观的合成孔径雷达系统帮助理解。见图4。
雷达信号的距离向高分辨通过发射宽带信号,方位向高分辨通过合成孔径处理。上图表示,一携带雷达的飞机,飞过目标区域(100m100m)附近,要找到目标区域中的一个目标点。 式中,Ls为合成孔径的有效长度,它是目标在天线波瓣宽度之内时飞机飞过的距离。
要实现雷达成像,需要经过以下几个步骤:
1.计算合成孔径长度Ls,并确定飞机在天线波瓣宽度之内飞行时间
2.回波采样
3.对回波进行距离压缩
4.方位压缩(获得图像)
整个设计的实验参数为:
载频f0, 发送带宽B1, 脉宽Tp, 飞机飞行速度Vp, 采样频率fs, 发射脉冲fp, 目标区域参数:R0 , x0 , y0 ,天线半波角θ0。
参数设置见代码。
3 线性调频信号(LFM)
脉冲压缩雷达能同时提高雷达的作用距离和距离分辨率。这种体制采用宽脉冲发射以提高发射的平均功率,保证足够大的作用距离;而接受时采用相应的脉冲压缩算法获得窄脉冲,以提高距离分辨率,较好的解决雷达作用距离与距离分辨率之间的矛盾。
脉冲压缩雷达最常见的调制信号是线性调频(Linear Frequency Modulation)信号,接收时采用匹配滤波器(Matched Filter)压缩脉冲。
这里,我们把雷达接收到的LFM信号的表达为:
5 脉冲压缩算法
5.1距离压缩信号
雷达的距离分辨率一般取决于信号带宽,普通脉冲雷达的时宽脉宽满足:
而在我们这次实验中,未压缩信号的脉宽为20微秒,经过匹配滤波器压缩信号后,回波信号脉宽为0.01微秒,大大减小了脉宽,使得回波信号在时域上的距离分辨率更高。
我们用Matlab实现雷达信号的距离压缩时,发现,信号在压缩处理后,变成了条状图样,见图5,这也表现了同一方位上的信号在距离上实现了压缩。
可以发现,经过匹配滤波器压缩处理后的图像分辨率大大提升,可以很好的区分物体,而且,考虑到徙动干扰较小,我们没有加入徙动矫正的步骤。
%%%雷达系统参数设置
f0 = 10e9;%载频
b1 = 100e6;%发送带宽
tp = 20e-6;%脉宽
vp = 100;%飞机飞行速度
fs = 250e6;%采样频率
fp = 1000;%发射脉冲
d = 1;%天线长度
c = 3e8;%电磁波传播速度
%%%目标区域参数设置
r = 10000;
x0 = 100;
y0 = 100;
r0 = [0 20 -2];
a0 = [0 0 -1];
nn = length(r0);
%%%衍生数据
j=sqrt(-1);
ppi=3.1415926;
piii= d/2;%方向分辨率
k= b1./tp;%调频速率
lanmada= c./f0;%波长
ka= -2(vp.^2)./(lanmada.r);%飞机飞行指数
sita= 0.886.lanmada/d;%雷达波束角
ls= rsita; %距离分辨率
%采样回波
ang=-ls/(2vp):1/fp:ls/(2vp);%飞机飞行时间
N=length(ang)
tao0= 2(r-x0/2)/c;%目标起始采集点t=tao0:1/fs:tao0+tp;%一次扫描时间
M=length(t);
rhui_matric= zeros(N,M);
h=exp(jppik*(t-tao0).^2);%匹配滤波信号符号
H=fft(h);
for i=1:N
rr=0;
for nj=1:nn
rp= vp.ang(i);%把飞机的运行位置放入一个矩阵
%rt_maxtric= sqrt((rp_maxtric-a0).2+(r+r0).2);
rt = sqrt((rp-a0(nj)).2+(r+r0(nj)).2);
tao= 2rt/c;
%%回波压缩 rr=rr+exp(-j4ppirt/lanmada).exp(jppik*(t-tao).^2).*(t>=tao&t<=tao+tp); end
R=fft(rr);
RC=R.*conj(H);
rc=ifft(RC);
rhui_matric(i,:)=rc;
end
为了保护版权,就没有给出完整代码啦!!!!!第一次写博客,哇哈哈哈。