统计信号处理基础 习题解答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估计量。

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

相关推荐
stbomei1 天前
基于 MATLAB 的信号处理实战:滤波、傅里叶变换与频谱分析
算法·matlab·信号处理
khystal2 天前
HUMS 2023齿轮箱数据分析
数据库·数据分析·信号处理
SKYDROID云卓小助手11 天前
三轴云台之控制信号解析与执行
运维·服务器·网络·人工智能·信号处理
khystal12 天前
ISTA为什么要加上软阈值激活函数?r若没有L1 正则化也要加其他激活函数吗?
神经网络·信号处理
2401_8238682212 天前
织构表面MATLAB仿真
人工智能·机器学习·matlab·信号处理
霖0012 天前
高级项目——基于FPGA的串行FIR滤波器
人工智能·经验分享·matlab·fpga开发·信息与通信·信号处理
刘小小_算法工程师13 天前
「ECG信号处理——(23)基于ECG和PPG信号的血压预测」2025年8月12日
信号处理
鸭鸭鸭进京赶烤14 天前
EI检索-学术会议 | 人工智能、虚拟现实、可视化
人工智能·物联网·5g·信息可视化·云计算·vr·信号处理
青草地溪水旁17 天前
Linux 信号处理标志sa_flags详解
linux·信号处理
青草地溪水旁17 天前
如何理解SA_RESTART”被信号中断的系统调用自动重启“?
linux·信号处理·sa_restart