【小问题】距离估计和频率估计的方差下界推导出距离估计的方差下界
【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;若不存在ÿ…...
XML Group端口详解
在XML数据映射过程中,经常需要对数据进行分组聚合操作。例如,当处理包含多个物料明细的XML文件时,可能需要将相同物料号的明细归为一组,或对相同物料号的数量进行求和计算。传统实现方式通常需要编写脚本代码,增加了开…...
基于距离变化能量开销动态调整的WSN低功耗拓扑控制开销算法matlab仿真
目录 1.程序功能描述 2.测试软件版本以及运行结果展示 3.核心程序 4.算法仿真参数 5.算法理论概述 6.参考文献 7.完整程序 1.程序功能描述 通过动态调整节点通信的能量开销,平衡网络负载,延长WSN生命周期。具体通过建立基于距离的能量消耗模型&am…...
基于Flask实现的医疗保险欺诈识别监测模型
基于Flask实现的医疗保险欺诈识别监测模型 项目截图 项目简介 社会医疗保险是国家通过立法形式强制实施,由雇主和个人按一定比例缴纳保险费,建立社会医疗保险基金,支付雇员医疗费用的一种医疗保险制度, 它是促进社会文明和进步的…...
解锁数据库简洁之道:FastAPI与SQLModel实战指南
在构建现代Web应用程序时,与数据库的交互无疑是核心环节。虽然传统的数据库操作方式(如直接编写SQL语句与psycopg2交互)赋予了我们精细的控制权,但在面对日益复杂的业务逻辑和快速迭代的需求时,这种方式的开发效率和可…...
Vue2 第一节_Vue2上手_插值表达式{{}}_访问数据和修改数据_Vue开发者工具
文章目录 1.Vue2上手-如何创建一个Vue实例,进行初始化渲染2. 插值表达式{{}}3. 访问数据和修改数据4. vue响应式5. Vue开发者工具--方便调试 1.Vue2上手-如何创建一个Vue实例,进行初始化渲染 准备容器引包创建Vue实例 new Vue()指定配置项 ->渲染数据 准备一个容器,例如: …...
vue3 字体颜色设置的多种方式
在Vue 3中设置字体颜色可以通过多种方式实现,这取决于你是想在组件内部直接设置,还是在CSS/SCSS/LESS等样式文件中定义。以下是几种常见的方法: 1. 内联样式 你可以直接在模板中使用style绑定来设置字体颜色。 <template><div :s…...
Unit 1 深度强化学习简介
Deep RL Course ——Unit 1 Introduction 从理论和实践层面深入学习深度强化学习。学会使用知名的深度强化学习库,例如 Stable Baselines3、RL Baselines3 Zoo、Sample Factory 和 CleanRL。在独特的环境中训练智能体,比如 SnowballFight、Huggy the Do…...
Android15默认授权浮窗权限
我们经常有那种需求,客户需要定制的apk集成在ROM中,并且默认授予其【显示在其他应用的上层】权限,也就是我们常说的浮窗权限,那么我们就可以通过以下方法在wms、ams等系统服务的systemReady()方法中调用即可实现预置应用默认授权浮…...
CMake控制VS2022项目文件分组
我们可以通过 CMake 控制源文件的组织结构,使它们在 VS 解决方案资源管理器中以“组”(Filter)的形式进行分类展示。 🎯 目标 通过 CMake 脚本将 .cpp、.h 等源文件分组显示在 Visual Studio 2022 的解决方案资源管理器中。 ✅ 支持的方法汇总(共4种) 方法描述是否推荐…...
安卓基础(aar)
重新设置java21的环境,临时设置 $env:JAVA_HOME "D:\Android Studio\jbr" 查看当前环境变量 JAVA_HOME 的值 echo $env:JAVA_HOME 构建ARR文件 ./gradlew :private-lib:assembleRelease 目录是这样的: MyApp/ ├── app/ …...
