使用Matlab仿真雷达信号处理系统.doc

毕业设计（论文）报告纸

摘　要
运用数字信号处理理论和Matlab 软件研究的脉冲压缩多普勒雷达的信号处理仿真问题,提出了一个仿真模型,该模型能够仿真雷达信号、系统噪声与杂波的产生和脉冲压缩多普勒雷达系统中信号的动态处理过程,最后结合MIMO雷达信号特点 ,显示了使用Matlab 仿真雷达信号处理系统方便快捷的特点。
关键词: MIMO 模糊图 脉冲压缩

Abstract
The use of digital signal processing theory and Matlab software research Doppler radar pulse compression signal processing simulation, a simulation model to simulation of radar signals, the system noise and clutter of the generation and pulse compression Doppler radar system Dynamic signal processing, the final combination of the characteristics of MIMO radar signal, indicating the use of Matlab simulation of the radar signal processing system characterized by convenient and efficient.
Key words: MIMO. Fuzzy Graph .pulse compressio

目 录
摘　要 2
ABSTRACT 3
第一章.雷达的基本原理 5
1.1雷达的原理 5
1.2雷达的用途 7
1.2.1双/多基地雷达 7
1.2.2相控阵雷达 7
1.2.3宽带／超宽带雷达 8
1.2.4合成孔径雷达 9
1.2.5毫米波雷达 9
1.2.6激光雷达 9
1.3相控阵雷达有多神 10
第二章．脉冲压缩基本原理 10
2.1引言 10
2.2线性调频（LFM）脉冲压缩雷达仿真 11
2.2.1雷达工作原理 11
2.2.2线性调频（LFM）信号 14
2.2.3LFM脉冲的匹配滤波 18
2.2.4 Matlab仿真结果 25
第三章 MIMO体制的原理 32
3.1引言 32
3.2 MIMO的几个概念 33
3.2.1 宽发窄收 33
3.2.2 综合脉冲和综合孔径 33
3.2.3 稀布阵 35
3.3 MIMO的性能及原理 36
3.3.1 正交编码的发射信号 36
3.3.2 波束综合 38
3.4 MIMO的组成和工作过程 43
第四章.结束语 44
参考文献 45
致 谢 46
附录：科技文献翻译 47
英文原文 47
中文翻译 57

