CN101950024B - 用于局域增强系统的码载一致性检测方法 - Google Patents

Abstract

本发明公开了一种应用于局域增强系统的码载一致性检测方法，主要解决目前码载一致性检测中码载偏离度平滑精度低、平滑时延大的问题。其检测步骤是：首先对获得的载波相位观测量进行周跳检测，获得消除周跳的载波相位观测量和码伪距观测量；接着，计算码载偏离度和码载偏离度变化率；然后，对码载偏离度变化率进行小波分析，获得平滑的码载偏离度变化率并从中分离出电离层延迟变化率、多径延迟变化率以及接收机噪声残留变化率分量；最后，根据三种分量数据分别进行电离层异常、多径异常和接收机噪声异常检验，生成码载一致性检验结果矩阵。本发明具有平滑精度、无拉伸过程和时间延迟小的优点，可用于卫星导航局域增强系统的完好性监测。

Description

用于局域增强系统的码载一致性检测方法

技术领域

本发明属于卫星导航领域，涉及卫星导航领域的局域增强系统，特别是一种应用于卫星导航局域增强系统的码载一致性检测方法。

背景技术

卫星导航局域增强系统是一个地基差分系统，通过地面差分站发播卫星差分数据，用户将差分数据用于解算来获得高精度的导航参数，如速度、位置和时间。该系统采用的差分技术是利用多个参考接收机测量的位置信息或距离信息及其他导航信息的相关性能消除大部分卫星钟误差、星历误差及大气延迟误差，大大提高了导航定位精度。普遍认为，它可以满足非精密进近以及I、II直至III类精密进近着陆要求。

在卫星信号传递的过程中存在电离层、对流层折射引起的信号时延，在信号接收的过程中还存在由于天线周围障碍物环境和地面环境造成的多路径信号干扰和其他射频信号干扰，在接收机解调导航信号的过程中还存在接收机的噪声干扰，由于为用户播发导航数据的卫星处于高空之中，其广播信号在传播到地面或近地面的时候信号强度已经非常微弱，更加剧了上述因素造成的传播延迟、观测精度下降等问题，这些对卫星导航局域增强系统完好性将产生非常大的影响。

局域增强系统的信号质量监测关注的就是卫星信号在传播和解调过程中存在的对定位精度有影响的各种误差源，其中最主要的一种方法就是码载一致性检测方法。码载一致性检测指的是通过对码载偏离度的变化率进行分析评估以检验接收机获取的码伪距观测量和载波相位观测量的一致性的方法，码载偏离度变化率包含了电离层延迟的变化率、天线多径残留的变化率以及接收机噪声残留，在对码载偏离度变化率进行评估的过程中，可以完成上述信号质量阶段最关键的三种误差来源的评估检测，从而为系统完好性服务提供保障。

目前国内外对码载一致性检测的研究并不多见，国内对码载一致性的研究处于起步阶段，没有成型的技术方法，而国外也只有少数研究机构有一定的深入研究，且在对码载偏离度进行估计的时候通常采用GMA(Geometric Moving Averaging)方法，这种方法得到的码载偏离度数据存在较为明显的时间延迟，而且精度也不太高，同时只能将之作为一个整体进行评估检验，无法分离电离层延迟、天线多径延迟和接收机噪声残留，因此在实际应用中无法满足更高条件的完好性要求。

发明内容

本发明的目的在于克服上述已有技术的不足，提供一种用于局域增强系统的码载一致性检测方法，该方法引入小波分析技术，通过小波分解与重构对码载偏离度变化率数据进行平滑，以有效减小检测延迟、提高检测精度，同时根据码载偏离度变化率数据三种误差组成的频率特性，将电离层延迟变化率、天线多径延迟变化率和接收机噪声残留变化率三种误差分量分离出来，并据此对卫星导航信号是否存在电离层异常、天线多径异常以及接收机噪声异常现象进行完整的检测，为局域增强系统提供更高等级的完好性服务。

