一种多基准时钟加权合成方法

摘要

本发明公开了一种多基准时钟加权合成方法，解决了组合时钟精度和稳定度差的问题，所述方法包括以下步骤：从多个时钟源获得时钟信号；对所述时钟信号进行去除噪声合成，得到去除噪声的时钟信号；对所述去除噪声的时钟信号，先利用KNN分类算法进行分类，再对分类后的信号进行加权合成，得到时钟合成信号。在对所述时钟信号进行去除噪声合成，得到去除噪声的时钟信号之前，对所述时钟信号进行异常信号检测，去除所述时钟信号中的奇异点。本发明所实现的频率基准信号频率准确度在n×E‑14量级，优于高品质单个铯钟1～2个数量级；本发明可以同时改善时钟源的短期和长期稳定性。

权利要求

1.一种多基准时钟加权合成方法，其特征在于，包括以下步骤：从多个时钟源获得时钟信号；

对所述时钟信号进行去除噪声合成，得到去除噪声的时钟信号；

对所述去除噪声的时钟信号，先利用KNN分类算法进行分类，再对分类后的信号进行加权合成，得到时钟合成信号，本步骤进一步包含：将所述去除噪声的时钟信号进行两两比对，生成比对向量，并根据所述比对向量计算得到频偏矩阵；

根据所述频偏矩阵生成待测样本；

对所述待测样本进行分类，得到分类系数；

根据所述频偏矩阵计算所述时钟源的归一化权值；

根据所述分类系数和所述归一化权值进行时钟合成，得到时钟合成信号。

2.如权利要求1所述的一种多基准时钟加权合成方法，其特征在于，在对所述时钟信号进行去除噪声合成，得到去除噪声的时钟信号之前，对所述时钟信号进行异常信号检测，去除所述时钟信号中的奇异点。

3.如权利要求2所述的一种多基准时钟加权合成方法，其特征在于，对所述时钟信号进行异常信号检测，去除所述时钟信号中的奇异点，进一步包含：选取所述异常信号检测的测试参考信号，所述测试参考信号为所述时钟信号的任一历史输出值；

将所述时钟信号与所述测试参考信号进行比较，得到比较差值；

对所述比较差值进行判断，当所述比较差值大于奇异点判断阈值时，所述时钟信号为奇异点。

4.如权利要求2所述的一种多基准时钟加权合成方法，其特征在于，对所述时钟信号进行异常信号检测，去除所述时钟信号中的奇异点，进一步包含：将所述时钟信号进行两两比对，并计算得到两两比对差值；

根据所述两两比对差值将所述时钟信号分为两类，当所述两两比对差值大于奇异点判断阈值时，判定两个时钟信号为不同类别，当所述两两比对差值不大于所述奇异点判断阈值时，判定两个时钟信号为相同类别；

根据所述时钟信号的分类情况，将信号个数少的一类时钟信号判定为奇异点。

5.如权利要求1所述的一种多基准时钟加权合成方法，其特征在于，对所述时钟信号进行去除噪声合成，得到去除噪声的时钟信号，进一步包含：对所述时钟信号进行小波函数分解，得到小波分解细节系数；

对所述小波分解细节系数进行硬阈值判断筛选，得到筛选后的小波分解细节系数，所述硬阈值判断筛选为，当小波分解细节系数小于小波分解硬阈值时，将该小波分解细节系数置0，当小波分解细节系数大于或等于所述小波分解硬阈值时，保留该小波分解细节系数；

根据所述筛选后的小波分解细节系数，对所述时钟信号进行重构。

6.如权利要求1所述的一种多基准时钟加权合成方法，其特征在于，对所述时钟信号进行去除噪声合成，得到去除噪声的时钟信号，进一步包含：对所述时钟信号进行小波函数分解，得到小波分解细节系数；

对所述小波分解细节系数根据设定的软阈值筛选比例进行软阈值判断筛选，得到筛选后的小波分解细节系数，所述软阈值判断筛选为，当小波分解细节系数小于小波分解硬阈值时，将该小波分解细节系数置0，当小波分解细节系数大于或等于所述小波分解硬阈值时，将该小波分解细节系数减去阈值得到新的小波分解细节系数小波分解细节系数小波分解细节系数小波分解细节系数；

根据所述筛选后的小波分解细节系数，对所述时钟信号进行重构。

