发明名称 卫星移动通信系统中对频率斜升和频偏估计的方法和装置
摘要 本发明公布了一种卫星移动通信系统中对频率斜升和频偏估计的方法和装置。所述方法包括:根据实际通信系统的参数,求出归一化频率斜升值a<sub>D</sub>的范围,并计算出频率斜升值的试探间隔Δa<sub>D</sub>;利用频偏估计方法对归一化频率斜升值a<sub>D</sub>进行多次试探,每次试探的间隔为Δa<sub>D</sub>,每次试探得到一个或多个(a<sub>D</sub>,f<sub>D</sub>)数对,当所有的试探结束后,全部(a<sub>D</sub>,f<sub>D</sub>)数对就组成了一个集合;将集合中所有的(a<sub>D</sub>,f<sub>D</sub>)数对代入似然函数G(a<sub>D</sub>,f<sub>D</sub>),找到使该似然函数取极大值时的那一组(a<sub>D</sub>,f<sub>D</sub>),作为频率斜升、频偏的估计值。本发明的具有信噪比门限低、运算复杂度低、估计精度高的优点。在实际应用中,频偏估计的方法可以沿用现成的模块,不需要改动太大便可以实现对频率斜升的估计,非常方便。
申请公布号 CN102437992B 申请公布日期 2014.04.30
申请号 CN201110414453.5 申请日期 2011.12.13
申请人 北京大学 发明人 蒋伟;乐天;甘淼;郑杨;徐国珍;罗武
分类号 H04L27/26(2006.01)I;H04L25/02(2006.01)I 主分类号 H04L27/26(2006.01)I
代理机构 北京万象新悦知识产权代理事务所(普通合伙) 11360 代理人 贾晓玲
主权项 1.一种卫星移动通信系统中对频率斜升和频偏估计的方法,包括如下步骤:1)根据实际通信系统的参数,求出归一化频率斜升值a<sub>D</sub>的范围,并计算出频率斜升值的试探间隔Δa<sub>D</sub>;2)以试探间隔Δa<sub>D</sub>对a<sub>D</sub>进行多次试探,在每个试探的a<sub>D</sub>值下,纠正a<sub>D</sub>带来的影响,接着利用频偏估计方法,每次试探得到一个或多个(a<sub>D</sub>,f<sub>D</sub>)数对,当所有的试探结束后,全部(a<sub>D</sub>,f<sub>D</sub>)数对就组成了一个集合;3)将第2)步的集合中所有的(a<sub>D</sub>,f<sub>D</sub>)数对代入似然函数G(a<sub>D</sub>,f<sub>D</sub>),找到使该似然函数取极大值时的那一组(a<sub>D</sub>,f<sub>D</sub>),作为频率斜升、频偏的估计值;其中,步骤1)中所述实际通信系统的参数为:实际通信系统的频率斜升最大值a、载波符号速率f<sub>s</sub>、发射信号的总符号数N;步骤1)中所述试探间隔Δa<sub>D</sub>的取值如下:<maths num="0001"><![CDATA[<math><mrow><mi>&Delta;</mi><msub><mi>a</mi><mi>D</mi></msub><mo>=</mo><mn>2</mn><mo>&CenterDot;</mo><mfrac><mi>&Delta;&phi;</mi><mrow><mi>&pi;</mi><mo>&CenterDot;</mo><msup><mrow><mo>(</mo><mi>N</mi><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac></mrow></math>]]></maths>其中,N为发射信号的总符号数,π为圆周率,Δφ是通信系统可容忍的最大相位偏差;步骤2)对归一化频率斜升值a<sub>D</sub>进行多次“试探”的方法为,以Δa<sub>D</sub>作试探间隔,对闭区间[-|a<sub>D</sub>|<sub>max</sub>,|a<sub>D</sub>|<sub>max</sub>]内的值进行逐一试探,所述<img file="FDA0000396763980000012.GIF" wi="410" he="191" />其中,|a<sub>D</sub>|<sub>max</sub>表示|a<sub>D</sub>|的最大值,a为实际通信系统的频率斜升最大值,f<sub>s</sub>为载波符号速率;步骤3)中所述似然函数G(a<sub>D</sub>,f<sub>D</sub>)如下得到:<maths num="0002"><![CDATA[<math><mrow><mi>G</mi><mrow><mo>(</mo><msub><mi>a</mi><mi>D</mi></msub><mo>,</mo><msub><mi>f</mi><mi>D</mi></msub><mo>)</mo></mrow><mo>=</mo><mo>|</mo><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>p</mi></msub></munderover><msub><msub><mi>z</mi><mi>k</mi></msub><mi>m</mi></msub><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><mrow><mo>(</mo><mi>&pi;</mi><msub><mi>a</mi><mi>D</mi></msub><msubsup><mi>k</mi><mi>m</mi><mn>2</mn></msubsup><mo>+</mo><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>D</mi></msub><msub><mi>k</mi><mi>m</mi></msub><mo>)</mo></mrow></mrow></msup><mo>|</mo></mrow></math>]]></maths>其中,<maths num="0003"><![CDATA[<math><mrow><msub><msub><mi>z</mi><mi>k</mi></msub><mi>m</mi></msub><mo>=</mo><msub><msub><mi>y</mi><mi>k</mi></msub><mi>m</mi></msub><mo>&CenterDot;</mo><mfrac><msubsup><mi>x</mi><msub><mi>k</mi><mi>m</mi></msub><mo>*</mo></msubsup><msqrt><msub><mi>E</mi><mi>s</mi></msub></msqrt></mfrac></mrow></math>]]></maths>其中,a<sub>D</sub>为归一化频率斜升值,f<sub>D</sub>为归一化多普勒频偏值,e为自然对数的底,j为虚数单位<img file="FDA0000396763980000022.GIF" wi="130" he="88" />E<sub>s</sub>为接收信号的能量,k<sub>m</sub>为导频符号的下标,1≤m≤N<sub>p</sub>,N<sub>p</sub>为导频符号数,<img file="FDA0000396763980000025.GIF" wi="76" he="67" />为第k<sub>m</sub>个导频符号,星号表示共轭,<img file="FDA0000396763980000026.GIF" wi="63" he="55" />为第k<sub>m</sub>个接收信号。
地址 100871 北京市海淀区颐和园路5号