为了实现上述目的，本发明提供一种用于卫星导航局域增强系统码载一致性检测的方法，包括以下步骤：

(1)获取历元k上的码伪距观测量和载波相位观测量利用周跳检验算法检验载波相位观测量中是否存在周跳，如果存在则消除周跳，得到无周跳的载波相位观测量；

(2)利用历元k上无周跳的载波相位观测量

和码伪距观测量

计算码载偏离度数值：

Z_CMC(k)＝2d_iono(k)+d_code-mp(k)+d_rn(k)-IA，

其中，d_iono(k)为历元k上码伪距中电离层延迟误差，d_code-mp(k)为历元k上码伪距中多径延迟误差，d_rn(k)为历元k上接收机噪声残留，IA表示整周模糊度的误差；

(3)利用历元k上的码载偏离度数值Z_CMC(k)和上一历元k-1上的码载偏离度数值Z_CMC(k-1)，计算码载偏离度变化率：

DEV_P-C(k)＝Z_CMC(k)-Z_CMC(k-1)＝2d_iono(k)+d_{code_mp}(k)+d_rn(k)，

其中，电离层延迟变化率； $d_{\mathrm{iono}}\left(k\right)=\frac{d_{\mathrm{iono}}\left(k\right)-d_{\mathrm{iono}}(k-1)}{T_s},$

多径延迟变化率： $d_{\mathrm{code}\_\mathrm{mp}}\left(k\right)=\frac{d_{\mathrm{code}\_\mathrm{mp}}\left(k\right)-d_{\mathrm{code}\_\mathrm{mp}}(k-1)}{T_s},$

接收机噪声残留变化率： $d_{\mathrm{rn}}\left(k\right)=\frac{d_{\mathrm{rn}}\left(k\right)-d_{\mathrm{rn}}(k-1)}{T_s};$

(4)采用小波分析方法对步骤(3)所得码载偏离度变化率DEV_P-C(k)数据，进行如下处理：

4a)对码载偏离度变化率数据进行小波分解： ${\mathrm{DEV}}_{P-C}\left(k\right)=a_l+\underset{i=1}{\overset{l}{\mathrm{\Σ}}}{d}_i,$

其中，l是小波分解的层次；a_l是在层次l上的平均部分分解系数，频率范围为0到(1/2^l)*(f_s/2)Hz；d_i是在层次i上的细节部分分解系数，频率范围为从(1/2^i-1)*(f_s/2)到(1/2ⁱ)*(f_s/2)Hz；f_s＝1/R_s是观测信号的采样频率[Hz]；R_s是数据采样间隔[秒]；

4b)根据局域增强系统覆盖范围内电离层频率特性、参考天线周围多径频率特性、接收机自身噪声频率特性以及步骤4a)中得到的小波分解结果，分别重构电离层延迟变化率时域分量DEV_iono(k)、天线多径延迟变化率DEV_mp(k)以及接收机噪声残留变化率DEV_rn(k)；

4c)分别将所述的DEV_iono(k)、DEV_mp(k)和DEV_rn(k)与通过统计分析所得的相应门限比对，如果大于相应的门限，则向系统提出异常告警，否则判断为正常状态；

(5)根据步骤(4)得到的三种异常检测判断结果，经过综合判定形成如下矩阵，该矩阵即为码载一致性检测结果，描述了系统当前码载一致性状态，用于局域增强系统完好性监视：

