[发明专利]机群链路通信/测量综合信道体制的信号处理系统有效
申请号: | 201010103897.2 | 申请日: | 2010-01-29 |
公开(公告)号: | CN101789809A | 公开(公告)日: | 2010-07-28 |
发明(设计)人: | 杨宜康;陈晓敏;安军社 | 申请(专利权)人: | 中国科学院空间科学与应用研究中心 |
主分类号: | H04B1/707 | 分类号: | H04B1/707;H04B17/00;H04L25/02;H04L27/18 |
代理公司: | 北京慧泉知识产权代理有限公司 11232 | 代理人: | 王顺荣;唐爱华 |
地址: | 100190 *** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明一种机群链路通信/测量综合信道体制的信号处理系统,可以在电路板的数字信号处理器DSP和FPGA上实现的单信道内非扩频通信数传/扩频测距体制的体系构架。本发明保留非扩频体制信息传输率高的优点,吸纳扩频体制方便测距、抗干扰的优点,通过同步时分复用信道的方式在单信道内综合实现非扩频通信数传/扩频测距综合体制;同时,本发明提供了基于载波多普勒频移的实时预报算法的通信/测量综合信道体制的的信号跟踪方法,解决了高动态条件下通信/测量综合信道体制的信号跟踪难题。本发明公开的方法能够广泛应用于基于抑制载波调制直接序列扩频体制的卫星导航接收机、测距系统、通信/测距功能综合化系统。 | ||
搜索关键词: | 机群 通信 测量 综合 信道 体制 信号 处理 系统 | ||
【主权项】:
1.一种机群链路通信/测量综合信道体制的信号处理系统,针对机群编队任务的自主性特点,构造了基于扩频测量/非扩频数传等时复用物理信道的机群链路通信与测量综合信道体制及其实现方法,其特征在于:该信号处理系统具体包括信号生成模块及信号跟踪模块;具体技术方案如下:(一)、机群链路通信/测量综合信道体制的信号生成模块该模块是在统一信道中实现扩频测距与高速通信数传的综合信道体制——扩频测量/非扩频数传等时复用信道体制,在单载波物理信道中利用等时复用原理,交替传输低码率扩频基带信号和高码率非扩频基带信号;扩频测量/非扩频数传等时复用信道体制满足以下性质:(1)扩频测量段与非扩频数传段同步复用物理信道,以相同的基带形式调制在同一载波上,扩频测量段与非扩频数传段之间无接缝和间隙;(2)对UTC定义的1PPS间隔时间段内,即UTC定义的1s间隔,物理信道被等分为F个等时复用段,周期为T=1/,第1段的段边沿与1PPS严格对齐,F决定了测量信息输出率;(3)在每个复用段内,定义两个虚拟信道:扩频信道和非扩频数传信道,这两个信道依次占用物理信道的各个等时复用段,每个等时复用段都具有相同的结构和占空比,扩频信道占用等时复用段周期的r%,则非扩频数传信道占用等时复用段周期的1-r%,在每个复用段内非扩频数传信道周期为扩频信道周期的整数倍数关系;(4)扩频信道和非扩频数传信道的传输时钟同源,来自机群链路通信终端的频标,非扩频数传信道数传码速率为扩频信道码片速率的整数倍数关系,比例为P;(5)在每个复用段内扩频信道内同步加载M个传输帧,在每个传输帧内同步加载N个数据位,每个数据位周期为扩频码周期,扩频码长为L位;(6)本体制在QPSK调制体制实现有两种方式:①I/Q两支路同时实施;②对I支路实施,Q支路为非扩频数传通道;(7)本体制启动时发送方采用纯扩频体制引导,待接收方跟踪锁定后对发送方传送遥控指令或者数据注入控制模式切换,发送方确认后转入扩频测量/非扩频数传等时复用体制;当接收方跟踪失锁时,对发送方传送遥控指令或者数据注入控制模式切换,发送方确认后转入纯扩频体制;(8)扩频测量/非扩频数传等时复用体制和纯扩频体制的整个模式确认、切换过程为双向反馈确认和主动控制过程,接收方根据跟踪锁定条件对发送方通过遥控指令或者数据注入进行模式控制,并且发送方在扩频信道和非扩频数传信道的传输帧内事先定义引导标记供接收方提取用于识别同步复用的边界;在机群链路扩频测量/非扩频数传等时复用体制中,机群链路终端采用伪码解扩/解调+非扩频数传解调组合环节,利用当前信息相互引导获得更高的性能;另外,机群编队成员节点之间的相对运动和节点本地频标漂移引起的视在多普勒频移将会带来下一段扩频信道捕获跟踪的困难,发生跟踪失锁甚至失捕;因此,各阶多普勒参数的高精度测量和计算对准确预报下一段扩频信道的多普勒频移非常重要,并直接决定是否能够准确保持跟踪锁定状态或者失捕时快速重捕:对于每秒复用段数F=5,扩频信道占空比r%=25%,扩频码长L=1023,扩频码速率=5.115Mcps,非扩频数传码速率:扩频码片速率倍率P=20,每个复用段的扩频信道内传输帧数M=2,则非扩频数传码速率=102.3Mbps,传输帧长N=125bit,每个复用段内扩频信道持续周期为T=F-1×r%=50ms,对于一般接收单元的捕获跟踪时间已经足够,若存在前几次的多普勒频移参数测量估计值,则能获得更快的捕获时间;本质上提出的扩频测量/非扩频数传等时复用体制对扩频信号的捕获跟踪类似于等时突发扩频通信模式,因此充分利用先验信息辅助捕获和跟踪能保证捕获概率、失锁概率和接收性能满足要求,使其与常规的连续扩频体制相同;设计原则如下:①扩频测量信道的占空比r%由信息数据传输所需码速率和最高传输码率的比值决定,结合秒内分段数F,确定每段内扩频信道持续周期T=F-1×r%,此周期必须满足接收单元完成捕获并转入粗跟踪状态——跟踪残余多普勒在±250Hz内、扩频码片对齐到1/4码片精度;②扩频测量信道的占空比r%受通信数传任务需求约束,秒内分段数F由测量输出率决定,这两者都有一定选择空间,需根据测量双方相对运动动态程度确定,F过小导致测量输出率不足,不满足测量要求,过大了将引起每段内扩频信道持续周期缩短,不利于接续捕获跟踪,需折中考虑;③测量输出率较低时每段内扩频信道持续周期T较长,有利于频域并行直接捕获,测量输出率较高时每秒内测量输出点数F较多,有利于跟踪保持;(二)、机群链路通信/测量综合信道体制的的信号跟踪模块对机群编队内成员飞机之间相对运动,引入统计估计模型,采用统计滤波算法,利用当前多普勒频移测量值实时估计载波参数,包括:载波相位、载波多普勒频移、多普勒频移的1阶/2阶导数、扩频码相位,码时钟则利用载波辅助码环跟踪算法实时计算;扩频信道结束并进入非扩频数传信道后,这些参数能够用于实时递推预报载波参数和码相位参数,接收方可利用这些参数实时计算载波跟踪环路开环控制所需的载波NCO频率控制字、码NCO时钟频率控制字、码相位实时估计值,控制载波跟踪环和码跟踪环NCO历元计数器,在非扩频数传信道周期内引导开环跟踪并预报进入下一段扩频测量信道的前边界时刻及其载波多普勒频率;给出以下三个参数预报公式:①载波多普勒频移的实时预报值计算公式f d ( t k + 1 ) = f d ( t k ) + f · d ( t k ) Δt + 1 2 f · · d ( t k ) Δ t 2 + 1 6 f · · · d ( ϵ k ) Δ t 3 + μ f Δt - - - ( 1 ) ]]>( k = 0,1,2 , . . . , N , Δt = t k - t k + 1 , ϵ k ∈ [ t k , t k + 1 ] ) ]]> 公式(1)中:fd(tk)、
分别为载波多普勒频移在tk时刻的0阶、1阶、2阶导数;公式(1)式右侧前4项为fd(t)(t∈[tk,tk+1])3阶连续可导成分的泰勒展开式(
为余项),第5项μfΔt为fd(t)(t∈[tk,tk+1])3阶连续可导成分之外的剩余成分,这两项为未知项,成为预报计算公式的模型误差;令:δfd(tk)、
分别为fd(tk)、
的估计误差,则tk+1时刻的载波多普勒频移的预报误差δfd(tk+1|tk)为:δ f d ( t k + 1 | t k ) = δ f d ( t k ) + δ f · d ( t k ) Δt + 1 2 δ f · · d ( t k ) Δ t 2 + 1 6 f · · · d ( ϵ k ) Δ t 3 + μΔt - - - ( 2 ) ]]>( k = 0,1,2 , . . . , N , Δt = t k - t k + 1 , , ϵ k ∈ [ t k , t k + 1 ] ) ]]> ②码时钟频率的实时预报值计算公式由于码相位跟踪过程依赖码时钟控制码NCO,采用载波辅助码环跟踪策略通过载波NCO频率控制字来计算码NCO频率控制字,设:fcode(t)为码时钟频率(fcode为标称值),fRF(t)载波频率(fRF为标称值),令:λ为码时钟(码片速率)与射频输出频点的比率,则:λ = f code f RF = f code ( t ) f RF ( t ) - - - ( 3 ) ]]> fRF(t)=fRF+fd(t) (4)根据公式(3)、公式(4)得:fcode(t)=λfRF(t) =fcode+λfd(t) (5)fcode(n)(t)=λfRF(n)(t) =λfd(n)(t) (n=1,2,3,...,) (6)公式(6)中:fcode(n)(t)为码时钟频率在t时刻的n阶导数(n=1,2,3,...);根据公式(1)、公式(5)得:f code ( t k + 1 ) = f code + λ f d ( t k + 1 ) ]]>= f code + λ f d ( t k ) + λ f · d ( t k ) Δt + 1 2 λ f · · d ( t k ) Δ t 2 + 1 6 λ f · · · d ( ϵ k ) Δ t 3 + λ μ f Δt - - - ( 7 ) ]]>= f code ( t k ) + f · code ( t k ) Δt + 1 2 f · · code ( t k ) Δ t 2 + 1 6 λ f · · · d ( ϵ k ) Δ t 3 + λ μ f Δf ]]>( k = 0,1,2 , . . . , N , Δt = t k - t k + 1 ) ]]> 根据公式(1)~公式(7),得到tk+1时刻码时钟频率的预报误差δfcode(tk+1|tk)为:δ f code ( t k + 1 | t k ) = λδ f d ( t k ) + λδ f · d ( t k ) Δt + 1 2 λδ f · · d ( t k ) Δ t 2 + 1 6 λ f · · · d ( ϵ k ) Δ t 3 + λ μ f Δt - - - ( 8 ) ]]>( k = 0,1,2 , . . . , N , Δt = t k - t k + 1 , ϵ k ∈ [ t k , t k + 1 ] ) ]]> 在非扩频数传信道周期内,利用公式(1)、公式(2)实时计算载波多普勒参数、码多普勒参数的预报值并生成频率控制字用于控制载波NCO、码NCO在非扩频数传信道开环工作,这种模式基于参数预报实施开环跟踪;载波NCO和码NCO计数器历元计数输出值达到设定值——由非扩频数传信道周期决定,则可判断到达了下一段扩频测量信道的前沿;在扩频测量信道周期内,以前一段非扩频数传信道周期的计数中止时刻外推预报值作为载波捕获/跟踪环路的初值转入对扩频信号的闭环跟踪,基于载波跟踪环路自动频率跟踪+锁相环与码跟踪环路数字延迟锁定环闭环控制载波NCO、码NCO工作,这种模式基于跟踪残差反馈控制实施闭环跟踪;③码相位的实时预报值计算公式给出机群编队内成员节点之间相对速度各阶导数与多普勒频移各阶导数的关系公式:v ( n ) ( t ) = c f d ( n ) ( t ) f RF - - - ( 9 ) ]]> 由公式(9),考虑运动学理论,得到:d ( t k + 1 ) = d ( t k ) + v ( t k ) Δt + 1 2 v · ( t k ) Δ t 2 + 1 6 v · · ( t k ) Δ t 3 + 1 12 v · · · ( ϵ k ) Δ t 4 + μ v Δt - - - ( 10 ) ]]>( k = 0,1,2 , . . . , N , Δt = t k - t k + 1 , ϵ k ∈ [ t k , t k + 1 ] ) ]]> 公式(10)中:d(tk)为tk时刻两机载天线相位中心之间的几何距离;公式(10)式右侧前4项为v(t)(t∈[tk,tk+1])3阶连续可导成分的泰勒展开式(
为余项),第5项μvΔt为v(t)(t∈[tk,tk+1])3阶连续可导成分之外的剩余成分,这两项为未知项,成为预报计算公式的模型误差;由公式(10)得:τ ( t k + 1 ) = τ ( t k ) + f d ( t k ) f RF Δt + 1 2 f · d ( t k ) f RF Δ t 2 + 1 6 f · · d ( t k ) f RF Δ t 3 + 1 12 f · · · d ( ϵ k ) f RF Δ t 4 + μ v c Δt ]]>( k = 0,1,2 , . . . , N , Δt = t k - t k + 1 , ϵ k ∈ [ t k , t k + 1 ] ) - - - ( 11 ) ]]> 公式(11)中:τ(tk)为tk时刻电波在机群编队内两成员节点之间的空间传播时延,τ(tk)=d(tk)/c(c为真空光速);令:δτ(tk)为τ(tk)的估计误差,则tk+1时刻的预报误差δτ(tk+1|tk)为:δτ ( t k + 1 | t k ) = δτ ( t k ) + δ f d ( t k ) f RF Δt + 1 2 δ f · d ( t k ) f RF Δ t 2 + 1 6 δ f · · d ( t k ) f RF Δ t 3 ]]>+ 1 12 f · · · d ( ϵ k ) f RF Δ t 4 + μ v c Δt , ( k = 0,1,2 , . . . , N , Δt = t k - t k + 1 , ϵ k ∈ [ t k , t k + 1 ] ) - - - ( 12 ) ]]> 定义码相位对齐误差为φcode(t),有:φ code ( t ) = δτ ( t ) T code = δτ ( t ) · f code - - - ( 13 ) ]]> 由公式(12)、公式(13)得tk+1时刻的码相位预报误差为:φ code ( t k + 1 | t k ) = [ δτ ( t k ) + δ f d ( t k ) f RF Δt + 1 2 δ f · d ( t k ) f RF Δ t 2 + 1 6 δ f · · d ( t k ) f RF Δ t 3 ]]>+ 1 12 f · · · d ( ϵ k ) f RF Δ t 4 + μ v c Δt ] · f code - - - ( 14 ) ]]>= φ code ( t k ) + λδ f d ( t k ) Δt + 1 2 λδ f · d ( t k ) Δ t 2 + 1 6 λδ f · · d ( t k ) Δ t 3 ]]>+ 1 12 λ f · · · d ( ϵ k ) Δ t 4 + μ v c T code Δt , ( k = 0,1 , . . . , N , Δt = t k - t k + 1 , ϵ k ∈ [ t k , t k + 1 ] ) ]]> 公式(14)中:φcode(tk)为tk时刻的码相位对齐误差;当前段非扩频数传信道结束时,码相位预报误差φcode(tk+1|tk)将成为下一段扩频测量信道的初始码相位对齐误差;因此,载波频率初始跟踪误差、码时钟频率初始跟踪误差、码相位初始跟踪误差由多普勒频移外推预报误差引起,因素有:①当前载波参数、码相位参数测量误差;②机群编队内成员节点之间相对运动状态在非扩频数传信道周期内发生变化,引起多普勒频率三阶导数增大,加速度发生突变或者不可导;③统计估计模型的准确性和滤波算法的性能引起多普勒频移各阶导数估计误差;公式(2)描述了进入下一段扩频测量信道时的载波频率初始跟踪误差,公式(8)描述了进入下一段扩频测量信道时的码时钟频率初始跟踪误差,公式(14)描述了进入下一段扩频测量信道时的扩频码相位初始跟踪误差;简单起见,仅考虑fd(t)三阶可导成分,忽略泰勒展开式余项;令:sup t | φ code ( t ) | = μ 0 - - - ( 16 ) ]]>sup t | δ f d ( t ) | = μ 1 - - - ( 17 ) ]]>sup t | δ f · d ( t ) | = μ 2 - - - ( 18 ) ]]>sup t | δ f · · d ( t ) | = μ 3 - - - ( 19 ) ]]> 在每段扩频测量信道末端,扩频码相位跟踪误差约为码片宽度的1%,载波跟踪误差小于1Hz;取较为恶劣的估计误差条件,令公式(16)~(19)中:μ0=0.01,μ1=1Hz,μ2=50Hz/s,μ3=500Hz/s2;若:码时钟fcode=5.115MHz,fRF=2250.6MHz(S频段),则λ=1/440,Ts=0.2s,利用公式(2)、公式(14)计算得:δfd(tk+1|tk)<21Hz、φcode(tk+1|tk)<0.0142,载波多普勒预报误差主要来源于载波多普勒频移一阶、二阶导数估计误差,而码相位预报误差主要来源于当前码相位估计误差,载波跟踪环和码跟踪环在比较恶劣的估计误差条件下仍然能够获得比较小的初始跟踪误差,几乎是保持锁定状态;可以认为,载波参数估计误差对载波多普勒频移预报精度影响远大于码相位预报精度;因此,在扩频测量/非扩频数传等时复用体制中,为保持对载波频率和码相位准确可靠地开环跟踪,必须获得准确的多普勒频移及其一阶、二阶导数估计值;在扩频信道内前边界配置引导信息,如:信息位重复多次“1”或者“0”以引导和确保捕获跟踪成功,并对扩频传输帧提供前方保护;在非扩频数传信道内前边界配置多次重复的同步码以引导和确保解调数据同步跟踪成功;若接收方确认已失捕失锁,则利用前向链路纯扩频信道向发送方发送指令使返向链路转入指定的纯扩频通信模式,直到重捕成功并跟踪锁定后再指令控制发送方返向信道转入扩频测量/非扩频数传等时复用体制。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国科学院空间科学与应用研究中心,未经中国科学院空间科学与应用研究中心许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201010103897.2/,转载请声明来源钻瓜专利网。