CN102778606A - 一种交流电的瞬时相位的数字化测量方法 - Google Patents

一种交流电的瞬时相位的数字化测量方法 Download PDF

Info

Publication number
CN102778606A
CN102778606A CN2012102610173A CN201210261017A CN102778606A CN 102778606 A CN102778606 A CN 102778606A CN 2012102610173 A CN2012102610173 A CN 2012102610173A CN 201210261017 A CN201210261017 A CN 201210261017A CN 102778606 A CN102778606 A CN 102778606A
Authority
CN
China
Prior art keywords
sampled
cycle
instantaneous phase
sampled point
time
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2012102610173A
Other languages
English (en)
Other versions
CN102778606B (zh
Inventor
戴尔晗
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Post and Telecommunication University
Nanjing University of Posts and Telecommunications
Original Assignee
Nanjing Post and Telecommunication University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing Post and Telecommunication University filed Critical Nanjing Post and Telecommunication University
Priority to CN201210261017.3A priority Critical patent/CN102778606B/zh
Publication of CN102778606A publication Critical patent/CN102778606A/zh
Application granted granted Critical
Publication of CN102778606B publication Critical patent/CN102778606B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明公开了一种交流电的瞬时相位的数字化测量方法,针对连续的数字采样信号,在信号过负峰值后选取一个采样点P1,在后续采样点中选取两个连续点,使P1采样值的负值P4在两点之间,利用插值法选取P4。利用P1的发生时间和P4的虚拟发生时间计算电信号零点,由一系列零点求得电信号周期。根据前一时段信号过零点时间、预测的当前周波周期以及时间,得到当前瞬时相位。本发明不需要将交流电转换成方波,克服了抗干扰能力差的缺点,尤其适合于交流电的瞬时相位测量,系统检测方便,硬件开销小;同时实现起来简单合理,数据准确、可靠;在测量的精度、抗干扰性等方面得到很大的提高。

Description

一种交流电的瞬时相位的数字化测量方法
技术领域
本发明涉及一种交流电的瞬时相位的数字化测量方法,也适用于低频率信号的瞬时相位测量。
背景技术
在日常生活、生产过程中常常需要测量交流电的瞬时相位,测量的瞬时相位越准确,应用的效果往往更好。在现有的瞬时相位测量方法中,大多数采用传统的零交法测量某一路信号的瞬时相位,这种算法易于实现,但是容易受到干扰的影响,使得测量的精度较低。
发明内容
本发明所要解决的技术问题是针对现有技术中交流电的瞬时相位测量精度不够准确,提出了一种交流电的瞬时相位的数字化测量方法,相比较传统的零交法而言,在测量的精度、抗干扰性等发明得到很大的提高。
本发明为解决上述技术问题采用以下技术方案:
一种交流电的瞬时相位的数字化测量方法,包括以下步骤:
步骤1),对被测电信号进行采样,得到被测电信号的采样值,然后生成采样信号;
步骤2),在被测电信号的采样信号中,当被测电信号的周波过负峰值后未到达零点时选择一个采样点作为第一采样点P1i,设该第一采样点P1i的采样时间为T1i,采样值为S1i
步骤3),当被测电信号的周波越过零点后未到达正峰值时,在该段采样信号中选取连续的两个采样点,即第二采样点P2i和第三采样点P3i,设这两个点的采样时间分别是T2i和T3i,采样值分别为S2i和S3i,其中采样值S2i和S3i需满足|S2i|≤|S1i|≤|S3i|且|S2i|<|S3i|;
步骤4),将第一采样点P1i采样值的相反值所对应的点作为虚拟采样点P4i,计算虚拟采样点P4i的虚拟发生时间 T 4 i = T 2 i + ( T 3 i - T 2 i ) × ( S 1 i + S 2 i ) S 2 i - S 3 i ;
步骤5),计算过零点时间:
Figure BDA00001935331200012
步骤6),重复步骤2)至步骤5)得到若干个过零点时间:Tz1,Tz2,Tz3...Tzi,Tzi+1,Tzi+2...Tzn,其中,i=1、2、…n,n为正整数;则被测电信号的周期:Tpi=Tzi+1-Tzi,因此获得的信号的周期依次为Tp1,Tp2,Tp3,...Tpi-1,Tpi,...Tpk,其中k为正整数且k=n-1;
步骤7),根据步骤6)计算得出的被测电信号的周期,预测当前周波的周期;
步骤8),根据前一个时段交流电信号的过零点的时间、步骤7)预测的当前周波的周期以及当前时间,计算得出当前的瞬时相位。
进一步的,本发明的一种交流电的瞬时相位的数字化测量方法,步骤1)所述的采样是等时间间隔采样或不等时间间隔采样。
进一步的,本发明的一种交流电的瞬时相位的数字化测量方法,步骤1)所述的采样为对信号的整周波进行采样。
进一步的,本发明的一种交流电的瞬时相位的数字化测量方法,所述步骤7)的具体步骤如下:
设由步骤6)计算得到的被测电信号前几个周波的周期依次为Tpk-3,Tpk-2,Tpk-1,Tpk;则预测的当前周波的周期为:Tp=K1×Tpk+K2×Tpk-1+K3×Tpk-2+K4×Tpk-3;其中,K1,K2,K3,K4为任意常数,但须满足条件K1+K2+K3+K4=1,k为大于3的正整数。
进一步的,本发明的一种交流电的瞬时相位的数字化测量方法,所述步骤8)的具体步骤如下:设被测电信号的前一个时段过零点的时间为Tz,当前时间为Tc,由步骤7)预测得到的被测电信号当前周波的周期为Tp,则被测电信号当前的瞬时相位以度数值表示为:
Figure BDA00001935331200021
或者以弧度值表示为:
Figure BDA00001935331200022
进一步的,本发明的一种交流电的瞬时相位的数字化测量方法,在所述步骤3)中,当选取的第一至第三采样点的采样值符合|S2i|=|S1i|=|S3i|条件时,采取以下两种处理方案中的任意一种:
a),重新选取第二采样点P2i和第三采样点P3i,使得第二采样点的采样值S2i和第三采样点的采样值S3i须满足|S2i|<|S3i|,然后继续步骤4)至步骤8);
b),在第二采样点P2i与第三采样点P3i的采样时间区间[T2i,T3i]上任取一点作为虚拟采样点P4i所对应的虚拟采样时间T4i,然后继续步骤5)至步骤8)。
进一步的,本发明的一种交流电的瞬时相位的数字化测量方法,在所述方案b)中,取第二采样点P2i与第三采样点P3i的采样时间区间[T2i,T3i]上的中点,计算虚拟采样点P4i的虚拟发生时间
Figure BDA00001935331200031
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
1)传统的零交法采用符号相反的两个连续点来确定零点并继而确定交流电的瞬时相位,虽然算法物理概念清晰,但是容易受谐波、测量误差等的干扰,测量精度低。只有准确定位零点,才能计算出精确的瞬时相位。本发明由选取的3个采样点可经过计算得出过零点。确定了信号的过零点之后,可以计算出交流信号的频率和周期。最后根据前一个时段交流电信号的过零点的时间、预测的当前周波的周期以及当前时间,计算得出当前的瞬时相位。相比较传统的零交法而言,运算量有所增加,但测量的精度、抗干扰性得到了很大的提高。
2)本发明所涉及的采样可以是等时间间隔采样,也可以是不等时间间隔采样,如果是等时间间隔采样的话,则运算量可以经过优化而进一步减少。
附图说明
图1是由P1i,P2i,P3i通过插值的方式获取P4i,再进一步由T1i和T4i计算过零点的示意图;
图2是计算交流电周期的示意图;
图3是计算某一个采样点瞬时相位的示意图。
具体实施方式
下面结合附图对本发明的技术方案做进一步的详细说明:
本发明的实质是在信号过负峰值后选取一个采样点P1,然后在信号过零点后选择连续的两个点P2和P3,P2和P3的选取条件是P1采样值的负值在区间[P2采样值,P3采样值]内。则可以通过插值的方法在P2和P3所组成的线段上求取一点P4,使得P4的数值等于P1采样值的负值,同时可获取P4点的虚拟发生时间。由P1的采样发生时间和P4的虚拟发生时间计算电信号的零点,由一系列电信号的过零点计算电信号的周期。最后根据前一个时段交流电信号的过零点的时间、预测的当前周波的周期以及当前时间,计算得出当前的瞬时相位。选取的P1不宜过于接近零点,如果接近零点的话容易受到噪声的干扰导致测量结果不够准确。
为获得较为精确的测量结果,建议的采样点P1的选取方案是:不宜过于接近零点。
本发明的特点和优点将通过实例结合附图进行详细说明。本发明的原理通过测量某一个采样点的瞬时相位来进行说明。
本发明的具体实施过程如下:
1.对被测信号进行采样,这里所述的采样为对整周波进行的采样。可以是等时间间隔采样,也可以是不等时间间隔采样。
2.在测电信号过负峰值后选择一个采样点P1i,设该采样点的采样时间为T1i,采样值为S1i;当采样信号越过零点后,在后续的采样值中选取连续的两个点P2i和P3i,设这两个点的采样时间分别是T2i和T3i,两个点的采样值分别为S2i和S3i,这两个点具有以下特点:S1i为负数,S2i和S3i同为正数,且满足|S2i|≤|S1i|≤|S3i|、|S2i|<|S3i|。
3.计算P4i点的虚拟发生时间T4i,由于P4i的数值等于P1i采样值的负值,则:T4i可表示为
Figure BDA00001935331200041
当选取的采样点的采样值符合|S2i|=|S1i|=|S3i|条件时,可以采取两种处理方法:一种方法是重新选取P2i和P3i并重新计算;另一种方法是T4i在区间[T2i,T3i]上任取一点,为方便计算,取该区间的中点,则根据计算获得的T4i,过零点时间的计算如下:
Figure BDA00001935331200043
4.具体操作可以参见图1。图1上空心园点P1i,P2i,P3i是根据要求选择出的三个点,空心三角形点P4i是经过插值后得到的虚拟点,点P4i的数值为点P1i采样值的负值。T1i是P1i的采样时间,T4i是P4i的虚拟发生时间,水平轴上的空心园点Tzi是计算后获得的过零点。
5.得到若干个过零点时间后可进行信号周期的计算,图2中的低频信号一共存在3个负峰值,因此在负峰值出现后一共进行了3次插值并计算,水平轴上时间分别为Tz1、Tz2、Tz3的空心园点是计算得出的3个过零点,由于被测信号的周期可表示为:Tpi=Tzi+1-Tzi,因此图2中信号周期可表示为:Tp2=Tz3-Tz2、Tp1=Tz2-Tz1
6.测量交流电瞬时相位的一个前提是获得当前周波的周期,但是在当前周波没有结束之前,当前周波的周期还无法直接获得,考虑到交流电的周期(频率)一般不会发生急剧的变化,本发明的实例采用了一个公式来预测当前周波的周期Tp。该公式可表示为:Tp=K1×Tpk+K2×Tpk-1+K3×Tpk-2+K4×Tpk-3,附图进行计算时采用了K1=2,K2=-1,K3=K4=0。
7.瞬时相位的计算可参见图3,计算的是图3中采样点S的瞬时相位,设图3中采样点S的采样时间是Tc,经计算得到的上一个交流电的过零点的时间为Tz,则采样点S的瞬时相位可表示为
Figure BDA00001935331200051
其中的Tp是先前预测得到的交流电当前周波的周期。计算得出的结果是以度数表示的相位值。图3中的其他采样点的瞬时相位也可以采用类似的方法获得。一个周波内任意点的瞬时相位的计算方法同上。
优选方案是所述的采样为等时间间隔采样。
综上所述,本发明所涉及的一种交流电的瞬时相位的数字化测量方法,不需要将交流电转换成方波,克服了传统的过零点检测抗干扰能力差的缺点,尤其适合于交流电的瞬时相位测量,系统检测方便,硬件开销小;系统的抗干扰能力较强;实现起来简单合理,数据准确、可靠;有较广泛的应用前景。当然,本发明也可以应用于低频正弦信号的瞬时相位测量。