7.如权利要求1所述的一种多基准时钟加权合成方法，其特征在于，根据所述频偏矩阵计算所述时钟源的归一化权值为：其中μi为第i个时钟源的归一化权值，N为所述时钟源的个数，μi′为第i个时钟源的权值，μi′表示为：其中fij为第i台时钟源时钟信号和第j台时钟源时钟信号比对计算得到的频偏矩阵。

8.如权利要求1所述的一种多基准时钟加权合成方法，其特征在于，所述时钟合成信号为：其中，TA(t)为所述时钟合成信号，ai为所述分类系数，μi为第i个时钟源的归一化权值，TAi(t)为第i个时钟源获得的时钟信号。

9.如权利要求1所述的一种多基准时钟加权合成方法，其特征在于，对所述待测样本进行分类，得到分类系数采用的分类器为KNN分类器。

说明书

一种多基准时钟加权合成方法

技术领域

[0001] 本发明涉及卫星导航系统，特别涉及一种多基准时钟加权合成方法。

背景技术

[0002] 在卫星导航系统中需要应用高精度时钟进行导航和定位，而时钟的时间/频率精度是由原子钟的物理特性决定的，目前大多采用时钟合成的方法来提高时间/频率精度。现有技术中采用的时间合成方法包括：加权平均算法，Kalman滤波算法和小波分解算法，加权平均算法得到的时钟合成信号短期稳定性改善效果较差，Kalman滤波算法和小波分解算法得到的时钟合成信号对时钟准确度方面改善不明显。

发明内容

[0003] 本发明提供一种多基准时钟加权合成方法，解决了组合时钟精度和稳定度差的问题。

[0004] 一种多基准时钟加权合成方法，包括以下步骤：从多个时钟源获得时钟信号；对所述时钟信号进行去除噪声合成，得到去除噪声的时钟信号；对所述去除噪声的时钟信号，先利用KNN分类算法进行分类，再对分类后的信号进行加权合成，得到时钟合成信号。在对所述时钟信号进行去除噪声合成，得到去除噪声的时钟信号之前，对所述时钟信号进行异常信号检测，去除所述时钟信号中的奇异点。

[0005] 优选地，对所述时钟信号进行异常信号检测，去除所述时钟信号中的奇异点，进一步包含：选取所述异常信号检测的测试参考信号，所述测试参考信号为所述时钟信号的任一历史输出值；将所述时钟信号与所述测试参考信号进行比较，得到比较差值；对所述比较差值进行判断，当所述比较差值大于奇异点判断阈值时，所述时钟信号为奇异点。

[0006] 优选地，对所述时钟信号进行异常信号检测，去除所述时钟信号中的奇异点，进一步包含：将所述时钟信号进行两两比对，并计算得到两两比对差值；根据所述两两比对差值将所述时钟信号分为两类，当所述两两比对差值大于奇异点判断阈值时，判定两个时钟信号为不同类别，当所述两两比对差值不大于所述奇异点判断阈值时，判定两个时钟信号为相同类别；根据所述时钟信号的分类情况，将信号个数少的一类时钟信号判定为奇异点。

[0007] 进一步地，对所述时钟信号进行去除噪声合成，得到去除噪声的时钟信号，进一步包含：对所述时钟信号进行小波函数分解，得到小波分解细节系数对所述小波分解细节系数进行硬阈值判断筛选，得到筛选后的小波分解细节系数，所述硬阈值判断筛选为，当小波分解细节系数小于小波分解硬阈值时，将该小波分解细节系数置0，当小波分解细节系数大于或等于所述小波分解硬阈值时，保留该小波分解细节系数，根据所述筛选后的小波分解细节系数，对所述时钟信号进行重构。

[0008] 优选地，对所述时钟信号进行去除噪声合成，得到去除噪声的时钟信号，进一步包含：对所述时钟信号进行小波函数分解，得到小波分解细节系数；

[0009] 对所述小波分解细节系数根据设定的软阈值筛选比例进行软阈值判断筛选，得到筛选后的小波分解细节系数，所述软阈值判断筛选为，当小波分解细节系数小于小波分解硬阈值时，将该小波分解细节系数置0，当小波分解细节系数大于或等于所述小波分解硬阈值时，将该小波分解细节系数减去阈值得到新的小波分解细节系数小波分解细节系数小波分解细节系数小波分解细节系数；

[0010] 根据所述筛选后小波分解细节系数，对所述时钟信号进行重构。

