发明名称 一种适用于并行传输系统符号检测的计算高效的MMSE-ⅡCPD算法
摘要 本发明实施例提供了一种适用于并行传输系统符号检测的计算高效的MMSE-IICPD算法,涉及通信技术领域,可以在检测性能维持不变的情况下,显著降低算法的复杂度。所述算法包括:针对并行传输的N个符号x=[x<sub>1</sub>,...,x<sub>N</sub>]<sup>T</sup>中的当前检测符号x<sub>i</sub>进行干扰抵消<img file="DSA0000112024150000011.GIF" wi="460" he="96" />对y<sub>i</sub>进行MMSE滤波,获得滤波结果:<img file="DSA0000112024150000012.GIF" wi="1550" he="88" /></maths>其中,f<sub>i,new</sub>是MMSE滤波矩阵F<sub>new</sub>的第i行,F<sub>new</sub>=DH<sup>H</sup>(HDH<sup>H</sup>+N<sub>0</sub>I<sub>M</sub>)<sup>-1</sup>,<img file="DSA0000112024150000013.GIF" wi="655" he="70" />,D=diag(d),d=[d<sub>1</sub>,...,d<sub>N</sub>]<sup>T</sup>为x=[x<sub>1</sub>,...,x<sub>N</sub>]<sup>T</sup>的剩余干扰方差,根据所述滤波结果计算当前检测符号x<sub>i</sub>在各调制星座点a<sub>k</sub>的度量值ω<sub>new</sub>(a<sub>k</sub>),根据所述x<sub>i</sub>在各星座点的度量值,计算得出x<sub>i</sub>的软判决和剩余干扰方差供下一次迭代使用。
申请公布号 CN104506470A 申请公布日期 2015.04.08
申请号 CN201410803041.4 申请日期 2014.12.12
申请人 西安电子科技大学 发明人 曹琲琲
分类号 H04L25/03(2006.01)I;H04L25/02(2006.01)I;H04B7/08(2006.01)I 主分类号 H04L25/03(2006.01)I
代理机构 代理人
主权项 一种适用于并行传输系统符号检测的计算高效的基于最小均方误差滤波的迭代干扰抵消并行检测MMSE‑IICPD算法,其特征在于,包括以下步骤:S1、针对并行传输的N个符号x=[x<sub>1</sub>,...,x<sub>N</sub>]<sup>T</sup>中的当前检测符号x<sub>i</sub>进行干扰抵消;其中,x<sub>i</sub>干扰抵消后的信号为:<img file="FSA0000112024170000011.GIF" wi="456" he="85" />公式中,r=[r<sub>1</sub>,...,r<sub>M</sub>]<sup>T</sup>为接收信号,h<sub>j</sub>为M行N列等效信道矩阵H的第j列向量,<img file="FSA0000112024170000012.GIF" wi="44" he="64" />为共道符号的N维软判决向量<img file="FSA0000112024170000013.GIF" wi="302" he="62" />的第i个元素,i=1,...,N,j=1,...,N;S2、对y<sub>i</sub>进行MMSE滤波,获得滤波结果:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>z</mi><mrow><mi>i</mi><mo>,</mo><mi>new</mi></mrow></msub><mo>=</mo><msub><mi>f</mi><mrow><mi>i</mi><mo>,</mo><mi>new</mi></mrow></msub><msub><mi>y</mi><mi>i</mi></msub><mo>=</mo><msub><mi>f</mi><mrow><mi>i</mi><mo>,</mo><mi>new</mi></mrow></msub><msub><mi>h</mi><mi>i</mi></msub><msub><mi>x</mi><mi>i</mi></msub><mo>+</mo><msubsup><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>&NotEqual;</mo><mn>1</mn></mrow><mi>N</mi></msubsup><msub><mi>f</mi><mrow><mi>i</mi><mo>,</mo><mi>new</mi></mrow></msub><msub><mi>h</mi><mi>j</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>j</mi></msub><mo>-</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mi>j</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>f</mi><mrow><mi>i</mi><mo>,</mo><mi>new</mi></mrow></msub><mi>w</mi><mo>=</mo><msub><mi>&mu;</mi><mrow><mi>i</mi><mo>,</mo><mi>new</mi></mrow></msub><msub><mi>x</mi><mi>i</mi></msub><mo>+</mo><msub><mi>n</mi><mrow><mi>i</mi><mo>,</mo><mi>new</mi></mrow></msub></mrow>]]></math><img file="FSA0000112024170000014.GIF" wi="1517" he="85" /></maths>其中,f<sub>i,new</sub>是MMSE滤波矩阵F<sub>new</sub>的第i行,F<sub>new</sub>=DH<sup>H</sup>(HDH<sup>H</sup>+N<sub>0</sub>I<sub>M</sub>)<sup>‑1</sup>,<img file="FSA0000112024170000015.GIF" wi="668" he="67" />D=diag(d),d=[d<sub>1</sub>,...,d<sub>N</sub>]<sup>T</sup>为x=[x<sub>1</sub>,...,x<sub>N</sub>]<sup>T</sup>的剩余干扰方差,diag(·)表示由向量构成的对角矩阵,接收信号r=Hx+w,其中w为M维加性复白高斯噪声列向量,w的方差为N<sub>0</sub>I<sub>M</sub>;当d<sub>i</sub><ζ时令d<sub>i</sub>=ζ;S3、根据所述滤波结果计算当前检测符号x<sub>i</sub>在各调制星座点a<sub>k</sub>的度量值ω<sub>new</sub>(a<sub>k</sub>),计算公式如下:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>&omega;</mi><mi>new</mi></msub><mrow><mo>(</mo><msub><mi>a</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mo>[</mo><mfrac><mrow><mo>-</mo><msup><mrow><mo>|</mo><msub><mi>z</mi><mrow><mi>i</mi><mo>,</mo><mi>new</mi></mrow></msub><mo>/</mo><msub><mi>&mu;</mi><mrow><mi>i</mi><mo>,</mo><mi>new</mi></mrow></msub><mo>-</mo><msub><mi>a</mi><mi>k</mi></msub><mo>|</mo></mrow><mn>2</mn></msup></mrow><mrow><msubsup><mi>&sigma;</mi><msub><mi>n</mi><mrow><mi>i</mi><mo>,</mo><mi>new</mi></mrow></msub><mn>2</mn></msubsup><mo>/</mo><msubsup><mi>&mu;</mi><mrow><mi>i</mi><mo>,</mo><mi>new</mi></mrow><mn>2</mn></msubsup></mrow></mfrac><mo>]</mo><mo>,</mo><msub><mi>a</mi><mi>k</mi></msub><mo>&Element;</mo><mi>&Omega;</mi></mrow>]]></math><img file="FSA0000112024170000016.GIF" wi="1303" he="231" /></maths>其中,滤波输出干扰加噪声能量<img file="FSA0000112024170000017.GIF" wi="578" he="67" />Ω是调制星座点集合,k=1,...,K,K是星座点数量或称星座规模;S4、根据所述x<sub>i</sub>在各星座点的度量值,计算得出x<sub>i</sub>的软判决和剩余干扰方差供下一次迭代使用。
地址 710071 陕西省西安市太白南路2号