【小问题】距离估计和频率估计的方差下界推导出距离估计的方差下界
【1】OFDM Radar Algorithms in Mobile Communication Networks pp34
文章目录
如何根据 d ^ = n ^ c 0 2 Δ f N P e r \hat{d} = \frac{\hat{n}c_0}{2\Delta f N_{\mathrm{Per}}} d^=2ΔfNPern^c0和 var [ ω ^ ] ≥ 6 σ N 2 ( N 2 − 1 ) N \operatorname{var}[\hat{\omega}] \geq \frac{6\sigma_N^2}{(N^2-1)N} var[ω^]≥(N2−1)N6σN2推出 var [ d ^ ] ≥ 6 σ N 2 ( N 2 − 1 ) N ( c 0 4 π Δ f ) 2 \operatorname{var}[\hat{d}] \geq \frac{6\sigma_N^2}{(N^2-1)N} \left( \frac{c_0}{4\pi \Delta f} \right)^2 var[d^]≥(N2−1)N6σN2(4πΔfc0)2
要从给定的距离估计公式和频率估计的方差下界推导出距离估计的方差下界,我们可以按照以下步骤进行:
1. 频率和距离之间的关系
已知距离估计的公式为:
d ^ = n ^ c 0 2 Δ f N P e r \hat{d} = \frac{\hat{n} c_0}{2 \Delta f N_{\mathrm{Per}}} d^=2ΔfNPern^c0
其中:
- n ^ \hat{n} n^ 是从周期图中检测到的频率索引。
- c 0 c_0 c0 是信号传播速度。
- Δ f \Delta f Δf 是频率分辨率。
- N P e r N_{\mathrm{Per}} NPer 是周期图的样本数量。
2. 计算 d ^ \hat{d} d^ 对 n ^ \hat{n} n^ 的导数
我们可以计算 d ^ \hat{d} d^ 对 n ^ \hat{n} n^ 的偏导数,以便将其用于方差的转换:
∂ d ^ ∂ n ^ = c 0 2 Δ f N P e r \frac{\partial \hat{d}}{\partial \hat{n}} = \frac{c_0}{2 \Delta f N_{\mathrm{Per}}} ∂n^∂d^=2ΔfNPerc0
3. 将频率的方差转化为距离的方差
根据 CRB 的转换关系:
var [ d ^ ] = ∣ ∂ d ^ ∂ n ^ ∣ 2 var [ n ^ ] \operatorname{var}[\hat{d}] = \left| \frac{\partial \hat{d}}{\partial \hat{n}} \right|^2 \operatorname{var}[\hat{n}] var[d^]= ∂n^∂d^ 2var[n^]
将上面计算的导数代入:
var [ d ^ ] = ( c 0 2 Δ f N P e r ) 2 var [ n ^ ] \operatorname{var}[\hat{d}] = \left( \frac{c_0}{2 \Delta f N_{\mathrm{Per}}} \right)^2 \operatorname{var}[\hat{n}] var[d^]=(2ΔfNPerc0)2var[n^]
4. 从频率的 CRB 获得 var [ n ^ ] \operatorname{var}[\hat{n}] var[n^]
根据频率估计的 Cramér-Rao Bound,我们有:
var [ ω ^ ] ≥ 6 σ N 2 ( N 2 − 1 ) N \operatorname{var}[\hat{\omega}] \geq \frac{6\sigma_N^2}{(N^2-1)N} var[ω^]≥(N2−1)N6σN2
注意到频率 ω ^ \hat{\omega} ω^ 与索引 n ^ \hat{n} n^ 的关系为:
ω ^ = 2 π n ^ N P e r \hat{\omega} = \frac{2\pi \hat{n}}{N_{\mathrm{Per}}} ω^=NPer2πn^
因此:
n ^ = N P e r ω ^ 2 π \hat{n} = \frac{N_{\mathrm{Per}} \hat{\omega}}{2\pi} n^=2πNPerω^
推导方差的关系:
var [ n ^ ] = ( N P e r 2 π ) 2 var [ ω ^ ] \operatorname{var}[\hat{n}] = \left( \frac{N_{\mathrm{Per}}}{2\pi} \right)^2 \operatorname{var}[\hat{\omega}] var[n^]=(2πNPer)2var[ω^]
5. 将 var [ n ^ ] \operatorname{var}[\hat{n}] var[n^] 代入距离方差公式
将 var [ ω ^ ] \operatorname{var}[\hat{\omega}] var[ω^] 的 CRB 代入:
var [ n ^ ] ≥ ( N P e r 2 π ) 2 6 σ N 2 ( N 2 − 1 ) N \operatorname{var}[\hat{n}] \geq \left( \frac{N_{\mathrm{Per}}}{2\pi} \right)^2 \frac{6\sigma_N^2}{(N^2-1)N} var[n^]≥(2πNPer)2(N2−1)N6σN2
将其代入到距离方差公式中:
var [ d ^ ] ≥ ( c 0 2 Δ f N P e r ) 2 ( N P e r 2 π ) 2 6 σ N 2 ( N 2 − 1 ) N \operatorname{var}[\hat{d}] \geq \left( \frac{c_0}{2 \Delta f N_{\mathrm{Per}}} \right)^2 \left( \frac{N_{\mathrm{Per}}}{2\pi} \right)^2 \frac{6\sigma_N^2}{(N^2-1)N} var[d^]≥(2ΔfNPerc0)2(2πNPer)2(N2−1)N6σN2
6. 简化公式
最终的方差表达式变为:
var [ d ^ ] ≥ 6 σ N 2 ( N 2 − 1 ) N ( c 0 4 π Δ f ) 2 \operatorname{var}[\hat{d}] \geq \frac{6\sigma_N^2}{(N^2-1)N} \left( \frac{c_0}{4\pi \Delta f} \right)^2 var[d^]≥(N2−1)N6σN2(4πΔfc0)2
这就是你提到的结果:
var [ d ^ ] ≥ 6 σ N 2 ( N 2 − 1 ) N ( c 0 4 π Δ f ) 2 \operatorname{var}[\hat{d}] \geq \frac{6\sigma_N^2}{(N^2-1)N} \left( \frac{c_0}{4\pi \Delta f} \right)^2 var[d^]≥(N2−1)N6σN2(4πΔfc0)2
总结来说,主要步骤是通过计算导数关系,将频率的方差转化为距离的方差,并结合频率的 CRB 得到最终结果。
相关文章:
【小问题】距离估计和频率估计的方差下界推导出距离估计的方差下界
【1】OFDM Radar Algorithms in Mobile Communication Networks pp34 文章目录 1. 频率和距离之间的关系2. 计算 d ^ \hat{d} d^ 对 n ^ \hat{n} n^ 的导数3. 将频率的方差转化为距离的方差4. 从频率的 CRB 获得 var [ n ^ ] \operatorname{var}[\hat{n}] var[n^]5. 将 …...

