CN102608415B - 基于加权双拟合的软件频率跟踪算法 - Google Patents
基于加权双拟合的软件频率跟踪算法 Download PDFInfo
- Publication number
- CN102608415B CN102608415B CN201210029921.1A CN201210029921A CN102608415B CN 102608415 B CN102608415 B CN 102608415B CN 201210029921 A CN201210029921 A CN 201210029921A CN 102608415 B CN102608415 B CN 102608415B
- Authority
- CN
- China
- Prior art keywords
- coordinate
- parabola
- value
- maximal value
- fit
- 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
Abstract
本发明提供一种基于加权双拟合的软件频率跟踪算法，包括以下步骤：提取正弦波最大值采样数据，并在同一波峰内提取所述最大值前后各两点采样数据，共获得5点采样数据；将最大值及其前两点数据代入拟合抛物线方程，求出第1条拟合抛物线最大值所在的x坐标值；将最大值及其后两点数据代入拟合抛物线方程，求出第2条拟合抛物线最大值所在的x坐标值；根据两条拟合抛物线所求出的x坐标值，按照加权因子计算出采样值波形最大值所在的x坐标值；计算出采样值波形两个相邻最大值的x坐标值差，获得采样值波形频率。本算法简单有效，响应速度快，可靠性高，对采样精度的要求较低，不受频率变化的影响，有较好的适用性。
Description
技术领域
本发明涉及电力技术领域，尤其涉及一种电网频率的软件测量方法。
背景技术
电网频率是指示电力系统工作状态的重要依据，也是电力系统运行的主要控制参数，频率质量直接影响着电力系统的安全、稳定和优质运行，微机保护许多算法都是建立在与采样频率有关的基础之上，因此需要实时监测系统的频率变化，根据系统频率的波动来获得测量和计算的准确性。
频率测量的方法主要有硬件测量和软件测量两种：硬件测量是将交流信号整形变成方波后，测量方波的跳变宽度，但在交流过零时，运放和光电隔离器存在线性区，容易产生方波的抖动，引起频率测量误差，并且需要额外增加硬件测频电路，占用单片机外部中断。软件测量是根据相似三角形的频率计算方法，利用交流信号过零时近似直线的特点计算两个过零点的宽度，这种方法对于A/D采样精度要求较高，并且采样值在过零时数值较小，离散性较大，频率测量结果不稳定。
发明内容
本发明提出一种新型的、有效的电力系统软件频率跟踪算法。如图1所示，通过采样取值拟合两个抛物线方程，并进行加权运算，获得抛物线极值之后，通过与相邻波峰极值求差从而测得电网频率。本发明提供的算法对硬件要求低，且测量结果稳定。
为了达到上述目的，本发明提供如下技术方案：
基于加权双拟合的软件频率跟踪算法，包括以下步骤：
（1）提取正弦波最大值采样数据，并在同一波峰内提取所述最大值前后各两点采样数据，共获得5点采样数据；
（2）将最大值及其前两点数据代入拟合抛物线方程，求出第1条拟合抛物线最大值所在的 坐标值；
（3）将最大值及其后两点数据代入拟合抛物线方程，求出第2条拟合抛物线最大值所在的坐标值；
（4）根据两条拟合抛物线所求出的坐标值，按照加权因子计算出采样值波形最大值所在的坐标值；
（5）计算出采样值波形两个相邻最大值的坐标值差，获得采样波形频率。
作为优选方案，所述步骤（2）和步骤（3）中拟合抛物线方程为。
作为优选，所述 5点采样数据，分别用点进行标识，其中为最大值点，为同一波峰内最大值点之前的两个点，即，为同一波峰内最大值点之后的两个点，即。
作为优选，所述步骤（2）为：将代入方程，计算出两个系数，然后求出第1条拟合抛物线最大值所在的坐标值。
作为优选，所述步骤（3）为：将代入方程，计算出两个系数，然后求出第2条拟合抛物线最大值所在的坐标值。
作为优选，所述步骤（4）为：根据公式计算出采样值波形最大值，其中，为加权因子。
作为优选，步骤（5）为：根据公式计算出采样值波形两个相邻最大值的坐标值差，并根据计算获得采样值波形频率。
该算法简单有效，响应速度快，可靠性高，对采样精度的要求较低，不受频率变化的影响，有较好的适用性。此外无需硬件开支，克服了硬件测量方法的不足，消除了硬件测频由于器件的因素而引入的误差，提高频率测量的精度。
附图说明
图1为本发明公开的算法中被采样跟踪的正弦波形及拟合抛物线示意图。
其中x为采样点的间距，相邻两个采样的间隔为1，y为采样值数据。
具体实施方法
以下将结合具体实施例对本发明提供的技术方案进行详细说明：
实施例一：
在本实施例中，被测试的交流信号频率为52Hz，峰值为11585，软件频率跟踪过程如下：
（1）采用数字化变电站通用的每周波80点采样速率，对该交流信号的正弦波形进行采样。取该正弦波最大值及其同一波峰内最大值的前后各两点，共获得5点采样值数据，分别用点标识，本例中采样点数据为，其中为最大值点，为同一波峰内之前的两个点，为同一波峰内之后的两个点。
（2）将最大值及其前两点数据代入拟合抛物线方程，计算出，然后求出第1条拟合抛物线最大值所在的坐标值；
（3）将最大值及其后两点数据代入拟合抛物线方程，计算出，然后求出第2条拟合抛物线最大值所在的坐标值；
（4）根据两条拟合抛物线求出的最大值所在坐标值，按照加权因子计算出优化后的采样值波形最大值所在的坐标值；
（5）同样提取下一个采样值正弦波最大值及其前后各两点，共5点采样值数据为，其中为最大值点；同样按照上述步骤（2）～（4）方法，求出下一个最大值，。由此计算出采样值波形两个相邻最大值的坐标值差，获得采样值波形频率，频率误差。
由上述实施例可知，本算法精度高，误差小，且由于采样值在峰值附近，数值较大，离散性较小，易于测量。
Claims (5)
1.一种基于加权双拟合的软件频率跟踪算法，其特征在于，包括以下步骤：
(1)提取正弦波最大值采样数据，并在同一波峰内提取所述最大值前后各两点采样数据，共获得5点采样数据；
(2)将最大值及其前两点数据代入拟合抛物线方程，求出第1条拟合抛物线最大值所在的x坐标值；
(3)将最大值及其后两点数据代入拟合抛物线方程，求出第2条拟合抛物线最大值所在的x坐标值；
(4)根据两条拟合抛物线所求出的x坐标值，按照加权因子计算出采样值波形最大值所在的x坐标值；
(5)计算出采样值波形两个相邻最大值的x坐标值差，获得采样波形频率。
2.根据权利要求1所述的基于加权双拟合的软件频率跟踪算法，其特征在于：所述步骤(2)和步骤(3)中拟合抛物线方程为y＝a1x2+b1x+c1。
3.根据权利要求1或2所述的基于加权双拟合的软件频率跟踪算法，其特征在于：步骤(1)中所述5点采样数据，分别用点p1(x1,y1),p2(x2,y2),p3(x3,y3),p4(x4,y4),p5(x5,y5)进行标识，其中p3(x3,y3)为最大值点，p1(x1,y1),p2(x2,y2)为同一波峰内最大值点之前的两个取值点，p4(x4,y4),p5(x5,y5)为同一波峰内最大值点之后的两个取值点。
4.根据权利要求3所述的基于加权双拟合的软件频率跟踪算法，其特征在于所述步骤(2)为：将p1(x1,y1),p2(x2,y2),p3(x3,y3)代入方程y＝a1x2+b1x+c1，计算出a1,b1两个系数，然后求出第1条拟合抛物线最大值所在的x坐标值xm1＝-b1/(2a1)。
5.根据权利要求4所述的基于加权双拟合的软件频率跟踪算法，其特征在于所述步骤(3)为：将p3(x3,y3),p4(x4,y4),p5(x5,y5)代入方程 y＝a2x2+b2x+c2，计算出a2,b2两个系数，然后求出第2条拟合抛物线最大值所在的x坐标值xm2＝-b2/(2a2)，
所述步骤(4)为：根据公式xmax＝k1*xm1+k2*xm2计算出采样值波形最大值xmax；
所述步骤(5)为：根据公式T＝|xmax1-xmax2|计算出采样值波形两个相邻最大值的x坐标值差T，并根据T计算获得采样值波形频率f。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210029921.1A CN102608415B (zh) | 2012-02-10 | 2012-02-10 | 基于加权双拟合的软件频率跟踪算法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210029921.1A CN102608415B (zh) | 2012-02-10 | 2012-02-10 | 基于加权双拟合的软件频率跟踪算法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102608415A CN102608415A (zh) | 2012-07-25 |
| CN102608415B true CN102608415B (zh) | 2015-06-17 |
Family
ID=46525942
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210029921.1A Active CN102608415B (zh) | 2012-02-10 | 2012-02-10 | 基于加权双拟合的软件频率跟踪算法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102608415B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102879639A (zh) * | 2012-09-13 | 2013-01-16 | 华中科技大学 | 一种电力系统中频率的实时测量方法 |
| CN103076487B (zh) * | 2013-01-10 | 2015-05-27 | 安徽省临泉县嘉柏列科技有限公司 | 大功率电池测试设备采样负载突变电压峰值算法 |
| CN104569581B (zh) * | 2015-01-30 | 2017-05-03 | 湖南大学 | 一种电网频率测量的多水平集单周期估计方法 |
| CN105403767A (zh) * | 2015-10-21 | 2016-03-16 | 广东美的制冷设备有限公司 | 输入空调器的交流电源的电压频率检测方法、系统和空调器 |
| CN105652082B (zh) * | 2015-12-30 | 2019-06-11 | 深圳市科陆电子科技股份有限公司 | 高精度电信号频率测量方法及其装置 |
| CN108226633B (zh) * | 2018-01-02 | 2020-12-11 | 京东方科技集团股份有限公司 | 频率检测方法和频率检测装置 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4761703A (en) * | 1987-08-31 | 1988-08-02 | Electric Power Research Institute, Inc. | Rotor fault detector for induction motors |
| CN101806832A (zh) * | 2010-04-15 | 2010-08-18 | 南京邮电大学 | 一种低频率信号的频率测量方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7684759B2 (en) * | 2006-05-19 | 2010-03-23 | Siemens Industry, Inc. | Method for determining the damped natural frequencies of a dynamic system |
-
2012
- 2012-02-10 CN CN201210029921.1A patent/CN102608415B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4761703A (en) * | 1987-08-31 | 1988-08-02 | Electric Power Research Institute, Inc. | Rotor fault detector for induction motors |
| CN101806832A (zh) * | 2010-04-15 | 2010-08-18 | 南京邮电大学 | 一种低频率信号的频率测量方法 |
Non-Patent Citations (4)
| Title |
|---|
| 张凤蛟.快速曲线拟合的方法.《延边大学学报(自然科学版)》.2006,第32卷(第3期),第208-211页. * |
| 杨晓琴.等高线生成算法的研究.《中国优秀博硕士学位论文全文数据库(硕士)基础科学辑》.2004,(第04期),正文第34-38页. * |
| 顾及光滑拟合精度的地形曲线采样点布设规则研究;李伟涛等;《资源开发与市场》;20111231;第27卷(第10期);第884-886页 * |
| 飞机AMU自动测试设备数据采集分析系统研究;黄静静;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20070228(第02期);正文第41页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102608415A (zh) | 2012-07-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102608415B (zh) | 基于加权双拟合的软件频率跟踪算法 | |
| CN106841778B (zh) | 基于pmu实现的次同步和超同步谐波参数的处理方法 | |
| CN103989462B (zh) | 一种脉搏波形第一特征点和第二特征点的提取方法 | |
| CN106405230B (zh) | 频率测量方法和装置 | |
| CN102818930B (zh) | 一种高精度快速计算电力谐波参数的方法 | |
| CN106018958B (zh) | 孤网变频系统中二次侧电流电压频率跟踪方法 | |
| CN102879639A (zh) | 一种电力系统中频率的实时测量方法 | |
| CN103412200A (zh) | 一种三相交流电相序检测方法及装置 | |
| CN104502707A (zh) | 一种基于三次样条插值的电力系统同步相量测量方法 | |
| CN102901856A (zh) | 基于相位搜索的电缆线路避雷器阻性电流检测方法 | |
| CN102809687B (zh) | 一种交流电频率的数字化测量方法 | |
| CN103018555A (zh) | 一种高精度的电力参数软件同步采样方法 | |
| CN108333426A (zh) | 基于傅氏算法的电力系统频率测量方法 | |
| WO2014008792A1 (zh) | 配网电压波动定位方法及定位装置 | |
| CN111625769A (zh) | 一种基于拉格朗日插值和三次指数平滑的pmu-scada数据对时与融合方法 | |
| CN105093163A (zh) | 一种电子式电能表的误差调试方法 | |
| CN104237633B (zh) | 一种频率快速测量方法 | |
| CN102122894B (zh) | 一种直流调速器同步信号检测方法 | |
| CN104330644B (zh) | 检测电力系统中正弦波信号幅值的方法和系统 | |
| CN103207307A (zh) | 金属氧化物避雷器泄漏电流组成成分测量方法 | |
| CN106291068B (zh) | 交流电机电流有效值的检测计算方法 | |
| CN103439569A (zh) | 一种宽范围量程自适应电压质量监测方法 | |
| CN112129989B (zh) | 一种基于自适应聚类的电压暂降分段刻画方法 | |
| CN102095936A (zh) | 一种快速测量交流电信号相位差的方法 | |
| 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 |