CN109030942B - 谐相角分析方法 - Google Patents

谐相角分析方法 Download PDF

Info

Publication number
CN109030942B
CN109030942B CN201811014781.4A CN201811014781A CN109030942B CN 109030942 B CN109030942 B CN 109030942B CN 201811014781 A CN201811014781 A CN 201811014781A CN 109030942 B CN109030942 B CN 109030942B
Authority
CN
China
Prior art keywords
sampling
analysis
harmonic
data
sampling point
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
Application number
CN201811014781.4A
Other languages
English (en)
Other versions
CN109030942A (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.)
Jiangsu University of Technology
Original Assignee
Jiangsu University of Technology
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 Jiangsu University of Technology filed Critical Jiangsu University of Technology
Priority to CN201811014781.4A priority Critical patent/CN109030942B/zh
Publication of CN109030942A publication Critical patent/CN109030942A/zh
Application granted granted Critical
Publication of CN109030942B publication Critical patent/CN109030942B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measuring Frequencies, Analyzing Spectra (AREA)
  • Measuring Phase Differences (AREA)
  • Complex Calculations (AREA)

Abstract

本发明涉及一种谐相角分析方法,包括如下几个步骤:(1)采样W+2个采样点数据(W由积分方法决定);(2)从采样点i=0应用准同步DFT公式
Figure 100004_DEST_PATH_IMAGE002
分析W+1个数据获得基波信息
Figure 100004_DEST_PATH_IMAGE004
Figure 100004_DEST_PATH_IMAGE006
;(3)从采样点i=1应用准同步DFT公式分析W+1个数据获得基波信息
Figure 100004_DEST_PATH_IMAGE008
Figure 100004_DEST_PATH_IMAGE010
;(4)应用公式
Figure 100004_DEST_PATH_IMAGE012
计算信号的频率漂移
Figure 100004_DEST_PATH_IMAGE014
;(5)从采样点i=0应用准同步DFT公式分析W+1个数据获得各次谐波信息
Figure 100004_DEST_PATH_IMAGE016
Figure 100004_DEST_PATH_IMAGE018
;(6)应用公式
Figure 100004_DEST_PATH_IMAGE020
计算各次谐波的幅相角;(7)应用公式
Figure 100004_DEST_PATH_IMAGE022
线性修正各次谐波的谐相角。本发明的方法,有助于电能质量监控、电子产品生产检验、电器设备监控等应用谐波分析的领域更加精确的获得各次谐波的幅相角和频率等信息。

Description

