CN103575979A - 一种数字化测量交流电频率的方法 - Google Patents

Abstract

本发明公开了一种数字化测量交流电频率的方法，采用数字积分并插值的方式处理连续的数字采样信号，利用信号过正负峰值后的积分开始点P₁点和根据插值的方式获取的积分结束点k的虚拟发生时间计算电信号的零点，有一系列电信号的过零点计算电信号的频率或周期则可以由P₁的采样发生时间和k的虚拟发生时间计算电信号的零点，由一系列电信号的过零点计算电信号的频率或周期。本发明不需要将低频信号转换成方波，尤其适合于低频正弦信号的频率测量，***检测方便，硬件开销小；***的抗干扰能力较强；实现起来简单合理，数据准确、可靠。

Description

一种数字化测量交流电频率的方法

技术领域

本发明涉及一种数字化测量交流电频率的方法，也可以用于固定频率或频率变化不大的低频正弦信号、三角波信号、锯齿波信号等等的频率测量。

背景技术

在日常生活、生产过程中常常需要监测信号的频率。在现有的低频信号的频率测量方法中，多是将低频信号先通过整形电路形成方波，检测方波相邻两个上升沿或者下降沿的时间间隔T，求倒数得出频率，但是这种方法的抗干扰能力较差。除此之外还有一些软件测量频率的方法，比如：零交法、解析法、误差最小化原理类算法、DFT类算法、正交去调制法等等。这些算法有些易于实现，但由于精度较低而不能投入实用；有些算法的测量效果不错，但由于运算量过大也不能投入实用。业界需要的是一种可以实用且测量精度较准确的方法。

发明内容

本发明所要解决的技术问题是针对现有技术中低频率电信号的频率测量精度不够准确或者是测量过程中运算量过大的缺点，提出一种数字化测量交流电频率的方法，该方法也可以用于固定频率或频率变化不大的低频正弦信号、三角波信号、锯齿波信号等等的频率测量。

本发明为解决上述技术问题采用以下技术方案：

该方法的原理是采用数字积分并插值的方式处理连续的数字采样信号，在信号过正负峰值之后选取一个采样点P₁作为积分开始点，在后续的采样点中会存在这样连续的两个点P_i和P_i+1，如果从P₁到P_i的数字积分数值和从P₁到P_i+1的数字积分数值的乘积小于等于零时,则可以在P_i和P_i+1之间通过插值的方式获取一个积分结束点k，k的获取条件是从P₁到k的数字积分为0。则可以由P₁的采样发生时间和k的虚拟发生时间计算电信号的零点，由一系列电信号的过零点计算电信号的频率或周期；具体步骤如下：

步骤1），对被测电信号进行采样，得到被测电信号的采样值，然后生成采样信号；

步骤2）,在被测电信号的采样信号中，当被测电信号的周波越过任意一个峰值后选择m个在采样时间上按次序排放的采样点，设这m个采样点的采样时间分别为t₁、t₂、...t_i、t_i+1、...t_m，采样值分别为x₁、x₂、...x_i、x_i+1...x_m，其中i、m均为自然数且1≤i<m；

步骤3），令S_i为从第1采样点(x₁,t₁)到第i采样点(x_i,t_i)的数字积分，则S_i+1为从第1采样点(x₁,t₁)到第i+1采样点(x_i+1,t_i+1)的数字积分；当S_i和S_i+1的乘积小于等于0时，在第i个采样点和第i+1个采样点之间通过插值的方式获取一个数字积分结束点k，使得从第1采样点(x₁,t₁)到积分结束点(x_k,t_k)的数字积分为零，计算得出积分结束点k的虚拟发生时间t_k；

步骤4），计算过零点时间：

步骤5），重复步骤2）至步骤4）得到若干个过零点时间：T_z1,T_z2,T_z3...T_zj,T_zj+1，T_zj+2...T_zn，其中，j=1、2、…n，则被测交流电信号的周期：T_pj=T_zj+2-T_zj或者T_pj=2×(T_zj+1-T_zj)，被测交流电信号的频率： $F_j=\frac{1}{T_{\mathrm{pj}}}=\frac{1}{T_{\mathrm{zj}+2}-T_{\mathrm{zj}}}$ 或者： $F_j=\frac{1}{T_{\mathrm{pj}}}=\frac{1}{2\&times;(T_{\mathrm{zj}+1}-T_{\mathrm{zj}})}.$

进一步的，本发明的一种数字化测量交流电频率的方法，步骤1）所述对被测电信号进行采样是等时间间隔采样或者是不等时间间隔采样。

进一步的，本发明的一种数字化测量交流电频率的方法，步骤1）所述的采样为对被测电信号的整周波进行采样。

进一步的，本发明的一种数字化测量交流电频率的方法，所述第1采样点的选取方案为：当被测电信号的周波过任意一个峰值后未到达零点时选择一个采样点作为第1采样点。

进一步的，本发明的一种数字化测量交流电频率的方法，还包括以下步骤：重复步骤1）到步骤5），测量出被测电信号的若干个周期或频率值，然后对得到的若干个周期或频率值采用平均法计算出最终频率。

