CN105337681B - 一种时钟同步规划中的prc接入网元选择方法 - Google Patents

Abstract

本发明公开了一种时钟同步规划中的PRC接入网元选择方法，首先根据网络中网元个数为参考，确定网络中所需规划的PRC个数；其次对所有网元按照其连接的网元个数多少进行降序排列，选择出备选网元；然后，根据所有网元到备选网元的时钟跟踪长度，按照约束的定量参数条件，即备选网元在有效跳数内覆盖网元个数尽可能多、覆盖网元的平均时钟跟踪长度尽可能小、PRC优先接入网元在网络中部署尽可能均匀这三条原则，通过迭代优化法来选择并确定PRC优先接入网元。对于大规模和复杂结构网络，能够通过计算机自动计算代替人工规划PRC接入网元，可有效降低时钟同步规划工作的困难性，提高规划效率，可保证时钟规划方案的整体性能。

Description

一种时钟同步规划中的PRC接入网元选择方法

技术领域

本发明涉及一种PRC接入网元的选择，特别是涉及一种时钟同步规划中的PRC接入网元选择方法，属于光通信技术领域。

背景技术

随着通信技术和业务的不断发展，光通信网的规模在不断扩大，光网络设备之间的信息交互对于时钟同步和时间性能的要求也越来越高。

同步网络结构上可分为准同步网络、主从同步网络以及混合同步网络三种。

准同步网络中所有时钟设备(Building Integrated Timing Supply,BITS)独立运行，不依赖于其他时钟，每个主参考时钟(Primary Reference Clock,PRC)负责同步各自区域内设备，因为主参考时钟PRC的精度较高，所以不同区域间设备亦可实现高精度的同步，从而达到整个网络的同步。

主从同步网络是一种分级结构的同步网，网络中的时钟设备分为多级，各级时钟由于配置的差异而具有不同的稳定性和准确度，下级设备严格同步上级设备，具有严格等级性。

混合同步网络指准同步和主从同步综合应用的组网方式，是最常用的组网方式。

对于光通信网络的时钟同步规划，一般会做以下技术参数的约束：2台时钟设备之间网元数量不超过N₁个，整个同步链不能超过N₂个网元，整个同步链不能超过N₃台时钟设备，同步时钟链路不可以出现环路。此外，还要综合权衡方案的技术性能与经济性。

考虑到时钟同步规划设计应该具有高可靠、高性能、易维护、可扩展等特点，当网络规模较大，网元节点个数过多时，为了提高同步网络的可靠性，进一步减少故障风险，可以通过划分时钟同步区的方法来优化设计。当前，时钟同步规划的一般步骤为：首先，获取客户需求，包括时钟数量要求，部署节点限制等；其次，根据网络拓扑结构及规模，确定时钟同步区数量及范围；再次，遵照规划规则，在同步区内进行时钟部署及时钟流规划；最后，检查规划时钟流是否存在环路，并评估时钟同步规划方案的性能。

目前，对于时钟同步规划，大多还依赖于人工规划，随着网络规模的扩大以及结构的复杂性提高，人工规划越来越困难，时钟规划方案正确性、可靠性很难控制。在确定时钟设备接入节点后，目前已有最小生成树等算法被用来规划时钟流链路，但是大规模和复杂结构网络下，同步区划分以及PRC时钟部署的仍主要依赖人工，无法实现智能规划。

发明内容

本发明的主要目的在于，克服现有技术中的不足，提供一种时钟同步规划中的PRC接入网元选择方法，特别适用于大规模和复杂结构网络的时钟同步规划。

本发明所要解决的技术问题是能够实现大规模和复杂结构网络中的PRC时钟接入网元的自动选择，可用到网络时钟同步规划算法中，可克服人工规划的复杂性及性能评估的不确定性等问题。

为了达到上述目的，本发明所采用的技术方案是：

一种时钟同步规划中的PRC接入网元选择方法，包括以下步骤：

1)输入网元信息参数，所述网元信息包括网元个数N、网元类型、网元ID以及网元连接关系；

将每个网元记为NE_n，其中n＝1,2,…,N，N为自然数；

