CN113466552B - 一种定间隔采样下频率跟踪方法 - Google Patents
一种定间隔采样下频率跟踪方法 Download PDFInfo
- Publication number
- CN113466552B CN113466552B CN202110793568.3A CN202110793568A CN113466552B CN 113466552 B CN113466552 B CN 113466552B CN 202110793568 A CN202110793568 A CN 202110793568A CN 113466552 B CN113466552 B CN 113466552B
- Authority
- CN
- China
- Prior art keywords
- index
- data
- sampling
- frequency
- points
- 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
- 238000005070 sampling Methods 0.000 title claims abstract description 43
- 238000000034 method Methods 0.000 title claims abstract description 20
- 238000004364 calculation method Methods 0.000 claims abstract description 16
- 238000005259 measurement Methods 0.000 claims abstract description 8
- 238000004422 calculation algorithm Methods 0.000 abstract description 7
- 230000001360 synchronised effect Effects 0.000 abstract description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000758 substrate Substances 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R25/00—Arrangements for measuring phase angle between a voltage and a current or between voltages or currents
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种定间隔采样下频率跟踪方法,包括如下步骤:(1)装置AD过采样;(2)根据装置AD采样频率、电网实时频率计算调整系数;(3)提取最新两周波数据,并将提取数据前后对倒;(4)将步骤(3)中的数据,基于拉格朗日插值函数推导FFT计算数据,推导数据用T表示;(5)将步骤(4)中推导出的数据对倒;(6)将步骤(5)中提取的数据进行傅里叶FFT处理即可得到系统的测量量。本发明基于拉格朗日插值算法的采样值调整方法,根据系统频率将采样到的数据通过拉格朗日插值算法调整到整周期,解决系统频率偏差时傅里叶计算准确性问题,对于多设备采样同步或者无法进行硬件频率跟踪的场景,具有良好的表现。
Description
技术领域
本发明涉及数字信号处理技术领域,尤其是一种定间隔采样下频率跟踪方法。
背景技术
电力系统中,继电保护装置通过AD芯片将系统中电压电流量转换成离散的数字信号采集到CPU中,数字量经过快速傅里叶(FFT)处理得到系统电压电流的基波幅值、相位、谐波含量等测量量,其测量精度的是衡量继电保护装置性能的重要指标。当电力系统频率偏离50Hz时,在固定采样间隔条件下,对被测信号采样量很难做到整周期截断,FFT算法的泄漏现象和栅栏效应将会导致测量量出现较大的测量误差。例如固定32点采样,由于系统频率偏差,采样到的32点数据多于或者少于一个周波,导致傅里叶计算出现较大误差。为解决此问题,现有设备中通常会应用一种频率跟踪技术,根据检测到的电网频率实时调整AD采样间隔,保证每周波AD采样点数相同,以应对系统频率飘动时傅里叶计算准确性问题。但是对于多设备(设备大于等于2)采样同步的场景中实时调整AD采样间隔将会降低装置间同步准确性或提高设备同步性能要求,同时有些设备无法实时调整AD采样间隔,例如多间隔DTU、不同间隔信号用同一AD芯片采集时。
发明内容
本发明所要解决的技术问题在于,提供一种定间隔采样下频率跟踪方法,基于拉格朗日插值算法的采样值调整方法,根据系统频率将采样到的数据通过拉格朗日插值算法调整到整周期,解决系统频率偏差时傅里叶计算准确性问题。
为解决上述技术问题,本发明提供一种,包括如下步骤:
(1)装置AD过采样;
(2)根据装置AD采样频率、电网实时频率计算调整系数;
(3)提取最新两周波数据,并将提取数据前后对倒;
(4)将步骤(3)中的数据,基于拉格朗日插值函数推导推导FFT计算数据,推导数据用T表示;
(5)将步骤(4)中推导出的数据对倒;
(6)将步骤(5)中提取的数据进行傅里叶FFT处理即可得到系统电压电流的基波幅值、相位、谐波含量的测量量。
优选的,步骤(1)中,装置AD过采样具体为:FFT采用32点计算时,AD采样速率采用其整数倍,倍数用M表示,选用32、64、96、128点等采样速率。理论上AD采样速率倍数越高,准确度越高。
优选的,步骤(2)中,根据装置采样频率、电网实时频率计算调整系数K具体为:
FCaiYang为装置采样频率,FXiTong为电网实时频率,通过装置频率采集回路计算得到。
优选的,步骤(3)中,提取最新2周波数据,并将提取数据前后对倒具体为:采样频率50HZ下固定128点采样,计采样间隔时间τ,提取256点数据,数据计为x1、x2、x3...x256,即计为(0,x1)(τ,x2)(2τ,x3)...(255τ,x256),256τ等于1/FCaiYang;
将前后对倒,即x1、x256对倒,x2、x255对倒,依次类推,计为(0,x256)(τ,x255)(2τ,x254)...(255τ,x1)。
优选的,步骤(4)中,将步骤(3)中的数据,基于拉格朗日插值函数推导FFT计算数据,推导数据用T表示具体为:32点傅里叶计算,即计算推导出T1、T2、T3...T32;推导点依次为(0,T1)(Mτ,T2)(2Mτ,T3)(3Mτ,T4)...(31Mτ,T32),31Mτ等于256τ;
推导数据第一点为提取数据第一点T1=x1,即(0,T1)=(0,x1);
对于第2到32点数据通过如下公式推导,n为对应的2到32点点号:
index_f=K*n*M,K为步骤(2)计算出调整系数,M为步骤(1)中倍数,若M等于4;
index_n=[index_f],对index_f取整;
将(index_n-1,x[index_n])、(index_n,x[index_n+1])、
(index_n+1,x[index_n+2])三点建立拉格朗日插值函数,即可计算出Tn:
Tn=x[index_n]*(index_n+1-index_f)*(index_n+2-index_f)/2+x[index_n+1]*(index_n-index_f)*(index_n+2-index_f)-x[index_n+2]*(index_n-index_f)*(index_n+1-index_f)/2。
优选的,步骤(5)中,将步骤(4)中推导出的数据对倒具体为:T1与T32对倒,T2与T31对倒,以此类推,即可完成整数据窗数据推导截取即频率跟踪处理,即(0,T32)(Mτ,T31)(2Mτ,T30)...(31Mτ,T1)。
本发明的有益效果为:本发明基于拉格朗日插值算法的采样值调整方法,根据系统频率将采样到的数据通过拉格朗日插值算法调整到整周期,解决系统频率偏差时傅里叶计算准确性问题,对于多设备采样同步、或者无法进行硬件频率跟踪的场景,具有良好的表现。
具体实施方式
本发明以系统频率50Hz,装置固定每周期采样128点,间隔4点提取32点数据,经过快速傅里叶(FFT)处理得到系统电压电流的基波幅值、相位、谐波含量等测量量为例,具体实施如下:
步骤1:装置AD过采样,即FFT采用32点计算,AD采样速率是其4倍128点采样;
步骤2:根据装置采样频率、电网实时频率计算系数K:
FCaiYang为装置采样频率,采样频率50HZ;
FXiTong为电网实时频率,通过装置频率采集回路计算得到。
步骤3:提取最新2周波数据,并将提取数据前后对倒。采样间隔时间τ,提取256点数据,数据计为x1、x2、x3...x256,即计为(0,x1)(τ,x2)(2τ,x3)...(255τ,x256),256τ等于1/50ms。
将前后数据对倒,即x1、x256对倒,x2、x255对倒,依次类推,计为(0,x256)(τ,x255)(2τ,x254)...(255τ,x1);
步骤4:将步骤3中数据应用本步骤中方法推导出FFT计算数据,推导数据用T表示,本发明以32点傅里叶计算为例,即计算推导出T1、T2、T3...T32;推导点依次为(0,T1)(Mτ,T2)(2Mτ,T3)(3Mτ,T4)...(31Mτ,T32),31Mτ等于256τ。
推导数据第一点为提取数据第一点T1=x1,即(0,T1)=(0,x1);
对于第2到32点数据通过如下公式推导,n为对应的2到32点数据:
index_f=K*n*M,K为步骤2计算系数,M等于4;
index_n=[index_f],对index_f取整;
将(index_n-1,x[index_n])、(index_n,x[index_n+1])、
(index_n+1,x[index_n+2])三点建立拉格朗日插值函数,求出截取点Tn:
Tn=x[index_n]*(index_n+1-index_f)*(index_n+2-index_f)/2+x[index_n+1]*(index_n-index_f)*(index_n+2-index_f)-x[index_n+2]*(index_n-index_f)*(index_n+1-index_f)/2;
步骤5:将步骤4推导出数据对倒,即T1与T32对倒,T2与T31对倒,以此类推,即可完成本发明专利整数据窗数据推导即频率跟踪处理,即(0,T32)(Mτ,T31)(2Mτ,T30)...(31Mτ,T1);
步骤6:将步骤5提取数据进行傅里叶(FFT)处理即可得到系统电压电流的基波幅值、相位、谐波含量等测量量。
Claims (2)
1.一种定间隔采样下频率跟踪方法,其特征在于,包括如下步骤:
(1)装置AD过采样;
(2)根据装置AD采样频率、电网实时频率计算调整系数;根据装置采样频率、电网实时频率计算调整系数K具体为:
FCaiYang为装置AD采样频率,FXiTong为电网实时频率,通过装置频率采集回路计算得到;
(3)提取最新两周波数据,并将提取数据前后对倒;具体为:采样频率50HZ下固定128点采样,计采样间隔时间τ,提取256点数据,数据计为x1、x2、x3...x256,即计为(0,x1)(τ,x2)(2τ,x3)...(255τ,x256),256τ等于1/FCaiYang;
将前后对倒,即x1、x256对倒,x2、x255对倒,依次类推,计为(0,x256)(τ,x255)(2τ,x254)...(255τ,x1);
(4)将步骤(3)中的数据,基于拉格朗日插值函数推导FFT计算数据,推导数据用T表示;具体为:32点傅里叶计算,即计算推导出T1、T2、T3...T32;推导点依次为(0,T1)(Mτ,T2)(2Mτ,T3)(3Mτ,T4)...(31Mτ,T32),31Mτ等于256τ;
推导数据第一点为提取数据第一点T1=x1,即(0,T1)=(0,x1);
对于第2到32点数据通过如下公式推导,n为对应的2到32点点号:
index_f=K*n*M,K为步骤(2)计算出调整系数,M为步骤(1)中倍数,
若M等于4;
index_n=[index_f],对index_f取整;
将(index_n-1,x[index_n])、(index_n,x[index_n+1])、
(index_n+1,x[index_n+2])三点建立拉格朗日插值函数,即可计算出Tn:
Tn=x[index_n]*(index_n+1-index_f)*(index_n+2-index_f)/2
+x[index_n+1]*(index_n-index_f)*(index_n+2-index_f)
-x[index_n+2]*(index_n-index_f)*(index_n+1-index_f)/2;
(5)将步骤(4)中推导出的数据对倒;具体为:T1与T32对倒,T2与T31对倒,以此类推,即可完成整数据窗数据推截取即频率跟踪处理,即(0,T32)(Mτ,T31)(2Mτ,T30)...(31Mτ,T1);
(6)将步骤(5)中提取的数据进行傅里叶FFT处理即可得到系统电压电流的基波幅值、相位、谐波含量的测量量。
2.如权利要求1所述的定间隔采样下频率跟踪方法,其特征在于,步骤(1)中,装置AD过采样具体为:FFT采用32点计算时,AD采样速率采用其整数倍,倍数用M表示,选用32、64、96、128点采样速率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110793568.3A CN113466552B (zh) | 2021-07-14 | 2021-07-14 | 一种定间隔采样下频率跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110793568.3A CN113466552B (zh) | 2021-07-14 | 2021-07-14 | 一种定间隔采样下频率跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113466552A CN113466552A (zh) | 2021-10-01 |
CN113466552B true CN113466552B (zh) | 2024-02-02 |
Family
ID=77880178
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110793568.3A Active CN113466552B (zh) | 2021-07-14 | 2021-07-14 | 一种定间隔采样下频率跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113466552B (zh) |
Citations (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1629873A1 (ru) * | 1988-02-15 | 1991-02-23 | Харьковский Институт Радиоэлектроники Им.Акад.М.К.Янгеля | Способ определени коэффициента гармоник |
JP2007171833A (ja) * | 2005-12-26 | 2007-07-05 | Kenwood Corp | 電子透かしエンコード装置、電子透かしデコード装置、電子透かしエンコード方法、電子透かしデコード方法及びプログラム |
KR20110046754A (ko) * | 2009-10-29 | 2011-05-06 | 한국전기연구원 | 재샘플링 기법이 적용된 시각 동기화 페이저 측정장치 |
CN102338827A (zh) * | 2011-06-10 | 2012-02-01 | 中国矿业大学 | 用于电网波形畸变分析和电力谐波参数自动监测方法 |
CN102435844A (zh) * | 2011-11-01 | 2012-05-02 | 南京磐能电力科技股份有限公司 | 一种频率无关的正弦信号相量计算方法 |
CN102801557A (zh) * | 2012-07-27 | 2012-11-28 | 中国电力科学研究院 | 一种提高变电站网络采样可靠性的系统方法 |
CN103196547A (zh) * | 2013-03-11 | 2013-07-10 | 安徽新力电业科技咨询有限责任公司 | 一种实现旋转机械振动信号同步阶比跟踪分析方法 |
CN103324603A (zh) * | 2013-06-07 | 2013-09-25 | 山东山大电力技术有限公司 | 一种不同采样速率的插值变换方法 |
CN103728619A (zh) * | 2014-01-16 | 2014-04-16 | 北京航空航天大学 | 基于变重频技术的机载大斜视条带sar成像方法 |
CN104934032A (zh) * | 2014-03-17 | 2015-09-23 | 华为技术有限公司 | 根据频域能量对语音信号进行处理的方法和装置 |
CN105004913A (zh) * | 2015-07-21 | 2015-10-28 | 黎文安 | 一种电力系统交流采样频率实时跟踪方法 |
CN105261202A (zh) * | 2015-11-05 | 2016-01-20 | 江苏省电力公司电力科学研究院 | 基于b样条插值小波突变点检测的传输延时测量方法 |
CN106597095A (zh) * | 2016-12-22 | 2017-04-26 | 南京因泰莱电器股份有限公司 | 一种双cpu等间隔采样的频率跟踪实现方法 |
CN108072868A (zh) * | 2016-11-11 | 2018-05-25 | 中国科学院沈阳自动化研究所 | 一种基于fmcw雷达信号频率细化的高精度测距方法 |
CN109633266A (zh) * | 2019-02-26 | 2019-04-16 | 重庆新世杰电气股份有限公司 | 一种频率测量方法、系统、装置及计算机可读存储介质 |
WO2020220804A1 (zh) * | 2019-04-28 | 2020-11-05 | 哈尔滨工业大学 | 一种可变采样率的傅里叶分析方法 |
CN112526209A (zh) * | 2020-10-20 | 2021-03-19 | 江苏宝亨新电气有限公司 | 电力系统同步相量测量方法 |
-
2021
- 2021-07-14 CN CN202110793568.3A patent/CN113466552B/zh active Active
Patent Citations (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1629873A1 (ru) * | 1988-02-15 | 1991-02-23 | Харьковский Институт Радиоэлектроники Им.Акад.М.К.Янгеля | Способ определени коэффициента гармоник |
JP2007171833A (ja) * | 2005-12-26 | 2007-07-05 | Kenwood Corp | 電子透かしエンコード装置、電子透かしデコード装置、電子透かしエンコード方法、電子透かしデコード方法及びプログラム |
KR20110046754A (ko) * | 2009-10-29 | 2011-05-06 | 한국전기연구원 | 재샘플링 기법이 적용된 시각 동기화 페이저 측정장치 |
CN102338827A (zh) * | 2011-06-10 | 2012-02-01 | 中国矿业大学 | 用于电网波形畸变分析和电力谐波参数自动监测方法 |
CN102435844A (zh) * | 2011-11-01 | 2012-05-02 | 南京磐能电力科技股份有限公司 | 一种频率无关的正弦信号相量计算方法 |
CN102801557A (zh) * | 2012-07-27 | 2012-11-28 | 中国电力科学研究院 | 一种提高变电站网络采样可靠性的系统方法 |
CN103196547A (zh) * | 2013-03-11 | 2013-07-10 | 安徽新力电业科技咨询有限责任公司 | 一种实现旋转机械振动信号同步阶比跟踪分析方法 |
CN103324603A (zh) * | 2013-06-07 | 2013-09-25 | 山东山大电力技术有限公司 | 一种不同采样速率的插值变换方法 |
CN103728619A (zh) * | 2014-01-16 | 2014-04-16 | 北京航空航天大学 | 基于变重频技术的机载大斜视条带sar成像方法 |
CN104934032A (zh) * | 2014-03-17 | 2015-09-23 | 华为技术有限公司 | 根据频域能量对语音信号进行处理的方法和装置 |
CN105004913A (zh) * | 2015-07-21 | 2015-10-28 | 黎文安 | 一种电力系统交流采样频率实时跟踪方法 |
CN105261202A (zh) * | 2015-11-05 | 2016-01-20 | 江苏省电力公司电力科学研究院 | 基于b样条插值小波突变点检测的传输延时测量方法 |
CN108072868A (zh) * | 2016-11-11 | 2018-05-25 | 中国科学院沈阳自动化研究所 | 一种基于fmcw雷达信号频率细化的高精度测距方法 |
CN106597095A (zh) * | 2016-12-22 | 2017-04-26 | 南京因泰莱电器股份有限公司 | 一种双cpu等间隔采样的频率跟踪实现方法 |
CN109633266A (zh) * | 2019-02-26 | 2019-04-16 | 重庆新世杰电气股份有限公司 | 一种频率测量方法、系统、装置及计算机可读存储介质 |
WO2020220804A1 (zh) * | 2019-04-28 | 2020-11-05 | 哈尔滨工业大学 | 一种可变采样率的傅里叶分析方法 |
CN112526209A (zh) * | 2020-10-20 | 2021-03-19 | 江苏宝亨新电气有限公司 | 电力系统同步相量测量方法 |
Non-Patent Citations (5)
Title |
---|
Closed-form reconstruction of images from irregular 2-D discrete Fourier samples using the Good-Thomas FFT;A.E. Yagle等;Proceedings 2000 International Conference on Image Processing (Cat. No.00CH37101);全文 * |
Sample Rate Converter As a Means of Reducing Measurment Error of the Voltage Spectrum by Application of FFT;Andrey N. Serov等;2019 29th International Conference Radioelektronika (RADIOELEKTRONIKA);全文 * |
一种高精度加窗插值FFT谐波分析方法;汪立;刘翠琳;;陕西电力(第05期);全文 * |
基于插值同步算法的频谱泄漏分析与仿真;曾泽昊;余有灵;许维胜;;系统仿真技术(第03期);全文 * |
用于船舶电站控制器的交流采样算法研究;朱翟;唐石青;王良秀;殷非;;船舶工程(第06期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113466552A (zh) | 2021-10-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Offelli et al. | The influence of windowing on the accuracy of multifrequency signal parameter estimation | |
US7480580B2 (en) | Apparatus and method for estimating synchronized phasors at predetermined times referenced to an absolute time standard in an electrical system | |
CN106209701B (zh) | Alpha稳定分布噪声环境下MFSK信号符号速率估计方法及装置 | |
CN102338827B (zh) | 用于电网波形畸变分析和电力谐波参数自动监测方法 | |
WO2003090327A2 (en) | Protective relay with synchronized phasor measurement capability for use in electric power systems | |
CN108614155B (zh) | 一种加入汉明窗的同步相量测量方法及系统 | |
US20110077885A1 (en) | apparatus and method for estimating synchronized phasors at predetermined times referenced to a common time standard in an electrical system | |
CN103995181B (zh) | 数字化变电站电能质量谐波分析方法 | |
CN102331526A (zh) | 汉宁窗函数连续频谱内插电力谐波参数获取方法 | |
Mayo Jr | Spectrum measurements with laser velocimeters | |
CN103543335A (zh) | 一种同步相量测量方法 | |
CN110837003A (zh) | 一种基于三角窗的双窗全相位dft同步相量测量方法及系统 | |
CN102253282A (zh) | 纳托尔窗函数连续频谱内插电力谐波参数获取方法 | |
CN105137181A (zh) | 基于Nuttall-Kaiser组合窗双谱线插值的谐波分析方法 | |
CN113466552B (zh) | 一种定间隔采样下频率跟踪方法 | |
Lesher et al. | SUNBURST/spl minus/a network of GIC monitoring systems | |
CN103605904B (zh) | 基于误差估算的自补偿电力系统幅值算法 | |
Hayakawa et al. | Recent findings on VLF/ELF sferics | |
Bhatti et al. | Implementation of low cost non-DFT based Phasor Measurement Unit for 50 Hz power system | |
Maharjan et al. | Improved off-nominal operation of phasor measurement units using discrete fourier transformation | |
CN111697541A (zh) | 一种配电网线路电流差动保护同步方法及系统 | |
CN116582400A (zh) | 基于短时脉冲干扰的l1载波恢复方法、装置及智能终端 | |
CN107870338A (zh) | 一种低更新频度的卫星导航载波跟踪方法 | |
De Waele et al. | A time domain error measure for resampled irregular data | |
CN109709422B (zh) | 一种时钟漂移的消除方法及其装置 |
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 |