${\mathrm{CCD}}_{\mathrm{matrix}}=\left[\begin{array}{ccccc}{\mathrm{PRN}}_1& {\mathrm{IAR}}_1& {\mathrm{AMR}}_1& {\mathrm{RNR}}_1& {\mathrm{TR}}_1\\ {\mathrm{PRN}}_2& {\mathrm{IAR}}_2& {\mathrm{AMR}}_2& {\mathrm{RNR}}_2& {\mathrm{TR}}_2\\ \·& \·& \·& \·& \·\\ \·& \·& \·& \·& \·\\ \·& \·& \·& \·& \·\\ {\mathrm{PRN}}_{\mathrm{SN}}& {\mathrm{IAR}}_{\mathrm{SN}}& {\mathrm{AMR}}_{\mathrm{SN}}& {\mathrm{RNR}}_{\mathrm{SN}}& {\mathrm{TR}}_{\mathrm{SN}}\end{array}\right]$

其中，SN为导航系统中卫星数量，PRN_k表示卫星号，IAR_k表示电离层异常检测结果，AMR_k表示天线多径异常检测结果，RNR_k表示接收机噪声异常检测结果，TR_k为综合检测结果，0表示码载一致性正常，1表示一种组成分量出现异常，2表示两种组成分量出现异常，3表示三种组成分量全部出现异常，1≤k≤SN。

本发明与现有技术相比具有如下优点：

1、本发明由于使用了小波平滑方法，提高了码载偏离度变化率数据的精确度，减小了平滑延迟时间，消除了平滑中拉伸的过程；

2、本发明由于了采用小波重构的方法，根据三种误差组成的频率特性从码载偏离度变化率中将电离层延迟变化率、天线多径延迟变化率和接收机噪声残留变化率孤立出来，实现了对电离层异常、天线多径异常以及接收机噪声异常的分项检测；

3、本发明步骤5)中得到的码载一致性检测结果矩阵，可以提供分等级的异常检测结果，能够为局域增强系统码载一致性检测提供更高级的完好性服务。

附图说明

图1为本发明的总流程图；

图2为本发明的周跳检验子流程图；

图3为本发明的小波分析子流程图；

图4为本发明中小波平滑所得结果与经典GMA方法所得结果精度对比图。

具体实施方式

下面将结合附图对本发明作进一步的详细说明。

参照图1，本发明包括以下步骤：

步骤1，获取历元k上的码伪距观测量和载波相位观测量

并对载波相位观测量进行周跳检测与消除。

参照图2，本步骤的具体实现如下：

(2A)计算历元k上的载波校正量

${\mathrm{\φ}}_c\left(k\right)={\mathrm{\ρ}}_{\mathrm{phase}}^{\mathrm{sv}}\left(k\right)-R\left(k\right)+\mathrm{\τ}\left(k\right)-\mathrm{\φ}\left(0\right),$

其中，

是载波相位观测量；R(k)是天线与卫星间的真距；τ(k)是卫星钟校正量；初始的载波相位校正量φ(0)等于φ(0)-R(0)+τ(0)，并且消除了整周模糊度；

(2B)根据载波校正量的二次模型：

利用连续三个时刻的数据，通过最小二乘法解算系数 和

(2C)设

Acc为载波相位加速度，如果Acc介于0到1之间，则认为没有发生周跳，否则，认为发生了周跳。

步骤2，利用历元k上无周跳的载波相位观测量

和码伪距观测量

计算码载偏离度：

$Z_{\mathrm{CMC}}={\mathrm{\ρ}}_{\mathrm{code}}^{\mathrm{sv}}-{\mathrm{\ρ}}_{\mathrm{phase}}^{\mathrm{sv}}=d_{\mathrm{iono}}+d_{\mathrm{code}-\mathrm{teack}}+d_{\mathrm{code}-\mathrm{mp}}+d_{\mathrm{code}-\mathrm{noise}}$

$-(-d_{\mathrm{iono}}+d_{\mathrm{phase}-\mathrm{track}}+d_{\mathrm{phase}-\mathrm{mp}}+d_{\mathrm{phase}-\mathrm{noise}}+\mathrm{IA})+(\mathrm{GroupDelay}-\mathrm{PhaseDelay}),$

$=2d_{\mathrm{iono}}+d_{\mathrm{code}-\mathrm{mp}}+d_{\mathrm{rn}}-\mathrm{IA}$

