【FMCW系统性能参数之测量精度公式推导】

news2025/2/24 20:20:19

本文编辑:调皮哥的小助理

连续多篇文章都在说FMCW雷达系统性能参数这个事儿,如:

(1)从奈奎斯特采样定理推导FMCW雷达系统性能参数
(2)从FMCW毫米波雷达系统的性能参数理解4D成像毫米波雷达的设计思路
(3)雷达入门 | FMCW毫米波雷达系统的性能参数分析

因为目标检测与参数估计这些是雷达系统的“第一性原理的第二个层次”,所以需要搞得清清楚楚,明明白白的。关于雷达的“第一性原理的第一个层次”,我认为是电磁场与电磁波(麦克斯韦方程组)、天线等,这才是雷达的真正源头与底层逻辑。无论是日常的科学研究,还是将来的找工作,我想这些基本的理论,是会经常被遇到的。

之前在(3)文章中提到过FMCW雷达系统性能参数的简要介绍,其中关于距离、速度和角度的测量精度直接给出了公式,缺乏详细的推导过程,正好之前在【雷达技术交流群】中有读者提到过,于是我觉得有必要分享。

下面来看看详细的过程,首先给出结论,需要注意的是本文仅针对的雷达体制是FMCW雷达,而不是脉冲雷达,虽然文章的末尾会给出不同波形的介绍。注意复信号和实信号的信噪比计算。

测距精度:

σ R = c 3.6 B 2 S N R \sigma_R=\frac{c}{3.6 B \sqrt{2 S N R}} σR=3.6B2SNR c

测速精度:

σ v = λ 3.6 N T c S N R \sigma_v=\frac{\lambda}{3.6 N T_c \sqrt{S N R}} σv=3.6NTcSNR λ

测角精度:

σ θ = θ 3 d B 1.6 2 S N R \sigma_\theta=\frac{\theta_{3 d B}}{1.6 \sqrt{2 S N R}} σθ=1.62SNR θ3dB

要得到公式的推导,需要参考比较老的文献,其中一本书是丁鹭飞的《雷达系统》,是1984年的出版的;另一本是李蕴滋的《雷达工程学》,这本书可以说比较适合入门了。其中,前一本书有详细的证明过程,后一本书直接给出了结论,所以我推荐先看后面一本,然后再去看详细的过程。

这里就直接看丁鹭飞的雷达系统》了,具体在第二章的雷达测量精度和分辨率。这本书可以在网上购买旧版的,价格也不贵,电子版在会员资料库中。

在这里插入图片描述

其中,在2.1小节就提到了最大似然估计法,这个方法用于雷达系统对目标进行参数估计,之前有读者提到距离估计所做的距离维FFT在某种程度上可以看作是最大似然估计,其实我现在觉得他说的对,不光是距离,就连速度、角度等估计都是。而今天所要讨论的雷达测量理论精度,就是建立在最大似然估计法上的。

1、测距精度

测距精度其实就是时延估计精度,在丁鹭飞的《雷达系统》一书中,先以脉冲雷达为例,按照最大似然估计法,测距时候首先计算相关积分(脉冲压缩),然后求出最大值出现的时间即为时延估计值。步骤如下:

(1)时延估计

(2)最大似然估计

(3)精度分析

最终经过一系列数学计算,就得到了脉冲雷达下的测距精度,如下所示。

σ R = c 2 1 B 2 S N R \sigma_R=\frac{c}{2} \frac{1}{B \sqrt{2 S N R}} σR=2cB2SNR 1

其中,由于该过程公式复杂,内容繁多,不便于在这里展示,请各位读者详细看书,这里仅是一个引导作用,还请谅解。另外距离估计精度时间估计精度的系数乘以c/2。

但是这个公式对不同的波形表达不同,比如矩形脉冲、高斯脉冲、辛克波形以及调频连续波。这里对调频连续波进行推导后,得到测距精度为:

σ R = 3 c 2 π B 2 E N 0 \sigma_R=\frac{\sqrt{3} c}{2 \pi B \sqrt{\frac{2 E}{N_0}}} σR=2πBN02E 3 c

