欢迎光临散文网 会员登陆 & 注册

高斯白噪声信道下后验概率的计算

2022-08-03 00:19 作者:乐吧的数学  | 我要投稿

这篇小笔记,是对文章和视频《LDPC 低密度奇偶校验码的软判决译码算法浅析》的一个补充。在留言区有人问,后验概率为什么是已知的,所以,我就写了这个小短文。

录制的视频在:https://www.bilibili.com/video/BV1SV4y1T7pe/


这篇文章也可以独立看,就是如何分析高斯白噪声信道下的后验概率。



我们知道高斯白噪声符合以下概率密度(这是均值为 0 ,方差为 %5Csigma%5E2 的高斯分布):

p(n)%20%3D%20%5Cfrac%7B1%7D%7B%5Csqrt%7B2%5Cpi%7D%5Csigma%7D%0A%20%20%20%20e%5E%7B-%5Cfrac%7Bn%5E2%7D%7B2%5Csigma%5E2%7D%7D系统框图如下:


那么,发送 c 的情况下,收到的是 r 的概率是多少呢?

r = c + n

则:

n = r -c

则:
p(r-c)%20%3D%20%5Cfrac%7B1%7D%7B%5Csqrt%7B2%5Cpi%7D%5Csigma%7D%0A%20%20%20%20e%5E%7B-%5Cfrac%7B(r-c)%5E2%7D%7B2%5Csigma%5E2%7D%7D


c 是已知,确定的,所以,上面的式子中隐含的一个条件是 c 已知,在 c 已知并且是确定某个值的情况下, r - c 与 r 是一一对应的,则

p(r%7Cc)%3Dp(r-c%7Cc)%20%3D%5Cfrac%7B1%7D%7B%5Csqrt%7B2%5Cpi%7D%5Csigma%7D%0A%20%20%20%20e%5E%7B-%5Cfrac%7B(r-c)%5E2%7D%7B2%5Csigma%5E2%7D%7D


这是先验概率公式,即确定发送的数据,计算收到 r 的概率。



下面,需要推导一个后验概率:

发送的是 c,收到的是 r,则,我们想要知道一个概率,即“ 收到  r 的情况下,发送端发送的 c 是 1 的概率”:

p(c|r)

我们用条件概率和全概率公式,做一个推导

p(c%7Cr)%3D%20%5Cfrac%7Bp(c%2Cr)%7D%7Bp(r)%7D%20%3D%20%5Cfrac%7Bp(r%7Cc)p(c)%7D%7B%5Csum_c%20p(r%7Cc)p(c)%7D


其中 p(c) 表示发送的是 c 的概率,我们一般假定,发送的数据,是等概率出现的,例如我们发送的是 0  和 1 两种数据,则 0 和 1 是等概率出现的。

所以,根据上面的公式,就可以计算出后验概率。



补充:

在视频中,我讲解的时候,写的公式有点不准确,现在补充如下:

p(c_i%7Cr_i)%3D%20%5Cfrac%7Bp(c_i%2Cr_i)%7D%7Bp(r_i)%7D%20%3D%20%5Cfrac%7Bp(r_i%7Cc_i)p(c_i)%7D%7B%5Csum_%7Bc_i%7D%20p(r_i%7Cc_i)p(c_i)%7D


例如,我们收到的数据是  0.23. 我们想知道发送的数据是 1 的概率,假设噪声是符合均值为 0 方差为 1 的高斯分布,则:

p(c_i%3D1%7Cr_i%3D0.23)%3D%20%5Cfrac%7Bp(c_i%3D1%2Cr_i%3D0.23)%7D%7Bp(r_i%3D0.23)%7D%20%3D%20%5Cfrac%7Bp(r_i%3D0.23%7Cc_i%3D1)p(c_i%3D1)%7D%7B%5Csum_%7Bc_i%7D%20p(r_i%3D0.23%7Cc_i)p(c_i)%7D%20%5C%5C%0A%3D%5Cfrac%7Bp(r_i%3D0.23%7Cc_i%3D1)p(c_i%3D1)%7D%0A%20%20%7B%20p(r_i%3D0.23%7Cc_i%3D0)p(c_i%3D0)%2Bp(r_i%3D0.23%7Cc_i%3D1)p(c_i%3D1)%7D%20%5C%5C%0A%20%20%3D%5Cfrac%7Bp(r_i%3D0.23%7Cc_i%3D1)%7D%0A%20%20%7B%20p(r_i%3D0.23%7Cc_i%3D0)%2Bp(r_i%3D0.23%7Cc_i%3D1)%7D%20%5C%5C%0A%20%20%3D%5Cfrac%7B0.2966%7D%7B0.2966%2B0.3885%7D%20%5C%5C%0A%20%20%3D0.4329


同理,可以计算出来收到的数据是  0.23 的情况下. 我发送的数据是 0 的概率
p(c_i%3D0%7Cr_i%3D0.23)%20%3D%5Cfrac%7Bp(r_i%3D0.23%7Cc_i%3D0)%7D%0A%20%20%7B%20p(r_i%3D0.23%7Cc_i%3D0)%2Bp(r_i%3D0.23%7Cc_i%3D1)%7D%20%5C%5C%0A%20%20%3D%5Cfrac%7B0.3885%7D%7B0.2966%2B0.3885%7D%20%5C%5C%0A%20%20%3D0.5671

我们看到,上面的计算还是挺复杂的,需要计算 e 的指数,那么,就有把

%5Cfrac%7Bp(c%3D1%7Cr)%7D%7Bp(c%3D0%7Cr)%7D%3D%5Cfrac%20%7B%20%5Cfrac%7B1%7D%7B%5Csqrt%7B2%5Cpi%7D%5Csigma%7D%0A%20%20%20%20e%5E%7B-%5Cfrac%7B(r-1)%5E2%7D%7B2%5Csigma%5E2%7D%7D%20%7D%0A%7B%20%5Cfrac%7B1%7D%7B%5Csqrt%7B2%5Cpi%7D%5Csigma%7D%0A%20%20%20%20e%5E%7B-%5Cfrac%7B(r-0)%5E2%7D%7B2%5Csigma%5E2%7D%7D%20%7D%0A%20%20%20%20%3D%0A%20%20%20%20e%5E%7B%5Cfrac%7B(r-0)%5E2-(r-1)%5E2%7D%7B2%5Csigma%5E2%7D%7D%20%20%20%3D%20e%5E%7B%5Cfrac%7B2r-1%7D%7B2%5Csigma%5E2%7D%7D


如果这个值大于 1, 则可以判决 发送的数据 c 是 1,如果小于 1,则认为发送的数据 c 是 0.  如果对上面再去自然对数,则可以把指数的计算拿掉:


ln%20%5Cfrac%7Bp(c%3D1%7Cr)%7D%7Bp(c%3D0%7Cr)%7D%20%3D%20%5Cfrac%7B2r-1%7D%7B2%5Csigma%5E2%7D

高斯白噪声信道下后验概率的计算的评论 (共 条)

分享到微博请遵守国家法律