发明名称 一种多带正交频分复用超宽带系统的信道估计方法
摘要 本发明公开了一种多带正交频分复用超宽带系统的信道估计方法,优点在于采用自相关特性较好的m序列作为时域训练序列并附加循环前缀,在接收端通过对去掉循环前缀的接收信号与训练序列作互相关运算和对各个训练序列作自相关运算来获得信道的冲激响应估计值,并利用m序列的自相关矩阵具有对角占优特性,首先分别通过对m序列的自相关矩阵进行一对角分解或三对角分解,然后采用一阶逆矩阵的逼近方法,有效的避免了复杂的求逆运算,从而使运算量降低了一个数量级,而性能逼近常规的时域信道估计方法,是一种超宽带系统的快速有效的信道估计方法,易于实现。
申请公布号 CN101447969A 申请公布日期 2009.06.03
申请号 CN200810164224.0 申请日期 2008.12.31
申请人 宁波大学 发明人 李有明;李新苗;徐铁锋
分类号 H04L27/26(2006.01)I;H04L25/03(2006.01)I 主分类号 H04L27/26(2006.01)I
代理机构 宁波奥圣专利代理事务所(普通合伙) 代理人 程晓明
主权项 1、一种多带正交频分复用超宽带系统的信道估计方法,包括以下步骤:①在发送端,首先对输入的数据信号进行正交相移调制处理得到调制信号;②然后对调制信号依次进行串并转换、傅里叶逆变换和并串转换处理,形成多个OFDM符号;③再在形成的多个OFDM符号中,每隔设定数量的OFDM符号插入一个长度为L<sub>P</sub>的m序列s,将m序列s作为一个训练序列,并根据信道特性在训练序列前附加一个长度为L<sub>C</sub>的循环前缀,得到附加循环前缀后的训练序列,用x表示,x=[x(0),x(1),…x(L<sub>P</sub>+L<sub>C</sub>-1)];④最后将附加循环前缀后的训练序列x和形成的OFDM符号一起经载波调制处理后通过超宽带信道传输至接收端,在传输过程中附加循环前缀后的训练序列x和OFDM符号受到信道衰落和高斯白噪声的影响;⑤在接收端,定义接收端接收到的经信道衰落和高斯白噪声影响后的附加循环前缀的训练序列x为第一接收信号,定义接收端接收到的经信道衰落和高斯白噪声影响后的OFDM符号为第二接收信号,将第一接收信号用抽头延迟线模型表示为<maths num="0001"><![CDATA[<math><mrow><mi>r</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>L</mi><mi>C</mi></msub><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>h</mi><mi>t</mi></msub><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>n</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>其中,k=0,1,…,L<sub>p</sub>+L<sub>c</sub>-1,r(k)为第k时刻的第一接收信号,h表示由信道的各个多径的系数构成的矩阵向量,<maths num="0002"><![CDATA[<math><mrow><mi>h</mi><mo>=</mo><msup><mrow><mo>[</mo><msub><mi>h</mi><mn>0</mn></msub><mo>,</mo><msub><mi>h</mi><mn>1</mn></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><msub><mi>h</mi><msub><mi>L</mi><mrow><mi>C</mi><mo>-</mo><mn>1</mn></mrow></msub></msub><mo>]</mo></mrow><mi>T</mi></msup><mo>,</mo></mrow></math>]]></maths>h<sub>t</sub>为信道的第t个多径系数,h应满足条件:{h<sub>t</sub>=0|L≤t≤L<sub>C</sub>-1},L为信道的阶数,x为附加循环前缀后的训练序列,x(k-t)为第k-t时刻的附加循环前缀后的训练序列,n为高斯白噪声,n(k)为第k时刻的高斯白噪声;⑥首先对第一接收信号r(k)进行去载波调制,并对去载波调制处理后的第一接收信号进行去循环前缀处理得到<maths num="0003"><![CDATA[<math><mrow><mover><mi>r</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>L</mi><mi>C</mi></msub><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>h</mi><mi>j</mi></msub><msup><mi>s</mi><mi>j</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mi>n</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>其中,k=0,1,…,L<sub>p</sub>+L<sub>c</sub>-1,<img file="A200810164224C0002111517QIETU.GIF" wi="111" he="71" />为去循环前缀后的第k时刻的第一接收信号,h表示由信道的各个多径系数构成的矩阵向量,<maths num="0004"><![CDATA[<math><mrow><mi>h</mi><mo>=</mo><msup><mrow><mo>[</mo><msub><mi>h</mi><mn>0</mn></msub><mo>,</mo><msub><mi>h</mi><mn>1</mn></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><msub><mi>h</mi><msub><mi>L</mi><mrow><mi>C</mi><mo>-</mo><mn>1</mn></mrow></msub></msub><mo>]</mo></mrow><mi>T</mi></msup><mo>,</mo></mrow></math>]]></maths>h<sub>j</sub>为信道的第j个多径系数,h应满足条件:{h<sub>j</sub>=0|L≤j≤L<sub>C</sub>-1},L为信道的阶数,n为高斯白噪声,n(k)为第k时刻的高斯白噪声,s<sup>j</sup>为m序列s循环右移j位后的m序列,s<sup>j</sup>(k)为m序列s循环右移j位后的第k时刻的序列;⑦然后计算去循环前缀后的第一接收信号<img file="A200810164224C00025.GIF" wi="19" he="32" />与m序列s循环右移i位后的m序列s<sup>i</sup>的互相关矩阵C和各个训练序列s的自相关矩阵C<sub>P</sub>,C=[C(i,j)],C(i,j)为去循环前缀后的第一接收信号<img file="A200810164224C00031.GIF" wi="20" he="32" />与m序列s循环右移i后的m序列s<sup>i</sup>的归一化互相关系数,<maths num="0005"><![CDATA[<math><mrow><mi>C</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>/</mo><msub><mi>L</mi><mi>P</mi></msub><mo>)</mo></mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>L</mi><mi>P</mi></msub><mo>-</mo><mn>1</mn></mrow></munderover><mover><mi>r</mi><mo>~</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msup><mi>s</mi><mi>i</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>C<sub>P</sub>=[C<sub>P</sub>(i,j)],C<sub>P</sub>(i,j)为m序列s循环右移j位后的m序列s<sup>j</sup>和m序列s循环右移i位后的m序列s<sup>i</sup>的归一化自相关系数,<maths num="0006"><![CDATA[<math><mrow><mrow><mfenced open='' close='' separators=' '><mtable><mtr><mtd><msub><mi>C</mi><mi>P</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>/</mo><msub><mi>L</mi><mi>P</mi></msub><mo>)</mo></mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>L</mi><mi>P</mi></msub><mo>-</mo><mn>1</mn></mrow></munderover><msup><mi>s</mi><mi>j</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msup><mi>s</mi><mi>i</mi></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>=</mo><mrow><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn><mo>,</mo></mtd><mtd><mi>i</mi><mo>=</mo><mi>j</mi></mtd></mtr><mtr><mtd><mo>-</mo><mn>1</mn><mo>/</mo><msub><mi>L</mi><mi>P</mi></msub><mo>,</mo></mtd><mtd><mi>i</mi><mo>&NotEqual;</mo><mi>j</mi></mtd></mtr></mtable></mfenced></mrow></mtd></mtr></mtable></mfenced></mrow><mo>,</mo></mrow></math>]]></maths>其中,i=0,1,…,L<sub>p</sub>,j=0,1,…,L<sub>p</sub>,k=0,1,…,L<sub>p</sub>+L<sub>c</sub>-1,<img file="A200810164224C00034.GIF" wi="80" he="43" />为去循环前缀后的第k时刻的第一接收信号,s<sup>j</sup>(k)为m序列s循环右移j位后的第k时刻的序列,s<sup>i</sup>(k)为m序列s循环右移i位后的第k时刻的序列。⑧再根据去循环前缀后的第一接收信号<img file="A200810164224C00035.GIF" wi="19" he="33" />与m序列s循环右移i位后的m序列s<sup>i</sup>的互相关矩阵C和各个训练序列s的自相关矩阵C<sub>P</sub>,计算信道的冲激响应估计值<img file="A200810164224C00036.GIF" wi="50" he="46" /><maths num="0007"><![CDATA[<math><mrow><mover><mi>h</mi><mo>~</mo></mover><mo>=</mo><msubsup><mi>C</mi><mi>p</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mi>C</mi><mo>,</mo></mrow></math>]]></maths>其中,<img file="A200810164224C00038.GIF" wi="48" he="52" />为自相关矩阵C<sub>P</sub>的逆矩阵;其特征在于根据所述的自相关矩阵C<sub>P</sub>的对角占优性,将所述的自相关矩阵C<sub>P</sub>分解为第一矩阵和第二矩阵之和,将所述的第一矩阵记为D,将所述的第二矩阵记为E,C<sub>P</sub>=D+E,在所述的第一矩阵D和所述的第二矩阵E满足‖D<sup>-1</sup>E‖&lt;1时,计算所述的自相关矩阵C<sub>P</sub>的逆矩阵<img file="A200810164224C00039.GIF" wi="92" he="53" /><maths num="0008"><![CDATA[<math><mrow><msubsup><mi>C</mi><mi>p</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>=</mo><mrow><mo>(</mo><mi>I</mi><mo>-</mo><msup><mi>D</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>E</mi><mo>+</mo><msup><mrow><mo>(</mo><msup><mi>D</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>E</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>+</mo><msup><mrow><mo>(</mo><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>m</mi></msup><msup><mrow><mo>(</mo><msup><mi>D</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>E</mi><mo>)</mo></mrow><mi>m</mi></msup><mo>+</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>)</mo></mrow><msup><mi>D</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>,</mo></mrow></math>]]></maths>其中,符号“‖‖”为范数符号,I为单位矩阵,D<sup>-1</sup>为第一矩阵D的逆矩阵,m=1,2,…,∞;再根据<maths num="0009"><![CDATA[<math><mrow><msubsup><mi>C</mi><mi>p</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>=</mo><mrow><mo>(</mo><mi>I</mi><mo>-</mo><msup><mi>D</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>E</mi><mo>+</mo><msup><mrow><mo>(</mo><msup><mi>D</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>E</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>+</mo><msup><mrow><mo>(</mo><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>m</mi></msup><msup><mrow><mo>(</mo><msup><mi>D</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>E</mi><mo>)</mo></mrow><mi>m</mi></msup><mo>+</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>)</mo></mrow><msup><mi>D</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup></mrow></math>]]></maths>计算<img file="A200810164224C000312.GIF" wi="47" he="51" />的一阶近似值,<maths num="0010"><![CDATA[<math><mrow><msubsup><mi>C</mi><mi>p</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>&ap;</mo><mrow><mo>(</mo><mi>I</mi><mo>-</mo><msup><mi>D</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>E</mi><mo>)</mo></mrow><msup><mi>D</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>&ap;</mo><msup><mi>D</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>-</mo><msup><mi>D</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>E</mi><msup><mi>D</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>.</mo></mrow></math>]]></maths>
地址 315211浙江省宁波市江北区风华路818号