简化之后近似表示为:

σ R = c 3.6 B 2 S N R \sigma_R=\frac{c}{3.6 B \sqrt{2 S N R}} σR=3.6B2SNR c

也可以表示为:

σ R = Δ R 1.8 2 S N R \sigma_R=\frac{\Delta R}{1.8 \sqrt{2 S N R}} σR=1.82SNR ΔR

也就是说,距离测量理论精度和发射信号波形、测量时得信噪比SNR、信号的带宽(理论距离分辨率、瑞利距离分辨率)有关,带宽越大、信噪比越高,测量精度越高。

2、测速精度

同理,速度测量精度主要是估计多普勒频移,也就是说速度测量精度,可以等效为多普勒频率的估计精度,多普勒频率估计精度与测速精度的系数是半波长。推导的方法和步骤与距离测量精度基本一样,具体过程也是非常繁杂,需要详细看书。

这里给出最终估计的速度误差(精度)表示为:

σ v = λ 2 π T e 2 E N 0 = λ 2 π N T c 2 S N R , T e = N T c \sigma_v=\frac{\lambda}{2 \pi T_e \sqrt{\frac{2 E}{N_0}}}=\frac{\lambda}{2 \pi N T_c \sqrt{2 S N R}}, T_e=N T_c σv=2πTeN02E λ=2πNTc2SNR λ,Te=NTc

其中Te是信号持续时间,Tc是单Chirp扫频时间,N是表示Chirp数。同样也是不同的波形,测速精度表达方式不一样, 如矩形脉冲、高斯脉冲、辛克波形、线性调频连续波以及相参脉冲串。其他波形看书,这里给出线性调频连续波的最终表达公式如下:

σ v = λ 3 2 π T e 2 E N 0 = λ 3 2 π N T c 2 S N R = λ 3.6 N T c 2 S N R \sigma_v=\frac{\lambda \sqrt{3}}{2 \pi T_e \sqrt{\frac{2 E}{N_0}}}=\frac{\lambda \sqrt{3}}{2 \pi N T_c \sqrt{2 S N R}}=\frac{\lambda}{3.6 N T_c \sqrt{2 S N R}} σv=2πTeN02E λ3 =2πNTc2SNR λ3 =3.6NTc2SNR λ

也就是说,速度测量理论精度和发射信号波形、测量时得信噪比SNR、信号的持续时间(理论速度分辨率、瑞利速度分辨率)有关,信号总持续时间越长、信噪比越高,测量精度越高。 同时波长越小,精度也越高,这也是77G雷达优于24G雷达的主要优势之一。

另外,还有距离速度联合估计精度,相关内容见书中。

3、测角精度

方法与前面二者大致相同,略有区别。当口面场为等幅度分布时,测角误差可以表示为:

σ θ = λ 3 π D 2 E N 0 \sigma_\theta=\frac{\lambda \sqrt{3}}{\pi D \sqrt{\frac{2 E}{N_0}}} σθ=πDN02E λ3

可知,测角精度与信噪比和天线的归一化口径D/λ有关,波长降低也能提高测角精度。 一般在3dB波束宽度或者称为半功率点波束宽度,可以得到D/λ的关系为:

θ 3 d B = 0.88 λ D \theta_{3 d B}=0.88 \frac{\lambda}{D} θ3dB=0.88Dλ

推导出:

σ θ = θ 3 d B 3 0.88 π 2 E N 0 \sigma_\theta=\frac{\theta_{3 d B} \sqrt{3}}{0.88 \pi \sqrt{\frac{2 E}{N_0}}} σθ=0.88πN02E θ3dB3

这个时候得到的测角精度为:

σ θ = θ 3 d B 3 0.88 π 2 E N 0 = θ 3 d B 1.6 2 S N R \sigma_\theta=\frac{\theta_{3 d B} \sqrt{3}}{0.88 \pi \sqrt{\frac{2 E}{N_0}}}=\frac{\theta_{3 d B}}{1.6 \sqrt{2 S N R}} σθ=0.88πN02E θ3dB3 =1.62SNR θ3dB

