发明名称 一种FFH系统中基于PN序列似然比的同步捕获方法
摘要 一种FFH系统中基于PN序列似然比的同步捕获方法,属于通信抗干扰技术领域,涉及快速跳频(Fast Frequency Hopping,FFH)通信系统,尤其是FFH通信系统的同步方法。本发明采用双图案的同步捕获头设计,充分利用了跳频频率和同步PN序列信息;采用早、迟门信号同时进行解调,可以在不增加系统实现复杂度的基础上,提高未完全同步时解调的可靠性;采用基于频率-PN序列图案似然比作为同步捕获的软判决量,减少了同步信息的损失。本发明可使FFH系统在干扰信道和衰落信道下具有更强的鲁棒性,提高了FFH系统的同步捕获性能。
申请公布号 CN101674104B 申请公布日期 2013.01.09
申请号 CN200910167726.3 申请日期 2009.09.23
申请人 电子科技大学 发明人 程郁凡;刘乔;陈亚丁;李少谦
分类号 H04B1/7156(2011.01)I 主分类号 H04B1/7156(2011.01)I
代理机构 电子科技大学专利中心 51203 代理人 葛启函
主权项 1.一种FFH系统中基于PN序列似然比的同步捕获方法,包括以下步骤:设定FFH系统发送方采用M个同步跳频频率调制长度为N的同步PN序列,M≤N,K=N/M,K为正整数,每个跳频频率调制一个同步PN-chip,M个同步跳频频率构成了一个同步跳频频率图案F={f<sub>0</sub>,f<sub>1</sub>,f<sub>i</sub>,...,f<sub>M-1</sub>},f<sub>i</sub>记为第i个同步频率,0≤i≤M-1,N个同步PN-chip构成一个同步PN序列图案S={PN<sub>0</sub>,PN<sub>1</sub>,PN<sub>i</sub>,...,PN<sub>N-1</sub>},同步跳频频率图案和同步PN序列图案共同构成跳频频率-PN序列双图案;步骤1:同步处理开始时,接收方采用与发送方相同的同步跳频频率图案,以该同步跳频频率图案中的任意一个频率作为起始频率进行解跳处理,同时同步跳数计数器的计数值C从1开始计数;步骤2:解跳时,分别对第C跳的早、迟门信号进行A/D转换,再经过快速傅氏变换和取模的平方处理,即|FFT|<sup>2</sup>处理,由于有K=N/M个相位依次相差为M的早、迟门本地同步PN序列值,故得到第C跳第n组本地同步PN序列值对应的早门基带信号的2FSK调制频点能量值<img file="FSB00000884217600011.GIF" wi="97" he="81" />和<img file="FSB00000884217600012.GIF" wi="126" he="77" />同时得到第C跳第n组本地同步PN序列值对应的迟门基带信号的2FSK调制频点能量值<img file="FSB00000884217600013.GIF" wi="94" he="84" />和<img file="FSB00000884217600014.GIF" wi="136" he="66" />其中,<img file="FSB00000884217600015.GIF" wi="93" he="95" />和<img file="FSB00000884217600016.GIF" wi="94" he="94" />分别表示第C跳第n组本地同步PN序列值对应的早门和迟门2FSK基带信号中一个调制频点,即数据信道调制频点的能量值,<img file="FSB00000884217600017.GIF" wi="102" he="90" />和<img file="FSB00000884217600018.GIF" wi="123" he="110" />分别表示第C跳第n组本地同步PN序列值对应的早门和迟门2FSK基带信号中另一个调制频点,即空闲信道调制频点的能量值;n=1,2,…,K;步骤3:根据步骤2所得到早门的<img file="FSB00000884217600019.GIF" wi="116" he="67" /><img file="FSB000008842176000110.GIF" wi="94" he="68" />与迟门的<img file="FSB000008842176000111.GIF" wi="118" he="110" /><img file="FSB000008842176000112.GIF" wi="118" he="90" />分别计算第C跳时,第n组本地同步PN序列值对应的早、迟门同步PN-chip似然比<img file="FSB000008842176000113.GIF" wi="75" he="73" />和<img file="FSB000008842176000114.GIF" wi="142" he="89" /><maths num="0001"><![CDATA[<math><mrow><msubsup><mi>v</mi><mi>z</mi><mi>nC</mi></msubsup><mo>=</mo><mi>ln</mi><mo>[</mo><mfrac><mrow><mi>f</mi><mrow><mo>(</mo><msubsup><mi>X</mi><mrow><mi>z</mi><mn>1</mn></mrow><mi>nC</mi></msubsup><mo>|</mo><msub><mi>H</mi><mn>1</mn></msub><mo>)</mo></mrow><mi>g</mi><mrow><mo>(</mo><msubsup><mi>X</mi><mrow><mi>z</mi><mn>0</mn></mrow><mi>nC</mi></msubsup><mo>|</mo><msub><mi>H</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow><mrow><mi>f</mi><mrow><mo>(</mo><msubsup><mi>X</mi><mrow><mi>z</mi><mn>0</mn></mrow><mi>nC</mi></msubsup><mo>|</mo><msub><mi>H</mi><mn>1</mn></msub><mo>)</mo></mrow><mi>g</mi><mrow><mo>(</mo><msubsup><mi>X</mi><mrow><mi>z</mi><mn>1</mn></mrow><mi>nC</mi></msubsup><mo>|</mo><msub><mi>H</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow></mfrac><mo>]</mo></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msubsup><mi>v</mi><mi>c</mi><mi>nC</mi></msubsup><mo>=</mo><mi>ln</mi><mo>[</mo><mfrac><mrow><mi>f</mi><mrow><mo>(</mo><msubsup><mi>X</mi><mrow><mi>c</mi><mn>1</mn></mrow><mi>nC</mi></msubsup><mo>|</mo><msub><mi>H</mi><mn>1</mn></msub><mo>)</mo></mrow><mi>g</mi><mrow><mo>(</mo><msubsup><mi>X</mi><mrow><mi>c</mi><mn>0</mn></mrow><mi>nC</mi></msubsup><mo>|</mo><msub><mi>H</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow><mrow><mi>f</mi><mrow><mo>(</mo><msubsup><mi>X</mi><mrow><mi>c</mi><mn>0</mn></mrow><mi>nC</mi></msubsup><mo>|</mo><msub><mi>H</mi><mn>1</mn></msub><mo>)</mo></mrow><mi>g</mi><mrow><mo>(</mo><msubsup><mi>X</mi><mrow><mi>c</mi><mn>1</mn></mrow><mi>nC</mi></msubsup><mo>|</mo><msub><mi>H</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow></mfrac><mo>]</mo></mrow></math>]]></maths>其中:f(x|H<sub>1</sub>)为理想情况时,即收发端跳频频率完全对准、且时间窗口也完全对准时,2FSK信号对应调制频点能量值的概率密度函数;g(x|H<sub>0</sub>)为2FSK信号另一个调制频点能量值的概率密度函数;当C≥N时,舍去前C-N跳同步PN-chip似然比值,由第n组本地同步PN序列值对应的当前最新N跳信号得到的N个早门同步PN-chip似然比,按顺序依次构成第n组早门PN似然比序列<img file="FSB00000884217600021.GIF" wi="520" he="96" />相应的N个迟门同步PN-chip似然比,按顺序依次构成第n组迟门PN似然比序列<img file="FSB00000884217600022.GIF" wi="580" he="133" />步骤4:分别计算K组早、迟门同步PN似然比序列<img file="FSB00000884217600023.GIF" wi="491" he="118" />和<img file="FSB00000884217600024.GIF" wi="570" he="138" />的代数和,得到K个早门同步信号似然比相关检测值<img file="FSB00000884217600025.GIF" wi="228" he="168" />和K个迟门同步信号似然比相关检测值<img file="FSB00000884217600026.GIF" wi="242" he="142" />n=1,2,…,K;步骤5:早、迟门分别进行K路似然比相关检测值门限判决;判断由步骤4所得的2K个同步信号似然比相关检测值中是否存在大于或等于同步捕获门限ζ<sub>Syn-th</sub>的同步信号似然比相关检测值;若存在,则转入步骤6;若不存在且C<mN,m为正整数,则在下一跳同步跳数计数器C增加1的同时返回步骤2;若不存在且C≥mN,则在下一跳将同步跳数计数器C清零后重新开始从1计数,并更换本地初始频率,返回步骤2;步骤6:进行同步PN序列定位;根据步骤1中采用的起始频率的频率号i,0≤i≤M-1,并根据出现大于或等于同步捕获门限ζ<sub>Syn-th</sub>的同步信号似然比相关检测值所在的同步PN似然比序列中的n值和同步跳数计数器的计数值C,计算接收方需要等待的跳数W;接收方等待W跳后,能使下一跳收发双方频率在第0号频率f<sub>0</sub>上对准,并且PN序列相位偏移为0的位置;其中,等待跳数W的计算公式为:W=(m×N-i-C-n×M)%N,其中符号%表示求模运算;步骤7:进行捕获验证;如果连续V个同步PN似然比序列都检测到大于或等于同步捕获门限ζ<sub>Syn-th</sub>的同步信号似然比相关检测值,则同步捕获成功,系统转入同步跟踪状态;如果连续V个同步PN似然比序列中任何一个未能检测到大于或等于同步捕获门限ζ<sub>Syn-th</sub>的同步信号似然比相关检测值,则下一跳更换本地初始频率,则转入步骤2重新进行同步捕获。
地址 611731 四川省成都市高新区(西区)西源大道2006号