[0011] 进一步地，本发明所述多基准时钟加权合成方法，对所述去除噪声的时钟信号，进行加权合成，得到时钟合成信号，进一步包含：将所述去除噪声的时钟信号进行两两比对，生成比对向量，并根据所述比对向量计算得到频偏矩阵；根据所述频偏矩阵生成待测样本；对所述待测样本进行分类，得到分类系数；

[0012] 根据所述频偏矩阵计算所述时钟源的归一化权值；根据所述分类系数和所述归一化权值进行时钟合成，得到时钟合成信号。

[0013] 优选地，作为本申请多基准时钟加权合成方法的实施例，根据所述频偏矩阵计算所述时钟源的归一化权值为

[0014]

[0015] 其中μi为第i个时钟源的归一化权值，N为所述时钟源的个数，μi′为第i个时钟源的权值，μi′表示为

[0016]

[0017] 其中fij为第i台时钟源时钟信号和第j台时钟源时钟信号比对计算得到的频偏矩阵，所述时钟合成信号为

[0018]

[0019] 其中，TA(t)为所述时钟合成信号，ai为所述分类系数，μi为第i个时钟源的归一化权值，TAi(t)为第i个时钟源获得的时钟信号。

[0020] 优选地，对所述待测样本进行分类，得到分类系数采用的分类器为KNN分类器。

[0021] 本发明有益效果包括：本发明所实现的频率基准信号频率准确度在n×E-14量级，优于高品质单个铯钟1～2个数量级；本发明可以同时改善时钟源的短期和长期稳定性，即频率准确度达到n×E-14量级的同时，观察窗口100s以下时的频率阿伦方差也比单个时钟源低2～3个数量级。

附图说明

[0022] 此处所说明的附图用来提供对本发明的进一步理解，构成本发明的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：

[0023] 图1为一种一种多基准时钟加权合成方法流程示意图；

[0024] 图2为一种包含异常信号检测的多基准时钟加权合成方法流程示意图。

具体实施方式

[0025] 为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明具体实施例及相应的附图对本发明技术方案进行清楚、完整地描述。显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

[0026] 以下结合附图，详细说明本发明各实施例提供的技术方案。

[0027] 图1为一种多基准时钟加权合成方法流程示意图。本申请实施例提出的一种多基准时钟加权合成方法，具体包括以下步骤：

[0028] 步骤101，从多个时钟源获得时钟信号。

[0029] 在步骤101中，所述时钟信号为所述时钟源的时间/频率信号。

[0030] 需要说明的是，选用多个时钟源的目的是组成组合时钟，所述时钟源的个数是预先设定的，要求不小于3个，具体个数这里不做限定。

[0031] 还需要说明的是，所述多个时钟源可以是原子钟，可以是GNSS接收机，这里不做具体限定。

[0032] 步骤102，对所述时钟信号进行去除噪声合成，得到去除噪声的时钟信号。

[0033] 在步骤102中，对所述时钟信号进行去除噪声合成，得到去除噪声的时钟信号，进一步包含：对所述时钟信号进行小波函数分解，得到小波分解细节系数；对所述小波分解细节系数进行硬阈值判断筛选，得到筛选后的小波分解细节系数，所述硬阈值判断筛选为，当小波分解细节系数小于小波分解硬阈值时，将该小波分解细节系数置0，当小波分解细节系数大于或等于所述小波分解硬阈值时，保留该小波分解细节系数小波分解细节系数小波分解细节系数；根据所述筛选后小波分解细节系数，对所述时钟信号进行重构。

[0034] 在步骤102中，M层小波函数分解可表示为

[0035]

[0036] 其中，TAi(t)表示第i台钟在t时刻的时钟信号，N为参加归算的时钟源总数，为小波分解的粗糙系数， 为小波分解的细节系数。保留所有的粗糙系数，对所述小波分解细节系数进行硬阈值判断筛选，得到筛选后的小波分解细节系数，所述硬阈值判断筛选为，当小波分解细节系数小于小波分解硬阈值(公式(2)中的Threshold)时，将该小波分解细节系数置0，当小波分解细节系数大于或等于所述小波分解硬阈值时，保留该小波分解细节系数，即：

[0037]

[0038] 根据阈值处理后的小波分解细节系数进行时钟信号的重构，即：

[0039]