所以更加进一步的结论是3dB波束越窄,也就是角度分辨率越低,测角精度越高。因此,无论是测距、测速还是测角,分辨率和信噪比都会影响测量精度的。

感谢大家的关注,博主会继续学习,争取分享更多有益于大家共同学习的文章,更多资料详情大家也可移步https://mp.weixin.同步查看!

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

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

相关文章

深度学习语义分割篇——FCN原理详解篇

🍊作者简介:秃头小苏,致力于用最通俗的语言描述问题 🍊往期回顾:目标检测系列——开山之作RCNN原理详解    目标检测系列——Fast R-CNN原理详解    目标检测系列——Faster R-CNN原理详解 🍊近期目标&…

Unity --- 3d数学 --- 坐标系统

1.世界坐标系是固定不动的 2.每一个游戏物体在世界坐标系中都有对应的坐标和方向 1.轴心点的位置不是固定的,是可以人为设定的 1.Screen Space --- 屏幕坐标 2.我们看到的屏幕其实就是相机所在的平面的位置 --- 而屏幕坐标系的Z其实就是游戏中的物体到相机平面的…

GDOUCTF

WEB hate eat snake 这是一个JS的题目,但是这个题目好像有点奇怪,不是很理解,当时我找到了我寝室JS的大哥,跟大哥说了一下我的思路,就是他根据这个time然后/1000转化为秒,就当作是我们玩游戏的一个分数&a…

速卖通韩国下载量再次登顶,7500万投资换来回报

韩国市场,还是一片蓝海。 速卖通终于等到了回报。 近日,数据平台 Data ai 显示,3 月 9 日以来,速卖通再次成为韩国购物 App 下载量第一名,超过当地电商平台 Coupang。 这或许和速卖通近日在韩国的布局有关。 前些天…

使用File System Access API 让浏览器可以操作文件

使用File System Access API 让浏览器可以操作文件 在早期我们开始学习前端三件套时,经常会听到这样的说法:浏览器是一个沙盒,它不允许我们操作本地文件,但是现在这个说法已经不再适用了,因为我们可以使用 File Syste…

2023接口自动化测试,完整入门篇

1. 什么是接口测试 顾名思义,接口测试是对系统或组件之间的接口进行测试,主要是校验数据的交换,传递和控制管理过程,以及相互逻辑依赖关系。其中接口协议分为HTTP,WebService,Dubbo,Thrift,Socket等类型,测试类型又主…

BBR算法

BBR算法 简述 bbr算法为google在2016年提出,用于改善tcp的性能,提升稳定性,降低延迟,更好地应对网络损伤。在整个算法调节周期中,bbr算法都在尽力维持最大bw和最小rtt。 对比传统的tcp算法 传统算法不能区分是拥塞导…

【RabbitMQ】Java操作RabbitMQ之入门Demo

目录 一、项目创建 二、生产者 三、消费者 一、项目创建 我们先在idea里创建两个Maven项目一个项目作为生产者&#xff0c;另一个作为消费者。创建完成后&#xff0c;在各自的pom.xml文件里引入Java使用RabbitMQ的依赖 <dependency><groupId>com.rabbitmq</g…

Android 单位和尺寸 以及ViewPager的使用

一 px 和 pt 1 px pixels 像素点 2 pt point 一个标准的长度单位&#xff0c;1pt 1/72英寸&#xff0c;用于印刷业&#xff0c;非常简单易用。 二 dp 和 sp 1 dp 就是 dip device independent pixels(设备独立像素) 不同设备有不同的显示效果。这个和设备硬件有关系&…

【数学推导找规律】全排列的价值【蓝桥杯13届】

欢迎观看我的博客&#xff0c;如有问题交流&#xff0c;欢迎评论区留言&#xff0c;一定尽快回复&#xff01;&#xff08;大家可以去看我的专栏&#xff0c;是所有文章的目录&#xff09;   文章字体风格&#xff1a; 红色文字表示&#xff1a;重难点★✔ 蓝色文字表示&#…

ChatGPT资讯—2023.4.14

