一种电力谐波参数软件同步采样方法技术

本发明专利技术公开一种电力谐波参数软件同步采样方法，包括如下步骤：a1：设置每周期采样点数N，N为每周期的采样点数；a2：测出前W个周期时间长度Xi，并存储，其中i＝1、2…W；a3：根据测量的前W个周期时间长度，计算平滑系数α，并储存α值；a4：根据三次指数平滑模型，预测下一采样周期的时间长度；a5：根据公式T

【技术实现步骤摘要】
一种新型电力谐波参数软件同步采样方法
本专利技术涉及电能质量参数测量领域，特别涉及一种软件同步采样的控制方法。
技术介绍
同步采样是指信号周期T，采样周期间隔Ts，及采样点数N严格满足等式T＝NTs的采样方法。理论上，只要满足T＝NTs并且N>2M(M为被测信号最高谐波次数)，同步采样不存在测量方法上的误差，然而由于信号周期T的测量并不总是准确的，并且电网信号的频率总是在一定范围内波动，因此在实际采样过程中，采样周期与被测信号周期不能实现同步，而产生频谱泄露现象和栅栏效应。同步采样的方法分为硬件同步采样方法和软件同步采样方法两种。其中软件同步采样方法较硬件同步采样方法更为简单，但是目前所应用的软件同步采样方法计算量较为复杂，或精确度不能满足要求。因此，为确保谐波、间谐波测量的准确性，需要根据之前测得的信号周期时间长度序列对当前测量周期的时间长度进行预测。因此，如何通过适当的算法，选用合适的数据进行预测，即本领域技术人员需要解决的技术问题。并且，针对上述问题，如何开发出新的同时满足精确度要求并且较为简单的同步采样方法，来满足谐波、间谐波测量的要求，即为本领域技术人员的研究方向所在。
技术实现思路
本专利技术的主要目的是为了更准确的预测当前测量周期的时间长度，并且满足准确性和实时性的要求，减少计算量，进而提出的一种新型电力谐波参数软件同步采样方法，该方法适用于50Hz电网内，为谐波、间谐波测量采集信号。为了达到上述目的，本专利技术提供一种新型电力谐波参数软件同步采样方法，包括如下步骤：步骤a1：设置每周期采样点数N，N为每周期的采样点数；步骤a2：测出前W个周期时间长度Xi，并存储，其中i＝1、2…W；步骤a3：根据测量的前W个周期时间长度，计算平滑系数α，并储存α值；步骤a4：根据三次指数平滑模型，预测下一采样周期的时间长度；步骤a5：根据公式Ts＝T/N，计算出采样时间间隔，进行下一周期采样，并储存采样值；步骤a6：判断已采样周期数是否达到W个，若是，进入步骤a7，若否，进入步骤a8；步骤a7：去掉最早一周期数据，加入最新一周期的采样值，通过离散傅里叶变换方法计算谐波参数，返回步骤a8；步骤a8：去掉最早一周期数据，通过过零点检测，存储最新一周期的周期时间长度，返回步骤a3。较佳的实施方式，在步骤a1中，设置一个测量周期内的采样点数为N，其中保证N>2M，M为被测信号最高谐波次数。较佳的实施方式，所述步骤a2的检测方式是通过过零点检测进行。较佳的实施方式，在步骤a3中，当误差平方和最小时，此时的平滑系数α为最准确，最小误差平方公式其中K表示预测的次数，表示第t次预测值，Xt表示第t次实际值。较佳的实施方式，在步骤a4中，三次指数平滑模型公式为：其中，St(1)为第t个周期的一次指数平滑值，St(2)为第t个周期的二次指数平滑值，St(3)为第t个周期的三次指数平滑值，α为平滑系数，Xt为时间序列{Xi}第t个周期的的实际时间长度。较佳的实施方式，在步骤a7中，采用离散傅里叶变换的方法来计算谐波参数。较佳的实施方式，最优平滑系数α的选择包括如下步骤：b1：选定t个连续时段，记录其测量值，并获得平滑初值；b2：设定步长λ＝0.01，将α∈(0,1)区间分成100个距离相同的区间，区间端点值依次记为αk，k＝1,2，…；b3：计算出每个αk对应的f值；b4：选择最小的f值，其对应的αk为最优平滑系数α；b5：预测并存储第t+1个时间段的数据以及最优平滑系数α；b6：去掉最早的数据，加入最新的t+1个周期的数据，得到新的平滑初值，并返回步骤b2。较佳的实施方式，对误差平方f值的计算采用地毯式遍历搜索的方法。较佳的实施方式，选取每测量十周期计算一次最优平滑系数，选用的数据同样为当前测量周期的前十个周期的预测时间长度和实际测量时间长度。较佳的实施方式，改变预测步长，依次预测当前周期的后m个周期的时间长度。与现有技术相比，本专利技术的有益效果在于：本专利技术可更准确的预测当前测量周期的时间长度，并且满足准确性和实时性的要求，同时减少计算量，适用于50Hz电网内，为谐波、间谐波测量采集信号附图说明图1为专利技术一种新型电力谐波参数软件同步采样方法流程图；图2最优平滑系数计算流程图。具体实施方式下面将参考附图并结合实施例来详细说明本专利技术。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。如图1所示，为本专利技术一种新型电力谐波参数软件同步采样方法流程图，其具体流程如下：a1：设置每周期采样点数N；在此步骤中，可设置一个测量周期内的采样点数为N，其中保证N>2M(M为被测信号最高谐波次数)。a2：测出前W个周期时间长度Xi(i＝1、2…W)，并存储，该检测方式可通过过零点检测进行，过零点检测采用比较器过零比较电路，可以更准确的进行过零检测。；在此实施例中，取W＝10为例，即通过测量前十个周期的周期时间长度，得到前十个周期的时间长度序列{Xi}，i＝1,2,…10，计算得到前十个周期的时间长度平均值为a3：根据步骤a2测量的前十个周期时间长度，计算平滑系数α，并储存α值；平滑系数α取值范围为α∈(0,1)，反映了不同时期历史数据在平滑预测时所占得比重，当α为定值时，会使得平滑模型难以跟踪周期时间长度序列自身的变化，这样会导致预测不够准确。所以本专利技术提出一种动态平滑系数的计算方法。如附图2所示，来选取平滑系数的值。根据经验当误差平方和最小时，此时的平滑系数α为最准确。最小误差平方公式其中K表示预测的次数，表示第t次预测值，Xt表示第t次实际值。a4：根据三次指数平滑模型，预测下一采样周期的时间长度；在此步骤中，是通过三次指数平滑模型，来对下一采样周期的时间长度进行预测，由于电网频率变化具有非线性和无规则性的特点，所以选择三次指数平滑模型进行预测。三次指数平滑模型公式为：其中St(1)为第t个周期的一次指数平滑值，St(2)为第t个周期的二次指数平滑值，St(3)为第t个周期的三次指数平滑值，α为平滑系数，Xt为时间序列{Xi}第t个周期的的实际时间长度。以十个周期的周期时间长度为例，可以直接取Xt＝X10。则第t+m期的预测参数为：其中m为预测步长，at、bt、ct为预测系数，at、bt、ct的计算公式为：在选取初值时，由于数据量仅为十个，所以取初值为前十个周期的平均数，即S1(1)＝S1(2)＝S1(3)＝Xav。选取预测步长m值为1，每次通过前十个周期预测出下一周期时间长度为令其中Ts为采样时间间隔，根据Ts进行下一周期的采样，根据下一周期采样得出的数据，加上前九个周期采样得到的数据，进行计算，通过离散傅里叶变换(DFT)，得出谐波、间谐波测量值。a5：根据公式Ts＝T/N，N为每周期的采样点数，计算出采样时间间隔，进行下一周期采样，并储存采样值；a6：判断已采样周期数是否达到十个，若是，进入步骤a7，若否，进入步骤a8；a7：去掉最早一周期数据，加入最新一周期的采样值，通过离散傅里叶变换方法计算谐波参数，返回步骤a8。步骤a7中，采用离散傅里叶变换(DFT)的方法来计算谐波参数，相比较于其它方法而言，离散傅里叶方法更为实用，有效并且较为稳定，并且可以通过快速傅里叶变换本文档来自技高网...

