【小问题】距离估计和频率估计的方差下界推导出距离估计的方差下界
【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;若不存在ÿ…...
AI编程协作实战:从提示工程到全周期开发工作流
1. 项目概述:从零开始的AI编程启蒙最近在GitHub上看到一个挺有意思的项目,叫jnMetaCode/ai-coding-101。光看名字,你大概能猜到它和AI编程有关,但“101”这个后缀又暗示着它的入门属性。作为一个在代码堆里摸爬滚打多年的老码农&a…...
ATF IronPython集成:如何在C应用中嵌入Python脚本引擎的完整指南
ATF IronPython集成:如何在C#应用中嵌入Python脚本引擎的完整指南 【免费下载链接】ATF Authoring Tools Framework (ATF) is a set of C#/.NET components for making tools on Windows. ATF has been in continuous development in Sony Computer Entertainments …...
基于MCP协议的智能文档处理工具simdoc-mcp:从RAG原理到Claude集成实战
1. 项目概述:从“文档理解”到“智能交互”的范式跃迁最近在折腾一个挺有意思的开源项目,叫simdoc-mcp。乍一看这个名字,可能有点摸不着头脑,svd-ai-lab是背后的团队,simdoc是核心,mcp是关键协议。简单来说…...
企业如何通过Taotoken实现API密钥的统一管理与审计
🚀 告别海外账号与网络限制!稳定直连全球优质大模型,限时半价接入中。 👉 点击领取海量免费额度 企业如何通过Taotoken实现API密钥的统一管理与审计 在将大模型能力集成到企业业务流程的过程中,一个常见的挑战是如何安…...
故障诊断创新算法之【先验知识+协同学习】基于故障特征掩码引导和潜在特征拆分的自编码器机械故障诊断(PyTorch)
小样本条件下,纯数据驱动方法很容易陷入过拟合和特征盲目提取,所以提出一种物理引导的深度诊断范式:将轴承内圈、外圈、滚动体的故障特征频率先验显式编码为故障特征掩码,并引入Huber函数构建先验引导损失,迫使网络学习…...
【日常小问】解决 Jenkins 部署 Spring Cloud 微服务到 Docker 容器启动失败的问题
一、问题出现在使用 Jenkins 进行 CI/CD 部署 Spring Cloud 微服务项目时,遇到了一个让人头疼的问题:所有通过 Jenkins 构建的 Docker 容器启动后立即退出,状态码为 Exited (1)。查看容器日志,报错信息如下:**********…...
DNS 服务器学习笔记:核心总结与实验指南
DNS 服务器学习笔记:核心总结与实验指南 📌 一、文章核心重点总结 1. DNS 基础知识 什么是 DNS? DNS(Domain Name System,域名系统)是互联网的“电话簿”,负责将人类易记的域名(如 w…...
JAVA学习之JAVASE基础
集合列表ListArrayList利用空参创建的集合,在底层创建一个默认长度为0的数组添加第一个元素时,底层会创建一个新的长度为10的数组存满时,会扩容1.5倍一次存多个元素,1.5倍还不够,则新创建的数组长度以实际为准LinkedLi…...
GitHub Actions 工作流中的输出处理
在现代软件开发中,CI/CD(持续集成和持续交付)是确保代码质量和自动化部署的关键环节。GitHub Actions 作为 GitHub 提供的 CI/CD 工具,支持通过工作流文件定义自动化任务。本文将结合一个实际的 GitHub Actions 工作流实例,探讨如何处理 Python 脚本的输出,并根据该输出决…...
AI助力船舶稳性计算:Gemini3.1Pro设计辅助新思路
在船舶设计工作中,稳性计算一直是非常核心、也非常严谨的环节。无论是新船方案设计、改装评估,还是载况校核,都需要围绕重心、浮心、横稳心、复原力臂、装载状态、自由液面影响等内容进行系统分析。过去这些资料往往分散在规范条文、设计手册…...