Selenium爬虫技术:如何模拟鼠标悬停抓取动态内容
介绍 在当今数据驱动的世界中,抓取动态网页内容变得越来越重要,尤其是像抖音这样的社交平台,动态加载的评论等内容需要通过特定的方式来获取。传统的静态爬虫方法难以处理这些由JavaScript生成的动态内容,Selenium爬虫技术则是一…...

Z-BlogPHP显示错误Undefined array key 0 (set_error_handler)的解决办法
今天打开博客的时候,意外发现页面,打开均显示错误:Undefined array key 0 (set_error_handler)。 博客程序采用的是Z-BlogPHP。百度了一圈没有找到解决办法,在官方论坛里也没找到解决办法。 于是开始自己排查原因。我服务器采用…...
java-实例化一个List,然后添加数据的方法详解
在Java中,实例化一个 List 并向其中添加数据非常简单。List 是一个接口,因此我们通常使用它的常见实现类 ArrayList 或 LinkedList。以下是一些常见的操作方法: ### 1. 使用 ArrayList 实例化并添加数据 java import java.util.ArrayList; …...

【Linux系统】Ubuntu的简单操作
什么是 Ubuntu? Ubuntu(乌帮图)是一个非洲词汇,它的意思是“人性对待他人”或“群在故我在”。Ubuntu发行版将Ubuntu精神带到软件世界之中。 目前已有大量各种各样基于GNU/Linux的操作系统,例如:Debian,SuSE,Gentoo,R…...