其中，

为用户至卫星的伪距测量值，为卫星至监测站的几何距离，Δt_ref为监测站接收机时钟相对于基准时间系统的时间偏差，d_iono为电离层延迟，d_tropo表示对流层延迟，URE为卫星星历引入的误差，d_code-track为接收机运动引入的码跟踪误差，d_code-mp为码伪距的多径误差，d_code-noise为码噪声误差，GroupDelay为码的群延迟，

为用户至卫星的伪距测量值，

为卫星至监测站的几何距离，Δt_ref为监测站接收机时钟相对于基准时间系统的时间偏差，Δt^sv为卫星时钟相对于基准时间系统的时间偏差，d_iono为电离层延迟，d_tropo为对流层延迟，URE为卫星星历引入的误差，d_phase-track为接收机运动引入的载波相位跟踪误差，d_phase-mp为载波的多径误差，d_phase-noise为载波噪声误差，IA为整周模糊度的误差，PhaseDelay为载波的相位延迟，d_rn＝d_code-noise-d_phase-noise为接收机噪声残留变化率。

步骤3，利用历元k上的码载偏离度数值Z_CMC(k)和上一历元k-1上的码载偏离度数值Z_CMC(k-1)，计算码载偏离度变化率：

DEV_P-C＝Z_CMC(k)-Z_CMC(k-1)

＝2(d_iono(k)-d_iono(k-1))+(d_{code_mp}(k)-d_{code_mp}(k-1))+(d_rn(k)-d_rn(k-1))，

＝2d_iono+d_{code_mp}+d_rn

其中，电离层延迟变化率

多径延迟变化率

接收机噪声残留变化率

步骤4，采用小波分析方法对步骤3所得码载偏离度变化率DEV_P-C(k)数据进行处理。

参照图3，本步骤的具体实现如下：

(4A)对码载偏离度变化率数据进行小波分解：

${\mathrm{DEV}}_{P-C}=a_l+\underset{i=1}{\overset{l}{\mathrm{\Σ}}}{d}_i,$

其中，l是小波分解的层次；a_l是在层次l上的平均部分分解系数，频率范围为0到(1/2^l)*(f_s/2)Hz；d_i是在层次i上的细节部分分解系数，频率范围为从(1/2^i-1)*(f_s/2)到(1/2ⁱ)*(f_s/2)Hz；f_s＝1/R_s是观测信号的采样频率[Hz]；R_s是数据采样间隔[秒]；

(4B)根据步骤(4A)所得分解系数重构码载偏离度变化率三种组成分量：

(4B1)根据局域增强系统覆盖范围内电离层频率特性，选择低通滤波器、高通滤波器以及重构所需分解系数，从所选最低层次开始，重构电离层延迟变化率时域分量DEV_iono(k)：

(4B11)对该最低层次平均部分分解系数进行上抽样，并求解与低通滤波器的低通卷积；

(4B12)对该最低层次细节部分系数进行上抽样，并求解与高通滤波器的高通卷积；

(4B13)将步骤(4B11)和(4B12)中得到的低通卷积和高通卷积求和，得到该层次重构结果，当重构层次小于等于选择的最大层次时，将重构结果作为下一层的平均部分系数与下一层细节部分系数一起重复步骤(4B11)和(4B12)；当重构层数大于选择的最大层数时，重构结束，得到的平均部分即为电离层延迟变化率时域分量DEV_iono(k)；

(4B2)根据参考天线周围多径频率特性，选择低通滤波器、高通滤波器以及重构所需分解系数，从所选最低层次开始，重构天线多径变化率时域分量DEV_mp(k)：

(4B21)对该最低层次平均部分分解系数进行上抽样，并求解与低通滤波器的低通卷积；

(4B22)对该最低层次细节部分系数进行上抽样，并求解与高通滤波器的高通卷积；

