CN101871965B - 电力正弦信号过零时间、频率、相位差的检测方法 - Google Patents
电力正弦信号过零时间、频率、相位差的检测方法 Download PDFInfo
- Publication number
- CN101871965B CN101871965B CN2010101980808A CN201010198080A CN101871965B CN 101871965 B CN101871965 B CN 101871965B CN 2010101980808 A CN2010101980808 A CN 2010101980808A CN 201010198080 A CN201010198080 A CN 201010198080A CN 101871965 B CN101871965 B CN 101871965B
- Authority
- CN
- China
- Prior art keywords
- signal
- zero
- zero crossing
- overbar
- frequency
- 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.)
- Active
Links
Images
Landscapes
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
本发明公开了一种电力正弦信号过零时间、频率、相位差的检测方法,所述过零时间的检测方法,是利用正弦信号在过零附近都可以近似等效为线性信号的基本原理,在粗略检测到信号过零位置后,通过一元线性回归理论统计信号过零附近电压幅值,从而获得正弦信号过零时间。对正弦信号过零点频率进行检测的方法,是用所得信号各过零点相位差,通过累积平均后得到信号每个过零点的周期值,然后求倒计算可得信号频率值。对相位差进行检测的方法,是用各不同信号的过零点相位值,算出不同信号之间的相位差。本发明能对实际信号的过零点的相对位置辨别不敏感,同时能在数字检测系统中实现较高检测准确度,而且对信号源的噪声与干扰也有很强的抑制作用,最终保证了该方法在高精度检测场所中的应用。
Description
技术领域
本发明涉及一种电力正弦信号的检测方法,尤其是正弦信号过零时间、频率、相位差的检测方法。
背景技术
频率是电力系统的主要性能参数,是各电力系统设备、仪器仪表的必检电气参数之一。目前,通过国内外学者的研究与发展,已经发明了许多基于各种原理与方法的检测手段,有基于简单软硬件实现的脉冲计数法、过零检测法,以及基于复杂算法的傅里叶变换法、卡尔曼滤波法等。
在众多检测法中,过零判别检测法是一种简单实用、应用广泛的电力正弦信号相位、频率检测方法。其基本原理是通过判别正弦电压过零点来确定信号相位,以及利用信号两过零点的时间间隔来计算频率。但在实际应用中由于电力系统谐波、噪声与干扰等原因,导致检测到的实际信号在过零点附近来回抖动,不但增加了软件过零判别、信号去抖的工作量,而且使得相位、频率检测误差较大,这限制了该方法在精密设备中的应用。另外,数字化检测设备处理的都是模拟信号经ADC采样后的数字信号,由于信号采样率的限制,数字信号相位有最小分辨率问题。因此,在数字化设备中,常规的过零检测法难以获得较高的信号频率、相位检测精度,这进一步限制了其在高精度检测场所的应用。
发明内容
针对现有技术在实际应用中过零点辨别困难以及数字信号处理中相位分辨率不足的问题,本发明提供了一种对信号过零检测不敏感,对噪声、干扰抑制能力强,并能有效提高数字正弦信号频率、相位检测准确度的电力正弦信号过零时间、频率、相位差的检测方法。
本发明提供的这种电力正弦信号过零时间的检测方法,利用正弦信号在过零附近都可以近似等效为线性信号的基本原理,在粗略检测到信号过零位置后,通过一元线性回归理论统计信号过零附近电压幅值,获得正弦信号过零时间,所述正弦信号过零时间适用下列公式:
其中,右式为信号在内存中的计算公式,式中:
所取分析信号的时间均值
u(k):所取分析信号采样幅值
k:信号过零点在内存中的地址偏移量
n:所取分析信号的点数
在信号过零左右基本对称地取n个数字采样值作为独立观察样本,为:(ti,ui),i=1,2,…,n。
本发明提供的对正弦信号过零点频率进行检测的方法,是用所得信号各过零点时间差,通过累积平均后得到信号每个过零点的周期值,然后求倒计算可得信号频率值。
检测频率适用如下公式:
f:被测信号频率
f0:标准信号频率值,中国电力系统中f0=50Hz;
T0:标准信号周期值T0=1/f0;
N:标准信号周期T0内的采样点数;
本发明提供的对正弦信号过零点相位差进行检测的方法,是用各不同信号的过零点时间差,换算出不同信号之间的相位差。
检测相位差适用如下公式:
Δω:信号间相位差
两信号多过零点时间差平均值
f:被测信号频率
f0:标准信号频率值,中国电力系统中f0=50Hz;
T0:标准信号周期值T0=1/f0;
N:标准信号周期T0内的采样点数。
本发明电力正弦信号过零时间、频率及相位差检测方法的工作原理是:首先利用正弦信号在过零附近都可以近似等效为线性信号的基本原理,在粗略检测到信号过零位置后,将过零点附近左右的信号值进行一元线性回归分析,从而计算得到该信号在横轴坐标上的截距。该截距即为该信号过零时间,利用该信号上不同过零时间之间的时间差即可计算得出该信号的周期,进而可精确计算出信号频率;利用不同信号上的过零点时间差,即可计算得出不同信号之间的相位关系。因为一元线性回归分析是满足最小二乘法原理的信号统计分析方法,它所计算出来的结果满足误差平方和最小关系,因此它不仅对信号噪声和干扰有很强的抑制作用,而且计算结果也具有很高的准确性。由此使本发明能对实际信号的过零点的相对位置辨别不敏感,同时能在数字检测系统中实现较高的相位分辨率,而且对信号源的噪声与干扰也有很强的抑制作用,最终保证了数字化设备在高精度检测场所的应用。
附图说明
图1是正弦信号过零点线性化示意图。
图2是在不同采样点数下的算法实现效果图。
图3是信号频率、相位检测算法流程图。
图4是信号过零检测程序流程图。
图5是过零检测算法示意图。
图6是算法过零查找与频率计算示意图。
图7是过零信号左右指标关系图。
图8是一元线性回归算法程序流程图。
具体实施方式
以下结合附图对本发明做详细说明:
参见图1设输入正弦信号为众所周知所有正弦信号在过零点附近都可以近似等效为线性信号。这是因为正弦函数的幂级数为在正弦信号过零处附近x取值一般都很小,因此忽略高次项式可近似为sinx≈x,如附图1所示。若近似式与真值相对误差要求≤0.1%,则x取值范围为x∈[-0.0774,0.0774]。
在实际问题中,研究相互联系的变量之间的相互关系通常有两类:一类是确定性关系,其特点是由一个或一组变量的值可唯一确定另一个变量的值;另一类是关系是非确定关系,其特点是变量之间相互联系,但又不能由一组变量的值唯一确定另一个变量的值。回归分析就是为了寻找不完全确定变量间的数量关系,通过统计方法进行推断的一种方法。当自变量只有一个的回归分析叫一元回归分析,同时数学关系式为线性时称为一元线性回归。
从电力系统的宏观来看,系统电压波形、幅值与时间变量t为正弦函数关系,大体上属于确定性关系,并且在过零附近可以局部线性化;但是由于电气噪声与系统干扰的影响,该正弦函数的参数:电压Um、频率f与相位都是服从一定分布函数的随机变量,导致系统在过零点相位这类局部特征呈现一定随机性。因此,我们可以把u(t)的过零点附近的观察结果看成由两部分叠加而成的:一部分由时间t的线性函数引起,记为a·t+b;另一部分是由随机因素引起的,记为ε,即
u(t)=a·t+b+ε (1)
式中参数a与b主要是由信号的电压、频率、相位三者随机变量的数学期望值EUm、Ef、决定;对于随机变量ε,它是引起信号过零点抖动、甚至重复过零的主要因素,它通常服从正态分布N(0,σ2)。于是可知,信号过零局部u(t)是随机变量,且
u(t)~N(a·t+b,σ2) (2)
而Eu(t)=a·t+b是时间t的线性函数。一般地,在数字检测系统中信号通常由模拟信号u(t)经过模数转换为离散数字信号。将输入按每周期采样点数为N进行数字化处理,则可得采样数字信号为注意式中电压Um、频率f与相位都是随机变量,f0是标准信号频率值,在中国电力系统中,该值为50Hz。在信号过零左右基本对称地取n个数字采样值作为独立观察样本,为:
(ti,ui),i=1,2,…,n (3)
那么,利用u、t之间近似线性关系,则可假定
式中各εi相互独立,i=1,2,…,n,这就是数字正弦信号在过零点附近的一元线性回归的数学模型。
现在求解式(8)关于参数的a、b的最小二乘估计,记令:
整理得:
联立求解可得:
最后,将式(11)代入式(6)即可计算当前信号过零时间的最佳估计值:
利用相同信号前后两过零时间差即为信号周期,求倒得到该信号的频率;利用不同信号的过零时间差,即可换算得到两信号间的相位差。
关于所需采样点数的问题,由于本发明是基于正弦电压过零处信号可以局部线性化的这一特性为基本前提的,因此所需采样点数应保证信号线性化所要求精度。例如,前面所述线性化误差要求≤0.1%,正弦信号过零点附近最大弧度应≤0.0774。设信号采样率为每周期256点,信号基频为50Hz,那么每个采样点间隔弧度最小分辨率为因此,在信号过零点左右可各取个点共6点进行回归分析。
事实上,由于工程上更关心的是信号过零点时间的准确位置(即信号与时间轴的截距),而并非信号拟合曲线精度。另外,由于正弦信号在零点左右为奇对称信号,因此实际算法选择点数条件可以放宽,而并非需要受线性化误差的限制。从理论上分析,设输入信号为N为信号周期内采样点数。以零点为中心左右对称的共取n个点时,由于正弦信号的对称性与时间轴选择的对称性计算可得 代入式(12)可得t=0。图2给出了算法在周波采样点数256点带噪声信号下不同点数的计算结果效果图。在绝大对多的实际检测系统中,6点分析已经可以得到非常理想的检测精度与效果。但在一些干扰特别大的特殊应用场合,更多的采样点数可以使算法对信号噪声有更好的抑制作用,这需要根据系统设计与实际应用环境而定。
图3给出了基于一元线性回归理论的数字正弦信号频率测量(a)、信号间相位差检测(b)的算法实现流程图。这两个算法的实现流程基本一致,主要都是由信号过零检测、信号回归分析与目标参数计算三步实现。下面对算法流程的实现过程进行详细介绍:
1.过零检测
图4给了信号过零检测流程图。因为正弦信号在一个周期里有正向穿越与负向穿越两次过零,为了防止程序重复判断过零点,流程中增加了一个过零检测使能位EnableFlg,EnableFlg=1时使能检测,EnableFlg=0时禁止检测。它只有在当前信号值小于检测使能负电平-Ven后才会开启信号的过零检测。其中,Ven可以设为系统最小输入交流电压峰值的80%。前面算法分析可知,回归算法是对信号过零点左右信号进行统计分析得到计算结果,对信号实际过零点位置不敏感。因此程序在过零检测使能情况下,只要简单辨断当前信号开始大于零即可,这样极大降低了信号过零的辨别条件。图4过零检测流程的过程示意图,如图5所示。
2.过零信号线性回归分析
在发现信号过零后,即要开展过零信号的回归分析。由于回归算法需要对过零点左右信号进行统计计算,但是检测设备在实时采样过程中当前发现过零点的右侧信号还未采得,所以算法只能对上一过零位置附近的信号进行分析,并结合前一过零信号计算结果进行参数计算,整个过程可参见图6所示。
另外,由于实际数字检测系统中,采样信号都是保存在系统内部的一片连续的缓冲区内。因此,在回归算法实际实现过程中,时间参量应为信号存储在系统内存中的信号偏移量。设当前信号过零点地址偏移量为k,在当前过零点左右对称的选择n个点进行回归分析,如图7所示。通常为了算法实现方便,一般取n为偶数,利用图7指标分布规律可对回归分析式(12)进行化简得:
将上述各式计算结果代回式(12)得到化简式:
图8为式(16)计算的程序设计流程图,算法在选定采样点数后可以将部分算式进行预先处理。例如取8点进行分析,即n=8,则代入式(16)可得:
3.目标参数计算
为了进一步提高算法对信号噪声与干扰的抑制能力,对于频率、相位差这类目标参数的计算采用了多周期累积求平均值的方法。另外,在大多数电力检测系统中,秒级左右的数据更新速度已经能够完全满足用户需求。在本发明的应用中,对频率采用了32周波(数据更新间隔时间约为32×0.02S=0.64S)累积数据求一次周期平均值:
式中Ti为前述算法每次发现过零后计算所得的信号周期。因为在数字系统中,是以信号采样时间间隔为单位时间的数字量,所以最后需要将该值转换为时间量,最后求倒得到信号频率(设采样率为T0为信号额定周期,为信号额定频率,中国电力系统中f0=50Hz):
对相位差计算采用了64周波(数据更新间隔时间约为64×0.02S=1.28S)累积值求一次平均值:
式中f为当前信号频率。
在MATLAB环境对本发明方法进行仿真,仿真中采样率取12.8kHz(基波50Hz每周期采样256点),在零点左右各取4点共8点进行分析。为验证本发明方法的有效性,以下仿真的被检信号取非整数频率,表1给出了仿真结果以及相对误差,可见在无噪声下本发明方法计算误差有效位可达6位数以上。
表1MATLAB仿真实验结果
实际现场环境中信号中夹杂了大量随机噪声,为能真实模拟现场环境,仿真中对信号加入方差为2的正态分布随机噪声,表2给出了仿真结果以及相对误差。可见在噪声环境下,算法也有非常高的准确度。
表2仿真实验结果(加入信号噪声)
将本发明方法设计成程序应用于0.1S高准确度等级电能表上,该表设计采用TI高精度24bits多通道ADC芯片ADS1278,配合ADI的高速DSP芯片BF533。实现本发明方法时仍采用8点分析,由于ADC分辨率非常高,因此用本发明方法在过零检测中使用的过零检测使能电平为5V,在含有大量谐波的应用场合,该使能电平可以适当提高。应用三晖0.01%高准确度台体进行实验验证,得到实验数据如表3所示。
表3实验结果
实验结果表明用本发明方法检测结果准确度高、各相之间数据一致性好、方法稳定、可靠。并且本发明方法能使编程中流程清晰、过程简洁、运算量很小。
Claims (5)
2.一种用权利要求1所述方法对正弦信号频率进行检测的方法,其特征是用所得信号各过零点时间差,通过累积平均后得到信号每个过零点的周期平均值,然后求倒计算可得信号频率值。
4.一种用权利要求1所述方法对两正弦信号相位差进行检测的方法,其特征是用各不同信号的过零点时间差,换算出不同信号之间的相位差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101980808A CN101871965B (zh) | 2010-06-11 | 2010-06-11 | 电力正弦信号过零时间、频率、相位差的检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101980808A CN101871965B (zh) | 2010-06-11 | 2010-06-11 | 电力正弦信号过零时间、频率、相位差的检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101871965A CN101871965A (zh) | 2010-10-27 |
CN101871965B true CN101871965B (zh) | 2012-01-04 |
Family
ID=42996951
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010101980808A Active CN101871965B (zh) | 2010-06-11 | 2010-06-11 | 电力正弦信号过零时间、频率、相位差的检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101871965B (zh) |
Families Citing this family (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102072988B (zh) * | 2010-11-19 | 2013-02-06 | 山东大学 | 光伏逆变器接入点电压多过零点相位检测方法 |
CN102478599A (zh) * | 2010-11-22 | 2012-05-30 | 浙江昱能光伏科技集成有限公司 | 电网交流电压及其过零点的检测方法及检测装置 |
CN102095929B (zh) * | 2010-12-17 | 2013-01-02 | 南京邮电大学 | 一种快速测量交流电信号频率的方法 |
CN102095934B (zh) * | 2010-12-17 | 2012-11-28 | 南京邮电大学 | 一种交流电信号的相位差测量方法 |
CN102095937B (zh) * | 2010-12-17 | 2013-01-02 | 南京邮电大学 | 一种快速测量交流电信号瞬时相位的方法 |
CN102095936B (zh) * | 2010-12-17 | 2013-01-02 | 南京邮电大学 | 一种快速测量交流电信号相位差的方法 |
CN103575979B (zh) * | 2012-07-26 | 2016-03-02 | 南京邮电大学 | 一种数字化测量交流电频率的方法 |
CN102879639A (zh) * | 2012-09-13 | 2013-01-16 | 华中科技大学 | 一种电力系统中频率的实时测量方法 |
CN103888312B (zh) * | 2014-03-04 | 2017-06-09 | 京信通信系统(广州)有限公司 | 一种预失真系统的报警方法及装置 |
CN104330622B (zh) * | 2014-10-29 | 2017-06-09 | 广东电网有限责任公司电力科学研究院 | 电力系统中正弦波信号的频率测量方法及系统 |
CN105445547B (zh) * | 2015-07-09 | 2018-05-25 | 深圳市科润宝实业有限公司 | 正弦信号的频率检测方法和系统 |
CN105044455B (zh) * | 2015-07-09 | 2017-11-24 | 宁波大学 | 交流信号的频率和相位的测量方法 |
CN105259409B (zh) * | 2015-10-13 | 2017-11-28 | 安徽节源环保科技有限公司 | 一种基于过零时域特性的电力信号频率计算方法 |
CN105652082B (zh) * | 2015-12-30 | 2019-06-11 | 深圳市科陆电子科技股份有限公司 | 高精度电信号频率测量方法及其装置 |
CN106950420A (zh) * | 2017-05-11 | 2017-07-14 | 国家电网公司 | 工频正弦交流电信号的过零时刻检测系统及其方法 |
CN107834678A (zh) * | 2017-08-17 | 2018-03-23 | 山东电工豪迈节能科技有限公司 | 一元线性回归测频算法在直流充电机中的应用 |
EP3508865B1 (en) * | 2018-01-08 | 2022-07-20 | Delta Electronics (Thailand) Public Co., Ltd. | Method for estimating a signal property |
CN108318737A (zh) * | 2018-01-22 | 2018-07-24 | 西安工业大学 | 一种高精度数字鉴相方法 |
CN108957133A (zh) * | 2018-08-13 | 2018-12-07 | 国网浙江省电力有限公司绍兴供电公司 | 变电站电容器采样频率动态评估和选取方法 |
CN109342814B (zh) * | 2018-09-05 | 2021-04-09 | 苏州光格科技股份有限公司 | 用于确定电缆局放工频相位频率的方法、装置及存储介质 |
CN110069984A (zh) * | 2019-03-11 | 2019-07-30 | 珠海兴诺能源技术有限公司 | 一元线性回归测频方法、电子设备及存储介质 |
CN110045399B (zh) * | 2019-04-30 | 2020-11-27 | 中国人民解放军国防科技大学 | 一种基于选择性非等量采样的过零点偏差抑制方法 |
CN112134593B (zh) | 2019-06-24 | 2021-10-22 | 天地融科技股份有限公司 | 过零点信号输出和电力线数据发送方法及设备 |
CN110596456A (zh) * | 2019-08-09 | 2019-12-20 | 广州汉光电气股份有限公司 | 一种电网频率测试方法 |
CN111458558B (zh) * | 2020-04-03 | 2022-03-29 | 宁波奥克斯电气股份有限公司 | 一种参数获取方法、装置、电子设备及存储介质 |
CN111900007B (zh) * | 2020-07-29 | 2022-08-23 | 浙江方圆电气设备检测有限公司 | 一种选相合闸角度精准控制的方法 |
CN111896807B (zh) * | 2020-08-05 | 2023-03-14 | 威胜集团有限公司 | 基波频率测量方法、测量终端及存储介质 |
CN114325094B (zh) * | 2021-12-29 | 2022-11-11 | 中国科学院上海微系统与信息技术研究所 | 一种相位信息测量装置及方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5566761A (en) * | 1978-06-19 | 1980-05-20 | Matsushita Electric Ind Co Ltd | Detecting method and its device for source voltage zero point |
US4959558A (en) * | 1988-03-31 | 1990-09-25 | U.S. Philips Corporation | Circuit arrangement for detecting cross-over by an alternating voltage of a fixed reference voltage level |
WO2006043949A1 (en) * | 2004-10-20 | 2006-04-27 | Semiconductor Components Industries, Llc. | High accuracy zero crossing detector and method therefor |
-
2010
- 2010-06-11 CN CN2010101980808A patent/CN101871965B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN101871965A (zh) | 2010-10-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101871965B (zh) | 电力正弦信号过零时间、频率、相位差的检测方法 | |
CN101806832B (zh) | 一种低频率信号的频率测量方法 | |
CN100480706C (zh) | 利用固定相移来测量同频信号相位差的方法及电路 | |
CN101813725B (zh) | 一种低频率信号的相位差测量方法 | |
CN102288821B (zh) | 三相电路相位差的测量方法、测量装置 | |
CN101833036B (zh) | 一种交流电的瞬时相位测量方法 | |
CN102435844A (zh) | 一种频率无关的正弦信号相量计算方法 | |
CN102508031B (zh) | 一种基于傅里叶级数的局部放电脉冲相角测量方法 | |
CN106324340B (zh) | 一种同步相量和频率测量动态性能的方法 | |
CN102135568B (zh) | 交流电信号的周期检测方法及相位差检测方法 | |
CN102809687B (zh) | 一种交流电频率的数字化测量方法 | |
CN106018958A (zh) | 孤网变频系统中二次侧电流电压频率跟踪方法 | |
CN101907656B (zh) | 利用固定相移来测量同频信号相位差的方法 | |
CN102095929A (zh) | 一种快速测量交流电信号频率的方法 | |
CN103575979A (zh) | 一种数字化测量交流电频率的方法 | |
CN105548711B (zh) | 一种多频信息滤波递推解调方法 | |
CN101702617B (zh) | 高精度±180°数字鉴相方法及其实施装置 | |
Toral et al. | Reactive power and energy measurement in the frequency domain using random pulse arithmetic | |
Muttakin et al. | Design and simulation of quadrature phase detection in electrical capacitance volume tomography | |
Yang et al. | A Hilbert Transform based method for dynamic phase difference measurement | |
CN111551785B (zh) | 基于无迹卡尔曼滤波的频率与谐波检测方法 | |
CN104914305B (zh) | 一种基于最小二乘法的高精度频率估计方法 | |
CN106353596A (zh) | 一种低频弦波信号相位差测量方法 | |
Zaitseva | Some problems of discretization and determination of basic parameters of harmonic signals | |
Tao | A Moving single station passive ranging method by interchangeable relationship between frequency shift and path difference |
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: 20101027 Assignee: Shenzhen Wismei Technology Co.,Ltd. Assignor: Wasion Group Co., Ltd. Contract record no.: 2013440020317 Denomination of invention: Method for detecting zero crossing time, frequency and phase difference of power sinusoidal signals Granted publication date: 20120104 License type: Exclusive License Record date: 20131112 |
|
LICC | Enforcement, change and cancellation of record of contracts on the licence for exploitation of a patent or utility model |