Claims (7)

1.一种交流电的瞬时相位的数字化测量方法,其特征在于包括以下步骤:
步骤1),对被测电信号进行采样,得到被测电信号的采样值,然后生成采样信号;
步骤2),在被测电信号的采样信号中,当被测电信号的周波过负峰值后未到达零点时选择一个采样点作为第一采样点P1i,设该第一采样点P1i的采样时间为T1i,采样值为S1i
步骤3),当被测电信号的周波越过零点后未到达正峰值时,在该段采样信号中选取连续的两个采样点,即第二采样点P2i和第三采样点P3i,设这两个点的采样时间分别是T2i和T3i,采样值分别为S2i和S3i,其中采样值S2i和S3i需满足|S2i|≤|S1i|≤|S3i|;
步骤4),将第一采样点P1i采样值的相反值所对应的点作为虚拟采样点P4i,计算虚拟采样点P4i的虚拟发生时间 T 4 i = T 2 i + ( T 3 i - T 2 i ) &times; ( S 1 i + S 2 i ) S 2 i - S 3 i ;
步骤5),计算过零点时间:
Figure FDA00001935331100012
步骤6),重复步骤2)至步骤5)得到若干个过零点时间:Tz1,Tz2,Tz3...Tzi,Tzi+1,Tzi+2...Tzn,其中,i=1、2、…n,n为正整数;则被测电信号的周期:Tpi=Tzi+1-Tzi,因此获得的信号的周期依次为Tp1,Tp2,Tp3,...Tpn-1
步骤7),根据步骤6)计算得出的被测电信号的周期,预测当前周波的周期;
步骤8),根据前一个时段交流电信号的过零点的时间、步骤7)预测的当前周波的周期以及当前时间,计算得出当前的瞬时相位。
2.根据权利要求1所述的一种交流电的瞬时相位的数字化测量方法,其特征在于:步骤1)所述的采样是等时间间隔采样或不等时间间隔采样。
3.根据权利要求1所述的一种交流电的瞬时相位的数字化测量方法,其特征在于:步骤1)所述的采样为对信号的整周波进行采样。
4.根据权利要求1所述的一种交流电的瞬时相位的数字化测量方法,其特征在于:所述步骤7)的具体步骤如下:
设由步骤6)计算得到的被测电信号前几个周波的周期依次为Tpk-3,Tpk-2,Tpk-1,Tpk;则预测的当前周波的周期为:Tp=K1×Tpk+K2×Tpk-1+K3×Tpk-2+K4×Tpk-3;其中,K1,K2,K3,K4为任意常数,但须满足条件K1+K2+K3+K4=1,k为大于3的正整数。
5.根据权利要求1所述的一种交流电的瞬时相位的数字化测量方法,其特征在于:所述步骤8)的具体步骤如下:设被测电信号的前一个时段过零点的时间为Tz,当前时间为Tc,由步骤7)预测得到的被测电信号当前周波的周期为Tp,则被测电信号当前的瞬时相位以度数值表示为:
Figure FDA00001935331100021
或者以弧度值表示为:
Figure FDA00001935331100022
6.根据权利要求1所述的一种交流电的瞬时相位的数字化测量方法,其特征在于:在所述步骤3)中,当选取的第一至第三采样点的采样值符合|S2i|=|S1i|=|S3i|条件时,采取以下两种处理方案中的任意一种:
a),重新选取第二采样点P2i和第三采样点P3i,使得第二采样点的采样值S2i和第三采样点的采样值S3i须满足|S2i|<|S3i|,然后继续步骤4)至步骤8);
b),在第二采样点P2i与第三采样点P3i的采样时间区间[T2i,T3i]上任取一点作为虚拟采样点P4i所对应的虚拟采样时间T4i,然后继续步骤5)至步骤8)。
7.根据权利要求6所述的一种交流电的瞬时相位的数字化测量方法,其特征在于:在所述方案b)中,取第二采样点P2i与第三采样点P3i的采样时间区间[T2i,T3i]上的中点,计算虚拟采样点P4i的虚拟发生时间
Figure FDA00001935331100023
CN201210261017.3A 2012-07-26 2012-07-26 一种交流电的瞬时相位的数字化测量方法 Expired - Fee Related CN102778606B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210261017.3A CN102778606B (zh) 2012-07-26 2012-07-26 一种交流电的瞬时相位的数字化测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210261017.3A CN102778606B (zh) 2012-07-26 2012-07-26 一种交流电的瞬时相位的数字化测量方法

