CN102033161B - 一种交流电信号的频率测量方法 - Google Patents

Abstract

本发明提供了一种交流电信号的频率测量方法。该方法包括以下步骤：对被测电信号进行采样；根据采样值，在信号零点附近选择若干个采样点，计算若干个采样点的采样值的平均值

和平均采样发生时间；由当前信号的特点和

经查表获得一个时间

，由

来表示该信号的过零点；根据过零点，计算被测信号的周期或频率。本发明不需要将交流电信号转换成方波，尤其适合于固定频率的低频正弦信号的频率测量，系统检测方便，硬件开销小；系统具有一定的抗干扰能力；实现简单合理，数据准确、可靠。

Description

一种交流电信号的频率测量方法

技术领域

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

背景技术

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

发明内容

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

为解决上述技术问题，本发明提出一种交流电信号的频率测量方法，包括以下步骤：

步骤1)，对被测电信号进行采样，得到被测电信号的采样值；在电信号的采样值的某一个零点附近选择任意组合的n个采样点；

步骤2)，假设每个采样点的采样时间为t_i，采样值为y_i，则n个采样点的采样值的平均值M和平均采样发生时间T表示为：

$M_j=\frac{\underset{1}{\overset{n}{\mathrm{\Σ}}}{y}_i}{n},$ $T_j=\frac{\underset{1}{\overset{n}{\mathrm{\Σ}}}{t}_i}{n};$

上式中，n≥1，n、i、j均为自然数，则(M_j，T_j)表示一个虚拟点；

步骤3)，根据被测电信号在过零点附近的交流电信号是上升或下降的过程、电信号有效值或者峰值，以及步骤2)所述n个采样点的采样值的平均值M_j，通过查表的方式来获得一个时间差T′_j，所述时间差T′_j代表纯正弦信号下，离虚拟点(M_j，T_j)最近的过零点与该虚拟点之间的时间差；对应的过零点时间为：T_zj＝T_j+T′_j；

其中上述查表所涉及的表格共有两列，其中第一列是

或者表格的第二列是

对应的时间差T′_j。

步骤4)，依次在电信号采样值的其它零点附近选择任意组合的n个采样点，重复步骤2)和步骤3)得到若干个过零点时间，根据这若干个过零点时间，计算得出被测电信号的周期或频率。

进一步的，前述交流电信号的频率测量方法中，步骤1)所述对电信号进行采样是等时间间隔采样或者是不等时间间隔采样。

进一步的，前述交流电信号的频率测量方法中，步骤1)所述的采样为对电信号的整周波进行采样。

进一步的，前述交流电信号的频率测量方法中，步骤4)中所述的计算被测电信号的周期或频率的具体步骤如下：重复步骤2)、步骤3)得到若干个过零点时间T_z1，T_z2，T_z3...T_zj，T_zj+1，T_zj+2...T_zn...，则被测交流电信号的周期：T_pj＝T_zj+2-T_zj，被测交流电信号的频率：

进一步的，前述交流电信号的频率测量方法中，在步骤2)中还包括对采样点的平均值M进行修正的步骤，所述修正的步骤如下：

a)设由计算得到的交流电信号前几个周波的周期依次为T_pk-2，T_pk-2，T_pk-1，T_pk，当前几个周波的周期不存在或得不到时，令其等于交流电标准周期；

b)预测的当前周波的周期表示为：

T_p＝K₁×T_pk+K₂×T_pk-1+K₃×T_pk-2+K₄×T_pk-3；其中K₁、K₂、K₃、K₄为任意常数，但须满足条件K₁+K₂+K₃+K₄＝1；

c)M_修正值＝M×T_p÷交流电标准周期。

进一步的，前述交流电信号的频率测量方法中，还包括以下步骤：重复步骤1)到步骤4)，测量出被测电信号的若干个周期或频率值，然后对得到的若干个周期或频率值采用平均法计算出最终频率。

本发明采用以上技术方案具有如下有益效果：

1)传统的零交法采用符号相反的两个连续点来确定过零点，虽然算法物理概念清晰，但是容易受谐波、测量误差等的干扰，测量精度低。随着数据采集速度的提高，可以在零点附近任意选择若干个采样点。由选取的若干个采样点可经过计算、查表，进而得出过零点。确定了信号的过零点之后，可以计算出交流信号的频率和周期。相比较传统的零交法而言，系统的抗干扰能力有所提高。

2)该方法在进行计算时可选取若干个采样点来参加计算，这能在一定程度上降低系统的随机误差；该方法的另外一个优点是算法的计算量和选取的采样点的个数只呈现线形关系，在选取的采样点数增加时，运算的精度可以大大提高但运算量不会呈现指数增长。

3)本发明还可以实现有限的频率预测，具体可以参见图3，图3中的四个采样点来源于图2中的S0，S1，S2，S3，且都是正电平，在S3采样结束之后(此时信号还没有到达零点)，可以根据S0，S1，S2，S3四个采样点计算查表后获取信号的下一个过零点Z₀，因此采用此种方法可以实现有限的频率预测。

