统计信号处理基础 习题解答11-4

题目

观测到数据:, 假定未知参数A具有先验 PDF

其中,, 是方差为的WGN,且与A独立,求A的MAP估计量。


解答

根据题目条件,得到条件概率:

那么对于N个观察的独立数据,有:

因此,得到

根据函数形式,具有相同的极值点,因此对,存在:

令上式等于零,得到此时A的MAP估计量为:

实际过程中,如果,那么可以取到MAP估计量的值,如下图所示。

但如,概率密度如下图。

那么概率密度只能保留的部分,此时概率密度的最大值就发生在A=0处,因此,综合来看,A的MAP估计为:

上图仿真代码:

Matlab 复制代码
clc, clear, close all;
N=10;
mu=0.3;
sigma=1;
x=normrnd(mu,sigma,N,1)';
lamda=1;
A_Value=-3:0.01:3;
p=zeros(1,length(A_Value));
for index=1:length(A_Value)
    A=A_Value(index);
    p(1,index)=((2*pi*sigma)^(-N/2))*exp(-1/(2*sigma)*sum((x-A).^2))*lamda.*exp(-lamda*A);
end
figure(1)
plot(A_Value,p)
hold on
xline(mean(x)-sigma*lamda/N);
xline(0,"--");
hold off
相关推荐
极度畅想4 小时前
脑电数据预处理十六:自动化阈值检测原理与实践
信号处理·脑电数据处理·eeg预处理·伪迹去除·脑电预处理·阈值检测·数据剔除
通信小呆呆10 小时前
5G/6G时代的智能超表面:如何重构无线传播环境?
5g·重构·信息与通信·信号处理·超表面
码农hbk20 小时前
Linux signal 图文详解(三)信号处理
linux·信号处理
山河君2 天前
webrtc之高通滤波——HighPassFilter源码及原理分析
算法·音视频·webrtc·信号处理
山河君4 天前
webrtc之语音活动上——VAD能量检测原理以及源码详解
算法·音视频·webrtc·信号处理
迎风打盹儿4 天前
均匀圆形阵抗干扰MATLAB仿真实录与特点解读
matlab·信号处理·抗干扰·均匀圆阵·波束合成
Ro Jace5 天前
科研笔记:数学建模启发的课题研究方法
数学建模·信号处理
T.i.s6 天前
毫米波雷达信号处理步骤顺序
信号处理·毫米波雷达
huangyuchi.7 天前
【Linux系统】万字解析,进程间的信号
linux·服务器·信号处理·信号产生·linux信号·信号保存·操作系统如何运行
半桔8 天前
【Linux手册】Unix/Linux 信号:原理、触发与响应机制实战
linux·运维·unix·信号处理