第一章.雷达的基本原理
1.1雷达的原理
雷达(radar)原是“无线电探测与定位”的英文缩写。雷达的基本任务是探测感兴趣的目标，测定有关目标的距离、方问、速度等状态参数。雷达主要由天线、发射机、接收机（包括信号处理机）和显示器等部分组成。
雷达发射机产生足够的电磁能量，经过收发转换开关传送给天线。天线将这些电磁能量辐射至大气中，集中在某一个很窄的方向上形成波束，向前传播。电磁波遇到波束内的目标后，将沿着各个方向产生反射，其中的一部分电磁能量反射回雷达的方向，被雷达天线获取。天线获取的能量经过收发转换开关送到接收机，形成雷达的回波信号。由于在传播过程中电磁波会随着传播距离而衰减，雷达回波信号非常微弱，几乎被噪声所淹没。接收机放大微弱的回波信号，经过信号处理机处理，提取出包含在回波中的信息，送到显示器，显示出目标的距离、方向、速度等。
为了测定目标的距离，雷达准确测量从电磁波发射时刻到接收到回波时刻的延迟时间，这个延迟时间是电磁波从发射机到目标，再由目标返回雷达接收机的传播时间。根据电磁波的传播速度，可以确定目标的距离为：S=CT/2
其中S：目标距离
T：电磁波从雷达到目标的往返传播时间
C：光速
雷达测定目标的方向是利用天线的方向性来实现的。通过机械和电气上的组合作用，雷达把天线的小事指向雷达要探测的方向，一旦发现目标，雷达读出些时天线小事的指向角，就是目标的方向角。两坐标雷达只能测定目标的方位角，三坐标雷达可以测定方位角和俯仰角。
测定目标的运动速度是雷达的一个重要功能，—雷达测速利用了物理学中的多普勒原理．当目标和雷达之间存在着相对位置运动时，目标回波的频率就会发生改变，频率的改变量称为多普勒频移，用于确定目标的相对径向速度，通常，具有测速能力的雷达，例如脉冲多普勒雷达，要比一般雷达复杂得多。
雷达的战术指标主要包括作用距离、威力范围、测距分辨力与精度、测角分辨力与精度、测速分辨力与精度、系统机动性等。
其中，作用距离是指雷达刚好能够可靠发现目标的距离。它取决于雷达的发射功率与天线口径的乘积，并与目标本身反射雷达电磁波的能力（雷达散射截面积的大小）等因素有关。威力范围指由最大作用距离、最小作用距离、最大仰角、最小仰角及方位角范围确定的区域。
雷达的技术指标与参数很多，而且与雷达的体制有关，这里仅仅讨论那些与电子对抗关系密切的主要参数。 根据波形来区分，雷达主要分为脉冲雷达和连续波雷达两大类。当前常用的雷达大多数是脉冲雷达。常规脉冲雷达周期性地发射高频脉冲。相关的参数为脉冲重复周期（脉冲重复频率）、脉冲宽度以及载波频率。载波频率是在一个脉冲内信号的高频振荡频率，也称为雷达的工作频率。
雷达天线对电磁能量在方向上的聚集能力用波束宽度来描述，波束越窄，天线的方向性越好。但是在设计和制造过程中，雷达天线不可能把所有能量全部集中在理想的波束之内，在其它方向上在在着泄漏能量的问题。能量集中在主波束中，我们常常形象地把主波束称为主瓣，其它方向上由泄漏形成旁瓣。为了覆盖宽广的空间，需要通过天线的机械转动或电子控制，使雷达波束在探测区域内扫描。
概括起来，雷达的技术参数主要包括工作频率（波长）、脉冲重复频率、脉冲宽度、发射功率、天线波束宽度、天线波束扫描方式、接收机灵敏度等。技术参数是根据雷达的战术性能与指标要求来选择和设计的，因此它们的数值在某种程度上反映了雷达具有的功能。例如，为提高远距离发现目标能力，预警雷达采用比较低的工作频率和脉冲重复频率，而机载雷达则为减小体积、重量等目的，使用比较高的工作频率和脉冲重复频率。这说明，如果知道了雷达的技术参数，就可在一定程度上识别出雷达的种类。
1.2雷达的用途
雷达的用途广泛，种类繁多，分类的方法也非常复杂。通常可以按照雷达的用途分类，如预警雷达、搜索警戒雷达、无线电测高雷达、气象雷达、航管雷达、引导雷达、炮瞄雷达、雷达引信、战场监视雷达、机载截击雷达、导航雷达以及防撞和敌我识别雷达等。除了按用途分，还可以从工作体制对雷达进行区分。这里就对一些新体制的雷达进行简单的介绍。（军事观察·）
1.2.1双/多基地雷达
普通雷达的发射机和接收机安装在同一地点，而双/多基地雷达是将发射机和接收机分别安装在相距很远的两个或多个地点上，地点可以设在地面、空中平台或空间平台上。由于隐身飞行器外形的设计主要是不让入射的雷达波直接反射回雷达，这对于单基地雷达很有效。但入射的雷达波会朝各个方向反射，总有部分反射波会被双/多基地雷达中的一个接收机接收到。美国国防部从七十年代就开始研制、试验双/多基地雷达，较著名的“圣殿”计划就是专门为研究双基地雷达而制定的，已完成了接收机和发射机都安装在地面上、发射机安装在飞机上而接收机安装在地面上、发射机和接收机都安装在空中平台上的试验。俄罗斯防空部队已应用双基地雷达探测具有一定隐身能力的飞机。英国已于70年代末80年代初开始研制双基地雷达，主要用于预警系统。
1.2.2相控阵雷达
我们知道，蜻蜓的每只眼睛由许许多多个小眼组成，每个小眼都能成完整的像，这样就使得蜻蜓所看到的范围要比人眼大得多。与此类似，相控阵雷达的天线阵面也由许多个辐射单元和接收单元（称为阵元）组成，单元数目和雷达的功能有关，可以从几百个到几万个。这些单元有规则地排列在平面上，构成阵列天线。利用电磁波相干原理，通过计算机控制馈往各辐射单元电流的相位，就可以改变波束的方向进行扫描，故称为电扫描。辐射单元把接收到的回波信号送入主机，完成雷达对目标的搜索、跟踪和测量。每个天线单元除了有天线振子之外，还有移相器等必须的器件。不同的振子通过移相器可以被馈入不同的相位的电流，从而在空间辐射出不同方向性的波束。天线的单元数目越多，则波束在空间可能的方位就越多。这种雷达的工作基础是相位可控的阵列天线，“相控阵”由此得名。
相控阵雷达的优点
（1）波束指向灵活，能实现无惯性快速扫描，数据率高；（2）一个雷达可同时形成多个独立波束，分别实现搜索、识别、跟踪、制导、无源探测等多种功能； （3）目标容量大，可在空域内同时监视、跟踪数百个目标；（4）对复杂目标环境的适应能力强；（5）抗干扰性能好。全固态相控阵雷达的可靠性高，即使少量组件失效仍能正常工作。但相控阵雷达设备复杂、造价昂贵，且波束扫描范围有限，最大扫描角为90°～120°。当需要进行全方位监视时，需配置3～4个天线阵面。
相控阵雷达与机械扫描雷达相比，扫描更灵活、性能更可靠、抗干扰能力更强，能快速适应战场条件的变化。多功能相控阵雷达已广泛用于地面远程预警系统、机载和舰载防空系统、机载和舰载系统、炮位测量、靶场测量等。美国“爱国者”防空系统的AN／MPQ-53雷达、舰载“宙斯盾”指挥控制系统中的雷达、B-1B轰炸机上的APQ-164雷达、俄罗斯C-300防空武器系统的多功能雷达等都是典型的相控阵雷达。随着微电子技术的发展，固体有源相控阵雷达得到了广泛应用，是新一代的战术防空、监视、火控雷达。
1.2.3宽带／超宽带雷达
工作频带很宽的雷达称为宽带／超宽带雷达。隐身兵器通常对付工作在某一波段的雷达是有效的，而面对覆盖波段很宽的雷达就无能为力了，它很可能被超宽带雷达波中的某一频率的电磁波探测到。另一方面，超宽带雷达发射的脉冲极窄，具有相当高的距离分辨率，可探测到小目标。目前美国正在研制、试验超宽带雷达，已完成动目标显示技术的研究，将要进行雷达波形的试验。
1.2.4合成孔径雷达
合成孔径雷达通常安装在移动的空中或空间平台上，利用雷达与目标间的相对运动，将雷达在每个不同位置上接收到的目标回波信号进行相干处理，就相当于在空中安装了一个“大个”的雷达，这样小孔径天线就能获得大孔径天线的探测效果，具有很高的目标方位分辨率，再加上应用脉冲压缩技术又能获得很高的距离分辨率，因而能探测到隐身目标。合成孔径雷达在军事上和民用领域都有广泛应用，如战场侦察、火控、制导、导航、资源勘测、地图测绘、海洋监视、环境遥感等。美国的联合监视与目标攻击雷达系统飞机新安装了一部AN／APY3型X波段多功能合成孔径雷达，英、德、意联合研制的“旋风”攻击机正在试飞合成孔径雷达。
1.2.5毫米波雷达
工作在毫米波段的雷达称为毫米波雷达。它具有天线波束窄、分辩率高、频带宽、抗干扰能力强等特点，同时它工作在目前隐身技术所能对抗的波段之外，因此它能探测隐身目标。毫米波雷达还具有能力，特别适用于防空、地面作战和灵巧武器，已获得了各国的调试重视。例如，美国的“爱国者”防空导弹已安装了毫米波雷达导引头，目前正在研制更先进的毫米波导引头；俄罗斯已拥有连续波输出功率为10千瓦的毫米波雷达；英、法等国家的一些防空系统也都将采用毫米波雷达。
1.2.6激光雷达
工作在红外和可见光波段的雷达称为激光雷达。它由激光发射机、光学接收机、转台和信息处理系统等组成，激光器将电脉冲变成光脉冲发射出去，光接收机再把从目标反射回来的光脉冲还原成电脉冲，送到显示器。隐身兵器通常是针对微波雷达的，因此激光雷达很容易“看穿”隐身目标所玩的“把戏”；再加上激光雷达波束窄、定向性好、测量精度高、分辨率高，因而它能有效地探测隐身目标。激光雷达在军事上主要用于靶场测量、空间目标交会测量、目标精密跟踪和瞄准、目标成像识别、导航、精确制导、综合火控、直升机防撞、化学战剂监测、局部风场测量、水下目标探测等。美国国防部正在开发用于目标探测和识别的激光雷达技术，已进行了前视／下视激光雷达的试验，主要探测伪装树丛中的目标。法国和德国正在积极进行使用激光雷达探测和识别直升机的联合研究工作.
1.3相控阵雷达有多神
“宙斯盾”系统的核心就是SPY—1D相控阵雷达，特别是它出众的预警搜索能力和识别能力，仿佛给妄图“独立”的台湾新领导人一根救命稻草，一把梦幻的保护伞，而相控阵雷达又再一次走进国人的视线中。说到相控阵雷达或技术，大家可能很陌生，但如果说起去年美国军方关于中国如何监测其隐型战斗机的报道，大家可能就清楚了。用一大串电视接收天线来监视天空，经济又有效，这就是最原始、最基础的雷达，相控阵雷达。
第二章．脉冲压缩基本原理
2.1引言
现代雷达是一个十分复杂的工程系统 ,雷达系统的设计研究也是一项十分复杂的技术 ,人们越来越重视采用计算机仿真技术进行系统的分析和设计。利用计算机仿真技术进行雷达系统的建模与仿真 ,可以高效地完成系统的方案论证和性能评估 ,将雷达系统设计人员从繁重的设计工作中解脱出来 ,使雷达系统的设计更加方便、高效和优化 ,能够大大提高设计的可靠性,并可缩短设计周期 ,降低开发成本 ,所以它是当前和未来雷达与电子对抗领域研究中的一种重要手段。在雷达信号处理系统中 ,系统级仿真占有极其重要的地位。经过系统级仿真 ,能够确保产品在高层次上的设计正确性。利用系统仿真这一方法 ,可以找出系统各参数的最佳值来保证所设计的产品获得最佳性能。本文提出了一个脉冲压缩多普勒雷达信号处理系统的仿真模型 ,并利用 Matlab 提供的强大仿真平台对该信号处理系统进行了仿真。
2.2线性调频（LFM）脉冲压缩雷达仿真
2.2.1雷达工作原理
雷达是Radar（RAdio Detection And Ranging）的音译词，意为“无线电检测和测距”，即利用无线电波来检测目标并测定目标的位置，这也是雷达设备在最初阶段的功能。典型的雷达系统如图1.1，它主要由发射机，天线，接收机，数据处理，定时控制，显示等设备组成。利用雷达可以获知目标的有无，目标斜距，目标角位置，目标相对速度等。现代高分辨雷达扩展了原始雷达概念，使它具有对运动目标(飞机，导弹等)和区域目标(地面等)成像和识别的能力。雷达的应用越来越广泛。

