CN108918964A - 一种稀疏性增强的谐波分析方法 - Google Patents

一种稀疏性增强的谐波分析方法 Download PDF

Info

Publication number
CN108918964A
CN108918964A CN201810763295.6A CN201810763295A CN108918964A CN 108918964 A CN108918964 A CN 108918964A CN 201810763295 A CN201810763295 A CN 201810763295A CN 108918964 A CN108918964 A CN 108918964A
Authority
CN
China
Prior art keywords
spectral line
spectral
harmonic
main lobe
frequency
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
Application number
CN201810763295.6A
Other languages
English (en)
Other versions
CN108918964B (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.)
Suqian Power Supply Branch Jiangsu Electric Power Co Ltd
State Grid Corp of China SGCC
State Grid Jiangsu Electric Power Co Ltd
Original Assignee
Suqian Power Supply Branch Jiangsu Electric Power Co Ltd
State Grid Corp of China SGCC
State Grid Jiangsu Electric Power Co Ltd
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 Suqian Power Supply Branch Jiangsu Electric Power Co Ltd, State Grid Corp of China SGCC, State Grid Jiangsu Electric Power Co Ltd filed Critical Suqian Power Supply Branch Jiangsu Electric Power Co Ltd
Priority to CN201810763295.6A priority Critical patent/CN108918964B/zh
Publication of CN108918964A publication Critical patent/CN108918964A/zh
Application granted granted Critical
Publication of CN108918964B publication Critical patent/CN108918964B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

Landscapes

  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Resistance Or Impedance (AREA)
  • Complex Calculations (AREA)

Abstract

本发明属于电能质量分析与控制领域，特别涉及一种稀疏性增强的谐波分析方法。本发明方法首先将复杂谐波进行离散化，其次对离散信号进行加窗截断，然后获取频域内多根谱线的幅值，根据谱线的幅值，对相邻幅值两两做差得到新的序列，新的序列稀疏性特征更显著，采用1范数优化模型，得到无泄漏的精确幅值，进而可以通过传统的基于幅值的参数估计方法进行谐波分析，该方法通过增强稀疏性，提高了谐波参数估计的精度。

Description

