CN112557751A - 一种基于dft迭代法的谐波参数估计方法 - Google Patents
一种基于dft迭代法的谐波参数估计方法 Download PDFInfo
- Publication number
- CN112557751A CN112557751A CN202011406474.8A CN202011406474A CN112557751A CN 112557751 A CN112557751 A CN 112557751A CN 202011406474 A CN202011406474 A CN 202011406474A CN 112557751 A CN112557751 A CN 112557751A
- Authority
- CN
- China
- Prior art keywords
- frequency
- dft
- harmonic
- iteration
- value
- 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.)
- Granted
Links
Images
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/16—Spectrum analysis; Fourier analysis
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/40—Arrangements for reducing harmonics
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Emergency Protection Circuit Devices (AREA)
- Measuring Phase Differences (AREA)
Abstract
本发明公开了一种基于DFT迭代法的谐波参数估计方法,首先对信号进行采集和预处理,然后构造频率估计的代价函数,通过DFT迭代法对谐波的真实频率进行搜索迭代运算,以代价函数最小化为目标进行优化,最终迭代出精确的谐波参数估计值。该方法能应用于包含多个阶次谐波的复杂信号,考虑并计算了偶次谐波的影响,能得到逼近真实值的估计结果。
Description
技术领域
本发明涉及一种基于离散傅里叶变换的谐波信号频率估计方法,使用了DFT迭代法。
背景技术
随着可再生能源的日益普及和非线性负荷的广泛应用,电力系统的电能质量面临着诸多挑战,电力电子设备的广泛应用会导致严重的谐波污染,威胁到电力系统的安全稳定运行。因此谐波分析是近年来的一个研究热点。
离散傅里叶变换(Discrete Fourier Transform,DFT)算法在静态条件下应用价值良好,拥有物理意义直观明确的优点,已经被广泛运用于谐波测量中。
DFT算法有许多的限制和缺陷,比如频谱泄露、栅栏效应。在信号处理之前先要对信号进行采样,这相当于让信号通过了一个有限长的矩形窗函数,时域的相乘在频域表现为谱的卷积。矩形窗的幅度谱是带有大量旁瓣的采样函数。在相干采样时,实正弦信号的DFT频谱只有两根镜像的谱线,此时频谱泄露不存在。但是实际应用中,相干采样几乎不存在,在非相干采样情况下,窗函数频谱中的旁瓣导致了频谱泄露。栅栏效应源自DFT变换是离散变换,频谱中的每两根谱线之间的内容是未知的,这样由于频谱分辨率的局限,直接由频谱峰值进行频率估计存在误差。经典的基于DFT的算法为了方便分析会忽略负频率的分量,导致长程频谱泄露的影响没被充分考虑,其改进算法也或多或少忽略了负频率分量,但是实际情况中负频率的频谱泄露影响是非常可观的,特别是正负频率谱线靠近的时候,此时算法的性能会劣化。在谐波测量中会导致谐波定位不准,给谐波干扰治理带来困难。
发明内容
针对上述问题,本发明提出一种基于DFT迭代法的谐波参数估计方法,该方法考虑并计算了正负频率的谱叠加,提高了估计性能。
本发明为解决其技术问题,所采用的技术方案为:一种基于DFT迭代法的谐波参数估计方法,其步骤为:
一种基于DFT迭代法的谐波参数估计方法,其步骤包括:
A、在一个测量周期内,采样得到谐波信号,对其进行离散傅里叶变换得到信号的频谱,对频谱的峰值进行定位得到频率的频率估计值;
B、构造频率估计代价函数;
C、以代价函数最小化为目标,以步骤A所得频率估计值为初始值,通过DFT迭代法求频率的估计值;
D、构建DFT表达式列矩阵方程,将步骤C所得频率的估计值代入DFT表达式列矩阵方程,计算谐波的幅度和相位。
优选的,步骤B的具体步骤包括:
用clark变换处理三相电压信号产生输入相量信号
由此将scos(n)的N点DFT表示为
令
当已知两个不同的DFT序列S(k1)和S(k2)时,可构造如下矩阵方程:
μ和v为参考变量,它们的值和k无关,当ω0的值已知的时候,μ的值可由上式求出,该计算过程记为
同样的,另外一个计算v的过程记为
优选的,步骤C的具体步骤包括:
设定迭代次数M,跳出条件TOL,步进距离ò,迭代结果下限δa和上限δb;
f≥TOL时,设定
当f1<f2时,δb=δc;否则δa=δc;
迭代完成后得到频率的估计值
优选的,步骤D的具体步骤包括:
由正弦信号DFT表达式
列出矩阵方程:
有益效果:利用本发明公开的方法可以通过方程得到DFT单元与阶跃变化频率的关系,使用六个不同的DFT单元来消除符号转换的影响,从而实现对例如FSK的频率阶跃变化的单频信号进行高精度的频率估计。本发明能应用于包含多个阶次谐波的复杂信号,考虑并计算了偶次谐波的影响,能得到逼近真实值的估计结果。根据仿真结果可以看出,本发明方法对谐波的参数估计效果良好,可以实现准确的谐波处理定位。
附图说明
图1是本发明所采用DFT迭代法的算法流程图;
图2是无噪情况下,L为2.14,本发明方法面对2、3、4、5、6、7、9、11阶谐波时对基频的估计效果图;
图3是无噪情况下,L为2.14,本发明方法面对11阶谐波时对各个谐波幅度的估计效果图;
图4是无噪情况下,L为2.14,本发明方法面对11阶谐波时对各个谐波相位的估计效果图。
具体实施方式
下面结合是实施例对本发明作进一步的说明。
本发明公开了一种基于DFT迭代法的谐波参数估计方法,首先对信号进行采集和预处理,然后构造频率估计的代价函数,通过DFT迭代法对谐波的真实频率进行搜索迭代运算,以代价函数最小化为目标进行优化,最终迭代出精确的谐波参数估计值。具体步骤如下:
A数据采集与预处理
在一个测量周期内,采样得到谐波信号,对其进行离散傅里叶变换得到信号的频谱,对频谱的峰值进行定位得到频率的粗略的频率估计值;
B构造代价函数
首先由谐波的DFT表达式:
为方便分析,令
当已知两个不同的DFT序列S(k1)和S(k2)时,可以构造下面这个矩阵方程:
μ和v被称为参考变量,它们的值和k无关,当ω0的值已知的时候,μ的值可以由上式求出。把这个计算过程记为
同样的,另外一个计算v的过程记为
C通过DFT迭代法求频率的值
设定迭代次数M,跳出条件TOL,步进距离ò,迭代结果下限δa和上限δb;
f≥TOL时,设定
当f1<f2时,δb=δc;否则δa=δc。
迭代完成后得到频率的估计值
D通过DFT方法计算谐波的幅度和相位
由正弦信号DFT表达式
列出矩阵方程:
由上式可以很容易求出各个谐波的幅度和相位。
本发明中的谐波次数通常为2、3、4、5、6、7、9、11次,每次计算时取其中一个数字;
为进一步说明该迭代方法,通过仿真实验来测试其性能;
对于频率估计仿真,设定DFT点数为128,L为2.14;迭代次数为100次,跳出条件为10^(-3),步进距离为10^(-5),图2展示了该发明在面对2、3、4、5、6、7、9、11次谐波时对基频的估计效果;图3展示了无噪情况下,本发明方法面对11阶谐波时对各个谐波幅度的估计效果图;图4展示了无噪情况下,本发明方法面对11阶谐波时对各个谐波相位的估计效果图。
根据仿真结果可以看出,本发明方法对谐波的参数估计效果良好,可以实现准确的谐波处理定位。
Claims (4)
1.一种基于DFT迭代法的谐波参数估计方法,其特征在于,其步骤为:
A、在一个测量周期内,采样得到谐波信号,对其进行离散傅里叶变换得到信号的频谱,对频谱的峰值进行定位得到频率的频率估计值;
B、构造频率估计代价函数;
C、以代价函数最小化为目标,以步骤A所得频率估计值为初始值,通过DFT迭代法求频率的估计值;
D、构建DFT表达式列矩阵方程,将步骤C所得频率的估计值代入DFT表达式列矩阵方程,计算谐波的幅度和相位。
2.根据权利要求1所述的一种基于DFT迭代法的谐波参数估计方法,其特征在于,步骤B的具体步骤包括:
用clark变换处理三相电压信号产生输入相量信号
由此将scos(n)的N点DFT表示为
令
当已知两个不同的DFT序列S(k1)和S(k2)时,可构造如下矩阵方程:
μ和v为参考变量,它们的值和k无关,当ω0的值已知的时候,μ的值可由上式求出,该计算过程记为
同样的,另外一个计算v的过程记为
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011406474.8A CN112557751B (zh) | 2020-12-03 | 2020-12-03 | 一种基于dft迭代法的谐波参数估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011406474.8A CN112557751B (zh) | 2020-12-03 | 2020-12-03 | 一种基于dft迭代法的谐波参数估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112557751A true CN112557751A (zh) | 2021-03-26 |
CN112557751B CN112557751B (zh) | 2023-07-18 |
Family
ID=75048253
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011406474.8A Active CN112557751B (zh) | 2020-12-03 | 2020-12-03 | 一种基于dft迭代法的谐波参数估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112557751B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113341224A (zh) * | 2021-06-08 | 2021-09-03 | 国网湖南省电力有限公司 | 一种电力系统低频振荡信号测量方法及装置 |
CN114252700A (zh) * | 2021-10-26 | 2022-03-29 | 深圳市锐风电子科技有限公司 | 一种基于正余弦算法的电力谐波检测方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012013883A1 (fr) * | 2010-07-30 | 2012-02-02 | Itron France | Determination de la frequence fondamentale d'un signal periodique incluant des composantes harmoniques |
CN103399203A (zh) * | 2013-08-09 | 2013-11-20 | 重庆大学 | 一种基于复合迭代算法的谐波参数高精度估计方法 |
CN103809023A (zh) * | 2014-01-26 | 2014-05-21 | 西南交通大学 | 基于二分搜索的电网同步谐波相量测量方法 |
CN104833853A (zh) * | 2015-05-14 | 2015-08-12 | 电子科技大学 | 一种频率自适应的滑窗dft谐波检测方法 |
CN107085140A (zh) * | 2017-04-25 | 2017-08-22 | 东南大学 | 基于改进的SmartDFT算法的非平衡系统频率估计方法 |
CN107102255A (zh) * | 2017-05-31 | 2017-08-29 | 太原科技大学 | 单一adc采集通道动态特性测试方法 |
CN107271768A (zh) * | 2017-05-26 | 2017-10-20 | 东南大学 | 一种最小二乘拟合动态频率测量方法 |
CN108037361A (zh) * | 2017-12-05 | 2018-05-15 | 南京福致通电气自动化有限公司 | 一种基于滑动窗dft的高精度谐波参数估计方法 |
CN109587094A (zh) * | 2019-01-03 | 2019-04-05 | 钟祥博谦信息科技有限公司 | 一种基于gd算法的非线性失真信号解调方法及系统 |
CN111398731A (zh) * | 2020-03-09 | 2020-07-10 | 西南交通大学 | 基于多频率-泰勒模型滤除衰减直流的动态相量测量方法 |
-
2020
- 2020-12-03 CN CN202011406474.8A patent/CN112557751B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012013883A1 (fr) * | 2010-07-30 | 2012-02-02 | Itron France | Determination de la frequence fondamentale d'un signal periodique incluant des composantes harmoniques |
CN103399203A (zh) * | 2013-08-09 | 2013-11-20 | 重庆大学 | 一种基于复合迭代算法的谐波参数高精度估计方法 |
CN103809023A (zh) * | 2014-01-26 | 2014-05-21 | 西南交通大学 | 基于二分搜索的电网同步谐波相量测量方法 |
CN104833853A (zh) * | 2015-05-14 | 2015-08-12 | 电子科技大学 | 一种频率自适应的滑窗dft谐波检测方法 |
CN107085140A (zh) * | 2017-04-25 | 2017-08-22 | 东南大学 | 基于改进的SmartDFT算法的非平衡系统频率估计方法 |
CN107271768A (zh) * | 2017-05-26 | 2017-10-20 | 东南大学 | 一种最小二乘拟合动态频率测量方法 |
CN107102255A (zh) * | 2017-05-31 | 2017-08-29 | 太原科技大学 | 单一adc采集通道动态特性测试方法 |
CN108037361A (zh) * | 2017-12-05 | 2018-05-15 | 南京福致通电气自动化有限公司 | 一种基于滑动窗dft的高精度谐波参数估计方法 |
CN109587094A (zh) * | 2019-01-03 | 2019-04-05 | 钟祥博谦信息科技有限公司 | 一种基于gd算法的非线性失真信号解调方法及系统 |
CN111398731A (zh) * | 2020-03-09 | 2020-07-10 | 西南交通大学 | 基于多频率-泰勒模型滤除衰减直流的动态相量测量方法 |
Non-Patent Citations (2)
Title |
---|
ASJA DERVIŠKADI´ 等: "Iterative-Interpolated DFT for Synchrophasor Estimation: A Single Algorithm for P- and M-Class Compliant PMUs", 《IEEE TRANSACTIONS ON INSTRUMENTATION AND MEASUREMENT》 * |
石立新等: "基于迭代法的谐波信号基音的精确估计", 《电子器件》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113341224A (zh) * | 2021-06-08 | 2021-09-03 | 国网湖南省电力有限公司 | 一种电力系统低频振荡信号测量方法及装置 |
CN113341224B (zh) * | 2021-06-08 | 2022-05-24 | 国网湖南省电力有限公司 | 一种电力系统低频振荡信号测量方法及装置 |
CN114252700A (zh) * | 2021-10-26 | 2022-03-29 | 深圳市锐风电子科技有限公司 | 一种基于正余弦算法的电力谐波检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112557751B (zh) | 2023-07-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108037361B (zh) | 一种基于滑动窗dft的高精度谐波参数估计方法 | |
Su et al. | Power harmonic and interharmonic detection method in renewable power based on Nuttall double‐window all‐phase FFT algorithm | |
CN107102255B (zh) | 单一adc采集通道动态特性测试方法 | |
CN112557751B (zh) | 一种基于dft迭代法的谐波参数估计方法 | |
Borkowski et al. | Frequency estimation in interpolated discrete fourier transform with generalized maximum sidelobe decay windows for the control of power | |
Belega et al. | Amplitude and phase estimation of real-valued sine wave via frequency-domain linear least-squares algorithms | |
CN110320400B (zh) | 准同步采样和改进能量算子的电压闪变包络参数提取方法 | |
CN112255457B (zh) | 适用于自动准同期装置的相角差测量方法 | |
CN109541312A (zh) | 一种新能源汇集地区次同步谐波检测方法 | |
Liu et al. | A morlet wavelet-based two-point FIR filter method for phasor estimation | |
CN109669072A (zh) | 一种配电网的自适应同步相量量测方法 | |
Koteswara Rao et al. | Accurate phasor and frequency estimation during power system oscillations using least squares | |
CN116047163A (zh) | 一种电力系统间谐波检测方法及装置 | |
CN112816779B (zh) | 一种解析信号生成的谐波实信号参数估计方法 | |
CN105372492B (zh) | 基于三条dft复数谱线的信号频率测量方法 | |
CN103543331A (zh) | 一种计算电信号谐波和间谐波的方法 | |
CN109541304A (zh) | 基于六项最小旁瓣窗插值的电网高次弱幅值谐波检测方法 | |
Petrović | Frequency and parameter estimation of multi-sinusoidal signal | |
CN112362968A (zh) | 一种基于预调制cdsc与sdft的单相谐波实时提取方法 | |
CN106324342A (zh) | 一种基于查表的谐波检测方法 | |
CN117234176A (zh) | 一种基于扫频试验的控制回路频域辨识方法 | |
Jiao et al. | An approach for electrical harmonic analysis based on interpolation DFT | |
CN112946374B (zh) | 基于卷积窗函数的三相不平衡度检测方法及装置 | |
Firouzjah et al. | A predictive current control method for shunt active filter with windowing based wavelet transform in harmonic detection | |
CN114487589A (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 |