合成孔径雷达

**

合成孔径雷达

**

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
为了保护版权，就没有给出完整代码啦！！！！！第一次写博客，哇哈哈哈。