[0040] 在步骤102中，对所述时钟信号进行去除噪声合成，得到去除噪声的时钟信号，进一步包含：对所述时钟信号进行小波函数分解，得到小波分解细节系数；对所述小波分解细节系数根据设定的软阈值筛选比例进行软阈值判断筛选，得到筛选后的小波分解细节系数，所述软阈值判断筛选为，当小波分解细节系数小于小波分解硬阈值时，将该小波分解细节系数置0，当小波分解细节系数大于或等于所述小波分解硬阈值时，将该小波分解细节系数减去阈值得到新的小波分解细节系数；根据所述筛选后的小波分解细节系数，对所述时钟信号进行重构。

[0041] 需要说明的是，对所述时钟信号进行小波函数分解，可以采用Haar小波、Daubechies小波、Coiflets小波等可以作为基函数，也可以采用其他小波作为基函数，这里不做限定。

[0042] 需要说明的是，理论上小波分解层数越多去噪精度越高，但是相应计算的复杂度也提升，计算时间越长，实际应用中在精度和计算复杂度上取一个折中数值。从纯理论的角度，M可以取任意正整数；从工程的角度考虑，一般我们将M取到3～10之间就可以达到比较好的效果。

[0043] 步骤103，对所述去除噪声的时钟信号，进行加权合成，得到时钟合成信号。

[0044] 在步骤103中，对所述去除噪声的时钟信号，进行加权合成，得到时钟合成信号，进一步包含：将所述去除噪声的时钟信号进行两两比对，生成比对向量Xij＝TAi'-TAj'，并根据所述比对向量计算得到频偏矩阵，偏移矩阵是一个N×N的矩阵，记为FN×N，这个矩阵中的每一个元素

[0045]

[0046] 其中i,j＝1,2,…,N；P是采样区间的时间序列长度，根据所述频偏矩阵生成待测样本(待测样本就是钟组中的每一台钟)；对所述待测样本进行分类，得到分类系数：按照每个时钟对系统贡献为正或负对TST进行分类，并赋予+1/-1的分类系数ai。本发明中定义两台钟对比频偏的绝对值为两台钟的距离，即dist＝|fij|)

[0047] 需要说明的是，KNN算法是一种经典的分类算法，核心思想是如果一个样本在特征空间中的k个最相邻(距离最小)的样本中的大多数属于某一个类别，则该样本也属于这个类别，并具有这个类别上样本的特性。

[0048] 例如，KNN(K-Nearest Neighbor)算法的实现步骤如下：

[0049] 步骤1：初始化距离为最大值；

[0050] 步骤2：计算待测样本和每个训练样本的距离dist；

[0051] 步骤3：得到目前K个最邻近样本中的最大距离maxdist；

[0052] 步骤4：如果dist小于maxdist，则将训练样本作为K-最近样本集合；

[0053] 步骤5：重复步骤2、3、4，直到所有待测样本和训练样本的距离都计算完毕；

[0054] 步骤6：统计K个最邻近样本中每个类别出现的次数；

[0055] 步骤7：选择出现频次最大的类别作为待测样本的类别

[0056] 针对本发明，训练样本选用几个典型的已知种类的时钟源(比如铷钟、铯钟、氢钟、星载钟各选择一个)，用训练样本训练分类器；而待测样本就是我们钟组中的N台钟，KNN算法中的距离就是频偏矩阵，放入分类器进行分类，并赋予不同的分类系数ai＝+1或者ai＝-1根据所述频偏矩阵计算所述时钟源的归一化权值；根据所述分类系数和所述归一化权值进行时钟合成，得到时钟合成信号。

[0057] 根据所述频偏矩阵计算所述时钟源的归一化权值为

[0058]

[0059] 其中μi为第i个时钟源的归一化权值，N为所述时钟源的个数，μi′为第i个时钟源的权值，μi′表示为

[0060]

[0061] 其中fij为第i台时钟源时钟信号和第j台时钟源时钟信号比对计算得到的频偏矩阵，所述时钟合成信号为

[0062]

[0063] 其中，TA(t)为所述时钟合成信号，ai为所述分类系数，μi为第i个时钟源的归一化权值，TAi(t)为第i个时钟源获得的时钟信号。

[0064] 在步骤103中，对所述待测样本进行分类，得到分类系数采用的分类器为KNN分类器，按照每个时钟对系统贡献为正或负对所述待测样本进行分类，产生+1/-1的分类系数。

