CN109428598B - 数据采样的校准方法和*** - Google Patents

Abstract

本发明公开了一种数据采样的时间同步、校准方法和***，其中，该方法包括：采用外部时钟源的实际时钟频率除以预设的目标采样率，得到实际采样间隔，实际采样间隔包括整数部分TIMERP_INT和小数部分TIMERP_FRAG；基于时间同步信号的固有时间延迟与计算实际采样间隔的时间延迟之和，计算得到固定延迟时间；采用固定延迟时间乘以时钟频率，得到固定延迟时间对应的时钟周期数；采用时钟周期数除以采样间隔的整数部分，得到商值和余数，记余数为FRAGO；将余数作为采样计数器在当前时间同步周期的第一次计数的初始计数值；基于初始计数值，控制AD采样器进行第一次采样。本发明能够减少后期校准的计算量。

Description

数据采样的校准方法和***

技术领域

本发明属于数据采样技术领域，尤其涉及一种数据采样的校准方法和***。

背景技术

全球电力骨干网络均是三相交流电力***，频率为50Hz或60Hz。为了保证电力***的稳定运行，需要对三相交流电的电压和电流进行准确的测量。目前，主要的测量技术是基于计算机的数字化测量技术，即通过电压和电流传感器(或者互感器)把高电压、大电流的信号(即电力***一次信号)转换为适合测量的幅值较小的电压或者电流信号(即电力***二次信号)，然后送入ADC(模数转换器)进行量化变成数字信号，并利用数字信号处理技术对采集到的信号进行相应的变换、计算，从而感知单个电力设备以及局部电网或者全局电网的工作状态，并进一步做故障判断与隔离、稳态与暂态控制、参数估计、故障预测等等应用。

由于电力***设备本身在物理空间上是分散的，因此对设备的数据采集也是分散进行的，然后通过有线或者无线通信把分散的数据汇集到主站，才能够做进一步计算。即智能电网的技术方向是通过布设大量的智能传感器，来感知各种电力设备(如变压器、线路、开关)的数据，并由单个设备的采样数据计算出电网的工作状态。

在对电力***进行数据采集以及工作状态计算时，由于电力信号在电力设备中的传输速度接近光速，即每秒30万公里，或者说每微秒(us)300米。为了准确感知电力***的整体状态，必须要获取同一时刻各个电力设备的数据，如果各个设备的采样数据在时间上差异太大，计算出的结果是没有意义的。比如在小电流接地配电网的单相接地故障定位这个应用领域，通过暂态零序电流进行故障定位是公认效果最好的故障定位方法，但是故障定位的前提是零序电流一定要准确，各相电流的采样时间偏差要控制在10us以内。在行波保护和雷电波定位这类应用领域，由于需要用时间差来计算空间距离，因此各个数据采集点的时间偏差要控制在1us以内。而由于电力***中各个设备在地理位置上是分散的，因此要保证它们的采样时间严格同步，或者说把采样时间的偏差控制在10us以内，甚至1us，这是非常有挑战的。为了满足电力信号测量精度的要求，现代电力***一般都采用较高的采样频率对电力信号进行高速采样，比如每个周期128、256、512个采样点。如果电力信号的频率是50Hz，则每秒钟的数据采样频率分别为6400,12800,25600次。

假如要求实现12800的数据采样频率，一般会用12800整数倍频的高频时钟进行分频，例如用25.6M(即25600000)的时钟进行2000倍分频，或者用12.8M的时钟进行1000倍分频，这样就生成了频率为12800的采样信号。

电子***上广泛使用低成本、低功耗的无源石英晶体振荡器(简称无源晶振，或者晶振)作为时钟源。石英晶体振荡器标称的理想频率和实际频率都是有一定偏差的，典型的频率偏差为±20PPM，这个偏差主要受晶体加工工艺的影响。在使用过程当中，石英晶体振荡器的频率还会受到温度的影响，随着时间的推移，石英晶体振荡器的频率还有老化现象。总之，常用的石英晶体振荡器的频率并非恒定，而是因个体差异、受工作环境影响、随时间缓慢变化的。

出于低成本和低功耗的考虑，电力***配电网的数据采集设备大量采用普通的石英晶体振荡器作为采集设备的时钟源。即使两个采集设备以一秒的周期进行严格的时间同步，即每秒初两个采集设备的时间是严格同步的，同步误差为0，但是由于±20PPM频率偏差的影响，在这一秒钟的末尾，每个设备的时间和标准时间相比会产生±20us的偏差，造成两个采集设备之间会累计产生最大40us的时间偏差。40us的时间偏差会造成合成的零序电流当中有比较大的合成误差，严重影响了故障判断的准确性。如果时间同步的周期大于一秒钟，累计的时间偏差会成倍增加。每次时间同步会都需要数据采集设备消耗额外的能量，对于功耗敏感的数据采集设备而言，频繁的时间同步是一个巨大的功耗负担。

为了满足电力信号测量精度的要求，现代电力***一般都采用较高的采样频率对电力信号进行高速采样，比如每个周期128、256、512个采样点。如果电力信号的频率是50Hz，则每秒钟的数据采样频率分别为6400,12800,25600次，如果电力信号的频率是60Hz，则每秒钟的数据采样频率分别为7680，15360，30720次。

假如要求实现12800的数据采样频率，一般会用12800整数倍频的高频时钟进行分频，例如用25.6M(即25600000)的时钟进行2000倍分频，或者用12.8M的时钟进行1000倍分频，这样就生成了频率为12800的采样信号。

电子***上广泛使用低成本、低功耗的无源石英晶体振荡器(简称无源晶振，或者晶振)作为时钟源。石英晶体振荡器标称的理想频率和实际频率都是有一定偏差的，典型的频率偏差为±20PPM，这个偏差主要受晶体加工工艺的影响。在使用过程当中，石英晶体振荡器的频率还会受到温度的影响，随着时间的推移，石英晶体振荡器的频率还有老化现象。总之，常用的石英晶体振荡器的频率并非恒定，而是因个体差异、受工作环境影响、随时间缓慢变化的。

