发明名称 一种基于共轭结构前导的联合定时与频率同步方法
摘要 本发明提供了一种基于共轭结构前导的联合定时与频率同步方法,该方法的步骤如下:1、计算接收信号r(n)和前导符号c(n)的相关信号r<sub>0</sub>(n,d);2、计算得到差分相关信号p(m,d);3、对差分相关信号p(m,d)进行加权累加得到加权相关函数P(d);4、利用滑动窗内接收信号的能量对加权相关函数P(d)进行归一化,得到归一化定时度量值M(d);5、对归一化定时度量M(d)求取最大值,得到定时偏移估计量<img file="DDA0000710000250000011.GIF" wi="432" he="88" />6、计算小数部分载波频偏估计量<img file="DDA0000710000250000012.GIF" wi="56" he="72" />7、利用小数部分的载波频偏估计量<img file="DDA0000710000250000013.GIF" wi="40" he="66" />对接收信号进行补偿得到补偿后的接收信号 <img file="DDA0000710000250000014.GIF" wi="176" he="64" />8、利用小数频偏补偿后的接收信号计算整数部分的频偏估计量<img file="DDA0000710000250000015.GIF" wi="56" he="56" />9、载波频偏估计量<img file="DDA0000710000250000016.GIF" wi="208" he="72" />该方法只采用一个具有共轭结构的前导符号进行联合符号定时与载波频率同步,符号定时同步不受频率偏差的影响。
申请公布号 CN104836769A 申请公布日期 2015.08.12
申请号 CN201510217786.7 申请日期 2015.04.30
申请人 北京遥测技术研究所;航天长征火箭技术有限公司 发明人 马云思;周三文;闫朝星
分类号 H04L27/26(2006.01)I;H04L25/02(2006.01)I 主分类号 H04L27/26(2006.01)I
代理机构 中国航天科技专利中心 11009 代理人 范晓毅
主权项 一种基于共轭结构前导的联合定时与频率同步方法,其特征在于包括如下步骤:(1)、在长度为N的滑动窗内,对OFDM接收机的接收信号r(n)和前导符号c(n)进行共轭相乘,得到相关信号r<sub>0</sub>(n,d):r<sub>0</sub>(n,d)=r(n+d)c<sup>*</sup>(n);其中,n=0、1、…、N‑1;d=0、1、…、M<sub>s</sub>×N<sub>s</sub>;N<sub>s</sub>=N+N<sub>g</sub>,N为设定的一个OFDM符号的样值数,M<sub>s</sub>为设定的每帧数据中OFDM符号个数,N<sub>g</sub>为OFDM符号的循环前缀个数;并且在所述前导符号c(n)具有共轭结构,即:[c(0) c(1) … c(N/2‑1)]=[c<sup>*</sup>(N/2) c<sup>*</sup>(N/2+1) … c<sup>*</sup>(N‑1)];(2)、根据步骤(1)得到的滑动相关信号r<sub>0</sub>(n,d)计算得到差分相关信号p(m,d):<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>p</mi><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>d</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mi>m</mi></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>r</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>d</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msup><msub><mi>r</mi><mn>0</mn></msub><mo>*</mo></msup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mi>m</mi><mo>,</mo><mi>d</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000710000220000011.GIF" wi="722" he="142" /></maths>其中,m=1、2、…、M<sub>0</sub>,M<sub>0</sub>为正整数且M<sub>0</sub>≤N‑1;(3)、对步骤(2)计算得到的差分相关信号p(m,d)进行加权累加,得到加权相关函数P(d);(4)、利用滑动窗内接收信号的能量对步骤(3)计算得到的加权相关函数P(d)进行归一化,得到归一化定时度量值M(d);(5)、对归一化定时度量M(d)求取最大值;差分相关定时偏移估计量<img file="FDA0000710000220000012.GIF" wi="46" he="62" />等于所述最大值对应的d的取值,即<img file="FDA0000710000220000013.GIF" wi="442" he="102" />(6)、计算小数部分载波频偏估计量<img file="FDA0000710000220000014.GIF" wi="70" he="80" /><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><mi>&xi;</mi><mo>^</mo></mover><mo>=</mo><mfrac><mn>1</mn><mi>&pi;</mi></mfrac><mi>angle</mi><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>/</mo><mn>2</mn><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>r</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mover><mi>&epsiv;</mi><mo>^</mo></mover><mo>)</mo></mrow><mo>&CenterDot;</mo><msubsup><mi>r</mi><mn>0</mn><mo>*</mo></msubsup><mrow><mo>(</mo><mi>n</mi><mo>+</mo><mfrac><mi>N</mi><mn>2</mn></mfrac><mo>,</mo><mover><mi>&epsiv;</mi><mo>^</mo></mover><mo>)</mo></mrow><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000710000220000015.GIF" wi="834" he="153" /></maths>其中,angle()代表取相位运算;(7)、利用步骤(6)计算得到的小数部分的载波频偏估计量<img file="FDA0000710000220000021.GIF" wi="46" he="89" />对接收信号进行补偿运算,得到小数频偏补偿后的接收信号<img file="FDA0000710000220000022.GIF" wi="175" he="74" /><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>r</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mover><mi>&epsiv;</mi><mo>^</mo></mover><mo>)</mo></mrow><mo>=</mo><msub><mi>r</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mover><mi>&epsiv;</mi><mo>^</mo></mover><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>j</mi><mn>2</mn><mi>&pi;n</mi><mover><mi>&xi;</mi><mo>^</mo></mover><mo>/</mo><mi>N</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000710000220000023.GIF" wi="738" he="92" /></maths>(8)、利用小数频偏补偿后的接收信号<img file="FDA0000710000220000024.GIF" wi="148" he="76" />求取目标函数Γ(q),然后对所述目标函数Γ(q)求取最大值,并将所述最大值对应的参数q作为整数部分的频偏估计量<img file="FDA0000710000220000025.GIF" wi="68" he="70" />具体实现过程如下:(8a)、计算目标函数Γ(q):<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mi>q</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>|</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>r</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mover><mi>&epsiv;</mi><mo>^</mo></mover><mo>)</mo></mrow><mo>&CenterDot;</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><mn>2</mn><mi>&pi;qn</mi><mo>/</mo><mi>N</mi></mrow></msup><mo>|</mo></mrow><mn>2</mn></msup><mo>,</mo><mi>q</mi><mo>=</mo><mo>-</mo><mfrac><mi>N</mi><mn>2</mn></mfrac><mo>,</mo><mo>-</mo><mfrac><mi>N</mi><mn>2</mn></mfrac><mo>+</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mfrac><mi>N</mi><mn>2</mn></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000710000220000026.GIF" wi="1312" he="162" /></maths>(8b)、然后对所述目标函数Γ(q)求取最大值,并将所述最大值对应的参数q作为整数部分的频偏<img file="FDA0000710000220000027.GIF" wi="56" he="56" />即<img file="FDA0000710000220000028.GIF" wi="416" he="100" />(9)、将小数部分的载波频偏估计量<img file="FDA0000710000220000029.GIF" wi="48" he="86" />和整数部分的频偏估计量<img file="FDA00007100002200000210.GIF" wi="44" he="75" />求和,得到载波频偏估计量<img file="FDA00007100002200000211.GIF" wi="64" he="68" />即<img file="FDA00007100002200000212.GIF" wi="219" he="84" />
地址 100076 北京市丰台区北京9200信箱74分箱