[0065] 需要说明的是，KNN算法是一种经典的分类算法，核心思想是如果一个样本在特征空间中的k个最相邻(距离最小)的样本中的大多数属于某一个类别，则该样本也属于这个类别，并具有这个类别上样本的特性。

[0066] 本发明使用机器学习分类器和小波变换相结合的方法加权生成高精度时钟基准，时钟源经小波变换后对特定小波分解细节系数进行修正，再将信号进行重构，从而减小时钟源信号中的噪声成分；同时，将时钟源之间两两比对，计算频偏作为加权系数，利用机器学习分类器对时钟进行分类，从而对加权系数做修正，用修正后的加权系数对降噪后的时钟信号进行加权输出，可以有效提高时钟精度和稳定度。

[0067] 图2为一种包含异常信号检测的多基准时钟加权合成方法流程示意图。本实施所述方法包含以下步骤：

[0068] 步骤101，从多个时钟源获得时钟信号。

[0069] 步骤201，对所述时钟信号进行异常信号检测，去除所述时钟信号中的奇异点。

[0070] 在步骤201中，对所述时钟信号进行异常信号检测，去除所述时钟信号中的奇异点，进一步包含：选取所述异常信号检测的测试参考信号，所述测试参考信号为所述时钟信号的任一历史输出值；将所述时钟信号与所述测试参考信号进行比较，得到比较差值；对所述比较差值进行判断，当所述比较差值大于奇异点判断阈值时，所述时钟信号为奇异点。

[0071] 优选地，在步骤201中，对所述时钟信号进行异常信号检测，去除所述时钟信号中的奇异点，进一步包含：将所述时钟信号进行两两比对，并计算得到两两比对差值；根据所述两两比对差值将所述时钟信号分为两类，当所述两两比对差值大于奇异点判断阈值时，判定两个时钟信号为不同类别，当所述两两比对差值不大于所述奇异点判断阈值时，判定两个时钟信号为相同类别；根据所述时钟信号的分类情况，将信号个数少的一类时钟信号判定为奇异点。

[0072] 步骤202，对所述时钟信号进行去除噪声合成，得到去除噪声的时钟信号。

[0073] 在步202中，对所述时钟信号进行去除噪声合成，得到去除噪声的时钟信号，进一步包含：对所述时钟信号进行小波函数分解，得到小波分解细节系数；对所述小波分解细节系数进行硬阈值判断筛选，得到筛选后的小波分解细节系数，所述硬阈值判断筛选为，当小波分解细节系数小于小波分解硬阈值时，将该小波分解细节系数置0，当小波分解细节系数大于或等于所述小波分解硬阈值时，保留该小波分解细节系数；根据所述筛选后的小波分解细节系数，对所述时钟信号进行重构。

[0074] 优选地，对所述时钟信号进行去除噪声合成，得到去除噪声的时钟信号，进一步包含：对所述时钟信号进行小波函数分解，得到小波分解细节系数；对所述小波分解细节系数根据设定的软阈值筛选比例进行软阈值判断筛选，得到筛选后的小波分解细节系数，所述软阈值判断筛选为，当小波分解细节系数小于小波分解硬阈值时，将该小波分解细节系数置0，当小波分解细节系数大于或等于所述小波分解硬阈值时，将该小波分解细节系数减去阈值得到新的小波分解细节系数；根据所述筛选后的小波分解细节系数，对所述时钟信号进行重构。

[0075] 步骤203，对所述去除噪声的时钟信号，进行加权合成，得到时钟合成信号。

[0076] 在步骤203中，对所述去除噪声的时钟信号，进行加权合成，得到时钟合成信号，进一步包含：将所述去除噪声的时钟信号进行两两比对，生成比对向量，并根据所述比对向量计算得到频偏矩阵；根据所述频偏矩阵生成待测样本；对所述待测样本进行分类，得到分类系数；根据所述频偏矩阵计算所述时钟源的归一化权值；根据所述分类系数和所述归一化权值进行时钟合成，得到时钟合成信号。

[0077] 本实施例所实现的频率基准信号频率准确度在n×E-14量级，优于高品质单个铯钟1～2个数量级；本发明可以同时改善时钟源的短期和长期稳定性，即频率准确度达到n×E-14量级的同时，观察窗口100s以下时的频率阿伦方差也比单个时钟源低2～3个数量级。

[0078] 需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。

[0079] 以上所述仅为本发明的实施例而已，并不用于限制本发明。对于本领域技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本发明的权利要求范围之内。