【技术保护点】
一种新型电力谐波参数软件同步采样方法，其特征在于，包括如下步骤：步骤a1：设置每周期采样点数N，N为每周期的采样点数；步骤a2：测出前W个周期时间长度Xi，并存储，其中i＝1、2…W；步骤a3：根据测量的前W个周期时间长度，计算平滑系数α，并储存α值；步骤a4：根据三次指数平滑模型，预测下一采样周期的时间长度；步骤a5：根据公式Ts＝T/N，计算出采样时间间隔，进行下一周期采样，并储存采样值；步骤a6：判断已采样周期数是否达到W个，若是，进入步骤a7，若否，进入步骤a8；步骤a7：去掉最早一周期数据，加入最新一周期的采样值，通过离散傅里叶变换方法计算谐波参数，返回步骤a8；步骤a8：去掉最早一周期数据，通过过零点检测，存储最新一周期的周期时间长度，返回步骤a3。

【技术特征摘要】
1.一种新型电力谐波参数软件同步采样方法，其特征在于，包括如下步骤：步骤a1：设置每周期采样点数N，N为每周期的采样点数；步骤a2：测出前W个周期时间长度Xi，并存储，其中i＝1、2…W；步骤a3：根据测量的前W个周期时间长度，计算平滑系数α，并储存α值；步骤a4：根据三次指数平滑模型，预测下一采样周期的时间长度；步骤a5：根据公式Ts＝T/N，计算出采样时间间隔，进行下一周期采样，并储存采样值；步骤a6：判断已采样周期数是否达到W个，若是，进入步骤a7，若否，进入步骤a8；步骤a7：去掉最早一周期数据，加入最新一周期的采样值，通过离散傅里叶变换方法计算谐波参数，返回步骤a8；步骤a8：去掉最早一周期数据，通过过零点检测，存储最新一周期的周期时间长度，返回步骤a3。2.根据权利要求1所述的一种新型电力谐波参数软件同步采样方法，其特征在于，在步骤a1中，设置一个测量周期内的采样点数为N，其中保证N>2M，M为被测信号最高谐波次数。3.根据权利要求1所述的一种新型电力谐波参数软件同步采样方法，其特征在于，所述步骤a2的检测方式是通过过零点检测进行。4.根据权利要求1所述的一种新型电力谐波参数软件同步采样方法，其特征在于，在步骤a3中，当误差平方和最小时，此时的平滑系数α为最准确，最小误差平方公式其中K表示预测的次数，表示第t次预测值，Xt表示第t次实际值。5.根据权利要求1所述的一种新型电力谐波参数软件同步采样方法，其特征在于，在步骤a4中，三次指...

【专利技术属性】
技术研发人员：佟为明，金显吉，李中伟，林景波，刘勇，赵志衡，陶耕宇，
申请(专利权)人：哈尔滨工业大学，
类型：发明
国别省市：黑龙江,23