一种稀疏性增强的谐波分析方法
技术领域
本发明属于电能质量分析与控制领域，特别涉及一种稀疏性增强的谐波分析方法。
背景技术
实际电网的复杂谐波包含丰富的间谐波成分，且间谐波有一些不同于整数次谐波的特征，如果简单的使用FFT方法将会出现很多问题，这也是复杂平稳谐波分析最大的挑战。相对于整数次谐波，间谐波有如下特点：
1、间谐波频率不是基波频率的整数倍，甚至很难确定其大概周期，所以要实现对间谐波分量的同步采样基本是不可能的。
2、间谐波和整数次谐波/间谐波的频谱可能靠的很近，这就要求复杂谐波的分析方法必须具有很高的分辨率。
3、间谐波的强度很弱，使其更容易受频谱泄露的影响，尤其是当间谐波和整数次谐波靠的很近的时候，这种影响就更为明显。
在IEC规定下，谱线间的频率间隔为△f＝fs/N＝5Hz，Hanning窗的主瓣宽度为20Hz，对于任意的相邻的整数次谐波而言，这已经足以将其分辨开来。再加之Hanning的加权作用，谐波分量的泄露的影响已经非常小，故高精度的估计整数次谐波分量的参数是较容易实现的。不过对于间谐波而言，考虑到间谐波上述的三个特点，FFT的局限性就非常严重，具体体现在以下三方面：
1、对于整数次谐波，选取合适的采样频率和采样点数，通过加窗减小频谱泄露的影响，再利用插值法克服栅栏效应，即可实现整数次谐波参数的高精度估计，但对于间谐波而言，由于间谐波的频率不是基波的整数倍，故即使在同步采样的时候(对整数次谐波而言)，也很难实现对间谐波的高精度分析。
2、由于间谐波的幅值仅为整数次谐波的百分之几或更小，于是临近的比较强的谐波分量的频谱泄露有可能导致出现虚假间谐波分量或者间谐波谱峰消失，当间谐波和谐波间隔比较小的时候，这种现象会更突出。
3、FFT分析时，为了将信号的两个比较接近频率分量分辨开来，必须提高频率分辨率，这可以通过增加采样时间或者采用主瓣比较窄的窗函数来实现。
增加采样时间确实可以达到高的分辨率，但对于间谐波而言，单纯的增加采样时间是不合适的，这是因为实际的间谐波频谱是随时间变化的，且具有一定随机性，比如由电弧产生的间谐波分量，所以增加时长的后果可能导致前后分析的间谐波频谱发生变化，使结果失去意义。
选取具有窄的主瓣的窗函数是也可以提高频率分辨率，但是提高的程度是比较受限的，不可能达到任意的分辨率，另一方面，提高分辨率和降低频谱泄露效应也是矛盾的：如果从降低频谱泄露的角度看，应该选择旁瓣电平低衰减快的窗函数；但旁瓣衰减快的窗函数同时其主瓣宽度也会增大，这时必然就降低了频率分辨率，同时主瓣间互相干扰的程度也增大，这会对参数估计造成更严重的误差。
在国内已经有利用稀疏性进行谐波分析的思路，但这些方法直接利用频域内的幅值，稀疏性不强。
发明内容
为了解决上述背景技术提出的技术问题，本发明旨在提供一种稀疏性增强的谐波分析方法方法，克服传统FFT法分析谐波所存在的弊端，提出一种基于增强稀疏性的谐波分析犯法，降低频谱泄露给频域稀疏性带来的影响，开拓了复杂谐波分析的思路。
本发明是采取以下技术方案实现的：
一种复杂电力谐波的参数估计方法，包括以下步骤：
1)稳态谐波信号x(t)包含间谐波，对信号x(t)进行离散化，得到数字序列x(n)；
2)选择窗函数，对步骤(1)中得到的数字序列x(n)进行加窗截断，得到加窗信号的频谱；
3)根据步骤(2)得到的加窗信号的频谱，设主瓣内最大谱线位于第k0根频谱处，计算第k0根频谱处相邻两根谱线的幅值之差得到相邻谱线差的序列，
h(m)＝y(k0-m)-y(k0-m-1)，
上式中，y(k0-m)与y(k0-m-1)为主瓣内谐波信号第k0-m根和k0-m-1谱线，下标k0-m根和k0-m-1表示主瓣内谱线的索引，m＝0表示主瓣内最高谱峰处谱线，最高谱峰处谱线左侧由近及远的各条谱线依次用下标m＝-1，m＝-2,…，m＝-I表示，最高谱峰处谱线左侧由近及远的各条谱线依次用下标，m＝1,，m＝2,…,，m＝I表示，I为谱峰左/右边相邻谱线数目之差的总个数，I为大于1的自然数；
4)根据相邻谱线差的序列，建立如下优化模型：
min||H||1
s.t H＝FXt
上式中，min表示最小值，H＝[h(-I),h(-I+1),…,h(I-1),h(I)]T，F为傅里叶差矩阵；Xt为时域信号序列；
求解步骤(4)的非线性规划问题，得到向量H；
5)根据主瓣内谱峰的位置和H，分别求主瓣内谱线的去泄露谱线
y(k0-1)＝y(k0)-h(0)
y(k0-2)＝y(k0)-h(0)-h(1)
6)对去泄露谱线Y＝[y(-k0-I),...y(k0-1),y(k0),y(k0+1),...y(k0+I)]，应用经典的已有谐波参数估方法得到谐波参数A、δ以及f。
步骤4)中，求得向量H的具体方法为：
4-1)根据傅里叶变换
当频谱两两做差之后，
即H＝FXt，则，
显然的稀疏性增强，谐波信号对应的谱峰更显著；
4-2)通过求解获得H，即进而通过主瓣内最高谱线和可重构频域内谱线，降低频谱泄露带来的影响，提高谐波参数估计的精度。
传统的谐波参数估计方法如比值法或者能量重心法认为其他泄露分量的影响有限，当两个谐波信号在频域内靠的很近时，这种假设显然不成立；而已有的基于稀疏的谐波参数估计方法认为频率稀疏，从而将该问题转化成1范数优化的压缩感知问题，然后这种方法的问题随着频域内谐波分量的靠近而存在适配现象，导致稀疏性和降低和求解精度的下降。
本发明为了提高频域内稀疏性，采用相邻谱线相减的形式进行求解，其有益之处在于：
在某谐波信号的频谱集中的主瓣或者旁瓣内，频域幅值衰减很快，以主瓣内最大谱峰右边的谱线为例，谱峰y(k0)为最大，[y(k0+1),y(k0+2),...,y(k0+I)]序列逐渐减小至0，且相邻两谱线差越来越小，那么谱线通过两两之差，显然可以进一步凸显频域幅值的稀疏性；
已有的压缩感知方法，由于适配现象，当存在多个靠的较近的谐波，导致估计精度的下降，而两两做差，使频域内谱峰的效应更明显，有助于参数估计。
具体实施方式
本发明一种复杂电力谐波的参数估计方法，包括以下步骤：
(1)稳态谐波信号x(t)包含间谐波，对信号x(t)进行离散化，得到数字序列x(n)；
(2)选择窗函数，对数字序列x(n)进行加窗截断，得到加窗信号的频谱；
(3)根据步骤(2)得到的加窗信号的频谱，设主瓣内最大谱线位于第k0根频谱处，计算处相邻两根谱线的幅值之差得到：
h(m)＝y(k0-m)-y(k0-m-1)
上式中，y(k0-m)与y(k0-m-1)为主瓣内谐波信号第k0-m根和k0-m-1谱线，下标k0-m根和k0-m-1表示主瓣内谱线的索引，m＝0表示主瓣内最高谱峰处谱线，最高谱峰处谱线左侧由近及远的各条谱线依次用下标m＝-1，m＝-2,…，m＝-I表示，最高谱峰处谱线左侧由近及远的各条谱线依次用下标，m＝1,，m＝2,…,，m＝I表示，I为谱峰左/右边相邻谱线数目之差的总个数；
(4)根据相邻谱线差的序列，建立如下优化模型：
min||H||1
s.t H＝FXt
上式中，min表示最小值，H＝[h(-I),h(-I+1),…,h(I-1),h(I)]T，F为傅里叶差矩阵，具体形式见下文；Xt为时域信号序列；
求解步骤(4)的非线性规划问题，得到向量H。
(5)根据主瓣内谱峰的位置和H，分别求主瓣内谱线的去泄露谱线
(6)对去泄露谱线Y＝[y(-k0-I),...y(k0-1),y(k0),y(k0+1),...y(k0+I)]，应用经典的已有谐波参数估方法得到谐波参数A、δ以及f。
以下将通过具体实施例对本发明的技术方案进行详细说明。
一种复杂电力谐波的参数估计方法，具体过程如下。
设一稳态谐波信号为的形式，以采样频率fs＝1/Ts对x(t)离散化，得数字序列：
其中，A为信号幅值，f0为信号频率，为信号初相，fs为离散化的采样频率，Ts为采样周期。
选择窗函数w(n)，对信号x(n)加窗截断，只考虑在正频率f0附近的频谱，则加窗信号x(n)w(n)的频谱为：
X(k)＝AejθW(k-k0')
其中，W(*)为窗函数的连续频谱函数，N为采样点数，k是频谱所对应的数字位置，k0'＝f0/△f，Δf为频率分辨率，设频谱最高谱峰点为k0，则k0＝k0'+δ，δ为由非同步采样引起的频率偏移量。
为叙述方便，这里以主瓣内谱峰处相邻的三根谱线为例进行分析，三根谱线的对应的幅值分别为y1＝|X(k1)|、y0＝|X(k0)|以及y-1＝|X(k-1)|，设δ为由非同步采样的引起的频率偏移量，如果所加的窗函数为Hanning窗，则无其他分量干扰时：
当无其他分量干扰时，上述三个等式严格相等，但由于间谐波的存在，三个等式是无法严格满足相等的，即三根谱线实际值应该是：
其中，ε1、εmax、ε2是由于其他分量泄露引起的未知干扰量。
根据傅里叶变换：
当频谱两两做差之后：
即H＝FXt，则：
显然的稀疏性增强，谐波信号对应的谱峰更显著。上述实施例仅为说明本发明的技术思想，不能以此限定本发明的保护范围，凡是按照本发明提出的技术思想，在技术方案基础上所做的任何改动，均落入本发明保护范围之内。

Claims (2)

1.一种复杂电力谐波的参数估计方法，其特征在于，包括以下步骤：
1)稳态谐波信号x(t)包含间谐波，对信号x(t)进行离散化，得到数字序列x(n)；
2)选择窗函数，对步骤(1)中得到的数字序列x(n)进行加窗截断，得到加窗信号的频谱；
3)根据步骤(2)得到的加窗信号的频谱，设主瓣内最大谱线位于第k0根频谱处，计算第k0根频谱处相邻两根谱线的幅值之差得到相邻谱线差的序列，
h(m)＝y(k0-m)-y(k0-m-1)，
上式中，y(k0-m)与y(k0-m-1)为主瓣内谐波信号第k0-m根和k0-m-1谱线，下标k0-m根和k0-m-1表示主瓣内谱线的索引，m＝0表示主瓣内最高谱峰处谱线，最高谱峰处谱线左侧由近及远的各条谱线依次用下标m＝-1，m＝-2,…，m＝-I表示，最高谱峰处谱线左侧由近及远的各条谱线依次用下标，m＝1,，m＝2,…,，m＝I表示，I为谱峰左/右边相邻谱线数目之差的总个数，I为大于1的自然数；
4)根据相邻谱线差的序列，建立如下优化模型，
min||H||1
s.t H＝FXt，
上式中，min表示最小值，H＝[h(-I),h(-I+1),…,h(I-1),h(I)]T，F为傅里叶差矩阵，Xt为时域信号序列；
求解步骤(4)的非线性规划问题，得到向量H；
5)根据主瓣内谱峰的位置和H，分别求主瓣内谱线的去泄露谱线
y(k0-1)＝y(k0)-h(0)
y(k0-2)＝y(k0)-h(0)-h(1)
6)对去泄露谱线Y＝[y(-k0-I),...y(k0-1),y(k0),y(k0+1),...y(k0+I)]，应用经典的已有谐波参数估方法得到谐波参数A、δ以及f。
2.根据权利要求1所述的复杂电力谐波的参数估计方法，其特征在于，步骤4)中，求得向量H的具体方法为：
4-1)根据傅里叶变换
当频谱两两做差之后，
即H＝FXt，则，
显然的稀疏性增强，谐波信号对应的谱峰更显著；
4-2)通过求解获得H，即进而通过主瓣内最高谱线和能重构频域内谱线，降低频谱泄露带来的影响，提高谐波参数估计的精度。
CN201810763295.6A 2018-07-11 2018-07-11 一种稀疏性增强的谐波分析方法 Active CN108918964B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810763295.6A CN108918964B (zh) 2018-07-11 2018-07-11 一种稀疏性增强的谐波分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810763295.6A CN108918964B (zh) 2018-07-11 2018-07-11 一种稀疏性增强的谐波分析方法