将网元连接关系通过网元连接对＜NE_i,NE_j＞来表征，其中i,j＝1,2,…,N；

2)根据网络中输入的网元个数，通过式(1)计算所需限定的主参考时钟PRC个数Nprc_min，

其中，N为网元个数，N₂为网络中需要约束的同步链所要限定的网元个数；

并确定该网络中所需规划的主参考时钟PRC个数Nprc_se，t要求Nprc_set≥Nprc_min；

3)遍历网络中所有网元，对每个网元所直接连接的网元个数进行统计，形成网元连接个数序列M，记为

M＝{M₁,M₂,…,M_i}，i＝1,2,…,N

其中，M_i表示与第i个网元直接连接的网元个数；

4)主参考时钟PRC接入备选网元集的选择；

对形成的网元连接个数序列M，根据每个网元直接连接的网元个数进行降序排列；

取排位于前Y个网元作为备选网元，取Y≥2*Nprc_set；

将排位于前Y个网元汇集成主参考时钟PRC接入备选网元集，记为

{NE_PRC_y},y＝1,2,…Y；

5)遍历Y个备选网元，计算出网络中所有网元对该Y个备选网元的时钟跟踪长度，并形成时钟跟踪长度矩阵D，记为

其中，d_i,j表示第i个网元到第j个备选网元的时钟跟踪长度值；

6)根据时钟跟踪长度矩阵D、并按照约束的定量参数条件，通过迭代优化法来选择主参考时钟PRC优先接入网元；

所述约束的定量参数条件包括备选网元在有效跳数内覆盖网元个数需取上限值，覆盖网元的平均时钟跟踪长度需取下限值，主参考时钟PRC优先接入网元在网络中部署需分布均匀；

所述主参考时钟PRC优先接入网元为Nprc个，是从主参考时钟PRC接入备选网元集中选择并确定的。

本发明进一步设置为：所述步骤6)中通过迭代优化法来选择主参考时钟PRC优先接入网元，具体步骤为，

6-1)初始化备选网元覆盖程度的权重矩阵；

首先，初始化备选网元在有效跳数内覆盖网元标识矩阵N_c；

建立一个N×Y维的全零新矩阵作为备选网元在有效跳数内覆盖网元标识矩阵N_c，令N_c(find(D≤N₁))＝1，即与时钟跟踪长度矩阵D中所有小于等于N₁的元素坐标相对应的备选网元在有效跳数内覆盖网元标识矩阵N_c的元素值均设为1；其中，find(D≤N₁)表示时钟跟踪长度矩阵D中满足d_i,j≤N₁的元素所对应的坐标(i,j)，N₁表示2台时钟设备之间所要限定的网元数量；

其次，初始化覆盖网元的有效时钟追踪距离矩阵D_c＝D·N_c，即为时钟跟踪长度矩阵D与备选网元在有效跳数内覆盖网元标识矩阵N_c的内积，获得覆盖网元的有效时钟追踪距离矩阵D_c中的元素值为有效跳数N₁内每个网元到备选网元的时钟追踪距离；

于是，根据备选网元在有效跳数内覆盖网元标识矩阵N_c和网络中所需规划的主参考时钟PRC个数Nprc_set，初始化备选网元覆盖程度的权重矩阵N_c′＝N_c·Nprc_set；

6-2)对备选网元，计算其在有效跳数N₁内备选网元覆盖程度的加权系数；

将备选网元覆盖程度的权重矩阵N_c′中每列元素相加，得到一个相加序列n_c，记为

n_c＝{n₁,n₂,…,n_i}，i＝1,2,…,Y

其中，n_i表示第i个备选网元在有效跳数N₁内备选网元覆盖程度的加权系数；

6-3)根据备选网元覆盖程度的加权系数值来判断并确定主参考时钟PRC优先接入网元；

将相加序列n_c中备选网元覆盖程度的加权系数的上限值记为加权系数顶值；

若相加序列n_c中加权系数顶值有且仅有一个、即唯一，则令主参考时钟PRC接入备选网元集{NE_PRC_y}中的y等于加权系数顶值所对应的备选网元序号，确定第y个备选网元作为主参考时钟PRC优先接入网元，并令已选主参考时钟PRC优先接入网元个数作+1累算、即N_prc＝N_prc+1，后跳转到步骤6-6)；