图1.1：简单脉冲雷达系统框图
雷达发射机的任务是产生符合要求的雷达波形（Radar Waveform），然后经馈线和收发开关由发射天线辐射出去，遇到目标后，电磁波一部分反射，经接收天线和收发开关由接收机接收，对雷达回波信号做适当的处理就可以获知目标的相关信息。
假设理想点目标与雷达的相对距离为R，为了探测这个目标，雷达发射信号,电磁波以光速C向四周传播，经过时间后电磁波到达目标，照射到目标上的电磁波可写成：。电磁波与目标相互作用，一部分电磁波被目标散射，被反射的电磁波为，其中为目标雷达散射截面（Radar Cross Section ,简称RCS），反映目标对电磁波的散射能力。再经过时间后，被雷达接收天线接收的信号为。
如果将雷达天线和目标看作一个系统，便得到如图1.2的等效，而且这是一个LTI（线性时不变）系统。

图1.2：雷达等效于LTI系统
等效LTI系统的冲击响应可写成:
(1.1)
M表示目标的个数，为目标散射特性，是光速在雷达与目标之间往返一次的时间：
(1.2)
式中，为第i个目标与雷达的相对距离。
雷达发射信号经过该LTI系统，得输出信号(即雷达的回波信号)：
(1.3)
那么，怎样从雷达回波信号提取出表征目标特性的(表征相对距离)和(表征目标反射特性)呢？常用的方法是让通过雷达发射信号的匹配滤波器，如图1.3。

