专利名称：一种多相关值定宽拟合伪码跟踪鉴相方法
技术领域：
本发明涉及的是一种信号处理方法，具体地说是一种在实际快时变低信噪比等 复杂信道下，实现高稳定性伪码跟踪环路伪码相位差提取的方法。(二) 背景技术扩频体制因具有抗干扰能力强、保密性强、可实现码分多址和高精度测量的优 点，广泛应用于各种通信、导航系统中，如CDMA、 GPS等。然而在无线电各个频段 内，存在着各式各样的人为或非人为的干扰，这将会严重的影响伪码跟踪环路的鉴 相特性曲线，使其不再呈现理想情况下的单峰值三角特性相关曲线，而是时变的存 在多处抖动的复杂相关曲线，非常不利于伪码跟踪环路接收信号伪码相位差的准确 提取。图1为在理想环境下伪码自相关峰，图2为在复杂快时变低信噪比环境下伪 码自相关峰，从两图中可以看出，在这种复杂信道下，常规伪码跟踪环路采用的超 前滞后鉴相方式，通过超前滞后两个支路做差来得到关于伪码相位差的观测量，得 到的相位差将抖动非常大，信号跟踪点可能锁定在存在较大误差的相位点上，容易 偏离真实伪码相位差，从而给伪码跟踪带来较大的跟踪相位误差。并且实际相关曲 线是关于相关峰值平方的曲线，由于噪声的存在，相关值总是非零的，因此采用双 曲线拟合的鉴相方式会引入负值估计的相关峰值，对相关曲线中非零相关长度为2 个码片长度未做限制，将会引入一定误差，因此需要研究复杂环境下的伪码鉴相方 式。扩频接收机数字信号处理部分一般均采用FPGA+DSP工作方式，这种方式用于 扩频导航信号的捕获跟踪是公认已知的，例如Artech House出版的，题为"GPS Principles and Applications (GPS原理与应用)"(ISBN 0-89006-793-7) —书 中第5章(GPS satellite signal acquisition and tracking "GPS卫星信号捕 获和跟踪")。接收机中各通道的数字相关器、积分清洗器、伪码发生器、载波发生 器及逻辑控制电路等由FPGA芯片來完成，DSP芯片读取各通道积分清洗器的值， 依据相应的控制算法控制各通道的硬件电路，并监控环路的锁定情况，同时完成数 据采集、载体位置解算等任务。采用多个通道并行工作以同时跟踪不同导航台信号。扩频导航信号的捕获跟踪环路采用公认的基于能量判决的同相lH交支路实现环路 的闭合，例如科学出版社出版的，题为"扩频通信数字基带信号处理算法及其VLSI 实现"(ISBN7-03-013158-4) —书中第7章(扩频码序列的捕获)和第8章(数字 延迟锁定跟踪环)。经过接收信号的]f.确捕获后，本地伪码发生器与接收信号的伪码相位差控制在 半个码片之内。接收机将转入伪码跟踪，通过FPGA得到各个支路的相关值平方输 入到DSP中，由DSP中的伪码鉴相算法得到接收信号与本地信号相位差，相位差 再经过环路滤波器输出控制量取整后输入至FPGA，控制本伪码相位，从而准确跟 踪接收伪码相位。
发明内容本发明的目的在于提供一种能提高复杂环境下的伪码跟踪精度和稳定度，消除 噪声的影响，能够得到本地信号与接收信号伪码相位差的精确估计值的一种多相关 值定宽拟合伪码跟踪鉴相方法。本发明的目的是这样实现的本地伪码发生器产生多路超甜和滞后的本地伪码并分别与接收信号进行相关， 经过积分清洗和平方环节后，运用数据拟合鉴相方式得到本地对准码相位与接收信 号相位差，再经过环路滤波器后调整本地码相位。 '本发明还可以包括1、 所采用的相关数据必须关于对称支路两边对称。2、 所述的数据拟合是采用具有广义高斯函数形式的式^ = /(x)^^T^、'"进行采 样点数据拟合，本地伪码发生器产生的多路超甜滞后伪码相位分别对应函数中的自 变量x，而经过积分清洗器和平方环后得到的不同相位相关峰值则对应着^值，即利用(x。，x,…xj和对应的0;。,乂…少J估计参数a,6,c ，参数c则对应接收码相位与 本地码相位之差。3、 在拟合算法中引入宽度限制条件，并考虑到噪声门限F"。,、e,引入与噪声门限 相交的x,，^之间要保持两个码片宽度的限制条件；由式;;=/(X) = "e—WH" (1)得"/W = ",c>《)w (2) 将式(2)展丌得x12 =￡:土、/_丄111"^ (3)引入宽度限制条件得:x,-x2 =2」」ln^^ = 2 (4)从而得到a与6的关系 "《/ (5)将式(5)带入式(1)可得 "/W = L/—(6)由于式(6)为非线性方程，因此对式(6)拟合甜首先需对改式进行线性化处 理，采用泰勒级数展丌对方程(6)进行线性化处理；取6。,c。为待估计量6，c近似值，估计量6,c则可表示为p "o + & lc = c0 +《式中A,《为估计量与所取近似值的偏差； 将式(6)在&，c。处展丌并取一次项，则有、&=[1 + & _ (x - c0 )2 & + 2 V《- 2V(A. ] Z 式(7)中，6。,c。，P;,为己知量，xj为跟踪环路测量值，A，《为未知量，式 (7)可以写成、欲'y = /(x) = /(&，c。) + (^) m,《+ (^)|『《(7)其中令/ = 乙、/"'( 2 ' "L、/"( 2[1-(X-02]，则式(8)可以表示为 / = ^ + (9)利用最小二乘法对《，《进行估计，得< A [厅][;^]-[T][rp]其中[]为高斯总和符号'[，]=|>,2 =i[2F_V"( 2(x, '/=0 /=0[W] = i>,/,，其余与之类似；得到《,《的估计量后令{"=&+& ，则下一次最小二乘法估计中的估计量&c则可表示为[c = " +《再将式(6)在"，^处泰勒级数展开，并利用最小二乘法估计出《,《，这样经过反复迭代运算稳定后，实现对6,c的准确估计，从而得到本地伪码与接收信号的相位差，将相位差送入环路滤波器，得到对伪码时钟的控制量，进一歩控制本地伪 码发生器，使本地伪码进一歩与接收信号相位对齐。本发明针对采用扩频体制的无线电通信、导航系统提出一种伪码广义高斯曲线 拟合鉴相方法，利用广义高斯曲线来拟合多支路不同延迟相位的输入相关曲线，通 过将非线性最小二乘算法线性化來估计本地伪码发生器与接收信号的伪码相位差， 从而实现伪码跟踪环路的高精度及稳定度跟踪。本发明针对扩频系统接收机伪码跟踪环提出了一种伪码相位差鉴相方法，通过本地伪码发生器产生的多路不同延迟的伪码信号与接收信号进行相关得到各路相 关值，然后利用广义高斯曲线去拟合相关值得到最大值对应的伪码相位差，从而完 成本地与接收伪码相位差的高精度和稳定度提取。发明中通过观察广义高斯曲线中参数6,c的估计值的收敛曲线来验证此算法的可行性，然后通过比较广义高斯曲线拟合鉴相方式与常规超甜滞后鉴相方式的伪码相位差收敛曲线来验证广义高斯曲 线拟合鉴相方式的高精度和稳定度的优点。图6和给出了这两种鉴相方式下的码跟 踪环路控制曲线，由图对比可知，采用本发明所述的多相关值定宽拟合鉴相方法的 环路跟踪控制更加稳定，并具有较小的抖动。图7和图8给出了两种方式下的伪码 跟踪环相位跟踪误差曲线对比图，从图中可以看出，采用本发明所述的多相关值定 宽拟合鉴相方法的环路跟踪误差较常规超甜滞后鉴相方式减小一半左右，明显提高 了相位跟踪精度。(四)
图1是理想环境下伪码信号自相关峰； 图2是复杂环境下伪码信号自相关峰；图3为不同时延相关值拟合鉴相方式的伪码跟踪原理框图； 图4是本地码片相位采样点示意图； 图5是相关峰值拟合宽度限制示意图；图6为传统超前滞后鉴相方式和多相关值定宽拟合鉴相方式下伪码跟踪环路控 制量曲线对比图；图7为传统超前滞后鉴相方式和多相关值定宽拟合鉴相方式下伪码跟踪误差曲 线对比图；图8为图7的局部放大图。
具体实施方式
下面结合附图举例对本发明做更详细地描述接收机扩频信号进行捕获，使本地产生的伪码相位与接收信号伪码相位差在半 个码片之内，之后转入伪码跟踪。转入伪码跟踪后，通过本地伪码发生器产生的多 路超前滞后对称的不同时延本地伪码信号与接收信号进行相关和平方运算得到各 路相关峰值，然后利用上述广义高斯曲线去拟合接收信号的相关曲线，通过最小二乘法不断迭代修改广义高斯曲线中参数估计值，经过多次迭代稳定后，将此时提取 的伪码相位差观测量经过环路滤波器送入到伪码时钟控制器中，调整本地伪码相位 使之与接收信号的伪码相位差趋向于零。 -参数估计值迭代稳定前需不断调整广义高斯曲线的参数值，并且此时得到的伪 码相位差观测量不反馈给伪码时钟控制器，只有参数估计值当迭代稳定后，才利用 提取的伪码相位差观测量控制整个伪码跟踪环路，并且不再修改广义高斯曲线的参 数值。
权利要求
1、一种多相关值定宽拟合伪码跟踪鉴相方法，其特征是本地伪码发生器产生多路超前和滞后的本地伪码并分别与接收信号进行相关，经过积分清洗和平方环节后，运用数据拟合鉴相方式得到本地对准码相位与接收信号相位差，再经过环路滤波器后调整本地码相位。
2、 根据权利要求l所述的一种多相关值定宽拟合伪码跟踪鉴相方法，其特 征是所采用的相关数据必须关于对称支路两边对称。
3、 根据权利要求1或2所述的一种多相关值定宽拟合伪码跟踪鉴相方法， 其特征是所述的数据拟合是采用具有广义高斯函数形式的式y = /(X) = "e—W—。2进行采样点数据拟合，本地伪码发生器产生的多路超附带后伪码相位分别对应函数中的自变量;c ，而经过积分清洗器和平方环后得到的不 同相位相关峰值则对应着y值，即利用^。,x,…xJ和对应的UvnJ估计参 数a，6,c ，参数c则对应接收码相位与本地码相位之差。
4、 根据权利要求3所述的一种多相关值定宽拟合伪码跟踪鉴相方法，其特征是在拟合算法中引入宽度限制条件，并考虑到噪声门限^_,，引入与噪声门 限相交的X，，A之间要保持两个码片宽度的限制条件；由式;;=/(" = "^"—c)2 (1)得y = /(X) = W—"K)、^股(2) 将式(2)展丌得引入宽度限制条件得:一 x， = 2j——ln从而得到"与6的关系:<formula>formula see original document page 3</formula>由于式(6)为非线性方程，因此对式(6)拟合前首先需对该式进行线性 化处理，采用泰勒级数展丌对方程(6)进行线性化处理；取^,c。为待估计量6，c近似值，估计量6，c则可表示为<formula>formula see original document page 3</formula>式中《,《为估计量与所取近似值的偏差;将式(6)在6。,C。处展开并取一次项，则有 "胁爆。)+ (盖"《+(芸)ICc<formula>formula see original document page 3</formula>式(7)中，6。,c。，F,为已知量，x,y为跟踪环路测量值，《，《为未知量， 式(7)可以写成<formula>formula see original document page 3</formula>其中令/ ="<formula>formula see original document page 3</formula>利用最小二乘法对A，《,进行估计，得其中[]为高斯总和符号，<formula>formula see original document page 4</formula>其余与之类似；得到&，《的估计量后令j^
5 <image>image see original document page 4</image>
6，c则可表示为 <formula>formula see original document page 4</formula>再将式(6)在^,4处泰勒级数展开，并利用最小二乘法估计出《,《，这样 经过反复迭代运算稳定后，实现对6,c的准确估计，从而得到本地伪码与接收信号的相位差，将相位差送入环路滤波器，得到对伪码时钟的控制量，进一歩控 制本地伪码发生器，使本地伪码进一歩与接收信号相位对齐。<formula>formula see original document page 4</formula>，则下一次最小二乘法估计中的估计j
全文摘要
本发明提供的是一种多相关值定宽拟合伪码跟踪鉴相方法。本地伪码发生器产生多路超前和滞后的本地伪码并分别与接收信号进行相关，经过积分清洗和平方环节后，运用数据拟合鉴相方式得到本地对准码相位与接收信号相位差，再经过环路滤波器后调整本地码相位。本发明针对采用扩频体制的无线电通信、导航系统提出一种伪码广义高斯曲线拟合鉴相方法，利用广义高斯曲线来拟合多支路不同延迟相位的输入相关曲线，通过将非线性最小二乘算法线性化来估计本地伪码发生器与接收信号的伪码相位差，从而实现伪码跟踪环路的高精度及稳定度跟踪。
文档编号H04B1/69GK101247138SQ200810064170
公开日2008年8月20日 申请日期2008年3月25日 优先权日2008年3月25日
发明者徐定杰, 锋 沈, 伟 王, 冰 薛, 郝燕玲, 韦金辰, 平 黄 申请人:哈尔滨工程大学