发明名称 一种多天线正交频分复用系统中的基于子带的检测方法
摘要 一种多天线正交频分复用系统中基于子带的检测方法,属于无线数据传输技术领域,其特征在于,接收端将全部子载波分为若干个由连续子载波组成的子带。对每个子带,检测器计算模块根据估计出的信道矩阵,计算出某些子载波的检测矩阵,对这些检测矩阵进行拟合,将拟合参数传给检测模块。检测模块重构出所有子载波上的检测矩阵,对接收数据进行检测。当发送端采用预编码时,对每个子带,接收端的预编码计算模块得到预编码矩阵的拟合参数,将其同时传给检测器计算模块和发送端的预编码模块。预编码模块重构出预编码矩阵,对发送信号预编码。本方法降低了检测矩阵的计算量及其在不同模块间的数据传输量。有预编码时,可降低预编码矩阵的计算量和反馈量。
申请公布号 CN101459645A 申请公布日期 2009.06.17
申请号 CN200910076389.7 申请日期 2009.01.15
申请人 清华大学 发明人 周世东;陈翔;张秀军;肖元章;肖立民;李云洲;丁国鹏
分类号 H04L27/26(2006.01)I;H04B7/04(2006.01)I 主分类号 H04L27/26(2006.01)I
代理机构 北京众合诚成知识产权代理有限公司 代理人 朱 琨
主权项 1、一种多天线正交频分复用系统中的基于子带的检测方法,其特征在于,在接收端依次执行以下步骤:步骤(1)把接收到的全部有效子载波任意地分为若干个子带,每个子带包含多个连续的子载波;步骤(2)对于每一个所述的子带,检测器计算模块按以下步骤得到检测矩阵的拟合参数,并将所述的拟合参数传输到一个检测模块:步骤(2.1)计算当前所述子带内属于第m根发送天线的共n个导频符号处的信道估计值序列h<sub>m</sub>[k<sub>1</sub>],...,h<sub>m</sub>[k<sub>i</sub>],...,h<sub>m</sub>[k<sub>n</sub>],其中m为所述发送天线的序号,m=1,2,...,M,i为所述导频符号的序号,i=1,2,...,n,其中,n不大于当前所述子带内的子载波总数B,k<sub>i</sub>为所述导频符号所在子载波的序号,k<sub>i</sub>=k<sub>1</sub>,k<sub>2</sub>,...,k<sub>n</sub>,用y[k<sub>i</sub>]除以x[k<sub>i</sub>]的第m个元素,即可得到信道矩阵H[k<sub>i</sub>]的第m列上的信道估计值h<sub>m</sub>[k<sub>i</sub>],所述y[k<sub>i</sub>]为第k<sub>i</sub>个子载波上的接收信号,x[k<sub>i</sub>]为对应于所述第k<sub>i</sub>个子载波的发送信号,x[k<sub>i</sub>]的第m个元素为预先设定的导频符号,其余元素为0,步骤(2.2)对于所述信道估计值h<sub>m</sub>[k<sub>i</sub>]的第l个元素hl<sub>m</sub>[k<sub>i</sub>],用下述线性函数进行拟合,条件是均方根误差<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>E</mi><mi>lm</mi><mi>h</mi></msubsup><mo>=</mo><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msup><mrow><mo>(</mo><msub><mi>h</mi><mi>lm</mi></msub><mrow><mo>[</mo><msub><mi>k</mi><mi>i</mi></msub><mo>]</mo></mrow><mo>-</mo><msubsup><mi>a</mi><mi>lm</mi><mi>h</mi></msubsup><mo>&CenterDot;</mo><msub><mi>k</mi><mi>i</mi></msub><mo>-</mo><msubsup><mi>b</mi><mi>lm</mi><mi>h</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup></mrow></math>]]></maths>最小,从而得到当前所述子带内所有子载波上的信道响应<maths num="0002"><![CDATA[<math><mrow><msub><mi>h</mi><mi>lm</mi></msub><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow><mo>=</mo><msubsup><mi>a</mi><mi>lm</mi><mi>h</mi></msubsup><mo>&CenterDot;</mo><mi>k</mi><mo>+</mo><msubsup><mi>b</mi><mi>lm</mi><mi>h</mi></msubsup></mrow></math>]]></maths>其中,k为当前所述子带内的任一子载波的序号,参数<img file="A200910076389C00023.GIF" wi="49" he="56" />和<img file="A200910076389C00024.GIF" wi="46" he="55" />的表达式如下:<maths num="0003"><![CDATA[<math><mrow><msubsup><mi>a</mi><mi>lm</mi><mi>h</mi></msubsup><mo>=</mo><mfrac><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msub><mi>k</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mi>h</mi><mi>lm</mi></msub><mrow><mo>[</mo><msub><mi>k</mi><mi>i</mi></msub><mo>]</mo></mrow><mo>-</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><mo>&CenterDot;</mo><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msub><mi>k</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msub><mi>h</mi><mi>lm</mi></msub><mrow><mo>[</mo><msub><mi>k</mi><mi>i</mi></msub><mo>]</mo></mrow></mrow><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msubsup><mi>k</mi><mi>i</mi><mn>2</mn></msubsup><mo>-</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><mo>&CenterDot;</mo><msup><mrow><mo>(</mo><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msub><mi>k</mi><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msubsup><mi>b</mi><mi>lm</mi><mi>h</mi></msubsup><mo>=</mo><mfrac><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msub><mi>h</mi><mi>lm</mi></msub><mrow><mo>[</mo><msub><mi>k</mi><mi>i</mi></msub><mo>]</mo></mrow><mo>-</mo><msubsup><mi>a</mi><mi>lm</mi><mi>h</mi></msubsup><mo>&CenterDot;</mo><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msub><mi>k</mi><mi>i</mi></msub></mrow><mi>n</mi></mfrac><mo>,</mo></mrow></math>]]></maths>步骤(2.3)按步骤(2.1)和步骤(2.2)所述的方法对其它子带进行同样操作,步骤(2.4)按步骤(2.1)到步骤(2.3)所述的方法对其它发送天线进行同样的操作,步骤(2.5)计算当前所述子带内的第k<sub>1</sub>,k<sub>2</sub>,...,k<sub>n</sub>个子载波上的n个检测矩阵R[k<sub>1</sub>],...,R[k<sub>n</sub>],使得diag(R[k<sub>i</sub>]<sup>H</sup>H[k<sub>i</sub>])=1其中,diag(·)表示一个矩阵的对角线元素组成的向量,1为一个每个元素均为1的向量,R[k<sub>i</sub>]=R[k<sub>i</sub>]Λ[k<sub>i</sub>],其中,Λ[k<sub>i</sub>]是一个对角矩阵,且<maths num="0005"><![CDATA[<math><mrow><mi>R</mi><mrow><mo>[</mo><msub><mi>k</mi><mi>i</mi></msub><mo>]</mo></mrow><mo>=</mo><mi>H</mi><mrow><mo>[</mo><msub><mi>k</mi><mi>i</mi></msub><mo>]</mo></mrow><msup><mrow><mo>(</mo><mi>H</mi><msup><mrow><mo>[</mo><msub><mi>k</mi><mi>i</mi></msub><mo>]</mo></mrow><mi>H</mi></msup><mi>H</mi><mrow><mo>[</mo><msub><mi>k</mi><mi>i</mi></msub><mo>]</mo></mrow><mo>+</mo><msubsup><mi>&sigma;</mi><mi>n</mi><mn>2</mn></msubsup><mi>I</mi><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>,</mo></mrow></math>]]></maths>其中,R[k<sub>i</sub>]为第k<sub>i</sub>个子载波上的检测矩阵,H[k<sub>i</sub>]为等效信道矩阵,<img file="A200910076389C00032.GIF" wi="47" he="55" />为高斯白噪声的功率,I为单位矩阵,由此得到第k<sub>1</sub>,k<sub>2</sub>,...,k<sub>n</sub>个子载波上的检测矩阵R[k<sub>1</sub>],...,R[k<sub>n</sub>],步骤(2.6)根据所述的n个检测矩阵R[k<sub>1</sub>],...,R[k<sub>n</sub>],通过多项式拟合得到当前所述子带内所有子载波上的检测矩阵,步骤(2.6.1)对于所述检测矩阵R[k<sub>i</sub>]的第l行第m列的元素r<sub>lm</sub>[k<sub>i</sub>],用下述线性函数进行拟合<maths num="0006"><![CDATA[<math><mrow><msub><mover><mi>r</mi><mo>&OverBar;</mo></mover><mi>lm</mi></msub><mo>=</mo><mrow><mo>[</mo><msub><mi>k</mi><mi>i</mi></msub><mo>]</mo></mrow><mo>=</mo><msubsup><mi>a</mi><mi>lm</mi><mi>r</mi></msubsup><mo>&CenterDot;</mo><msub><mi>k</mi><mi>i</mi></msub><mo>+</mo><msubsup><mi>b</mi><mi>lm</mi><mi>r</mi></msubsup></mrow></math>]]></maths>其中,<img file="A200910076389C00034.GIF" wi="48" he="48" />和<img file="A200910076389C00035.GIF" wi="46" he="49" />为拟合参数,在均方误差<maths num="0007"><![CDATA[<math><mrow><msubsup><mi>E</mi><mi>lm</mi><mi>r</mi></msubsup><mo>=</mo><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msup><mrow><mo>(</mo><msub><mover><mi>r</mi><mo>&OverBar;</mo></mover><mi>lm</mi></msub><mrow><mo>[</mo><msub><mi>k</mi><mi>i</mi></msub><mo>]</mo></mrow><mo>-</mo><msubsup><mi>a</mi><mi>lm</mi><mi>r</mi></msubsup><mo>&CenterDot;</mo><msub><mi>k</mi><mi>i</mi></msub><mo>-</mo><msubsup><mi>b</mi><mi>lm</mi><mi>r</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup></mrow></math>]]></maths>最小的条件下,有<maths num="0008"><![CDATA[<math><mrow><msubsup><mi>a</mi><mi>lm</mi><mi>r</mi></msubsup><mo>=</mo><mfrac><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msub><mi>k</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mover><mi>r</mi><mo>&OverBar;</mo></mover><mi>lm</mi></msub><mrow><mo>[</mo><msub><mi>k</mi><mi>i</mi></msub><mo>]</mo></mrow><mo>-</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><mo>&CenterDot;</mo><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msub><mi>k</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msub><mover><mi>r</mi><mo>&OverBar;</mo></mover><mi>lm</mi></msub><mrow><mo>[</mo><msub><mi>k</mi><mi>i</mi></msub><mo>]</mo></mrow></mrow><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msubsup><mi>k</mi><mi>i</mi><mn>2</mn></msubsup><mo>-</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><mo>&CenterDot;</mo><msup><mrow><mo>(</mo><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msub><mi>k</mi><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac></mrow></math>]]></maths><maths num="0009"><![CDATA[<math><mrow><msubsup><mi>b</mi><mi>lm</mi><mi>h</mi></msubsup><mo>=</mo><mfrac><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msub><mover><mi>r</mi><mo>&OverBar;</mo></mover><mi>lm</mi></msub><mrow><mo>[</mo><msub><mi>k</mi><mi>i</mi></msub><mo>]</mo></mrow><mo>-</mo><msubsup><mi>a</mi><mi>lm</mi><mi>h</mi></msubsup><mo>&CenterDot;</mo><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msub><mi>k</mi><mi>i</mi></msub></mrow><mi>n</mi></mfrac><mo>,</mo></mrow></math>]]></maths>步骤(2.6.2)用步骤(2.6.1)所述方法得到当前所述子带内的检测矩阵的所有元素的拟合系数,步骤(2.7)所述检测器计算模块把步骤(2.6.2)所述的拟合系数传输给检测模块;步骤(3)所述检测模块根据所述步骤(2.7)传送来的检测矩阵拟合系数,按下式计算当前所述子带内第k个子载波上的检测矩阵的第l行第m列的元素:<maths num="0010"><![CDATA[<math><mrow><msub><mover><mi>r</mi><mo>&OverBar;</mo></mover><mi>lm</mi></msub><mo>=</mo><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow><mo>=</mo><msubsup><mi>a</mi><mi>lm</mi><mi>r</mi></msubsup><mo>&CenterDot;</mo><mi>k</mi><mo>+</mo><msubsup><mi>b</mi><mi>lm</mi><mi>r</mi></msubsup><mo>;</mo></mrow></math>]]></maths>步骤(4)所述检测模块根据步骤(3)得到的检测矩阵对接收信号进行检测。
地址 100084北京市100084-82信箱