发明名称 一种MIMO系统中酉空时码的迭代检测方法
摘要 一种MIMO系统中酉空时码的迭代检测方法,属无线通信中的检测技术领域,其步骤如下:(1)非相干检测;(2)估计信道信息;(3)相干检测;(4)判断是否结束迭代。本发明方法解决了MIMO系统用于信道估计的导频开销问题和酉空时码的非相干检测性能较差的问题,提出了一种MIMO系统中酉空时码的迭代检测方法,目的在于提供MIMO系统中酉空时码的一种迭代检测方法,采用非相干检测和相干检测相结合的迭代方法,既不需要预先获知信道状态信息,从而避免了发送序列中的导频开销,又解决了酉空时码的单纯非相干检测性能较差的问题。
申请公布号 CN102098242B 申请公布日期 2013.08.14
申请号 CN201010597734.4 申请日期 2010.12.21
申请人 山东大学 发明人 杜正锋;袁东风;张彭;周晓天
分类号 H04L25/02(2006.01)I;H04B7/08(2006.01)I;H04L1/06(2006.01)I 主分类号 H04L25/02(2006.01)I
代理机构 济南金迪知识产权代理有限公司 37219 代理人 许德山
主权项 1.一种MIMO系统中酉空时码的迭代检测方法,应用于无线MIMO系统中,该检测方法步骤如下:(1)利用传统的酉空时码系统的非相干检测方法进行非相干检测:在第k个码块持续时间,接收到的信号为<maths num="0001"><![CDATA[<math><mrow><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msqrt><mfrac><mi>&rho;T</mi><mi>M</mi></mfrac></msqrt><mi>&Psi;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mi>h</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mi>w</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></math>]]></maths>式中,x(k)为接收信号向量,ρ为每个接收天线上的平均信噪比,T为每个码块持续时间包含的符号周期数,M为发送天线数,Ψ(k)为由酉空时码构成的发送信号矩阵,<img file="FSA00000392468300012.GIF" wi="400" he="61" />其中I<sub>N</sub>为N×N的单位矩阵,<img file="FSA00000392468300013.GIF" wi="38" he="38" />表示Kronecker乘积,Φ(k)为在第k个码块持续时间取自星座集<img file="FSA00000392468300014.GIF" wi="548" he="79" />中的T×M的酉空时码矩阵,上角标<img file="FSA00000392468300015.GIF" wi="84" he="28" />表示矩阵或向量的共轭转置,h(k)为复值信道增益向量,w(k)为接收到的噪声向量;在未知信道状态信息的条件下,假定在第k个码块持续时间内,<img file="FSA00000392468300016.GIF" wi="109" he="63" />是等概率发送的,进行基于最大似然准则的非相干检测<maths num="0002"><![CDATA[<math><mrow><mover><mi>&Psi;</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mi>arg</mi><munder><mi>max</mi><mrow><msub><mi>&Psi;</mi><mi>l</mi></msub><mo>&Element;</mo><mo>{</mo><msub><mi>&Psi;</mi><mn>1</mn></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>&Psi;</mi><mi>L</mi></msub><mo>}</mo></mrow></munder><msubsup><mrow><mo>|</mo><mo>|</mo><msubsup><mi>&Psi;</mi><mi>l</mi><mo>+</mo></msubsup><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo><mo>|</mo></mrow><mi>F</mi><mn>2</mn></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中<img file="FSA00000392468300018.GIF" wi="62" he="73" />表示求Frobenius范数,<img file="FSA00000392468300019.GIF" wi="283" he="59" />检测结果<img file="FSA000003924683000110.GIF" wi="114" he="69" />作为本发明迭代方法的初始值;(2)利用上一步非相干检测的结果,根据最小均方误差准则估计信道信息<img file="FSA000003924683000111.GIF" wi="125" he="70" />基于最小均方误差(MMSE)准则进行信道估计,定义第<img file="FSA000003924683000112.GIF" wi="213" he="73" />个码块持续时间的信道观测值如下<maths num="0003"><![CDATA[<math><mrow><mover><mi>h</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>=</mo><msqrt><mfrac><mi>M</mi><mi>&rho;T</mi></mfrac></msqrt><msup><mover><mi>&Psi;</mi><mo>^</mo></mover><mo>+</mo></msup><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>&tau;</mi><mo>)</mo></mrow><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><mo>=</mo><mover><mi>h</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>h</mi><mi>e</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>&tau;</mi><mo>)</mo></mrow></mrow></math>]]></maths>其中<maths num="0005"><![CDATA[<math><mrow><mover><mi>h</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>=</mo><mi>h</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>+</mo><msqrt><mfrac><mi>M</mi><mi>&rho;T</mi></mfrac></msqrt><msup><mover><mi>&Psi;</mi><mo>^</mo></mover><mo>+</mo></msup><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>&tau;</mi><mo>)</mo></mrow><mi>w</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>为仅包含噪声的信道观测值,而<maths num="0006"><![CDATA[<math><mrow><msub><mi>h</mi><mi>e</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>=</mo><mo>[</mo><msup><mover><mi>&Psi;</mi><mo>^</mo></mover><mo>+</mo></msup><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>-</mo><msup><mi>&Psi;</mi><mo>+</mo></msup><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>]</mo><mi>&Psi;</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>&tau;</mi><mo>)</mo></mrow><mi>h</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>则为估计误差分量;由仅包含噪声的观测值<img file="FSA00000392468300021.GIF" wi="172" he="66" />来求维纳滤波系数,将求得的维纳滤波系数乘以信道观测值<img file="FSA00000392468300022.GIF" wi="200" he="63" />得到信道的MMSE估计;定义<img file="FSA00000392468300023.GIF" wi="193" he="69" />为向量<img file="FSA00000392468300024.GIF" wi="175" he="66" />的第j个元素,j=1,...,MN,即<img file="FSA00000392468300025.GIF" wi="188" he="68" />为第(k+τ)个码块持续时间内第j个收发天线对之间的信道衰落系数;定义<img file="FSA00000392468300026.GIF" wi="1049" he="86" />其中上角标“T”表示矩阵或向量的转置,则h<sub>j</sub>(k)的基于<img file="FSA00000392468300027.GIF" wi="116" he="67" />的MMSE估计值为<maths num="0007"><![CDATA[<math><mrow><msub><mover><mi>h</mi><mo>^</mo></mover><mi>j</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msup><mi>c</mi><mo>+</mo></msup><mover><mi>h</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(5)中向量c为维纳滤波系数,可由如下Wiener-Hopf方程得到<maths num="0008"><![CDATA[<math><mrow><mi>c</mi><mo>=</mo><msubsup><mi>R</mi><mrow><mover><mi>h</mi><mo>&OverBar;</mo></mover><mover><mi>h</mi><mo>&OverBar;</mo></mover></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mi>r</mi><mrow><mover><mi>h</mi><mo>&OverBar;</mo></mover><mi>h</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(6)中<img file="FSA000003924683000210.GIF" wi="452" he="77" />为<img file="FSA000003924683000211.GIF" wi="117" he="69" />的自相关矩阵<img file="FSA000003924683000212.GIF" wi="445" he="77" />其中上角标“*”表示标量的共轭,即<img file="FSA000003924683000213.GIF" wi="46" he="43" />为第j个收发天线的信道衰落系数在第(k-τ<sub>1</sub>),(k-τ<sub>1</sub>+1),...,(k+τ<sub>2</sub>)个码块持续时间的观测值(仅包含噪声分量)和其在第k个码块持续时间的值之间的互相关向量;由于<img file="FSA000003924683000214.GIF" wi="737" he="75" />噪声采样点之间以及噪声与信道衰落系数之间相互独立,于是可以得到<img file="FSA000003924683000215.GIF" wi="67" he="54" />的第(s,t)个元素为<maths num="0009"><![CDATA[<math><mrow><msub><mi>r</mi><mi>st</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn><mo>+</mo><msqrt><mi>M</mi><mo>/</mo><mi>&rho;T</mi></msqrt><mo>,</mo></mtd><mtd><mi>s</mi><mo>=</mo><mi>t</mi></mtd></mtr><mtr><mtd><msub><mi>R</mi><mi>h</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>-</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo></mtd><mtd><mi>s</mi><mo>&NotEqual;</mo><mi>t</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>互相关向量<img file="FSA000003924683000217.GIF" wi="46" he="43" />可由下式得到<maths num="0010"><![CDATA[<math><mrow><msub><mi>r</mi><mrow><mover><mi>h</mi><mo>&OverBar;</mo></mover><mi>h</mi></mrow></msub><mo>=</mo><msup><mrow><mo>[</mo><msub><mi>R</mi><mi>h</mi></msub><mrow><mo>(</mo><mo>-</mo><msub><mi>&tau;</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>,</mo><msub><mi>R</mi><mi>h</mi></msub><mrow><mo>(</mo><mo>-</mo><msub><mi>&tau;</mi><mn>1</mn></msub><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>R</mi><mi>h</mi></msub><mrow><mo>(</mo><msub><mi>&tau;</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>]</mo></mrow><mi>T</mi></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>从式(6)-(8)可见,滤波系数c与通道编号j和码块编号k无关,能预先计算出来,将式(5)写为向量的形式,则有<maths num="0011"><![CDATA[<math><mrow><mover><mi>h</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mover><mi>H</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mi>c</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中<img file="FSA000003924683000220.GIF" wi="983" he="73" />为在第(k-τ<sub>1</sub>),(k-τ<sub>1</sub>+1),...,(k+τ<sub>2</sub>)个码块持续时间内由信道的观测值所构成的矩阵;(3)利用估计出的信道<img file="FSA000003924683000221.GIF" wi="98" he="69" />进行相干检测:估计出信道状态信息<img file="FSA000003924683000222.GIF" wi="99" he="70" />以后,利用此信道信息进行相干检测,即:<maths num="0012"><![CDATA[<math><mrow><mover><mi>&Psi;</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mi>arg</mi><munder><mi>max</mi><mrow><msub><mi>&Psi;</mi><mi>l</mi></msub><mo>&Element;</mo><mo>{</mo><msub><mi>&Psi;</mi><mn>1</mn></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>&Psi;</mi><mi>L</mi></msub><mo>}</mo></mrow></munder><msup><mrow><mo>|</mo><mo>|</mo><mi></mi><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msqrt><mfrac><mi>&rho;T</mi><mi>M</mi></mfrac></msqrt><msub><mi>&Psi;</mi><mi>l</mi></msub><mover><mi>h</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow></math>]]></maths>通过式(10)更新对于<img file="FSA00000392468300032.GIF" wi="114" he="67" />的估计;(4)判断是否结束迭代:若达到预先设定的迭代次数,结束迭代,输出检测结果<img file="FSA00000392468300033.GIF" wi="143" he="69" />否则,返回第二步,迭代过程继续进行。
地址 250100 山东省济南市历城区山大南路27号