图1.3：雷达回波信号处理
的匹配滤波器为：
(1.4)
于是， （1.5）
对上式进行傅立叶变换：
(1.6)
如果选取合适的，使它的幅频特性为常数，那么1.6式可写为：
(1.7)
其傅立叶反变换为： (1.8)
中包含目标的特征信息和。从 中可以得到目标的个数M和每个目标相对雷达的距离： （1.9）

这也是线性调频（LFM）脉冲压缩雷达的工作原理。
2.2.2线性调频（LFM）信号
脉冲压缩雷达能同时提高雷达的作用距离和距离分辨率。这种体制采用宽脉冲发射以提高发射的平均功率，保证足够大的作用距离；而接受时采用相应的脉冲压缩算法获得窄脉冲，以提高距离分辨率，较好的解决雷达作用距离与距离分辨率之间的矛盾。
脉冲压缩雷达最常见的调制信号是线性调频（Linear Frequency Modulation）信号,接收时采用匹配滤波器（Matched Filter）压缩脉冲。
LFM信号(也称Chirp 信号)的数学表达式为：
(2.1)
式中为载波频率，为矩形信号，
（2.2）
，是调频斜率，于是，信号的瞬时频率为，如图2.1

图2.1 典型的chirp信号（a）up-chirp（b）down-chirp
将2.1式中的up-chirp信号重写为：
（2.3）
式中，
（2.4）
是信号s(t)的复包络。由傅立叶变换性质，S(t)与s(t)具有相同的幅频特性，只是中心频率不同而以，因此，Matlab仿真时，只需产生S(t)。以下Matlab程序产生2.4式的chirp信号，并作出其时域波形和幅频特性，如图2.2。