谐相角分析方法
本申请是申请号为:201510258020.3,发明创造名称为《一种谐相角分析方法》,申请日为:2015年05月19日的发明专利申请的分案申请。
技术领域
本发明涉及一种高精度的谐相角分析方法。
背景技术
谐波分析技术在电能质量监控、电子产品生产检验、电器设备监控等众多领域应用广泛,是进行电网监控、质量检验、设备监控的重要技术手段。目前谐波分析应用最广泛的技术是离散傅里叶变换(DFT)和快速傅里叶变换(FFT)。准同步采样技术和DFT技术相结合的谐波分析技术能够提高谐波分析的精度,其算式为:
Figure DEST_PATH_IMAGE002
式中:k为需要获得的谐波的次数(如基波k=1,3次谐波k=3);sin和cos分别为正弦和余弦函数;而a k b k 分别为k次谐波的实部和虚部;n为迭代次数;W由积分方法决定,采用复化梯形积分方法时,W=nN
Figure DEST_PATH_IMAGE004
为一次加权系数;
Figure DEST_PATH_IMAGE006
,为所有加权系数之和;f(i)为分析波形的第i个采样值;N为周期内采样次数。
在工程应用中,谐波分析总是进行有限点的采样和难以做到严格意义的同步采样。这样,在应用准同步DFT进行谐波分析时,就会存在由于截断效应导致的长范围泄漏和由于栅栏效应导致的短范围泄漏,使得分析结果精度不高,甚至不可信。
图1给出了应用准同步DFT对于任一给定实例进行谐波分析的误差图。从图中可以发现,准同步DFT算法的谐相角除了50Hz时其余均误差极大,基本不可信。
发明内容
本发明要解决的技术问题是提供一种高精度的谐相角分析方法,以有效改进准同步DFT谐波分析技术的分析误差,获得高精度的谐波分析结果,从而提高基于谐波分析理论的电能质量监控、电子产品生产检验、电器设备监控等领域仪器设备的质量和状态判断的有效性。
实现本发明目的的技术方案是提供一种谐相角分析方法,包括以下步骤:
(1)等间距采样W+2个采样点数据{f(i),i=0,1,…,W+1}(W由所选择的积分方法决定,本发明并不指定某一种积分方法,常用的积分方法有复化梯形积分方法W=nN、复化矩形积分方法W=n(N-1)、复化辛普森积分方法W=n(N-1)/2等,可以根据本发明应用的实际情况来选择合适的积分方法。一般以复化梯形积分方法效果较理想。N为一个理想周期内的采样点数。
(2)从采样点i=0开始应用准同步DFT公式
Figure DEST_PATH_IMAGE008
,分析W+1个数据获得基波信息
Figure DEST_PATH_IMAGE010
Figure DEST_PATH_IMAGE012
(3)从采样点i=1应用准同步DFT公式分析W+1个数据获得基波信息
Figure DEST_PATH_IMAGE014
Figure DEST_PATH_IMAGE016
(4)应用公式:
Figure DEST_PATH_IMAGE018
,计算信号的频率漂移
Figure DEST_PATH_IMAGE020
(5)从采样点i=0开始应用准同步DFT公式分析W+1个数据获得各次谐波信息
Figure DEST_PATH_IMAGE022
Figure DEST_PATH_IMAGE024
(6)应用公式
Figure DEST_PATH_IMAGE026
计算各次谐波的幅相角;
(7)应用公式
Figure DEST_PATH_IMAGE028
线性修正各次谐波的谐相角。
等间隔采样是根据进行谐波分析的理想信号的周期T和频率f,在一个周期内采样N点,即采样频率为f s =Nf,且N≥64
所述的采样W+2个采样点数据是根据所选择的积分方法而作相应选择,若采用复化梯形积分方法,则W=nN;若采用复化矩形积分方法,则W=n(N-1);若采用复化辛普森积分方法,则W=n(N-1)/2。然后根据采样频率f s =Nf,获得采样点数据序列{f(i),i=0,1,…,W+ 1}n为迭代次数,一般n≥3;最后对该数据序列进行谐波分析。
一次迭代系数
Figure 151241DEST_PATH_IMAGE004
由积分方法、理想周期采样点N和迭代次数n决定,具体推导过程参见文献【戴先中.准同步采样应用中的若干问题[J].电测与仪表, 1988, (2): 2-7.】。
Figure 271644DEST_PATH_IMAGE006
为所有加权系数之和。
信号频率的漂移
Figure DEST_PATH_IMAGE030
是根据相邻采样点基波相角差与理想周期内采样点数N的固定关系而获得的,信号频率的漂移
Figure 989065DEST_PATH_IMAGE030
也可用于修正基波和高次谐波的频率f 1 和高次谐波的频率f k
Figure DEST_PATH_IMAGE032
)。
本发明具有积极的效果:(1)高精度的谐相角分析结果。如对于图1给定的分析实例,本发明获得的分析精度提高到10-8级(图2)。
(2)本发明所述的方法从根本上解决了准同步DFT谐相角分析精度低的问题,而无需进行复杂的反演和修正,算法简单。
(3)相对于准同步DFT,本发明所述的谐波分析技术只需要增加一个采样点就解决了准同步DFT分析误差大的问题,易于实现。
(4)应用本发明来改进现有的仪器设备,技术上是可行,并且不需要增加任何的硬件开销就可使分析结果可以提高到10-8级。
(5)本方法也同样也适用于进行多次迭代而非一次迭代的谐波分析过程,此时只需要把一次迭代分解成多次迭代实现就可以了。一次迭代和多次迭代本质上是一样的,只是在计算时多次迭代进行分步计算,而一次迭代是把多次迭代的过程合并到迭代系数
Figure 587536DEST_PATH_IMAGE004
中一次计算完成,所以本发明同样适用于多次迭代过程。
附图说明
图1为准同步DFT的谐相角分析误差图。
图2为本发明的谐相角分析误差图。
具体实施方式
(实施例1)
本实施例的一种谐相角分析方法,包括以下步骤:
首先,等间隔采样W+2个采样点,以获得被分析信号的离散序列{f(i),i=0,1,…, Wq+1}W的值由积分方法、迭代次数n和理想周期内采样点数N共同决定。
等间隔采样是指:根据进行谐波分析的理想信号的频率(如工频信号频率f为50Hz,周期为20mS)确定采样频率f S =Nf,在采样频率f S 的作用下在一个周期内均匀地采样N点。一般地,周期采样点N=64或以上就能获得较好的谐波分析结果,而迭代次数n=3~5就能获得较理想的谐波分析结果。
积分方法有复化梯形积分方法、复化矩形积分方法、辛普森方法等多种,可以根据实际情况进行选择。若采用复化梯形积分方法,则W=nN;若采用复化矩形积分方法,则W=n (N-1);若采用复化辛普森积分方法,则W=n(N-1)/2
其次,从采样点i=0开始应用准同步DFT公式
Figure 870750DEST_PATH_IMAGE008
分析W+1个数据获得基波信息
Figure 159124DEST_PATH_IMAGE010
Figure 894999DEST_PATH_IMAGE012
再次,从采样点i=1应用准同步DFT公式分析W+1个数据获得基波信息
Figure 562741DEST_PATH_IMAGE014
Figure 700461DEST_PATH_IMAGE016
再次,应用公式:
Figure 897087DEST_PATH_IMAGE018
,计算信号的频率漂移
Figure 120258DEST_PATH_IMAGE020
再次,从采样点i=0开始应用准同步DFT公式分析W+1个数据获得各次谐波信息
Figure 326112DEST_PATH_IMAGE022
Figure 318338DEST_PATH_IMAGE024
然后,应用公式
Figure DEST_PATH_IMAGE034
计算各次谐波的幅相角;
最后,应用公式
Figure DEST_PATH_IMAGE036
线性修正各次谐波的谐相角。
本技术领域的普通技术人员应当认识到,以上的实施例仅是用来说明本发明,而并非作为对本发明的限定,本发明还可以变化成更多的方式,只要在本发明的实质精神范围内,对以上所述实施例的变化、变型都将落在本发明的权利要求书范围内。

