发明名称 非对数域最大后验概率Turbo译码方法
摘要 本发明公布了一种数字通信差错控制并行级联卷积特博(Turbo)码的非对数域最大后验概率(MAP)译码方法。本译码方法使用单极性二进制信号{0,1}代替传统的双极性二进制信号{+1,-1}来实现状态转移概率测度的计算,从而得到一种非对数域的软输入软输出(SISO)迭代译码算法。
申请公布号 CN100454768C 申请公布日期 2009.01.21
申请号 CN200510021016.1 申请日期 2005.06.01
申请人 西南交通大学 发明人 赵宏宇;范平志
分类号 H03M13/29(2006.01);H03M13/23(2006.01);H03M13/00(2006.01);H04L1/00(2006.01) 主分类号 H03M13/29(2006.01)
代理机构 成都博通专利事务所 代理人 陈树明
主权项 1、一种非对数域最大后验概率Turbo译码方法,该方法采用两个完全相同的软输入软输出(SISO)译码流程对通信系统接收端解调器输出的信息比特实数观测值序列和校验比特实数观测值序列进行译码,每个软输入软输出(SISO)译码流程内部包括前向迭代、后向迭代和外信息计算三个运算步骤,其特征在于:(1)第j个(j=0,1)软输入软输出(SISO)译码流程有三个实数输入信号序列,分别是指数信息比特实数观测值<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>x</mi><mi>k</mi><mo>/</mo></msubsup><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><msub><mi>L</mi><mi>c</mi></msub><msub><mi>x</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>指数校验比特实数观测值<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>y</mi><mrow><mi>j</mi><mo>,</mo><mi>k</mi></mrow><mo>/</mo></msubsup><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><msub><mi>L</mi><mi>c</mi></msub><msub><mi>y</mi><mrow><mi>j</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>)</mo></mrow></mrow></math>]]></maths>和信息比特先验概率似然比l<sub>(j+1)mod2,k</sub>,其中,x<sub>k</sub>是通信系统接收端解调器输出的第k个信息比特X<sub>k</sub>的实数观测值,y<sub>j,k</sub>是对应于第j(j=0,1)个系统递归卷积(RSC)分量编码器输入信息比特X<sub>k</sub>时输出的校验比特Y<sub>j,k</sub>的实数观测值,L<sub>c</sub>为信道可靠性测度,L<sub>c</sub>=4/N<sub>0</sub>,N<sub>0</sub>为高斯白噪声(AWGN)的单边功率谱密度;(2)第j个(j=0,1)软输入软输出(SISO)译码流程有两个实数输出信号序列,分别是送到另一个软输入软输出(SISO)译码流程的先验概率似然比1j,k和用作译码判决的后验概率似然比Λ(X<sub>k</sub>),其中,<maths num="0003"><![CDATA[<math><mrow><mi>&Lambda;</mi><mrow><mo>(</mo><msub><mi>X</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>=</mo><msubsup><mi>x</mi><mi>k</mi><mo>/</mo></msubsup><msub><mi>l</mi><mrow><mn>0</mn><mo>,</mo><mi>k</mi></mrow></msub><msub><mi>l</mi><mrow><mn>1</mn><mo>,</mo><mi>k</mi></mrow></msub><mo>,</mo></mrow></math>]]></maths>当Λ(X<sub>k</sub>)≥1时,将X<sub>k</sub>译为1,否则将X<sub>k</sub>译为0;(3)第j个(j=0,1)软输入软输出(SISO)译码流程输出的关于第k个信息比特X<sub>k</sub>的先验概率似然比为<maths num="0004"><![CDATA[<math><mrow><msub><mi>l</mi><mrow><mi>j</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>=</mo><mfrac><mrow><munder><mi>&Sigma;</mi><msup><mi>s</mi><mo>&prime;</mo></msup></munder><msub><mi>&alpha;</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><msup><mi>s</mi><mo>&prime;</mo></msup><mo>)</mo></mrow><msubsup><mi>&gamma;</mi><mi>k</mi><mrow><mo>(</mo><mi>e</mi><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><msubsup><mi>Y</mi><msup><mi>s</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><msub><mi>&beta;</mi><mi>k</mi></msub><mrow><mo>(</mo><msubsup><mi>t</mi><msup><mi>s</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mo>)</mo></mrow></mrow><mrow><munder><mi>&Sigma;</mi><msup><mi>s</mi><mo>&prime;</mo></msup></munder><msub><mi>&alpha;</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><msup><mi>s</mi><mo>&prime;</mo></msup><mo>)</mo></mrow><msubsup><mi>&gamma;</mi><mi>k</mi><mrow><mo>(</mo><mi>e</mi><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><msubsup><mi>Y</mi><msup><mi>s</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><msub><mi>&beta;</mi><mi>k</mi></msub><mrow><mo>(</mo><msubsup><mi>t</mi><msup><mi>s</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>其中,当<maths num="0005"><![CDATA[<math><mrow><msubsup><mi>Y</mi><msup><mi>s</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>=</mo><mn>0</mn><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>0,1</mn><mo>)</mo></mrow></mrow></math>]]></maths>时,信息比特X<sub>k</sub>的先验概率似然比的状态转移概率测度<maths num="0006"><![CDATA[<math><mrow><msubsup><mi>&gamma;</mi><mi>k</mi><mrow><mo>(</mo><mi>e</mi><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><msubsup><mi>Y</mi><msup><mi>s</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo>=</mo><mn>1</mn><mo>,</mo></mrow></math>]]></maths>否则,当<maths num="0007"><![CDATA[<math><mrow><msubsup><mi>Y</mi><msup><mi>s</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>=</mo><mn>1</mn></mrow></math>]]></maths>时,<maths num="0008"><![CDATA[<math><mrow><msubsup><mi>&gamma;</mi><mi>k</mi><mrow><mo>(</mo><mi>e</mi><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><msubsup><mi>Y</mi><msup><mi>s</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo>=</mo><msubsup><mi>y</mi><mrow><mi>j</mi><mo>,</mo><mi>k</mi></mrow><mo>/</mo></msubsup><mo>;</mo></mrow></math>]]></maths>而前向迭代计算公式为<maths num="0009"><![CDATA[<math><mrow><msub><mi>&alpha;</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munder><mi>&Sigma;</mi><msup><mi>s</mi><mo>&prime;</mo></msup></munder><munder><mi>&Sigma;</mi><mi>i</mi></munder><msub><mi>&alpha;</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><msup><mi>s</mi><mo>&prime;</mo></msup><mo>)</mo></mrow><msub><mi>&gamma;</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><msubsup><mi>Y</mi><msup><mi>s</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow></mrow><mrow><munder><mi>&Sigma;</mi><mi>s</mi></munder><munder><mi>&Sigma;</mi><msup><mi>s</mi><mo>&prime;</mo></msup></munder><munder><mi>&Sigma;</mi><mi>i</mi></munder><msub><mi>&alpha;</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><msup><mi>s</mi><mo>&prime;</mo></msup><mo>)</mo></mrow><msub><mi>&gamma;</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><msubsup><mi>Y</mi><msup><mi>s</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>(i=0,1),后向迭代计算公式为<maths num="0010"><![CDATA[<math><mrow><msub><mi>&beta;</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munder><mi>&Sigma;</mi><msup><mi>s</mi><mo>&prime;</mo></msup></munder><munder><mi>&Sigma;</mi><mi>i</mi></munder><msub><mi>&beta;</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><msup><mi>s</mi><mo>&prime;</mo></msup><mo>)</mo></mrow><msub><mi>&gamma;</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><msubsup><mi>Y</mi><mi>t</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow></mrow><mrow><munder><mi>&Sigma;</mi><mi>s</mi></munder><munder><mi>&Sigma;</mi><msup><mi>s</mi><mo>&prime;</mo></msup></munder><munder><mi>&Sigma;</mi><mi>i</mi></munder><msub><mi>&alpha;</mi><mi>k</mi></msub><mrow><mo>(</mo><msup><mi>s</mi><mo>&prime;</mo></msup><mo>)</mo></mrow><msub><mi>&gamma;</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><msubsup><mi>Y</mi><msup><mi>s</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>(i=0,1)计算,其中,状态转移概率测度γ<sub>k</sub>(i,Y<sub>m</sub><sup>(i)</sup>)的全部4种可能值按下表计算。<img file="C2005100210160003C1.GIF" wi="1186" he="526" />
地址 610031四川省成都市二环路北一段111号