PMCW体制雷达系列文章(4) – PMCW雷达之抗干扰

说明

    本文作为PMCW体制雷达系列文章之一,主要聊聊FMCW&PMCW两种体制雷达的干扰问题。事实上不管是通信领域还是雷达领域,对于一切以电磁波作为媒介的信息传递活动,干扰是无处不在的。近年来,随着雷达装车率的提高,车载毫米波雷达之间的干扰问题受到了越来越多的重视,雷达在受到干扰后可能会出现漏检、虚警等问题,这些对于辅助/自动驾驶而言会有极大的风险,不可不察

    规划本话题的初衷是为了说明PMCW雷达相较于FMCW雷达在抗干扰上的优势,不过这句话的意思并不是说PMCW雷达就不会受到干扰,而是PMCW雷达在同等前提下受到干扰的影响要小一些,且其抑制干扰的手段也要多一些,关于这两种体制雷达的干扰影响细节及其区别是本文重点论述的方面。

    本文算是对我之前所写两篇博文:车载毫米波雷达信号处理中的干扰问题-CSDN博客PMCW体制雷达系列文章(1) – PMCW体制雷达综述-CSDN博客的补充,所以建议读者在阅读本文之前,可以先看前面两篇。本文内容主要来源[1],在其基础上加上一些自己的理解所成,参考文献[1]及其references是很好的研究资料,建议感兴趣的读者可以深入研读

Blog

20231119 博文第一次写作

目录

说明

目录

一、雷达间干扰问题概述

二、FMCW雷达之间的干扰

三、PMCW雷达之间的干扰

四、雷达间干扰问题的解决

五、总结

六、参考资料


一、雷达间干扰问题概述

    试想一种简单的行车场景,如图1.1左图所示,自车安装了前雷达,目标车在同一车道的前方,并装有后向雷达,该雷达与自车雷达类似:比如两者的发射功率P、天线增益G、工作频段及调制方式一致。此时对于自车雷达而言,目标车是其需要检测的目标物,而目标车上的雷达所发射的电磁波对自车雷达而言是干扰,该场景下干扰源与目标同距(距离都为R),右图所示的场景类似,只是此时目标(自车前方车辆)与干扰源不同距。