Claims (2)

1.一种谐相角分析方法,其特征在于包括:包括以下步骤:
(1)等间距采样W+2个采样点数据{f(i),i=0,1,…,W+1};所述的采样W+2个采样点数据采用复化梯形积分方法,则W=nN
(2)从采样点i=0开始应用准同步DFT公式
Figure DEST_PATH_IMAGE001
分析W+1个数据获得基波信息
Figure 753338DEST_PATH_IMAGE002
Figure DEST_PATH_IMAGE003
(3)从采样点i=1应用准同步DFT公式分析W+1个数据获得基波信息
Figure 873741DEST_PATH_IMAGE004
Figure DEST_PATH_IMAGE005
(4)应用公式:
Figure 591161DEST_PATH_IMAGE006
,计算信号的频率漂移
Figure DEST_PATH_IMAGE007
(5)从采样点i=0开始应用准同步DFT公式分析W+1个数据获得各次谐波信息
Figure 924053DEST_PATH_IMAGE008
Figure DEST_PATH_IMAGE009
(6)应用公式
Figure 679038DEST_PATH_IMAGE010
计算各次谐波的幅相角;
(7)应用公式
Figure DEST_PATH_IMAGE011
线性修正各次谐波的谐相角;
式中:k为需要获得的谐波的次数;sin和cos分别为正弦和余弦函数;而a k b k 分别为k次谐波的实部和虚部;n为迭代次数;W由积分方法决定;
Figure 439184DEST_PATH_IMAGE012
为一次加权系数;
Figure DEST_PATH_IMAGE013
,为所有加权系数之和;f(i)为分析波形的第i个采样值;N为周期内采样次数;
所述步骤(1)中,等间距采样是根据进行谐波分析的理想信号的周期T和频率f,在一个周期内采样N点,即采样频率为f s =Nf,且N≥64
2.根据权利要求1所述的谐相角分析方法,其特征在于:所述步骤(1)中,所述的采样W+ 2个采样点数据是根据所选择的积分方法而作相应选择,然后根据采样频率f s =Nf,获得采样点数据序列{f(i),i=0,1,…,W+1}n为迭代次数,n≥3;最后对该数据序列进行谐波分析。
CN201811014781.4A 2015-05-19 2015-05-19 谐相角分析方法 Active CN109030942B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811014781.4A CN109030942B (zh) 2015-05-19 2015-05-19 谐相角分析方法

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201510258020.3A CN104849551B (zh) 2015-05-19 2015-05-19 一种谐相角分析方法
CN201811014781.4A CN109030942B (zh) 2015-05-19 2015-05-19 谐相角分析方法

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
CN201510258020.3A Division CN104849551B (zh) 2015-05-19 2015-05-19 一种谐相角分析方法

Publications (2)

Publication Number Publication Date
CN109030942A CN109030942A (zh) 2018-12-18
CN109030942B true CN109030942B (zh) 2020-12-04

Family

ID=53849342

Family Applications (4)