标准日志插件项目【C/C++】
博客主页:花果山~程序猿-CSDN博客 文章分栏:项目日记_花果山~程序猿的博客-CSDN博客 关注我一起学习,一起进步,一起探索编程的无限可能吧!让我们一起努力,一起成长! 目录 一,项目介…...

SpingBoot原理
SpingBoot原理 在前面十多天的课程当中,我们学习的都是web开发的技术使用,都是面向应用层面的,我们学会了怎 么样去用。而我们今天所要学习的是web后端开发的最后一个篇章springboot原理篇,主要偏向于底 层原理。 我们今天的课程…...

Cout输出应用举例
Cout输出应用 在main.cpp里输入程序如下: #include <iostream> //使能cin(),cout(); #include <stdlib.h> //使能exit(); #include <sstream> #include <iomanip> //使能setbase(),setfill(),setw(),setprecision(),setiosflags()和res…...
java的无锁编程和锁机制
Java 的并发编程中,为了保证线程安全和高性能,采用了两种主要的同步手段:锁机制和无锁编程。以下是对锁机制、无锁编程、死锁及其避免的详细讲解。 一、无锁编程 无锁编程通过原子操作来避免传统锁,从而减少线程的上下文切换&am…...
vue实现富文本编辑器上传(粘贴)图片 + 文字
vue实现富文本编辑器上传(粘贴)图片 文字 1.安装插件 npm install vue-quill-editor -s2.在使用vue-quill-editor富文本的时候,对于图片的处理经常是将图片转换成base64,再上传数据库,但是base64不好存储。 原理&a…...

子集和全排列(深度优先遍历)问题
欢迎访问杀马特主页:小小杀马特主页呀! 目录 前言: 例题一全排列: 1.题目介绍: 2.思路汇总: 3.代码解答: 例题二子集: 题目叙述: 解法一: 1.思路汇总…...
判断检测框是否在感兴趣区域(ROI)内
判断检测框是否在感兴趣区域(ROI)内 在计算机视觉和图像处理中,我们经常需要确定一个矩形检测框是否位于一个特定的感兴趣区域(Region of Interest, ROI)内。这个ROI可以是一个多边形,而检测框则是一个矩形…...