Publications (2)

Publication Number Publication Date
CN108918964A true CN108918964A (zh) 2018-11-30
CN108918964B CN108918964B (zh) 2020-06-26

Family

ID=64410491

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810763295.6A Active CN108918964B (zh) 2018-07-11 2018-07-11 一种稀疏性增强的谐波分析方法

Country Status (1)

Country Link
CN (1) CN108918964B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109709378A (zh) * 2019-01-21 2019-05-03 国网安徽省电力有限公司六安供电公司 瞬变电信号的频率和幅值自适应算法
CN110045184A (zh) * 2019-04-02 2019-07-23 国网上海市电力公司 一种基于压缩感知macsmp的超次谐波测量方法
CN110716088A (zh) * 2019-10-30 2020-01-21 国网上海市电力公司 一种基于压缩感知macsmp的超高次谐波测量方法
CN111024650A (zh) * 2019-12-10 2020-04-17 中南大学 一种基于信号稀疏重构的玻璃药瓶内气体浓度检测方法及装置
CN111458562A (zh) * 2020-04-15 2020-07-28 国网江苏省电力有限公司宿迁供电分公司 一种基于差商的谐波分析中谱峰可信度判断方法
CN112269962A (zh) * 2020-09-30 2021-01-26 中国船舶重工集团公司七五0试验场 一种辐射噪声线谱稀疏反演方法
CN113971258A (zh) * 2021-09-16 2022-01-25 清华大学 基于深度学习的信号频谱增强方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4542336A (en) * 1983-05-11 1985-09-17 Simmonds Precision Products, Inc. Method and apparatus for sampling broad band spectra in fuel quantity measurement systems
CN106526569A (zh) * 2016-10-28 2017-03-22 中国人民解放军信息工程大学 基于交替迭代的宽带mimo雷达稀疏频谱波形设计方法
CN106788449A (zh) * 2017-01-13 2017-05-31 重庆邮电大学 一种实时电能质量数据参数化压缩方法
CN107576866A (zh) * 2017-08-22 2018-01-12 西南交通大学 一种基于光滑逼近稀疏重构与内点法的多谐波源辨识方法
CN108107269A (zh) * 2017-12-07 2018-06-01 中国矿业大学 一种谐波分析中幅值求解方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4542336A (en) * 1983-05-11 1985-09-17 Simmonds Precision Products, Inc. Method and apparatus for sampling broad band spectra in fuel quantity measurement systems
CN106526569A (zh) * 2016-10-28 2017-03-22 中国人民解放军信息工程大学 基于交替迭代的宽带mimo雷达稀疏频谱波形设计方法
CN106788449A (zh) * 2017-01-13 2017-05-31 重庆邮电大学 一种实时电能质量数据参数化压缩方法
CN107576866A (zh) * 2017-08-22 2018-01-12 西南交通大学 一种基于光滑逼近稀疏重构与内点法的多谐波源辨识方法
CN108107269A (zh) * 2017-12-07 2018-06-01 中国矿业大学 一种谐波分析中幅值求解方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张凤珍等: "基于差分的稀疏度自适应重构算法", 《计算机辅助设计与图形学学报》 *
陈允杰等: "一种基于协同稀疏和全变差的高光谱线性解混方法", 《自动化学报》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109709378A (zh) * 2019-01-21 2019-05-03 国网安徽省电力有限公司六安供电公司 瞬变电信号的频率和幅值自适应算法
CN109709378B (zh) * 2019-01-21 2021-02-19 国网安徽省电力有限公司六安供电公司 瞬变电信号的频率和幅值自适应算法
CN110045184A (zh) * 2019-04-02 2019-07-23 国网上海市电力公司 一种基于压缩感知macsmp的超次谐波测量方法
CN110716088A (zh) * 2019-10-30 2020-01-21 国网上海市电力公司 一种基于压缩感知macsmp的超高次谐波测量方法
CN111024650A (zh) * 2019-12-10 2020-04-17 中南大学 一种基于信号稀疏重构的玻璃药瓶内气体浓度检测方法及装置
CN111458562A (zh) * 2020-04-15 2020-07-28 国网江苏省电力有限公司宿迁供电分公司 一种基于差商的谐波分析中谱峰可信度判断方法
CN112269962A (zh) * 2020-09-30 2021-01-26 中国船舶重工集团公司七五0试验场 一种辐射噪声线谱稀疏反演方法
CN112269962B (zh) * 2020-09-30 2024-05-31 中国船舶重工集团公司七五0试验场 一种辐射噪声线谱稀疏反演方法
CN113971258A (zh) * 2021-09-16 2022-01-25 清华大学 基于深度学习的信号频谱增强方法及系统

