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

题目

对于后验PDF

画出时的PDF。然后,对于同样的值,求 MMSE和 MAP 的估计量。


解答

根据题目条件,通过Matlab,画出两个函数的概率密度函数图像,其中的x设定为2.

复制代码
clc, clear, close all;
x=2;
e=0.5;
theat=-10:0.01:10;
p_theat=e/sqrt(2*pi)*exp(-0.5*(theat-x).^2)+(1-e)/sqrt(2*pi)*exp(-0.5*(theat+x).^2);
figure(1)
plot(theat,p_theat)
hold on
xline(x);
xline(-x)
hold off

e=0.75;
p_theat=e/sqrt(2*pi)*exp(-0.5*(theat-x).^2)+(1-e)/sqrt(2*pi)*exp(-0.5*(theat+x).^2);
figure(2)
plot(theat,p_theat)
hold on
xline(x);
xline(-x)
hold off

时,做出图案:

此时可以发现,是个偶函数,因此:

也就是此时的的MMSE估计量为0。

而对的观察可知,该函数存在两个极值点,所以的MMSE估计量不唯一:


计算MMSE:

上式的积分,分别是的数学期望,因此可以得到:

因此:当

此时代入也验证了上半题目的结论。

而观察图像的特性,可以发现该函数存在一个极值点,也就是的MAP估计量为:

经过这个例子分析,可以发现,MAP估计量不一定等于MMSE估计量。

书上举得联合分布是高斯,先验也是高斯情况下两者相等,只是特例。

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