CN104849551B - 一种谐相角分析方法 - Google Patents
一种谐相角分析方法 Download PDFInfo
- Publication number
- CN104849551B CN104849551B CN201510258020.3A CN201510258020A CN104849551B CN 104849551 B CN104849551 B CN 104849551B CN 201510258020 A CN201510258020 A CN 201510258020A CN 104849551 B CN104849551 B CN 104849551B
- Authority
- CN
- China
- Prior art keywords
- formula
- analysis
- data
- phase angle
- harmonic
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 42
- 238000000034 method Methods 0.000 claims abstract description 37
- 230000010354 integration Effects 0.000 claims abstract description 24
- 238000005070 sampling Methods 0.000 claims abstract description 19
- 238000012544 monitoring process Methods 0.000 abstract description 7
- 238000004519 manufacturing process Methods 0.000 abstract description 3
- 238000007689 inspection Methods 0.000 abstract description 2
- 230000009286 beneficial effect Effects 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 2
- 241001515806 Stictis Species 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
Landscapes
- Measuring Frequencies, Analyzing Spectra (AREA)
- Measuring Phase Differences (AREA)
- Complex Calculations (AREA)
Abstract
本发明涉及一种谐相角分析方法,包括如下几个步骤:(1)采样W+2个采样点数据(W由积分方法决定);(2)从采样点i=0应用准同步DFT公式分析W+1个数据获得基波信息和(3)从采样点i=1应用准同步DFT公式分析W+1个数据获得基波信息和(4)应用公式计算信号的频率漂移μ;(5)从采样点i=O应用准同步DFT公式分析W+1个数据获得各次谐波信息和(6)应用公式汁算各次谐波的幅相角;(7)应用公式线性修正各次谐波的谐相角。本发明的方法,有助于电能质量监控、电子产品生产检验、电器设备监控等应用谐波分析的领域更加精确的获得各次谐波的幅相角和频率等信息。
Description
技术领域
本发明涉及一种高精度的谐相角分析方法。
背景技术
谐波分析技术在电能质量监控、电子产品生产检验、电器设备监控等众多领域应用广泛,是进行电网监控、质量检验、设备监控的重要技术手段。目前谐波分析应用最广泛的技术是离散傅里叶变换(DFT)和快速傅里叶变换(FFT)。准同步采样技术和DFT技术相结合的谐波分析技术能够提高谐波分析的精度,其算式为:
式中:k为需要获得的谐波的次数(如基波k=1,3次谐波k=3);sin和cos分别为正弦和余弦函数;而ak和bk分别为k次谐波的实部和虚部;n为迭代次数;W由积分方法决定,采用复化梯形积分方法时,W=nN;γi为一次加权系数;为所有加权系数之和;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公式
分析W+1个数据获得基波信息和
(3)从采样点i=1应用准同步DFT公式分析W+1个数据获得基波信息和
(4)应用公式:计算信号的频率漂移μ;
(5)从采样点i=0开始应用准同步DFT公式分析W+1个数据获得各次谐波信息和
(6)应用公式计算各次谐波的幅相角;
(7)应用公式线性修正各次谐波的谐相角。
等间隔采样是根据进行谐波分析的理想信号的周期T和频率f,在一个周期内采样N点,即采样频率为fs=Nf,且N≥64。
所述的采样W+2个采样点数据是根据所选择的积分方法而作相应选择,若采用复化梯形积分方法,则W=nN;若采用复化矩形积分方法,则W=n(N-1);若采用复化辛普森积分方法,则W=n(N-1)/2。然后根据采样频率fs=Nf,获得采样点数据序列{f(i),i=0,1,...,W+1},n为迭代次数,一般n≥3;最后对该数据序列进行谐波分析。
一次迭代系数γi由积分方法、理想周期采样点N和迭代次数n决定,具体推导过程参见文献【戴先中.准同步采样应用中的若干问题[J].电测与仪表,1988,(2):2-7.】。
为所有加权系数之和。
信号频率的漂移μ是根据相邻采样点基波相角差与理想周期内采样点数N的固定关系而获得的,信号频率的漂移μ也可用于修正基波和高次谐波的频率f1和高次谐波的频率fk(fk=kμfs/N)。
本发明具有积极的效果:(1)高精度的谐相角分析结果。如对于图1给定的分析实例,本发明获得的分析精度提高到10-8级(图2)。
(2)本发明所述的方法从根本上解决了准同步DFT谐相角分析精度低的问题,而无需进行复杂的反演和修正,算法简单。
(3)相对于准同步DFT,本发明所述的谐波分析技术只需要增加一个采样点就解决了准同步DFT分析误差大的问题,易于实现。
(4)应用本发明来改进现有的仪器设备,技术上是可行,并且不需要增加任何的硬件开销就可使分析结果可以提高到10-8级。
(5)本方法也同样也适用于进行多次迭代而非一次迭代的谐波分析过程,此时只需要把一次迭代分解成多次迭代实现就可以了。一次迭代和多次迭代本质上是一样的,只是在计算时多次迭代进行分步计算,而一次迭代是把多次迭代的过程合并到迭代系数γi中一次计算完成,所以本发明同样适用于多次迭代过程。
附图说明
图1为准同步DFT的谐相角分析误差图。
图2为本发明的谐相角分析误差图。
具体实施方式
(实施例1)
本实施例的一种谐相角分析方法,包括以下步骤:
首先,等间隔采样W+2个采样点,以获得被分析信号的离散序列{f(i),i=0,1,...,Wq+1}。W的值由积分方法、迭代次数n和理想周期内采样点数N共同决定。
等间隔采样是指:根据进行谐波分析的理想信号的频率(如工频信号频率f为50Hz,周期为20mS)确定采样频率fS=Nf,在采样频率fS的作用下在一个周期内均匀地采样N点。一般地,周期采样点N=64或以上就能获得较好的谐波分析结果,而迭代次数n=3~5就能获得较理想的谐波分析结果。
积分方法有复化梯形积分方法、复化矩形积分方法、辛普森方法等多种,可以根据实际情况进行选择。若采用复化梯形积分方法,则W=nN;若采用复化矩形积分方法,则W=n(N-1);若采用复化辛普森积分方法,则W=n(N-1)/2。
其次,从采样点i=0开始应用准同步DFT公式
分析W+1个数据获得基波信息和
再次,从采样点i=1应用准同步DFT公式分析W+1个数据获得基波信息和
再次,应用公式:计算信号的频率漂移μ;
再次,从采样点i=0开始应用准同步DFT公式分析W+1个数据获得各次谐波信息和
然后,应用公式计算各次谐波的幅相角;
最后,应用公式线性修正各次谐波的谐相角。
本技术领域的普通技术人员应当认识到,以上的实施例仅是用来说明本发明,而并非作为对本发明的限定,本发明还可以变化成更多的方式,只要在本发明的实质精神范围内,对以上所述实施例的变化、变型都将落在本发明的权利要求书范围内。
Claims (1)
1.一种谐相角分析方法,其特征在于包括:包括以下步骤:
(1)等间距采样W+2个采样点数据{f(i),i=0,1,...,W+1};
(2)从采样点i=0开始应用准同步DFT公式
分析W+1个数据获得基波信息和
(3)从采样点i=1应用准同步DFT公式分析W+1个数据获得基波信息和
(4)应用公式:计算信号的频率漂移μ;
(5)从采样点i=0开始应用准同步DFT公式分析W+1个数据获得各次谐波信息和
(6)应用公式计算各次谐波的幅相角;
(7)应用公式线性修正各次谐波的谐相角;
式中:k为需要获得的谐波的次数;sin和cos分别为正弦和余弦函数;而ak和bk分别为k次谐波的实部和虚部;n为迭代次数;W由积分方法决定,采用复化梯形积分方法时,W=nN;γi为一次加权系数;为所有加权系数之和;f(i)为分析波形的第i个采样值;N为周期内采样次数。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811014772.5A CN109142866B (zh) | 2015-05-19 | 2015-05-19 | 一种基于线性修正算法的谐相角分析方法 |
CN201510258020.3A CN104849551B (zh) | 2015-05-19 | 2015-05-19 | 一种谐相角分析方法 |
CN201811015823.6A CN108982966B (zh) | 2015-05-19 | 2015-05-19 | 基于线性修正算法的谐相角分析方法 |
CN201811014781.4A CN109030942B (zh) | 2015-05-19 | 2015-05-19 | 谐相角分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510258020.3A CN104849551B (zh) | 2015-05-19 | 2015-05-19 | 一种谐相角分析方法 |
Related Child Applications (3)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811014772.5A Division CN109142866B (zh) | 2015-05-19 | 2015-05-19 | 一种基于线性修正算法的谐相角分析方法 |
CN201811014781.4A Division CN109030942B (zh) | 2015-05-19 | 2015-05-19 | 谐相角分析方法 |
CN201811015823.6A Division CN108982966B (zh) | 2015-05-19 | 2015-05-19 | 基于线性修正算法的谐相角分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104849551A CN104849551A (zh) | 2015-08-19 |
CN104849551B true CN104849551B (zh) | 2018-11-23 |
Family
ID=53849342
Family Applications (4)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811015823.6A Active CN108982966B (zh) | 2015-05-19 | 2015-05-19 | 基于线性修正算法的谐相角分析方法 |
CN201510258020.3A Active CN104849551B (zh) | 2015-05-19 | 2015-05-19 | 一种谐相角分析方法 |
CN201811014772.5A Active CN109142866B (zh) | 2015-05-19 | 2015-05-19 | 一种基于线性修正算法的谐相角分析方法 |
CN201811014781.4A Active CN109030942B (zh) | 2015-05-19 | 2015-05-19 | 谐相角分析方法 |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811015823.6A Active CN108982966B (zh) | 2015-05-19 | 2015-05-19 | 基于线性修正算法的谐相角分析方法 |
Family Applications After (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811014772.5A Active CN109142866B (zh) | 2015-05-19 | 2015-05-19 | 一种基于线性修正算法的谐相角分析方法 |
CN201811014781.4A Active CN109030942B (zh) | 2015-05-19 | 2015-05-19 | 谐相角分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (4) | CN108982966B (zh) |
Families Citing this family (2)
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 (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102393488A (zh) * | 2011-08-24 | 2012-03-28 | 江苏技术师范学院 | 一种谐波分析方法 |
CN102508026A (zh) * | 2011-10-21 | 2012-06-20 | 江苏技术师范学院 | 一种电能质量谐波分析仪的谐波分析方法 |
DE102012204963A1 (de) * | 2012-03-28 | 2013-10-02 | Siemens Aktiengesellschaft | Verfahren und Vorrichtung zur Reduzierung von Strom-Oberschwingungen |
CN103439566A (zh) * | 2011-10-21 | 2013-12-11 | 常州顺创电气科技有限公司 | 精度较高的moa阻性电流测试仪的工作方法 |
CN103604989A (zh) * | 2011-10-21 | 2014-02-26 | 蒋春花 | 一种电能质量谐波分析仪的谐波分析方法 |
Family Cites Families (16)
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 |
US5400366A (en) * | 1992-07-09 | 1995-03-21 | Fujitsu Limited | Quasi-synchronous detection and demodulation circuit and frequency discriminator used for the same |
EP1628452A3 (en) * | 1998-01-30 | 2006-05-10 | Matsushita Electric Industrial Co., Ltd. | Modulation method and radio communication system |
JP2006276006A (ja) * | 2005-03-01 | 2006-10-12 | Nagoya Institute Of Technology | 電力系統における高調波解析法 |
CN101261293A (zh) * | 2007-03-08 | 2008-09-10 | 国网南京自动化研究院 | 基于自适应滤波器的电力稳态信号跟踪测量法 |
CN102135567A (zh) * | 2011-02-25 | 2011-07-27 | 深圳市业通达实业有限公司 | 一种电力***交流采样实时频率跟踪及谐波测量方法 |
CN102156250A (zh) * | 2011-03-17 | 2011-08-17 | 华北电力大学(保定) | 一种基于等效模型的介质损耗因数测量方法 |
CN102331526B (zh) * | 2011-06-10 | 2013-07-17 | 中国矿业大学 | 汉宁窗函数连续频谱内插电力谐波参数获取方法 |
CN103969507B (zh) * | 2011-10-21 | 2016-10-19 | 江苏理工学院 | 一种电能质量谐波分析方法 |
CN102435860B (zh) * | 2011-10-21 | 2013-11-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 | 安徽理工大学 | 电网的电压和电流信号的频率偏移情况下同步相量修正的方法 |
CN103576002B (zh) * | 2013-11-11 | 2016-01-20 | 华北电力大学(保定) | 一种容性绝缘设备介质损耗角的计算方法 |
CN103760425A (zh) * | 2014-01-22 | 2014-04-30 | 湖南大学 | 一种基于时域准同步的介损角快速测量方法及其测量装置 |
CN104181391A (zh) * | 2014-08-21 | 2014-12-03 | 艾德克斯电子(南京)有限公司 | 数字功率计谐波检测的方法 |
-
2015
- 2015-05-19 CN CN201811015823.6A patent/CN108982966B/zh active Active
- 2015-05-19 CN CN201510258020.3A patent/CN104849551B/zh active Active
- 2015-05-19 CN CN201811014772.5A patent/CN109142866B/zh active Active
- 2015-05-19 CN CN201811014781.4A patent/CN109030942B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102393488A (zh) * | 2011-08-24 | 2012-03-28 | 江苏技术师范学院 | 一种谐波分析方法 |
CN102508026A (zh) * | 2011-10-21 | 2012-06-20 | 江苏技术师范学院 | 一种电能质量谐波分析仪的谐波分析方法 |
CN103439566A (zh) * | 2011-10-21 | 2013-12-11 | 常州顺创电气科技有限公司 | 精度较高的moa阻性电流测试仪的工作方法 |
CN103604989A (zh) * | 2011-10-21 | 2014-02-26 | 蒋春花 | 一种电能质量谐波分析仪的谐波分析方法 |
DE102012204963A1 (de) * | 2012-03-28 | 2013-10-02 | Siemens Aktiengesellschaft | Verfahren und Vorrichtung zur Reduzierung von Strom-Oberschwingungen |
Also Published As
Publication number | Publication date |
---|---|
CN104849551A (zh) | 2015-08-19 |
CN109030942A (zh) | 2018-12-18 |
CN109142866A (zh) | 2019-01-04 |
CN109142866B (zh) | 2020-10-09 |
CN108982966A (zh) | 2018-12-11 |
CN109030942B (zh) | 2020-12-04 |
CN108982966B (zh) | 2020-10-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103454497B (zh) | 基于改进加窗离散傅立叶变换的相位差测量方法 | |
CN103869162B (zh) | 一种基于时域准同步的动态信号相量测量方法 | |
CN103454494B (zh) | 一种高精度的谐波分析方法 | |
CN109752584B (zh) | 一种周期信号有效值测量方法 | |
CN102508026B (zh) | 一种电能质量谐波分析仪的谐波分析方法 | |
CN103983849B (zh) | 一种实时高精度的电力谐波分析方法 | |
CN102435860B (zh) | 一种介质损耗电流测试仪的工作方法 | |
CN102809687A (zh) | 一种交流电频率的数字化测量方法 | |
CN104849569B (zh) | 一种介质损耗测量方法 | |
CN103969508A (zh) | 一种实时高精密的电力谐波分析方法及装置 | |
CN104849551B (zh) | 一种谐相角分析方法 | |
CN103543331B (zh) | 一种计算电信号谐波和间谐波的方法 | |
CN104849530B (zh) | 一种moa阻性电流基波的测量方法 | |
CN105067882B (zh) | 一种电量幅值的测定方法 | |
Salcic et al. | An improved Taylor method for frequency measurement in power systems | |
CN104407197A (zh) | 一种基于三角函数迭代的信号相量测量的方法 | |
CN103983852B (zh) | 电能质量谐波分析仪的谐波分析方法 | |
CN105277780A (zh) | 电压过零点的获取方法和装置 | |
CN110083891A (zh) | 一种准同步dft幅值线性修正系数的计算方法 | |
CN103592513A (zh) | 电力信号谐波分析方法和装置 | |
CN104764930B (zh) | 基于相关特性的谐波发射水平估计方法 | |
CN103969507B (zh) | 一种电能质量谐波分析方法 | |
CN103604989A (zh) | 一种电能质量谐波分析仪的谐波分析方法 | |
CN103592512A (zh) | 一种电能质量谐波分析仪的谐波分析方法 | |
CN103267896B (zh) | 周期信号初相角的测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CP02 | Change in the address of a patent holder | ||
CP02 | Change in the address of a patent holder |
Address after: No. 1801 Zhong Wu Avenue, Changzhou, Jiangsu Province, Jiangsu Patentee after: Jiangsu University of Technology Address before: 213001 1801 Zhong Wu Avenue, Zhong Lou District, Changzhou, Jiangsu Patentee before: Jiangsu University of Technology |