附图说明：

图1是正弦信号有从低向高过零点和从高向低过零点的两种可能的示意图。

图2是在零点附近进行等时间间隔采样而获取的8个采样点的示意图(采样点的电平由高到低)。

图3是由4个正电平采样点进行计算后查表所获得的信号过零点的示意图。

图4是由3个正电平和1个负电平采样点进行计算后查表所获得的信号过零点的示意图。

图5是由2个正电平和2个负电平采样点进行计算后查表所获得的信号过零点的示意图。

图6是在零点附近进行等时间间隔进行采样而获取的8个采样点的示意图(采样点的电平由低到高)。

图7是在不等时间间隔采样情形下，由3个正电平和1个负电平采样点进行计算后查表所获得的信号过零点的示意图。

图8是对交流电信号进行采样后在每一个零点附近进行计算并查表后获得过零点的示意图。

具体实施方式

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

本发明的实质是在零点附近选择若干个采样点，计算若干个采样点的平均值M和平均采样发生时间T，然后根据当前信号的特点和平均值M查表获取T′，由T+T′来表示信号的过零点，从而计算出信号的频率或周期。选取的采样点不宜过多，选取的采样点如果较多的话，计算复杂度会高一些；选取得采样点也不宜过少，选取的采样点如果较少的话，容易受到干扰的影响。因此用户可以根据需要选取合理的采样点组合和个数来进行计算。

为获得较为精确的测量结果，建议的采样点的选取方案是：正的采样点的数量和负的采样点的数量相等或者大致相等，采样点的个数则可根据具体实际情况来确定。

本发明的原理通过测量交流电信号的两个周波的频率来进行说明，随着采样速度的提高，在每个周波内的采样点也越来越多，当采样次数足够多时，可以在零点附近选取若干个采样点，对采样点计算并查表后可获取交流电信号的过零点。后续的具体实施方式中都采用了四个采样点来计算过零点，实际应用中可以根据需要适当增加或者减少采样点的个数。

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

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

2.在信号的每一个零点附近选择任意组合的若干个采样点。选出的采样点的采样值可以全部为正；也可以是全部为负；也可以根据需要选择若干个采样值为正、若干个采样值为负的采样点。

3.得到过零点的过程如下：假设选取的采样点数为n，每个采样点的采样时间为t_i，采样值为y_i，则n个采样点的采样值的平均值M和平均采样发生时间T可表示为：

$M_j=\frac{\underset{1}{\overset{n}{\mathrm{\Σ}}}{y}_i}{n},$ $T_j=\frac{\underset{1}{\overset{n}{\mathrm{\Σ}}}{t}_i}{n};$

上式中，n≥1，i、j均为自然数，(M_j，T_j)可看作是一个虚拟点；

获取(M_j，T_j)后，可由当前信号的特点、信号的有效值(或者是峰值)和M_j通过查表来获得T′_j，T′_j的物理意义是纯正弦信号下，离虚拟点(M_j，T_j)最近的过零点与该虚拟点之间的时间差。过零点附近的交流电信号可能是上升的(图1中的2)，也有可能是下降的(图1中的1)。针对过零点附近是上升的交流电信号(频率为50Hz)设计一个表格，表格的第一列是

表格的第二列是

对应的T′_j。通过该表格，可以根据

的数值来获取对应的T′_j，具体的表格如下(表格的第一列也可以设计成

)：

针对过零点附近是下降的交流电信号(频率为50Hz)设计一个表格，表格的第一列是

表格的第二列是

对应的T′_j。通过该表格，可以根据

的数值来获取对应的T′_j，具体的表格如下(表格的第一列也可以设计成

)：

上述两张表格的粒度还可以设计的更细一些，这样可以提高T′_j的准确性。需要注意的是，上述表格的

的范围是[-0.5，0.5]，这个范围可以根据应用的需要做适当的增加或者减少。如果需要测量其他低频率信号的频率时，需要重新设计上述表格。

查表后，则对应的过零点时间可表示为：T_zj＝T_j+T′_j。

4.在被测信号频率偏离交流电标准频率时，可以对采样点的平均值M进行修正来提高测量的准确性。修正的方法如下：设由计算得到的交流电信号前几个周波的周期依次为T_pk-3，T_pk-2，T_pk-1，T_pk(在前面周波的周期没有计算得到时，可令部分或者全部的T_pk-3，T_pk-2，T_pk-1，T_pk等于交流电标准周期)；则预测的当前周波的周期可表示为：T_p＝K₁×T_pk+K₂×T_pk-1+K₃×T_pk-2+K₄×T_pk-3；公式中的K₁，K₂，K₃，K₄可为任意常数，但须满足条件K₁+K₂+K₃+K₄＝1。则令M_修正值＝M×T_p÷交流电标准周期。在被测信号频率偏移交流电标准频率时，采用M_修正值来查表能得到更加准确的测量结果。在后续的实验中，令K₁＝1，K₂＝K₃＝K₄＝0，则T_p＝T_pk，即认为当前周波的周期与上一周波的周期相同。

