【小问题】距离估计和频率估计的方差下界推导出距离估计的方差下界
【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;若不存在ÿ…...
未来机器人的大脑:如何用神经网络模拟器实现更智能的决策?
编辑:陈萍萍的公主一点人工一点智能 未来机器人的大脑:如何用神经网络模拟器实现更智能的决策?RWM通过双自回归机制有效解决了复合误差、部分可观测性和随机动力学等关键挑战,在不依赖领域特定归纳偏见的条件下实现了卓越的预测准…...
java调用dll出现unsatisfiedLinkError以及JNA和JNI的区别
UnsatisfiedLinkError 在对接硬件设备中,我们会遇到使用 java 调用 dll文件 的情况,此时大概率出现UnsatisfiedLinkError链接错误,原因可能有如下几种 类名错误包名错误方法名参数错误使用 JNI 协议调用,结果 dll 未实现 JNI 协…...
linux arm系统烧录
1、打开瑞芯微程序 2、按住linux arm 的 recover按键 插入电源 3、当瑞芯微检测到有设备 4、松开recover按键 5、选择升级固件 6、点击固件选择本地刷机的linux arm 镜像 7、点击升级 (忘了有没有这步了 估计有) 刷机程序 和 镜像 就不提供了。要刷的时…...
从零实现STL哈希容器:unordered_map/unordered_set封装详解
本篇文章是对C学习的STL哈希容器自主实现部分的学习分享 希望也能为你带来些帮助~ 那咱们废话不多说,直接开始吧! 一、源码结构分析 1. SGISTL30实现剖析 // hash_set核心结构 template <class Value, class HashFcn, ...> class hash_set {ty…...
Spring AI 入门:Java 开发者的生成式 AI 实践之路
一、Spring AI 简介 在人工智能技术快速迭代的今天,Spring AI 作为 Spring 生态系统的新生力量,正在成为 Java 开发者拥抱生成式 AI 的最佳选择。该框架通过模块化设计实现了与主流 AI 服务(如 OpenAI、Anthropic)的无缝对接&…...
拉力测试cuda pytorch 把 4070显卡拉满
import torch import timedef stress_test_gpu(matrix_size16384, duration300):"""对GPU进行压力测试,通过持续的矩阵乘法来最大化GPU利用率参数:matrix_size: 矩阵维度大小,增大可提高计算复杂度duration: 测试持续时间(秒&…...
SpringCloudGateway 自定义局部过滤器
场景: 将所有请求转化为同一路径请求(方便穿网配置)在请求头内标识原来路径,然后在将请求分发给不同服务 AllToOneGatewayFilterFactory import lombok.Getter; import lombok.Setter; import lombok.extern.slf4j.Slf4j; impor…...
SQL慢可能是触发了ring buffer
简介 最近在进行 postgresql 性能排查的时候,发现 PG 在某一个时间并行执行的 SQL 变得特别慢。最后通过监控监观察到并行发起得时间 buffers_alloc 就急速上升,且低水位伴随在整个慢 SQL,一直是 buferIO 的等待事件,此时也没有其他会话的争抢。SQL 虽然不是高效 SQL ,但…...
Linux nano命令的基本使用
参考资料 GNU nanoを使いこなすnano基础 目录 一. 简介二. 文件打开2.1 普通方式打开文件2.2 只读方式打开文件 三. 文件查看3.1 打开文件时,显示行号3.2 翻页查看 四. 文件编辑4.1 Ctrl K 复制 和 Ctrl U 粘贴4.2 Alt/Esc U 撤回 五. 文件保存与退出5.1 Ctrl …...
FFmpeg:Windows系统小白安装及其使用
一、安装 1.访问官网 Download FFmpeg 2.点击版本目录 3.选择版本点击安装 注意这里选择的是【release buids】,注意左上角标题 例如我安装在目录 F:\FFmpeg 4.解压 5.添加环境变量 把你解压后的bin目录(即exe所在文件夹)加入系统变量…...