图1.1  雷达干扰示意图[1]

    后续内容分析同距离场景。由雷达方程可得,自车雷达所接收到的目标反射回波功率Pr和干扰源所发射的电磁波功率PI分别为:

                                                (1-1)       

                                                  (1-2)      

    式中,Pt为雷达发射功率,G为天线增益,λ为载波波长,σ为目标车雷达散射截面积(RCS)。(对雷达方程不是很熟悉的读者可以参考我之前的博文:毫米波雷达的系统设计细节(1) - 从雷达方程到雷达威力图_雷达系统设计-CSDN博客)。由两者之比可以得到信号干扰比(signal-to-interference, SIR)为:

                                                     (1-3)

    引入Gp为雷达对目标回波信号的处理增益,该值可以提高雷达的SIR,则上式变成:

                                                   (1-4)

    从目标检测的角度来说,SIR值应该越大越好,但由于干扰相较于正常目标回波由正比于距离平方的能量附加(干扰源发射的电磁波少走了一半的路程),干扰的存在会严重影响目标检测。以式(1-4)为参考,举个例子,假设目标的RCS值为10dBsm(小车的典型RCS值),距离为100m,雷达对目标的信号处理增益为50dB(该值已经很大了,该处理增益主要来源于对ADC数据的二维压缩、相干/非相干积累等),带入上式可得SIR等于9dB,当距离为200m时,SIR为3dB,考虑到小车的RCS值在是实际运动过程中会随入射角、距离等因素发生较大变化(甚至低至0dBsm),而且这里的讨论假定信号处理过程中干扰无增益,总之,这两个例子下得到的SIR值表明:干扰功率已经很接近信号功率,所以目标的正常检测会受到较大影响

    前文的讨论以接近理想情况的场景结合雷达方程对干扰影响做出了大致的评估。对干扰影响的量化,更一般的方法是以干扰噪声比(interference-noise-ratio, INR)来评估,为简化讨论:假定干扰信号的能量像噪声一样均匀分布在被干扰雷达的有效带宽内,这种情况下干扰的影响是最小的(被噪声化了)。我们假设干扰源的发射带宽为B,则干扰源所发射的电磁波信号在被干扰雷达(Victim radar)上的功率谱密度PSDint为:

                            (1-5)

    同样,该式来源于雷达方程。式中,Pt为干扰雷达的发射功率;Gt(Gr)分别是干扰源雷达的发射天线增益以及被干扰雷达接收天线增益;NTx(NRx)对应干扰源的发射天线数量以及被干扰雷达的接收天线数量;Llos_T(Llos_R)为传输损耗,该损耗包含了在雷达内部收发硬件链路上的损耗、电磁波在空气中传播的损耗以及电磁波穿透汽车保险盖和雷达天线罩时的损耗;R为干扰源与被干扰雷达之间的相对距离;DF是干扰源的发射脉冲在被干扰雷达脉冲接收时间段的占比,该值取值范围是[0 1],K是为FMCW体制雷达设定的参数(因为FMCW雷达的下变频以及低通滤波的过程中会将干扰信号与发射信号差频中超过滤波带宽的部分滤除掉,而PMCW雷达是发射单频信号,且我们假定干扰与发射信号载波频率相同,对于PMCW雷达,该K值为1),它的定义为:被干扰雷达混频解调后在基带中的干扰信号的功率谱密度与在射频接收端接收到的干扰信号的功率谱密度之比:

                                                 (1-6)

   式中,是基带中干扰信号的功率谱密度,是射频接收端接收到的干扰信号的功率谱密度。功率谱密度与带宽成反比,所以也可以将之转换成频带宽度之比:是干扰雷达的射频发射带宽,是混频滤波后在基带中的干扰带宽,干扰带宽与干扰源和被干扰雷达两者的调频斜率等有关,关于K值的理解和分析放在了第二章。

    上述内容对干扰的功率谱密度做了定义,被干扰雷达的干扰噪声比INR为:

                                                    (1-7)

    式中,对应噪声的功率谱密度,且:

                                                (1-8)

    式中,k为玻尔兹曼常数,To为环境温度,Fn为接收机噪声系数,接收机的输入噪声功率(热噪)为: Ni = kBTo,B为输入带宽,Fn定义为输入信噪比除以输出信噪比,也即该值正比于输出噪声功率No除以输入噪声功率Ni,于是输出噪声功率No = Fn*Ni,带入Ni = kBTo即可得式(1-8)。

    结合前述各式可知,INR的大小与诸如距离、带宽、干扰源的发射功率等因素有关,下图为基于前述各式的仿真结果[1],(设置K=1):

图1.2  INR仿真结果

    图1.2表征的是被干扰雷达正常波束指向(下图,波束指向0°)以及波束指向2°(上图)时不同距离下的INR值。可以看到INR的值整体上是很大的!且随着距离的减小而增加,小于150米时该值大于10dB,高于噪声10dB的干扰功率会极大地影响雷达的目标检测:包括淹没回波中能量较低的弱目标或较远距离的目标,这会造成雷达漏检、降低雷达的最远检测距离等问题。

    以上算是建立了干扰的数学模型和分析框架,并对干扰问题做了基础探讨,后文分别细化FMCW与PMCW两种体制雷达下的干扰问题。

二、FMCW雷达之间的干扰

    FMCW雷达发射调频连续波(chirp)信号对目标进行探测,目标回波在接收端与本振信号进行混频、滤波并被ADC采集,对采集的时域信号进行傅里叶变换后得到目标对应的频率信息并完成对目标距离的检测,正常情况且不考虑噪声时前述处理过程如下图所示:

图2.1  FMCW体制雷达正常目标回波信号及其处理示意图

    当FMCW雷达间存在干扰时,干扰的波形和本振之间会有不同的“位置关系”,如下图所示:

图2.2  FMCW体制雷达间的不同干扰示意图

    上图例举了FMCW雷达间可能的三种干扰情况,第3种情况下干扰源与目标回波完全一致,其结果时产生无法与正常目标区分的虚假目标,这是雷达虚警的来源之一,不过在实际应用中,考虑到发射信号与干扰信号发射周期、帧周期等可能不一样,首先这种情况发生的概率极低,此外,随着自车或干扰源所在车辆的运动,干扰信号与发射信号之间的相对时序关系等会发生变化,从而使得这种情况下检出的虚假目标不具备连续性,该虚假目标会被后续的跟踪等算法很容易地滤除掉,本节主要考虑前两种更具一般性的干扰情况。

    前两种情况下干扰信号的调频斜率与发射信号不一样,其中情况1所示干扰信号的斜率为负,这两种情况下因为两信号混频后的差频频率会从某个较大值逐渐减小至0,随后逐渐增加,所以干扰信号的能量会分布在较大带宽(当然,小于滤波带宽)内,这两种情况都会导致局部或全部滤波带宽频域内的噪底抬升,如下图所示:

