发明名称 一种修正的模判决多模盲均衡方法
摘要 本发明公开了一种修正的模判决多模盲均衡方法MMDMMA,重新定义了代价函数,由所定义的代价函数按随机梯度下降法获得前向均衡器和反馈滤波器的权向量更新公式;使用具有时变特性的泄漏因子替现有模判决多模盲均衡方法MDMMA的固定判决因子,能够准确跟踪时变信号。这种替代的效果是:在收敛前,判决误差大,时变泄漏因子也大,收敛速度快;在收敛后,判决误差小,时变泄漏因子也很小,收敛速度慢;采用时变泄漏因子避免了MDMMA方法盲均衡过程的发散又避免了误调。本发明方法均方误差小、收敛速度快、抑制码间干扰和纠正相位旋转能力强,相较模判决多模盲均衡方法MDMMA和常模盲均衡方法CMA有明显的提高。
申请公布号 CN103516648A 申请公布日期 2014.01.15
申请号 CN201310493509.X 申请日期 2013.10.18
申请人 南京信息工程大学 发明人 郭业才;王南南;张鑫凤;黄友锐
分类号 H04L25/03(2006.01)I 主分类号 H04L25/03(2006.01)I
代理机构 南京众联专利代理有限公司 32206 代理人 顾进;叶涓涓
主权项 1.一种修正的模判决多模盲均衡方法,包括如下步骤:步骤A,发射信号a(k)经过脉冲响应信道h(k)并加入信道高斯白噪声n(k)后,得到前向均衡器输入信号y(k):y(k)=a(k)h(k)+n(k),其中k为正整数且表示时间序列;步骤B,步骤A所述的前向均衡器输入信号y(k)经前向均衡器得到前向均衡器输出信号z<sub>F</sub>(k):<img file="FDA0000398417560000011.GIF" wi="424" he="77" />其中,f<sub>F</sub>(k)为前向均衡器的权向量;上标T表示转置;下标F表示前向均衡器;步骤C,基于步骤B所述的前向均衡器输出信号z<sub>F</sub>(k)和反馈滤波器输出信号z<sub>B</sub>(k)得到判决器Qu的输入信号z(k):z(k)=z<sub>F</sub>(k)-z<sub>B</sub>(k),其中<img file="FDA0000398417560000012.GIF" wi="417" he="75" />f<sub>B</sub>(k)为反馈滤波器的权向量;下标B表示反馈滤波器;步骤D,步骤C所述的判决器输入信号z(k)经判决器Qu判决得到发射信号序列a(k)的估计<img file="FDA0000398417560000013.GIF" wi="136" he="75" />其特征在于,所述前向均衡器的权向量f<sub>F</sub>(k)通过以下公式进行更新:f<sub>F</sub>(k+1)=f<sub>F</sub>(k)-μ<sub>F</sub>[ρ(k)(e<sub>MM_R</sub>(k)+je<sub>MM_I</sub>(k))+(1-ρ(k))(e<sub>R</sub>(k)+je<sub>I</sub>(k))]y<sup>*</sup>(k);其中,μ<sub>F</sub>为前馈滤波器权向量f<sub>F</sub>(k)的迭代步长,e<sub>MM_R</sub>(k)与e<sub>MM_I</sub>(k)为多模误差信号e<sub>MM</sub>(k)的实部和虚部,而e<sub>MM</sub>(k)=z(k)[|z(k)|<sup>2</sup>-r<sub>MM</sub>(k)],r<sub>MM</sub>(k)是k时刻判决器输入信号z(k)的星座图中离z(k)最近的圆上所有的发射信号统计的平方,下标MM表示多模;e<sub>R</sub>(k)与e<sub>I</sub>(k)为z(k)通过误差生成函数生成的误差信号e(k)的实部和虚部,且e(k)=z(k)[|z(k)|<sup>2</sup>-r<sub>CM</sub>],r<sub>CM</sub>=E[a<sup>4</sup>(k)]/E[a<sup>2</sup>(k)]为发射信号a(k)的常模值,下标CM表示常模;<img file="FDA0000398417560000019.GIF" wi="188" he="83" />为虚数单位;所述反馈滤波器的权向量f<sub>B</sub>(k)通过以下公式进行更新:<maths num="0001"><![CDATA[<math><mrow><mrow><msub><mi>f</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msub><mi>f</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>&mu;</mi><mi>B</mi></msub><mo>[</mo><mi>&rho;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>e</mi><mrow><mi>MM</mi><mo>_</mo><mi>R</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>je</mi><mrow><mi>MM</mi><mo>_</mo><mi>I</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&rho;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>e</mi><mi>R</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>je</mi><mi>I</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mo>]</mo></mrow><msup><mover><mi>a</mi><mo>^</mo></mover><mo>*</mo></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>式中,μ<sub>B</sub>为反馈滤波器权向量f<sub>B</sub>(k)的迭代步长;ρ(k)为新的时变泄漏因子,ρ(k)的数学表达式为:<img file="FDA00003984175600000110.GIF" wi="1616" he="240" />其中,|x|表示x的绝对值;r<sub>CM_R</sub>与r<sub>CM_I</sub>是发射信号a(k)常数模r<sub>CM</sub>=E[a<sup>4</sup>(k)]/E[a<sup>2</sup>(k)]的实部与虚部,<img file="FDA0000398417560000017.GIF" wi="587" he="84" /><img file="FDA0000398417560000018.GIF" wi="578" he="84" />a<sub>R</sub>(k)与a<sub>I</sub>(k)是发射信号a(k)的实部与虚部;且r<sub>CM_R</sub>=r<sub>CM_I</sub>;r<sub>MM_R</sub>(k)、r<sub>MM_I</sub>(k)分别是r<sub>MM</sub>(k)的实部和虚部,且满足<img file="FDA0000398417560000021.GIF" wi="796" he="110" />min表示取最小值;r<sub>MM_i</sub>表示第i个圆半径的平方;<img file="FDA0000398417560000027.GIF" wi="97" he="63" />为判决误差函数,且<maths num="0002"><![CDATA[<math><mrow><mover><mi>e</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>z</mi><mi>R</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mrow><mo>(</mo><msup><mrow><mo>|</mo><msub><mi>z</mi><mi>R</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><msup><mrow><mo>|</mo><msub><mover><mi>a</mi><mo>^</mo></mover><mi>R</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><mo>+</mo><msub><mi>jz</mi><mi>I</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mrow><mo>(</mo><msup><mrow><mo>|</mo><msub><mi>z</mi><mi>I</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><msup><mrow><mo>|</mo><msub><mover><mi>a</mi><mo>^</mo></mover><mi>I</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>其中<img file="FDA0000398417560000024.GIF" wi="133" he="73" />与<img file="FDA0000398417560000025.GIF" wi="120" he="72" />分别是判决输出<img file="FDA0000398417560000026.GIF" wi="106" he="70" />的实部和虚部;z<sub>R</sub>(k)与z<sub>I</sub>(k)分别是z(k)的实部和虚部;α&gt;0、β&gt;0。
地址 210044 江苏省南京市宁六路219号