T=10e-6; %发射脉宽10us
B=30e6; %调频带宽30MHz
K=B/T; %频率调制斜率
Fs=2.5*B;Ts=1/Fs; %计算机仿真的采样频率和采样周期
N=T/Ts; %采样点数
t=linspace(-T/2,T/2,N);
St=exp(j*pi*K*t.^2); %产生线性调频信号
subplot(211)
plot(t*1e6,real(St));
xlabel('\fontsize{9}时间(us)');title('\fontsize{9}LFM脉冲的时域波形');
grid on;axis tight;
subplot(212)
freq=linspace(-Fs/2,Fs/2,N);
plot(freq*1e-6,fftshift(abs(fft(St))));
xlabel('\fontsize{9}频率(MHz)');title('\fontsize{9}LFM脉冲的幅频特性');
grid on;axis tight;

仿真结果显示：

图2.2：LFM信号的时域波形和幅频特性

2.2.3LFM脉冲的匹配滤波
信号的匹配滤波器的时域脉冲响应为：
（3.1）
是使滤波器物理可实现所附加的时延。理论分析时，可令＝0，重写3.1式，
（3.2）
将2.1式代入3.2式得:
(3.3 )

图3.1：LFM信号的匹配滤波
如图3.1,经过系统得输出信号，

当时,
(3.4)
当时,
(3.5)
合并3.4和3.5两式：
（3.6）
3.6式即为LFM脉冲信号经匹配滤波器得输出,它是一固定载频的信号。当时，包络近似为辛克（sinc）函数。
（3.7）

