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

卷积码的 BCJR 译码算法 (二)--计算 γ

2023-01-05 12:33 作者:乐吧的数学  | 我要投稿

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

前面文章已经推导出来下面这个状态转移的条件概率,这个概率是进一步计算发送比特后验概率的基础。

P(%5Cpsi_t%3Dp%2C%5Cpsi_%7Bt%2B1%7D%3Dq%7Cr)%20%20%3D%20%5Cfrac%7B1%7D%7Bp(r)%7D%20%5Ctimes%20%20%20%20p(%20%5Cpsi_t%3Dp%20%2C%20r_%7Br%3Ct%7D)%20%5Ctimes%20p(%5Cpsi_%7Bt%2B1%7D%3Dq%2C%20r_t%20%7C%20%20%5Cpsi_t%3Dp)%20%5Ctimes%20%20p(r_%7Br%3Et%7D%20%7C%20%5Cpsi_%7Bt%2B1%7D%3Dq%20)%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5Ctag%7B1%7D


为了后面表达方便,我们把 (1)  中等号右边的三个部分,分别记为:

%5Cbegin%7Baligned%7D%20%5Calpha_t(p)%20%26%3D%20%20p(%20%5Cpsi_t%3Dp%20%2C%20r_%7Br%3Ct%7D)%20%20%20%5C%5C%20%5Cgamma_t(p%2Cq)%20%26%3D%20%20p(%5Cpsi_%7Bt%2B1%7D%3Dq%2C%20r_t%20%7C%20%20%5Cpsi_t%3Dp)%20%5C%5C%20%5Cbeta_%7Bt%2B1%7D(q)%20%26%3D%20p(r_%7Br%3Et%7D%20%7C%20%5Cpsi_%7Bt%2B1%7D%3Dq%20)%20%5Cend%7Baligned%7D%20%20%20%20%20%20%5Ctag%7B2%7D


则公式 (1) 就可以简写为

P(%5Cpsi_t%3Dp%2C%5Cpsi_%7Bt%2B1%7D%3Dq%7Cr)%20%3D%20%5Calpha_t(p)%20%5Cgamma_t(p%2Cq)%20%5Cbeta_%7Bt%2B1%7D(q)



公式 (2) 中的 第二个,实际上是可以计算的,我们做一下推导:

%5Cbegin%7Baligned%7D%20%5Cgamma_t(p%2Cq)%20%26%3D%20%20p(%5Cpsi_%7Bt%2B1%7D%3Dq%2C%20r_t%20%7C%20%20%5Cpsi_t%3Dp)%20%5C%5C%20%26%3D%20p(%20r_t%20%7C%20%5Cpsi_%7Bt%2B1%7D%3Dq%2C%20%20%5Cpsi_t%3Dp)%20p(%5Cpsi_%7Bt%2B1%7D%3Dq%7C%5Cpsi_t%3Dp)%20%5Cend%7Baligned%7D%20%20%20%20%20%20%5Ctag%7B3%7D


其中

p(%5Cpsi_%7Bt%2B1%7D%3Dq%7C%5Cpsi_t%3Dp)%20%3D%20p(X_t%20%3D%20x_t%5E%7Bp%2Cq%7D)%20%20%20%20%5Ctag%7B4%7D


公式 (4) 的含义,就是在 t 时刻状态为 p 时, t+1 时刻转到状态 q 的概率,也就是在 t 时刻,输入的比特是让状态从 p 转到 q 的值,例如,在第 6 时刻,状态是 1,那么转到 第 7 时刻状态为 2 的概率,就是在 6 时刻输入比特是 0 的概率,因为输入 0 ,能让状态从 1 转到 2 .  一般是等概率的假定,所以,这个一般就是 1/2.

再来看公式 (3) 中的另外一块:

p(%20r_t%20%7C%20%5Cpsi_%7Bt%2B1%7D%3Dq%2C%20%20%5Cpsi_t%3Dp)%20%3D%20p(%20r_t%20%7C%20X_t%20%3D%20x_t%5E%7Bp%2Cq%7D%20)%20%20%3Dp(%20r_t%20%7C%20X_t%20%3D%20v_t%5E%7Bp%2Cq%7D%20)%20%3Dp(%20r_t%20%7C%20a%20)%20%20%5Ctag%7B5%7D


上面的含义,是从 t 时刻的状态 p 转到 t+1  时刻的状态 q 的条件下,收到是 r_t 的概率。"从 t 时刻的状态 p 转到 t+1  时刻的状态 q "  对应的是一个输入比特,我们记为 x_t%5E%7Bp%2Cq%7D, 此时,对应一个输出 v_t%5E%7Bp%2Cq%7D,经过调制后得到数据为 a .  数据 a 通过高斯高斯白噪声信道送出,多个数据依次通过信道送出:

r_t%5E%7B(i)%7D%20%3D%20a%5E%7B(i)%7D%20%2B%20n_t


各自都是符合均值为 0 方差为 %5Csigma%5E2%20 的高斯分布:

p(r_t%5E%7B(i)%7D%7Ca%5E%7B(i)%7D%20)%20%3D%20%5Cfrac%7B1%7D%7B%5Csqrt%7B2%5Cpi%7D%20%5Csigma%7D%20e%5E%7B-%5Cfrac%7B%20(r_t%5E%7B(i)%7D-a%5E%7B(i)%7D)%5E2%20%7D%7B2%5Csigma%5E2%7D%7D

 