Publications (2)

Publication Number Publication Date
CN102778606A true CN102778606A (zh) 2012-11-14
CN102778606B CN102778606B (zh) 2014-11-19

Family

ID=47123577

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210261017.3A Expired - Fee Related CN102778606B (zh) 2012-07-26 2012-07-26 一种交流电的瞬时相位的数字化测量方法

Country Status (1)

Country Link
CN (1) CN102778606B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104330630A (zh) * 2014-10-31 2015-02-04 许继电气股份有限公司 一种高精度模拟量相位差测量方法
CN104614586A (zh) * 2015-01-13 2015-05-13 上海惠安系统控制有限公司 一种基于gps的电气量信号自适应采样方法
CN105954589A (zh) * 2016-04-19 2016-09-21 万高(杭州)科技有限公司 一种多芯片系统的相位测量方法及装置
CN113049872A (zh) * 2021-02-02 2021-06-29 厦门大学 一种用于多次低频电流纹波抑制的电信号零点检测方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1815248A (zh) * 2005-02-02 2006-08-09 艾默生网络能源系统有限公司 交流电频率监测方法
CN101833036A (zh) * 2010-04-15 2010-09-15 南京邮电大学 一种交流电的瞬时相位测量方法
CN102095937A (zh) * 2010-12-17 2011-06-15 南京邮电大学 一种快速测量交流电信号瞬时相位的方法
CN102095935A (zh) * 2010-12-17 2011-06-15 南京邮电大学 一种交流电信号的瞬时相位测量方法
JP2012026836A (ja) * 2010-07-22 2012-02-09 Shindengen Electric Mfg Co Ltd 分散型電源の周波数検出方法及び系統連系保護装置
CN102435844A (zh) * 2011-11-01 2012-05-02 南京磐能电力科技股份有限公司 一种频率无关的正弦信号相量计算方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1815248A (zh) * 2005-02-02 2006-08-09 艾默生网络能源系统有限公司 交流电频率监测方法
CN101833036A (zh) * 2010-04-15 2010-09-15 南京邮电大学 一种交流电的瞬时相位测量方法
JP2012026836A (ja) * 2010-07-22 2012-02-09 Shindengen Electric Mfg Co Ltd 分散型電源の周波数検出方法及び系統連系保護装置
CN102095937A (zh) * 2010-12-17 2011-06-15 南京邮电大学 一种快速测量交流电信号瞬时相位的方法
CN102095935A (zh) * 2010-12-17 2011-06-15 南京邮电大学 一种交流电信号的瞬时相位测量方法
CN102435844A (zh) * 2011-11-01 2012-05-02 南京磐能电力科技股份有限公司 一种频率无关的正弦信号相量计算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
GERSON EDUARDO MOD等: "Zero Crossing Determination by Linear Interpolation of Sampled Sinusoidal Signals", 《20041EEEIPES TRANSMISSION & DISTRIBUTION CONFERENCE & EXPOSITION: LATIN AMERICA》 *
尚秋峰等: "基于虚拟仪器技术的信号瞬时频率提取", 《电测与仪表》 *
李辉等: "基于经验模态分解的瞬时相位分析方法的应用", 《振动、测试与诊断》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104330630A (zh) * 2014-10-31 2015-02-04 许继电气股份有限公司 一种高精度模拟量相位差测量方法
CN104614586A (zh) * 2015-01-13 2015-05-13 上海惠安系统控制有限公司 一种基于gps的电气量信号自适应采样方法
CN104614586B (zh) * 2015-01-13 2019-01-18 上海惠安系统控制有限公司 一种基于gps的电气量信号自适应采样方法
CN105954589A (zh) * 2016-04-19 2016-09-21 万高(杭州)科技有限公司 一种多芯片系统的相位测量方法及装置
CN105954589B (zh) * 2016-04-19 2019-03-22 杭州万高科技股份有限公司 一种多芯片系统的相位测量方法及装置
CN113049872A (zh) * 2021-02-02 2021-06-29 厦门大学 一种用于多次低频电流纹波抑制的电信号零点检测方法
CN113049872B (zh) * 2021-02-02 2021-11-16 厦门大学 一种用于多次低频电流纹波抑制的电信号零点检测方法