图3.2：匹配滤波的输出信号
如图3.2，当时，为其第一零点坐标；当时，，习惯上，将此时的脉冲宽度定义为压缩脉冲宽度。
(3.8)
LFM信号的压缩前脉冲宽度T和压缩后的脉冲宽度之比通常称为压缩比D，
（3.9）
3.9式表明，压缩比也就是LFM信号的时宽频宽积。
由2.1,3.3,3.6式，s(t),h(t),so(t)均为复信号形式，Matab仿真时，只需考虑它们的复包络
S(t),H(t),So(t)。以下Matlab程序段仿真了图3.1所示的过程，并将仿真结果和理论进行对照。
T=10e-6; %发射脉宽10us
B=30e6; %调频带宽30MHz
K=B/T; %频率调制斜率
Fs=10*B;Ts=1/Fs; %计算机仿真的采样频率和采样周期
N=T/Ts;
t=linspace(-T/2,T/2,N);
St=exp(j*pi*K*t.^2); %产生线性调频信号
Ht=exp(-j*pi*K*t.^2); %匹配滤波器单位冲激响应
Sot=conv(St,Ht); %匹配滤波输出
figure(1)
L=2*N-1;
t1=linspace(-T,T,L);
Z=abs(Sot);Z=Z/max(Z); %仿真信号对数归一化
Z=20*log10(Z+1e-6);
Z1=abs(sinc(B.*t1)); %产生理伦输出信号并归一化
Z1=20*log10(Z1+1e-6);
t1=t1*B; %时间轴与1/B归一化
plot(t1,Z,t1,Z1,'r.');
axis([-15,15,-50,inf]);grid on;
legend('\fontsize{9}仿真结果','理论结果');
xlabel('\fontsize{9}时间\times\itB (s)');
ylabel('\fontsize{9}相对幅度dB');
title('\fontsize{9}脉冲压缩后的回波(归一化后)');
figure(2) %仿真结果局部放大
N0=3*Fs/B;
t2=-N0*Ts:Ts:N0*Ts;
t2=B*t2;
plot(t2,Z(N-N0:N+N0),t2,Z1(N-N0:N+N0),'r.');
axis([-inf,inf,-50,inf]);grid on;
set(gca,'Ytick',[-13.4,-4,0],'Xtick',[-3,-2,-1,-0.5,0,0.5,1,2,3]);
legend('\fontsize{9}仿真结果','理论结果');
xlabel('\fontsize{9}时间\times\itB (s)');
ylabel('\fontsize{9}相对幅度dB');
title('\fontsize{9}脉冲压缩后的回波(归一化后局部图)');
仿真结果如图3.3，3 .4。

图3.3：脉冲压缩后的回波

图3.4：脉冲压缩后的回波(局部图)

图3.4中，时间轴进行了归一化，（）。图中反映出理论与仿真结果吻合良好。第一零点出现在（即）处，此时相对幅度-13.4dB。压缩后的脉冲宽度近似为（），此时相对幅度-4dB,这理论分析（图3.2）一致。

上面只是对各个信号复包络的仿真，实际雷达系统中，LFM脉冲的处理过程如图2.4。

图3.5： LFM信号的接收处理过程
雷达回波信号（1.4式）经过正交解调后，得到基带信号，再经过匹配滤波脉冲压缩后就可以作出判决。正交解调原理如图3.6，雷达回波信号经正交解调后得两路相互正交的信号I(t)和Q(t)。一种数字方法处理的的匹配滤波原理如图3.7。

图3.6：正交解调原理

图3.7：一种脉冲压缩雷达的数字处理方式
2.2.4 Matlab仿真结果
（1）任务：对以下雷达系统仿真。
雷达发射信号参数：
幅度：1.0
信号波形：线性调频信号
频带宽度：30兆赫兹（30MHz）
脉冲宽度：10微妙（20us）
中心频率：1GHz（109Hz）
雷达接收方式：
正交解调接收
距离门：10Km~15Km
目标：
Tar1：10.5Km
Tar2：11Km
Tar3：12Km
Tar4：12Km＋5m
Tar5：13Km
Tar6：13Km＋2m
（2）系统模型：
结合以上分析，用Matlab仿真雷达发射信号，回波信号，和压缩后的信号的复包络特性，其载频不予考虑（实际中需加调制和正交解调环节），仿真信号与系统模型如图4.1。

