CN117054737B - 一种自相关滤波计算电力供电频率的方法和装置 - Google Patents
一种自相关滤波计算电力供电频率的方法和装置 Download PDFInfo
- Publication number
- CN117054737B CN117054737B CN202311026967.2A CN202311026967A CN117054737B CN 117054737 B CN117054737 B CN 117054737B CN 202311026967 A CN202311026967 A CN 202311026967A CN 117054737 B CN117054737 B CN 117054737B
- Authority
- CN
- China
- Prior art keywords
- frequency
- filtering
- period
- voltage
- voltage signal
- 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
- 238000001914 filtration Methods 0.000 title claims abstract description 21
- 238000000034 method Methods 0.000 title claims abstract description 17
- 238000005070 sampling Methods 0.000 claims description 19
- 238000007493 shaping process Methods 0.000 claims description 5
- 230000003111 delayed effect Effects 0.000 claims description 4
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 5
- 238000001514 detection method Methods 0.000 abstract description 3
- 238000007796 conventional method Methods 0.000 abstract description 2
- 238000012544 monitoring process Methods 0.000 abstract description 2
- 230000001934 delay Effects 0.000 abstract 1
- 235000019800 disodium phosphate Nutrition 0.000 abstract 1
- 238000005259 measurement Methods 0.000 description 5
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/02—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
- G01R19/25—Arrangements for measuring currents or voltages or for indicating presence or sign thereof using digital measurement techniques
- G01R19/2503—Arrangements for measuring currents or voltages or for indicating presence or sign thereof using digital measurement techniques for measuring voltage only, e.g. digital volt meters (DVM's)
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/40—Arrangements for reducing harmonics
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
本发明提出一种自相关滤波计算电力供电频率的方法和装置,属于用于各种电能质量检测设备和电力录波装置的自动监测方法。首先对AD采样的电压信号滤波,而后分别用延时m和2m的信号与自身相乘,乘得的数据经过低通滤波器后依公式求得当前的瞬时频率值。多个工频周期的瞬时频率值求平均可得到频率的精确值。本发明的优点是延时大大小于常规方法,计算精度高,适于各种计算机、DSP或FPGA实现。
Description
技术领域
本发明涉及电力***技术领域,具体的说是一种自相关滤波计算电力供电频率的方法和装置,用于电能质量检测设备的频率测量。
背景技术
随着电力电子技术的发展,电能质量监测在电力***中的应用越来越广泛,而频率测量是电能质量的重要组成部分。对电力频率的高精度测量是对电能质量其余很多参数测量的基础。
常用的频率测量方法有闸门法、周期法以及等精度法。但是这些方法有如下缺点:1、需要使用调解电路将电压转换为方波,然后使用定时器/计数器计算频率;2、不能提供随AD采样率逐点的频率信息;3、需要测量多个工频周期才可以得到足够的精度,瞬时测量的分辨率和精度不够理想。
随着计算机硬件水平的发展,计算复杂度对于算法的限制越来越小,因此可以考虑使用较复杂算法,将AD变换后的数值计算成为频率信息。
发明内容
针对现有技术中存在的上述不足之处,本发明要解决的技术问题是提供一种自相关滤波计算电力供电频率的方法和装置,用于各种电能质量检测设备和电力录波装置。
本发明为实现上述目的所采用的技术方案是:一种自相关滤波计算电力供电频率的方法,包括以下步骤:
以采样周期Ts采样被分析的电压信号得到V(n);
滤除V(n)的谐波和直流分量,得到V`(n);
将V`(n)分别延时mTs和2mTs后与V`(n)相乘,其中,m为整形常量;
滤除相乘结果中的高频分量;
根据滤除高频分量的相乘结果,获得瞬时频率;
根据瞬时频率获得一段时间内的电压频率。
所述将V`(n)分别延时mTs和2mTs后与V`(n)相乘,具体为:所述相乘S1(n)和S2(n)结果表达为:
S1(n)=V`(n)V`(n+m)
S2(n)=V`(n)V`(n+2m)
所述电压信号V(n)表示如下:
V(n)=Ucos(2π×n×Ts×f(n))
式中,U为基波电压,f(n)为瞬时频率;
将V(n)代入S1(n)和S2(n)中并三角分解,得:
式中,ρ=2π×Ts×f(n)。
所述滤除高频分量的相乘结果为:
所述瞬时频率为:
所述一段时间内的电压频率为:
其中,N为该段时间内对应的采样点数量。
一种自相关滤波计算电力供电频率的装置,包括:
采样模块,用于以采样周期Ts采样被分析的电压信号V(n);带通滤波器,用于滤除V(n)的谐波和直流分量,得到V`(n);
延时模块,用于将V`(n)分别延时mTs和2mTs,其中,m为整形常量;
乘法模块,用于将延时后的结果与V`(n)相乘;
低通滤波器,用于滤除相乘结果中的高频分量;
处理模块,用于根据滤除高频分量的相乘结果,获得瞬时频率;
确定模块,用于根据瞬时频率获得一段时间内的电压频率。
本发明具有以下优点及有益效果:
1、本发明频率计算的延时取决于AD采样周期和滤波器的阶数,一般在1到2个工频周期就可以得到很高的精度,计算延时大大小于常规方法。
2、计算精度高。
3、适于各种计算机、DSP或FPGA实现。
附图说明
图1为本发明实施例的方法流程图;
图2为本发明实施例的计算原理框图;
图3为本发明实施例的装置结构图。
具体实施方式
下面结合附图及实施例对本发明做进一步的详细说明。
如图1和图2所示,一种自相关滤波计算电力供电频率的方法,包括以下步骤:
步骤a.以采样周期为Ts采样被分析的电压信号得到V(n),n为电压采样序号,直接用分压器或从电压互感器PT二次侧取得电网的电压信号,经过抗混叠滤波器送达AD采样处。用带通滤波器滤除谐波和直流分量,得到V`(n),然后分别延时mTs和2mTs后和自己相乘,得:
S1(n)=V′(n)V′(n+m)
S2(n)=V′(n)V′(n+2m) (1)
其中,m为整形常量,一般2mTS对应1/8周期左右,实际上m>=1算法都成立。电压信号V(n)可表示如下:
V(n)=Ucos(2π×n×Ts×f(n)) (2)
式中U为基波电压,f(n)为瞬时频率,代入式(1)并三角分解,得:
式中,ρ=2π×Ts×f(n)。
步骤b.将S1(n)、S2(n)通过低通滤波器滤除高频分量,得:
步骤c.解该方程,得
在延时小于π/2时,加号成立,则可解得频率的瞬时值:
式中,f(n)为瞬时频率。
步骤d.瞬时频率在一段时间的算数平均就是频率的精确值。在一段时间对瞬时频率求平均,得到频率的精确值f:
其中,N为该段时间内对应的采样点数量,n=0~N-1。
如图3所示,一种自相关滤波计算电力供电频率的装置,包括:
采样模块,用于以采样周期Ts采样被分析的电压信号V(n);
带通滤波器,用于滤除V(n)的谐波和直流分量,得到V`(n),该带通滤波器截止频率在90Hz;
延时模块,用于将V`(n)分别延时mTs和2mTs,其中,m为整形常量;
乘法模块,用于将延时后的结果与V`(n)相乘;
低通滤波器,用于滤除相乘结果中的高频分量,该低通滤波器的截止频率为45Hz;
处理模块,用于根据滤除高频分量的相乘结果,获得瞬时频率;
确定模块,用于根据瞬时频率获得一段时间内的电压频率。
Claims (2)
1.一种自相关滤波计算电力供电频率的方法,其特征在于,包括以下步骤:
以采样周期Ts采样被分析的电压信号得到V(n);
滤除V(n)的谐波和直流分量,得到V`(n);
将V`(n)分别延时mTs和2mTs后与V`(n)相乘,其中,m为整形常量;
滤除相乘结果中的高频分量;
根据滤除高频分量的相乘结果,获得瞬时频率;
根据瞬时频率获得一段时间内的电压频率;
所述将V`(n)分别延时mTs和2mTs后与V`(n)相乘,具体为:
所述相乘S1(n)和S2(n)结果表达为:
S1(n)=V`(n)V`(n+m)
S2(n)=V`(n)V`(n+2m)
所述电压信号V(n)表示如下:
V(n)=Ucos(2π×n×Ts×f(n))
式中,U为基波电压,f(n)为瞬时频率;
将V(n)代入S1(n)和S2(n)中并三角分解,得:
式中,ρ=2π×Ts×f(n);
所述滤除高频分量的相乘结果为:
所述瞬时频率为:
2.根据权利要求1所述的一种自相关滤波计算电力供电频率的方法,其特征在于,所述一段时间内的电压频率为:
其中,N为该段时间内对应的采样点数量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311026967.2A CN117054737B (zh) | 2023-08-15 | 2023-08-15 | 一种自相关滤波计算电力供电频率的方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311026967.2A CN117054737B (zh) | 2023-08-15 | 2023-08-15 | 一种自相关滤波计算电力供电频率的方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117054737A CN117054737A (zh) | 2023-11-14 |
CN117054737B true CN117054737B (zh) | 2024-03-29 |
Family
ID=88668728
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311026967.2A Active CN117054737B (zh) | 2023-08-15 | 2023-08-15 | 一种自相关滤波计算电力供电频率的方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117054737B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4433422A (en) * | 1982-03-29 | 1984-02-21 | Sperry Corporation | Frequency measuring system for alternating frequency signals |
US5471133A (en) * | 1992-04-06 | 1995-11-28 | Siemens Aktiengesellschaft | Digital device for measuring the frequency of an electrical signal |
JPH09119947A (ja) * | 1995-10-26 | 1997-05-06 | Railway Technical Res Inst | 特定周波数成分信号の予測方法および装置並びにアクティブフィルタのデジタル制御システム |
US5812737A (en) * | 1995-01-09 | 1998-09-22 | The Board Of Trustees Of The Leland Stanford Junior University | Harmonic and frequency-locked loop pitch tracker and sound separation system |
CN102075269A (zh) * | 2010-12-24 | 2011-05-25 | 北京航空航天大学 | 高速跳频信号的跳速和频率转换时间的测量方法 |
CN102347926A (zh) * | 2011-09-26 | 2012-02-08 | 豪威科技(上海)有限公司 | 载波频率捕获方法及装置 |
JP2012163543A (ja) * | 2010-09-30 | 2012-08-30 | Daihen Corp | 周波数検出装置 |
CN102983829A (zh) * | 2012-11-02 | 2013-03-20 | 南昌航空大学 | 基于电域延迟自相关的激光器频率差锁定方法 |
CN105044456A (zh) * | 2015-07-21 | 2015-11-11 | 电子科技大学 | 一种基于正交子带的电网瞬时频率测量与跟踪方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8138969B2 (en) * | 2008-10-22 | 2012-03-20 | Bae Systems Information And Electronic Systems Integration Inc. | Monobit based low cost high performance radar warning receiver |
US9553620B2 (en) * | 2014-07-16 | 2017-01-24 | Raytheon Company | Signal detection and characterization |
-
2023
- 2023-08-15 CN CN202311026967.2A patent/CN117054737B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4433422A (en) * | 1982-03-29 | 1984-02-21 | Sperry Corporation | Frequency measuring system for alternating frequency signals |
US5471133A (en) * | 1992-04-06 | 1995-11-28 | Siemens Aktiengesellschaft | Digital device for measuring the frequency of an electrical signal |
US5812737A (en) * | 1995-01-09 | 1998-09-22 | The Board Of Trustees Of The Leland Stanford Junior University | Harmonic and frequency-locked loop pitch tracker and sound separation system |
JPH09119947A (ja) * | 1995-10-26 | 1997-05-06 | Railway Technical Res Inst | 特定周波数成分信号の予測方法および装置並びにアクティブフィルタのデジタル制御システム |
JP2012163543A (ja) * | 2010-09-30 | 2012-08-30 | Daihen Corp | 周波数検出装置 |
CN102075269A (zh) * | 2010-12-24 | 2011-05-25 | 北京航空航天大学 | 高速跳频信号的跳速和频率转换时间的测量方法 |
CN102347926A (zh) * | 2011-09-26 | 2012-02-08 | 豪威科技(上海)有限公司 | 载波频率捕获方法及装置 |
CN102983829A (zh) * | 2012-11-02 | 2013-03-20 | 南昌航空大学 | 基于电域延迟自相关的激光器频率差锁定方法 |
CN105044456A (zh) * | 2015-07-21 | 2015-11-11 | 电子科技大学 | 一种基于正交子带的电网瞬时频率测量与跟踪方法 |
Non-Patent Citations (3)
Title |
---|
Transmission Gates Combined With Level-Restoring CMOS Gates Reduce Glitches in Low-Power Low-Frequency Multipliers;Flavio Carbognani, et al.;IEEE Transactions on Very Large Scale Integration (VLSI) Systems;20080731;第16卷(第7期);第830-836页 * |
一种新的载波频率估计方法;钟志 等;哈尔滨商业大学学报(自然科学版);20070615;第23卷(第03期);第340-342页 * |
两种简单实用的电网频率实时测量方法;林的方 等;电网与清洁能源;20090425;第25卷(第04期);第37-39页 * |
Also Published As
Publication number | Publication date |
---|---|
CN117054737A (zh) | 2023-11-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zygarlicki et al. | A reduced Prony's method in power-quality analysis—parameters selection | |
CN102116798A (zh) | 一种电网频率测量方法及装置 | |
CN102879639A (zh) | 一种电力***中频率的实时测量方法 | |
CN108318852B (zh) | 一种智能电能表方波影响试验方法 | |
CN103941088A (zh) | 一种基于三相信号的电力***频率快速测量方法 | |
CN104122444A (zh) | 全数字中频频谱分析仪及频谱分析方法 | |
CN103983849A (zh) | 一种实时高精度的电力谐波分析方法 | |
CN102495250A (zh) | 一种基于Hilbert算法的准同步宽频无功电能表及其采样方法 | |
CN110837000B (zh) | 基于fpga的频率测量*** | |
CN104155521A (zh) | 相位差的确定方法和装置 | |
CN114252837B (zh) | 用于数字功率计谐波测量的高精度均等取样***及方法 | |
RU2406094C2 (ru) | Способ мгновенного определения коэффициента искажения сигналов в электрической сети переменного тока и соответствующее устройство | |
CN103969508A (zh) | 一种实时高精密的电力谐波分析方法及装置 | |
CN117169590B (zh) | 一种基于软件变采样率的电力谐波分析的方法和装置 | |
CN106154034A (zh) | 一种功率模块的谐波测量装置和方法 | |
CN117054737B (zh) | 一种自相关滤波计算电力供电频率的方法和装置 | |
CN104330644B (zh) | 检测电力***中正弦波信号幅值的方法和*** | |
Tiwari et al. | Measurement of instantaneous power quality parameters using UWPT and Hilbert transform and its FPGA implementation | |
Valenzuela et al. | Real-time interharmonics detection and measurement based on FFT algorithm | |
US8433743B2 (en) | Root mean square (RMS) metering devices and methods for generating RMS current level to both high or low frequency within signal | |
CN107515332B (zh) | 一种基于频谱分析和同步采样的直流电能计量装置及方法 | |
CN109975598B (zh) | 一种周期信号有效值测量电路 | |
CN112557781A (zh) | 一种适用于校准器的宽频域信号测量方法 | |
CN109884388B (zh) | 一种基于半周期移相法的电网频率计量、测量装置及方法 | |
Bhatti et al. | Harmonic detection using microcontroller |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |