CN108982966B - 基于线性修正算法的谐相角分析方法 - Google Patents
基于线性修正算法的谐相角分析方法 Download PDFInfo
- Publication number
- CN108982966B CN108982966B CN201811015823.6A CN201811015823A CN108982966B CN 108982966 B CN108982966 B CN 108982966B CN 201811015823 A CN201811015823 A CN 201811015823A CN 108982966 B CN108982966 B CN 108982966B
- 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
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 55
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 10
- 238000005070 sampling Methods 0.000 claims abstract description 49
- 238000000034 method Methods 0.000 claims abstract description 37
- 230000010354 integration Effects 0.000 claims abstract description 23
- 238000012544 monitoring process Methods 0.000 abstract description 8
- 238000007689 inspection Methods 0.000 abstract description 4
- 238000004519 manufacturing process Methods 0.000 abstract description 3
- 230000009286 beneficial effect Effects 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000004888 barrier function Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000012804 iterative process 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
- 230000008092 positive effect Effects 0.000 description 1
Images
Landscapes
- Measuring Frequencies, Analyzing Spectra (AREA)
- Measuring Phase Differences (AREA)
- Complex Calculations (AREA)
Abstract
Description
本申请是申请号为:201510258020.3,发明创造名称为《一种谐相角分析方法》,申请日为:2015年05月19日的发明专利申请的分案申请。
技术领域
本发明涉及一种高精度的谐相角分析方法。
背景技术
谐波分析技术在电能质量监控、电子产品生产检验、电器设备监控等众多领域应用广泛,是进行电网监控、质量检验、设备监控的重要技术手段。目前谐波分析应用最广泛的技术是离散傅里叶变换(DFT)和快速傅里叶变换(FFT)。准同步采样技术和DFT技术相结合的谐波分析技术能够提高谐波分析的精度,其算式为:
式中:k为需要获得的谐波的次数(如基波k=1,3次谐波k=3);sin和cos分别为正弦和余弦函数;而a k 和b k 分别为k次谐波的实部和虚部;n为迭代次数;W由积分方法决定,采用复化梯形积分方法时,W=nN;为一次加权系数;,为所有加权系数之和;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公式
等间隔采样是根据进行谐波分析的理想信号的周期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;最后对该数据序列进行谐波分析。
本发明具有积极的效果:(1)高精度的谐相角分析结果。如对于图1给定的分析实例,本发明获得的分析精度提高到10-8级(图2)。
(2)本发明所述的方法从根本上解决了准同步DFT谐相角分析精度低的问题,而无需进行复杂的反演和修正,算法简单。
(3)相对于准同步DFT,本发明所述的谐波分析技术只需要增加一个采样点就解决了准同步DFT分析误差大的问题,易于实现。
(4)应用本发明来改进现有的仪器设备,技术上是可行,并且不需要增加任何的硬件开销就可使分析结果可以提高到10-8级。
(5)本方法也同样也适用于进行多次迭代而非一次迭代的谐波分析过程,此时只需要把一次迭代分解成多次迭代实现就可以了。一次迭代和多次迭代本质上是一样的,只是在计算时多次迭代进行分步计算,而一次迭代是把多次迭代的过程合并到迭代系数中一次计算完成,所以本发明同样适用于多次迭代过程。
附图说明
图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公式
本技术领域的普通技术人员应当认识到,以上的实施例仅是用来说明本发明,而并非作为对本发明的限定,本发明还可以变化成更多的方式,只要在本发明的实质精神范围内,对以上所述实施例的变化、变型都将落在本发明的权利要求书范围内。
Claims (2)
1.一种基于线性修正算法的谐相角分析方法,其特征在于:包括以下步骤:
(1)等间距采样W+2个采样点数据{f(i),i=0,1,…,W+1};所述的采样W+2个采样点数据采用复化辛普森积分方法,则W=n(N-1)/2;
(2)从采样点i=0开始应用准同步DFT公式
式中:k为需要获得的谐波的次数;sin和cos分别为正弦和余弦函数;而a k 和b k 分别为k次谐波的实部和虚部;n为迭代次数;W由积分方法决定;为一次加权系数;为所有加权系数之和;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;最后对该数据序列进行谐波分析。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811015823.6A CN108982966B (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 | 一种谐相角分析方法 |
CN201811015823.6A CN108982966B (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 |
---|---|
CN108982966A CN108982966A (zh) | 2018-12-11 |
CN108982966B true CN108982966B (zh) | 2020-10-09 |
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 After (3)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
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 | 谐相角分析方法 |
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幅值线性修正系数的计算方法 |
Family Cites Families (21)
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 | 中国矿业大学 | 汉宁窗函数连续频谱内插电力谐波参数获取方法 |
CN102393488B (zh) * | 2011-08-24 | 2013-12-11 | 江苏理工学院 | 一种谐波分析方法 |
CN103439566B (zh) * | 2011-10-21 | 2014-12-24 | 常州顺创电气科技有限公司 | 一种moa阻性电流测试仪的工作方法 |
CN103969507B (zh) * | 2011-10-21 | 2016-10-19 | 江苏理工学院 | 一种电能质量谐波分析方法 |
CN102435860B (zh) * | 2011-10-21 | 2013-11-06 | 江苏理工学院 | 一种介质损耗电流测试仪的工作方法 |
CN103592512A (zh) * | 2011-10-21 | 2014-02-19 | 蒋春花 | 一种电能质量谐波分析仪的谐波分析方法 |
CN102508026B (zh) * | 2011-10-21 | 2014-08-06 | 江苏理工学院 | 一种电能质量谐波分析仪的谐波分析方法 |
CN103604989A (zh) * | 2011-10-21 | 2014-02-26 | 蒋春花 | 一种电能质量谐波分析仪的谐波分析方法 |
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 | 华北电力大学(保定) | 一种容性绝缘设备介质损耗角的计算方法 |
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
Also Published As
Publication number | Publication date |
---|---|
CN104849551A (zh) | 2015-08-19 |
CN104849551B (zh) | 2018-11-23 |
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 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103454494B (zh) | 一种高精度的谐波分析方法 | |
CN102508026B (zh) | 一种电能质量谐波分析仪的谐波分析方法 | |
CN110837001B (zh) | 一种电力***中谐波和间谐波的分析方法与装置 | |
CN104897960B (zh) | 基于加窗四谱线插值fft的谐波快速分析方法及*** | |
CN103869162A (zh) | 一种基于时域准同步的动态信号相量测量方法 | |
CN106018956B (zh) | 一种加窗谱线插值的电力***频率计算方法 | |
CN102435860B (zh) | 一种介质损耗电流测试仪的工作方法 | |
CN105137180B (zh) | 基于六项余弦窗四谱线插值的高精度谐波分析方法 | |
CN103018555B (zh) | 一种高精度的电力参数软件同步采样方法 | |
CN105307095B (zh) | 一种基于fft的高分辨率音频频率测量方法 | |
CN103983849B (zh) | 一种实时高精度的电力谐波分析方法 | |
CN108896944A (zh) | 一种同步测量装置实验室校准仪及其同步相量测量方法 | |
CN109239463B (zh) | 一种基于线性修正算法的介质损耗测量方法 | |
CN108982966B (zh) | 基于线性修正算法的谐相角分析方法 | |
CN108957097B (zh) | 金属氧化物避雷器的阻性电流基波的测量方法 | |
CN104931777B (zh) | 一种基于两条dft复数谱线的信号频率测量方法 | |
CN105372492B (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 |