若相加序列n_c中加权系数顶值至少为两个、即存在有相同值，则进入步骤6-4)；

6-4)根据覆盖网元的平均时钟跟踪长度来判断并确定主参考时钟PRC优先接入网元；

对于相加序列n_c中加权系数顶值所对应的备选网元序号，计算备选网元序号对应的备选网元所覆盖网元的平均时钟跟踪长度；将相加序列n_c中加权系数顶值在覆盖网元的有效时钟追踪距离矩阵D_c中所对应的列元素相加并除以加权系数顶值，即为加权系数顶值对应的备选网元所覆盖网元的平均时钟跟踪长度；

按照约束的定量参数条件，将备选网元所覆盖网元的平均时钟跟踪长度需取下限值，备选网元设为共有S个；

若S+N_prc≤Nprc_set，则将对应的S个备选网元全部作为主参考时钟PRC优先接入网元，并令已选的主参考时钟PRC优先接入网元个数作+S累算、即N_prc＝N_prc+S，后跳转到步骤6-6)；

若S+N_prc＞Nprc_set，则进入步骤6-5)；

6-5)计算S个备选网元之间的相互距离，根据备选网元到其他备选网元距离之和降序排列，取前Nprc_set-N_prc个备选网元作为主参考时钟PRC优先接入网元，并令N_prc＝Nprc_set；

6-6)若N_prc＜Nprc_set，则更新备选网元覆盖程度的权重矩阵N_c′，再对备选网元覆盖程度的权重矩阵N_c′中的每一列减去已确定的备选网元序号所对应的备选网元在有效跳数内覆盖网元标识矩阵N_c中列的值，从而得到新标识矩阵N_c′；

令新标识矩阵N_c′中小于0的元素值为0，并令新标识矩阵N_c′和备选网元在有效跳数内覆盖网元标识矩阵N_c中与本次迭代已确定的备选网元序号所对应的列的值均为0，重复步骤6-2)～步骤6-6)；

直至N_prc＝Nprc_set，则迭代结束。

与现有技术相比，本发明具有的有益效果是：

本发明首先根据网络中网元个数为参考，确定网络中所需规划的PRC个数Nprc_set；其次对所有网元按照其连接的网元个数多少进行降序排列，选择前Y(Y≥2*Nprc_set)个网元作为备选网元；然后，根据所有网元到备选网元的时钟跟踪长度，按照约束的定量参数条件，即备选网元在有效跳数内覆盖网元个数尽可能多、覆盖网元的平均时钟跟踪长度尽可能小、PRC优先接入网元在网络中部署尽可能均匀这三条原则，通过迭代优化法来选择并确定Nprc_set个备选网元作为PRC优先接入网元。本发明给出的算法和其中采用迭代寻优的思路，对于大规模和复杂结构网络，能够通过计算机自动计算代替人工规划PRC接入网元，可有效的降低时钟同步规划工作的困难性，并提高规划效率；而且能够较常规方法更好地保证时钟规划方案的整体性能，可克服人工规划的复杂性及性能评估的不确定性。

上述内容仅是本发明技术方案的概述，为了更清楚的了解本发明的技术手段，下面结合附图对本发明作进一步的描述。

附图说明

图1为本发明一种时钟同步规划中的PRC接入网元选择方法的流程图；

图2为本发明实施例中网络拓扑实例示意图；

图3为本发明实施例中网络节点时钟同步规划实例示意图。

具体实施方式

下面结合说明书附图，对本发明作进一步的说明。

如图1所示，本发明提供一种时钟同步规划中的PRC接入网元选择方法，包括以下步骤：

1)输入网元信息参数，所述网元信息包括网元个数N、网元类型、网元ID以及网元连接关系；将每个网元记为NE_n，其中n＝1,2,…,N，N为自然数；将网元连接关系通过网元连接对＜NE_i,NE_j＞来表征，其中i,j＝1,2,…,N。

2)根据网络中输入的网元个数，通过式(1)计算所需限定的主参考时钟PRC个数Nprc_min，

