CN107632200B - 基于离散傅里叶变换的谐波检测方法 - Google Patents
基于离散傅里叶变换的谐波检测方法 Download PDFInfo
- Publication number
- CN107632200B CN107632200B CN201710773955.4A CN201710773955A CN107632200B CN 107632200 B CN107632200 B CN 107632200B CN 201710773955 A CN201710773955 A CN 201710773955A CN 107632200 B CN107632200 B CN 107632200B
- Authority
- CN
- China
- Prior art keywords
- current
- phase
- amplitude
- active
- fundamental wave
- 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
Images
Landscapes
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于离散傅里叶变换的谐波检测方法,其包括基于离散傅里叶变换的谐波检测方法,其特征在于,包括在设定周期内采集下,采集电网侧的三相电流,并对采集的三相电流做离散处理;构建三相电流的基波有功电流幅值和基波无功电流幅值的计算模型;通过滑窗迭代方法对采集电网侧的三相电流进行离散傅里叶变换正变换,分别得到三相电流的基波有功电流幅值和基波无功电流幅值;采用坐标变换计算三相电流的基波有功电流正序分量;根据采集的三相电流及基波有功电流正序分量,计算检测的三相谐波电流。
Description
技术领域
本发明涉及数字信号处理技术领域,具体涉及一种基于离散傅里叶变换的谐波检测方法。
背景技术
在并联型有源电力滤波***(APF)的设计中,需要对电网侧电流进行补偿,实现功率因数的校正,通常的做法是检测网侧电流基波分量,再将网侧电流与基波电流相减,得到谐波电流量作为电流指令信号产生PWM开关信号。在实际的工程研制中谐波检测技术是解决问题的关键。
目前传统的谐波检测方法(ip-iq法)是建立在瞬时无功理论基础上,采样三相网侧电流ia、ib、ic,通过park变换得到有功电流分量ip和无功电流分量iq,设计两个数字低通滤波环节(LPF),得到其直流分量。再通过2r-3s反变换得到网侧电流基波分量iaf、ibf、icf,最后与采样值做差得到三相谐波电流。
该方法控制原理简单明了,但谐波检测效果与LPF性能有直接关系。LPF截止频率越低,谐波检测精度越高,但动态响应特性过慢;LPF截止频率过高,对低次谐波检测效果变差,又会影响谐波检测精度,且该方法还无法实现单次谐波的检测。
发明内容
针对现有技术中的上述不足,本发明提供了一种快速、准确检测谐波的基于离散傅里叶变换的谐波检测方法。
为了达到上述发明目的,本发明采用的技术方案为:
提供一种基于离散傅里叶变换的谐波检测方法,其包括:
在设定周期内采集下,采集电网侧的三相电流,并对采集的三相电流做离散处理;三相电流分别为a相采集电流、b相采集电流和c相采集电流;
构建三相电流的基波有功电流幅值和基波无功电流幅值的计算模型:
其中,XA为采样电流的基波有功电流幅值;XB为采样电流的基波无功电流幅值;n为采样点的序号;N为一个采样周期内的采样点个数;i为采样电流;
通过滑窗迭代方法对采集电网侧的三相电流进行离散傅里叶变换正变换,分别得到三相电流的基波有功电流幅值和基波无功电流幅值:
当采样点的个数小于采样周期内的采样点个数时,三相电流的基波有功电流幅值和基波无功电流幅值的迭代公式:
当采样点的个数大于采样周期内的采样点个数时,三相电流的基波有功电流幅值和基波无功电流幅值的迭代公式:
采用坐标变换计算三相电流的基波有功电流正序分量:
其中,iA为a相基波有功电流正序分量;iB为b相基波有功电流正序分量;iC为c相基波有功电流正序分量;XaA为a相电流的有功幅值;XaB为a相电流的无功幅值;XbA为b相电流的有功幅值;XbB为b相电流的无功幅值;XcA为c相电流的有功幅值;XcB为c相电流的无功幅值;
根据采集的三相电流及基波有功电流正序分量,计算检测的三相谐波电流:
其中,iah为a相的谐波电流;ibh为b相的谐波电流;ich为c相的谐波电流;ia为a相采集电流;ib为b相采集电流;ic为c相采集电流。
本发明的有益效果为:本方案在负载发生变化后的第一个周期内和一个周期后采用不同的方式进行三相电流基波无功电流幅值,通过这种方式再结合基波有功电流正序分量算例求解出的谐波电流,能够在一个采样周期内将电网侧的三相电流调整至稳定状态。
本方案在进行谐波电流计算过程中每一步迭代的乘法因子都是常数,同时无需重构锁相环,因此运算量大大减少,从而大幅度提高了***的相应速率,该方法适用于并联型有源电力滤波***的谐波电流检测时,可以实现快速响应,并且改善THD(谐波失真)指标,从而实现指定次数谐波的补偿。
采用本方法能实现快速准确地完成谐波电流的检测,计算量小,抗干扰能力强,易于工程应用实现,其能够保证检测速度的同时保证检测精度。
附图说明
图1为基于离散傅里叶变换的谐波检测方法的流程图。
图2为滑窗迭代DFT算法原理图。
图3为现有技术的方法与本方案的方法进行谐波检测时的响应能力对比图。
图4为三相四线的APF的***模型的a相仿真波形图。
图5为三相四线的APF的***模型补偿后的三相电流波形图。
图6a为补偿前的谐波失真的仿真图。
图6b为补偿后的谐波失真的仿真图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
参考图1,图1示出了基于离散傅里叶变换的谐波检测方法的流程图;如图1所示,该谐波检测方法100包括步骤101至步骤105。
在步骤101中,在设定周期内采集下,采集三相四线制并联型有源电力滤波***电网侧的三相电流,并采用AD芯片对对采集的三相电流做做离散处理;其中,三相电流分别为a相采集电流、b相采集电流和c相采集电流。
在步骤102中,构建三相电流的基波有功电流幅值和基波无功电流幅值的计算模型:
其中,XA为采样电流的基波有功电流幅值;XB为采样电流的基波无功电流幅值;n为采样点的序号;N为一个采样周期内的采样点个数;i为采样电流;
在步骤103中,通过滑窗迭代方法(参考图2)对采集电网侧的三相电流进行离散傅里叶变换正变换,分别得到三相电流的基波有功电流幅值和基波无功电流幅值:
当采样点的个数小于采样周期内的采样点个数时,三相电流的基波有功电流幅值和基波无功电流幅值的迭代公式:
当采样点的个数大于采样周期内的采样点个数时,三相电流的基波有功电流幅值和基波无功电流幅值的迭代公式:
在步骤104中,采用坐标变换计算三相电流的基波有功电流正序分量:
其中,iA为a相基波有功电流正序分量;iB为b相基波有功电流正序分量;iC为c相基波有功电流正序分量;XaA为a相电流的有功幅值;XaB为a相电流的无功幅值;XbA为b相电流的有功幅值;XbB为b相电流的无功幅值;XcA为c相电流的有功幅值;XcB为c相电流的无功幅值。
在本发明的一个实施例中,采用坐标变换计算三相电流的基波有功电流正序分量进一步包括:
计算三相电流的正序分量:
提取三相电流正序分量的实部,舍掉虚部得到三相电流的基波有功电流正序分量。
在步骤105中,根据采集的三相电流及基波有功电流正序分量,计算检测的三相谐波电流:
其中,iah为a相的谐波电流;ibh为b相的谐波电流;ich为c相的谐波电流;ia为a相采集电流;ib为b相采集电流;ic为c相采集电流。
下面结合具体的实例,对本方案的谐波检测方法的效果进行仿真:
设定采样频率6.4K,待采样的电流信号包含5次、7次谐波,并且幅值与相位偏移均不同,采用Matlab Simulink仿真平台采用现有技术的ip-iq法和本方案的谐波检测方法进行仿真,两种谐波检测方式的响应能力如图3所示,由图3可知,本方案在0s和0.2s出现负载变换时,本方案在一个周期后电流就趋于稳定,而现有技术需要3个周期的时间才能使电流趋于稳定,由此可见采用本方案进行谐波检测能够大幅度提高相应速度。
随后搭建三相四线的APF的***模型,仿真参数设置如下:
电网侧220V/50Hz三相交流电;负载侧经过不可控整流桥接阻感负载,R=10Ω,L=2mH;逆变电路输出滤波电感La=Lb=Lc=1mH;直流侧储能电容C1=C2=20000uH,采样频率设置为24kHz,一个周期有480个采样点,仿真时间设置为1s。***稳定运行,并且保持了较好的补偿电流跟踪特性。使用本方案的谐波检测方法实现了***的控制,仿真结构如图4至图6b。
通过仿真波形可以看出,本方案的谐波检测方法的检测效果可以满足IEEE929-2000的标准。
Claims (3)
1.基于离散傅里叶变换的谐波检测方法,其特征在于,包括:
在设定周期内,采集电网侧的三相电流,并对采集的三相电流做离散处理;三相电流分别为a相采集电流、b相采集电流和c相采集电流;
构建三相电流的基波有功电流幅值和基波无功电流幅值的计算模型:
其中,XA为采样电流的基波有功电流幅值;XB为采样电流的基波无功电流幅值;n为采样点的序号;N为一个采样周期内的采样点个数;i为采样电流;
通过滑窗迭代方法对采集电网侧的三相电流进行离散傅里叶变换正变换,分别得到三相电流的基波有功电流幅值和基波无功电流幅值:
当采样点的个数小于采样周期内的采样点个数时,三相电流的基波有功电流幅值和基波无功电流幅值的迭代公式:
当采样点的个数大于采样周期内的采样点个数时,三相电流的基波有功电流幅值和基波无功电流幅值的迭代公式:
采用坐标变换计算三相电流的基波有功电流正序分量:
其中,iA为a相基波有功电流正序分量;iB为b相基波有功电流正序分量;iC为c相基波有功电流正序分量;XaA为a相电流的有功幅值;XaB为a相电流的无功幅值;XbA为b相电流的有功幅值;XbB为b相电流的无功幅值;XcA为c相电流的有功幅值;XcB为c相电流的无功幅值;
根据采集的三相电流及基波有功电流正序分量,计算检测的三相谐波电流:
其中,iah为a相的谐波电流;ibh为b相的谐波电流;ich为c相的谐波电流;ia为a相采集电流;ib为b相采集电流;ic为c相采集电流。
3.根据权利要求1所述的基于离散傅里叶变换的谐波检测方法,其特征在于,在设定周期内采集下,采集三相四线制并联型有源电力滤波***中电网侧的三相电流。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710773955.4A CN107632200B (zh) | 2017-08-31 | 2017-08-31 | 基于离散傅里叶变换的谐波检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710773955.4A CN107632200B (zh) | 2017-08-31 | 2017-08-31 | 基于离散傅里叶变换的谐波检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107632200A CN107632200A (zh) | 2018-01-26 |
CN107632200B true CN107632200B (zh) | 2020-05-08 |
Family
ID=61100726
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710773955.4A Active CN107632200B (zh) | 2017-08-31 | 2017-08-31 | 基于离散傅里叶变换的谐波检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107632200B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109358223A (zh) * | 2018-09-29 | 2019-02-19 | 海特尔机电工程技术(马鞍山)有限公司 | 一种滑动窗dft谐波电流检测方法及装置、存储介质 |
CN110011311A (zh) * | 2019-05-16 | 2019-07-12 | 中国测试技术研究院电子研究所 | 一种谐波控制装置 |
CN110208593B (zh) * | 2019-06-10 | 2022-04-08 | 北京利达英杰联合技术有限公司 | 一种用于电气防火限流式保护器的检测电流的方法 |
CN111077370A (zh) * | 2020-01-02 | 2020-04-28 | 哈尔滨理工大学 | 一种改进的递归离散傅里叶变换检测法 |
CN111208340A (zh) * | 2020-02-26 | 2020-05-29 | 泰州学院 | 一种基于傅里叶变换的单相基波无功电流精确检测方法 |
CN111896802B (zh) * | 2020-06-30 | 2023-03-31 | 成都四威功率电子科技有限公司 | 一种频率自适应的采样方法 |
CN114137283A (zh) * | 2021-11-05 | 2022-03-04 | 江苏安方电力科技有限公司 | 一种无锁相环的单相基波有功电流精确检测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1194887A (ja) * | 1997-09-18 | 1999-04-09 | Toenec Corp | リアルタイム高調波解析装置 |
JP2000180484A (ja) * | 1998-12-18 | 2000-06-30 | Nissin Electric Co Ltd | 高調波測定装置 |
CN102393488A (zh) * | 2011-08-24 | 2012-03-28 | 江苏技术师范学院 | 一种谐波分析方法 |
CN104749435A (zh) * | 2015-04-01 | 2015-07-01 | 四川厚天科技股份有限公司 | 一种无累积误差滑窗的dft谐波检测方法 |
CN105203853A (zh) * | 2015-09-11 | 2015-12-30 | 国网智能电网研究院 | 一种大容量高频变压器寄生电容的测量方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101893652B (zh) * | 2010-06-30 | 2012-07-04 | 中南大学 | 一种基于电压矢量空间变换的谐波和无功电流检测方法 |
CN103267895B (zh) * | 2013-04-25 | 2015-04-22 | 国家电网公司 | 一种电网电流的谐波电流检测方法 |
CN104833853B (zh) * | 2015-05-14 | 2017-08-11 | 电子科技大学 | 一种频率自适应的滑窗dft谐波检测方法 |
-
2017
- 2017-08-31 CN CN201710773955.4A patent/CN107632200B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1194887A (ja) * | 1997-09-18 | 1999-04-09 | Toenec Corp | リアルタイム高調波解析装置 |
JP2000180484A (ja) * | 1998-12-18 | 2000-06-30 | Nissin Electric Co Ltd | 高調波測定装置 |
CN102393488A (zh) * | 2011-08-24 | 2012-03-28 | 江苏技术师范学院 | 一种谐波分析方法 |
CN104749435A (zh) * | 2015-04-01 | 2015-07-01 | 四川厚天科技股份有限公司 | 一种无累积误差滑窗的dft谐波检测方法 |
CN105203853A (zh) * | 2015-09-11 | 2015-12-30 | 国网智能电网研究院 | 一种大容量高频变压器寄生电容的测量方法 |
Non-Patent Citations (2)
Title |
---|
基于DFT的周期信号谐波特性测试与仿真研究;马令坤 等;《计算机仿真》;20151130;第32卷(第11期);第207-211页 * |
电力谐波滑窗迭代DFT检测算法的研究与仿真;陆秀令 等;《***仿真学报》;20080731;第20卷(第14期);第3652-3655页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107632200A (zh) | 2018-01-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107632200B (zh) | 基于离散傅里叶变换的谐波检测方法 | |
CN103257271B (zh) | 一种基于stm32f107vct6的微电网谐波与间谐波检测装置及检测方法 | |
CN102621382B (zh) | 一种检测电力***中电气信号频率、相位和幅值的方法 | |
CN106501574A (zh) | 一种新型有源电力滤波器谐波电流检测方法 | |
CN102590646A (zh) | 一种基于相序滤波器的正序、负序与谐波电流检测方法 | |
CN203287435U (zh) | 一种基于stm32f107vct6的微电网谐波与间谐波检测装置 | |
CN104698254B (zh) | 一种电网基波正序电压提取方法及锁相方法 | |
CN103095296B (zh) | 用于svc控制***的新型软件锁相环的实现方法 | |
CN102664545B (zh) | 逆变器电压前馈直流分量的抑制方法及逆变器的控制方法 | |
CN108365617A (zh) | 一种应用于电网电压不平衡及畸变条件下的锁相方法 | |
CN104833853A (zh) | 一种频率自适应的滑窗dft谐波检测方法 | |
CN107863774B (zh) | 谐波指令电流的获取方法、装置、***及可读存储介质 | |
CN103647550A (zh) | 一种用于动态电压无功补偿的锁相环方法 | |
CN110596455B (zh) | 一种工频电参数提取方法、***及计算机可读存储介质 | |
CN113014250A (zh) | 一种可消除直流偏移电压的锁相环及其锁相控制方法 | |
CN110661263B (zh) | 含有自适应延时滤波器的锁频环及基于该锁频环的并网逆变器控制方法 | |
CN107221931A (zh) | 基于有源电力滤波器的z源逆变器并网控制方法 | |
CN205753451U (zh) | 一种基于电网同步的锁相装置 | |
CN113346899B (zh) | 一种基于级联型滤波器的三相并网软件锁相环 | |
CN103176030A (zh) | 配电***的间谐波检测方法 | |
CN112485522B (zh) | 基于电能数据感知的平顶窗函数同步相量测量方法及装置 | |
CN104820129B (zh) | 一种基波正序有功电流的检测方法 | |
CN108957130B (zh) | 一种基于瞬时无功功率理论的电力***谐波检测方法 | |
CN106771597B (zh) | 一种省略Park变换的无锁相环式SAI谐波电流检测方法 | |
CN104076204A (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 |