参数估计理论

估计理论的主要任务是在某种信号假设下,估算该信号中某个参数(比如幅度、相位、达到时间)的具体取值。
参数估计:先假定研究的问题具有某种数学模型, 如正态分布,二项分布,再用已知类别的学习样本估计里面的参数。例如:ARMA模型功率谱估计。
非参数估计:不假定数学模型,用已知类别的学习样本的先验知识直接估计数学模型。例如:经典
功率谱估计。

一、估计量的性能

1.1 估计的基本概念

由随机信号的N个样本{xi (n) i=1,2,...,N}获得的真实参数θ1 , θ2 ,… θp的估计是一个将N维样本空间X^{N}映射为p维参数空间Θ的函数T,记作T: X^{N} → Θ。

估计偏差:

如果估计偏差等于零,则称为无偏估计。
渐进无偏估计:

无偏估计一定是渐进无偏的,反之不一定成立。

1.2 常用的估计量

(1)均值估计
(2)方差估计
(3)自相关估计
无偏估计
渐进无偏估计
(4)互相关估计

二、Fisher信息与Cramer-Rao下界

2.1 品质函数

在真实参数θ已知的条件下,样本x获得估计量是否最优,可由品质函数来评估。
品质函数:

2.2 fisher信息

品质函数V(x)的方差称为fisher信息。

2.3 Cramer-Rao下界

三、Bayes估计

3.1 代价函数

代价函数:衡量估计值与参数真实值之间误差的函数。
代价函数函数需要满足非负性和具有极小值两个条件。
常用代价函数:
(1)误差绝对值代价函数
(2)误差平方代价函数
(3)均匀代价函数

3.2 Bayes估计

假设被估计量θ的先验概率密度函数为p(θ),代价函数是随机变量θ和观测信号x的函数,因此平均
代价函数为

其中, p(x, θ)是θx的联合概率密度函数。
Bayes估计:在p(θ)已知,代价函数确定条件下, 使平均代价函数C最小化的参数估计。
        贝叶斯估计是把待估的参数θ作为具有某种先验分布p(θ)的随机变量,通过随机样本xθ的分布进行修正,由样本x进行修正后的概率密度p(θ|x)称为后验概率。贝叶斯估计的实质是利用后验概率p(θ|x)对θ进行推断。

3.3 估计步骤

① 确定θ的先验概率密度p(θ)。
② 确定样本x的条件概率密度p(x|θ),它是θ的函数。
③ 利用贝叶斯公式,求θ的后验概率

④ 计算参数θ的贝叶斯估计

四、最大似然估计

最大似然估计常用来估计未知的非随机参量或者概率密度函数未知的随机参量。
        假设被估计的参数θ为未知常数,给定样本x1 , x2 ,…. xN,则将样本的联合概率密度p(x1 , x2 ,…. xN |θ)称为似然函数。为便于计算,似然函数通常写为对数形式:ln p(x1 , x2 ,…. xN |θ)
参数θ的最大似然估计在似然函数达到极大值时求得。

为得到θ的最大似然估计,需对似然函数求导,并令导数等于0。

如果N个样本相互独立,而参数θ为向量θ=(θ1 , θ2 ,…,θp)^T ,则p个未知参量可以通过下式求解。

五、线性均方估计

        线性均方估计的规则,就是把估计量构造成观测量的线性函数,同时要求估计量的均方误差最小。
线性估计:
已知观测样本为xi (i=1,…,N),则参数θ的估计值可以写为

估计量的均方误差为
线性均方估计通过选择最佳系数aib,使得估计量的均方误差最小。
均方误差分别对aib求偏导,并令结果等于0

整理可得
估计的均值
估计的均方误差

六、最小二乘估计

        Bayes估计、线性均方估计需要知道被估计量的先验概率密度; 最大似然估计需要知道似然函数;最小二乘估计不需要先验统计特性,适用范围更广。
未知参数向量为:
观测信号模型为:
其中
观测向量
噪声向量
观测矩阵
最小二乘估计的代价函数为:
所以,最小二乘估计量,是满足下述方程的解:
由于
所以
若观测噪声n的均值矢量为零,则最小二乘估计的均方误差为