出于低成本和低功耗的考虑，电力***配电网的数据采集设备大量采用普通的石英晶体振荡器作为采集设备的时钟源。即使两个采集设备以一秒的周期进行严格的时间同步，即每秒初两个采集设备的时间是严格同步的，同步误差为0，但是由于±20PPM频率偏差的影响，在这一秒钟的末尾，每个设备的时间和标准时间相比会产生±20us的偏差，造成两个采集设备之间会累计产生最大40us的时间偏差。40us的时间偏差会造成合成的零序电流当中有比较大的合成误差，严重影响了故障判断的准确性。如果时间同步的周期大于一秒钟，累计的时间偏差会成倍增加。每次时间同步会都需要数据采集设备消耗额外的能量，对于功耗敏感的数据采集设备而言，频繁的时间同步是一个巨大的功耗负担。

发明内容

(一)发明目的

本发明的目的是提供一种能够减少后期数据校准的计算量的数据采样的时间同步、校准方法和***。

(二)技术方案

为解决上述问题，本发明的第一方面提供了一种数据采样的时间同步、校准方法，包括：采用外部时钟源的实际时钟频率除以预设的目标采样率，得到实际采样间隔，所述实际采样间隔包括整数部分TIMERP_INT和小数部分TIMERP_FRAG；基于时间同步信号的固有时间延迟与计算所述实际采样间隔的时间延迟之和，计算得到固定延迟时间；采用所述固定延迟时间乘以所述时钟频率，得到所述固定延迟时间对应的时钟周期数；采用所述时钟周期数除以所述采样间隔的整数部分，得到商值和第一余数，记所述第一余数为FRAGO；将所述第一余数作为采样计数器在当前时间同步周期的第一次计数的初始计数值；基于所述初始计数值，控制AD采样器进行第一次采样。

根据本发明的另一个方面，提供一种数据采样的时间同步、校准***，包括：采样间隔计算模块，用于采用外部时钟源实际的时钟频率除以预设的目标采样率，得到实际采样间隔，所述实际采样间隔包括整数部分TIMERP_INT和小数部分TIMERP_FRAG；固定延迟时间计算模块，用于基于时间同步信号的固有时间延迟与计算所述实际采样间隔的时间延迟之和，计算得到固定延迟时间；时钟周期数计算模块，用于采用所述固定延迟时间乘以所述时钟频率，得到所述固定延迟时间对应的时钟周期数；采样预设值计算模块，用于采用所述时钟周期数除以所述采样间隔的整数部分，得到商值和第一余数，记所述第一余数为FRAGO，以及将基于所述第一余数转换后的数值作为采样计数器的预设值；控制模块，用于每当计数器达到预设值或预设值的整数倍时，控制AD采样器进行采样。

本发明通过首先计算固定延迟时间，并计算得到固定时间延迟对应的时钟周期数，且将固定时间延迟对应的时钟周期数处以实际采样间隔的整数部分得到商值和第一余数，将第一余数作为采样计数器的第一次计数的初始计数值，从而基于所述初始计数值，控制AD采样器进行第一次采样。

(三)有益效果

本发明的上述技术方案具有如下有益的技术效果：由于误差是由固定时间延迟和守时误差两部分构成，因此在将固定时间延迟在前期采样过程中就消除掉之后，就能够减小采样过程中理想采样点的采样数值误差，从而在保证数据采样的时间同步精度的前提下，能够减少插值计算的计算量，使得需要校准的数据减少，并节省了功耗。

附图说明

图1是本发明实施例一的数据采样的时间同步、校准方法流程图；

图2是本发明实施例二介绍的关于实施例一中步骤S2的一种具体实施方式的流程图；

图3是本发明实施例三介绍的关于实施例一中步骤S6的一种具体实施方式的流程图；

图4是本发明实施例四介绍的关于实施例一中步骤S1的一种具体实施方式的流程图；

图5是本发明实施例五的数据采样的时间同步、校准方法的流程图；

图6是本发明实施例五中关于时间同步周期的时间轴示意图；

图7是本发明前述实施例中的实际采样间隔对应的时钟周期数的一种具体实施方式的流程图；

图8是本发明实施例六的数据采样的时间同步、校准方法流程图；

图9是关于本发明实施例六中步骤S801的原理性示意图；

图10是关于本发明实施例六中步骤S805的一种具体实施方式的流程图；

图11是本发明一个示例的示意图；

图12是本发明实施例七的一种数据采样的时间同步、校准***的结构示意图；

图13是本发明实施例八的一种数据采样的时间同步、校准***的结构示意图；

图14是本发明实施例九的一种数据采样的时间同步、校准***的结构示意图；

图15是本发明实施例十的一种数据采样的时间同步、校准***的结构示意图；

图16是关于本发明实施例十中关于校准模块的一种具体实施例方式的结构示意图；

图17是关于本发明实施例十的校准模块中插值子模块的结构示意图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚明了，下面结合具体实施方式并参照附图，对本发明进一步详细说明。应该理解，这些描述只是示例性的，而并非要限制本发明的范围。此外，在以下说明中，省略了对公知结构和技术的描述，以避免不必要地混淆本发明的概念。

请参阅图1，是本发明实施例一的数据采样的校准方法流程图，包括如下步骤S1-S6：

S1，采用外部时钟源的实际时钟频率除以预设的目标采样率，得到实际采样间隔，所述实际采样间隔包括整数部分TIMERP_INT和小数部分TIMERP_FRAG；

具体地，外部时钟源的时钟频率是指每一秒钟之内外部时钟源的振荡次数。预设的目标采样率可以由用户自行设置。在大多数情况下，使用时钟频率除以目标采样率，都是无法整除的，因而得到的实际采样间隔就会存在整数部分和小数部分。例如：时钟频率为25600000，实际频率偏差是10PPM，则测量得到的实际的时钟频率为25600256，假如数据采样频率为12800，则实际采样间隔是2000.02个时钟周期。