正点原子阿尔法ARM开发板-IMX6ULL(九)——关于SecureCRT连接板子上的ubuntu
文章目录 一、拨码器二、SecureCRT 一、拨码器 emmm,也是好久没学IMX6ULL了,也是忘了拨码器决定了主板的启动方式 一种是直接从TF卡中读取文件(注意这里是通过imdownload软件编译好了之后,通过指令放入TF卡) 一种是现在这种用串口…...
微信支付Java+uniapp微信小程序
JS: request.post(/vip/pay, {//这是自己写的java支付接口id: this.vipInfo.id,payWay: wechat-mini}).then((res) > {let success (res2) > {//前端的支付成功回调函数this.$refs.popup.close();// 支付成功刷新当前页面setTimeout(() > {this.doGetVipI…...

【NOIP提高组】加分二叉树
【NOIP提高组】加分二叉树 💐The Begin💐点点关注,收藏不迷路💐 设一个n个节点的二叉树tree的中序遍历为(l,2,3,…,n),其中数字1,2,3,…,n为节点编号。每个节点都有一个分数(均为正整…...

HarmonyOS 相对布局(RelativeContainer)
1. HarmonyOS 相对布局(RelativeContainer) 文档中心:https://developer.huawei.com/consumer/cn/doc/harmonyos-guides-V5/arkts-layout-development-relative-layout-V5 RelativeContainer为采用相对布局的容器,支持容器内部的子元素设…...
webpack5搭建react脚手架详细步骤
1. 初始化项目 首先,创建一个新目录并初始化项目: bash mkdir create-react cd create-react pnpm init --y git init 这里使用pnpm作为包管理工具,因为它在处理依赖和速度上表现更好。 2. 安装React和TypeScript 安装React和React-DOM…...
速盾:高防cdn怎么拦截恶意ip?
高防CDN(Content Delivery Network)是一种用于防御网络攻击和提供高可用性的服务。它通过分发网络流量,将用户的请求导向最近的服务器,从而提高网站的加载速度和稳定性。然而,不可避免地,有些恶意IP地址会试…...

太阳能面板分割系统:训练自动化
太阳能面板分割系统源码&数据集分享 [yolov8-seg-EfficientHead&yolov8-seg-vanillanet等50全套改进创新点发刊_一键训练教程_Web前端展示] 1.研究背景与意义 项目参考ILSVRC ImageNet Large Scale Visual Recognition Challenge 项目来源AAAI Globa…...

C++笔记---位图
1. 位图的概念 位图(Bitmap)是一种基于位操作的数据结构,用于表示一组元素的集合信息。它通常是一个仅包含0和1的数组,每个元素对应一个二进制位,若该元素存在,则对应的位为1;若不存在ÿ…...

多云管理“拦路虎”:深入解析网络互联、身份同步与成本可视化的技术复杂度
一、引言:多云环境的技术复杂性本质 企业采用多云策略已从技术选型升维至生存刚需。当业务系统分散部署在多个云平台时,基础设施的技术债呈现指数级积累。网络连接、身份认证、成本管理这三大核心挑战相互嵌套:跨云网络构建数据…...
【Linux】C语言执行shell指令
在C语言中执行Shell指令 在C语言中,有几种方法可以执行Shell指令: 1. 使用system()函数 这是最简单的方法,包含在stdlib.h头文件中: #include <stdlib.h>int main() {system("ls -l"); // 执行ls -l命令retu…...
QMC5883L的驱动
简介 本篇文章的代码已经上传到了github上面,开源代码 作为一个电子罗盘模块,我们可以通过I2C从中获取偏航角yaw,相对于六轴陀螺仪的yaw,qmc5883l几乎不会零飘并且成本较低。 参考资料 QMC5883L磁场传感器驱动 QMC5883L磁力计…...

对WWDC 2025 Keynote 内容的预测
借助我们以往对苹果公司发展路径的深入研究经验,以及大语言模型的分析能力,我们系统梳理了多年来苹果 WWDC 主题演讲的规律。在 WWDC 2025 即将揭幕之际,我们让 ChatGPT 对今年的 Keynote 内容进行了一个初步预测,聊作存档。等到明…...
生成 Git SSH 证书
🔑 1. 生成 SSH 密钥对 在终端(Windows 使用 Git Bash,Mac/Linux 使用 Terminal)执行命令: ssh-keygen -t rsa -b 4096 -C "your_emailexample.com" 参数说明: -t rsa&#x…...

Module Federation 和 Native Federation 的比较
前言 Module Federation 是 Webpack 5 引入的微前端架构方案,允许不同独立构建的应用在运行时动态共享模块。 Native Federation 是 Angular 官方基于 Module Federation 理念实现的专为 Angular 优化的微前端方案。 概念解析 Module Federation (模块联邦) Modul…...
css的定位(position)详解:相对定位 绝对定位 固定定位
在 CSS 中,元素的定位通过 position 属性控制,共有 5 种定位模式:static(静态定位)、relative(相对定位)、absolute(绝对定位)、fixed(固定定位)和…...
OpenLayers 分屏对比(地图联动)
注:当前使用的是 ol 5.3.0 版本,天地图使用的key请到天地图官网申请,并替换为自己的key 地图分屏对比在WebGIS开发中是很常见的功能,和卷帘图层不一样的是,分屏对比是在各个地图中添加相同或者不同的图层进行对比查看。…...

Redis数据倾斜问题解决
Redis 数据倾斜问题解析与解决方案 什么是 Redis 数据倾斜 Redis 数据倾斜指的是在 Redis 集群中,部分节点存储的数据量或访问量远高于其他节点,导致这些节点负载过高,影响整体性能。 数据倾斜的主要表现 部分节点内存使用率远高于其他节…...
全面解析各类VPN技术:GRE、IPsec、L2TP、SSL与MPLS VPN对比
目录 引言 VPN技术概述 GRE VPN 3.1 GRE封装结构 3.2 GRE的应用场景 GRE over IPsec 4.1 GRE over IPsec封装结构 4.2 为什么使用GRE over IPsec? IPsec VPN 5.1 IPsec传输模式(Transport Mode) 5.2 IPsec隧道模式(Tunne…...