(4B23)将步骤(4B21)和(4B22)中得到的低通卷积和高通卷积求和，得到该层次重构结果，当重构层次小于等于选择的最大层次时，将重构结果作为下一层的平均部分系数与下一层细节部分系数一起重复步骤(4B21)和(4B22)；当重构层数大于选择的最大层数时，重构结束，得到的平均部分即为天线多径变化率时域分量DEV_mp(k)；

(4B3)根据接收机自身噪声频率特性，选择低通滤波器、高通滤波器以及重构所需分解系数，从所选最低层次开始，重构接收机噪声变化率时域分量DEV_rn(k)：

(4B31)对该最低层次平均部分分解系数进行上抽样，并求解与低通滤波器的低通卷积；

(4B32)对该最低层次细节部分系数进行上抽样，并求解与高通滤波器的高通卷积；

(4B33)将步骤(4B31)和(4B32)中得到的低通卷积和高通卷积求和，得到该层次重构结果，当重构层次小于等于选择的最大层次时，将重构结果作为下一层的平均部分系数与下一层细节部分系数一起重复步骤(4B31)和(4B32)；当重构层数大于选择的最大层数时，重构结束，得到的平均部分即为接收机噪声变化率时域分量DEV_rn(k)；

(4C)分别将所述的DEV_iono(k)、DEV_mp(k)和DEV_rn(k)与统计分析所得相应门限进行比对，如果大于相应的门限，则向系统提出异常告警，否则判断为正常状态。

步骤5，根据步骤4得到的三种异常检测判断结果，经过综合判定形成如下矩阵，该矩阵即为码载一致性检测结果：

${\mathrm{CCD}}_{\mathrm{matrix}}=\left[\begin{array}{ccccc}{\mathrm{PRN}}_1& {\mathrm{IAR}}_1& {\mathrm{AMR}}_1& {\mathrm{RNR}}_1& {\mathrm{TR}}_1\\ {\mathrm{PRN}}_2& {\mathrm{IAR}}_2& {\mathrm{AMR}}_2& {\mathrm{RNR}}_2& {\mathrm{TR}}_2\\ \·& \·& \·& \·& \·\\ \·& \·& \·& \·& \·\\ \·& \·& \·& \·& \·\\ {\mathrm{PRN}}_{\mathrm{SN}}& {\mathrm{IAR}}_{\mathrm{SN}}& {\mathrm{AMR}}_{\mathrm{SN}}& {\mathrm{RNR}}_{\mathrm{SN}}& {\mathrm{TR}}_{\mathrm{SN}}\end{array}\right],$

其中，SN为导航系统中卫星数量，PRN_k表示卫星号，IAR_k表示电离层异常检测结果，AMR_k表示天线多径异常检测结果，RNR_k表示接收机噪声异常检测结果，TR_k为综合检测结果，0表示码载一致性正常，1表示一种组成分量出现异常，2表示两种组成分量出现异常，3表示三种组成分量全部出现异常，1≤k≤SN。

本发明的效果可由如下实例验证：

选择局域增强系统地面站参考接收机与配套天线作为实测对象，天线位于建筑物楼顶，距楼顶平面高度为约1.5米，周围无明显遮挡，根据长期统计分析结果显示，三种主要误差组成电离层延迟误差、多径延迟误差和接收机噪声中，电离层误差拥有最低的频谱，变化范围为0-1.4e-4Hz；在采样频率为1Hz的情况下，多径误差的频率谱组成范围从大约0.003到0.02Hz；接收机噪声频谱大体上是高斯分布，可以通过一个1-2-1分解因子在细节层次“1”上将接收机噪声分解并隔离出来，将三个重构的组成分量组合得到码载偏离度变化率平滑结果。采用局域增强系统中参考接收机对GPS卫星信号进行持续检测，以第2号卫星为例，通过计算码载偏离度变化率数据，并对该数据分别进行GMA平滑和小波平滑，对结果进行绘图显示。