其中，N为网元个数，N₂为网络中需要约束的同步链所要限定的网元个数；并确定该网络中所需规划的主参考时钟PRC个数Nprc_set，要求Nprc_set≥Nprc_min。

3)遍历网络中所有网元，对每个网元所直接连接的网元个数进行统计，形成网元连接个数序列M，记为

M＝{M₁,M₂,…,M_i}，i＝1,2,…,N

其中，M_i表示与第i个网元直接连接的网元个数。

4)主参考时钟PRC接入备选网元集的选择；

对形成的网元连接个数序列M，根据每个网元直接连接的网元个数进行降序排列；取排位于前Y个网元作为备选网元，取Y≥2*Nprc_set；将排位于前Y个网元汇集成主参考时钟PRC接入备选网元集，记为{NE_PRC_y},y＝1,2,…Y。

5)遍历Y个备选网元，计算出网络中所有网元对该Y个备选网元的时钟跟踪长度，并形成时钟跟踪长度矩阵D，记为

其中，d_i,j表示第i个网元到第j个备选网元的时钟跟踪长度值。

6)根据时钟跟踪长度矩阵D，并考虑主时钟部署网元要尽可能在有效跳数N₁内覆盖尽量多的网元、覆盖网元的平均时钟跟踪长度尽可能小、同时要尽量满足主参考时钟PRC优先接入网元在整个网络中分布的均匀性；通过如下迭代优化法来选择主参考时钟PRC优先接入网元，从主参考时钟PRC接入备选网元集中选择并确定的Nprc个备选网元作为主参考时钟PRC优先接入网元。

6-1)初始化备选网元覆盖程度的权重矩阵；

首先，初始化备选网元在有效跳数内覆盖网元标识矩阵N_c；

建立一个N×Y维的全零新矩阵作为备选网元在有效跳数内覆盖网元标识矩阵N_c，令N_c(find(D≤N₁))＝1，即与时钟跟踪长度矩阵D中所有小于等于N₁的元素坐标相对应的备选网元在有效跳数内覆盖网元标识矩阵N_c的元素值均设为1；其中，find(D≤N₁)表示时钟跟踪长度矩阵D中满足d_i,j≤N₁的元素所对应的坐标(i,j)，N₁表示2台时钟设备之间所要限定的网元数量；

其次，初始化覆盖网元的有效时钟追踪距离矩阵D_c＝D·N_c，即为时钟跟踪长度矩阵D与备选网元在有效跳数内覆盖网元标识矩阵N_c的内积，获得覆盖网元的有效时钟追踪距离矩阵D_c中的元素值为有效跳数N₁内每个网元到备选网元的时钟追踪距离；

于是，根据备选网元在有效跳数内覆盖网元标识矩阵N_c和网络中所需规划的主参考时钟PRC个数Nprc_set，初始化备选网元覆盖程度的权重矩阵N_c′＝N_c·Nprc_set。

6-2)对备选网元，计算其在有效跳数N₁内备选网元覆盖程度的加权系数；

将备选网元覆盖程度的权重矩阵N_c′中每列元素相加，得到一个相加序列n_c，记为

n_c＝{n₁,n₂,…,n_i}，i＝1,2,…,Y

其中，n_i表示第i个备选网元在有效跳数N₁内备选网元覆盖程度的加权系数。

6-3)根据备选网元覆盖程度的加权系数值来判断并确定主参考时钟PRC优先接入网元；

将相加序列n_c中备选网元覆盖程度的加权系数的上限值记为加权系数顶值；

若相加序列n_c中加权系数顶值有且仅有一个、即唯一，则令主参考时钟PRC接入备选网元集{NE_PRC_y}中的y等于加权系数顶值所对应的备选网元序号，确定第y个备选网元作为主参考时钟PRC优先接入网元，并令已选主参考时钟PRC优先接入网元个数作+1累算、即N_prc＝N_prc+1，后跳转到步骤6-6)；

若相加序列n_c中加权系数顶值至少为两个、即存在有相同值，则进入步骤6-4)。

6-4)根据覆盖网元的平均时钟跟踪长度来判断并确定主参考时钟PRC优先接入网元；