S2，基于时间同步信号的固有时间延迟与计算实际采样间隔的时间延迟之和，计算得到固定延迟时间；

具体地，时间同步信号的作用是使各个采集器进行时间同步。时间同步信号可以是直接来自GPS/北斗接收模块的秒脉冲，也可以是通过有线方式传送的时间同步信息，还可以是通过无线方式传送的时间同步信息。上述三种不同的时间同步方式会有不同的时间同步延时，且延时的偏差也不同。这些方式当中，来自全球定位***的同步秒脉冲的精确度最高，可以达到广域10ns级别的同步精度。而采用有线方式可以达到1us～10us的同步精度，以及采用无线方式可以达到几十us～几百us的同步精度。具体地，根据时间同步延时的平均值可计算得到固有时间延时。

其中，计算实际采样间隔的时间延迟是指从接收到时间同步信号到完成采样间隔计算的时间延迟也可以固定下来。而当实际采样间隔的计算方法确定之后，其所使用的硬件电路也就定下来了，因为每个硬件电路的计算延时是和电路设计有关，因此计算实际采样间隔的时间延迟也就确定了。

S3，采用固定延迟时间乘以时钟频率，得到固定延迟时间对应的时钟周期数，记为FIXED_DLY_CLK；

S4，采用时钟周期数除以实际采样间隔的整数部分，得到商值和第一余数，记第一余数为FRAGO；

具体地，采用FIXED_DLY_CLK除以实际采样间隔的整数部分TIMERP_INT，得到商为N0,含义为固定延迟时间对应了多少个数据采样点。第一余数FRAG0,含义是除了N0个采样点以外，仍然有多少个时钟周期的延时。

S5，将第一余数作为采样计数器在当前时间同步周期的第一次计数的初始计数值；

S6，基于初始计数值，控制AD采样器进行第一次采样。

针对上述步骤S5和步骤S6，在具体实施过程中，由于实际的时钟频率是变化的，计算出的采样间隔时钟周期数不可能正好是整数，而是带有小数部分。因此，无法实现计数。需要通过以下方式对第一余数进行转换为能够用于计数的预设值，例如：小数部分数值为0.95，则通过以下方式计算得到经过处理后的小数部分数值：int(0.95*(1<<N_FRAG))，通过该计算方式能够得到能够用以小数计数器的计数值。具体地，如图2所示，本发明是采用如下方式计算实际采样间隔时钟周期数，得到整数部分以及经过处理能够用于计数器计数的小数部分：

S201，采用前一时间同步周期的实际时钟频率CLK_FREQ除以目标采样率SPS_TARGET，其中，目标采样率是预先设定好的，得到第一商值(即实际采样间隔的整数部分)，记为TIMERP_INT，第一余数(即实际采样间隔的余数部分)，记为CLK_REMAIN；

S202，采用前一时间同步周期的实际时钟频率CLK_FREQ除以第一商值TIMERP_INT，得到第二商值，记为SPS0；

S203，采用第一余数CLK_REMAIN除以第二商值SPS0，得到第三商值，记为TIMERP_FRAG。

具体地，在步骤S203中，当第一余数小于第二商值时，将第一余数左移N位，再采用左移N位后的第一余数除以第二商值，得到的结果的整数作为经过处理后的小数部分数值用以小数计数器进行计数。

在步骤S203中，由于第一余数CLK_REMAIN小于第二商值SPS0，因此如果直接做整数除法商是0。为了保证实际采样间隔时钟周期数足够的计算精度，在计算除法以前需要把第一余数CLK_REMAIN左移N_FRAG位，比如N_FRAG为8、16、20、24、或者32位，第一余数左移以后的结果再除以第二商值SPS0，得到的计算结果为第三商值TIMERP_FRAG的最高N_FRAG个小数位。

采用本发明的上述步骤S201-S203的方法计算得到的经过处理后的小数部分数值在用于小数计数器计数时，能够将多个采集器之间的时间偏差尽可能减小。

下面通过举例对上述实施方式进行解释说明：

例如：以标称频率为25600000，实际频率偏差是10PPM，则实际时钟频率为25600256。假如数据的目标采样率为12800，则实际采样间隔时钟周期数是2000.02个时钟周期。则计算方法如下：

S110：25600256/12800，第一商值TIMERP_INT是2000，第一余数CLK_REMAIN为256；

S111：25600256/2000,第二商值SPS0为12800；

S112：第一余数CLK_REMAIN除以第二商值SPS0，即256/12800，在做除法以前，假如左移的位数是24位，则实际的运算是(256<<24)/12800＝(2^32)/(2^8*50)＝(2^24/50)＝335544＝0x51EB8，即TIMERP_FRAG为335544，小数计数器的预设值就是(1<<24)即2^24。

通过小数计数器能够达到实际补偿的时钟周期数为：SPS0*TIMERP_FRAG/小数计数器的预设值,即335544*12800/2^24＝0xFFFFF00/2^24＝0xFFFFF00>>24＝0xFF,即255个时钟。而最初需要补偿的CLK_REMAIN为256，二者仅相差1个时钟周期。

需要说明的是：在S110，S111，S112的计算过程当中，所有的除法都是整数除法，每次除法都会得到相应的商和余数。在步骤S12中，可选地，可以将第一整数计数器的计数值减去整数部分数值以恢复初始状态重新开始计数，还可以是每当第一整数计数器的计数值等于整数部分的整数倍时，则触发小数计数器开始累加小数部分进行累加计数。

其中，时间同步周期是指自接收到一次全球定位***发送的时间同步信号开始，将当前时间作为时间同步周期的起始时间，将下一次全球定位***发送的时间同步信号的时间作为时间同步周期的结束时间，以此构成一个时间同步周期。具体地，除去第一个时间同步信号之外，其他时间同步信号的接收时间都是既作为前一时间同步周期的结束时间，同时也作为当前时间同步周期的开始时间。

