1.一种锁相环百皮秒级时钟调控方法，其特征在于，所述方法包括：
根据振荡器时钟的历元间接收机钟差得到：
式中，“[·]”表示四舍五入取整，y(i)表示历元i的时钟调控量，tepoch表示相邻历元之间的时间间隔，z为压控晶体振荡器的最小数字调节量，初始化期间，DPT时间传递估计的接收机钟差表示为：
dtr(i)＝ttran(i)+dt′r(i)                   (2)
其中，dtr(i)是初始化期间估计的接收机钟差，ttran(i)是DPT初始化期间的误差，dt′r(i)是接收机钟差；
在第一次估计接收机钟差之后，接收机晶体振荡器通过频率控制进行调整，在振荡过程中增加一个预期的时间长度T0，通过设置频率补偿的预期时间长度T0来获得接收机的预期钟差变化量，计算如式(3)所示：
对应的时钟调控量为：
在本地晶体振荡器实时控制的过程中，对每个历元进行迭代计算，以确保DPT收敛后的晶振时钟利用式(4)锁定时频参考；
获取参考时间与被控晶振的频率偏差，利用数模转换器将数字量频率偏差转换成对应的电压调节值，并利用所述电压调节值对晶振进行实时调整，频率调节范围为：
fA＝f0±Δf                        (5)
其中f0为OCXO晶振的标称频率值；±Δf为OCXO晶振的实际工作频率与标称频率偏差的可调节的上下界；
估计最小频率偏差的调节值为：
式中，Vmax为可调节的最大电压值，Vmin为可调节的最小电压值，V为可调节电压区间，fmax为可调节的最大频率值，fmin为可调节的最小频率值，通过式(6)可实现频率偏差向电压调节量的转换。
2.根据权利要求1所述的方法，其特征在于，所述方法还包括：
通过如下公式计算稳态相差：
其中 为稳态相差，E(s)为误差传递函数，θi(s)为输入时钟信号；
二阶环由一个积分器和一个常数项组成，即
式中，τ1表示积分器的系数，τ2表示常数项的系数；
环路传递函数H(s)和误差传递函数E(s)分别表示为：
式中，K1表示二阶锁相环电路中的数模转换器的系数，K2表示二阶锁相环电路中的ocxo的系数，s表示自动控制原理里的拉普拉斯变换的变量；
将式(8)和(9)中分母写成归一表达，得到：
式中ωn表示环路的特征频率，ζ表示阻尼系数；通过式(10)，环路的传递函数H(s)和误差传递函数E(s)表示为：
根据各项噪声的叠加，二阶锁相环的输出噪声的功率谱密度So(f)表示为：
So(f)＝Sref(f)+SDAC(f)+SOCXO(f)                                  (12)其中，Sref(f)和SOCXO(f)分别表示为精密时间传递噪声功率谱密度和振荡器自由振荡功率谱密度，SDAC(f)表示为数模转换模块引入的相位抖动噪声功率谱密度，通过数模转换量化噪声建模，实现二阶锁相环路参数的确定；
根据对输入噪声的抑制效果确定阻尼系数ζ；
根据总输出相位方差确定特征频率ωn。
3.根据权利要求2所述的方法，其特征在于，确定的阻尼系数ζ＝0.707。
4.根据权利要求2所述的方法，其特征在于，根据总输出相位方差确定特征频率ωn，包括：
基于系统总输出相位方差最小原则，确定二阶锁相环的等效噪声带宽Bn为：
将带宽Bn取为交点频率fc，通过式(14)获得特征频率ωn：
5.一种锁相环百皮秒级时钟调控装置，其特征在于，所述装置包括时钟精密调控模块，所述时钟精密调控模块被配置为：
根据振荡器时钟的历元间接收机钟差得到：
式中，“[·]”表示四舍五入取整，y(i)表示历元i的时钟调控量，tepoch表示相邻历元之间的时间间隔，z为压控晶体振荡器的最小数字调节量，初始化期间，DPT时间传递估计的接收机钟差表示为：
dtr(i)＝ttran (i)+dtr′(i)                                           (2)其中，dtr(i)是初始化期间估计的接收机钟差，ttran(i)是DPT初始化期间的误差，dtr′(i)是接收机钟差；
在第一次估计接收机钟差之后，接收机晶体振荡器通过频率控制进行调整，在振荡过程中增加一个预期的时间长度T0，通过设置频率补偿的预期时间长度T0来获得接收机的预期钟差变化量，计算如式(3)所示：
对应的时钟调控量为：
在本地晶体振荡器实时控制的过程中，对每个历元进行迭代计算，以确保DPT收敛后的晶振时钟利用式(4)锁定时频参考；
获取参考时间与被控晶振的频率偏差，利用数模转换器将数字量频率偏差转换成对应的电压调节值，并利用所述电压调节值对晶振进行实时调整，频率调节范围为：
fA＝f0±Δf                                                       (5)其中f0为OCXO晶振的标称频率值；±Δf为OCXO晶振的实际工作频率与标称频率偏差的可调节的上下界；
估计最小频率偏差的调节值为：
式中，Vmax为可调节的最大电压值，Vmin为可调节的最小电压值，V为可调节电压区间，fmax为可调节的最大频率值，fmin为可调节的最小频率值，通过式(6)可实现频率偏差向电压调节量的转换。
6.根据权利要求5所述的装置，其特征在于，所述装置还包括晶振时钟调控模块，所述晶振时钟调控模块被配置为：
通过如下公式计算稳态相差：
其中 为稳态相差，E(s)为误差传递函数，θi(s)为输入信号；
二阶环由一个积分器和一个常数项组成，即
式中，τ1表示积分器的系数，τ2表示常数项的系数；
环路传递函数H(s)和误差传递函数E(s)分别表示为：
式中，K1表示二阶锁相环电路中的数模转换器的系数，K2表示二阶锁相环电路中的ocxo的系数，s表示自动控制原理里的拉普拉斯变换的变量；
将式(8)和(9)中分母写成归一表达，得到：
式中ωn表示环路的特征频率，ζ表示阻尼系数；通过式(10)，环路的传递函数H(s)和误差传递函数E(s)表示为：
根据各项噪声的叠加，二阶锁相环的输出噪声的功率谱密度So(f)表示为：
So(f)＝Sref(f)+SDAC(f)+SOCXO(f)                                  (12)其中，Sref(f)和SOCXO(f)分别表示为精密时间传递噪声功率谱密度和振荡器自由振荡功率谱密度，SDAC(f)表示为数模转换模块引入的相位抖动噪声功率谱密度，通过数模转换量化噪声建模，实现二阶锁相环路参数的确定；
根据对输入噪声的抑制效果确定阻尼系数ζ；
根据总输出相位方差确定特征频率ωn。
7.根据权利要求6所述的装置，其特征在于，确定的阻尼系数ζ＝0.707。
8.根据权利要求6所述的装置，其特征在于，所述晶振时钟调控模块被进一步配置为：
基于系统总输出相位方差最小原则，确定二阶锁相环的等效噪声带宽Bn为：
将带宽Bn取为交点频率fc，通过式(14)获得特征频率ωn：
9.一种锁相环百皮秒级时钟调控系统，其特征在于：所述系统包括：
存储器，用于存储计算机程序；
处理器，用于执行所述计算机程序以实现如权利要求1至4中任一项所述的方法。
10.一种存储有指令的非暂时性计算机可读存储介质，当所述指令由处理器执行时，执行根据权利要求1至4中任一项所述的方法。