发明名称 一种抗多径干扰LDPC解码器
摘要 本发明公开了一种抗多径干扰LDPC解码器,所提出的解码器包括解码模块和时变FIR滤波模块,时变FIR滤波模块用于对多径干扰进行抑制。本发明的解码器中的时变FIR滤波模块与LDPC解码模块独立,时变FIR滤波模块作为LDPC解码器的预处理器,在LDPC解码之前,消除大部分多径干扰;从而,LDPC解码器可获得很好的BER解码性能,能够对多径干扰环境下的信号正确译码。
申请公布号 CN101931417B 申请公布日期 2013.01.23
申请号 CN200910237861.0 申请日期 2009.11.12
申请人 北京交通大学 发明人 肖扬
分类号 H03M13/11(2006.01)I 主分类号 H03M13/11(2006.01)I
代理机构 北京市商泰律师事务所 11255 代理人 麻吉凤;毛燕生
主权项 1.一种抗多径干扰LDPC解码器,包括解码模块和滤波模块,其特征在于,所述解码模块与滤波模块级联,所述滤波模块用于对多径干扰进行抑制;该滤波模块为时变FIR滤波模块,该模块接收的信道信号为x(k)被划分为:<img file="FSB00000950112800011.GIF" wi="616" he="131" />其中,x<sub>d</sub>(k)为信道在数据时隙输出的比特序列;x<sub>t</sub>(k)为信道在训练时隙输出的比特序列;所述解码器还包括有系数调整模块,与所述时变FIR滤波模块相连接,用于在训练时隙调整所述时变FIR滤波模块的滤波系数以抑制多径干扰;所述系数调整模块包括:本地参考信号获取单元,用于依据LDPC编码的奇偶校验矩阵H,预先存储所述训练时隙比特序列的LDPC编码样本;最小均方算法调整单元,用于依据如下规则,对所述时变FIR滤波模块的系数进行调整:时变FIR滤波模块的系数向量w(k)=[w<sub>0</sub>(k),w<sub>1</sub>(k),w<sub>2</sub>(k),...,w<sub>N</sub>(k)]<sup>T</sup>、误差信号e(k)、训练时隙比特序列x<sub>t</sub>(k)以及收敛因子μ(k)满足的关系为,w(k+1)=w(k)+2μ(k)e(k)x<sub>t</sub>(k);其中,收敛因子为<img file="FSB00000950112800012.GIF" wi="419" he="96" />误差信号为<img file="FSB00000950112800013.GIF" wi="686" he="91" />d(k)为本地参考信号,本地参考信号d(k)是优先存储的训练时序比特序列的LDPC编码样本;其中,y(k)为时变FIR滤波模块的输出信号,x<sub>t</sub><sup>T</sup>(k)为x<sub>t</sub>(k)的转置;所述解码模块包括:初始化单元,用于依据<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>f</mi><mi>n</mi><mn>1</mn></msubsup><mo>=</mo><mn>1</mn><mo>/</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mn>2</mn><msub><mi>ay</mi><mi>n</mi></msub><mo>/</mo><msup><mi>&sigma;</mi><mn>2</mn></msup><mo>)</mo></mrow></mrow><mo>,</mo></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msubsup><mi>f</mi><mi>n</mi><mn>0</mn></msubsup><mo>=</mo><mn>1</mn><mo>-</mo><msubsup><mi>f</mi><mi>n</mi><mn>1</mn></msubsup></mrow></math>]]></maths>计算参数<img file="FSB00000950112800023.GIF" wi="205" he="57" />其中,σ<sup>2</sup>=N<sub>0</sub>/2为信道噪声的方差;a为信道衰落因子,y<sub>n</sub>为所述时变FIR滤波模块软输出;并且,令<img file="FSB00000950112800024.GIF" wi="178" he="52" /><maths num="0003"><![CDATA[<math><mrow><msubsup><mi>q</mi><mi>mm</mi><mn>0</mn></msubsup><mo>=</mo><msubsup><mi>f</mi><mi>n</mi><mn>0</mn></msubsup><mo>;</mo></mrow></math>]]></maths>水平迭代单元,用于对m=1,....,M,n∈N(m),依据<maths num="0004"><![CDATA[<math><mrow><msub><mi>dq</mi><mi>mn</mi></msub><mo>=</mo><msubsup><mi>q</mi><mi>mn</mi><mn>0</mn></msubsup><mo>-</mo><msubsup><mi>q</mi><mi>mn</mi><mn>1</mn></msubsup><mo>,</mo></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msub><mi>dr</mi><mi>mn</mi></msub><mo>=</mo><msubsup><mi>r</mi><mi>mn</mi><mn>0</mn></msubsup><mo>-</mo><msubsup><mi>r</mi><mi>mn</mi><mn>1</mn></msubsup><mo>=</mo><munder><mi>&Pi;</mi><mrow><msup><mi>n</mi><mo>&prime;</mo></msup><mo>&Element;</mo><mi>N</mi><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>/</mo><mi>n</mi></mrow></munder><msub><mi>dq</mi><msup><mi>mn</mi><mo>&prime;</mo></msup></msub><mo>,</mo></mrow></math>]]></maths>计算dq<sub>mn</sub>、dr<sub>mn</sub>;进而依据<maths num="0006"><![CDATA[<math><mrow><msubsup><mi>r</mi><mi>mn</mi><mn>0</mn></msubsup><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>dr</mi><mi>mn</mi></msub><mo>)</mo></mrow><mo>/</mo><mn>2</mn><mo>,</mo></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><msubsup><mi>r</mi><mi>mn</mi><mn>1</mn></msubsup><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>dr</mi><mi>mn</mi></msub><mo>)</mo></mrow><mo>/</mo><mn>2</mn></mrow></math>]]></maths>获取<img file="FSB000009501128000210.GIF" wi="206" he="59" />垂直迭代单元,用于依据获取的<img file="FSB000009501128000211.GIF" wi="51" he="58" />和<img file="FSB000009501128000212.GIF" wi="51" he="58" />更新概率值<img file="FSB000009501128000213.GIF" wi="61" he="59" />和<img file="FSB000009501128000214.GIF" wi="86" he="58" />其中,<maths num="0008"><![CDATA[<math><mrow><msubsup><mi>q</mi><mi>mn</mi><mn>0</mn></msubsup><mo>=</mo><msub><mi>&alpha;</mi><mi>mn</mi></msub><msubsup><mi>f</mi><mi>n</mi><mn>0</mn></msubsup><munder><mi>&Pi;</mi><mrow><msup><mi>m</mi><mo>&prime;</mo></msup><mo>&Element;</mo><mi>M</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>/</mo><mi>m</mi></mrow></munder><msubsup><mi>r</mi><mrow><msup><mi>m</mi><mo>&prime;</mo></msup><mi>n</mi></mrow><mn>0</mn></msubsup><mo>,</mo></mrow></math>]]></maths><maths num="0009"><![CDATA[<math><mrow><msubsup><mi>q</mi><mi>mn</mi><mn>1</mn></msubsup><mo>=</mo><msub><mi>&alpha;</mi><mi>mn</mi></msub><msubsup><mi>f</mi><mi>n</mi><mn>1</mn></msubsup><munder><mi>&Pi;</mi><mrow><msup><mi>m</mi><mo>&prime;</mo></msup><mo>&Element;</mo><mi>M</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>/</mo><mi>m</mi></mrow></munder><msubsup><mi>r</mi><mrow><msup><mi>m</mi><mo>&prime;</mo></msup><mi>n</mi></mrow><mn>1</mn></msubsup><mo>;</mo></mrow></math>]]></maths>α<sub>mn</sub>为归一化系数,使得<maths num="0010"><![CDATA[<math><mrow><msubsup><mi>q</mi><mi>mn</mi><mn>0</mn></msubsup><mo>+</mo><msubsup><mi>q</mi><mi>mn</mi><mn>1</mn></msubsup><mo>=</mo><mn>1</mn><mo>;</mo></mrow></math>]]></maths>并且,依据<maths num="0011"><![CDATA[<math><mrow><msubsup><mi>q</mi><mi>n</mi><mn>0</mn></msubsup><mo>=</mo><msub><mi>&alpha;</mi><mi>n</mi></msub><msubsup><mi>f</mi><mi>n</mi><mn>0</mn></msubsup><munder><mi>&Pi;</mi><mrow><mi>m</mi><mo>&Element;</mo><mi>M</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></munder><msubsup><mi>r</mi><mi>mn</mi><mn>0</mn></msubsup><mo>,</mo></mrow></math>]]></maths><maths num="0012"><![CDATA[<math><mrow><msubsup><mi>q</mi><mi>n</mi><mn>1</mn></msubsup><mo>=</mo><msub><mi>&alpha;</mi><mi>n</mi></msub><msubsup><mi>f</mi><mi>n</mi><mn>1</mn></msubsup><munder><mi>&Pi;</mi><mrow><mi>m</mi><mo>&Element;</mo><mi>M</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></munder><msubsup><mi>r</mi><mi>mn</mi><mn>0</mn></msubsup></mrow></math>]]></maths>计算伪后验概率<img file="FSB000009501128000220.GIF" wi="41" he="57" />和<img file="FSB000009501128000221.GIF" wi="67" he="57" />其中,α<sub>n</sub>为归一化系数,使得<img file="FSB000009501128000222.GIF" wi="220" he="56" />其中,N(m)为校验矩阵第m行的列号集合,该集合的总数目为N,M(n)为校验矩阵第n列的行号集合,该集合的总数目为M;译码单元,用于根据预置规则进行译码,所述预置规则为,当<img file="FSB000009501128000223.GIF" wi="110" he="58" />时令<img file="FSB000009501128000224.GIF" wi="135" he="58" />反之<img file="FSB000009501128000225.GIF" wi="140" he="58" />m=0;如果校验方程rH<sup>T</sup>=0,则译码成功;其中r为LDPC译码输出比特序列,用向量表示。
地址 100044 北京市海淀区西直门外上园村3号