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.svg 和 5c7ed096-0736-eb11-8da9-e4434bdf6706.svg 分别表示信号载频和光速),

  •  617ed096-0736-eb11-8da9-e4434bdf6706.svg 和 627ed096-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的模范审稿人。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/489131.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

人类历史十大经典物理实验,有些实验就在你身边为何你没注意到?

来源:数学职业家实验是科学研究的基本方法之一,也是人类认识自然现象、自然性质、自然规律的途径。现代实验技术的发展,不断地揭示和发现各种新的物理现象,日益加深人们对客观世界规律的正确认识,从而推动物理学的向前…

一、linux搭建jenkins+github详细步骤

事情缘由: 现在在做的主要工作是通过jenkinspostman实现api的自动化测试,想要达到的效果是,api自动化测试定时跑脚本的同时,github有新的代码提交,jenkins会自动检测部署新提交的代码,并且再次将api自动化脚…

matlab插值与拟合例题_MATLAB中数据插值和数据拟合的用法

一、数据插值:插值是在一组已知数据点的范围内添加新数据点的技术。可以使用插值来填充缺失的数据、对现有数据进行平滑处理以及进行预测等。MATLAB 中的插值技术可分为适用于网格上的数据点和散点数据点。从数学上来说,数据插值是一种函数逼近的方法。数…

python中with的用法_python 中 with的用法

with 语句使用于对资源进行访问的场合,确保不管使用过程是否发生异常都会执行必要的"清理"操作,释放资源,比如文件使用后自动关闭/线程中锁的自动获取和释放with open("1.txt") as file:data file.read()使用with,能够减少冗长,还能自动处理上…

半导体行业深度报告:从应用到行业的全面复苏

来源:国金证券一、2020-2021年全球半导体市场投资展望多种因素导致全球半导体市场于 2019 年同比下跌近 13%到 4,102 亿美元,而存储器行业同比下跌超过 30%,逻辑半导体同比下跌近 2%。存储器市场占全球半导体市场达到近三年低点的 27%。但受惠…

ES6精华: 解构运算符 扩展运算符 剩余运算符

http://www.tuicool.com/articles/26bAzmm转载于:https://www.cnblogs.com/xuezizhenchengxuyuan/p/6378382.html

mysql 8 修改root密码忘记_忘记mysql8或者mariadb5及以上 的root密码如何更改

实验环境操作系统:lsb_release -aNo LSB modules are available.Distributor ID:UbuntuDescription:Ubuntu 20.04.1 LTSRelease:20.04Codename:focalmysql数据库版本:mysql基本信息具体步骤1. 停止mysql (sudo service mysqld stop 或者 kill -TERM mysq…

配置tomcat_Tomcat 配置必备的 10 个小技巧,让你轻松玩转Tomcat

现在开发Java Web应用,建立和部署Web内容是一件很简单的工作。使用Jakarta Tomcat作为Servlet和JSP容器的人已经遍及全世界。Tomcat具有免费、跨平台等诸多特性,并且更新得很快,现在非常的流行。你所需要做的就是:按照你的需求配置…

卫星移动通信现状与未来发展

来源:电子万花筒摘要:卫星移动通信系统具有覆盖范围广,对地面情况不敏感等优势,已经成为地面移动通信领域重要的组成部分,尤其是在空中、海洋、荒漠戈壁等地面无线网络难以覆盖的地方。随着科学技术的不断进步&#xf…

【SpringCloud】Netflix源码解析之Ribbon:负载均衡策略的定义和实现

Ribbon负载均衡策略定义 IRule其实就只做了一件事情Server choose(Object key),可以看到这个功能是在LB中定义(要求)的,LB把这个功能委托给IRule来实现。不同的IRule可以向LB提供不同的负载均衡算法。 public interface IRule{ pu…

mysql load data infile 重写_mysql load data infile 命令的数据导入

文章介绍了mysql load data infile 命令的数据导入,该方式比直接的insert的效率要高,按照官方的说法是要比insert语句快上20倍,有需要的朋友可参考。使用方式如下:代码如下复制代码mysql>load data local infile "D:/ab.txt" in…

DARPA发布战略框架文件旨在 为美国国家安全创建突破性新技术能力

来源:ARPA网站/图片来自互联网12月4日,DARPA公开了于今年8月份完成的2019年度战略框架文件《面向国家安全创建技术突破和新能力》,阐述了该机构面临的威胁环境,并提出了未来几年的发展计划。1应关注的重点该框架文件指出&#xff…

python数学公式pdf文件的转换_python转换文件 多种文件转换为pdf

1.office文件这里用的是win32com, 需要注意已经安装的python版本是32位还是64位。安装后导入from win32com.client import Dispatch, constants, gencache, DispatchEx(1) word转PDFgencache.EnsureModule({00020905-0000-0000-C000-000000000046}, 0, 8, 4)w DispatchEx(&quo…

rust矿洞绳子怎么爬下_rust矿洞绳子怎么爬下_打工小伙爬冰救人,每动一下都能听到冰面碎裂声...

近日,在辽宁沈阳鲁迅公园,一男子掉进湖中冰窟窿,露出脑袋不停地呼救。生死时刻,一位正在附近找零活的小伙,不顾冰面碎裂发出的咔咔声,一步步爬冰救人,最终在岸上众人帮助下,成功救出…

Python爬虫入门四之Urllib库的高级用法

1.设置Headers 有些网站不会同意程序直接用上面的方式进行访问,如果识别有问题,那么站点根本不会响应,所以为了完全模拟浏览器的工作,我们需要设置一些Headers 的属性。 首先,打开我们的浏览器,调试浏览器F…

Magic Leap 价值 20 多亿美元的 AR 幻梦,现在醒了

来源:传感器技术那个用特效来吹出 AR 神话的 Magic Leap,终于有了新动静。 当地时间 12 月 10 日,增强现实(AR)领域的“知名”公司 Magic Leap 宣布,将 Magic Leap One Creator Edition 改名为 Magic Leap …

mysql约束_不是吧,阿Sir,MySQL约束你竟然还不懂!

以前写的太乱了,翻出来重新整理下系列目录:MySQL入门,问题不大【增删改查极速上手】(一) 引入约束(1) 约束出现在哪里? 想要讲解约束,就要知道约束用在哪里,用来干嘛?SQL 语言通过定义一个关系所…

dc持久内存与mysql_Calypso Systems推出测试软件和服务器测试傲腾数据中心级持久内存...

Calypso Systems测试Optane DCPMM(DC Persistent Memory Module 数据中心级持久性内存模块)的性能,带宽超过36500MB/s(顺序读写速度128000MB/s),每秒读写速度为3240万IOPS(顺序读取速度5000),平均响应时间达到0.11微秒。Calypso PM(持久性内存…

计算机行业研究及2020年策略:聚焦主赛道,投资真成长(84页)

来源:国元证券1、行业涨幅位居前列,个股表现精彩纷呈……2、IT 产业大变革,新机遇不断涌现近年来,以 5G、物联网、云计算、大数据、人工智能、区块链等技术为代表的信息 科技行业快速发展,新兴技术不断涌现&#xff0c…

最长回文串_第78天——第78题(最长回文串 )

今天又是阴天,不过阴天凉快,我喜欢。第78天——第78题(最长回文串)看题目!给定一个包含大写字母和小写字母的字符串,找到通过这些字母构造成的最长的回文串。在构造过程中,请注意区分大小写。比…