发明名称 一种无人机数据链路频偏检测方法
摘要 本发明提供了一种无人机数据链路频偏检测方法,包括数据链路系统功能初始化、小数倍频偏值计算、整数倍频偏值计算和频偏值有效性判别的步骤。本发明提出的频偏检测方法,采用PN码作为同步信息序列,并利用通信系统原有模块,进行同步检测,不再占用新的资源。
申请公布号 CN106549893A 申请公布日期 2017.03.29
申请号 CN201510589532.8 申请日期 2015.09.16
申请人 中国航空工业第六一八研究所 发明人 任继军
分类号 H04L27/00(2006.01)I 主分类号 H04L27/00(2006.01)I
代理机构 中国航空专利中心 11008 代理人 杜永保
主权项 一种无人机数据链路频偏检测方法,针对以QPSK调制方式为主的数据链进行频偏检测,包括以下步骤:步骤一:数据链路系统功能初始化;在数据链路系统的发射端,将要发射的信号分成I,Q两路,然后将I,Q两路信号分别进行组帧,信号帧的结构为:同步信息序列,之后依次重复为帧起始、数据区、校验信息、帧结束;I,Q两路的同步信息序列分别选择N位的PN序列,N为2的幂次且N≥10;步骤二:小数倍频偏值计算;设发射端的同步信息时域序列为m<sub>t</sub>(n),n为整数且n≥1,周期为<img file="FDA0000804112960000011.GIF" wi="214" he="117" />K为大于等于2的整数,[ ]表示按四舍五入准则取整操作。根据公式<img file="FDA0000804112960000012.GIF" wi="119" he="71" /><img file="FDA0000804112960000013.GIF" wi="510" he="85" />计算在发射端同步信息序列之间的互相关值,并把最大互相关值所对应的相位作为发射相位θ<sub>t</sub>;对接收到的信息数据序列进行解调和译码得到接收端的同步信息时域序列m<sub>r</sub>(n),n为整数且n≥1。设同步信息序列的起始帧头位置为n<sub>0</sub>,根据公式<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>R</mi><mi>r</mi></msub><mo>=</mo><msubsup><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>T</mi><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mi>m</mi><mi>r</mi></msub><mrow><mo>(</mo><msub><mi>n</mi><mn>0</mn></msub><mo>+</mo><mi>i</mi><mo>)</mo></mrow><msup><msub><mi>m</mi><mi>r</mi></msub><mo>*</mo></msup><mrow><mo>(</mo><msub><mi>n</mi><mn>0</mn></msub><mo>+</mo><mi>i</mi><mo>+</mo><mi>T</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000804112960000014.GIF" wi="936" he="87" /></maths>计算在接收端同步信息序列之间的互相关值,并把最大互相关值所对应的相位作为接收相位θ<sub>r</sub>;通过比较接收相位θ<sub>r</sub>和发射相位θ<sub>t</sub>之间的差别,根据公式<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><mi>&epsiv;</mi><mo>^</mo></mover><mo>=</mo><mfrac><mi>K</mi><mrow><mn>2</mn><mi>&pi;</mi></mrow></mfrac><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>r</mi></msub><mo>-</mo><msub><mi>&theta;</mi><mi>t</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000804112960000015.GIF" wi="421" he="111" /></maths>来计算小数倍频偏值<img file="FDA0000804112960000016.GIF" wi="54" he="70" />步骤三:整数倍频偏值计算;根据计算出来的频点小数倍频偏值对接收信号y(n)进行频偏补偿,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msup><msub><mi>m</mi><mi>r</mi></msub><mo>*</mo></msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>m</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msup><mi>e</mi><mrow><mi>j</mi><mover><mi>&epsiv;</mi><mo>^</mo></mover></mrow></msup></mrow>]]></math><img file="FDA0000804112960000017.GIF" wi="493" he="86" /></maths>然后对频偏补偿后的信号m<sub>t</sub><sup>*</sup>(n)进行FFT变换,转换成频域序列M<sub>r</sub><sup>*</sup>(n);M<sub>r</sub><sup>*</sup>(n)=FFT(m<sub>r</sub><sup>*</sup>(n))对发射端同步信息序列m<sub>t</sub>(n)进行FFT变换,转换成频域序列M<sub>t</sub>(n);M<sub>t</sub>(n)=FFT(m<sub>t</sub>(n))FFT点数为N,N选取128、512、1024、2048这些2的幂次值;将M<sub>r</sub><sup>*</sup>(n)与M<sub>t</sub>(n)进行互相关运算:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>C</mi><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>M</mi><mi>t</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>+</mo><mi>f</mi><mo>)</mo></mrow><msup><msub><mi>M</mi><mi>r</mi></msub><mo>*</mo></msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><mo>-</mo><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>&le;</mo><mi>f</mi><mo>&le;</mo><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000804112960000021.GIF" wi="1318" he="222" /></maths>根据计算得到的最大互相关值所对应的频率作为整数倍频偏;步骤四:频偏值有效性判别;将接收端的I,Q两路数据分别进行各自的小数倍和整数倍频偏值补偿,根据得到的频偏补偿后的I<sup>*</sup>、Q<sup>*</sup>两路数据进行能量计算,I<sup>*2</sup>+Q<sup>*2</sup>,以此作为频偏值有效性的判别,如果(I<sup>*2</sup>+Q<sup>*2</sup>)&gt;(I<sup>2</sup>+Q<sup>2</sup>),则此频偏值有效,将以此作为接收信号的频偏补偿值;反之,则判为无效,重新进行迭代计算。
地址 710065 陕西省西安市雁塔区电子一路92号