其中，在进行计数时，是通过两个整数计数器和一个小数计数器来控制AD采样器的采样，两个整数计数器分别记为第一整数计数器和第二整数计数器，步骤S6的一种具体实施方式，如图3所示，包括：

S301，计算前一时间同步周期的实际采样间隔时钟周期数，实际采样间隔时钟周期数包括整数部分和小数部分；

S302，在当前时间同步周期内，当外部时钟源的第一个时钟周期信号到来时，分别触发第一、第二整数计数器在所述初始计数值的基础上进行加1计数；

S303，当第一整数计数器的计数值等于前一时间同步周期的实际采样间隔时钟周期数的整数部分数值时，触发小数计数器的计数值加上前一时间同步周期的实际采样间隔时钟周期数的小数部分转换得到的整数值，将所述第一整数计数器的计数值减去整数部分数值；

S304，当小数计数器的累加计数值大于或等于预设值时，触发第二整数计数器的计数值减1，将所述小数计数器的计数值减去所述预设值；

S305，当第二整数计数器的计数值等于前一时间同步周期的实际采样间隔时钟周期数的整数部分数值时，发送采样脉冲至采样器，以进行采样，且将第一、第二整数计数器和小数计数器清零作为各自的初始计数值重新开始计数，将所述第二整数计数器的计数值减去整数部分数值。

本发明实施例通过基于时间同步信号的固有时间延迟与计算实际采样间隔的时间延迟之和，计算得到固定延迟时间；进而采用固定延迟时间乘以时钟频率，得到所述固定延迟时间对应的时钟周期数；采用时钟周期数除以采样间隔的整数部分，得到商值和余数，记余数为FRAGO；将基于余数转换后的数值作为采样计数器的预设值，以控制AD采样器进行采样。由于将固定时间延迟提前在时间同步的过程中就进行了消除，从而不需要在数据校准的过程中再考虑此部分固定时间延迟造成的误差，这样可以在保证数据采样同步精度的前提下，进一步减少插值计算的计算量，从而节省了功耗。

请参阅图4，是本发明是实例四的方法流程图，本发明实施例四是本发明实施例一的步骤S1的一种具体实施方式，步骤S1包括步骤S401-S403：

S401，采用时钟计数器对外部时钟源的时钟周期信号进行加1计数；

在本发明该实施例中，需要设置一时钟计数器对外部时钟源的振荡周期进行计数，时钟周期信号则是指外部时钟源的振荡次数，时钟源每振荡一次则相当于一个时钟周期的结束点以及另一个时钟周期的起始点，相邻两次时钟周期信号构成一个时钟周期，每当时钟源发出一个时钟周期信号，则时钟计数器加1进行计数。

S402，每当接收到外部的时间同步信号时，读取所述时钟计数器的计数值；

优选地，时间同步信号来自GPS/北斗接收模块的秒脉冲、或通过有线方式传送的时间同步信号、或通过无线方式传送的时间同步信号。由于不同的时间同步方式会有不同的时间同步延时，并且延时的偏差也不同。这些方式当中，来自全球定位***的同步脉冲的精确度是最高的，可以达到广域10ns级别的同步精度。采用有线方式可以达到1us以内的同步精度。采用无线方式可以达到10us以内的同步精度。时间同步的固定时间延迟可以通过数据校准算法进行纠正，但是时间延迟的的随机性偏差是无法进行消除的。上述时间偏差均为固定时间延迟，对于数据采样来说，时间同步误差是由固定时间延迟和守时误差共同决定的，守时误差是由晶体振荡器的一些原因造成的误差，可以通过小数计数器来减小，当时间同步误差能够控制得足够小时，且守时误差也可以控制的足够小时，则时间同步精度就可以控制在1us以内。

S403，采用相邻两次读取到的计数值的差值除以接收到的相邻两次时间同步信号的时间间隔，计算得到所述外部时钟源实际的时钟频率。

由于计算得到的实际采样间隔时钟周期数不是正好为整数，因此其包含整数部分和小数部分。

请参阅图5，为本发明实施例五的方法流程图，本发明实施例五是在实施例一的基础上，还包括如下步骤S501-S505：

S501，采用时钟计数器对外部时钟源的时钟周期信号进行加1计数；

具体地，是采用时钟计数器对外部时钟源进行计数。每当外部时钟源有效，时钟计数器的计数就增加一。每一秒钟之内时钟计数器增加的数值就是准确的时钟频率。外部时钟源可以是采用晶体振荡器。由于两次时间同步信号的接收时间间隔比较短(一般为一秒钟到一分钟)，在此期间外部环境温度的变化不大，并且可以不考虑外部晶体振荡器的老化，因此可以认为晶体振荡器的频率基本是稳定不变的。

S502，每当接收到来自外部的时间同步信号时，读取时钟计数器的计数值；

优选地，时间同步信号来自GPS/北斗接收模块的秒脉冲、或通过有线方式传送的时间同步信号、或通过无线方式传送的时间同步信号。由于不同的时间同步方式会有不同的时间同步延时，并且延时的偏差也不同。这些方式当中，来自全球定位***的同步脉冲的精确度是最高的，可以达到广域10ns级别的同步精度。采用有线方式可以达到1us以内的同步精度。采用无线方式可以达到10us以内的同步精度。时间同步的固定时间延迟可以通过数据校准算法进行纠正，但是时间延迟的的随机性偏差是无法进行消除的。上述时间偏差均为固定时间延迟，对于数据采样来说，时间同步误差是由固定时间延迟和守时误差共同决定的，守时误差是由晶体振荡器的一些原因造成的误差，可以通过小数计数器来减小，当时间同步误差能够控制得足够小时，且守时误差也可以控制的足够小时，则时间同步精度就可以控制在1us以内。

S503，基于当前读取到的计数值和前一次读取到计数值，计算相邻两次计数值的差值绝对值；

