发明名称 一种相干通信机用的比特-符号信号处理方法
摘要 本发明公开了一种适用于自适应判决反馈均衡器与译码器级连、迭代信号处理的比特-符号信号处理方法,该方法包括:1)确定与通信机接收的符号的后验概率最大对应的表达式;2)根据SOVA算法,使步骤1)的表达式值最大,从而使通信机接收的符号的后验概率最大,得到符号的软判决;3)根据SOVA算法,使步骤1)的表达式值最大,从而使通信机接收的符号的后验概率最大,得到符号的硬判决。
申请公布号 CN101119177B 申请公布日期 2010.05.12
申请号 CN200610089137.4 申请日期 2006.08.04
申请人 中国科学院声学研究所 发明人 朱维庆;朱敏
分类号 H04L1/00(2006.01)I;H04B13/02(2006.01)I 主分类号 H04L1/00(2006.01)I
代理机构 北京泛华伟业知识产权代理有限公司 11280 代理人 高存秀
主权项 1.一种相干通信机用的比特-符号信号处理方法,包括如下步骤:1)确定与通信机接收的符号的后验概率最大对应的表达式如下:<maths num="0001"><![CDATA[<math><mrow><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><msub><mrow><mi>log</mi><mi>p</mi></mrow><mi>t</mi></msub><mrow><mo>(</mo><msub><mi>s</mi><mi>t</mi></msub><mo>)</mo></mrow><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><mo>{</mo><msup><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>t</mi><mo>,</mo><mi>I</mi></mrow></msub><mo>-</mo><msubsup><mi>s</mi><mrow><mi>t</mi><mo>,</mo><mi>I</mi></mrow><mi>i</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>t</mi><mo>,</mo><mi>Q</mi></mrow></msub><mo>-</mo><msubsup><mi>s</mi><mrow><mi>t</mi><mo>,</mo><mi>Q</mi></mrow><mi>i</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>}</mo></mrow></math>]]></maths>其中p<sub>t</sub>(s<sub>t</sub>)是先验概率,r<sub>t,I</sub>和r<sub>t,Q</sub>分别为t时刻接收符号的虚部和实部,s<sub>t,I</sub><sup>i</sup>和s<sub>t,Q</sub><sup>i</sup>分别为星座图中对应信号的虚部和实部;2)根据SOVA算法,确定分量卷积码网格图中转移路经的分支度量为:<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>v</mi><mi>t</mi><msub><mi>s</mi><mi>t</mi></msub></msubsup><mo>=</mo><msup><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>t</mi><mo>,</mo><mi>I</mi></mrow></msub><mo>-</mo><msubsup><mi>s</mi><mrow><mi>t</mi><mo>,</mo><mi>I</mi></mrow><mi>i</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>t</mi><mo>,</mo><mi>Q</mi></mrow></msub><mo>-</mo><msubsup><mi>s</mi><mrow><mi>t</mi><mo>,</mo><mi>Q</mi></mrow><mi>i</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><msub><mrow><mi>log</mi><mi>p</mi></mrow><mi>t</mi></msub><mrow><mo>(</mo><msub><mi>s</mi><mi>t</mi></msub><mo>)</mo></mrow></mrow></math>]]></maths>3)对于转移路径X,确定路径度量为:<maths num="0003"><![CDATA[<math><mrow><msubsup><mi>&mu;</mi><mi>t</mi><mi>X</mi></msubsup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><msup><mi>t</mi><mo>&prime;</mo></msup><mo>=</mo><mn>1</mn></mrow><mi>t</mi></munderover><msubsup><mi>v</mi><msup><mi>t</mi><mo>&prime;</mo></msup><msub><mi>s</mi><msup><mi>t</mi><mo>&prime;</mo></msup></msub></msubsup><mo>=</mo><msubsup><mi>&mu;</mi><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mi>X</mi></msubsup><mo>+</mo><msubsup><mi>v</mi><mi>t</mi><msub><mi>s</mi><mi>t</mi></msub></msubsup></mrow></math>]]></maths>4)对于MPSK调制方式,在任一时刻生存路径都对应M-1条竞争路径,确定生存路径的度量值为:<maths num="0004"><![CDATA[<math><mrow><msub><mi>&mu;</mi><mrow><mi>T</mi><mo>,</mo><mi>min</mi></mrow></msub><mo>=</mo><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><msup><mi>t</mi><mo>&prime;</mo></msup><mo>=</mo><mn>1</mn></mrow><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow></munderover><msubsup><mi>v</mi><msup><mi>t</mi><mo>&prime;</mo></msup><msubsup><mi>s</mi><msup><mi>t</mi><mo>&prime;</mo></msup><mn>0</mn></msubsup></msubsup><mo>)</mo></mrow><mo>+</mo><msubsup><mi>v</mi><mi>t</mi><msubsup><mi>s</mi><mi>t</mi><mn>0</mn></msubsup></msubsup><mo>+</mo><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><msup><mi>t</mi><mrow><mo>&prime;</mo><mo>=</mo><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msup><mi>T</mi></munderover><msubsup><mi>v</mi><mi>t</mi><msubsup><mi>s</mi><msup><mi>t</mi><mo>&prime;</mo></msup><mn>0</mn></msubsup></msubsup><mo>)</mo></mrow><mo>=</mo><msup><mi>&Lambda;</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><msubsup><mi>s</mi><mi>t</mi><mn>0</mn></msubsup><mo>)</mo></mrow></mrow></math>]]></maths>确定与信号s<sup>1</sup>,s<sup>2</sup>,……s<sup>M-1</sup>对应的竞争路径的度量值为:<maths num="0005"><![CDATA[<math><mrow><msub><mi>&mu;</mi><msup><mrow><mi>t</mi><mo>,</mo><mi>s</mi></mrow><mi>i</mi></msup></msub><mo>=</mo><mi>min</mi><mo>{</mo><msubsup><mi>&mu;</mi><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow><mi>f</mi></msubsup><mrow><mo>(</mo><msubsup><mi>l</mi><msup><mi>s</mi><mi>i</mi></msup><mo>&prime;</mo></msubsup><mo>)</mo></mrow><mo>+</mo><msubsup><mi>v</mi><mi>t</mi><msup><mi>s</mi><mi>i</mi></msup></msubsup><mrow><mo>(</mo><msubsup><mi>l</mi><msup><mi>s</mi><mi>i</mi></msup><mo>&prime;</mo></msubsup><mo>,</mo><msub><mi>l</mi><msup><mi>s</mi><mi>i</mi></msup></msub><mo>)</mo></mrow><mo>+</mo><msubsup><mi>&mu;</mi><mi>t</mi><mi>b</mi></msubsup><mrow><mo>(</mo><msub><mi>l</mi><msup><mi>s</mi><mi>i</mi></msup></msub><mo>)</mo></mrow><mo>}</mo><mo>=</mo><msup><mi>&Lambda;</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><msubsup><mi>s</mi><mi>t</mi><mi>i</mi></msubsup><mo>)</mo></mrow></mrow></math>]]></maths>其中,t时刻与生存路径对应的是星座图中的信号s<sup>0</sup>,M-1条竞争路径此时分别对应信号s<sup>1</sup>、s<sup>2</sup>、……s<sup>M-1</sup>,<img file="F2006100891374C00016.GIF" wi="454" he="63" />为与信号s<sup>1</sup>、s<sup>2</sup>、……s<sup>M-1</sup>对应的竞争路径的度量值;且<img file="F2006100891374C00017.GIF" wi="444" he="66" />代表相邻时刻的状态,M<sub>S</sub>为网格图中的状态数,<img file="F2006100891374C00018.GIF" wi="153" he="65" />是t-1时刻前向生存路径度量值,<img file="F2006100891374C00019.GIF" wi="131" he="65" />为t时刻反向生存路径的度量值,<img file="F2006100891374C000110.GIF" wi="192" he="71" />为t时刻生成信号s<sup>i</sup>的状态转移分支度量;5)确定星座图中各信号的概率为:<maths num="0006"><![CDATA[<math><mrow><mi>p</mi><mrow><mo>(</mo><msubsup><mover><mi>s</mi><mo>^</mo></mover><mi>t</mi><mi>i</mi></msubsup><mo>=</mo><msubsup><mi>s</mi><mi>t</mi><mi>i</mi></msubsup><mo>)</mo></mrow><mo>=</mo><msup><mi>&Lambda;</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><msubsup><mi>s</mi><mi>t</mi><mi>i</mi></msubsup><mo>)</mo></mrow><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>M</mi><mo>-</mo><mn>1</mn></mrow></munderover><msup><mi>&Lambda;</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><msubsup><mi>s</mi><mi>t</mi><mi>i</mi></msubsup><mo>)</mo></mrow><mi>i</mi><mo>=</mo><mn>0,1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mi>M</mi><mo>-</mo><mn>1</mn></mrow></math>]]></maths>6)确定系统比特c<sub>t</sub>分别为0和1的先验概率为:<maths num="0007"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mi>p</mi><mrow><mo>(</mo><msub><mi>c</mi><mi>t</mi></msub><mo>=</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mfrac><msup><mi>e</mi><mrow><mi>&Lambda;</mi><mrow><mo>(</mo><msub><mi>c</mi><mi>t</mi></msub><mo>)</mo></mrow></mrow></msup><mrow><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><mi>&Lambda;</mi><mrow><mo>(</mo><msub><mi>C</mi><mi>t</mi></msub><mo>)</mo></mrow></mrow></msup></mrow></mfrac></mtd></mtr><mtr><mtd><mi>p</mi><mrow><mo>(</mo><msub><mi>c</mi><mi>t</mi></msub><mo>=</mo><mn>0</mn><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><mi>&Lambda;</mi><mrow><mo>(</mo><msub><mi>c</mi><mi>t</mi></msub><mo>)</mo></mrow></mrow></msup></mrow></mfrac></mtd></mtr></mtable></mfenced></math>]]></maths>其中,先验概率估计Λ(c<sub>t</sub>)为译码器输出的系统比特概率软信息;7)确定时刻t接收符号r<sub>t</sub>与星座图中各信号的欧氏距离分别为:<maths num="0008"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mi>L</mi><mrow><mo>(</mo><msubsup><mi>s</mi><mi>t</mi><mn>0</mn></msubsup><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>t</mi><mo>,</mo><mi>I</mi></mrow></msub><mo>-</mo><msubsup><mi>s</mi><mrow><mi>t</mi><mo>,</mo><mi>I</mi></mrow><mn>0</mn></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>t</mi><mo>,</mo><mi>Q</mi></mrow></msub><mo>-</mo><msubsup><mi>s</mi><mrow><mi>t</mi><mo>,</mo><mi>Q</mi></mrow><mn>0</mn></msubsup><mo>)</mo></mrow><mn>2</mn></msup></mtd></mtr><mtr><mtd><mi>L</mi><mrow><mo>(</mo><msubsup><mi>s</mi><mi>t</mi><mi>I</mi></msubsup><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>t</mi><mo>,</mo><mi>I</mi></mrow></msub><mo>-</mo><msubsup><mi>s</mi><mrow><mi>t</mi><mo>,</mo><mi>I</mi></mrow><mn>1</mn></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>t</mi><mo>,</mo><mi>Q</mi></mrow></msub><mo>-</mo><msubsup><mi>s</mi><mrow><mi>t</mi><mo>,</mo><mi>Q</mi></mrow><mn>1</mn></msubsup><mo>)</mo></mrow><mn>2</mn></msup></mtd></mtr><mtr><mtd><mi>L</mi><mrow><mo>(</mo><msubsup><mi>s</mi><mi>t</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>t</mi><mo>,</mo><mi>I</mi></mrow></msub><mo>-</mo><msubsup><mi>s</mi><mrow><mi>t</mi><mo>,</mo><mi>I</mi></mrow><mn>2</mn></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>t</mi><mo>,</mo><mi>Q</mi></mrow></msub><mo>-</mo><msubsup><mi>s</mi><mrow><mi>t</mi><mo>,</mo><mi>Q</mi></mrow><mn>2</mn></msubsup><mo>)</mo></mrow><mn>2</mn></msup></mtd></mtr><mtr><mtd><mfenced open='' close=''><mtable><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr></mtable></mfenced></mtd></mtr><mtr><mtd><mfenced open='' close=''><mtable><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr></mtable></mfenced></mtd></mtr><mtr><mtd><mfenced open='' close=''><mtable><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr></mtable></mfenced></mtd></mtr><mtr><mtd><mi>L</mi><mrow><mo>(</mo><msubsup><mi>s</mi><mi>t</mi><mrow><mi>M</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>t</mi><mo>,</mo><mi>I</mi></mrow></msub><mo>-</mo><msubsup><mi>s</mi><mrow><mi>t</mi><mo>,</mo><mi>I</mi></mrow><mrow><mi>M</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>t</mi><mo>,</mo><mi>Q</mi></mrow></msub><mo>-</mo><msubsup><mi>s</mi><mrow><mi>t</mi><mo>,</mo><mi>Q</mi></mrow><mrow><mi>M</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mn>2</mn></msup></mtd></mtr></mtable></mfenced></math>]]></maths>8)由步骤7),根据星座图,得t时刻校验位的概率:<maths num="0009"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mi>p</mi><mrow><mo>(</mo><mi>Check</mi><mo>=</mo><mn>0</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>p</mi><mrow><mo>(</mo><mi>Check</mi><mo>=</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mn>1</mn><mo>-</mo><mi>p</mi><mrow><mo>(</mo><mi>Check</mi><mo>=</mo><mn>0</mn><mo>)</mo></mrow></mtd></mtr></mtable></mfenced></math>]]></maths>9)由步骤6)、步骤7)和步骤8)的结果,根据星座图,得到p<sub>t</sub>(s<sub>t</sub>);10)由步骤4)、步骤5)和步骤9)的结果得到<img file="F2006100891374C00024.GIF" wi="222" he="74" />成为符号的软判决,用于判决反馈自适应均衡器和译码器的软迭代信号处理。
地址 100080 北京市海淀区北四环西路21号
您可能感兴趣的专利