一、 最新资讯 1. 国内大模型&#xff1a; 4月8日&#xff0c;华为云人工智能领域首席科学家田奇首次介绍了华为“盘古大模型”的进展及其应用。 紧接着4月10日&#xff0c;一天内三个大模型争相“官宣”&#xff0c;包括搜狗创始人王小川宣布即将投入创业的大模型产品名为“…

机器学习 day04(梯度下降算法,学习率,偏导数)

梯度下降 我们可以用一种更系统的方法&#xff0c;来找到一组w&#xff0c;b&#xff0c;使成本函数的值最小。这个方法叫梯度下降算法&#xff0c;它可用于最小化任何函数&#xff0c;不仅仅包括线性回归的成本函数&#xff0c;也包括两个以上参数的其他成本函数在线性回归中…

2023最新面试题-Java-4

Date相关 1. java8的Date相关API&#xff1a; 常用 api 1、 获取当前日期 LocalDate.now() 2、创建日期 LocalDate date LocalDate.of(2020, 9, 21) 3、获取年份 date.getYear()//通过 TemporalField 接口的实现枚举类 ChronoField.YEAR 获取年份 date.get(ChronoFie…

车载通信——J1939 DM1

诊断故障代码&#xff08;DTC&#xff09;由4 个独立域构成&#xff1a; a. 可疑参数的编号(SPN) 19位 b. 故障模式标志(FMI) 5位 c. 发生次数(OC) 7位 d. 可疑参数编号的转化方式(CM) 1 位 Byte1 – Lamp Status Bits1 – 2 Protection Lamp 保护灯状态 …

基于 NIOSII 软核的流水灯实验

文章目录一、硬件部分设计1、进行 Qsys 系统设计二、逻辑连接三、软件设计四、总结一、硬件部分设计 1、进行 Qsys 系统设计 新建kernel.qsys文件&#xff1a; 设置clk时钟&#xff1a; 添加 Nios II 其余选项保持默认 添加 jtag uart 接口 添加片上存储器 On-Chip Memory(R…

【C++高级】手写线程池项目-经典死锁问题分析-简历项目输出指导

作为五大池之一&#xff0c; 线程池的应用非常广 泛&#xff0c;不管是客户端程序&#xff0c;还是后台服务程序&#xff0c;掌握线程池&#xff0c;是提高业务处理能力的必备模块 本课程将带你从零开始&#xff0c;设计一个支持fixed和cached模式的线程池&#xff0c;玩转C11、…

使用FRP(快速反向代理)实现内网穿透——以腾讯云服务器为例

一、FRP简介 FRP&#xff0c;即快速反向代理技术&#xff08;fast reverse proxy&#xff09;。本文的FRP程序是基于github开源项目GitHub - fatedier/frp。当前&#xff0c;该程序可实现&#xff1a;“将位于 NAT 或防火墙后面的本地服务器暴露给互联网”。它目前支持 TCP 和…

金山衍生新软件,wps,excel用户已在用,Access用户:以后就它了

我们常会用到微软的办公软件Word、Excel、PPT、Outlook&#xff0c;但在数据处理方面还是得看Access。Access用简短的表述来说就是微软开发的一个关系数据库管理系统。★好用&#xff0c;可门槛高&#xff0c;够不着 为啥说数据处理得看Access呢&#xff1f;举个例子&#xff0…

如何实现视觉识别形状

1. 功能说明 通过摄像头识别圆形及矩形两种形状。 2. 电子硬件 本实验中采用了以下硬件&#xff1a; 主控板 Basra主控板&#xff08;兼容Arduino Uno&#xff09; 扩展板 Bigfish2.1 电池7.4V锂电池通信2510通信转接板WiFi路由器 其它 摄像头 配置OpenCV的Visual Studio 2015.…

深入了解网络通信原理

同一个 vlan 同一个网段的 ip 地址能够直接互通&#xff0c;那么同一个 vlan 不同的网段能不能互通呢&#xff1f;海翎光电的小编整理了一篇文章&#xff0c;可以解决大家在网络中遇到的一些奇怪的问题&#xff0c;也是可以加深对网络的理解。 一、同 VLAN 不同网段能否 Ping 通…