七、加权最小二乘估计

通过加权的方式,提高估计的精确度—— 加权最小二乘估计,其代价函数为:
W称为加权矩阵,是N×N对称正定阵。
加权最小二乘估计量满足下述方程:
所以
加权最小二乘估计矢量的均方误差为:
若观测噪声矢量n的均值矢量为零,协方差矩阵为Cn,则最优的加权矩阵为:
估计量为:
此时,加权最小二乘估计矢量的均方误差为:

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.rhkb.cn/news/468440.html

如若内容造成侵权/违法违规/事实不符,请联系长河编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

[vulnhub] DarkHole: 2

https://www.vulnhub.com/entry/darkhole-2,740/ 端口扫描主机发现 探测存活主机,185是靶机 # nmap -sP 192.168.75.0/24 Starting Nmap 7.94SVN ( https://nmap.org ) at 2024-11-08 18:02 CST Nmap scan report for 192.168.75.1 Host is up (0.…

【温度表达转化】

【温度表达转化】 C语言代码C代码Java代码Python代码 💐The Begin💐点点关注,收藏不迷路💐 利用公式 C5∗(F−32)/9 (其中C表示摄氏温度,F表示华氏温度) 进行计算转化。 输出 输出一行&#x…

【Promise】JS 异步之宏队列与微队列

文章目录 1 原理图2 说明3 相关面试题3.1 面试题13.2 面试题23.3 面试题33.4 面试题4 1 原理图 2 说明 JS 中用来存储待执行回调函数的队列包含 2 个不同特定的队列:宏队列和微队列。宏队列:用来保存待执行的宏任务(回调),比如:定…

【Linux】Linux入门实操——vim、目录结构、远程登录、重启注销

一、Linux 概述 1. 应用领域 服务器领域 linux在服务器领域是最强的,因为它免费、开源、稳定。 嵌入式领域 它的内核最小可以达到几百KB, 可根据需求对软件剪裁,近些年在嵌入式领域得到了很大的应用。 主要应用:机顶盒、数字电视、网络…

ubuntu下aarch64-linux-gnu(交叉编译) gdb/gdbserver(二)

ubuntu下aarch64-linux-gnu(交叉编译) gdb/gdbserver(二) 本教程作为gdb/gdbserver编译安装教程的一个补充,教会大家如何使用gdb/gdbserver进行远程调试。 如上图所示,我们需要将编译后的gdbserver上传至目标设备,其上…

Flutter错误: uses-sdk:minSdkVersion 16 cannot be smaller than version 21 declared

前言 今天要做蓝牙通信的功能,我使用了flutter_reactive_ble这个库,但是在运行的时候发现一下错误 Launching lib/main.dart on AQM AL10 in debug mode... /Users/macbook/Desktop/test/flutter/my_app/android/app/src/debug/AndroidManifest.xml Err…

c中柔性数组

c99中,结构中最后一个元素允许是未知大小的数组,这就叫柔性数组成员。 柔性数组的特点 1.结构中柔性数组前必须至少有一个其他成员 2.sizeof返回的这种结构大小不包括柔性数组的内存 3.包含柔性数组成员的结构用malloc函数进行动态分配,并…

WPS 默认模板修改

重装系统把word自定义样式搞没了,安装office时间太长,转战wps 解决方案 打开wps 点击【新建】word空白文档 设置修改你自己的样式 点击文件–另存为–Microsoft Word 带宏的模板文件(*.dotm) 另存路径为如下: 查…

Ubuntu24.04网络异常与应对方案记录

PS: 参加过408改卷的ZJU ghsongzju.edu.cn 开启嘲讽: 你们知道408有多简单吗,操作系统真实水平自己知道就行~~ Requested credits of master in UWSC30,in ZJU24,domestic master is too simple ubuntu安全软件 在 U…

[C++11] Lambda 表达式