对于相加序列n_c中加权系数顶值所对应的备选网元序号，计算备选网元序号对应的备选网元所覆盖网元的平均时钟跟踪长度；将相加序列n_c中加权系数顶值在覆盖网元的有效时钟追踪距离矩阵D_c中所对应的列元素相加并除以加权系数顶值，即为加权系数顶值对应的备选网元所覆盖网元的平均时钟跟踪长度；

按照约束的定量参数条件，将备选网元所覆盖网元的平均时钟跟踪长度需取下限值，备选网元设为共有S个；判断S+N_prc是否大于规定个数；

若S+N_prc≤Nprc_set，则将对应的S个备选网元全部作为主参考时钟PRC优先接入网元，并令已选的主参考时钟PRC优先接入网元个数作+S累算、即N_prc＝N_prc+S，后跳转到步骤6-6)；

若S+N_prc＞Nprc_set，则进入步骤6-5)。

6-5)计算S个备选网元之间的相互距离，根据备选网元到其他备选网元距离之和降序排列，取前Nprc_set-N_prc个备选网元作为主参考时钟PRC优先接入网元，并令N_prc＝Nprc_set。

6-6)若N_prc＜Nprc_set，则更新备选网元覆盖程度的权重矩阵N_c′，再对备选网元覆盖程度的权重矩阵N_c′中的每一列减去已确定的备选网元序号所对应的备选网元在有效跳数内覆盖网元标识矩阵N_c中列的值，从而得到新标识矩阵N_c′；

令新标识矩阵N_c′中小于0的元素值为0，并令新标识矩阵N_c′和备选网元在有效跳数内覆盖网元标识矩阵N_c中与本次迭代已确定的备选网元序号所对应的列的值均为0，重复步骤6-2)～步骤6-6)；

直至N_prc＝Nprc_set，则迭代结束。

下面通过一个网络实施例的PRC接入网元规划来说明本发明算法的可行性，为了更好的评估PRC接入网元选择对于时钟规划性能的影响，给出以下量化参数定义。

网元时钟跟踪长度NTTL(全称为NE TIMER TRACK LENGTH)：网元时钟在时钟链中相对于基准时钟而言被SDH网元传递的次数。时钟每经过一个网元的传递，该网元的时钟跟踪长度就相对上游的网元值加1；第一个接入外部时钟的网元的时钟跟踪长度为0。

全网平均时钟跟踪长度NATTL(全称为Ne Average Time Track Length)：全网的网元时钟跟踪长度的平均值，计算公式为式(2)

其中，NTTL(NE_i)表示第i个网元NE_i的网元时钟跟踪长度。

全网最小平均时钟跟踪长度NATTL_min(全称为Min-Ne Average Time TrackLength)：在时钟接入点确定的情况下，各网元不同的时钟跟踪路径，导致可得到不同的全网平均跟踪长度。全网平均时钟跟踪长度NATTL中的最小值即是全网最小平均时钟跟踪长度，即为

NATTL_min＝min(NATTL)

如图2所示，构建了一个实例网络，图2中方框表示网元，方框内数字表示网元编号，无箭头直线表示2个网元间存在光纤链路连接、即网元间有可传输时钟信号的链路。

假设光通信网络的时钟同步规划中以下约束为条件：2台时钟设备之间网元数量不超过N₁(N₁＝10)个，整个同步链不能超过N₂(N₂＝60)个网元，同步时钟链路不可以出现环路。根据本发明的一种时钟同步规划中的PRC接入网元选择方法，针对实例网络进行PRC接入网元选取。

由图2可知，网络中网元总数为N＝32个，网元标记为NE_n(n＝1,2,…,32)；根据网络中网元个数，计算所需限定的主参考时钟PRC个数Nprc_min＝2，进一步确定该网络所需规划的主参考时钟PRC个数为Nprc_set＝2。根据每个网元的连接网元个数情况取前Y(Y≥4)个网元作为主参考时钟PRC接入备选网元集，主参考时钟PRC备选网元集如下，包含了连接网元个数为4和3的所有网元共14个，即为