Application Number Title Priority Date Filing Date
CN201811014781.4A Active CN109030942B (zh) 2015-05-19 2015-05-19 谐相角分析方法
CN201510258020.3A Active CN104849551B (zh) 2015-05-19 2015-05-19 一种谐相角分析方法
CN201811015823.6A Active CN108982966B (zh) 2015-05-19 2015-05-19 基于线性修正算法的谐相角分析方法
CN201811014772.5A Active CN109142866B (zh) 2015-05-19 2015-05-19 一种基于线性修正算法的谐相角分析方法

Family Applications After (3)

Application Number Title Priority Date Filing Date
CN201510258020.3A Active CN104849551B (zh) 2015-05-19 2015-05-19 一种谐相角分析方法
CN201811015823.6A Active CN108982966B (zh) 2015-05-19 2015-05-19 基于线性修正算法的谐相角分析方法
CN201811014772.5A Active CN109142866B (zh) 2015-05-19 2015-05-19 一种基于线性修正算法的谐相角分析方法

Country Status (1)

Country Link
CN (4) CN109030942B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109683016A (zh) * 2019-01-18 2019-04-26 江苏理工学院 一种谐波幅值分析方法
CN110083891A (zh) * 2019-04-10 2019-08-02 江苏理工学院 一种准同步dft幅值线性修正系数的计算方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5661761A (en) * 1992-07-09 1997-08-26 Fujitsu Limited Quasi-synchronous detection and demodulation circuit and frequency discriminator used for the same
EP1363438A3 (en) * 1998-01-30 2004-02-11 Matsushita Electric Industrial Co., Ltd. Modulation method and radio communication system
CN101261293A (zh) * 2007-03-08 2008-09-10 国网南京自动化研究院 基于自适应滤波器的电力稳态信号跟踪测量法
CN102135567A (zh) * 2011-02-25 2011-07-27 深圳市业通达实业有限公司 一种电力系统交流采样实时频率跟踪及谐波测量方法
CN102156250A (zh) * 2011-03-17 2011-08-17 华北电力大学(保定) 一种基于等效模型的介质损耗因数测量方法
CN102393488A (zh) * 2011-08-24 2012-03-28 江苏技术师范学院 一种谐波分析方法
CN102435860A (zh) * 2011-10-21 2012-05-02 江苏技术师范学院 一种介质损耗电流测试仪的工作方法
CN103604989A (zh) * 2011-10-21 2014-02-26 蒋春花 一种电能质量谐波分析仪的谐波分析方法
CN103760425A (zh) * 2014-01-22 2014-04-30 湖南大学 一种基于时域准同步的介损角快速测量方法及其测量装置

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS56148181A (en) * 1980-04-21 1981-11-17 Toshiba Corp Method of controlling single phase multistage rectifier
JP2006276006A (ja) * 2005-03-01 2006-10-12 Nagoya Institute Of Technology 電力系統における高調波解析法
CN102331526B (zh) * 2011-06-10 2013-07-17 中国矿业大学 汉宁窗函数连续频谱内插电力谐波参数获取方法
CN103439566B (zh) * 2011-10-21 2014-12-24 常州顺创电气科技有限公司 一种moa阻性电流测试仪的工作方法
CN103969507B (zh) * 2011-10-21 2016-10-19 江苏理工学院 一种电能质量谐波分析方法
CN102508026B (zh) * 2011-10-21 2014-08-06 江苏理工学院 一种电能质量谐波分析仪的谐波分析方法
CN103592512A (zh) * 2011-10-21 2014-02-19 蒋春花 一种电能质量谐波分析仪的谐波分析方法
JP5705102B2 (ja) * 2011-12-21 2015-04-22 三菱電機株式会社 絶縁劣化診断装置
CN102590615A (zh) * 2012-02-16 2012-07-18 安徽理工大学 电网的电压和电流信号的频率偏移情况下同步相量修正的方法
DE102012204963A1 (de) * 2012-03-28 2013-10-02 Siemens Aktiengesellschaft Verfahren und Vorrichtung zur Reduzierung von Strom-Oberschwingungen
CN103576002B (zh) * 2013-11-11 2016-01-20 华北电力大学(保定) 一种容性绝缘设备介质损耗角的计算方法
CN104181391A (zh) * 2014-08-21 2014-12-03 艾德克斯电子(南京)有限公司 数字功率计谐波检测的方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5661761A (en) * 1992-07-09 1997-08-26 Fujitsu Limited Quasi-synchronous detection and demodulation circuit and frequency discriminator used for the same
EP1363438A3 (en) * 1998-01-30 2004-02-11 Matsushita Electric Industrial Co., Ltd. Modulation method and radio communication system
CN101261293A (zh) * 2007-03-08 2008-09-10 国网南京自动化研究院 基于自适应滤波器的电力稳态信号跟踪测量法
CN102135567A (zh) * 2011-02-25 2011-07-27 深圳市业通达实业有限公司 一种电力系统交流采样实时频率跟踪及谐波测量方法
CN102156250A (zh) * 2011-03-17 2011-08-17 华北电力大学(保定) 一种基于等效模型的介质损耗因数测量方法
CN102393488A (zh) * 2011-08-24 2012-03-28 江苏技术师范学院 一种谐波分析方法
CN103454494A (zh) * 2011-08-24 2013-12-18 常州顺创电气科技有限公司 一种高精度的谐波分析方法
CN102435860A (zh) * 2011-10-21 2012-05-02 江苏技术师范学院 一种介质损耗电流测试仪的工作方法
CN103604989A (zh) * 2011-10-21 2014-02-26 蒋春花 一种电能质量谐波分析仪的谐波分析方法
CN103760425A (zh) * 2014-01-22 2014-04-30 湖南大学 一种基于时域准同步的介损角快速测量方法及其测量装置

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
High Precision Phase Measurement Using Reduced Sine and Cosine Tables;SYED MASUD MAHMUD 等;《IEEE TRANSACTIONS ON INSTRUMENTATION AND MEASUREMENT》;19900228;第39卷(第1期);第56-60页 *
QUASI-SYNCHRONOUS SAMPLING ALGORITHM AND ITS APPLICATIONS——3. HIGH ACCURATE MEASUREMENT OF FREQUENCY,FREQUENCY DEVIATION AND PHASE ANGLE DIFFERENCE IN POWER SYSTEMS;Xianzhong Dai 等;《Proceedings of the 1993 IEEE Instrumentation and Measurement Technology Conference》;19931231;第726-729页 *
基于准同步DFT的非整数谐波分析算法;傅中君 等;《仪器仪表学报》;20120131;第33卷(第1期);第235-240页 *
频谱校正理论的发展;丁康 等;《振动工程学报》;20000331;第13卷(第1期);第14-22页 *

