发明名称 高速移动环境下莱斯信道中多普勒频偏估计与补偿方法
摘要 本发明公开了一种高速移动环境下莱斯信道中多普勒频偏估计与补偿方法,主要解决现有技术中多普勒频偏估计精度不高和非最佳倍数补偿问题。其基本步骤为:首先利用接收信号通过频域信道估计算法或时域信道估计算法得到时域信道估计;再从时域信道估计中提取含有LOS分量的径;然后利用该含有LOS分量的径进行频偏估计得到高精度的频偏估计值,并通过最大化接收信号的信干比来确定最佳的频偏补偿倍数;最后对接收信号进行最佳频偏补偿,消除接收信号中多普勒频偏的影响。本发明有效地减小了多普勒频偏带来的系统性能损失,可用于提高OFDM/OFDMA系统对多普勒频偏的鲁棒性。
申请公布号 CN102404268B 申请公布日期 2014.02.12
申请号 CN201110366374.1 申请日期 2011.11.17
申请人 西安电子科技大学 发明人 任光亮;杨丽花;翟万涛;张会宁
分类号 H04L27/26(2006.01)I;H04L25/02(2006.01)I 主分类号 H04L27/26(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 1.一种高速移动环境下莱斯信道的多普勒频偏估计与补偿方法,包括:(1)频偏估计步骤:(1a)提取带有直射分量LOS的径首先,提取频域接收信号中的正交频分复用/正交频分多址OFDM/OFDMA导频符号,并利用该接收信号中的导频符号和发送信号中的导频符号进行信道估计,得到导频位置频域信道响应:<img file="FDA0000406857160000011.GIF" wi="567" he="173" />其中m<sub>p</sub>为带有导频符号的标识,k<sub>p</sub>为导频子载波位置的标识,Y(m<sub>p</sub>,k<sub>p</sub>)为位置(m<sub>p</sub>,k<sub>p</sub>)处的频域接收信号,X(m<sub>p</sub>,k<sub>p</sub>)为位置(m<sub>p</sub>,k<sub>p</sub>)处的频域发送信号;其次,对导频位置频域信道响应<img file="FDA0000406857160000012.GIF" wi="198" he="94" />进行内插,得到整个导频符号上其他数据子载波上的信道频域响应<img file="FDA0000406857160000013.GIF" wi="226" he="92" />并将当前带有导频符号的上的频域信道响应估计值记为<img file="FDA0000406857160000014.GIF" wi="208" he="94" />其中k<sub>d</sub>为数据子载波位置的标识,k为系统占用的子载波位置的标识;再次,对频域信道响应<img file="FDA0000406857160000015.GIF" wi="180" he="95" />进行IFFT变换处理,得到时域信道估计值<img file="FDA0000406857160000016.GIF" wi="638" he="150" />式中,<img file="FDA0000406857160000017.GIF" wi="184" he="100" />是时域信道第一径系数的估计值,<img file="FDA00004068571600000113.GIF" wi="220" he="86" />是时域信道第L-1径系数的估计值,<img file="FDA0000406857160000018.GIF" wi="78" he="89" />为整个OFDM/OFDMA符号上信道响应的平均值,且位于该OFDM/OFDMA符号的中间位置;最后,从时域信道估计值<img file="FDA0000406857160000019.GIF" wi="76" he="88" />中提取带有LOS分量的径,即第一径:<maths num="0001"><![CDATA[<math><mrow><msub><mover><mi>h</mi><mo>^</mo></mover><mi>los</mi></msub><mrow><mo>(</mo><msub><mi>m</mi><mi>p</mi></msub><mo>)</mo></mrow><mo>=</mo><mover><mi>h</mi><mo>^</mo></mover><mrow><mo>(</mo><msub><mi>m</mi><mi>p</mi></msub><mo>,</mo><mn>0</mn><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>(1b)利用估计得到的带有LOS分量的径<img file="FDA00004068571600000111.GIF" wi="185" he="105" />进行频偏估计,得到多普勒频偏估计值:<maths num="0002"><![CDATA[<math><mrow><mover><mi>&epsiv;</mi><mo>^</mo></mover><mo>=</mo><mfrac><mi>N</mi><mrow><mn>2</mn><mi>&pi;</mi><msub><mi>N</mi><mi>s</mi></msub><msub><mi>DN</mi><mi>m</mi></msub></mrow></mfrac><munder><mi>&Sigma;</mi><mrow><msub><mi>m</mi><mi>p</mi></msub><mo>&Element;</mo><msub><mi>&kappa;</mi><mi>p</mi></msub></mrow></munder><mi>arg</mi><mo>{</mo><msub><mover><mi>h</mi><mo>^</mo></mover><mi>los</mi></msub><mrow><mo>(</mo><msub><mi>m</mi><mi>p</mi></msub><mo>+</mo><mi>D</mi><mo>)</mo></mrow><msup><msub><mover><mi>h</mi><mo>^</mo></mover><mi>los</mi></msub><mo>*</mo></msup><mrow><mo>(</mo><msub><mi>m</mi><mi>p</mi></msub><mo>)</mo></mrow><mo>}</mo></mrow></math>]]></maths>式中,κ<sub>p</sub>为带有导频的OFDM/OFDMA符号位置集合,m<sub>p</sub>为带有导频的OFDM/OFDMA符号的标识,N<sub>s</sub>=N+N<sub>g</sub>,N和N<sub>g</sub>分别为FFT和循环前缀的长度,arg{}表示取角度符号,D为两个进行相关处理的OFDM/OFDMA符号间隔的符号数,N<sub>m</sub>为进行相关处理的OFDM/OFDMA组数,<img file="FDA0000406857160000021.GIF" wi="269" he="105" />为位置m<sub>p</sub>+D处带有LOS分量的径系数估计值,<img file="FDA0000406857160000022.GIF" wi="202" he="122" />表示<img file="FDA0000406857160000023.GIF" wi="185" he="105" />的共轭值;(2)最佳多普勒频偏补偿步骤:(2a)计算频偏补偿后信号的信干比SIR假设用于频偏补偿的频偏值为<img file="FDA0000406857160000024.GIF" wi="92" he="93" />其中ρ∈[0 1]为补偿倍数,<img file="FDA0000406857160000025.GIF" wi="43" he="79" />为频偏的估计值,则频偏补偿后信号的信干比SIR为:<maths num="0003"><![CDATA[<math><mrow><mi>SIR</mi><mo>=</mo><mo>[</mo><msub><mi>P</mi><mi>L</mi></msub><msup><mrow><mo>|</mo><mfrac><mrow><mi>sin</mi><mi>&pi;</mi><munder><mi>&epsiv;</mi><mo>&OverBar;</mo></munder></mrow><mrow><mi>sin</mi><mi>&pi;</mi><munder><mi>&epsiv;</mi><mo>&OverBar;</mo></munder><mo>/</mo><mi>N</mi></mrow></mfrac><mo>|</mo></mrow><mn>2</mn></msup><mo>+</mo><msub><mi>P</mi><mi>S</mi></msub><mrow><mo>(</mo><mi>N</mi><mo>+</mo><mn>2</mn><munderover><mi>&Sigma;</mi><mrow><mi>r</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>J</mi><mn>0</mn></msub><mrow><mo>(</mo><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>max</mi></msub><mi>r</mi><msub><mi>T</mi><mi>s</mi></msub><mo>)</mo></mrow><mi>cos</mi><mrow><mo>(</mo><mn>2</mn><mi>&pi;&rho;</mi><mover><mi>&epsiv;</mi><mo>^</mo></mover><mi>r</mi><mo>/</mo><mi>N</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>]</mo><mo>/</mo><msub><mi>P</mi><mi>ICI</mi></msub></mrow></math>]]></maths>式中,f<sub>max</sub>为最大多普勒频移,T<sub>s</sub>为采样间隔,J<sub>0</sub>(·)为第一类零阶贝塞尔函数,<img file="FDA0000406857160000027.GIF" wi="195" he="86" />ε为真实归一化频偏,P<sub>L</sub>为莱斯信道LOS分量的功率,P<sub>L</sub>=|c|<sup>2</sup>,c为信道第一径LOS分量系数,P<sub>S</sub>为莱斯信道所有散射分量的功率和,<img file="FDA0000406857160000028.GIF" wi="311" he="142" /><img file="FDA00004068571600000210.GIF" wi="51" he="66" />为信道第l<sub>p</sub>径的系数,var[]为取方差符号,P<sub>ICI</sub>为平均干扰功率:<maths num="0004"><![CDATA[<math><mfenced open='' close=''><mtable><mtr><mtd><msub><mi>P</mi><mi>ICI</mi></msub><mo>=</mo><mfrac><msub><mi>P</mi><mi>L</mi></msub><mrow><msub><mi>N</mi><mi>u</mi></msub><msup><mi>N</mi><mn>2</mn></msup></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>p</mi><mo>=</mo><mo>-</mo><msub><mi>N</mi><mi>u</mi></msub><mo>+</mo><mn>1</mn></mrow><mrow><msub><mi>N</mi><mi>u</mi></msub><mo>-</mo><mn>1</mn></mrow></munderover><mrow><mo>(</mo><msub><mi>N</mi><mi>u</mi></msub><mo>-</mo><mo>|</mo><mi>p</mi><mo>|</mo><mo>)</mo></mrow><msup><mrow><mo>|</mo><mfrac><mrow><mi>sin</mi><mi>&pi;</mi><mrow><mo>(</mo><mi>p</mi><mo>+</mo><munder><mi>&epsiv;</mi><mo>&OverBar;</mo></munder><mo>)</mo></mrow></mrow><mrow><mi>sin</mi><mo>[</mo><mi>&pi;</mi><mrow><mo>(</mo><mi>p</mi><mo>+</mo><munder><mi>&epsiv;</mi><mo>&OverBar;</mo></munder><mo>)</mo></mrow><mo>/</mo><mi>N</mi><mo>]</mo></mrow></mfrac><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><mfrac><msub><mi>P</mi><mi>L</mi></msub><msup><mi>N</mi><mn>2</mn></msup></mfrac><msup><mrow><mo>|</mo><mfrac><mrow><mi>sin</mi><mi>&pi;</mi><munder><mi>&epsiv;</mi><mo>&OverBar;</mo></munder></mrow><mrow><mi>sin</mi><mo>[</mo><mi>&pi;</mi><munder><mi>&epsiv;</mi><mo>&OverBar;</mo></munder><mo>/</mo><mi>N</mi><mo>]</mo></mrow></mfrac><mo>|</mo></mrow><mn>2</mn></msup><mo>+</mo><mfrac><mrow><msub><mi>N</mi><mi>u</mi></msub><mo>-</mo><mn>1</mn></mrow><mi>N</mi></mfrac><msub><mi>P</mi><mi>S</mi></msub></mtd></mtr><mtr><mtd><mo>+</mo><mfrac><mrow><mn>2</mn><msub><mi>P</mi><mi>S</mi></msub></mrow><mrow><msub><mi>N</mi><mi>u</mi></msub><msup><mi>N</mi><mn>2</mn></msup></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>q</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mi>q</mi><mo>)</mo></mrow><msub><mi>J</mi><mn>0</mn></msub><mrow><mo>(</mo><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>max</mi></msub><msub><mi>T</mi><mi>s</mi></msub><mi>q</mi><mo>)</mo></mrow><mo>[</mo><munderover><mi>&Sigma;</mi><mrow><mi>p</mi><mo>=</mo><mo>-</mo><msub><mi>N</mi><mi>u</mi></msub><mo>+</mo><mn>1</mn></mrow><mrow><msub><mi>N</mi><mi>u</mi></msub><mo>-</mo><mn>1</mn></mrow></munderover><mrow><mo>(</mo><msub><mi>N</mi><mi>u</mi></msub><mo>-</mo><mo>|</mo><mi>p</mi><mo>|</mo><mo>)</mo></mrow><mi>cos</mi><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>&pi;q</mi><mrow><mo>(</mo><mi>p</mi><mo>-</mo><mover><mi>&epsiv;</mi><mo>^</mo></mover><mo>)</mo></mrow></mrow><mi>N</mi></mfrac><mo>)</mo></mrow><mo>-</mo><msub><mi>N</mi><mi>u</mi></msub><mi>cos</mi><mrow><mo>(</mo><mfrac><mrow><mo>-</mo><mn>2</mn><mi>&pi;q</mi><mover><mi>&epsiv;</mi><mo>^</mo></mover></mrow><mi>N</mi></mfrac><mo>)</mo></mrow><mo>]</mo></mtd></mtr></mtable></mfenced></math>]]></maths>式中,N<sub>u</sub>为OFDM/OFDMA符号中有用子载波的个数;(2b)根据补偿后信号的信干比SIR的大小确定最佳的频偏补偿倍数ρ,即将SIR取最大值时对应的频偏补偿倍数ρ作为最佳的补偿倍数。
地址 710071 陕西省西安市西安市太白南路2号