进一步的，本发明的一种数字化测量交流电频率的方法，所述数字积分方式包括：梯形积分或者矩形积分；所述插值方式为：矩形插值或者梯形插值。

本发明采用以上技术方案与现有技术相比，具有以下技术效果：

1）传统的零交法采用符号相反的两个连续点来确定过零点，虽然算法物理概念清晰，但是容易受谐波、噪声等的干扰，测量精度低。针对电信号大多数是对称的特点，对被测电信号过峰值后进行数字积分，通过插值的方式获取一个积分结束点，使得从积分开始点到积分结束点的数字积分为零，由积分开始点的采样发生时间和积分结束点的虚拟发生时间计算电信号的零点。确定了信号的过零点之后，可以计算出电信号的频率和周期。相比较传统的零交法而言，***的抗干扰能力有极大的提高，频率测量的准确度也大大提高。实验表明即使电信号中混入了较大的谐波、白噪声，测量的精确度也非常之高。

2）该方法在进行计算时非常方便简捷，适合嵌入式***下使用。

附图说明

图1是采用矩形积分、矩形插值的方式计算得到积分结束点的虚拟发生时间，并获取过零点的示意图；

图2是采用梯形积分、梯形插值的方式计算得到积分结束点的虚拟发生时间，并获取过零点的示意图；

图3是采用梯形积分、矩形插值的方式计算得到积分结束点的虚拟发生时间，并获取过零点的示意图；

图4是对交流电信号进行采样后，在每一个峰值后进行矩形积分、矩形插值，获得过零点并计算周期的示意图；

图5是对交流电信号进行采样后，在每一个峰值后进行梯形积分、梯形插值，获得过零点并计算周期的示意图；

图6是对交流电信号进行采样后，在每一个峰值后进行梯形积分、矩形插值，获得过零点并计算周期的示意图。

具体实施方式

下面结合附图对本发明的技术方案做进一步的详细说明：

本发明的实质是采用数字积分并插值的方式处理连续的数字采样信号，在信号过正负峰值之后选取一个采样点P₁作为积分开始点，在后续的采样点中会存在这样连续的两个点P_i和P_i+1，如果从P₁到P_i的数字积分数值和从P₁到P_i+1的数字积分数值的乘积小于等于0时,则可以在P_i和P_i+1之间通过插值的方式获取一个积分结束点k，k的获取条件是从P₁到k的数字积分为0。则可以由P₁的采样发生时间和k的虚拟发生时间计算电信号的零点，由一系列电信号的过零点计算电信号的频率或周期。这里所谓的数字积分有如下几种方式：梯形积分、矩形积分方式；这里所谓的插值有矩形插值、梯形插值方式。

本发明的原理通过测量交流电信号的两个周波的频率来进行说明。

本发明的具体实施过程如下：

1.对被测电信号进行采样,这里所述的采样为对整周波进行的采样。可以是等时间间隔采样，也可以是不等时间间隔采样。

2.在被测电信号过正负峰值后选择若干个在采样时间上按次序排放的采样点，假设每个采样点的采样时间为t_i，采样值为x_i，i≥1，令S_i为从点(x₁,t₁)到第i点的数字积分；

3.如果S_i和S_i+1的乘积小于等于0，则可在第i点和第i+1点之间通过插值的方式获取一个积分结束点k，使得从点(x₁,t₁)到点k的数字积分为零，可计算得出点k的虚拟发生时间t_k；

4.所选择的峰值之后的过零点时间的计算公式如下：

5.计算过零点时间的具体操作可以参见图1、图2、图3。图1是采用矩形积分、矩形插值的方式计算得到积分结束点的虚拟发生时间，并获取过零点的示意图。图1中从t₁开始数字积分，t_k是插值获取的积分结束点的虚拟发生时间，图中从t₁到t_k的数字积分为0。T_z是计算获取的过零点。图2是采用梯形积分、梯形插值的方式计算得到积分结束点的虚拟发生时间，并获取过零点的示意图。图2中从t₁开始数字积分，t_k是插值获取的积分结束点的虚拟发生时间，图中从t₁到t_k的数字积分为0。T_z是计算获取的过零点。图3是采用梯形积分、矩形插值的方式计算得到积分结束点的虚拟发生时间，并获取过零点的示意图。图3中从t₁开始数字积分，t_k是插值获取的积分结束点的虚拟发生时间，图中从t₁到t_k的数字积分为0。T_z是计算获取的过零点。如果纯粹采用梯形积分法，在进行插值计算时存在计算量偏大的情况。图3在开始积分的时候采用梯形积分法，但是在进行插值运算时，采用了矩形插值的方法，这样做的目的是降低运算的复杂度。为取得较为准确的测量结果,在附图的计算中设该矩形的高度是邻近的两个采样点的平均值。