Also Published As

Publication number Publication date
CN109142866A (zh) 2019-01-04
CN108982966B (zh) 2020-10-09
CN108982966A (zh) 2018-12-11
CN109030942A (zh) 2018-12-18
CN104849551A (zh) 2015-08-19
CN109142866B (zh) 2020-10-09
CN104849551B (zh) 2018-11-23

Similar Documents

Publication Publication Date Title
CN103454494B (zh) 一种高精度的谐波分析方法
Belega et al. Frequency estimation by two-or three-point interpolated Fourier algorithms based on cosine windows
CN102508026B (zh) 一种电能质量谐波分析仪的谐波分析方法
CN110837001B (zh) 一种电力系统中谐波和间谐波的分析方法与装置
CN104897960B (zh) 基于加窗四谱线插值fft的谐波快速分析方法及系统
CN103869162A (zh) 一种基于时域准同步的动态信号相量测量方法
CN105137180B (zh) 基于六项余弦窗四谱线插值的高精度谐波分析方法
CN102435860B (zh) 一种介质损耗电流测试仪的工作方法
CN103983849B (zh) 一种实时高精度的电力谐波分析方法
CN109239463B (zh) 一种基于线性修正算法的介质损耗测量方法
CN109030942B (zh) 谐相角分析方法
CN108896944A (zh) 一种同步测量装置实验室校准仪及其同步相量测量方法
CN109490603B (zh) 阻性电流基波的测量方法
CN104931777B (zh) 一种基于两条dft复数谱线的信号频率测量方法
CN104407197B (zh) 一种基于三角函数迭代的信号相量测量的方法
CN103543331A (zh) 一种计算电信号谐波和间谐波的方法
CN105372492A (zh) 基于三条dft复数谱线的信号频率测量方法
CN104678170A (zh) 一种基于谐波分析仪的电力谐波分析方法和谐波分析仪
CN103983852B (zh) 电能质量谐波分析仪的谐波分析方法
CN110083891A (zh) 一种准同步dft幅值线性修正系数的计算方法
CN103592512A (zh) 一种电能质量谐波分析仪的谐波分析方法
CN103604989A (zh) 一种电能质量谐波分析仪的谐波分析方法
CN105372493B (zh) 基于三条dft复数谱线的信号幅值和相位测量方法
CN103969507B (zh) 一种电能质量谐波分析方法
CN109490625B (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