统计信号处理基础 习题解答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
相关推荐
JINX的诅咒17 小时前
CORDIC算法:三角函数的硬件加速革命——从数学原理到FPGA实现的超高效计算方案
算法·数学建模·fpga开发·架构·信号处理·硬件加速器
西北大程序猿21 小时前
linux进程信号 ─── linux第27课
linux·运维·服务器·信号处理
山河君1 天前
音频进阶学习二十四——IIR滤波器设计方法
学习·算法·音视频·信号处理
vonchenchen11 天前
nara wpe去混响学习笔记
机器学习·音视频·音频·信息与通信·信号处理
全栈开发圈1 天前
新书速览|深入理解DSP:基于TMS320F28379D的开发与实践
人工智能·信号处理·dsp开发
山河君1 天前
音频进阶学习二十五——脉冲响应不变法实现低通滤波器
学习·算法·音视频·信号处理
卡尔曼的BD SLAMer1 天前
期刊 | 不收版面费与审稿费的电子通信类期刊
算法·信息与通信·信号处理
北京青翼科技1 天前
【PCIE711-214】基于PCIe总线架构的4路HD-SDI/3G-SDI视频图像模拟源
图像处理·人工智能·fpga开发·信号处理
m0_555762902 天前
信号处理中的窗
信号处理
信号小探4 天前
信号与系统(郑君里)第一章-绪论 1-23 课后习题解答
算法·信息与通信·信号处理·抽象代数·傅里叶分析