Advertisement

mimo雷达信号处理_雷达学术入门脉冲雷达信号处理概述

阅读量:
复制代码

Reviewed by :@甜草莓 @Robert Zhou;

前置知识:概率论与统计学。面向人群:本科生、研究生/信号处理博士。 ``

编者:对于信号处理来说,雷达和通信一直是一体两面,从MIMO通信到MIMO雷达,从OFDM通信到Multicarrier雷达,很多通信和雷达领域的前沿都在事实上的相互影响。本篇文章将会介绍Fan Liu博士即将在IEEE Communications Letters发表的系列约稿里的第一篇,雷达基本原理。在此感谢Fan Liu博士不辞劳苦的中文翻译工作,并且愿意在公众号首发学术论文,来做雷达领域科普的精神和努力(:

本文旨在面向通信背景的读者简明扼要地概括脉冲雷达信号处理中的基本问题和方法

由于篇幅所限,本文将不会对相关问题做深入讨论。有兴趣的读者可以进一步阅读本文后面的参考文献。
4477993641f1c0aa154ac8c3475a58a7.png脉冲雷达基本模型

1. 基本原理

如上图所示,脉冲式雷达的基本工作流程,我们可以做如下阐述: 雷达首先发射一个探测脉冲,脉冲信号到达目标后被目标反射至雷达端。这一回波中包含了目标距离、速度、角度等参数信息。雷达的任务是对回波进行处理提取出目标信息,从而对目标进行定位或跟踪。 由于雷达周期性地发射信号并接收回波,这一发射-接收循环通常被称为脉冲重复周期(Pulse Repetition Interval, PRI) 。类似地,一秒钟发射的脉冲个数被定义为脉冲重复频率(Pulse Repetition Frequency, PRF)
a8956f5090fcb6e5eeadb11f93182205.png

配图仅为示意

那么延续前文思路,我们如果给出简洁的雷达建模? 如果我们考虑一个具有 417ed096-0736-eb11-8da9-e4434bdf6706.svg 发射天线和 437ed096-0736-eb11-8da9-e4434bdf6706.svg 接收天线的雷达,且雷达正在探测位于远场(Far Field)的一个点状目标。该目标相对于雷达的距离、速度以及角度分别表示为 457ed096-0736-eb11-8da9-e4434bdf6706.svg 。进一步地,将雷达发射的脉冲信号表示为 497ed096-0736-eb11-8da9-e4434bdf6706.svg ,则回波信号可以表示为 4c7ed096-0736-eb11-8da9-e4434bdf6706.svg

  • 527ed096-0736-eb11-8da9-e4434bdf6706.svg 包含了双程路径损耗(Round-trip Pathloss)及目标的雷达散射截面积(Radar Cross-Section, RCS),

  • 547ed096-0736-eb11-8da9-e4434bdf6706.svg 表示多普勒频率(577ed096-0736-eb11-8da9-e4434bdf6706.svg5c7ed096-0736-eb11-8da9-e4434bdf6706.svg 分别表示信号载频和光速),

  • 617ed096-0736-eb11-8da9-e4434bdf6706.svg627ed096-0736-eb11-8da9-e4434bdf6706.svg 分别是发射和接收方向矢量(Steering Vector),由发射以及接收天线阵列几何来决定,

  • 657ed096-0736-eb11-8da9-e4434bdf6706.svg 代表双程时延(Round-trip Delay),

  • 697ed096-0736-eb11-8da9-e4434bdf6706.svg 表示加性高斯白噪声,其方差为 6c7ed096-0736-eb11-8da9-e4434bdf6706.svg

注意:除了噪声以外,雷达还会收到来自其他方位的干扰信号。这些干扰可以来自其他信号源,也可以是不感兴趣的目标或者障碍物反射至雷达的回波。一般将后者称为杂波(Clutter)。通常需要对杂波进行抑制,以免对感兴趣的目标造成干扰。 受限于篇幅,我们将不在此文中讨论杂波抑制问题。

2.目标检测问题

通常我们认为,雷达信号处理会关心两种基本问题:检测(Detection)和估计(Estimation)。

其中,最简单的检测问题一般只考虑目标是否存在,通常可以建模为如下二元假设检验(Binary Hypothesis Testing):
6f7ed096-0736-eb11-8da9-e4434bdf6706.svg
717ed096-0736-eb11-8da9-e4434bdf6706.svg 假设(Null Hypothesis)代表雷达仅仅收到噪声的情况,727ed096-0736-eb11-8da9-e4434bdf6706.svg 假设则代表雷达同时收到了回波和噪声。

那么如何求解呢?

需要设计一个检测器(Detector)来对 737ed096-0736-eb11-8da9-e4434bdf6706.svg 进行检验。检测器 777ed096-0736-eb11-8da9-e4434bdf6706.svg 表征一种映射关系,将信号737ed096-0736-eb11-8da9-e4434bdf6706.svg映射为一个正实数,再将这一正实数与一个预先给定的门限 7e7ed096-0736-eb11-8da9-e4434bdf6706.svg 进行比较,从而决定选取717ed096-0736-eb11-8da9-e4434bdf6706.svg 还是727ed096-0736-eb11-8da9-e4434bdf6706.svg 假设。

这一过程可以表示为
857ed096-0736-eb11-8da9-e4434bdf6706.svg

上式中,如果 887ed096-0736-eb11-8da9-e4434bdf6706.svg则选取727ed096-0736-eb11-8da9-e4434bdf6706.svg ,反之选取717ed096-0736-eb11-8da9-e4434bdf6706.svg。根据场景和雷达所具有的不同先验信息,我们可以设计多种检测器来对回波信号进行检验。

常用的检测器包括似然比检验(Likelihood Ratio Test, LRT),广义似然比检验(Generalized Likelihood Ratio Test, GLRT),Rao检验(Rao Test)和Wald检验(Wald Test)等。

此外,为衡量目标检测的好坏,人们提出了多种评价指标来描述检测器的性能,最为常用的是检测概率(Detection Probability)和虚警概率(False-Alarm Probability)。

检测概率定义为目标存在且雷达判决727ed096-0736-eb11-8da9-e4434bdf6706.svg 假设为真的概率(即目标存在,且雷达也判断目标存在),虚警概率则定义为目标不存在且雷达判决727ed096-0736-eb11-8da9-e4434bdf6706.svg 假设为真的概率(即目标不存在,但雷达判断目标存在)。这两种概率可以用公式表示为:
917ed096-0736-eb11-8da9-e4434bdf6706.svg

除以上概率指标外,通常还关心所谓的漏警概率(False-Dismissal Probability) ,即存在目标但雷达判断目标不存在的概率,可记为 947ed096-0736-eb11-8da9-e4434bdf6706.svg

注意相比于漏警,虚警对雷达造成的伤害更大。因为一旦判断目标存在,雷达就要动用硬件和信号处理资源来对目标进行进一步探测与跟踪,虚警情形下,这将造成雷达资源的严重浪费。

有鉴于此,目标检测器的设计通常需要遵循所谓的奈曼-皮尔逊准则(Neyman-Pearson Criterion)。简言之,即在给定最小可容忍恒虚警概率的条件下最大化检测概率。

3. 参数估计问题

一旦确定目标存在,雷达就需要处理被噪声污染过的回波信号737ed096-0736-eb11-8da9-e4434bdf6706.svg,从而对目标的参数进行估计。如上文所述,对于点目标模型,感兴趣的参数通常包括距离,速度和角度,这些参数需要通过设计一个估计器(Estimator)来进行估计。

与检测器类似,估计器同样可以看做一种特殊的映射,即,将回波信号737ed096-0736-eb11-8da9-e4434bdf6706.svg从信号空间(Signal Space)映射到参数空间(Parameter Space)。这一过程可以表示为:
997ed096-0736-eb11-8da9-e4434bdf6706.svg

其中 9c7ed096-0736-eb11-8da9-e4434bdf6706.svg 是真实距离、速度和角度的估计值。

在实际情况中,一般需要针对这些参数设计独立的估计器,我们分别对此进行阐述。

  • 对于距离估计,一般采用所谓脉冲压缩方法来估计信号的时延 9d7ed096-0736-eb11-8da9-e4434bdf6706.svg ,然后算出对应的距离,其实质是利用延迟的发射信号副本对回波进行匹配滤波,匹配滤波器输出响应最大时对应的延迟即为目标的时延。这部分一般称为快时间(Fast-Time)信号处理,是对单个雷达脉冲内的信号样本(又称为“快拍”(Snapshot))的处理。

  • 对于速度估计,由于单个脉冲持续时间较短,其多普勒相移难以识别,通常需要接收多个脉冲的回波确保多普勒相移积累得足够大,再通过快速傅里叶变换(Fast Fourier Transform, FFT)来估计多普勒频率。这部分称为慢时间(Slow-Time)信号处理,即以一个脉冲对应一个慢时间单位的脉冲间信号处理。

  • 对于角度估计,可以采用经典的基于子空间的到达角估计算法,例如MUSIC或者ESPIRIT。

估计器的性能通常可以用均方误差(Mean Squared Error, MSE)进行描述。以角度估计为例,其均方误差定义为
a17ed096-0736-eb11-8da9-e4434bdf6706.svg

然而,参数估计的MSE一般很难求出闭式解,这是因为真值 a37ed096-0736-eb11-8da9-e4434bdf6706.svg 和估计值 a67ed096-0736-eb11-8da9-e4434bdf6706.svg 的分布在大多数情况下都难以得到。

作为一种替代手段,经常考虑的一个经典的性能指标是所谓的克拉美-罗下界(Cramér–Rao Lower Bound, CRLB)。简而言之,克拉美-罗界是所有无偏估计器的方差的下界(无偏估计定义为估计值的数学期望为真值的估计器,因此其方差等于MSE)。

这一下界指出,无偏估计器的估计方差应至少与Fisher信息(Fisher Information)的逆相当,即
a87ed096-0736-eb11-8da9-e4434bdf6706.svg

上式中,分母定义为关于 a37ed096-0736-eb11-8da9-e4434bdf6706.svg 的Fisher信息, ae7ed096-0736-eb11-8da9-e4434bdf6706.svg 是关于 a37ed096-0736-eb11-8da9-e4434bdf6706.svg 的似然函数(Likelihood Function)。

4. 总结

本文对脉冲式雷达的基本原理进行了简介,并重点介绍了雷达信号处理的两个基本问题:目标检测与参数估计。

由于篇幅所限,本文所考虑的只是最基本的点目标模型,对于雷达系统中的其他重要概念,如杂波干扰、扩展目标(Extended Target)、波形设计等均未作介绍。

感兴趣的读者请参考雷达和统计信号处理等教材,并欢迎与我单独交流。

参考文献

[1] M. A. Richards, J. Scheer, and W. A. Holm, Principles of Modern Radar: Basic Principles. New York, NY, USA: Scitech, 2010.[2] S. M. Kay, Fundamentals of Statistical Signal Processing. Englewood Cliffs, NJ, USA: Prentice Hall, 1998.

作者:刘凡,男,1992年出生,工学博士。现为英国伦敦大学学院电子与电气工程系博士后,“玛丽·居里”研究员。研究方向包括雷达通信一体化,车联网通信,5G毫米波通信和大规模MIMO等。先后于2013和2018年获得北京理工大学学士与博士学位,博士导师为王小谟院士。曾于2016-2018年赴伦敦大学学院访问研究。截止目前已写作学术论文40余篇,专利和专著章节各1篇。曾获得2019年中国电子学会优秀博士学位论文以及2018年欧盟“玛丽·居里”学者等多项奖励及荣誉。目前担任2020年IEEE ICC Workshop on Communication and Radar Spectrum Sharing论坛共同主席 ,以及多种IEEE期刊及国际学术会议的审稿人。曾于2017-2018被评为IEEE TWC, IEEE TCOM和IEEE Comms Letters的模范审稿人。

全部评论 (0)

还没有任何评论哟~