图4给出了2号卫星GPS 1743周内一段时间的码载偏离度变化率平滑结果轨迹，数据区间长度为2号星一个可见周期时长，图4A)为小波平滑结果轨迹，图4B)为GMA平滑结果，其平滑区间长度为200秒，由图4A)和图4B)可知，GMA平滑结果具有一个较明显的拉伸过程，且精度不高，无法精确体现实际码载偏离度变化率的波形，而小波平滑的结果与之相比没有拉伸过程且精度高，能够比较准确地体现码载偏离度变化率的变化情况。图4C)为小波平滑结果的局部放大轨迹，图4D)为相对应GMA平滑结果的局部放大轨迹，由图4C)和图4D)可知，GMA平滑结果具有较明显的时间延迟，而小波平滑所得结果延时几乎可以忽略。此结果标明，本方法可以提高码载一致性检测的精度，增强系统对码载偏离度异常检测的准确星，降低局域增强系统中码载一致性产生的完好性风险。

Claims (5)

1.一种用于卫星导航局域增强系统的码载一致性检测方法，包括如下步骤： 

(1)获取历元k上的码伪距观测量

和载波相位观测量

利用周跳检验算法检验载波相位观测量中是否存在周跳，如果存在则消除周跳，得到无周跳的载波相位观测量； 

(2)利用历元k上无周跳的载波相位观测量

和码伪距观测量计算码载偏离度数值： 

Z_CMC(k)＝2d_iono(k)+d_code-mp(k)+d_rn(k)-IA， 

其中，d_iono(k)为历元k上码伪距中电离层延迟误差，d_code-mp(k)为历元k上码伪距中多径延迟误差，d_rn(k)为历元k上接收机噪声残留，IA表示整周模糊度的误差； 

(3)利用历元k上的码载偏离度数值Z_CMC(k)和上一历元k-1上的码载偏离度数值Z_CMC(k-1)，计算码载偏离度变化率： 

其中，T_s为历元时刻k与相邻历元时刻(k-1)之间的时间间隔，单位为秒， 

电离层延迟变化率：

多径延迟变化率：

接收机噪声残留变化率：

(4)采用小波分析方法对步骤(3)所得码载偏离度变化率DEV_P-C(k)数据，进行如下处理： 

4a)对码载偏离度变化率数据进行小波分解：

其中，l是小波分解的层次；a_l是在层次l上的平均部分分解系数，频率范围为0到(1/2^l)*(f_s/2)Hz；d_i是在层次i上的细节部分分解系数，频率范围为从 (1/2^i-1)*(f_s/2)到(1/2ⁱ)*(f_s/2)Hz；f_s＝1/R_s是观测信号的采样频率[Hz]；R_s是数据采样间隔[秒]； 

4b)根据局域增强系统覆盖范围内电离层频率特性、参考天线周围多径频率特性、接收机自身噪声频率特性以及步骤4a)中得到的小波分解结果，分别重构电离层延迟变化率时域分量DEV_iono(k)、天线多径变化率时域分量DEV_mp(k)以及接收机噪声变化率时域分量DEV_rn(k)； 

4c)分别将所述的DEV_iono(k)、DEV_mp(k)和DEV_rn(k)与通过统计分析所得的相应门限比对，如果大于相应的门限，则向系统提出异常告警，否则判断为正常状态； 

(5)根据步骤(4)得到的三种异常检测判断结果，经过综合判定形成如下矩阵，该矩阵即为码载一致性检测结果： 

其中，SN为导航系统中卫星数量，PRN_m表示卫星号，IAR_m表示电离层异常检测结果，AMR_m表示天线多径异常检测结果，RNR_m表示接收机噪声异常检测结果，TR_m为综合检测结果，0表示码载一致性正常，1表示一种组成分量出现异常，2表示两种组成分量出现异常，3表示三种组成分量全部出现异常，1≤m≤SN。 