{NE_PRC_y}_Y＝{NE₃,NE₄,NE₅,NE₆,NE₇,NE₈,NE₉,NE₁₀,NE₁₁,NE₁₂,NE₁₃,NE₁₈,NE₂₃,NE₂₈}

按照本发明的迭代优化法来选择主参考时钟PRC优先接入网元，最终确定2个主参考时钟PRC优先接入网元为：{NE₃,NE₄}。

如图3所示，给出了本发明实施例中确定PRC接入网元后规划网络的时钟流链路示意图。图3中，圆圈表示PRC，带箭头的实线表示优先级为1级的时钟流链路，带箭头的虚线表示优先级为2级的时钟流链路，箭头方向为时钟流方向。

对于主参考时钟PRC接入备选网元集{NE_PRC_y}_Y中的所有备选网元，分别计算出该网络的全网最小平均时钟跟踪长度(NATTL_min)，如下表1所示。

网元	NE3	NE4	NE5	NE6	NE7	NE8	NE9
								NATTLmin	3.6	3.6	3.7	3.7	3.9	3.9	3.7
网元	NE10	NE11	NE12	NE13	NE18	NE23	NE28
								NATTLmin	3.7	3.9	3.9	4.5	4.5	4.5	4.5

表1

从时钟性能角度定量分析，本发明的方法计算得到的PRC优先接入网元具有最小的全网最小平均时钟跟踪长度，对于大规模和复杂结构网络，能够通过计算机自动计算代替人工规划PRC接入网元，可有效的降低时钟同步规划工作的困难性，并提高规划效率，同时能够更好地保证时钟规划方案的性能。

以上显示和描述了本发明的基本原理、主要特征及优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。

Claims (2)

1.一种时钟同步规划中的PRC接入网元选择方法，其特征在于，包括以下步骤：

1)输入网元信息参数，所述网元信息包括网元个数N、网元类型、网元ID以及网元连接关系；

将每个网元记为NE_n，其中n＝1,2,…,N，N为自然数；

将网元连接关系通过网元连接对＜NE_i,NE_j＞来表征，其中i,j＝1,2,…,N；

2)根据网络中输入的网元个数，通过式(1)计算所需限定的主参考时钟PRC个数Nprc_min，

其中，N为网元个数，N₂为网络中需要约束的同步链所要限定的网元个数；

并确定该网络中所需规划的主参考时钟PRC个数Nprc_se，t要求Nprc_set≥Nprc_min；

3)遍历网络中所有网元，对每个网元所直接连接的网元个数进行统计，形成网元连接个数序列M，记为

M＝{M₁,M₂,…,M_i}，i＝1,2,…,N

其中，M_i表示与第i个网元直接连接的网元个数；

4)主参考时钟PRC接入备选网元集的选择；

对形成的网元连接个数序列M，根据每个网元直接连接的网元个数进行降序排列；

取排位于前Y个网元作为备选网元，取Y≥2*Nprc_set；

将排位于前Y个网元汇集成主参考时钟PRC接入备选网元集，记为{NE_PRC_y},y＝1,2,…Y；

5)遍历Y个备选网元，计算出网络中所有网元对该Y个备选网元的时钟跟踪长度，并形成时钟跟踪长度矩阵D，记为

其中，d_i,j表示第i个网元到第j个备选网元的时钟跟踪长度值；

6)根据时钟跟踪长度矩阵D、并按照约束的定量参数条件，通过迭代优化法来选择主参考时钟PRC优先接入网元；

所述约束的定量参数条件包括备选网元在有效跳数内覆盖网元个数需取上限值，覆盖网元的平均时钟跟踪长度需取下限值，主参考时钟PRC优先接入网元在网络中部署需分布均匀；

所述主参考时钟PRC优先接入网元为Nprc个，是从主参考时钟PRC接入备选网元集中选择并确定的。

2.根据权利要求1所述的一种时钟同步规划中的PRC接入网元选择方法，其特征在于：所述步骤6)中通过迭代优化法来选择主参考时钟PRC优先接入网元，具体步骤为，

6-1)初始化备选网元覆盖程度的权重矩阵；

首先，初始化备选网元在有效跳数内覆盖网元标识矩阵N_c；

