1.一种低轨卫星的星载原子钟向UTC(k)的溯源方法，其特征在于，包括：
S1：低轨卫星从多个地面GNSS观测站获取GNSS导航卫星的精密星历，所述精密星历包括GNSS导航卫星的精密轨道坐标、相对于参考时间UTC(k)的精密钟差；通过GNSS接收模块获取GNSS导航卫星的广播星历和观测数据，所述观测数据包括码伪距和载波相位观测值；
S101，多个地面GNSS观测站观测GNSS导航卫星，计算GNSS导航卫星的精密轨道坐标和相对于参考时间UTC(k)的精密钟差；
S102，多个地面GNSS观测站向低轨卫星发送GNSS导航卫星的精密轨道坐标和相对于参考时间UTC(k)的精密钟差；
S103，在101步骤的同时，低轨卫星的星载GNSS接收模块接收GNSS导航卫星的广播星历、码伪距和载波相位观测值；
S104，低轨卫星的星载数据处理模块接收GNSS导航卫星的广播星历、码伪距和载波相位观测值以及GNSS导航卫星的精密轨道坐标和相对于参考时间UTC(k)的精密钟差；
S2：根据所述GNSS导航卫星的精密轨道坐标、相对于参考时间UTC(k)的精密钟差、所述广播星历以及所述码伪距和载波相位观测值解算低轨卫星的本地时间和参考时间UTC(k)的钟差Ts‑UTC(k)；
S201，对GNSS导航卫星的精密轨道坐标数据、精密钟差数据进行预处理；所述预处理包括误差项处理、插值处理；
S202，计算观测数据从GNSS导航卫星到低轨卫星的GNSS接收模块的传播路径中的各项误差修正项；
S203，将码伪距观测值、广播星历、各项误差修正项带入码伪距观测方程，经过伪距单点定位方法估算初始时刻低轨卫星相对于参考时间UTC(k)的钟差，估算的结果作为后续卡尔曼滤波计算的初始值；将码伪距、精密星历代入码伪距观测方程，将载波相位观测值、精密星历代入载波相位观测方程，联立码伪距观测方程和载波相位观测方程采用卡尔曼滤波方式进行精密单点定位技术的参数估计；对估计得到的参数进行残差检验得到钟差Ts‑UTC(k)；具体为：
由S201和S202步骤得到的低轨卫星的GNSS接收模块数据、广播星历、精密星历和各项误差修正项，形成码伪距观测方程和载波相位观测方程；其中，码伪距观测方程表示为式(1)：
(1)
式中，i代表GNSS导航卫星编号，j代表观测值频段号，代表光速， 代表观测的码伪距值， 代表GNSS导航卫星坐标和接收坐标间的物理距离， 和 分别代表GNSS接收模块参考时间以及GNSS导航卫星时间和卫星导航系统参考时的钟差， 代表GNSS接收模块和GNSS导航卫星硬件时延， 代表步骤S202中所述各项误差修正项，代表GNSS接收模块的码伪距观测噪声；
载波相位观测方程表示为式(2)：
(2)
式中， 代表载波相位观测值， 代表载波波长， 是以周围单位的原始载波观测值，从载波相位观测文件中获得， 代表整数模糊度， 代表频率相关的GNSS接收模块和GNSS导航卫星硬件延迟， 代表载波观测噪声；
先将码伪距观测值和广播星历的卫星坐标和钟差信息带入码伪距观测方程，经过伪距单点定位方法估算初始时刻低轨卫星的坐标和GNSS接收模块相对于参考时间UTC(k)的钟差，估算的结果作为后续卡尔曼滤波计算的初始值；然后将码伪距和载波相位观测值、精密卫星坐标和钟差代入码伪距观测方程和载波相位观测方程，采用卡尔曼滤波的方式进行PPP参数估计；某一时刻观测到n颗卫星，采用双频码伪距和载波相位观测数据，可建立4n个观测方程；误差方程由式(3)表示：
(3)
式中，n为GNSS导航卫星数，m为待估参数量，此处m=2n+4， 为观测值残差向量， 为系数矩阵， 为状态向量， 为观测值和计算值的差；状态向量的待估参数包括低轨卫星的坐标向量 ，本地钟差，第i个GNSS导航卫星的双频模糊度 和 ，状态向量由式(4)表示：
(4)
式中， 、 和 为低轨卫星的坐标向量 的三维坐标， 、 代表第n个GNSS导航卫星的双频模糊度；
由式(1)(2)组合，系数矩阵表示为式(5)：
(5)
式中， 、 、 分别为第1个GNSS导航卫星到低轨卫星连线的方向余弦， 、 、 分别为第n个GNSS导航卫星到低轨卫星连线的方向余弦；对应的系统噪声方差矩阵如式(6)所示：
(6)
式中，Q为系统噪声矩阵， 为低轨卫星的三维坐标的协方差矩阵， 为钟差的协方差矩阵， 和 为双频模糊度的协方差矩阵；
卡尔曼滤波过程表示为式(7)：
(7)
式中， 和 为更新后的状态矩阵和系统噪声方差矩阵，R为观测噪声的方差矩阵，可根据卫星高度角或信噪比随机模型确定，I为单位矩阵， 为矩阵H的转置；经过参数估计后，对计算得到的结果进行残差检验，通过检验的低轨卫星的本地时间和参考时间UTC(k)的钟差Ts‑UTC(k)是本步骤的输出结果；
S3：根据所述低轨卫星的本地时间和参考时间UTC(k)的钟差Ts‑UTC(k)控制星载原子钟，调整星载原子钟输出信号频率，使输出的卫星标准时间Ts与UTC(k)同步，实现星载原子钟到UTC(k)的溯源；
S301，低轨卫星的星钟控制模块对获得的钟差Ts‑UTC(k)进行预报；
S302，低轨卫星的星钟控制模块根据钟差预报值产生频率驾驭量；
S303，低轨卫星的星钟控制模块接收星载原子钟的原始信号；
S304，低轨卫星的星钟控制模块根据频率驾驭量驾驭控制星载原子钟，实现向UTC(k)的溯源并输出卫星标准时间Ts。
2.根据权利要求1所述的一种低轨卫星的星载原子钟向UTC(k)的溯源方法，其特征在于，所述精密轨道坐标的数据间隔为15分钟一组，所述精密钟差的数据为30秒一组。
3.根据权利要求1所述的一种低轨卫星的星载原子钟向UTC(k)的溯源方法，其特征在于，所述插值处理包括：所述GNSS导航卫星的精密轨道坐标数据采用高阶拉格朗日插值，所述精密钟差数据采用线性插值。
4.一种实现权利要求1所述的低轨卫星的星载原子钟向UTC(k)的溯源方法的系统，其特征在于，包括星载原子钟、GNSS接收模块、星钟控制模块和星载数据处理模块；
所述星载原子钟输出10MHz信号分别到星钟控制模块和星载数据处理模块，经过星钟控制模块进行放大分路，一路输入到GNSS接收模块作为参考源，另一路输入到星载数据处理模块用于比对计算；
所述星钟控制模块根据钟差数据调整星载原子钟输出时频信号的频率和相位，并输出调整后的时频信号；
所述GNSS接收模块通过天线接收GNSS导航卫星的广播星历和观测数据并输入到星载数据处理模块；同时接收GNSS观测站发送的精密星历并输入到星载数据处理模块；
所述星载数据处理模块根据输入的广播星历、观测数据和精密星历，计算得到本地频率源和参考频率源的钟差Ts‑UTC(k)；根据计算的钟差Ts‑UTC(k)、星载原子钟输入的信号以及星钟控制模块调整后的时频信号，计算星钟控制量并输入到星钟控制模块，实现星钟向UTC(k)的溯源。
5.根据权利要求4所述的一种低轨卫星的星载原子钟向UTC(k)的溯源系统，其特征在于，还包括上注通信链路，所述上注通信链路用于地面GNSS观测站与低轨卫星之间的数据通信。