发明名称 一种通用的非相干直接序列扩频信号跟踪方法
摘要 一种通用的非相干直接序列扩频信号跟踪方法,步骤如下:一:根据待跟踪的非相干直接序列扩频信号特征确定相干积分时间,并生成相干积分脉冲:二:计算需要配置的参数,计有载波中心频率、数据速率、扩频码速率和环路滤波器系数;三:数字下变频:输入的模拟中频信号经采样后,变为数字中频信号并与本地复现的数字中频载波相乘即可实现数字下变频;四:数字下变频结果与扩频码做相关运算:五:相干积分结果鉴频鉴相及滤波:相干积分结果I<sub>P</sub>、Q<sub>P</sub>用于载波环鉴频鉴相,鉴频及鉴相结果输出给二阶辅助三阶的载波环环路滤波器,滤波结果反馈控制数据NCO;通过上述五个步骤达到对扩频码速率和数据速率可变的非相干直接序列扩频信号的通用跟踪。
申请公布号 CN103414493A 申请公布日期 2013.11.27
申请号 CN201310334670.2 申请日期 2013.08.02
申请人 北京航空航天大学 发明人 常青;李育龙;徐勇;辜晓波
分类号 H04B1/7085(2011.01)I 主分类号 H04B1/7085(2011.01)I
代理机构 北京慧泉知识产权代理有限公司 11232 代理人 王顺荣;唐爱华
主权项 1.一种通用的非相干直接序列扩频信号跟踪方法,其特征在于:该方法包含如下步骤:步骤一:根据待跟踪的非相干直接序列扩频信号特征确定相干积分时间,并生成相干积分脉冲:相干积分脉冲有两类,一类用于确定载波环及码环中的相关器的相干积分起始,一类用于确定位同步环的相关器的相干积分起始;由于数据跳变的影响,载波环及码环的相关器相干积分时间应小于二分之一数据比特持续时间,而为了提高信噪比,相干积分时间应尽量长,综合考虑,在数据速率较高时,确定相干积分时间为半个数据比特持续时间;捕获提供给跟踪的多普勒范围为±250Hz,由于数据跳变,鉴频器鉴相范围为±1/(4*T),所以相干积分时间T最大取值为0.5ms;位同步环相干积分起始分为同相和中相脉冲,同相脉冲在位同步环锁定时,理论上与数据跳变沿对齐,而中相脉冲处于两同相脉冲正中间,即数据中点;数据速率为1Kbps时,半个数据比特持续时间是0.5ms,所以当数据速率大于等于1Kbps时,相关器相干积分脉冲等于同相脉冲和中相脉冲的逻辑或;而数据速率小于1Kbps时,不能简单的把载波环及码环中的相关器相干积分时间定为0.5ms,因为在数据速率大于250bps时,以固定的0.5ms为相干积分时间,必定会有些相干积分结果是跨越数据比特的,这对之后的锁频环的影响很大,很可能会失锁;所以设定数据速率大于等于500bps且小于1Kbps时,相干积分时间取为四分之一数据比特持续时间,积分脉冲与位同步的同相脉冲和中相脉冲有关系,即在同相脉冲和中相脉冲正中间又加入一个脉冲;同理数据速率大于等于250bps且小于500bps时,相干积分时间取为八分之一数据比特持续时间;数据速率小于250bps时,相干积分时间取为0.5ms,这种情况下,积分脉冲与位定时脉冲完全独立,即可确定载波环及码环的相关器相干积分时间如下:载波环及码环的相关器相干积分时间<maths num="0001"><![CDATA[<math><mrow><mi>T</mi><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>T</mi><mi>b</mi></msub></mtd><mtd><msub><mi>R</mi><mi>b</mi></msub><mo>></mo><mo>=</mo><mn>1</mn><mi>Kbps</mi></mtd></mtr><mtr><mtd><mfrac><mn>1</mn><mn>4</mn></mfrac><msub><mi>T</mi><mi>b</mi></msub></mtd><mtd><mn>500</mn><mi>bps</mi><mo>=</mo><mo>&lt;</mo><msub><mi>R</mi><mi>b</mi></msub><mo>&lt;</mo><mn>1</mn><mi>Kbps</mi></mtd></mtr><mtr><mtd><mfrac><mn>1</mn><mn>8</mn></mfrac><msub><mi>T</mi><mi>b</mi></msub></mtd><mtd><mn>250</mn><mi>bps</mi><mo>=</mo><mo>&lt;</mo><msub><mi>R</mi><mi>b</mi></msub><mo>&lt;</mo><mn>500</mn><mi>bps</mi></mtd></mtr><mtr><mtd><mn>0.5</mn><mi>ms</mi></mtd><mtd><msub><mi>R</mi><mi>b</mi></msub><mo>&lt;</mo><mn>250</mn><mi>bps</mi></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其中T<sub>b</sub>为数据比特持续时间,R<sub>b</sub>为数据速率,且<img file="FDA00003613785900022.GIF" wi="194" he="145" />综上,相干积分脉冲生成方法如下,根据数据速率,生成同相脉冲和中相脉冲,之后再由数据速率,确定载波环及码环的相关器相干积分时间T,即确定了载波环及码环的相关器的相干积分脉冲周期为T,并注意到当T小于0.5ms时,载波环及码环的相关器相干积分脉冲与同相脉冲有对应的相位关系,由此生成载波环及码环的相关器相干积分脉冲;步骤二:计算需要配置的参数为了能够跟踪锁定不同扩频码速率、不同数据速率的非相干直接序列扩频信号,需要根据相干积分时间计算出跟踪所需的配置参数;需要配置的参数有载波中心频率、数据速率、扩频码速率、环路滤波器系数;载波中心频率要换算为载波频率控制字,计算公式如下:<maths num="0002"><![CDATA[<math><mrow><msub><mi>FW</mi><mi>carr</mi></msub><mo>=</mo><msub><mi>f</mi><mi>carr</mi></msub><mo>*</mo><mfrac><msup><mn>2</mn><msub><mi>W</mi><mi>f</mi></msub></msup><msub><mi>f</mi><mi>s</mi></msub></mfrac></mrow></math>]]></maths>其中FW<sub>carr</sub>为载波频率控制字,f<sub>carr</sub>为载波频率,f<sub>s</sub>为系统采样率,W<sub>f</sub>为载波NCO位宽;数据速率要换算为数据NCO控制字,计算公式如下:<maths num="0003"><![CDATA[<math><mrow><msub><mi>FW</mi><mi>data</mi></msub><mo>=</mo><msub><mi>f</mi><mi>data</mi></msub><mo>*</mo><mfrac><msup><mn>2</mn><msub><mi>W</mi><mi>d</mi></msub></msup><msub><mi>f</mi><mi>s</mi></msub></mfrac></mrow></math>]]></maths>其中FW<sub>data</sub>为数据NCO控制字,f<sub>data</sub>为数据速率,f<sub>s</sub>为系统采样率,W<sub>d</sub>为数据NCO位宽;扩频码速率要换算为扩频码NCO控制字,计算公式如下:<maths num="0004"><![CDATA[<math><mrow><msub><mi>FW</mi><mi>code</mi></msub><mo>=</mo><msub><mi>f</mi><mi>code</mi></msub><mo>*</mo><mfrac><msup><mn>2</mn><msub><mi>W</mi><mi>c</mi></msub></msup><msub><mi>f</mi><mi>s</mi></msub></mfrac></mrow></math>]]></maths>其中FW<sub>code</sub>为扩频码NCO控制字,f<sub>code</sub>为扩频码速率,f<sub>s</sub>为系统采样率,W<sub>c</sub>为扩频码NCO位宽;环路滤波器系数由环路自然谐振角频率、相干积分时间、系统采样率和量化倍数共同决定;载波环环路滤波器系数计算公式如下:锁频环系数:<maths num="0005"><![CDATA[<math><mrow><msub><mi>C</mi><mrow><mi>f</mi><mn>1</mn></mrow></msub><mo>=</mo><msubsup><mi>&omega;</mi><mrow><mn>0</mn><mi>f</mi></mrow><mn>2</mn></msubsup><mo>*</mo><mi>T</mi><mo>*</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>&pi;</mi></mrow></mfrac><mo>*</mo><mfrac><msup><mn>2</mn><msub><mi>W</mi><mi>f</mi></msub></msup><msub><mi>f</mi><mi>s</mi></msub></mfrac><mo>*</mo><msup><mn>2</mn><msub><mi>N</mi><mi>f</mi></msub></msup></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><msub><mi>C</mi><mrow><mi>f</mi><mn>2</mn></mrow></msub><mo>=</mo><msub><mi>a</mi><mn>2</mn></msub><msub><mi>&omega;</mi><mrow><mn>0</mn><mi>f</mi></mrow></msub><mo>*</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>&pi;</mi></mrow></mfrac><mo>*</mo><mfrac><msup><mn>2</mn><msub><mi>W</mi><mi>f</mi></msub></msup><msub><mi>f</mi><mi>s</mi></msub></mfrac><mo>*</mo><msup><mn>2</mn><msub><mi>N</mi><mi>f</mi></msub></msup></mrow></math>]]></maths>其中a<sub>2</sub>=1.414,B<sub>nf</sub>=0.53ω<sub>0f</sub>,ω<sub>0f</sub>为环路自然谐振角频率,T为相干积分时间,f<sub>s</sub>为系统采样率,W<sub>f</sub>为载波NCO位宽,N<sub>f</sub>为系数量化位数;锁相环系数:<maths num="0007"><![CDATA[<math><mrow><msub><mi>C</mi><mrow><mi>p</mi><mn>1</mn></mrow></msub><mo>=</mo><msubsup><mi>&omega;</mi><mrow><mn>0</mn><mi>p</mi></mrow><mn>3</mn></msubsup><mo>*</mo><msup><mi>T</mi><mn>2</mn></msup><mo>*</mo><mfrac><msup><mn>2</mn><msub><mi>W</mi><mi>p</mi></msub></msup><msub><mi>f</mi><mi>s</mi></msub></mfrac><mo>*</mo><msup><mn>2</mn><msub><mi>N</mi><mi>p</mi></msub></msup></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><msub><mi>C</mi><mrow><mi>p</mi><mn>2</mn></mrow></msub><mo>=</mo><msub><mi>a</mi><mn>3</mn></msub><msubsup><mi>&omega;</mi><mrow><mn>0</mn><mi>p</mi></mrow><mn>2</mn></msubsup><mo>*</mo><mi>T</mi><mo>*</mo><mfrac><msup><mn>2</mn><msub><mi>W</mi><mi>p</mi></msub></msup><msub><mi>f</mi><mi>s</mi></msub></mfrac><mo>*</mo><msup><mn>2</mn><msub><mi>N</mi><mi>p</mi></msub></msup></mrow></math>]]></maths><maths num="0009"><![CDATA[<math><mrow><msub><mi>C</mi><mrow><mi>p</mi><mn>3</mn></mrow></msub><mo>=</mo><msub><mi>b</mi><mn>3</mn></msub><msub><mi>&omega;</mi><mrow><mn>0</mn><mi>p</mi></mrow></msub><mo>*</mo><mfrac><msup><mn>2</mn><msub><mi>W</mi><mi>p</mi></msub></msup><msub><mi>f</mi><mi>s</mi></msub></mfrac><mo>*</mo><msup><mn>2</mn><msub><mi>N</mi><mi>p</mi></msub></msup></mrow></math>]]></maths>其中a<sub>3</sub>=1.1,b<sub>3</sub>=2.4,B<sub>np</sub>=0.7845ω<sub>0p</sub>,ω<sub>0p</sub>为环路自然谐振角频率,T为相干积分时间,f<sub>s</sub>为系统采样率,W<sub>p</sub>为载波NCO位宽,N<sub>p</sub>为系数量化位数;码环环路滤波器系数计算公式如下:<maths num="0010"><![CDATA[<math><mrow><msub><mi>C</mi><mrow><mi>c</mi><mn>1</mn></mrow></msub><mo>=</mo><mfrac><mi>D</mi><mn>2</mn></mfrac><msubsup><mi>&omega;</mi><mrow><mn>0</mn><mi>c</mi></mrow><mn>2</mn></msubsup><mo>*</mo><mi>T</mi><mo>*</mo><mfrac><msup><mn>2</mn><msub><mi>W</mi><mi>c</mi></msub></msup><msub><mi>f</mi><mi>s</mi></msub></mfrac><mo>*</mo><msup><mn>2</mn><msub><mi>N</mi><mi>c</mi></msub></msup></mrow></math>]]></maths><maths num="0011"><![CDATA[<math><mrow><msub><mi>C</mi><mrow><mi>c</mi><mn>2</mn></mrow></msub><mo>=</mo><mfrac><mi>D</mi><mn>2</mn></mfrac><msub><mi>a</mi><mn>2</mn></msub><msub><mi>&omega;</mi><mrow><mn>0</mn><mi>c</mi></mrow></msub><mo>*</mo><mfrac><msup><mn>2</mn><msub><mi>W</mi><mi>c</mi></msub></msup><msub><mi>f</mi><mi>s</mi></msub></mfrac><mo>*</mo><msup><mn>2</mn><msub><mi>N</mi><mi>c</mi></msub></msup></mrow></math>]]></maths>其中a<sub>2</sub>=1.414,B<sub>nc</sub>=0.53ω<sub>0c</sub>,ω<sub>0c</sub>为环路自然谐振角频率,T为相干积分时间,f<sub>s</sub>为系统采样率,W<sub>c</sub>为扩频码NCO位宽,N<sub>c</sub>为系数量化位数,D为码间距;位同步环环路滤波器系数计算公式如下:<maths num="0012"><![CDATA[<math><mrow><msub><mi>C</mi><mrow><mi>d</mi><mn>1</mn></mrow></msub><mo>=</mo><msubsup><mi>&omega;</mi><mrow><mn>0</mn><mi>d</mi></mrow><mn>2</mn></msubsup><mo>*</mo><mi>T</mi><mo>*</mo><mfrac><msup><mn>2</mn><msub><mi>W</mi><mi>d</mi></msub></msup><msub><mi>f</mi><mi>s</mi></msub></mfrac><mo>*</mo><msup><mn>2</mn><msub><mi>N</mi><mi>d</mi></msub></msup></mrow></math>]]></maths><maths num="0013"><![CDATA[<math><mrow><msub><mi>C</mi><mrow><mi>d</mi><mn>2</mn></mrow></msub><mo>=</mo><msub><mi>&alpha;</mi><mn>2</mn></msub><msub><mi>&omega;</mi><mrow><mn>0</mn><mi>d</mi></mrow></msub><mo>*</mo><mfrac><msup><mn>2</mn><msub><mi>W</mi><mi>d</mi></msub></msup><msub><mi>f</mi><mi>s</mi></msub></mfrac><mo>*</mo><msup><mn>2</mn><msub><mi>N</mi><mi>d</mi></msub></msup></mrow></math>]]></maths>其中a<sub>2</sub>=1.414,B<sub>nd</sub>=0.53ω<sub>0d</sub>,ω<sub>0d</sub>为环路自然谐振角频率,T为相干积分时间,f<sub>s</sub>为系统采样率,W<sub>d</sub>为数据NCO位宽,N<sub>d</sub>为系数量化位数;步骤三:数字下变频输入的模拟中频信号经模数转换器ADC采样后,变为数字中频信号,数字中频信号与本地复现的数字中频载波相乘即可实现数字下变频;数字中频信号与本地载波NCO控制产生的sin和cos相乘,完成数字下变频,得到两路下变频结果,分别是同相I路的I<sub>down</sub>和正交Q路Q<sub>down</sub>;步骤四:数字下变频结果与扩频码做相关运算扩频码频率控制字控制扩频码NCO生成扩频码时钟,扩频码时钟驱动扩频码发生器生成超前E支路扩频码,E支路扩频码经延时后得到即时P支路扩频码和滞后L支路扩频码;数字下变频结果I<sub>down</sub>和Q<sub>down</sub>,分别根据步骤一生成的载波环及码环的相关器相干积分脉冲与本地扩频码发生器的超前E、即时P、滞后L支路扩频码做相关运算,输出结果为I<sub>E</sub>、I<sub>P</sub>、I<sub>L</sub>、Q<sub>E</sub>、Q<sub>P</sub>、Q<sub>L</sub>;I路下变频结果根据同相脉冲与即时P支路扩频码做相关运算,输出结果为I<sub>edge</sub>;I路下变频结果根据中相脉冲与即时P支路扩频码做相关运算,输出结果为I<sub>half</sub>;其中,I<sub>E</sub>为I支路超前相干积分结果;I<sub>P</sub>为I支路即时相干积分结果;I<sub>L</sub>为I支路滞后相干积分结果;Q<sub>E</sub>为Q支路超前相干积分结果;Q<sub>P</sub>为Q支路即时相干积分结果;Q<sub>L</sub>为Q支路滞后相干积分结果;I<sub>edge</sub>为同相脉冲支路相干积分果;I<sub>half</sub>为中相脉冲支路相干积分果;步骤五:相干积分结果鉴频鉴相及滤波相干积分结果I<sub>P</sub>、Q<sub>P</sub>用于载波环鉴频鉴相,鉴频及鉴相结果输出给二阶辅助三阶的载波环环路滤波器,滤波结果反馈控制步骤三中的本地载波NCO;相干积分结果I<sub>E</sub>、I<sub>L</sub>、Q<sub>E</sub>、Q<sub>L</sub>用于码环鉴相,鉴相结果输出给码环环路滤波器,滤波结果反馈控制步骤四的本地扩频码NCO,相干积分结果I<sub>edge</sub>、I<sub>half</sub>用于位同步环鉴相,鉴相结果输出给位同步环环路滤波器,位同步环环路滤波器与码环环路滤波器相同,滤波结果反馈控制数据NCO;通过上述五个步骤达到对扩频码速率和数据速率可变的非相干直接序列扩频信号的通用跟踪。
地址 100191 北京市海淀区学院路37号