例如：如图6所示，各个时间同步周期(T1，T2，T3)的时钟计数值分别是CLK_CNT1，CLK_CNT2，CLK_CNT3。每当接收到时间同步信号时，就读取时钟计数值，分别计算ABS(CLK_CNT2-CLK_CNT1),ABS(CLK_CNT3-CLK_CNT2),即计算出相邻两个时间同步周期的时钟计数值的差值的绝对值DELTA_CNT，该差值反映了相邻两个时间同步周期时钟频率的变化。

S504，判断相邻两次时钟源的计数值的差值绝对值是否大于或等于预设的时钟计数偏差值；

具体地，预设的时钟计数偏差值MAX_CNT的计算公式是：允许的最大采样误差(us)*时钟频率(MHz)，比如最大采样误差是1us，时钟频率为25.6MHz，则预设的时钟计数偏差值MAX_CNT＝1*25.6＝25.6。

S505，若相邻两次时钟源的计数值的差值绝对值大于或等于预设的时钟计数偏差值，则对采样的数据进行校准。

具体地，若相邻两次时钟源的计数值的差值绝对值大于或等于预设的时钟计数偏差值，则该时间同步周期内最后一个采样点的时间偏差超出了允许的范围，则需要进行采样数据的校准。

在本发明实施例的基础上，S504还存在其他的实施方式，包括：S506，若相邻两次时钟源的计数值的差值绝对值小于预设的时钟计数偏差值，则不对采样数据进行校准。

本发明上述实施例通过判断相邻两次时钟源的计数值的差值绝对值是否大于或等于预设的时钟计数偏差值，来决定是否需要数据校准，从而使得数据校准方法是按需进行的。由于电力***大部分时刻都处于正常工作状态，采样到的数据不会真正去进行计算，这些采样数据在一段时间之内，会自动被覆盖。出于低功耗的考虑，只有当采样数据需要被使用时，才会进行数据的校准，并且把校准以后的采样数据送给计算平台进行计算。那些没有被使用的采样数据，则不会进行采样数据的校准。这样大大减少了数据校准所消耗的功耗，同时也保证了真正关心的采样值的时间同步和数据准确度。

请参阅图7，如图7所示，实际采样间隔对应的时钟周期数是按照以下方法计算得到的：

S71，每当接收到时钟周期信号时，采用时钟计数器进行加1计数；

在本发明该实施例中，需要设置一时钟计数器对外部时钟源的振荡周期进行计数，时钟周期信号则是指外部时钟源的振荡次数，每当时钟源发出一个时钟周期信号，则时钟计数器加1进行计数。

S72，每当接收到时间同步信号时，读取时钟计数器的当前计数值；

优选地，时间同步信号来自GPS/北斗接收模块的秒脉冲、或通过有线方式传送的时间同步信号、或通过无线方式传送的时间同步信号。

S73，基于读取到的当前计数值和前一次读取到计数值之差，以及读取两次计数值的时间间隔，计算实际时钟频率；

S74，采用实际时钟频率乘以目标采样率，计算得到实际采样间隔对应的时钟周期数。

由于计算得到的实际采样间隔对应的时钟周期数不是正好为整数，因此其包含整数部分和小数部分。

请参阅图8，为本发明实施例六的方法流程图，本发明实施例六是在实施例五的基础上，还包括如下步骤S801-S805：

S801，在当前时间同步周期，基于前一时间同步周期的时间同步信号到来时读取的时钟计数器的计数值和目标数据采样率，计算得到实际采样时间间隔；

具体地，时间同步周期是分别以接收到的相邻两次时间同步信号为该周期的起始时间点和结束时间点。

如图9所示，S[0]、S[1]、S[2]为实际的采样时刻，其时间间隔是CLK_CNT1除以目标数据采样率。计算公式为：timerp1＝CLK_CNT1/SPS_TARGET；式中，SPS_TARGET为目标数据采样率，CLK_CNT1为第二个时间同步信号到来时时钟计数器的计数值，timerp1为实际采样时间间隔。

S802，在当前时间同步周期，基于当前时间同步周期的时间同步信号到来时读取的时钟计数器的计数值和目标数据采样率，计算得到理想采样时间间隔；

如图9所示，RS[0]、RS[1]、RS[2]、RS[3]、RS[4]是理想采样时刻的采样值，RS[0]、RS[1]、RS[2]、RS[3]、RS[4]之间的时间间隔是CLK_CNT2除以目标数据采样率SPS_TARGET。计算公式为：timerp2＝CLK_CNT2/SPS_TARGET；式中，SPS_TARGET为目标数据采样率，CLK_CNT2为第三个时间同步信号到来时时钟计数器的计数值，timerp2为理想采样时间间隔。

S803，采用理想采样时间间隔除以目标采样率，计算得到各个理想采样点对应的时钟时刻；

具体地，仍然以实施例一中所举例为示例，各个理想采样点对应的时钟时刻CLK_I的计算方法是：CLK_I＝timerp2*I；式中，CLK_I为T2时间段中第I个理想采样点对应的时钟时刻。

S804，采用各个理想采样点对应的时钟时刻除以所述实际采样时间间隔，计算得到各个理想采样点对应的时钟时刻对应的实际采样点的位置，记为J0；

具体地，计算方法是：J0＝CLK_I/timerp1；式中，J0为第I个理想采样点对应的时钟时刻对应的实际采样点。

S805，从当前采样周期1个理想采样点开始，基于所述各个理想采样点对应的实际采样点J0的采样值和与实际采样点J0相邻的实际采样点J1的采样值，以及相邻两个实际采样点的位置分别与当前理想采样点I之间的权重系数，利用插值算法计算得到各理想采样点的校准采样值。其中，插值算法可以采用线性插值、正弦插值、或者其它的非线性插值算法。但一般情况下，线性插值的插值精度就可以满足校准精度的要求。

请参阅图10，为本发明实施例七的方法流程图，本发明实施例七是在实施例六的基础上，步骤S805中相邻两个实际采样点的位置分别与当前理想采样点I之间的权重系数的计算方式的一种具体实施方式，包括S1001-S1003：

S1001，将当前理想采样点I对应的实际采样点J0位置取整；