Also Published As

Publication number Publication date
CN102778606B (zh) 2014-11-19

Similar Documents

Publication Publication Date Title
CN101833036B (zh) 一种交流电的瞬时相位测量方法
CN101806832B (zh) 一种低频率信号的频率测量方法
CN101813725B (zh) 一种低频率信号的相位差测量方法
CN100432681C (zh) 交流电频率监测方法
CN102435844B (zh) 一种频率无关的正弦信号相量计算方法
CN102809687B (zh) 一种交流电频率的数字化测量方法
CN102095934B (zh) 一种交流电信号的相位差测量方法
CN102818930B (zh) 一种高精度快速计算电力谐波参数的方法
CN102879639A (zh) 一种电力系统中频率的实时测量方法
CN102778606B (zh) 一种交流电的瞬时相位的数字化测量方法
CN102095929A (zh) 一种快速测量交流电信号频率的方法
CN102928666B (zh) 一种交流电的相位差的数字化测量方法
CN103018555A (zh) 一种高精度的电力参数软件同步采样方法
CN102095935B (zh) 一种交流电信号的瞬时相位测量方法
CN105486921A (zh) 凯撒三阶互卷积窗三谱线插值的谐波与间谐波检测方法
CN103575979B (zh) 一种数字化测量交流电频率的方法
CN102508022A (zh) 采用最优乘子牛顿算法检测电网频率的方法
CN102095936B (zh) 一种快速测量交流电信号相位差的方法
CN103543331A (zh) 一种计算电信号谐波和间谐波的方法
CN103575981A (zh) 一种交流电频率的精确测量方法
RU2540443C1 (ru) Способ определения места обрыва на воздушной линии электропередачи
CN104808060A (zh) 一种电信号相位差的数字化测量方法
CN103575991B (zh) 一种交流电的瞬时相位的精确测量方法
RU2482517C1 (ru) Трассопоисковый приемник
CN104808055A (zh) 一种电信号频率的数字化测量方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20121114

Assignee: Jiangsu Nanyou IOT Technology Park Ltd.

Assignor: Nanjing Post & Telecommunication Univ.

Contract record no.: 2016320000215

Denomination of invention: Digital measurement method for instantaneous phase of alternating current

Granted publication date: 20141119

License type: Common License

Record date: 20161118

LICC Enforcement, change and cancellation of record of contracts on the licence for exploitation of a patent or utility model
EC01 Cancellation of recordation of patent licensing contract

Assignee: Jiangsu Nanyou IOT Technology Park Ltd.

Assignor: Nanjing Post & Telecommunication Univ.

Contract record no.: 2016320000215

Date of cancellation: 20180116

EC01 Cancellation of recordation of patent licensing contract
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20141119

Termination date: 20170726

CF01 Termination of patent right due to non-payment of annual fee