统计信号处理-简单看看克拉美罗界

  各种研究领域(包括无线定位方向)都会碰到参数估计的问题,这时常常会看到克拉美罗界 (Cramér–Rao bound)  这个东西。很多随机信号的书都会介绍什么是克拉美罗界,但初学者学起来往往很吃力,本文从直观上简单讨论一下克拉美罗界的各个方面。

什么是参数估计问题

  假设一种最简单的情况:

    一个物理量为,我们使用某种方式去观测它,观测值为,由于存在噪声,此时为高斯噪声,

这种情况下,我们自然会直接使用观测值去估计,这时就会存在估计的误差,直观地理解,噪声的方差越大,估计就可能越不准确。

为什么要讨论克拉美罗界

  讨论克拉美罗界就是为了使用这个标准来衡量无偏估计量的性能。

  采用上面的方式,使用去估计,这个估计值会在真实值附近波动(看作随机变量)。我们需要使用一些标准来衡量这种估计的好坏,一个标准是估计值的平均,这里的这个估计量是无偏估计量。另一标准是这个估计值波动的剧烈程度,也就是方差。上面这个问题中,克拉美罗界就等于这个方差。

  可是为什么不直接讨论方差而要去计算克拉美罗界呢,因为方差是针对某一种特定的估计量(或者理解为估计方式)而言的,在上面的例子中,方差是估计量的方差()。对于稍微复杂一点点的问题,对的可以有各种不同的估计量,它们分别的方差是不同的。显然,对于无偏估计量而言,方差越小的估计方式性能越好,但是这个方差有一个下界,就是我们的克拉美罗界。

直观地理解克拉美罗界

  克拉美罗界本身不关心具体的估计方式,只是去反映:利用已有信息所能估计参数的最好效果。

  还是上面那个参数估计问题,当我们观察到的时候,我们可以知道真实值的概率密度分布是以为均值,为方差的正态分布,即:

                  

上图给出了两个似然函数的例子,直观地看,似然函数的“尖锐”性决定了我们估计位置参数的精度。这个“尖锐”性可以用对数似然函数峰值处的负的二阶导数来度量,即对数似然函数的曲率(对数似然函数就是在似然函数的基础山加一个自然对数,这样有利于计算)。计算过程我就不写了,有兴趣的可以自己算算,算完之后结果为:,是噪声的方差的倒数,也就是噪声越小,对数似然函数越尖锐。

  所以,可以这样理解,似然函数的“尖锐”程度,或者,符合似然函数分布的这组数据的方差,就是克拉美罗界。

不同的估计量(估计方式)是什么意思

  让我们来分析一个稍微复杂一点点的参数估计问题:

    一个物理量为,我们使用某种方式去观测它,观测值为,这是两个不同时刻的观测结果,一样的高斯噪声

  这种情况下,我们要估计,正常人可能会采用估计量,即前后两个观测的平均,也有人可能觉得这样计算量有点大,于是总是直接使用去估计,也有人觉得第二个观测值可能会受到系统影响而不准确,他更相信前面的观察值,于是总采取这样的估计量。这三个估计量都是无偏的:

  估计量的方差为:

  估计量的方差为:

  估计量的方差为:

  比较上面的三种估计量,第一种的方差最小,它的估计效果较好。实际上,如果第二个观测值真的不太准确,也就是后一个高斯噪声较大,这样的话也许第二个估计量就比较合适了。

  因此,不同的考虑方式可以产生各种不同的估计算法,这些不同的估计量都是在真实值附近波动的随机变量(有的有偏,有的无偏),它们分别的方差也是不一样的,但是数学家们证明了:任何无偏估计量的方差必定大于等于克拉美罗界。

克拉美罗界的基本计算

  我们假设这两次观察互相独立,仅受相同的高斯白噪声影响,那么根据已有的信息,真实值的似然函数为两个正态的概率密度分布相乘:(注意:pdf实际上应该再进行归一化处理,但是我们之后使用对数似然函数,乘不乘归一化系数都无所谓,对数之后变成了常数,求导的时候就没了)

与之前一样,可以计算出对数似然函数的二阶导数,得到结果为:。实际上,当观测数目为的时候,这个值将会是。也就是说,使用多个观测值的信息时,对数似然函数越“尖锐”。这个二阶导数(曲率)更一般的度量是(下面用来表示要估计的参数):

它度量了对数似然函数的平均曲率(很多情况下曲率与的值有关,取数学期望使得它仅为的函数),被称为数据的Fisher信息,直观地理解,信息越多,下限越低,它具有信息测度的基本性质(非负的、独立观测的可加性)。一般来说,Fisher信息的倒数就是克拉美罗界了,任何无偏估计量的方差满足:

大多情况下,这个不等式的右边(克拉美罗界)是的函数。

克拉美罗界的标准定义

  (定理:Cramer-Rao下限----标量参数)

  假定PDF 满足“正则”条件(对于所有的):

其中数学期望是对 求取的。那么,任何无偏估计量的方差必定满足:

其中导数是在的真值处计算的,数学期望是对求取的。而且,对于某个函数,当且仅当

时,对所有达到下限的无偏估计量就可以求得。这个估计量是,它是MVU估计量(最小方差无偏估计),最小方差是