Also Published As

Publication number Publication date
CN108918964B (zh) 2020-06-26

Similar Documents

Publication Publication Date Title
CN108918964A (zh) 一种稀疏性增强的谐波分析方法
CN106771591B (zh) 一种复杂电力谐波的参数估计方法
CN106483374B (zh) 一种基于Nuttall双窗全相位FFT的谐波间谐波检测方法
CN103454497B (zh) 基于改进加窗离散傅立叶变换的相位差测量方法
CN107085140B (zh) 基于改进的SmartDFT算法的非平衡系统频率估计方法
CN101915874B (zh) 一种基于傅立叶变换的谐波检测方法
CN109521275B (zh) 一种同步相量确定方法、系统、装置及可读存储介质
CN113378661B (zh) 一种基于改进小波阈值和相关检测的直流电能信号去噪方法
CN107102255B (zh) 单一adc采集通道动态特性测试方法
CN102288807A (zh) 一种测量电网电压闪变的方法
CN103353550A (zh) 一种测量电力系统信号频率及谐波参数的方法
CN103308804A (zh) 基于快速k-s变换电能质量扰动信号时频参数提取方法
CN107085144A (zh) 一种快速测量电力系统谐波的方法
CN103393435B (zh) 一种胎心音信号包络的取得方法及装置
Zeng et al. Parameter estimation of power system signals based on cosine self-convolution window with desirable side-lobe behaviors
CN109900959B (zh) 一种动态正弦畸变信号中谐波成分的提取方法
CN111693775A (zh) 一种输电网的谐波检测方法、装置和介质
CN105137181A (zh) 基于Nuttall-Kaiser组合窗双谱线插值的谐波分析方法
CN110940877A (zh) 一种基于加窗fft插值算法的电压闪变检测方法
CN107315714B (zh) 一种去卷积功率谱估计方法
CN109655665A (zh) 基于布莱克曼窗的全相位傅里叶谐波分析方法
CN110320400A (zh) 准同步采样和改进能量算子的电压闪变包络参数提取方法
CN105486921A (zh) 凯撒三阶互卷积窗三谱线插值的谐波与间谐波检测方法
CN109884390A (zh) 一种电力系统电压谐波次同步/超同步相量测量方法
CN108390698A (zh) 一种基于插值fft算法的电力线载波参数测量方法

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