图2.3  干扰情况下的时域&频域示意图

    图中显示干扰信号在频带内的功率要高于底噪15~20dB,在目标检测中,等价于噪底抬升了15~20dB,噪底的抬高会极大地影响雷达的目标检测:包括淹没回波能量较低的弱目标或远距处的目标,这会造成漏检、降低雷达的最远检测距离等问题。

   在前文第一章被干扰雷达接收到的干扰信号功率谱密度公式(1-5)中,我们在雷达方程的基础上增加了参数DF表征干扰源的发射脉冲在被干扰雷达脉冲接收时间段的占比,该值取值范围是[0 1],这从上图2.2中很容易理解(对于PMCW雷达也是这样)。此外,对于FMCW雷达,我们增加了参数K,并将之定义为干扰信号的发射带宽和基带信号中受干扰带宽的比值,不难理解的是,在发射带宽一定的前提下,基带信号中的受干扰带宽越大,对应干扰信号与雷达发射信号的差频“快速地贯穿”整个滤波频带,K越小,不同K值下的干扰对比如下图所示[1]:

图2.4  不同K值下的干扰影响对比示意图

    如上图所示,由于有带通滤波器的存在,两信号的差频只有在滤波器带宽内才会被ADC采集到,当K=1时,其差频快速地从滤波器的截止频率(负)变化为0,然后再增长至截止频率,K=10时,两信号的最大差频仍处在滤波器的截止频率内,所以干扰信号的存在时间甚至可以贯穿整个采样时宽,这种情况下对被干扰雷达的影响更大。

    以上分析只针对单个干扰信号对雷达的影响,实际应用中问题会更复杂:不同车辆上的雷达所发射的信号都可能对本车雷达造成干扰,且这些干扰对应的K值也不尽相同。

三、PMCW雷达之间的干扰

    PMCW雷达间的干扰相较于FMCW要简单一些,假设干扰源与受干扰雷达的载频、相位调制速率、以及带宽都完全一样,只是两者的编码方式不同。此时干扰信号再被干扰雷达内经过混频、低通滤波以及匹配滤波(因为编码方式不一样,所以干扰信号并不会产生正常目标回波那种匹配滤波后的“尖峰”)处理后在时域和频域都类似噪声分布,如下图所示:

图3.1  PMCW雷达间干扰影响

    还是如前式(1-5)所给出的,不同于FMCW雷达间的干扰,PMCW下干扰信号的能量被受干扰雷达接收的大小只与DF值有关,而FMCW雷达由于去斜接收会将频带进一步压缩至滤波器的窄带宽内,由此导致雷达基带中的干扰信号功率谱密度增加,进而恶化干扰噪声比INR,两种情况的对比如下述两式所示:

                                              (3-1)

                                          (3-2)

       

    式中,是受干扰雷达接收到的干扰信号功率,Bi是该接收到的干扰信号的频宽。对于FMCW雷达,接收信号会被进一步压缩至窄带宽内,对应K>1,而导致干扰的功率谱密度增加,进而提高INR,这也从底层原理上解释了为什么PMCW体制雷达相较于FMCW体制雷达更不易受干扰

    此外,不难理解的是,PMCW体制雷达也更容易避免干扰:FMCW雷达的射频信号需要占用很大的带宽(所用带宽越大其距离分辨率越高),而车载雷达可用频段是受限定的:76-81GHz,于是随着雷达数量的增加雷达间的相互干扰会不可避免地恶化,而PMCW雷达采用单频作载频,车载雷达的带宽有5GHz,所以不同雷达间有极大的空间可以错开载频,其次,得益于通信领域伪随机编码技术的使用,雷达间可以使用不同的编码方式、编码长度、相位调制速率等方式很好地相互隔绝,规避干扰。