总结

  估计一个参数,根据已有信息得到了似然函数(或者pdf),这个pdf的“尖锐”性,或者,符合似然函数分布的这组数据的方差,就是克拉美罗界,它可以通过对对数似然函数求二阶导再取倒数得到。克拉美罗界的计算不依赖具体的估计方式,它可以用来作为一个衡量估计方式好坏的标准,及估计量的方差越靠近克拉美罗界,效果越好。

(注:本文主要参考《统计信号处理基础-估计与检测理论》-国外电子与通信教材系列)

作者: rubbninja
出处: http://www.cnblogs.com/rubbninja/
关于作者:目前是学生,主要研究领域为机器学习与无线定位技术,欢迎讨论与指正!
版权声明:本文版权归作者和博客园共有,转载请注明出处。

时间: 2024-08-06 09:40:32

统计信号处理-简单看看克拉美罗界的相关文章

统计难题(简单字典树)

字典树(讲解+模板) 又称单词查找树,Trie树,是一种树形结构,是一种哈希树的变种.典型应用是用于统计,排序和保存大量的字符串(但不仅限于字符串),所以经常被搜索引擎系统用于文本词频统计.它的优点是:利用字符串的公共前缀来节约存储空间,最大限度地减少无谓的字符串比较,查询效率比哈希表高. 字典树跟字典很相似,当我们要查询一个单词是不是在字典中时,我们首先查询单词的第一个字母,之后确定了一个字母后,查找下一个字母,反复这种操作,直到找到单词,但是这种效率很低,在庞大的单词系统面前,则显得心塞了.

统计语言模型——《数学之美》

华为训练题:初级——输入一行字符,分别统计(简单题)

下面是AC代码 1 #include<iostream> 2 #include<string> 3 #include<vector> 4 using namespace std; 5 6 7 vector<int> jishu(string& str) 8 { 9 int len=str.size(); 10 vector<int> vec; 11 int re1=0,re2=0,re3=0,re4=0; 12 for(int i=0;i

简单爬虫一部美剧(二:下载进度条)

上篇获取到了每一集的下载url,这篇就用获取到的url下载视频 在下载文件时希望可以给出进度条之类的提示,在网上搜索了一波,发现有一个库可以实现:tqdm库,具体用法参考这篇博客:https://www.jianshu.com/p/1ed2a8b2c77b 在原来的类下面,再加一个方法,用来下载文件,如下 def download_file(url, name): """下载文件""" try: response = requests.get(u

估计理论—从经典到贝叶斯

本文内容主要参考Steven M.Kay的<统计信号处理基础——估计与检测理论>,该书中译本分类为“国外电子与通信教材系列”,应该会有一定局限性.本文是我看过该书后的一点点总结. 1.从最大似然估计看经典估计理论 最大似然估计(Maximum Likelihood Estimation,MLE)是一种很经典的估计方法.顾名思义,就是最大化似然函数的一种估计方法.似然函数(Likelihood function),即(其中$\theta$代表估计量,$X$代表观察值). 对于不同的问题,似然函数

Cramer-Rao Bounds (CRB)

克拉美-罗界.又称Cramer-Rao lower bounds(CRLB),克拉美-罗下界. 克拉美罗界是对于参数估计问题提出的,为任何无偏估计量的方差确定了一个下限.无偏估计量的方差只能无限制的逼近CRB,而不会低于CRB,因此这个界也可以称为CRLB,意为克拉美罗下界. CRLB可以用于计算无偏估计中能够获得的最佳估计精度,因此经常用于计算理论能达到的最佳估计精度,和评估参数估计方法的性能(是否接近CRLB下界). 克拉美罗界本身不关心具体的估计方式,只是去反映:利用已有信息所能估计参数的

[翻译和注解]Kafka Streams简介: 让流处理变得更简单

Introducing Kafka Streams: Stream Processing Made Simple 这是Jay Kreps在三月写的一篇文章,用来介绍Kafka Streams.当时Kafka Streams还没有正式发布,所以具体的API和功能和0.10.0.0版(2016年6月发布)有所区别.但是Jay Krpes在这简文章里介绍了很多Kafka Streams在设计方面的考虑,还是很值得一看的. 以下的并不会完全按照原文翻译,因为那么搞太累了……这篇文件的确很长,而且Jay

MATLAB信号处理工具箱函数列表分类

现将MATLAB信号处理工具箱函数进行分组,便于记忆查询和长期回顾.(只解释基本用途,具体用法请在help目录下查询) Waveform Generation(波形产生) chairp: 产生扫频余弦函数: diric: 产生Dirichlet或周期sinc函数: gauspuls: 产生高斯调制地正弦曲线脉冲: pulstran: 产生一个脉冲序列: rectpuls: 产生一个非周期的抽样方波: sawtooth: 产生锯齿波或三角波: sinc: 产生sinc函数,即sin(πt)/πt;

美学心得(第一百九十五集) 罗国正

  美学心得(第一百九十五集) 罗国正 (2018年12月) 2847.美的起点是主体本身,是感官和大脑,核心是大脑,最终在大脑落实.客体可以是诱因,但有时主体可以不要主体外的事物而产生美感.主体的大脑可以将主体的意象,或主体的某部分作为审美对象,诱发出美感.弄清这个根本的问题,就弄明白了美学的根源和逻辑起点.不然的话,就会渐行渐远,离题万里. 2848.狄德罗说:"从单一的关系的感觉得来的美,往往小于从多种关系的感觉得来的美."其实,单一的关系的感觉有强的,有弱的.多种关系的感觉也有