多个数据,即卷积码一个比特输出产生多少个比特的输出,我们记为 Q 个,则公式 (5) 为:

p(%20r_t%20%7C%20a%20)%20%3D%20%20%5Cfrac%7B1%7D%7B%20(2%5Cpi%20%5Csigma%5E2)%5E%7BQ%2F2%7D%7D%20e%5E%7B-%5Cfrac%7B%20%5Csum_%7Bi%3D1%7D%5EQ%20(r_t%5E%7B(i)%7D-a%5E%7B(i)%7D)%5E2%20%7D%7B2%5Csigma%5E2%7D%7D


我们举个例子说明一下,例如 t=6 , 即时刻 6. 状态 从 1 转到 2,则我们计算的概率为:

%5Cbegin%7Baligned%7D%20%5Cgamma_6(1%2C2)%20%26%3D%20%20p(%5Cpsi_7%3D2%2C%20r_6%20%7C%20%20%5Cpsi_6%3D1)%20%5C%5C%20%26%3D%20p(%20r_6%20%7C%20%5Cpsi_7%3D2%2C%20%20%5Cpsi_6%3D1)%20p(%5Cpsi_7%3D2%7C%5Cpsi_6%3D1)%20%20%5C%5C%20%26%3D%20p(%20r_6%20%7C%20%5Cpsi_7%3D2%2C%20%20%5Cpsi_6%3D1)%20p(X_6%3D0)%20%20%5C%5C%20%26%3D%20p(%20r_6%20%7C%20%5Cpsi_7%3D2%2C%20%20%5Cpsi_6%3D1)%20%5Ctimes%20%5Cfrac%7B1%7D%7B2%7D%20%5C%5C%20%5Cend%7Baligned%7D%20%20%20%20%20%20%5Ctag%7B6%7D%0A


其中:

%5Cbegin%7Baligned%7D%20p(%20r_6%20%7C%20%5Cpsi_7%3D2%2C%20%20%5Cpsi_6%3D1)%20%26%3D%20p(r_6%5E%7B(0)%7D%2C%20r_6%5E%7B(1)%7D%20%7C%20X_6%20%3D%200)%20%5C%5C%20%26%3D%20p(r_6%5E%7B(0)%7D%2C%20r_6%5E%7B(1)%7D%20%7C%20v_6%5E%7B(0)%7D%20%3D%200%2C%20v_6%5E%7B(0)%7D%20%3D%201)%20%20%5Cquad%20(%E5%8D%B7%E7%A7%AF%E7%BC%96%E7%A0%81)%5C%5C%20%26%3D%20p(r_6%5E%7B(0)%7D%2C%20r_6%5E%7B(1)%7D%20%7C%20a_6%5E%7B(0)%7D%20%3D%20-1%2C%20a_6%5E%7B(0)%7D%20%3D%20%2B1)%20%20%5Cquad%20(%E8%B0%83%E5%88%B6)%5C%5C%20%26%3D%20p(r_6%5E%7B(0)%7D%7C%20a_6%5E%7B(0)%7D%20%3D%20-1)%20%5Ctimes%20p(%20r_6%5E%7B(1)%7D%20%7C%20a_6%5E%7B(1)%7D%20%3D%20%2B1)%20%20%5Cquad%20(%E7%9B%B8%E4%BA%92%E7%8B%AC%E7%AB%8B)%5C%5C%20%5C%5C%20%26%3D%20%20%5Cfrac%7B1%7D%7B%5Csqrt%7B2%5Cpi%7D%20%5Csigma%7D%20e%5E%7B-%5Cfrac%7B%20(r_6%5E%7B(0)%7D-(-1))%5E2%20%7D%7B2%5Csigma%5E2%7D%7D%20%20%20%5Ctimes%20%20%20%20%20%20%20%20%20%5Cfrac%7B1%7D%7B%5Csqrt%7B2%5Cpi%7D%20%5Csigma%7D%20e%5E%7B-%5Cfrac%7B%20(r_6%5E%7B(1)%7D-(%2B1))%5E2%20%7D%7B2%5Csigma%5E2%7D%7D%20%20%20%5C%5C%20%5C%5C%20%26%3D%20%5Cfrac%7B1%7D%7B2%5Cpi%20%5Csigma%5E2%7D%20%20e%5E%7B-%5Cfrac%7B%20(r_6%5E%7B(0)%7D-(-1))%5E2%20%2B(r_6%5E%7B(1)%7D-(%2B1))%5E2%20%7D%7B2%5Csigma%5E2%7D%7D%20%20%5C%5C%20%5C%5C%20%26%20%5Cpropto%20%20e%5E%7B-%5Cfrac%7B%20(r_6%5E%7B(0)%7D%2B1)%5E2%20%2B(r_6%5E%7B(1)%7D-1)%5E2%20%7D%7B2%5Csigma%5E2%7D%7D%20%5Cend%7Baligned%7D%20%20%20%20%20%20%5Ctag%7B7%7D


把 (7) 代回 (6) 即可计算。我们在计算过程中,把公式 (6) 中的 1/2 和公式 (7) 中的 %5Cfrac%7B1%7D%7B2%5Cpi%20%5Csigma%5E2%7D%20都忽略掉不参与计算,因为这些在计算过程中都是不变的,最后在计算概率归一化的时候,能把他们消除掉。


卷积码的 BCJR 译码算法 (二)--计算 γ的评论 (共 条)

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