lambda 表达式(Lambda Expressions)作为一种匿名函数,为开发者提供了简洁、灵活的函数定义方式。相比传统的函数指针和仿函数,lambda 表达式在简化代码结构、提升代码可读性和编程效率方面表现出色。 Lambda 表达式的基本语法 在…

Docker平台搭建方法

Docker平台搭建方法 1.1在VMware中创建两个虚拟机,只需要1个网卡,连接192.168.200.0网络。 虚拟机分配2个CPU,2G内存,60G硬盘,主机名分别为server和client,IP地址分别为192.168.200.137和192.168.200.138。server节点还兼做regis…

【学习笔记】Kylin-Desktop-V10-SP1 麒麟系统知识4——设备设置

提示:学习麒麟Kylin-Desktop-V10-SP1系统设备设置相关知识,包含设备设置进入方法、配置打印机、设置鼠标、键盘相关参数(包含输入法的配置)、以及管理快捷键组合、和多屏协同相关配置 一、前期准备 成功安装麒麟系统&#xff08…

Linux应用项目之量产工具(一)——显示系统

目录 前言 项目特点及介绍 ① 简单易用 ② 软件可配置、易扩展 ③ 纯 C 语言编程 软件总框架 显示系统 1.数据结构抽象 disp_manager.h 2.Framebuffer编程 framebuffer.c 3.显示管理 disp_manager.c 4.单元测试 disp_test.c 顶层目录Makefile 顶层目录Makefil…

企微SCRM价格解析及其性价比分析

内容概要 在如今的数字化时代,企业对于客户关系管理的需求日益增长,而企微SCRM(Social Customer Relationship Management)作为一款新兴的客户管理工具,正好满足了这一需求。本文旨在为大家深入解析企微SCRM的价格体系…

leetcode92:反转链表||

给你单链表的头指针 head 和两个整数 left 和 right &#xff0c;其中 left < right 。请你反转从位置 left 到位置 right 的链表节点&#xff0c;返回 反转后的链表 。 示例 1&#xff1a; 输入&#xff1a;head [1,2,3,4,5], left 2, right 4 输出&#xff1a;[1,4,3,2…

Python——数列1/2,2/3,3/4,···,n/(n+1)···的一般项为Xn=n/(n+1),当n—>∞时,判断数列{Xn}是否收敛

没注释的源代码 from sympy import * n symbols(n) s n/(n1) print(数列的极限为&#xff1a;,limit(s,n,oo))

多线程的创建方式以及及Thread类详解

目录 一.线程的创建方法&#xff1a;&#xff08;重点&#xff09; 一&#xff1a;继承Thread类 写法一&#xff1a;正常写法 写法二&#xff1a;匿名内部类 二.实现Runnable接口 写法一&#xff1a;正常写法 写法二&#xff1a;匿名内部类 三. 实现 Callable 接口 ​…

成功解决WSL2上的Ubuntu22.04执行sudo apt-get update指令报错问题

问题&#xff1a;输入sudo apt-get update指令会显示如下报错 问题所在&#xff1a;Temporary failure in name resolution 显然是系统无法解析域名。这可能是 DNS 配置问题。 解决方案&#xff1a; 临时修改 DNS 配置 尝试手动修改 /etc/resolv.conf 文件来使用公共 DNS 服务…

【英特尔IA-32架构软件开发者开发手册第3卷:系统编程指南】2001年版翻译,2-19

文件下载与邀请翻译者 学习英特尔开发手册&#xff0c;最好手里这个手册文件。原版是PDF文件。点击下方链接了解下载方法。 讲解下载英特尔开发手册的文章 翻译英特尔开发手册&#xff0c;会是一件耗时费力的工作。如果有愿意和我一起来做这件事的&#xff0c;那么&#xff…

ffmpeg 视频滤镜:屏蔽边框杂色- fillborders

滤镜描述 fillborders 官网链接 > FFmpeg Filters Documentation fillborders滤镜有几种方式帮你屏蔽边框的杂色、不好的图案。 滤镜使用 参数 left <int> ..FV.....T. set the left fill border (from 0 to INT_MAX) (default 0)right …