四、雷达间干扰问题的解决

    随着雷达在车上的装配率越来越高,雷达间的干扰问题也越来越多地引起关注,关于干扰影响的减弱或消除的研究也逐渐增加,关于干扰的消除/减弱我在博文车载毫米波雷达信号处理中的干扰问题-CSDN博客中有过比较详细的说明,本章内容与之相互补充,总结来说,主要包括以下几个解决方向(更细节的一些研究、仿真和实测数据验证等工作读者可以参考[1]的参考文献):

    一是通过各种手段来规避干扰:1、通过改变雷达天线极化方向的方式使得雷达无法接收(或只能接收极少能量)与之不同极化方向干扰源的干扰信号;2、保持雷达发射脉冲的周期或帧周期不变,改变(随机/伪随机化)雷达发射的单脉冲或帧在其周期中的位置来避免干扰;3、因为雷达的带通滤波带宽恒定,可以在可用频带范围内改变发射信号的载频使得干扰信号与本雷达的差频超过滤波器频带范围而被滤掉,从而避免干扰。

    二是检测到干扰,并采取一定的手段来减弱干扰的影响:1、在所采集的ADC时域数据中,通过分析并找到干扰对应的位置,随后对这些数据进行修复(如置0、插值等手段)以减小干扰的影响;2、可以在后端的数据处理中通过算法获取干扰源所在的方向信息,随后在发射端使用波束成形等方法将雷达的波束指向远离干扰的方向以减弱干扰的影响。

    三是系统性地解决该问题:通过推动相关法规从源头上规避雷达间的相互干扰,比如可以给角雷达、前雷达分配不同频带,规定其使用不同的调制方式、编码方式等手段来系统地解决该问题。

五、总结

    本文论述雷达间的干扰问题,首先针对最简单的应用场景对干扰问题进行了数学建模,随后在干扰功率谱密度公式的指导下,比较细节地对FMCW雷达间、PMCW雷达间的干扰情况进行了论述,并分别从更不易受干扰影响以及更容易规避干扰两个角度探讨了PMCW体制雷达相较于FMCW体制雷达在抗干扰上的优势。本文的最后给出了关于抗干扰手段的总结。

六、参考资料

[1]Allond,Stephen,Stark,et al.Interference in Automotive Radar Systems: Characteristics, mitigation techniques, and current and future research[J].[2023-11-14].

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

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈qq邮箱809451989@qq.com,一经查实,立即删除!

相关文章

Python---练习:封装一个函数,用于生成指定长度的验证码

练习涉及相关链接:Python---练习:编写一段Python代码,生成一个随机的4位验证码-CSDN博客 Python----函数中的说明文档-CSDN博客Python---return返回值-CSDN博客 代码: # 定义一个generate_code()函数 def generate_code(num): …

SDUT OJ《算法分析与设计》搜索算法

A - 子集和问题 Description 子集和问题的一个实例为〈S,t〉。其中,S{ x1 , x2 ,…,xn }是一个正整数的集合,c是一个正整数。子集和问题判定是否存在S的一个子集S1,使得: 。 试设计一个解子…

【Web】Ctfshow SSTI刷题记录1

目录 ①web361 362-无过滤 ②web363-过滤单双引号 ③web364-过滤单双引号和args ④web365-过滤中括号[]、单双引号、args ⑤web366-过滤单双引号、args、中括号[]、下划线 ⑦web367-过滤单双引号、args、中括号[]、下划线、os ⑧web368-过滤单双引号、args、中括号[]、下…

【嵌入式 – GD32开发实战指南(ARM版本)】第2部分 外设篇 - 第3章 温度传感器DS18B20

1 理论分析 1.1 DS18B20概述 DS18B20 是 DALLAS 最新单线数字温度传感器,新的"一线器件"体积更小、适用电压更宽、更经济。Dallas 半导体公司的数字化温度传感器 DS1820 是世界上第一片支持 "一线总线"接口的温度传感器。 DS18B20采用的单总线协议,也…

Centos7 重置 Root 密码

Centos7 重置 Root 密码 1.启动服务器2.编辑启动项3.修改密码4.重新登陆 1.启动服务器 启动服务器后,不要直接进入系统,在开机页面按键盘【E】 2.编辑启动项 按【E】后进入如下页面,并按向下箭头,找到如图位置,添加如…

数据结构:红黑树的插入实现(C++)

个人主页 : 个人主页 个人专栏 : 《数据结构》 《C语言》《C》《Linux》 文章目录 一、红黑树二、红黑树的插入三、代码实现总结 一、红黑树 红黑树的概念: 红黑树是一颗二叉搜索树,但在每个节点上增加一个存储位表示节点的颜色&…

慢日志查询

概述 MySQL的慢查询日志是MySQL提供的一种日志记录,它用来记录在MySQL中响应时间超过阀值的语句,具体指运行时间超过 long_query_time 值的SQL,则会被记录到慢查询日志中,ong_query_time 的默认值为 10,意思是运行10S…

redis三种集群方式

