DSP多普勒雷达测速测距(精)

DSP实验课大作业设计

一实验目的

在DSP上实现线性调频信号的脉冲压缩、动目标显示（MTI）和动目标检测(MTD)，并将结果与MATLAB上的结果进行误差仿真。

二实验内容

2.1 MATLAB仿真

设定带宽、脉宽、采样率、脉冲重复频率，用MATLAB产生16个脉冲的LFM，每个脉冲有4个目标（静止，低速，高速），依次做

2.1.1 脉压

2.1.2 相邻2脉冲做MTI，产生15个脉冲

2.1.3 16个脉冲到齐后，做MTD，输出16个多普勒通道

2.2 DSP实现

将MATLAB产生的信号，在visual dsp中做脉压，MTI、MTD，并将结果与MATLAB作比较。

三实验原理

3.1 脉冲压缩原理及线性调频信号

雷达中的显著矛盾是：雷达作用距离和距离分辨率之间的矛盾以及距离分辨率和速度分辨率之间的矛盾。雷达的距离分辨率取决于信号带宽。在普通脉冲雷达中，雷达信号的时宽带宽积为一常量（约为1），因此不能兼顾距离分辨率和速度分辨力两项指标。脉冲压缩（PC）采用宽脉冲发射以提高发射的平均功率，保证足够的最大作用距离，而在接收时则采用相应的脉冲压缩法获得窄脉冲，以提高距离分辨率，因而能较好地解决作用距离和分辨能力之间的矛盾。

一个理想的脉冲压缩系统，应该是一个匹配滤波系统。它要求发射信号具有非线性的相位谱，并使其包络接近矩形；要求压缩网络的频率特性（包括幅频特性和相频特性）与发射脉冲信号频谱（包括幅度谱和相位谱）实现完全的匹配。

脉冲压缩按信号的调制规律（调频或调相）分类，可分为以下四种：

（1）线性调频脉冲压缩

（2）非线性调频脉冲压缩

（3）相位编码脉冲压缩

（4）时间频率编码脉冲压缩

本实验采用的是线性调频脉冲压缩。

线性调频信号是指频率随时间的变化而线性改变的信号。线性调频可以同时保留连续信号和脉冲的特性，并且可以获得较大的压缩比，有着良好的距离分辨率和径向速度分辨率，所以将线性调频信号作为雷达系统中一种常用的脉冲压缩信号。

接收机输入端的回波信号是经过调制的宽脉冲，所以在接收机中应该设置一个与发射信号频率匹配的滤波器，使回波信号变成窄脉冲，同时实现了宽脉冲的能量和窄脉冲的分辨能力。解决了雷达发射能量及分辨率之间的矛盾。

匹配滤波器是指输出信噪比最大准则下的最佳线性滤波器。根据匹配理论，匹配滤波器的传输特性：

*S其中，K为幅度归一化常数，为信号的复共轭。传输特性还可用它

。 h(t)的冲激响应来表示（时域表示）：

3.2 MTI（动目标显示）原理

在雷达终端的显示器，当杂波与运动的目标在显示器上同时出现时，观察将会变得困难，特别是在强杂波时，弱的目标信号将会完全被淹没。由于目标回波的多普勒频移远大于杂波的多普勒频移，从而可在频域上区分目标与杂波。MTI的设计主要是合理的选择滤波器，动目标显示滤波器（MTI）利用运动目标回波和杂波在频谱上的区别，有效地抑制杂波而提取信号。由下图可以看出回波信号（含杂波与目标）：

图3.1 回波信号

杂波主要集中在0频附近，在频率为0处，所设计的滤波器频率响应应有凹口。以使地物杂波在通过MTI滤波器后将受到很大的抑制，对于每一个脉冲，其回波信号中的杂波分量都基本相同，而运动目标产生多普勒频移。对连续脉冲对的回波进行相减，就可以完全对消杂波分量，而目标信号由于其相位的改变，通常不会完全被对消掉。

常用的MTI滤波器为一次相消器（二脉冲对消）。在同一距离上等效一个有效的采样间隔的离散时间序列x(n)。其时域方程为：，传输函数

为：，在Z=1处有单零点，频率响应为：

。其在零频有一凹口，如图所示：

图3.3 多普勒滤波器频响

此滤波器能够用来抑制噪声。但同时把静态目标也给对消掉，也使得运动目标的谱分量有部分的衰减，因此用MTI一次相消器检测不出静止目标。

3.3 MTD（动目标检测）原理

MTD与MTI同属雷达信号处理的频域处理范畴，但一般意义上来说，MTD是MTI的改进或更有效的频域处理技术。广义上地讲，MTD处理又是多普勒处理的一种特殊形式与传统的MTI相比，MTD改善了滤波器的频率特性，使之更接近于最佳匹配线性滤波，以提高改善因子；能够检测强地物杂波的低速目标甚至切向飞行的大目标。MTD的核心是线性MTI加窄带

多普勒滤波器组。根据最佳线性滤波理论，在杂波背景下检测运动目标回波，除了杂波抑制滤波器外，还应串接有对脉冲串信号匹配的滤波器，实际工作中，多普勒频移fd不能预知，需要采用一组相邻且部分重叠的滤波器组覆盖整个多普勒频率范围，这就是窄带多普勒滤波器组。输入N个脉冲的横向滤波器组的实现是由N个输出的横向滤波器（N个脉冲和N-1根迟延线）经过各脉冲不同的加权并求和后形成的。其实现方法下图所示：

图3.4 横向滤波器图3.5 N=16时滤波器频响

该滤波器的频率覆盖范围为0到fr，fr为雷达工作重复频率。气象杂波(如云雨等)和箔条杂波受气流和风力的影响，会相对雷达而动，其频谱中心不是固定在0频附近，而是在某个频率区间内变化的，抑制此类杂波用普通的MTI滤波器是不行的，而MTD滤波器则可以抑制此类杂波。

四实验步骤

以下是该实验中设定的几个参数

BandWidth=2.0e6----------------------------------------------带宽

TimeWidth=42.0e-6--------------------------------------------脉宽

Fs=2.0e6---------------------------------------------------采样率

PRT=240e-6-------------------------------------------脉冲重复周期

TargetDistance=[3000 8025 8025 25600] ------------------目标距离

TargetVelocity=[50 -100 0 283]---------------------------目标速度

%其中X=0; Y=8; Z=3;对应学号1302121083的后三位

假设接收到的回波数是16个，噪声为高斯随机噪声。

4.1 在MATLAB中产生线性调频信号。

4.2 根据目标距离得出目标回波在时域的延迟量，根据目标速度得出多普勒相

移，从而在MATLAB中产生4个目标的16个回波串（接收到的回波含噪声）。

4.3 由匹配滤波理论产生对应于目标回波的滤波系数（脉压系数），在时域中做

线性卷积，实现匹配滤波（时域脉压）；在频域中做回波信号和脉压系数的FFT，点乘后作逆FFT，实现频域脉压。两者进行比较，讨论其差别。

4.4 对16个去暂态点后的脉冲串按接收顺序进行排列，用一次相消器（一种滤波方式）实现MTI。

4.5 做16通道的FFT，实现MTD。

4.6 在DSP中对MATLAB产生的回波数据和脉压系数进行处理，实现频域脉

压。导入DSP的回波数据为时域数据，而脉压系数为频域数据。将导入DSP的时域回波数据进行一次FFT变换到频域，然后将其与频域脉压系数进行点积，得到频域脉压结果。对该结果再做一次逆FFT，将频域转换成时域。在这一步中需要调用库函数fft_flp32.asm。该子程序可实现8192