2.根据权利要求1所述的码载一致性检测方法，其中步骤(1)所述的利用周跳检验算法检验载波相位观测量中是否存在周跳，按如下步骤进行： 

2a)计算历元k上的载波校正量

(k)： 

其中，

是载波相位观测量，R(k)是天线与卫星间的真距，τ(k)是卫星钟校正量，初始的载波相位校正量φ(0)等于φ(0)-R(0)+τ(0)，并且消除了整周模糊度； 

2b)根据载波校正量

的二次模型：利 用连续三个时刻的数据，通过最小二乘法解算系数

和

2c)设

Acc为载波相位加速度，如果Acc介于0到1之间，则认为没有发生周跳，否则，认为发生了周跳。 

3.根据权利要求1所述的码载一致性检测方法，其中步骤4b)所述的重构电离层延迟变化率时域分量DEV_iono(k)，是根据局域增强系统覆盖范围内电离层频率特性，选择低通滤波器和高通滤波器，从步骤4a)中小波分解所得到的各层次平均部分分解系数和细节部分分解系数选择与电离层频率特性相吻合的分解系数，从选择的最低层次开始，按如下步骤进行： 

3a)对该最低层次平均部分分解系数进行上抽样，并求解与低通滤波器的低通卷积； 

3b)对该最低层次细节部分分解系数进行上抽样，并求解与高通滤波器的高通卷积； 

3c)将步骤3a)和3b)中得到的低通卷积和高通卷积求和得到该层次重构结果，当重构层次小于等于选择的最大层次时，将重构结果作为下一层的平均部分分解系数与所选下一层细节部分分解系数一起再重复步骤3a)和3b)，直到重构层数大于选择的最大层数时，重构结束，得到的平均部分即为电离层延迟变化率时域分量DEV_iono(k)。 

4.根据权利要求1所述的码载一致性检测方法，其中步骤4b)所述的重构天线多径变化率时域分量DEV_mp(k)，是根据参考天线周围多径频率特性，选择低通滤波器和高通滤波器，从步骤4a)中小波分解所得到的各层次平均部分分解系数和细节部分分解系数选择与多径频率特性相吻合的分解系数，从所选择最低层次开始，按如下步骤进行： 

4a)对该最低层次平均部分分解系数进行上抽样，并求解与低通滤波器的低通卷积； 

4b)对该最低层次细节部分分解系数进行上抽样，并求解与高通滤波器的高通卷积； 

4c)将步骤4a)和4b)中得到的低通卷积和高通卷积求和得到该层次重构结果， 当重构层次小于等于选择的最大层次时，将重构结果作为下一层的平均部分分解系数与所选下一层细节部分分解系数一起再重复步骤4a)和4b)，直到重构层数大于选择的最大层数时，重构结束，得到的平均部分即为天线多径变化率时域分量DEV_mp(k)。 

5.根据权利要求1所述的码载一致性检测方法，其中步骤4b)所述的重构接收机噪声变化率时域分量DEV_rn(k)，是根据接收机自身噪声频率特性，选择低通滤波器和高通滤波器，从步骤4a)中小波分解所得到的各层次平均部分分解系数和细节部分分解系数选择与噪声频率特性相吻合的分解系数，从所选择最低层次开始，按如下步骤进行 

5a)对该最低层次平均部分分解系数进行上抽样，并求解与低通滤波器的低通卷积； 

5b)对该最低层次细节部分分解系数进行上抽样，并求解与高通滤波器的高通卷积； 

5c)将步骤5a)和5b)中得到的低通卷积和高通卷积求和得到该层次重构结果，当重构层次小于等于选择的最大层次时，将重构结果作为下一层的平均部分分解系数与所选下一层细节部分分解系数一起再重复步骤5a)和5b)，直到重构层数大于选择的最大层数时，重构结束，得到的平均部分即为接收机噪声变化率时域分量DEV_rn(k)。 

Images