redis有三种集群方式:主从复制,哨兵模式和集群。 1.主从复制 主从复制原理: 从服务器连接主服务器,发送SYNC命令; 主服务器接收到SYNC命名后,开始执行BGSAVE命令生成RDB文件并使用缓冲区记录此后执行的所…

python-opencv 培训课程笔记(2)

python-opencv 培训课程笔记(2) 1.图像格式转换 先看一下cvtColor函数的例子 #默认加载彩图 pathrD:\learn\photo\cv\cat.jpg# imread(path,way) #way0 灰度图。way1 彩图 #默认彩图 imgcv2.imread(path) img_dogcv2.imread(path_dog) #图片格式的转化…

【开发流程】持续集成、持续交付、持续部署

一、开发工作流程 假设把开发流程分为以下几个阶段: 编码 -> 构建 -> 集成 -> 测试 -> 交付 -> 部署 如上图所示,持续集成、持续交付、持续部署有着不同的软件自动交付周期。 二、持续集成、持续交付、持续部署 1、持续集成 持续集成…

FPGA设计时序约束八、others类约束之Set_Case_Analysis

目录 一、序言 二、Set Case Analysis 2.1 基本概念 2.2 设置界面 2.3 命令语法 2.4 命令示例 三、工程示例 四、参考资料 一、序言 在Vivado的时序约束窗口中,存在一类特殊的约束,划分在others目录下,可用于设置忽略或修改默认的时序…

关于DBMS_STATS.GATHER_DATABASE_STATS_JOB_PROC的一些发现

任务在哪 这个是11g以后的自动收集统计信息的后台任务,10g之前是在dba_scheduler_jobs里查看 SQL> SELECT CLIENT_NAME ,STATUS ,MEAN_INCOMING_TASKS_7_DAYS,MEAN_INCOMING_TASKS_30_DAYS FROM DBA_AUTOTASK_CLIENT WHERE…

mac控制台命令小技巧

shigen日更文章的博客写手,擅长Java、python、vue、shell等编程语言和各种应用程序、脚本的开发。记录成长,分享认知,留住感动。 hello伙伴们,作为忠实的mac骨灰级别的粉丝,它真的给我带来了很多效率上的提升。那作为接…

算法通关村——数字中的统计、溢出、进制转换处理模板

数字与数学基础问题 1、数字统计 1.1、符号统计 LeetCode1822. 给定一个数组,求所有元素的乘积的符号,如果最终答案是负的返回-1,如果最终答案是正的返回1,如果答案是0返回0. 这题其实只用看数组中0和负数的个数就好了&#x…

SpirngBoot + Vue 前后端分离开发工具代码

✅作者简介:大家好,我是Leo,热爱Java后端开发者,一个想要与大家共同进步的男人😉😉 🍎个人主页:Leo的博客 💞当前专栏: Java从入门到精通 ✨特色专栏&#xf…

股东入股可用的出资形式主要有哪些

股东入股,可用的出资形式主要包括货币以及实物、知识产权、土地使用权等可以用货币估价并可以依法转让的非货币财产。 第一,货币。设立公司必然需要一定数量的流动资金。以支付创建公司时的开支和启动公司运营。因此,股东可以用货币出资。 第…

Kafka性能测试初探

相信大家对Kafka不会陌生,但首先还是要简单介绍一下。 Kafka是一种高性能的分布式消息系统,由LinkedIn公司开发,用于处理海量的实时数据流。它采用了发布/订阅模式,可以将数据流分发到多个消费者端,同时提供了高可靠性…

uniapp 微信小程序登录 新手专用 引入即可

预览 第一步导入插件 在引入的页面的登录按钮下拷贝一下代码 <template><view class"content"><button type"primary" click"login">微信登录</button></view><TC-WXlogin :wxloginwxlogin /> </templ…

阿坤老师的独特瓷器(Java详解)

一、题目描述 示例&#xff1a; 输入&#xff1a; 5 3 4 5 6 2 5 3 7 6 5 输出&#xff1a; 3 二、题解 思路分析&#xff1a; 题目要求我们计算出“独特瓷器”的个数&#xff0c;而“独特瓷器”是指对于瓷器A&#xff0c;没有另一个瓷器B&#xff0c;直径和高度都大于A。则…

Idea远程debug

Idea远程debug 一、方法 1.1、启动时加参数 java -agentlib:jdwptransportdt_socket,servery,suspendn,address5005 -jar remote-debug-0.0.1-SNAPSHOT.jar1.2、运行时加参数 在没有显式-agentlib:jdwp参数的情况下进行远程调试&#xff1a; 首先&#xff0c;确保您的Sprin…