建立一个N×Y维的全零新矩阵作为备选网元在有效跳数内覆盖网元标识矩阵N_c，令N_c(find(D≤N₁))＝1，即与时钟跟踪长度矩阵D中所有小于等于N₁的元素坐标相对应的备选网元在有效跳数内覆盖网元标识矩阵N_c的元素值均设为1；其中，find(D≤N₁)表示时钟跟踪长度矩阵D中满足d_i,j≤N₁的元素所对应的坐标(i,j)，N₁表示2台时钟设备之间所要限定的网元数量；

其次，初始化覆盖网元的有效时钟追踪距离矩阵D_c＝D·N_c，即为时钟跟踪长度矩阵D与备选网元在有效跳数内覆盖网元标识矩阵N_c的内积，获得覆盖网元的有效时钟追踪距离矩阵D_c中的元素值为有效跳数N₁内每个网元到备选网元的时钟追踪距离；

于是，根据备选网元在有效跳数内覆盖网元标识矩阵N_c和网络中所需规划的主参考时钟PRC个数Nprc_set，初始化备选网元覆盖程度的权重矩阵N_c′＝N_c·Nprc_set；

6-2)对备选网元，计算其在有效跳数N₁内备选网元覆盖程度的加权系数；

将备选网元覆盖程度的权重矩阵N_c′中每列元素相加，得到一个相加序列n_c，记为

n_c＝{n₁,n₂,…,n_i}，i＝1,2,…,Y

其中，n_i表示第i个备选网元在有效跳数N₁内备选网元覆盖程度的加权系数；

6-3)根据备选网元覆盖程度的加权系数值来判断并确定主参考时钟PRC优先接入网元；

将相加序列n_c中备选网元覆盖程度的加权系数的上限值记为加权系数顶值；

若相加序列n_c中加权系数顶值有且仅有一个、即唯一，则令主参考时钟PRC接入备选网元集{NE_PRC_y}中的y等于加权系数顶值所对应的备选网元序号，确定第y个备选网元作为主参考时钟PRC优先接入网元，并令已选主参考时钟PRC优先接入网元个数作+1累算、即N_prc＝N_prc+1，后跳转到步骤6-6)；

若相加序列n_c中加权系数顶值至少为两个、即存在有相同值，则进入步骤6-4)；

6-4)根据覆盖网元的平均时钟跟踪长度来判断并确定主参考时钟PRC优先接入网元；

对于相加序列n_c中加权系数顶值所对应的备选网元序号，计算备选网元序号对应的备选网元所覆盖网元的平均时钟跟踪长度；将相加序列n_c中加权系数顶值在覆盖网元的有效时钟追踪距离矩阵D_c中所对应的列元素相加并除以加权系数顶值，即为加权系数顶值对应的备选网元所覆盖网元的平均时钟跟踪长度；

按照约束的定量参数条件，将备选网元所覆盖网元的平均时钟跟踪长度需取下限值，备选网元设为共有S个；

若S+N_prc≤Nprc_set，则将对应的S个备选网元全部作为主参考时钟PRC优先接入网元，并令已选的主参考时钟PRC优先接入网元个数作+S累算、即N_prc＝N_prc+S，后跳转到步骤6-6)；

若S+N_prc＞Nprc_set，则进入步骤6-5)；

6-5)计算S个备选网元之间的相互距离，根据备选网元到其他备选网元距离之和降序排列，取前Nprc_set-N_prc个备选网元作为主参考时钟PRC优先接入网元，并令N_prc＝Nprc_set；

6-6)若N_prc＜Nprc_set，则更新备选网元覆盖程度的权重矩阵N_c′，再对备选网元覆盖程度的权重矩阵N_c′中的每一列减去已确定的备选网元序号所对应的备选网元在有效跳数内覆盖网元标识矩阵N_c中列的值，从而得到新标识矩阵N_c′；

令新标识矩阵N_c′中小于0的元素值为0，并令新标识矩阵N_c′和备选网元在有效跳数内覆盖网元标识矩阵N_c中与本次迭代已确定的备选网元序号所对应的列的值均为0，重复步骤6-2)～步骤6-6)；

直至N_prc＝Nprc_set，则迭代结束。