5.具体操作可以参见图2、图3、图4和图5。图2上的S0，S1，S2，S3，S4，S5，S6，S7是在零点附近进行等时间间隔采样而获取的8个采样点。图3中的空心三角形所在的点是由图2中的采样点S0，S1，S2，S3进行计算后获得的虚拟点(M₀，T₀)，Z₀是经查表后获得的过零点；图4中的空心三角形所在的点是由图2中的采样点S1，S2，S3，S4进行计算后获得的虚拟点(M₁，T₁)，Z₁是经查表后获得的过零点；图5中的空心三角形所在的点是由图2中的采样点S2，S3，S4，S5进行计算后获得的虚拟点(M₂，T₂)，Z₂是经查表后获得的过零点。图3、图4、图5说明了信号由高到低穿越零点时的情形。当信号由低向高穿越零点时(此种情况可以参见图6)计算的方法相同。图3、图4、图5的例子是等时间间隔采样，图7是不等时间间隔采样情形下经过计算获取虚拟点(M₃，T₃)，Z₃是经查表后获得的过零点。

6.重复上述步骤，可以得到若干个过零点时间T₁，T₂，T₃...T_j，T_j+1，T_j+2...T_n...，则被测交流电信号的周期：T_pj＝T_j+2-T_j，被测交流电信号的频率：

对交流电信号进行采样后在每一个零点附近进行计算并查表的图形如图8。图8中的交流电信号存在4个过零点，因此进行了四次计算并查表。计算时采用的是零点附近的2个正电平、2个负电平采样点，图中空心三角形所在的点为计算得到的虚拟点，图中的空心圆点是得到的过零点，过零点依次是：T_i-1，T_i，T_i+1，T_i+2。则该信号的周期可以表示为T_p＝T_i+1-T_i-1或者T_p＝T_i+2-T_i。频率F_p则是T_p的倒数。

优选方案是步骤1)所述的采样为等时间间隔采样。

采用本方法可以在一个信号周期内还可以进行多次频率测量。

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

Claims (5)

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

步骤1），对被测交流电信号进行采样，得到被测交流电信号的采样值；在被测交流电信号的采样值的某一个过零点附近选择任意组合的n个采样点；

步骤2），假设每个采样点的采样时间为t_i，采样值为y_i，则n个采样点的采样值的平均值M_j和平均采样发生时间T_j表示为：

$M_j=\frac{\underset{1}{\overset{n}{\mathrm{\Σ}}}{y}_i}{n},$ $T_j=\frac{\underset{1}{\overset{n}{\mathrm{\Σ}}}{t}_i}{n};$

上式中，n≥1，n、i、j均为自然数，则（M_j，T_j）表示一个虚拟点；

对n个采样点的采样值的平均值M_j进行修正，所述修正的步骤如下：

a）设由计算得到的被测交流电信号前几个周波的周期依次为T_pk-3,T_pk-2，T_pk-1,T_pk，当前几个周波的周期不存在或得不到时，令其等于交流电标准周期；

b）预测的当前周波的周期表示为：

T_p=K₁×T_pk+K₂×T_pk-1+K₃×T_pk-2+K₄×T_pk-3；其中K₁、K₂、K₃、K₄为任意常数，但须满足条件K₁+K₂+K₃+K₄=1；

c）M_修正值=M_j×T_p÷交流电标准周期；

步骤3），根据被测交流电信号在过零点附近的交流电信号是上升或下降的过程、被测交流电信号有效值或者峰值，以及步骤2）所述n个采样点的采样值的平均值的修正值M_修正值，通过查表的方式来获得一个时间差

所述时间差

代表纯正弦信号下，离虚拟点（M_修正值，T_j）最近的过零点与该虚拟点之间的时间差；对应的过零点时间为： $T_{\mathrm{zj}}=T_j+T_j^{\′};$

其中上述查表所涉及的表格共有两列，其中第一列是

或者

表格的第二列是

对应的时间差

步骤4），依次在被测交流电信号采样值的其它过零点附近选择任意组合的n个采样点，重复步骤2）和步骤3）得到若干个过零点时间，根据这若干个过零点时间，计算得出被测交流电信号的周期或频率。

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

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

4.根据权利要求1所述的交流电信号的频率测量方法，其特征在于：步骤4）中所述的计算被测交流电信号的周期或频率的具体步骤如下：重复步骤2）、步骤3）得到若干个过零点时间T_z1,T_z2,T_z3...T_zj,T_zj+1,T_zj+2...T_zn...，则被测交流电信号的周期：T_pj=T_zj+2-T_zj，被测交流电信号的频率：

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

Images