6.重复上述步骤，得到若干个过零点时间，T_z1,T_z2,T_z3...T_zj,T_zj+1，T_zj+2...T_zn，则被测交流电信号的周期：T_pj=T_zj+2-T_zj或者T_pj=2×(T_zj+1-T_zj)，被测交流电信号的频率： $F_j=\frac{1}{T_{\mathrm{pj}}}=\frac{1}{T_{\mathrm{zj}+2}-T_{\mathrm{zj}}}$ 或者： $F_j=\frac{1}{T_{\mathrm{pj}}}=\frac{1}{2\&times;(T_{\mathrm{zj}+1}-T_{\mathrm{zj}})}.$

7.对电信号进行采样后在每一个零点附近进行计算并测量频率的图形如图4、图5、图6。图4采用的是矩形积分、矩形插值的方法，图5采用的是梯形积分、梯形插值的方法，图6采用的是梯形积分、矩形插值的方法。图4、图5、图6中的交流电信号存在4个过零点，因此进行了四次积分并计算。图4、图5、图6中水平轴上的空心圆点是计算得到的过零点，过零点依次是：T_z1，T_z2，T_z3，T_z4。则该信号的周期可以表示为T_p=T_z3-T_z1或者T_p=T_z4-T_z2或者T_p=2×(T_z2-T_z1)或者T_p=2×(T_z3-T₂)或者T_p=2×(T_z4-T_z3)，频率F_p则是T_p的倒数。

8.重复上述步骤，测量出被测电信号的若干个周期或频率值，然后对得到的若干个周期或频率值采用平均法计算出最终频率。

从附图中可以看到，采用梯形积分、梯形插值法所获得的零点和实际信号的零点最为接近，采用梯形积分、矩形插值法所获得的零点和实际信号的零点有一定的误差，采用矩形积分、矩形插值法所获得的零点和实际信号的零点的误差最大。当然随着采样点数的增加，这种误差将会越来越小。考虑到嵌入式***的特点，优选方案是采用梯形积分，矩形插值法计算信号的过零点。

综上所述，本发明所涉及的一种数字化测量交流电频率的方法，不需要将交流电信号转换成方波，尤其适合于低频正弦信号的频率测量，***检测方便，硬件开销小；***的抗干扰能力很强；实现起来简单合理，数据准确、可靠；有较广泛的应用前景。

Claims (6)

1.一种数字化测量交流电频率的方法，其特征在于：包括以下步骤：

步骤1），对被测电信号进行采样，得到被测电信号的采样值，然后生成采样信号；

步骤2）,在被测电信号的采样信号中，当被测电信号的周波越过任意一个峰值后选择m个在采样时间上按次序排放的采样点，设这m个采样点的采样时间分别为t₁、t₂、...t_i、t_i+1、...t_m，采样值分别为x₁、x₂、...x_i、x_i+1...x_m，其中i、m均为自然数且1≤i<m；

步骤3），令S_i为从第1采样点(x₁,t₁)到第i采样点(x_i,t_i)的数字积分，则S_i+1为从第1采样点(x₁,t₁)到第i+1采样点(x_i+1,t_i+1)的数字积分；当S_i和S_i+1的乘积小于等于0时，在第i个采样点和第i+1个采样点之间通过插值的方式获取一个数字积分结束点k，使得从第1采样点(x₁,t₁)到积分结束点(x_k,t_k)的数字积分为零，计算得出积分结束点k的虚拟发生时间t_k；

步骤4），计算过零点时间：

步骤5），重复步骤2）至步骤4）依次得到若干个过零点时间：T_z1,T_z2,T_z3...T_zj,T_zj+1，T_zj+2...T_zn，其中，j=1、2、…n，则被测交流电信号的周期：T_pj=T_zj+2-T_zj或者T_pj=2×(T_zj+1-T_zj)，被测交流电信号的频率： $F_j=\frac{1}{T_{\mathrm{pj}}}=\frac{1}{T_{\mathrm{zj}+2}-T_{\mathrm{zj}}}$ 或者： $F_j=\frac{1}{T_{\mathrm{pj}}}=\frac{1}{2\&times;(T_{\mathrm{zj}+1}-T_{\mathrm{zj}})}.$

2.根据权利要求1所述的一种数字化测量交流电频率的方法，其特征在于：步骤1）所述对被测电信号进行采样是等时间间隔采样或者是不等时间间隔采样。

3.根据权利要求1所述的一种数字化测量交流电频率的方法，其特征在于：步骤1）所述的采样为对被测电信号的整周波进行采样。

4.根据权利要求1所述的一种数字化测量交流电频率的方法，其特征在于：所述第1采样点的选取方案为：当被测电信号的周波过任意一个峰值后未到达零点时选择一个采样点作为第1采样点。

5.根据权利要求1至4中任一项所述的一种数字化测量交流电频率的方法，其特征在于：还包括以下步骤：重复步骤1）到步骤5），测量出被测电信号的若干个周期或频率值，然后对得到的若干个周期或频率值采用平均法计算出最终频率。

6.根据权利要求1至4中任一项所述的一种数字化测量交流电频率的方法，其特征在于：所述数字积分方式包括：梯形积分或者矩形积分；所述插值方式包括：矩形插值或者梯形插值。

Images