[发明专利]用于时钟驯服的频率偏移量分析方法

权利要求书

1.用于时钟驯服的频率偏移量分析方法，其特征在于，包括：

S1、获取当前观测数据和历史观测数据；

S2、修正历史观测数据，具体包括：

S21、假设当前的本地时钟频率控制值为C₀、时间为t₀，相位为ph₀，之前一秒的数据为C₁、t₁、ph₁，……，数据n的时钟频率控制值为C_n、时间为t_n、相位为ph_n，本地时钟的频率值feq与执行的频率控制值C之间的对应关系为feq=f(C)；

S22、数据n的相位修正值计算方法为 ;

S23、数据n的相位经过修正后的值为；

S3、处理当前观测数据和修正后历史观测数据的异常数据；具体包括：

S31、用拉依达准则和接收机报文内容判断两种方法共同来筛选异常数据，当拉依达准则和接收机报文内容任一判断确定为异常数据时，对异常数据做异常丢弃标记；

S32、从异常丢弃标记序列开头Mark₀开始，依次向后寻找第一个未被标记异常丢弃的标志位，作为Mark_i；从Mark_i开始，依次向后寻找第一个被标记异常丢弃的标志位Mark_m+1；

S33、判断是否寻找到Mark_m+1，若是，则令m+1=m，并返回S32；若否，则Mark_m=Mark_n；Mark_i到Mark_m的数据是一段未被“断点”分割的连续观测数据，并进入S34；

S34、重复S33与S34，直到当前观测数据中数据Mark₀、Mark₁、……Mark_n全部被处理，得到若干组类似Mark_i到Mark_m的数据；

S35、对每一组Mark_i到Mark_m的数据进行处理，具体为：

对应的修正后相位序列ph-fix_i到ph-fix_m中每个元素值分别减去ph-fix_i，得到数据组ph-fix_i-ph-fix_i、ph-fix_i+1-ph-fix_i、……、ph-fix_m-ph-fix_i，存入数据集合Y；

将对应的时间序列t_i到t_m中每个元素值分别减去t_i，得到数据组t_i-t_i、t_i+1-t_i、……、t_m-t_i，存入数据集合X；

S4、采用加权最小二乘法对处理后的数据进行直线拟合计算频率偏移量。

2.根据权利要求1所述的用于时钟驯服的频率偏移量分析方法，其特征在于，在S4中具体包括：

S41、建立一个数量为n的数据集合Weight，用来代表数据集合X、Y中对应元素的加权值,Weight具体加权值由最初的观测数据中时间序列t₀、t₁、……、t_n酌情拟定：t_n为最近的时刻，对应的权重设置为100%；t_n-1、t_n-2、……、t₀依次减少权重；

S42、带入权重采用最小二乘法直线拟合公式进行直线拟合计算频率偏移量，表示为。

3.根据权利要求1所述的用于时钟驯服的频率偏移量分析方法，其特征在于，在S31中，接收机报文内容包括授时有效标志位和可用卫星数量，当卫星数量大于4和授时有效标志位位于NMEA0183协议GGA报文中时，该观测数据则为正常数据，反之观测数据则为异常数据。