S1002，将当前理想采样点I对应的实际采样点J0位置减去取整部分的数值，得到当前理想采样点I和另一实际采样点J1之间的权重系数，记为Q1；

S1003，将1减去Q1得到当前理想采样点I与计算得到的与当前理想采样点I对应的实际采样点J0之间的权重系数。

具体地，RS[I]为当前时间同步周期的第I个理想采样点的采样值，当前理想采样点I对应的时钟时刻对应的实际采样点位置为J0，根据J0的位置选取与其相邻的另一个实际采样点J1，其中，J0和J1分别位于当前理想采样点的两端；分别计算实际采样点J1与理想采样点N0+I之间的权重J0-int(J0)和实际采样点J0与理想采样点N0+I之间的权重1-(J0-int(J0))进行插值计算得到的。计算方法是：RS[I]＝插值算法(S[int(J0)],S[int(J0)+1],J0-int(J0)),式中，RS[I]为T2时间段中第I个理想采样点的采样值；S[int(J0)]为T2时间段第J0个实际采样点的采样值，S[int(J0)+1]为T2时间段第J1个实际采样点的采样值，J0-int(J0)为第J1个实际采样点和当前理想采样点之间的权重系数，1-(J0-int(J0))为与J0相邻的另一实际采样点J1与当前理想采样点之间的权重系数。

下面以线性差值为例对上述插值计算过程进行说明：

假如当前理想采样点I对应的实际采样点位置是10.3，要利用S[10]和S[11]进行插值，J0-int(J0)＝10.3-10＝0.3。也就是说理想采样点I和S[11]的权重系数是0.3，理想采样点N0+I和S[10]的权重系数就是(1-0.3)＝0.7。那么线性差值的结果就是0.7*S[10]+0.3*S[11]＝S[10]-0.3*S[10]+0.3*S[11]＝S[10]+0.3*(S[11]-S[10])。

本发明上述实施例的方法可以自动适应石英晶体振荡器的频率偏差，降低了数据采样的同步误差，把采样的同步误差控制在1us之内，提高了采样数据的可用性。

上述实施例的采样数据校准过程，可以在数据采集设备端进行，也可以在数据汇集设备端进行，也可以在云平台进行。如果不在数据采集设备进行数据校准，则需要数据采集设备提供用于数据校准的参数和原始采样数据。

如图11所示，通过下述示例对本发明的技术效果进行说明：

图11中，实心箭头指向的0us为理想采样点，

假设时间同步周期为50us，则实心箭头指向的50us、100us、150us和200us均为实际采样点，假设固定时间延迟为15us，则15us、65us、115us和165us均为现有技术的实际采样点，通过本发明的技术方案能够使得后期在校准时不再需要考虑固定时间延迟带来的误差，只需针对守时误差(指各个采样器在进行时间同步以后是否能够守住该同步后的时间)进行校准即可，能够减少后期校准的计算量，通过将固定时间延迟从触发采样器采样的第二整数计数器的计数值中减掉，使得在前期采样时就将固定时间延迟带来的采样误差消除掉。

如图12所示，一种数据采样的时间同步、校准***，包括：

采样间隔计算模块1，用于采用外部时钟源实际的时钟频率除以预设的目标采样率，得到实际采样间隔，所述实际采样间隔包括整数部分TIMERP_INT和小数部分TIMERP_FRAG；

固定延迟时间计算模块2，用于基于时间同步信号的固有时间延迟与计算所述实际采样间隔的时间延迟之和，计算得到固定延迟时间；

时钟周期数计算模块3，用于采用所述固定延迟时间乘以所述时钟频率，得到所述固定延迟时间对应的时钟周期数；

采样预设值计算模块4，用于采用所述时钟周期数除以所述采样间隔的整数部分，得到商值和余数，记所述余数为FRAGO，以及将基于所述余数转换后的数值作为采样计数器6的预设值；

控制模块5，用于每当所述计数器达到所述预设值或所述预设值的整数倍时，控制AD采样器进行采样。

如图13所示，所述采样计数器6包括第一整数计数器61、第二整数计数器62和小数计数器63；

所述***还包括外部时钟源7：

时钟源7，当自身第一个时钟周期信号发出时，分别触发第一、第二整数计数器在所述初始计数值的基础上进行加1计数；

第一整数计数器61，用于当其计数值等于前一时间同步周期的实际采样间隔时钟周期数的整数部分数值时，触发小数计数器63的计数值加上前一时间同步周期的实际采样间隔时钟周期数的小数部分转换得到的整数值，且将所述第一整数计数器的计数值减去整数部分数值；

小数计数器63，当其累加计数值大于或等于预设值时，触发第二整数计数器62的计数值减1，且将所述小数计数器的计数值减去所述预设值；

第二整数计数器62，用于当其计数值等于前一时间同步周期的实际采样间隔时钟周期数的整数部分数值时，发送采样脉冲至采样器，以进行采样，且将所述第二整数计数器的计数值减去整数部分数值。

如图14所示，一种数据采样的时间同步、校准***还包括：

时钟计数器8，用于对外部时钟源的时钟周期信号进行加1计数；

计数读取模块9，用于每当接收到外部的时间同步信号时，读取所述时钟计数器的计数值；

时钟频率计算模块10，用于采用相邻两次读取到的计数值的差值除以接收到的相邻两次时间同步信号的时间间隔，计算得到所述时钟频率。

如图15所示，一种数据采样的时间同步、校准***还包括：

时钟计数器8，用于对外部时钟源的时钟周期信号进行加1计数；

计数读取模块9，用于每当接收到外部的时间同步信号时，读取所述时钟计数器8的计数值；

差值绝对值计算模块11，用于基于当前读取到的计数值和前一次读取到计数值，计算相邻两次计数值的差值绝对值；

判断模块12，用于判断所述相邻两次时钟源的计数值的差值绝对值是否大于或等于预设的时钟计数偏差值；

校准模块13，用于当所述相邻两次时钟源的计数值的差值绝对值大于或等于预设的时钟计数偏差值时，对采样的数据进行校准。

