发明名称 正交频分复用移动通信系统中的信道估计方法
摘要 正交频分复用移动通信系统中时频最优的信道估计方法包括以下步骤:步骤一,在发送端构造时域、频域二维导频数据结构;步骤二,在接收端获取导频符号处信道参数的最小二乘估计;步骤三,将上述导频符号处信道参数在二维离散余弦变换域进行单点滤波,并通过扩展的二维逆离散余弦变换获取所有数据符号处的信道参数。本发明克服了现有OFDM信道估计方法在系统存在保护带以及终端高速移动时性能恶化的问题,能以较低的复杂度实现逼近最优的估计性能。
申请公布号 CN101252555B 申请公布日期 2011.02.16
申请号 CN200810020077.X 申请日期 2008.03.28
申请人 东南大学 发明人 高西奇;江彬;王闻今;吴大焰;尤肖虎
分类号 H04L25/02(2006.01)I;H04L27/26(2006.01)I 主分类号 H04L25/02(2006.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 叶连生
主权项 1.一种正交频分复用移动通信系统中的信道估计方法,其特征在于该方法包括以下步骤:步骤一、在发送端构造时域、频域二维导频数据结构,在频率方向放置L<sub>F</sub>个间隔为D<sub>f</sub>的导频符号,在时间方向放置的L<sub>T</sub>个间隔为D<sub>t</sub>的导频符号,频率方向的导频符号由长度为L<sub>F</sub>的Zadoff-Chu序列构成,时间方向的导频符号是频率方向导频符号的重复;步骤二、在接收端,按如下公式获取导频符号处信道参数的最小二乘估计:<maths num="0001"><![CDATA[<math><mrow><msub><mover><mi>H</mi><mo>^</mo></mover><mi>p</mi></msub><mo>=</mo><mi>diag</mi><msup><mrow><mo>{</mo><msub><mi>X</mi><mi>p</mi></msub><mo>}</mo></mrow><mi>H</mi></msup><msub><mi>Y</mi><mi>p</mi></msub></mrow></math>]]></maths>其中,X<sub>p</sub>是由所有导频符号组成的向量,Y<sub>p</sub>是由所有接收导频符号组成的向量,<img file="FSB00000125107600012.GIF" wi="62" he="71" />是由所有导频符号处的频域信道参数的最小二乘估计值所组成的向量,diag{X<sub>p</sub>}表示主对角元为X<sub>p</sub>的对角阵,上标(·)<sup>H</sup>表示矩阵的共轭转置运算;步骤三、在接收端,对步骤二中得到的导频符号处的信道参数<img file="FSB00000125107600013.GIF" wi="73" he="69" />进行二维离散余弦变换,得到信号d,即<maths num="0002"><![CDATA[<math><mrow><mi>d</mi><mo>=</mo><mrow><mo>(</mo><msubsup><mi>C</mi><msub><mi>L</mi><mi>T</mi></msub><mi>II</mi></msubsup><mo>&CircleTimes;</mo><msubsup><mi>C</mi><msub><mi>L</mi><mi>F</mi></msub><mi>II</mi></msubsup><mo>)</mo></mrow><msub><mover><mi>H</mi><mo>^</mo></mover><mi>p</mi></msub></mrow></math>]]></maths>其中:<img file="FSB00000125107600015.GIF" wi="39" he="48" />表示矩阵的Kronecker乘积,<img file="FSB00000125107600016.GIF" wi="81" he="88" />和<img file="FSB00000125107600017.GIF" wi="68" he="67" />分别是长度为L<sub>T</sub>和L<sub>F</sub>的第II类离散余弦变换矩阵,其定义式为:<maths num="0003"><![CDATA[<math><mrow><msub><mrow><mo>[</mo><msubsup><mi>C</mi><msub><mi>L</mi><mi>T</mi></msub><mi>II</mi></msubsup><mo>]</mo></mrow><mrow><mi>k</mi><mo>,</mo><mi>l</mi></mrow></msub><mo>=</mo><msub><mi>&kappa;</mi><mi>k</mi></msub><mi>cos</mi><mfrac><mrow><mi>&pi;k</mi><mrow><mo>(</mo><mi>l</mi><mo>+</mo><mn>0.5</mn><mo>)</mo></mrow></mrow><msub><mi>L</mi><mi>T</mi></msub></mfrac></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msub><mrow><mo>[</mo><msubsup><mi>C</mi><msub><mi>L</mi><mi>F</mi></msub><mi>II</mi></msubsup><mo>]</mo></mrow><mrow><mi>k</mi><mo>,</mo><mi>l</mi></mrow></msub><mo>=</mo><msub><mi>&mu;</mi><mi>k</mi></msub><mi>cos</mi><mfrac><mrow><mi>&pi;k</mi><mrow><mo>(</mo><mi>l</mi><mo>+</mo><mn>0.5</mn><mo>)</mo></mrow></mrow><msub><mi>L</mi><mi>F</mi></msub></mfrac></mrow></math>]]></maths>其中:<maths num="0005"><![CDATA[<math><mrow><msub><mi>&kappa;</mi><mi>k</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn><mo>/</mo><msqrt><msub><mi>L</mi><mi>T</mi></msub></msqrt><mo>,</mo></mtd><mtd><mi>k</mi><mo>=</mo><mn>0</mn><mo>,</mo></mtd></mtr><mtr><mtd><msqrt><mn>2</mn></msqrt><mo>/</mo><msqrt><msub><mi>L</mi><mi>T</mi></msub></msqrt><mo>,</mo></mtd><mtd><mi>k</mi><mo>&NotEqual;</mo><mn>0</mn><mo>.</mo></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><msub><mi>&mu;</mi><mi>k</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn><mo>/</mo><msqrt><msub><mi>L</mi><mi>F</mi></msub></msqrt><mo>,</mo></mtd><mtd><mi>k</mi><mo>=</mo><mn>0</mn><mo>,</mo></mtd></mtr><mtr><mtd><msqrt><mn>2</mn></msqrt><mo>/</mo><msqrt><msub><mi>L</mi><mi>F</mi></msub></msqrt><mo>,</mo></mtd><mtd><mi>k</mi><mo>&NotEqual;</mo><mn>0</mn><mo>.</mo></mtd></mtr></mtable></mfenced><mo>;</mo></mrow></math>]]></maths>步骤四、在接收端,按如下公式对步骤三得到的信号d进行单点滤波:<maths num="0007"><![CDATA[<math><mrow><mover><mi>d</mi><mo>~</mo></mover><mo>=</mo><mi>&Gamma;d</mi></mrow></math>]]></maths>其中,<img file="FSB000001251076000113.GIF" wi="26" he="48" />是单点滤波后的信号,Γ是一个对角阵,其第i个主对角元素为E{|d<sub>i</sub>|<sup>2</sup>}-σ<sub>z</sub><sup>2</sup>,d<sub>i</sub>表示d的第i个元素,E{·}表示期望运算,σ<sub>z</sub><sup>2</sup>是噪声方差;步骤五、在接收端,对步骤四中得到的信号<img file="FSB00000125107600021.GIF" wi="26" he="47" />进行扩展的二维逆离散余弦变换,获取所有数据符号处的信道参数,即<maths num="0008"><![CDATA[<math><mrow><mover><mi>H</mi><mo>^</mo></mover><mo>=</mo><msup><mrow><mo>(</mo><msubsup><munder><mi>C</mi><mo>&OverBar;</mo></munder><msub><mi>L</mi><mi>T</mi></msub><mi>II</mi></msubsup><mo>&CircleTimes;</mo><msubsup><munder><mi>C</mi><mo>&OverBar;</mo></munder><msub><mi>L</mi><mi>F</mi></msub><mi>II</mi></msubsup><mo>)</mo></mrow><mi>T</mi></msup><mover><mi>d</mi><mo>~</mo></mover></mrow></math>]]></maths>其中<img file="FSB00000125107600023.GIF" wi="66" he="75" />和<img file="FSB00000125107600024.GIF" wi="69" he="67" />都是扩展的第II类离散余弦变换矩阵。<img file="FSB00000125107600025.GIF" wi="67" he="67" />的定义式为:<maths num="0009"><![CDATA[<math><mrow><msub><mrow><mo>[</mo><msubsup><munder><mi>C</mi><mo>&OverBar;</mo></munder><msub><mi>L</mi><mi>T</mi></msub><mi>II</mi></msubsup><mo>]</mo></mrow><mrow><mi>k</mi><mo>,</mo><mi>l</mi></mrow></msub><mo>=</mo><msub><mi>&kappa;</mi><mi>k</mi></msub><mi>cos</mi><mfrac><mrow><mi>&pi;k</mi><mrow><mo>(</mo><mi>l</mi><mo>/</mo><msub><mi>D</mi><mi>t</mi></msub><mo>+</mo><mn>0.5</mn><mo>)</mo></mrow></mrow><msub><mi>L</mi><mi>T</mi></msub></mfrac></mrow></math>]]></maths>其中:0≤k≤<sub>L</sub>T-1,0≤l≤N-1。<img file="FSB00000125107600027.GIF" wi="69" he="67" />的定义式为:<maths num="0010"><![CDATA[<math><mrow><msub><mrow><mo>[</mo><msubsup><munder><mi>C</mi><mo>&OverBar;</mo></munder><msub><mi>L</mi><mi>F</mi></msub><mi>II</mi></msubsup><mo>]</mo></mrow><mrow><mi>k</mi><mo>,</mo><mi>l</mi></mrow></msub><mo>=</mo><msub><mi>&mu;</mi><mi>k</mi></msub><mo>=</mo><mi>cos</mi><mfrac><mrow><mi>&pi;k</mi><mrow><mo>(</mo><mi>l</mi><mo>/</mo><msub><mi>D</mi><mi>f</mi></msub><mo>+</mo><mn>0.5</mn><mo>)</mo></mrow></mrow><msub><mi>L</mi><mi>F</mi></msub></mfrac></mrow></math>]]></maths>其中:0≤k≤L<sub>F</sub>-1,0≤l≤K-1。
地址 211109 江苏省南京市江宁开发区东南大学路2号