主权项 |
一种适用于并行传输系统符号检测的计算高效的基于最小均方误差滤波的迭代干扰抵消并行检测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>Σ</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>≠</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>μ</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>ω</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>μ</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>σ</mi><msub><mi>n</mi><mrow><mi>i</mi><mo>,</mo><mi>new</mi></mrow></msub><mn>2</mn></msubsup><mo>/</mo><msubsup><mi>μ</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>∈</mo><mi>Ω</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>的软判决和剩余干扰方差供下一次迭代使用。 |