如图16所示，校准模块13还包括：

实际采样时间间隔计算子模块131，用于在当前时间同步周期，基于前一时间同步周期的时间同步信号到来时读取的时钟计数器的计数值和目标数据采样率，计算得到实际采样时间间隔；

理想采样时间间隔计算子模块132，用于在当前时间同步周期，基于当前时间同步周期的时间同步信号到来时读取的时钟计数器的计数值和目标数据采样率，计算得到理想采样时间间隔；

理想采样点时间计算子模块133，用于采用所述理想采样时间间隔除以目标采样率，计算得到各个理想采样点对应的时钟时刻；

理想采样点对应实际采样点的计算子模块134，用于采用所述各个理想采样点对应的时钟时刻除以所述实际采样时间间隔，计算得到各个理想采样点对应的时钟时刻对应的实际采样点的位置，记为J0；

插值子模块135，用于从当前时间同步周期的第1个理想采样点开始，基于所述各个理想采样点对应的实际采样点J0的采样值和与实际采样点J0相邻的实际采样点J1的采样值，以及相邻两个实际采样点的位置分别与当前理想采样点I之间的权重系数，利用插值算法计算得到各理想采样点的校准采样值。

如图17所示，插值子模块135包括：

取整单元1350，用于将当前理想采样点I对应的实际采样点J0位置取整；

第一权重系数计算单元1351，用于将当前理想采样点I对应的实际采样点J0位置减去取整部分的数值，得到当前理想采样点I和另一实际采样点J1之间的权重系数，记为Q1；

第二权重系数计算单元1352，用于将1减去Q1得到当前理想采样点I与计算得到的与当前理想采样点I对应的实际采样点J0之间的权重系数。

