CN117614786A - 基于卡尔曼滤波的载波和伪码的同步方法 - Google Patents
基于卡尔曼滤波的载波和伪码的同步方法 Download PDFInfo
- Publication number
- CN117614786A CN117614786A CN202311562977.8A CN202311562977A CN117614786A CN 117614786 A CN117614786 A CN 117614786A CN 202311562977 A CN202311562977 A CN 202311562977A CN 117614786 A CN117614786 A CN 117614786A
- Authority
- CN
- China
- Prior art keywords
- pseudo code
- carrier
- phase
- kalman filtering
- pseudo
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000001914 filtration Methods 0.000 title claims abstract description 101
- 238000000034 method Methods 0.000 title claims abstract description 47
- 239000011159 matrix material Substances 0.000 claims abstract description 63
- 230000004927 fusion Effects 0.000 claims abstract description 37
- 238000012545 processing Methods 0.000 claims description 30
- 230000008859 change Effects 0.000 claims description 24
- 230000000875 corresponding effect Effects 0.000 claims description 17
- 238000009825 accumulation Methods 0.000 claims description 14
- 238000013507 mapping Methods 0.000 claims description 10
- 238000004422 calculation algorithm Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000003111 delayed effect Effects 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 8
- 230000002596 correlated effect Effects 0.000 claims description 6
- 230000001360 synchronised effect Effects 0.000 claims description 6
- 230000007704 transition Effects 0.000 claims description 6
- 230000010354 integration Effects 0.000 claims description 5
- 238000005311 autocorrelation function Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 230000010355 oscillation Effects 0.000 claims description 3
- 230000010363 phase shift Effects 0.000 claims description 3
- 230000007480 spreading Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000004891 communication Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000008094 contradictory effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000010587 phase diagram Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 239000012050 conventional carrier Substances 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000026676 system process Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/0014—Carrier regulation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/0014—Carrier regulation
- H04L2027/0044—Control loops for carrier regulation
- H04L2027/0063—Elements of loops
- H04L2027/0069—Loop filters
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
本发明公开了一种基于卡尔曼滤波的载波和伪码的同步方法,该方法首先将载波鉴相输出进行卡尔曼滤波,一次迭代之后更新输出状态变量与协方差矩阵,与此同时将伪码也进行超前滞后伪码鉴相,鉴相输出结果为码偏移,将鉴相结果转化为码相位差,并进行卡尔曼滤波,一次迭代之后也更新输出状态变量与协方差矩阵;然后将伪码输出状态变量转换为载波多普勒偏移,对两组状态输出进行线性最小方差融合估计,求出新的状态变量值以及对应的协方差矩阵;最后将新的状态变量值与协方差矩阵反馈回卡尔曼滤波的步骤中,进行新一轮运算。本发明基于卡尔曼滤波,结合伪码与载波同步结果,提升了载波同步精度,提高了接收机性能。
Description
技术领域
本发明属于数字通信同步领域,具体为一种基于卡尔曼滤波的载波和伪码的同步方法。
背景技术
同步技术广泛应用于数字通信技术之中载波跟踪的目的是对接收本振精确同步,使得载波与本地本振之间的残余相位误差最小,传统的载波跟踪技术是在70年代建立的。传统的锁相环由三个基本组成模块:一个相位检测器,负责提供输出相位误差测量值;一个环路滤波器,平滑由相位检测输出的热噪声引起的变异值;最后是压控振荡器,用于根据环路滤波器输出产生本地载波校正接收信号。
但是这种传统的载波跟踪技术无法满足现阶段的要求,主要原因是这种结构在消除噪声与跟踪高动态方面的能力是相互矛盾的。在传统架构中,这两种需求是高度耦合的,因此在不引起一种需求急剧下降的情况下,很难为这两个方面提供更好的方案。例如,专注于高动态的设计容易遭受严重的抖动,这将引入严重的噪声到估计的载波相位。另一方面,专注于噪声抑制的设计能够使输出相位的抖动最小化,但要以承受动态应力误差为代价(即当跟踪环路无法跟随输入的高阶动态)。在这些情况下,最终的结果可能是跟踪环路丢失。
发明内容
本发明的目的在于提供了一种基于卡尔曼滤波的载波和伪码的同步方法,以解决传统环路的环路带宽与跟踪能力矛盾的问题的同时,提高在低信噪比情况下载波跟踪能力与收敛速度。
实现本发明目的的技术方案为:一种基于卡尔曼滤波的载波和伪码的同步方法,具体步骤为:
步骤1:对接收机收到的中频信号进行数字正交下变频,得到两路正交IQ信号;
步骤2:将两路正交IQ信号分别与超前支路伪码、滞后支路伪码以及中间支路伪码进行解扩处理得到六路解扩信号,对超前、滞后四路解扩信号进行伪码鉴相处理,对中间支路两路解扩信号进行载波鉴相处理;
步骤3:将伪码鉴相处理得到的信号作为伪码卡尔曼滤波的输入,将作为伪码卡尔曼滤波器的状态变量进行卡尔曼滤波,其中,/>为伪码相位差,ωτe为伪码频率差,/>为伪码频率变化率;
将载波鉴相处理得到的信号作为载波卡尔曼滤波的输入,将作为载波卡尔曼滤波器状态变量进行卡尔曼滤波,其中,/>为载波相位差,ωe为载波频率差,/>载波频率变化率;
步骤4:将伪码卡尔曼滤波输出与载波卡尔曼滤波输出分别反馈回本地伪码发生器与数字NCO进行补偿处理,获得同步结果;
步骤5:将步骤3中的两个卡尔曼滤波的状态变量值与协方差矩阵进行线性最小方差融合估计,将融合估计后的状态变量值与融合协方差矩阵作为新的卡尔曼滤波器的状态变量值与协方差矩阵,返回步骤3,对信号持续同步。
优选地,两路正交IQ信号分别与超前伪码、滞后伪码进行解扩处理得到的四路超前、滞后四路解扩信号具体为:
式中,IE为I路与超前支路伪码相关累加后的值,IL为I路与滞后支路伪码相关累加后的值,QE为Q路与超前支路伪码相关累加后的值,QL为Q路与滞后支路伪码相关累加后的值,A为信号幅值,Dk表示第k个符号数据,ε为本地中间支路参考伪码与接收信号的伪码相差值,d为本地生成的各路参考伪码之间的相差,sinc(z)=sin(z)/z,z为变量,fd为多普勒频偏估计残差,Tb为预计相关积分时间,k表示第k个符号,表示载波初始相差值,R(ε)为伪码的自相关函数。
优选地,伪码鉴相采用的鉴相算法为超前滞后能量鉴相算法:
D(ε)=(IL(k)2+QL(k)2)-(IE(k)2+QE(k)2)=A2(R2(ε+d)-R2(ε-d))
式中,D(ε)为伪码鉴相结果,IE为I路与超前支路伪码相关累加后的值,IL为I路与滞后支路伪码相关累加后的值,QE为Q路与超前支路伪码相关累加后的值,QL为Q路与滞后支路伪码相关累加后的值,A为信号幅值,ε为本地中间支路参考伪码与接收信号的伪码相差值,d为本地生成的各路参考伪码之间的相差。
优选地,对中间支路两路解扩信号进行载波鉴相处理的具体公式为:
其中,sign表示符号函数,I(k)与Q(k)分别为同相和正交支路与中间支路伪码进行解扩输出的解扩信号,err为鉴相输出表示本地载波与接收载波信号的相位差,fd为多普勒频偏估计残差,Tb为预计相关积分时间,k表示第k个符号,表示载波初始相差值。
优选地,将载波鉴相处理得到的信号作为载波卡尔曼滤波的输入,将作为载波卡尔曼滤波器状态变量进行卡尔曼滤波的具体方法为:
根据载波多普勒频率增量确定Δt时间内载波多普勒频率产生的载波相位增量,具体为:
确定载波相位的变化值,即Δt时间内本地振荡与接收信号相位的差值与多普勒频率在Δt时间内的增量之和,具体为:
其中/>其中,/>为载波多普勒相位差,ωe为载波多普勒频率差,/>为载波多普勒频率变化率;
将载波相位状态转移方程表示为:
其中,Wn为均值为0的高斯白噪声;
中间支路两路解扩信号与正交支路经过载波鉴相器,在鉴相更新周期内得到相位差的平均值:
其中,为相位差的平均值,Ts为载波鉴相更新周期;
确定观测矩阵中的各个元素,观测矩阵为Hf:
确定载波卡尔曼滤波状态矩阵φf与观测矩阵Hf分别为:
将计算得到的状态矩阵φf与观测矩阵Hf代入载波卡尔曼滤波的动态模型进行滤波处理,其中,载波卡尔曼滤波的动态模型具体为:
Xf(k+1)=φfXf(k)+Wf(k)
Yf(k)=HfXf(k)+Vf(k)
其中
优选地,将伪码鉴相处理得到的信号作为伪码卡尔曼滤波的输入,将作为伪码卡尔曼滤波器的状态变量进行卡尔曼滤波的具体方法为:
确定Δt时间内伪码的多普勒频率产生的伪码相位增量为:
Δt时间内产生的伪码相位变化值为本地伪码与接收伪码的相位的差值加上伪码多普勒频率在该时间内的增量:
其中/>
式中,ωτe为伪码频率差,为伪码频率变化率;
根据Δt时间内产生的伪码相位变化值,确定伪码相位状态转移方程为:
其中Wn为均值为0的高斯白噪声;
解扩后超前支路和同相支路经过伪码鉴相器,其输出为鉴相更新周期内得到的码相位差的平均值:
式中,Ts为更新周期;
将伪码的状态矩阵φτ与观测矩阵Hτ代入卡尔曼滤波公式,获得滤波结果,其中,伪码的状态矩阵φτ与观测矩阵Hτ分别为:
。
优选地,将步骤3中的两个卡尔曼滤波的状态变量值与协方差矩阵进行线性最小方差融合估计的具体过程为:
步骤5-1、将伪码卡尔曼滤波输出的伪码相位偏移映射为载波对应的多普勒相位偏移,具体过程如下:
伪码卡尔曼滤波输出为根据鉴相算法/>其中A为信号幅度增益,计算出实际的伪码偏移大小为/>
确定伪码频率偏移映射出的载波的多普勒频率偏移为:
其中,fc表示载波频率,Rc为PN码的码元速率,fτd表示伪码多普勒偏移频率,Rc1为实际的伪码速率,Nc为扩频码的长度,ε为实际的伪码偏移,Tc为伪码的码元周期;
根据载波的鉴相输出为sin(2πfdkTb),计算出由伪码映射出的载波多普勒相位偏移近似于2πfd1kTb,即将伪码卡尔曼滤波输出映射为载波卡尔曼滤波输出:
其中,为伪码卡尔曼滤波状态变量映射为载波卡尔曼滤波对应的状态变量,其协方差矩阵为/>
步骤5-2、将两个卡尔曼滤波输出结果进行线性最小方差融合估计,具体计算步骤如下:
线性最小方差融合估计:
对一个任意随机向量X∈Rn,假设有L个无偏独立估计对应的协方差矩阵为/>则/>就是X的线性融合估计,其中/>融合估计的协方差矩阵:
将伪码卡尔曼滤波输出映射为载波卡尔曼滤波的状态变量的值与载波卡尔曼滤波的输出Xf作为线性最小方差融合估计的输入变量,/>与Pf为对应的协方差矩阵;
融合估计步骤如下:
其中为线性最小方差融合估计值,/>为对应的协方差矩阵。将融合估计输出作为载波卡尔曼滤波预测步的输入,并通过步骤5-1的逆变换计算出伪码卡尔曼滤波预测步的状态变量更新值与协方差矩阵更新值。
本发明与现有技术相比其显著优点为:(1)本发明弥补了传统环路滤波器环路带宽与跟踪精度相矛盾的问题;(2)本发明通过对两个卡尔曼滤波输出结果进行线性最小方差融合估计,提高了同步精度,降低了其协方差,计算结果更准确;(3)在低信噪比情况下,能提高收敛速度,减少锁定时间。(4)本发明计算方法,计算复杂度适中,能满足信号同步跟踪的要求。
下面结合附图对本发明作进一步的详细描述
附图说明
图1是本发明的流程图。
图2单独使用卡尔曼滤波实现载波同步的鉴相图。
图3使用线性最小方差融合估计实现载波同步的鉴相图。
图4为卡尔曼滤波具体步骤流程图。
具体实施方式
一种基于卡尔曼滤波的载波和伪码的同步方法,包括以下步骤:
步骤1:首先对接收机收到的中频信号进行数字正交下变频,得到两路正交信号I,Q;
步骤2:将两路正交信号分别与超前伪码,滞后伪码以及中间支路伪码进行解扩处理得到六路解扩信号,对超前滞后四路解扩信号进行伪码鉴相处理,并对中间支路解扩输出两路信号进行载波鉴相处理;
步骤3:将伪码鉴相输出作为伪码卡尔曼滤波的输入,载波鉴相输出作为载波卡尔曼滤波的输入,将作为载波卡尔曼滤波器状态变量,其中/>为载波相位差,ωe为载波频率差,/>载波频率变化率,进行卡尔曼滤波,同时将/>作为伪码卡尔曼滤波器的状态变量进行卡尔曼滤波,其中/>为伪码相位差,ωτe为伪码频率差,/>为伪码频率变化率;
步骤4:将伪码卡尔曼滤波输出与载波卡尔曼滤波输出分别反馈回本地伪码发生器与数字NCO进行补偿处理,获得同步结果;
步骤5:将步骤3中的两个卡尔曼滤波的状态变量值与协方差矩阵进行线性最小方差融合估计,将融合估计后的状态变量值与融合协方差矩阵作为新的卡尔曼滤波器的状态变量值与协方差矩阵,返回步骤3,对信号持续同步。
进一步的实施例中,步骤2中分别对伪码与载波进行鉴相处理,具体步骤为:
步骤2-1、计算伪码鉴相的输出,具体为:
输入两路正交信号,将其分别与超前支路伪码与滞后支路伪码进行相关运算,计算出超前滞后四路解扩信号,其表达式为:
其中:
IE:I路与超前支路伪码相关累加后的值
IL:I路与滞后支路伪码相关累加后的值
QE:Q路与超前支路伪码相关累加后的值
QL:Q路与滞后支路伪码相关累加后的值
A:信号幅值
Dk:表示第k个符号数据
ε:本地中间支路参考伪码与接收信号的伪码相差(ε≤0.5)
d:本地生成的各路参考伪码之间的相差(本系统取d=0.5,即半个码长,则超前和滞后两个支路的伪码相差为一个码片宽度)
sinc(z)=sin(z)/z
fd:多普勒频偏估计残差
Tb:预计相关积分时间
k:表示第k个符号
表示载波初始相差
R(ε):伪码的自相关函数
这里认为sinc(fdTb)≈1,忽略噪声带来的影响。
伪码鉴相采用的鉴相算法为超前滞后能量鉴相算法:
D(ε)=(IL(k)2+QL(k)2)-(IE(k)2+QE(k)2)=A2(R2(ε+d)-R2(ε-d))
利用伪码序列的自相关函数计算出:
Dg(ε)表示归一化后伪码鉴相输出结果
步骤2-2、计算载波鉴相输出,具体为:
输入两路正交信号,将其与中间支路伪码进行解扩处理,计算出中间支路的解扩信号,采用硬判决鉴相算法:
其中sign表示符号函数,I(k)与Q(k)分别为同相和正交支路与中间支路伪码进行解扩输出的结果。err为鉴相输出表示本地载波与接收载波信号的相位差。
进一步的实施例中,步骤3中分别对载波与伪码鉴相输出进行卡尔曼滤波,具体步骤为:
步骤3-1、载波卡尔曼滤波状态矩阵与观测矩阵计算过程如下:
载波多普勒频率增量表达式为fd(t+τ)=fd(t)+adτ
其中fd(t)为t时刻载波多普勒频率值,ad为载波多普勒频率变化率,τ为时间间隔
因此在Δt时间内载波多普勒频率产生的载波相位增量为:
在Δt时间内产生的载波相位变化值为本地振荡与接收信号相位的差值加上多普勒频率在该时间内的增量
其中/>
ωe为载波多普勒频率差,为载波多普勒频率变化率,故将载波相位状态转移方程表示为:
其中Wn为均值为0的高斯白噪声Wn=[w0,w1,w2],其协方差为Q
其中Q0,Qf是由接收机基准时钟相位和频率引起的随机误差,Qa是由动态模型决定的。
本发明认为接收机基准时钟引起的随机误差为0,认为噪声主要是由动态模型决定的。
观测值为解扩后同相支路与正交支路经过载波鉴相器的输出,鉴相输出为鉴相更新周期内得到的相位差的平均值:
其中Ts为载波鉴相更新周期,为一个符号周期
因此载波卡尔曼滤波状态矩阵φf与观测矩阵Hf分别为:
步骤3-2、伪码卡尔曼滤波状态矩阵与观测矩阵计算过程如下:
伪码多普勒偏移频率为:
Rc:伪码速率
fc:载波频率
fd:载波多普勒频偏
伪码多普勒频率增量表达式为fτd(t+σ)=fτd(t)+aτdτ
上式中fτd(t)为t时刻伪码多普勒频率值,aτd为伪码多普勒频率变化率,τ为时间间隔,因此在Δt时间内伪码的多普勒频率产生的伪码相位增量为:
在Δt时间内产生的伪码相位变化值为本地伪码与接收伪码的相位的差值加上伪码多普勒频率在该时间内的增量
其中/>
ωτe为伪码频率差,为伪码频率变化率,故将伪码相位状态转移方程为:
其中Wn为均值为0的高斯白噪声Wn=[w0,w1,w2]
观测值为解扩后超前支路和同相支路经过鉴相器的输出,鉴相输出为鉴相更新周期(这里更新周期为Ts)内得到的码相位差的平均值:
因此伪码的状态矩阵φτ与观测矩阵Hτ分别为:
步骤3-3、对伪码鉴相输出与载波鉴相输出同时进行卡尔曼滤波
卡尔曼滤波是一种时域设计方法,他把状态空间的概念引入随机估计理论中,把信号的处理过程看成是白噪声作用下的线性系统的输出,估计过程利用状态方程,观测方程以及噪声误差(系统过程噪声与观测噪声),对这种输入输出关系进行表述。
使用如下的状态方程与观测方程描述该系统:
X(k+1)=φX(k)+W(k)
Y(k)=HX(k)+V(k)
其中k为离散时间变量,X(k)表示系统在时刻k的状态,Y(k)为对应状态的观测信号,φ为状态转移矩阵,H为观测矩阵,W(k)为系统噪声,V(k)为观测噪声。
下面为具体的卡尔曼滤波估计算法的计算过程:
(1)状态一步预测:
(2)状态更新:
(3)滤波增益值计算:K(k+1)=P(k+1|k)HT[HP(k+1|k)HT+R]-1
(4)一步预测协方差阵:P(k+1|k)=φP(k|k)φT+Q
(5)协方差阵更新:P(k+1|k+1)=[I-K(k+1|k)HP(k+1|k)]
上式中K为卡尔曼增益,I为单位阵,R为观测噪声的协方差矩阵,Q为系统噪声对应的协方差矩阵。
载波卡尔曼滤波的动态模型为:
Xf(k+1)=φfXf(k)+Wf(k)
Yf(k)=HfXf(k)+Vf(k)
其中预测方程组为:
Pf(k+1|k)=φfPf(k|k)φf T+Qf
观测方程组为:
Kf(k+1)=Pf(k+1|k)Hf T[HfPf(k+1|k)Hf T+Rf]-1
Pf(k+1|k+1)=[I-Kf(k+1|k)HfPf(k+1|k)]
在载波卡尔曼滤波中表示上一时刻(k时刻)的最佳估计结果,表示当前时刻的预测值,Pf(k|k)表示上一时刻估计协方差,Pf(k+1|k)表示当前时刻的协方差预测值,Kf(k+1)表示当前时刻的卡尔曼滤波增益,/>表示当前时刻最佳估计结果,Pf(k+1|k+1)表示当前时刻协方差最佳估计结果,Qf表示系统噪声协方差矩阵,Rf表示观测噪声协方差矩阵,/>表示当前时刻状态变量更新值,Yf(k+1)表示当前时刻观测变量观测值,φf表示状态矩阵,Hf表示观测矩阵。
伪码卡尔曼滤波的动态模型为:
其中预测方程组为:
观测方程组为:
/>
伪码卡尔曼滤波具体参数不再赘述,与以上载波卡尔曼滤波一致。
进一步的实施例中,步骤5中将载波与伪码卡尔曼滤波输出的状态变量与协方差矩阵进行融合的过程,具体步骤为:
步骤5-1、将伪码卡尔曼滤波输出的伪码相位偏移映射为载波对应的多普勒相位偏移;伪码频率映射为载波多普勒频率偏移,这里伪码频率变化率与载波频率变化率均为0,具体计算过程如下:
伪码卡尔曼滤波输出为根据鉴相算法/>其中A为信号幅度增益,即可计算出实际的伪码偏移大小为/>根据步骤3-2,伪码频率偏移映射出的载波的多普勒频率偏移为:
其中,fc表示载波频率,Rc为PN码的码元速率,fτd表示伪码多普勒偏移频率,Rc1为实际的伪码速率,Nc为扩频码的长度,ε为实际的伪码偏移,Tc为伪码的码元周期
载波的鉴相输出为sin(2πfdkTb),可以计算出由伪码映射出的载波多普勒相位偏移近似于2πfd1kTb,即可将伪码卡尔曼滤波输出映射为载波卡尔曼滤波输出:
其中为伪码卡尔曼滤波状态变量映射为载波卡尔曼滤波对应的状态变量。步骤5-2、将两个卡尔曼滤波输出结果进行线性最小方差融合估计,具体计算步骤如下:
线性最小方差融合估计:
对一个任意随机向量X∈Rn,假设有L个无偏独立估计对应的协方差矩阵为/>则/>就是X的线性融合估计,其中/>融合估计的协方差矩阵:
将伪码卡尔曼滤波输出映射为载波卡尔曼滤波的状态变量的值与载波卡尔曼滤波的输出Xf作为线性最小方差融合估计的输入变量,/>与Pf为对应的协方差矩阵。
融合估计步骤如下:
其中为线性最小方差融合估计值,/>为对应的协方差矩阵。将融合估计输出作为载波卡尔曼滤波预测步的输入,并通过步骤5-1的逆变换计算出伪码卡尔曼滤波预测步的状态变量更新值与协方差矩阵更新值,进行下一步卡尔曼滤波
下面结合附图与实例1对本发明作详细说明
实施例1
如图1所示,基于卡尔曼滤波的载波和伪码的同步方法,步骤如下:
仿真条件:本系统是基于扩频通信系统进行仿真,载波频率fc=150MHz,时钟速率为200MHz,PN码速率为Rc=40MHz,扩频码长为Nc=80,多普勒频率设置为8k,卡尔曼滤波迭代更新时间与鉴相算法更新时间保持一致为500KHz为一个符号的更新时间,信噪比为0dB。
步骤1:首先对接收机收到的中频信号进行数字正交下变频,得到两路正交信号I,Q。;
步骤2:将两路正交信号分别与超前伪码,滞后伪码以及中间支路伪码进行解扩处理得到六路解扩信号,对超前滞后四路解扩信号进行伪码鉴相处理,并对中间支路解扩输出两路信号进行载波鉴相处理;
步骤3:将伪码鉴相输出作为伪码卡尔曼滤波的输入,载波鉴相输出作为载波卡尔曼滤波的输入,将作为载波卡尔曼滤波器状态变量,其中/>为载波相位差,ωe为载波频率差,/>载波频率变化率,进行卡尔曼滤波,同时将/>作为伪码卡尔曼滤波器的状态变量进行卡尔曼滤波,其中/>为伪码相位差,ωτe为伪码频率差,/>为伪码频率变化率;
步骤4:将步骤3中的两个卡尔曼滤波的状态变量值与协方差矩阵进行线性最小方差融合估计,将更新后的状态变量值与融合协方差矩阵作为新的卡尔曼滤波器的状态变量值与协方差矩阵,进行下一步卡尔曼滤波。
步骤5:将伪码卡尔曼滤波输出与载波卡尔曼滤波输出分别反馈回本地伪码发生器与数字NCO进行补偿处理。
前260个符号并未进行同步,260个符号之后才开始进行同步。通过图2与图3对比看出,在信噪比为0dB时,使用线性最小方差融合估计卡尔曼滤波比单独使用卡尔曼滤波进行载波同步收敛速度更快,单独使用卡尔曼滤波,在第350个符号处完成了收敛,在使用最小方差融合估计后,在第314个符号处就完成了收敛,在低信噪比下,使用最小方差融合估计载波环收敛的速度远快于单独使用卡尔曼滤波器。大多数情况下信号的捕获需要捕获帧头,因此收敛速度提高对提高接收机性能有重要作用。牺牲了少量初始伪码精度,加快了环路收敛速率。
Claims (7)
1.一种基于卡尔曼滤波的载波和伪码的同步方法,其特征在于,具体步骤为:
步骤1:对接收机收到的中频信号进行数字正交下变频,得到两路正交IQ信号;
步骤2:将两路正交IQ信号分别与超前支路伪码、滞后支路伪码以及中间支路伪码进行解扩处理得到六路解扩信号,对超前、滞后四路解扩信号进行伪码鉴相处理,对中间支路两路解扩信号进行载波鉴相处理;
步骤3:将伪码鉴相处理得到的信号作为伪码卡尔曼滤波的输入,将作为伪码卡尔曼滤波器的状态变量进行卡尔曼滤波,其中,/>为伪码相位差,ωτe为伪码频率差,/>为伪码频率变化率;
将载波鉴相处理得到的信号作为载波卡尔曼滤波的输入,将作为载波卡尔曼滤波器状态变量进行卡尔曼滤波,其中,/>为载波相位差,ωe为载波频率差,/>载波频率变化率;
步骤4:将伪码卡尔曼滤波输出与载波卡尔曼滤波输出分别反馈回本地伪码发生器与数字NCO进行补偿处理,获得同步结果;
步骤5:将步骤3中的两个卡尔曼滤波的状态变量值与协方差矩阵进行线性最小方差融合估计,将融合估计后的状态变量值与融合协方差矩阵作为新的卡尔曼滤波器的状态变量值与协方差矩阵,返回步骤3,对信号持续同步。
2.根据权利要求1所述的基于卡尔曼滤波的载波和伪码的同步方法,其特征在于,两路正交IQ信号分别与超前伪码、滞后伪码进行解扩处理得到的四路超前、滞后四路解扩信号具体为:
式中,IE为I路与超前支路伪码相关累加后的值,IL为I路与滞后支路伪码相关累加后的值,QE为Q路与超前支路伪码相关累加后的值,QL为Q路与滞后支路伪码相关累加后的值,A为信号幅值,Dk表示第k个符号数据,ε为本地中间支路参考伪码与接收信号的伪码相差值,d为本地生成的各路参考伪码之间的相差,sinc(z)=sin(z)/z,z为变量,fd为多普勒频偏估计残差,Tb为预计相关积分时间,k表示第k个符号,表示载波初始相差值,R(ε)为伪码的自相关函数。
3.根据权利要求2所述的基于卡尔曼滤波的载波和伪码的同步方法,其特征在于,伪码鉴相采用的鉴相算法为超前滞后能量鉴相算法:
D(ω)=(IL(k)2+QL(k)2)-(IE(k)2+QE(k)2)=A2(R2(ε+d)-R2(ε-d))
式中,D(ε)为伪码鉴相结果,IE为I路与超前支路伪码相关累加后的值,IL为I路与滞后支路伪码相关累加后的值,QE为Q路与超前支路伪码相关累加后的值,QL为Q路与滞后支路伪码相关累加后的值,A为信号幅值,ε为本地中间支路参考伪码与接收信号的伪码相差值,d为本地生成的各路参考伪码之间的相差。
4.根据权利要求1所述的基于卡尔曼滤波的载波和伪码的同步方法,其特征在于,对中间支路两路解扩信号进行载波鉴相处理的具体公式为:
其中,sign表示符号函数,I(k)与Q(k)分别为同相和正交支路与中间支路伪码进行解扩输出的解扩信号,err为鉴相输出表示本地载波与接收载波信号的相位差,fd为多普勒频偏估计残差,Tb为预计相关积分时间,k表示第k个符号,表示载波初始相差值。
5.根据权利要求1所述的基于卡尔曼滤波的载波和伪码的同步方法,其特征在于,将载波鉴相处理得到的信号作为载波卡尔曼滤波的输入,将作为载波卡尔曼滤波器状态变量进行卡尔曼滤波的具体方法为:
根据载波多普勒频率增量确定Δt时间内载波多普勒频率产生的载波相位增量,具体为:
确定载波相位的变化值,即Δt时间内本地振荡与接收信号相位的差值与多普勒频率在Δt时间内的增量之和,具体为:
其中/>
其中,为载波多普勒相位差,ωe为载波多普勒频率差,/>为载波多普勒频率变化率;
将载波相位状态转移方程表示为:
其中,Wn为均值为0的高斯白噪声;
中间支路两路解扩信号与正交支路经过载波鉴相器,在鉴相更新周期内得到相位差的平均值:
其中,为相位差的平均值,Ts为载波鉴相更新周期;
确定观测矩阵中的各个元素,观测矩阵为Hf:
确定载波卡尔曼滤波状态矩阵φf与观测矩阵Hf分别为:
将计算得到的状态矩阵φf与观测矩阵Hf代入载波卡尔曼滤波的动态模型进行滤波处理,其中,载波卡尔曼滤波的动态模型具体为:
Xf(k+1)=φfXf(k)+Wf(k)
Yf(k)=HfXf(k)+Vf(k)
其中
6.根据权利要求1所述的基于卡尔曼滤波的载波和伪码的同步方法,其特征在于,将伪码鉴相处理得到的信号作为伪码卡尔曼滤波的输入,将作为伪码卡尔曼滤波器的状态变量进行卡尔曼滤波的具体方法为:
确定Δt时间内伪码的多普勒频率产生的伪码相位增量为:
Δt时间内产生的伪码相位变化值为本地伪码与接收伪码的相位的差值加上伪码多普勒频率在该时间内的增量:
其中/>
式中,ωτe为伪码频率差,为伪码频率变化率;
根据Δt时间内产生的伪码相位变化值,确定伪码相位状态转移方程为:
其中Wn为均值为0的高斯白噪声;
解扩后超前支路和同相支路经过伪码鉴相器,其输出为鉴相更新周期内得到的码相位差的平均值:
式中,Ts为更新周期;
将伪码的状态矩阵φτ与观测矩阵Hτ代入卡尔曼滤波公式,获得滤波结果,其中,伪码的状态矩阵φτ与观测矩阵Hτ分别为:
7.根据权利要求1所述的基于卡尔曼滤波的载波和伪码的同步方法,其特征在于,将步骤3中的两个卡尔曼滤波的状态变量值与协方差矩阵进行线性最小方差融合估计的具体过程为:
步骤5-1、将伪码卡尔曼滤波输出的伪码相位偏移映射为载波对应的多普勒相位偏移,具体过程如下:
伪码卡尔曼滤波输出为根据鉴相算法/>其中A为信号幅度增益,计算出实际的伪码偏移大小为/>
确定伪码频率偏移映射出的载波的多普勒频率偏移为:
其中,fc表示载波频率,Rc为PN码的码元速率,fτd表示伪码多普勒偏移频率,Rc1为实际的伪码速率,Nc为扩频码的长度,ε为实际的伪码偏移,Tc为伪码的码元周期;
根据载波的鉴相输出为sin(2πfdkTb),计算出由伪码映射出的载波多普勒相位偏移近似于2πfd1kTb,即将伪码卡尔曼滤波输出映射为载波卡尔曼滤波输出:
其中,为伪码卡尔曼滤波状态变量映射为载波卡尔曼滤波对应的状态变量,其协方差矩阵为/>
步骤5-2、将两个卡尔曼滤波输出结果进行线性最小方差融合估计,具体计算步骤如下:
线性最小方差融合估计:
对一个任意随机向量X∈Rn,假设有L个无偏独立估计对应的协方差矩阵为则/>就是X的线性融合估计,其中/>融合估计的协方差矩阵:
将伪码卡尔曼滤波输出映射为载波卡尔曼滤波的状态变量的值与载波卡尔曼滤波的输出Xf作为线性最小方差融合估计的输入变量,/>与Pf为对应的协方差矩阵;
融合估计步骤如下:
其中为线性最小方差融合估计值,/>为对应的协方差矩阵;
将融合估计输出作为载波卡尔曼滤波预测步的输入,并通过步骤4-1的逆变换计算出伪码卡尔曼滤波预测步的状态变量更新值与协方差矩阵更新值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311562977.8A CN117614786A (zh) | 2023-11-21 | 2023-11-21 | 基于卡尔曼滤波的载波和伪码的同步方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311562977.8A CN117614786A (zh) | 2023-11-21 | 2023-11-21 | 基于卡尔曼滤波的载波和伪码的同步方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117614786A true CN117614786A (zh) | 2024-02-27 |
Family
ID=89947325
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311562977.8A Pending CN117614786A (zh) | 2023-11-21 | 2023-11-21 | 基于卡尔曼滤波的载波和伪码的同步方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117614786A (zh) |
-
2023
- 2023-11-21 CN CN202311562977.8A patent/CN117614786A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113238261B (zh) | 低轨卫星扩频通信体制信号捕获跟踪系统 | |
JP3013763B2 (ja) | キャリア同期ユニット | |
US5898665A (en) | Coherent tracking apparatus and method for CDMA receiver | |
CN110912847B (zh) | 一种gmsk信号解调方法 | |
CN108768604B (zh) | 一种用于pcm/fm多符号检测的低复杂度位同步方法 | |
US6768780B1 (en) | Non-data aided maximum likelihood based feedforward timing synchronization method | |
CN108270715B (zh) | 适合高阶4096-qam的载波恢复系统和方法 | |
KR102341877B1 (ko) | 송신기와 수신기 및 해당 방법들 | |
CN109756968B (zh) | 单载波扩频系统的精确同步定时方法及精确同步定时装置 | |
US6430235B1 (en) | Non-data-aided feedforward timing synchronization method | |
CN113037671B (zh) | 一种低复杂度的高效soqpsk符号定时与相位联合同步算法 | |
WO1993022861A1 (en) | Receiver for digital communication system | |
JPH0131821B2 (zh) | ||
EP0838111B1 (en) | Carrier recovery for digitally phase modulated signals, using a known sequence | |
US7277509B2 (en) | Low complexity frequency-offset correction method | |
TW201412062A (zh) | 數位資料輔助頻率偏移估計 | |
JP3859903B2 (ja) | 周波数誤差推定装置およびその方法 | |
US5337331A (en) | Method of coherent demodulation for phase shift keying and a device for implementing the method | |
CN117614786A (zh) | 基于卡尔曼滤波的载波和伪码的同步方法 | |
CN113644934B (zh) | 一种星地异源扩跳频载波捕获频率补偿方法及系统 | |
CN114465691A (zh) | 一种低复杂度的恒包络相位调制信号采样偏差估计和补偿方法及系统 | |
US4263671A (en) | Sampling clock correction circuit | |
Shammakhi et al. | GPU-Based Parallel Algorithm for Wideband Signal Timing Recovery | |
CN115865301B (zh) | 一种卫星通信中抗大码偏、低功耗的定时同步方法 | |
CN114244673B (zh) | 一种基于数据辅助的频率跟踪方法、系统及存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |