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

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

相关推荐
Aspiresky1 天前
浅析Linux进程信号处理机制:基本原理及应用
linux·运维·信号处理
太爱学习了4 天前
FPGA雷达信号处理之:自适应门限阈值
fpga开发·信号处理
达不溜的日记5 天前
ADC模数转换器详解(基于STM32)
stm32·单片机·嵌入式硬件·信息与通信·信号处理
极度畅想7 天前
脑电数据预处理十六:自动化阈值检测原理与实践
信号处理·脑电数据处理·eeg预处理·伪迹去除·脑电预处理·阈值检测·数据剔除
通信小呆呆7 天前
5G/6G时代的智能超表面:如何重构无线传播环境?
5g·重构·信息与通信·信号处理·超表面
码农hbk7 天前
Linux signal 图文详解(三)信号处理
linux·信号处理
山河君9 天前
webrtc之高通滤波——HighPassFilter源码及原理分析
算法·音视频·webrtc·信号处理
山河君10 天前
webrtc之语音活动上——VAD能量检测原理以及源码详解
算法·音视频·webrtc·信号处理
迎风打盹儿11 天前
均匀圆形阵抗干扰MATLAB仿真实录与特点解读
matlab·信号处理·抗干扰·均匀圆阵·波束合成
Ro Jace11 天前
科研笔记:数学建模启发的课题研究方法
数学建模·信号处理