图4.1：雷达仿真等效信号与系统模型
function LFM_radar(T,B,Rmin,Rmax,R,RCS)
if nargin==0
T=10e-6; %发射脉宽
B=30e6; %调频带宽
Rmin=10000;Rmax=15000; %距离门
R=[10500,11000,12000,12005,13000,13002]; %目标位置向量，相对于雷达
RCS=[1 1 1 1 1 1]; %目标相对有效反射面
end
%参数设置
R=R-Rmin; %目标位置向量，相对距离门下限Rmin
C=3e8; %光速
K=B/T; %调频斜率
Rrec=Rmax-Rmin; %接收门宽/m
Trec=2*Rrec/C; %接收门宽/s
Fs=5*B;Ts=1/Fs; %计算机仿真的采样频率和周期
N0=ceil(T/Ts); %发射脉宽对应的采样点数
N=ceil(Trec/Ts); %接收窗对应的采样点数
%产生高采样率chirp信号
Rate=10;
t1=linspace(0,T,Rate*N0);
Chirp=exp(j*pi*K*t1.^2);
%产生回波信号
M=length(R);
Srt(1:N)=0; %初始化接收信号
for i=1:1:M
temp(1:N)=0;
tau=2*R(i)/C;
K0=ceil(tau/Ts) %第i个目标回波的相对位置
M0=fix((Ts*K0-tau)*Rate/Ts)
temp(K0:K0+N0-1)=RCS(i)*Chirp(M0+1:Rate:M0+Rate*(N0-1)+1);
Srt=Srt+temp; %将所有的目标回波叠加
end
%作出相关图形
subplot(211)
t=linspace(0,Trec,N);
plot(t*10^6,real(Srt));
xlabel('\fontsize{9}时间(us)');title('\fontsize{9}脉冲压缩前的回波');
subplot(212)
t1=linspace(-T/2,T/2,N0);
h=exp(-j*pi*K*t1.^2); %匹配滤波器单位冲激响应
Sot=conv(Srt,h); %匹配滤波
Z=abs(Sot);Z=Z/max(Z);
Z=20*log10(Z+1e-6);
dist=linspace(0,Rrec,N);
plot(dist+Rmin,Z(N0:N+N0-1));
axis([Rmin,Rmax,-60,0]);
xlabel('\fontsize{9}目标相对了雷达距离(m)');title('\fontsize{9}脉冲压缩后的回波');

函数LFM_radar的参数意义如下：
T：chirp信号的持续脉宽；
B：chirp信号的调频带宽；
Rmin：观测目标距雷达的最近位置；
Rmax：观测目标距雷达的最远位置；
R：一维数组，数组值表示每个目标相对雷达的斜距；
RCS：一维数组，数组值表示每个目标的雷达散射截面。

在Matlab指令窗中键入：

LFM_radar(10e-6,30e6,10000,15000,[10500,11000,12000,12008,13000,13005],[1,1,1,1,1,1])

得到的仿真结果如下：

图4.2：仿真结果
改变目标的相对距离，便可得到如图4.3的组图，当T=10us，B=30MHz时，雷达的距离分辨率为：
(4.1)

(a) 图中，两目标相距2m，显然区分不开；
(b) 图中，两目标相距5m，实际上是两目标的输出sinc包络叠加，可以看到他们的副瓣相互抵消；
(c) ，(d)两图中，两目标距离大于雷达的距离分辨率，因而能分辨出，而且，随着目标距离越大，雷达越容易区分。

图4.3：雷达分辨率仿真

第三章 MIMO体制的原理
3.1引言
本章将对MIMO雷达（MIMO)的原理性质加以介绍。给出MIMO雷达的框图，并介绍稀布阵、全向反射发射、综合脉冲和综合孔径、波束形成等基本理论，简要介绍MIMO处理的基本思想，说明MIMO是一种具有优越性能的新体制雷达。
3.2 MIMO的几个概念
3.2.1 宽发窄收
为了确定目标的方向，雷达天线定向发射和接收信号，发射信号越是集中在一个角落里，越是有利于增加测量距离。发射信号对于作用距离的影响，不仅仅是发射功率，更主要的是它的能量，例如脉冲雷达，脉冲功率并不大，但是宽度长，通过脉冲压缩得到相干积累而使得回波信号峰值提高。
监视雷达的威力范围不仅指作用距离，同时也包括覆盖的空域，常规雷达的工作方式使得它的波束空间角小，驻留时间短。为了保持一定的角分辨力，波束空间角不要太窄，而扫描时间又不允许改变，解决的方法之一是采取“宽发窄收”，接收波束为窄波束，以保持一定的角分辨力，同时将发射波束加宽。在宽发射波束覆盖的区域里，有多个接收波束同时接收，显然接收波束的驻留时间的延长与发射波束的加宽成正比，不过向一定方向的发射功率与发射波束的加宽成反比。在驻留时间内接收的回波能量是相同的，只要能有效的进行相干积累（相当于匹配滤波），两种方式的作用距离相同。但宽发窄收还能进一步延长驻留时间，将可以获得更大的作用距离。虽然采用“宽发”的方式，能量在空间分散了，不能有效地照射目标而不利于检测，实际上搜索雷达地任务和跟踪雷达不同，它要对规定地空域，进行搜索，它必须将能量辐射到所有需要搜索地空间，而何时出现目标是不可预知的，这两种方式并无差别。更本