发明名称 一种基于博弈论的MIMO信道跟踪方法
摘要 本发明公开了一种基于博弈论的MIMO信道跟踪方法,其步骤为:发射导频信号,通过导频信号进行信道估计,并将信道估计值作为跟踪的初始值;建模MIMO信道噪声为博弈对手;分解跟踪误差,确定跟踪目标函数;通过博弈论求解目标函数鞍点实现信道跟踪。本方法结合了博弈论理论和MIMO信道跟踪方法,跟踪精度高。
申请公布号 CN101697495A 申请公布日期 2010.04.21
申请号 CN200910044697.1 申请日期 2009.11.05
申请人 湖南大学 发明人 曾凡仔;刘洁;李仁发;罗娟;肖玲
分类号 H04B7/02(2006.01)I 主分类号 H04B7/02(2006.01)I
代理机构 长沙市融智专利事务所 43114 代理人 颜勇
主权项 1.一种基于博弈论的MIMO信道跟踪算法,其特征在于,包括以下步骤:1)发射导频信号,通过导频信号进行信道估计,并将信道估计值作为跟踪的初始值:所述发射导频信号为发射符号向量为S:S=[s<sub>1</sub>,…,s<sub>m</sub>]<sup>T</sup>,m为发送端天线数;接收符号向量等式为Y=HS+η;其中,H为信道传输矩阵,<img file="F2009100446971C0000011.GIF" wi="745" he="264" />其中n为接收端天线数,h<sub>ij</sub>,i=1,…,n;j=1,…,m表示第j根发射天线到第i根接收天线的传输系数,h<sub>i</sub>=[h<sub>i1</sub>  …  h<sub>im</sub>]<sup>T</sup>表示第i根接收天线的传输向量,同时假定h<sub>i</sub>相互独立;接收符号向量为Y:Y=(y<sub>1</sub>,…,y<sub>n</sub>),其中y<sub>i</sub>表示第i个接收单元接收的符号;η=(η<sub>1</sub>,…,η<sub>n</sub>)为各分量相互独立的、均值为0、实部与虚部方差相等的复高斯随机序列;根据接收符号向量等式Y=HS+η,估计出信道参数h<sub>ij</sub>以作为信道跟踪的初始值<img file="F2009100446971C0000012.GIF" wi="439" he="254" />2)建模MIMO信道跟踪的博弈对手:信道估计和信道噪声为一对博弈对手,信道估计的效用函数为:<maths num="0001"><![CDATA[<math><mrow><msub><mi>J</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>K</mi><mi>k</mi></msub><mo>,</mo><msub><mi>L</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mi>N</mi></munderover><mi>E</mi><mrow><mo>(</mo><msub><mi>e</mi><mrow><mn>1</mn><mo>,</mo><mi>k</mi></mrow></msub><msubsup><mi>e</mi><mrow><mn>1</mn><mo>,</mo><mi>k</mi></mrow><mi>T</mi></msubsup><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>信道噪声的效用函数为:<maths num="0002"><![CDATA[<math><mrow><msub><mi>J</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>K</mi><mi>k</mi></msub><mo>,</mo><msub><mi>L</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mi>N</mi></munderover><mi>E</mi><mrow><mo>(</mo><msub><mi>e</mi><mrow><mn>2</mn><mo>,</mo><mi>k</mi></mrow></msub><msubsup><mi>e</mi><mrow><mn>2</mn><mo>,</mo><mi>k</mi></mrow><mi>T</mi></msubsup><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>其中:<maths num="0003"><![CDATA[<math><mrow><msub><mi>e</mi><mrow><mn>1</mn><mo>,</mo><mi>k</mi></mrow></msub><mo>=</mo><mrow><mo>(</mo><mi>A</mi><mo>-</mo><msub><mi>K</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><msubsup><mi>S</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>T</mi></msubsup><mo>+</mo><msub><mi>L</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><msub><mi>G</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><msub><mi>e</mi><mrow><mn>1</mn><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>K</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><msub><mi>v</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>e</mi><mn>1,0</mn></msub><mo>=</mo><msub><mi>H</mi><mn>0</mn></msub><mo>;</mo></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msub><mi>e</mi><mrow><mn>2</mn><mo>,</mo><mi>k</mi></mrow></msub><mo>=</mo><mrow><mo>(</mo><mi>A</mi><mo>-</mo><msub><mi>K</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><msubsup><mi>S</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>T</mi></msubsup><mo>+</mo><msub><mi>L</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><msub><mi>G</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><msub><mi>e</mi><mrow><mn>2</mn><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>L</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><msub><mi>n</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>e</mi><mn>2,0</mn></msub><mo>=</mo><mn>0</mn><mo>;</mo></mrow></math>]]></maths>E(e<sub>1,k</sub>e<sub>1,k</sub><sup>T</sup>)表示对误差e<sub>1,k</sub>的平方求期望、E(e<sub>2,k</sub>e<sub>2,k</sub><sup>T</sup>)表示对误差e<sub>2,k</sub>的平方求期望;其中,矩阵A=diag(α<sub>1×n</sub>),其中α<sub>1×n</sub>为n维行向量,向量的每个分量为α,α=J<sub>0</sub>(2πf<sub>d</sub>T<sub>s</sub>)exp(2πf<sub>0</sub>T<sub>s</sub>),其中J<sub>0</sub>()表示第一类零阶贝塞尔函数,f<sub>d</sub>为由于移动台与基站的相对位移引起的多普勒频移,f<sub>0</sub>为由于移动台与基站的晶震失配引起的载波频率偏移,T<sub>s</sub>为信道的时延扩展比,K<sub>k</sub>为信道估计增益<img file="F2009100446971C0000021.GIF" wi="314" he="66" />其中S<sub>k</sub>为k时刻的发送符号向量;L<sub>k</sub>是需要求解的信道噪声增益,G<sub>k</sub>为根据先验知识给定的矩阵,v<sub>k</sub>单位方差的白噪声序列且γ<sub>k</sub>互不相关;γ<sub>k</sub>为信道状态噪声序列;n<sub>k</sub>为白噪声序列,且与v<sub>k</sub>不相关;H<sub>0</sub>为用初始时刻的导频采用最小二乘算法估计得到的信道传输系数;3)确定跟踪目标函数:跟踪目标函数为<img file="F2009100446971C0000022.GIF" wi="749" he="124" />上式中,trace表示矩阵求迹运算;4)通过微分博弈理论,求解目标函数,实现MIMO信道跟踪目标函数的解为:<img file="F2009100446971C0000023.GIF" wi="800" he="81" />Y<sub>k</sub>表示k时刻Y接收符号向量,<img file="F2009100446971C0000024.GIF" wi="57" he="61" />表示H的k时刻的估计值,初始估计值<img file="F2009100446971C0000025.GIF" wi="176" he="61" />H<sub>k</sub>表示k时刻的信道传输矩阵,H<sub>k</sub>=AH<sub>k-1</sub>+γ<sub>k</sub>。
地址 410082 湖南省长沙市岳麓区麓山南路2号
您可能感兴趣的专利