需要说明的是，本发明数据采样的校准***是与涉及计算机程序流程的数据采样的校准方法一一对应的装置/***，由于在前已经对数据采样的校准方法的步骤流程进行了详细描述，在此不再对数据采样的校准***的实施过程进行赘述。

应当理解的是，本发明的上述具体实施方式仅仅用于示例性说明或解释本发明的原理，而不构成对本发明的限制。因此，在不偏离本发明的精神和范围的情况下所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。此外，本发明所附权利要求旨在涵盖落入所附权利要求范围和边界、或者这种范围和边界的等同形式内的全部变化和修改例。

Claims (12)

1.一种数据采样的时间同步、校准方法，其特征在于，包括：

采用外部时钟源的实际时钟频率除以预设的目标采样率，得到实际采样间隔，所述实际采样间隔包括整数部分TIMERP_INT和小数部分TIMERP_FRAG；

基于时间同步信号的固有时间延迟与计算所述实际采样间隔的时间延迟之和，计算得到固定延迟时间；

采用所述固定延迟时间乘以所述时钟频率，得到所述固定延迟时间对应的时钟周期数；

采用所述时钟周期数除以所述实际采样间隔的整数部分，得到商值和第一余数，记所述第一余数为FRAGO；

将所述第一余数作为采样计数器在当前时间同步周期的第一次计数的初始计数值；

基于所述初始计数值，控制AD采样器进行第一次采样。

2.根据权利要求1所述的时间同步、校准方法，其特征在于，所述采样计数器包括第一、第二整数计数器和小数计数器；

其中，基于所述初始计数值，控制AD采样器进行采样的步骤包括：

在当前时间同步周期内，当外部时钟源的第一个时钟周期信号到来时，分别触发第一、第二整数计数器在所述初始计数值的基础上进行加1计数；

当第一整数计数器的计数值等于前一时间同步周期的实际采样间隔时钟周期数的整数部分数值时，触发小数计数器的计数值加上前一时间同步周期的实际采样间隔时钟周期数的小数部分转换得到的整数值，将所述第一整数计数器的计数值减去整数部分数值；

当小数计数器的累加计数值大于或等于预设值时，触发第二整数计数器的计数值减1，将所述小数计数器的计数值减去所述预设值；

当第二整数计数器的计数值等于前一时间同步周期的实际采样间隔时钟周期数的整数部分数值时，发送采样脉冲至采样器，以进行采样，将所述第二整数计数器的计数值减去整数部分数值。

3.根据权利要求1所述的时间同步、校准方法，其特征在于，所述外部时钟源的实际时钟频率根据以下方法计算：

采用时钟计数器对外部时钟源的时钟周期信号进行加1计数；

每当接收到外部的时间同步信号时，读取所述时钟计数器的计数值；

采用相邻两次读取到的计数值的差值除以接收到的相邻两次时间同步信号的时间间隔，计算得到所述外部时钟源实际的时钟频率。

4.根据权利要求1-3任一项所述的时间同步、校准方法，其特征在于，所述方法还包括：

采用时钟计数器对外部时钟源的时钟周期信号进行加1计数；

每当接收到来自外部的时间同步信号时，读取所述时钟计数器的计数值；

基于当前读取到的计数值和前一次读取到计数值，计算相邻两次计数值的差值绝对值；

判断相邻两次时钟源的计数值的差值绝对值是否大于或等于预设的时钟计数偏差值；

若相邻两次时钟源的计数值的差值绝对值大于或等于预设的时钟计数偏差值，则对采样的数据进行校准。

5.根据权利要求4所述的时间同步、校准方法，其特征在于，对所述采样数据进行校准的步骤包括：

在当前时间同步周期，基于前一时间同步周期的时间同步信号到来时读取的时钟计数器的计数值和目标数据采样率，计算得到实际采样时间间隔；

在当前时间同步周期，基于当前时间同步周期的时间同步信号到来时读取的时钟计数器的计数值和目标数据采样率，计算得到理想采样时间间隔；

采用所述理想采样时间间隔除以目标采样率，计算得到各个理想采样点I对应的时钟时刻；

采用所述各个理想采样点对应的时钟时刻除以所述实际采样时间间隔，计算得到各个理想采样点对应的时钟时刻对应的实际采样点的位置，记为J0；

从当前采样周期的第1个理想采样点开始，基于所述各个理想采样点对应的实际采样点J0的采样值和与实际采样点J0相邻的实际采样点J1的采样值，以及相邻两个实际采样点的位置分别与当前理想采样点I之间的权重系数，利用插值算法计算得到各理想采样点的校准采样值。

6.根据权利要求5所述的时间同步、校准方法，其特征在于，所述相邻两个实际采样点的位置分别与当前理想采样点I之间的权重系数按照以下方式计算：

将当前理想采样点I对应的实际采样点J0位置取整；

将当前理想采样点I对应的实际采样点J0位置减去取整部分的数值，得到当前理想采样点I和另一实际采样点J1之间的权重系数，记为Q1；

将1减去Q1得到当前理想采样点I与计算得到的与当前理想采样点I对应的实际采样点J0之间的权重系数。

7.一种数据采样的时间同步、校准***，其特征在于，包括：

采样间隔计算模块(1)，用于采用外部时钟源实际的时钟频率除以预设的目标采样率，得到实际采样间隔，所述实际采样间隔包括整数部分TIMERP_INT和小数部分TIMERP_FRAG；

固定延迟时间计算模块(2)，用于基于时间同步信号的固有时间延迟与计算所述实际采样间隔的时间延迟之和，计算得到固定延迟时间；

时钟周期数计算模块(3)，用于采用所述固定延迟时间乘以所述时钟频率，得到所述固定延迟时间对应的时钟周期数；

采样预设值计算模块(4)，用于采用所述时钟周期数除以所述采样间隔的整数部分，得到商值和第一余数，记所述第一余数为FRAGO，以及将基于所述第一余数转换后的数值作为采样计数器(6)的预设值；

控制模块(5)，用于每当所述计数器达到所述预设值或所述预设值的整数倍时，控制AD采样器进行采样。

8.根据权利要求7所述的时间同步、校准***，其特征在于，所述采样计数器(6)包括第一整数计数器(61)、第二整数计数器(62)和小数计数器(63)；

所述***还包括：

外部时钟源(7)，当自身第一个时钟周期信号发出时，分别触发第一、第二整数计数器在初始计数值的基础上进行加1计数；

第一整数计数器(61)，用于当其计数值等于前一时间同步周期的实际采样间隔时钟周期数的整数部分数值时，触发小数计数器(63)的计数值加上前一时间同步周期的实际采样间隔时钟周期数的小数部分转换得到的整数值，且将所述第一整数计数器的计数值减去整数部分数值；

小数计数器(63)，当其累加计数值大于或等于所述预设值时，触发第二整数计数器(62)的计数值减1，且将所述小数计数器的计数值减去所述预设值；

第二整数计数器(62)，用于当其计数值等于前一时间同步周期的实际采样间隔时钟周期数的整数部分数值时，发送采样脉冲至采样器，以进行采样，且将所述第二整数计数器的计数值减去整数部分数值。

9.根据权利要求7所述的时间同步、校准***，其特征在于，所述***还包括：

时钟计数器(8)，用于对外部时钟源的时钟周期信号进行加1计数；

计数读取模块(9)，用于每当接收到外部的时间同步信号时，读取所述时钟计数器的计数值；

时钟频率计算模块(10)，用于采用相邻两次读取到的计数值的差值除以接收到的相邻两次时间同步信号的时间间隔，计算得到所述时钟频率。

10.根据权利要求7-9任一项所述的时间同步、校准***，其特征在于，所述***还包括：

时钟计数器(8)，用于对外部时钟源的时钟周期信号进行加1计数；

计数读取模块(9)，用于每当接收到外部的时间同步信号时，读取所述时钟计数器(8)的计数值；

差值绝对值计算模块(11)，用于基于当前读取到的计数值和前一次读取到计数值，计算相邻两次计数值的差值绝对值；

判断模块(12)，用于判断相邻两次时钟源的计数值的差值绝对值是否大于或等于预设的时钟计数偏差值；

校准模块(13)，用于当相邻两次时钟源的计数值的差值绝对值大于或等于预设的时钟计数偏差值时，对采样的数据进行校准。

11.根据权利要求10所述的时间同步、校准***，其特征在于，所述校准模块(13)还包括：

实际采样时间间隔计算子模块(131)，用于在当前时间同步周期，基于前一时间同步周期的时间同步信号到来时读取的时钟计数器的计数值和目标数据采样率，计算得到实际采样时间间隔；

理想采样时间间隔计算子模块(132)，用于在当前时间同步周期，基于当前时间同步周期的时间同步信号到来时读取的时钟计数器的计数值和目标数据采样率，计算得到理想采样时间间隔；

理想采样点时间计算子模块(133)，用于采用所述理想采样时间间隔除以目标采样率，计算得到各个理想采样点对应的时钟时刻；

理想采样点对应实际采样点的计算子模块(134)，用于采用所述各个理想采样点对应的时钟时刻除以所述实际采样时间间隔，计算得到各个理想采样点对应的时钟时刻对应的实际采样点的位置，记为J0；

插值子模块(135)，用于从当前时间同步周期的第1个理想采样点开始，基于所述各个理想采样点对应的实际采样点J0的采样值和与实际采样点J0相邻的实际采样点J1的采样值，以及相邻两个实际采样点的位置分别与当前理想采样点I之间的权重系数，利用插值算法计算得到各理想采样点的校准采样值。

12.根据权利要求11所述的时间同步、校准***，其特征在于，所述插值子模块(135)包括：

取整单元(1350)，用于将当前理想采样点I对应的实际采样点J0位置取整；

第一权重系数计算单元(1351)，用于将当前理想采样点I对应的实际采样点J0位置减去取整部分的数值，得到当前理想采样点I和另一实际采样点J1之间的权重系数，记为Q1；

第二权重系数计算单元(1352)，用于将1减去Q1得到当前理想采样点I与计算得到的与当前理想采样点I对应的实际采样点J0之间的权重系数。

Images