CN103941090B - 基于谱线能量插值的谐波测量方法 - Google Patents
基于谱线能量插值的谐波测量方法 Download PDFInfo
- Publication number
- CN103941090B CN103941090B CN201410161877.9A CN201410161877A CN103941090B CN 103941090 B CN103941090 B CN 103941090B CN 201410161877 A CN201410161877 A CN 201410161877A CN 103941090 B CN103941090 B CN 103941090B
- Authority
- CN
- China
- Prior art keywords
- frequency
- window
- spectral
- spectral line
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 230000003595 spectral effect Effects 0.000 claims abstract description 101
- 238000001228 spectrum Methods 0.000 claims description 54
- 238000005070 sampling Methods 0.000 claims description 41
- 238000000691 measurement method Methods 0.000 claims description 16
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 6
- 238000010183 spectrum analysis Methods 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 230000001360 synchronised effect Effects 0.000 claims description 4
- 229910000831 Steel Inorganic materials 0.000 claims description 3
- 238000011084 recovery Methods 0.000 claims description 3
- 239000010959 steel Substances 0.000 claims description 3
- 238000001514 detection method Methods 0.000 abstract description 22
- 238000004422 calculation algorithm Methods 0.000 abstract description 19
- 238000006243 chemical reaction Methods 0.000 abstract description 5
- 230000000694 effects Effects 0.000 abstract description 5
- 230000033228 biological regulation Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Landscapes
- Complex Calculations (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
本发明公开了一种基于谱线能量插值的谐波测量方法,满足IEC标准规定的10个周波检测环境,精确检测出某分量的具体参数,而IEC方法只能粗略地检测出信号某频段内的幅值;首次从信号频域谱线能量的角度,通过乘以旋转相位因子,进行能量谱线的变换和插值,对干扰分量能量的相互抵消来消除其旁瓣叠加泄漏的影响,提高频率、幅值和相位的检测精度;当信号中的分量可能会淹没在噪声信号中,通过谱线能量变换将主要信号突出,从而有效地从噪声信号中提取有用信号,源于IEC标准推荐的DFT/FFT算法,原理简单,效果及稳定性较好,具有一定的工程运用价值和良好的应用前景。
Description
技术领域
本发明涉及一种基于谱线能量插值的谐波测量方法,属于电能质量分析技术领域。
背景技术
电能质量的优劣对电网的稳定坚强运行有着至关重要的影响,而谐波、间谐波参数指标是衡量供电质量的一项重要参数,随着国内蓬勃发展的高铁、电动汽车、风电和光伏发电等新能源的接入,能否精确地检测谐波、间谐波参数对其评估影响较大。国际电工委员会(International Electrotechnical Commission,IEC)出台的标准IEC61000-4-7和IEC61000-4-30中规定采用DFT/FFT进行频谱分析,并规定了分析窗口长度为10个周波(50Hz***),IEC标准规定了同步采样时进行10周波矩形窗的频谱分组检测方法,优点是可以快速地进行谐波、间谐波含量分频段检测,缺点是不能测出某一频率分量的具体参数,且在非同步采样时有较大的测量误差。
针对IEC标准中的缺陷,目前有很多的基于DFT/FFT改进算法,常用的两类算法是加窗插值算法和相位差算法,加窗插值算法一般是利用分量频率附近的最大和次大两根谱线的幅值进行插值计算,其精度与所加窗的类型有很大的关系,窗函数的旁瓣抑制性能越好其检测精度越高,但同时主瓣宽度越大,目前的方法中采用加窗插值算法会增加数据的长度,使分析窗口长度达到几十个周波,不能满足IEC标准的要求;相位差算法采用的是基于两段不同时刻的数据进行频谱分析,通过对应谱线的相位差求取分量的参数,因此采样长度至少为20个周波才能满足IEC标准的5Hz分辨率标准,同时,非同步采样带来的旁瓣泄漏对检测结果影响较大,上述的两种算法在计算对应谱线频谱值的时候,忽略了其它分量的旁瓣泄漏干扰,当频率分量相隔较近或者泄漏较严重时,会产生较大的误差。
发明内容
本发明的目的是克服现有技术中针对IEC标准中的缺陷,常用的两类算法是加窗插值算法和相位差算法存在的不足。本发明的基于谱线能量插值的谐波测量方法,能够解决IEC标准下10个周波数据长度的谐波、间谐波精确检测问题,同时使其能够消除其余分量的旁瓣干扰,以进一步提高谐波间谐波的检测精确度,源于IEC标准推荐的DFT/FFT算法,效果及稳定性较好,具有一定的工程运用价值和良好的应用前景。
为了达到上述的目的,本发明所采用的技术方案是:
一种基于谱线能量插值的谐波测量方法,其特征在于:包括以下步骤,
步骤(1),对于电网信号的离散采样,采样值g(nΔt)N为公式(1),
其中,Δt为单位采样间隔,M为采样分量个数,fi、Ai、为采样分量的频率、幅值和相位,n=0、1、…、N-1,N为采样窗口长度;
步骤(2),对步骤(1)的采样值g(nΔt)N按照IEC标准规定的10周波加窗,并根据DFT/FFT频谱分析得到频谱G(k),
其中,βi=fi/Δf为单位化频率值,Δf为采样窗长为N时对应的频率分辨率,W(k-βi)表示对应谱线(k-βi)处的窗谱值,W(k+βi)表示对应谱线(k+βi)处的窗谱值,G(k)对应谱线k处的频谱值;
步骤(3),将步骤(2)得到的频谱G(k)进行乘旋转相位因子处理得到新频谱
步骤(4),通过对新频谱相邻谱线相加以消除部分旁瓣叠加泄漏,求取对应的频率偏移值得到采样分量的频率fi;
步骤(5),根据新频谱作频谱平方处理,转化为谱线的能量,通过对采样分量的频率fi处的三根谱线能量相加来消除旁瓣叠加能量干扰,插值求得信号的幅值Ai和相位值
前述的基于谱线能量插值的谐波测量方法,其特征在于:步骤(2)对步骤(1)的采样值g(nΔt)N按照IEC标准规定的10周波加窗,当基频同步时加10周波的矩形窗,当非同步时加10周波的Hanning窗,当采样窗口长度N>>1时,对应窗谱函数的表达式如下:矩形窗为,
Hanning窗为,
前述的基于谱线能量插值的谐波测量方法,其特征在于:步骤(3)将频谱G(k)进行乘旋转相位因子处理得到新频谱的方法为,
(1)将公式(2)表示的频谱G(k)转换为公式(5),
其中,频谱G(k)分量的相位因子中含两项,其随着谱线位置的不同而变化,当k值的增大影响频谱G(k)分量的相位因子增大;
(2)乘以旋转相位因子以抵消公式(5)中的相位值随谱线位置k的变化而变化,乘旋转相位因子为代入公式(5)得到新频谱
(3)将公式(6)中的单位化频率值βi写成如公式(7)的形式,
其中,ki、为对应的采样分量的频率fi的最大谱线序号和频率偏移值,新频谱表示为,
即为,
其中,为乘旋转相位因子后频谱上其余分量对采样分量频率fi的旁瓣干扰叠加值。
前述的基于谱线能量插值的谐波测量方法,其特征在于:步骤(4),通过对新频谱相邻谱线相加以消除部分旁瓣叠加泄漏,求取对应的频率偏移值得到采样分量的频率fi的方法为,
(1)根据公式(9),通过对新频谱相邻谱线相加以消除部分旁瓣叠加泄漏,令Δ(ki)+Δ(ki+1)≈0,则定义比值因子α为公式(10),
(2)对于矩形窗,根据公式(3)有公式(10)化简为,
其中,αR为矩形窗对应的比值因子α;
(3)对于Hanning窗,根据公式(4)有公式(10)化简为,
其中,αH为Hanning窗对应的比值因子α;
(4)根据公式(13),反解出频率偏移值
其中,cw为对应窗函数的幅值恢复系数,加矩形窗时cw=1,加Hanning窗时cw=2;为对应谱线(ki-1)、ki、(ki+1)处的新频谱值;
(5)根据公式(14),计算采样分量的频率fi为,
前述的基于谱线能量插值的谐波测量方法,其特征在于:步骤(5),插值求得信号的幅值Ai和相位值的方法为,
(1)根据帕塞瓦尔定理,利用对应分量最大能量谱线处的三根谱线对应的频谱值进行平方处理转化为谱线能量,
其中,公式(15)、(16)、(17)中为对应分量最大能量谱线处的三根谱线对应的频谱值;
(2)为了消除其它分量旁瓣能量的干扰,对公式(15)、(16)、(17)的三根谱线能量相加为,
(3)根据乘以相位旋转因子后的旁瓣特性,做近似处理,公式(18)化简为,
(4)根据公式(19),得到加矩形窗或Hanning窗的幅值、相位信息的表达式,
对于矩形窗为,
对于Hanning窗为,
其中,
(5)根据公式(20)或(21),得到幅值Ai和相位值
本发明的有益效果是:本发明的基于谱线能量插值的谐波测量方法,与现有技术相比,具有以下优点:
1)本发明严格满足IEC标准规定的10个周波检测环境,可以精确检测出某分量的具体参数,而IEC方法只能粗略地检测出信号某频段内的幅值;
2)本发明首次从信号频域谱线能量的角度,通过乘以旋转相位因子,进行能量谱线的变换和插值,对干扰分量能量的相互抵消来消除其旁瓣叠加泄漏的影响,提高频率、幅值和相位的检测精度;
(3)当信号中的分量可能会淹没在噪声信号中,本发明通过谱线能量变换将主要信号突出,从而有效地从噪声信号中提取有用信号;
(4)本发明源于IEC标准推荐的DFT/FFT算法,原理简单,效果及稳定性较好,具有一定的工程运用价值和良好的应用前景。
附图说明
图1是本发明的基于谱线能量插值的谐波测量方法的流程图。
图2是本发明在IEC标准下加矩形窗时对应的频谱图。
图3是本发明在IEC标准下加Hanning窗时对应的频谱图。
图4是本发明的IEC标准下不同信噪比时的频率检测误差图。
图5是本发明的IEC标准下不同信噪比时的幅值检测误差图。
图6是本发明的IEC标准下不同信噪比时的相位检测误差图。
具体实施方式
下面将结合说明书附图,对本发明作进一步说明。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
本发明的基于谱线能量插值的谐波测量方法,能够解决IEC标准下10个周波数据长度的谐波、间谐波精确检测问题,同时使其能够消除其余分量的旁瓣干扰,通过谱线能量变换将主要信号突出,从而有效地从噪声信号中提取有用信号,以进一步提高谐波间谐波的检测精确度,源于IEC标准推荐的DFT/FFT算法,效果及稳定性较好,如图1所示,包括以下步骤,
步骤(1),对于电网信号的离散采样,采样值g(nΔt)N为公式(1),
其中,Δt为单位采样间隔,M为采样分量个数,fi、Ai、为采样分量的频率、幅值和相位,n=0、1、…、N-1,N为采样窗口长度;
步骤(2),对步骤(1)的采样值g(nΔt)N按照IEC标准规定的10周波加窗,并根据DFT/FFT频谱分析得到频谱G(k),
其中,βi=fi/Δf为单位化频率值,Δf为采样窗长为N时对应的频率分辨率,W(k-βi)表示对应谱线(k-βi)处的窗谱值,W(k+βi)表示对应谱线(k+βi)处的窗谱值,G(k)对应谱线k处的频谱值;
当基频同步时加10周波的矩形窗,当非同步时加10周波的Hanning窗,当采样窗口长度N>>1时,对应窗谱函数的表达式如下:矩形窗为,
Hanning窗为,
步骤(3),将步骤(2)得到的频谱G(k)进行乘旋转相位因子处理得到新频谱使得各分量主瓣内的相位严格相同,旁瓣内的相位严格相反、幅值相近,得到新频谱的方法为,
(1)将公式(2)表示的频谱G(k)转换为公式(5),
其中,频谱G(k)分量的相位因子中含两项,其随着谱线位置的不同而变化,当k值的增大影响频谱G(k)分量的相位因子增大;
(2)乘以旋转相位因子以抵消公式(5)中的相位值随谱线位置k的变化而变化,乘旋转相位因子为代入公式(5)得到新频谱
(3)将公式(6)中的单位化频率值βi写成如公式(7)的形式,
其中,ki、为对应的采样分量的频率fi的最大谱线序号和频率偏移值,新频谱表示为,
即为,
其中,为乘旋转相位因子后频谱上其余分量对采样分量频率fi的旁瓣干扰叠加值,对单个值Δ(km)来说,其旁瓣幅值相近、相位方向严格相反;
步骤(4),通过对新频谱相邻谱线相加以消除部分旁瓣叠加泄漏,求取对应的频率偏移值得到采样分量的频率fi,计算过程如下,
(1)根据公式(9),通过对新频谱相邻谱线相加以消除部分旁瓣叠加泄漏,令Δ(ki)+Δ(ki+1)≈0,则定义比值因子α为公式(10),
(2)对于矩形窗,根据公式(3)有公式(10)化简为,
其中,αR为矩形窗对应的比值因子α;
(3)对于Hanning窗,根据公式(4)有公式(10)化简为,
其中,αH为Hanning窗对应的比值因子α;
(4)根据公式(13),反解出频率偏移值
其中,cw为对应窗函数的幅值恢复系数,加矩形窗时cw=1,加Hanning窗时cw=2;为对应谱线(ki-1)、ki、(ki+1)处的新频谱值;
(5)根据公式(14),计算采样分量的频率fi为,
步骤(5),根据新频谱作频谱平方处理,转化为谱线的能量,信号频域上的每根谱线是蕴含能量信息的,考虑每根谱线所含能量不同,通过对采样分量的频率fi处的三根谱线能量相加来消除旁瓣叠加能量干扰,插值求得信号的幅值Ai和相位值计算过程如下,
(1)根据帕塞瓦尔定理,利用对应分量最大能量谱线处的三根谱线对应的频谱值进行平方处理转化为谱线能量,
其中,公式(15)、(16)、(17)中为对应分量最大能量谱线处的三根谱线对应的频谱值;
(2)为了消除其它分量旁瓣能量的干扰,对公式(15)、(16)、(17)的三根谱线能量相加为,
(3)根据乘以相位旋转因子后的旁瓣特性,做近似处理,公式(18)化简为,
(4)根据公式(19),得到加矩形窗或Hanning窗的幅值、相位信息的表达式,
对于矩形窗为,
对于Hanning窗为,
其中,
(5)根据公式(20)或(21),得到幅值Ai和相位值
以上步骤(1)-步骤(5),即实现了在IEC标准规定检测要求下,10个周波内的从信号谱线能量的角度进行谐波间谐波参数检测,下面根据一具体实施例,对本发明进行算法仿真验证,具体如下,
主动配电网进行仿真分析,IEC 61000-4-7标准规定谐波间谐波检测时的数据长度为10个周波,采样DFT/FFT方法进行分析,设含有多个频率分量的信号(表1),采样频率为fs=10240Hz,采样周期为10个周波,即2048点,对应的频率分辨率Δf=5Hz,
表1各分量参数理论值
对上述的信号分别进行本发明前述的加矩形窗和Hanning窗分析,并进行频谱转换,其谱线分布图如图2和图3所示,运用谱线能量插值方法检测的结果如表2所示,
表2本发明算法的检测结果
从表2中可以看出,本发明均可以在IEC标准规定下实现谐波、间谐波参数的较为准确的检测;在基波频率非同步采样时,基于Hanning窗的谱线能量插值方法(方法2)比基于矩形窗的谱线能量插值方法(方法1)具有更高的检测精度,方法2更适合非同步采样时的谐波间谐波参数具体测量。
以加Hanning窗的信号为例,为了验证本发明所提方法的有效性,对背景技术中两种主流方案:基于Hanning窗的双谱线插值法(方法3)、基于Hanning窗的相位差法(方法4)的检测结果进行对比,同时限定各方法的采样数据长度均为10个周波,方法3和方法4的检测结果如表3所示。
表3方法3和方法4的检测结果
为验证本发明方法的抗噪性能,同时与方法3和方法4的测量误差进行对比,模拟噪声环境为叠加10~100dB信噪比(SNR)的高斯白噪声信号,分别采用上述4种方法的检测结果进行对比,得到不同噪声强度下各算法的检测误差变化情况(以82Hz间谐波为例)如图4、图5和图6所示,均表明非同步采样时,采用加Hanning窗的测量效果要好于加矩形窗,本发明方法中加Hanning窗的谱线能量插值方法(方法2)比加Hanning窗的双谱线插值方法(方法3)、加Hanning窗的相位差方法(方法4)测量精度高1~2个数量级,由此验证了本发明中所提方法在IEC标准下的准确性和有效性,具有较强的抗干扰性能。
以上显示和描述了本发明的基本原理、主要特征及优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (5)
1.基于谱线能量插值的谐波测量方法,其特征在于:包括以下步骤,
步骤(1),对于电网信号的离散采样,采样值g(nΔt)N为公式(1),
其中,Δt为单位采样间隔,M为采样分量个数,fi、Ai、为采样分量的频率、幅值和相位,n=0、1、…、N-1,N为采样窗口长度;
步骤(2),对步骤(1)的采样值g(nΔt)N按照IEC标准规定的10周波加窗,并根据DFT/FFT频谱分析得到频谱G(k),
其中,βi=fi/Δf为单位化频率值,Δf为采样窗长为N时对应的频率分辨率,W(k-βi)表示对应谱线(k-βi)处的窗谱值,W(k+βi)表示对应谱线(k+βi)处的窗谱值,G(k)对应谱线k处的频谱值;
步骤(3),将步骤(2)得到的频谱G(k)进行乘旋转相位因子处理得到新频谱
步骤(4),通过对新频谱相邻谱线相加以消除部分旁瓣叠加泄漏,求取对应的频率偏移值得到采样分量的频率fi;
步骤(5),根据新频谱作频谱平方处理,转化为谱线的能量,通过对采样分量的频率fi处的三根谱线能量相加来消除旁瓣叠加能量干扰,插值求得信号的幅值Ai和相位值
2.根据权利要求1所述的基于谱线能量插值的谐波测量方法,其特征在于:步骤(2)对步骤(1)的采样值g(nΔt)N按照IEC标准规定的10周波加窗,当基频同步时加10周波的矩形窗,当非同步时加10周波的Hanning窗,当采样窗口长度N>>1时,对应窗谱函数的表达式如下:
矩形窗为,
Hanning窗为,
3.根据权利要求1所述的基于谱线能量插值的谐波测量方法,其特征在于:步骤(3)将频谱G(k)进行乘旋转相位因子处理得到新频谱的方法为,
(1)将公式(2)表示的频谱G(k)转换为公式(5),
其中,频谱G(k)分量的相位因子中含 两项,其随着谱线位置的不同而变化,当k值的增大影响频谱G(k)分量的相位因子增大;
(2)乘以旋转相位因子以抵消公式(5)中的相位值随谱线位置k的变化而变化,乘旋转相位因子为代入公式(5)得到新频谱
(3)将公式(6)中的单位化频率值βi写成如公式(7)的形式,
其中,ki、为对应的采样分量的频率fi的最大谱线序号和频率偏移值,新频谱表示为,
即为,
其中,为乘旋转相位因子后频谱上其余分量对采样分量频率fi的旁瓣干扰叠加值。
4.根据权利要求1所述的基于谱线能量插值的谐波测量方法,其特征在于:步骤(4),通过对新频谱相邻谱线相加以消除部分旁瓣叠加泄漏,求取对应的频率偏移值得到采样分量的频率fi的方法为,
(1)根据公式(9),通过对新频谱相邻谱线相加以消除部分旁瓣叠加泄漏,令Δ(ki)+Δ(ki+1)≈0,则定义比值因子α为公式(10),
(2)对于矩形窗,根据公式(3)有公式(10)化简为,
其中,αR为矩形窗对应的比值因子α;
(3)对于Hanning窗,根据公式(4)有公式(10)化简为,
其中,αH为Hanning窗对应的比值因子α;
(4)根据公式(13),反解出频率偏移值
其中,cw为对应窗函数的幅值恢复系数,加矩形窗时cw=1,加Hanning窗时cw=2; 为对应谱线(ki-1)、ki、(ki+1)处的新频谱值;
(5)根据公式(14),计算采样分量的频率fi为,
5.根据权利要求1所述的基于谱线能量插值的谐波测量方法,其特征在于:步骤(5),插值求得信号的幅值Ai和相位值的方法为,
(1)根据帕塞瓦尔定理,利用对应分量最大能量谱线处的三根谱线对应的频谱值进行平方处理转化为谱线能量,
其中,公式(15)、(16)、(17)中 为对应分量最大能量谱线处的三根谱线对应的频谱值;
(2)为了消除其它分量旁瓣能量的干扰,对公式(15)、(16)、(17)的三根谱线能量相加为,
(3)根据乘以相位旋转因子后的旁瓣特性,做近似处理,公式(18)化简为,
(4)根据公式(19),得到加矩形窗或Hanning窗的幅值、相位信息的表达式,
对于矩形窗为,
对于Hanning窗为,
其中,
(5)根据公式(20)或(21),得到幅值Ai和相位值
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410161877.9A CN103941090B (zh) | 2014-04-22 | 2014-04-22 | 基于谱线能量插值的谐波测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410161877.9A CN103941090B (zh) | 2014-04-22 | 2014-04-22 | 基于谱线能量插值的谐波测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103941090A CN103941090A (zh) | 2014-07-23 |
CN103941090B true CN103941090B (zh) | 2016-11-23 |
Family
ID=51188839
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410161877.9A Active CN103941090B (zh) | 2014-04-22 | 2014-04-22 | 基于谱线能量插值的谐波测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103941090B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104122443B (zh) * | 2014-08-04 | 2017-02-15 | 国家电网公司 | Iec框架下的邻近谐波间谐波分离测量方法 |
CN105372493B (zh) * | 2014-08-31 | 2018-05-08 | 常州昊云工控科技有限公司 | 基于三条dft复数谱线的信号幅值和相位测量方法 |
CN105372492B (zh) * | 2014-08-31 | 2018-06-22 | 常州昊云工控科技有限公司 | 基于三条dft复数谱线的信号频率测量方法 |
CN105588985A (zh) * | 2016-03-18 | 2016-05-18 | 国网上海市电力公司 | 一种基于相位因子变换插值的介损在线精确测量方法 |
CN108241106A (zh) * | 2017-12-29 | 2018-07-03 | 国网北京市电力公司 | 时延估计方法及装置 |
CN108802180B (zh) * | 2018-05-19 | 2021-01-15 | 南京理工大学 | 一种通过非线性混频声波来评估微裂纹宽度的方法 |
CN113030540B (zh) * | 2021-03-01 | 2022-07-26 | 湖南大学 | 一种分布式新能源并网的基波和谐波电能双向计量方法 |
CN115598416B (zh) * | 2022-09-16 | 2024-01-30 | 珠海多创科技有限公司 | 台区采样信号的处理方法、***、存储介质及计算机设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101261292A (zh) * | 2008-04-14 | 2008-09-10 | 湖南大学 | 基于5项Rife-Vincent(I)窗双谱线插值FFT的基波与谐波检测方法 |
CN101441233A (zh) * | 2008-12-02 | 2009-05-27 | 湖南海兴电器有限责任公司 | 基于Kaiser窗双谱线插值FFT的基波与谐波检测方法 |
CN101825660A (zh) * | 2010-05-05 | 2010-09-08 | 天津大学 | 欠采样下的正弦信号频率的高效测量方法及实施装置 |
CN103207319A (zh) * | 2013-03-12 | 2013-07-17 | 广东电网公司电力科学研究院 | 数字化变电站电力信号非同步采样条件下的谐波测量方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6915224B2 (en) * | 2002-10-25 | 2005-07-05 | Jung-Ching Wu | Method for optimum spectrum analysis |
-
2014
- 2014-04-22 CN CN201410161877.9A patent/CN103941090B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101261292A (zh) * | 2008-04-14 | 2008-09-10 | 湖南大学 | 基于5项Rife-Vincent(I)窗双谱线插值FFT的基波与谐波检测方法 |
CN101441233A (zh) * | 2008-12-02 | 2009-05-27 | 湖南海兴电器有限责任公司 | 基于Kaiser窗双谱线插值FFT的基波与谐波检测方法 |
CN101825660A (zh) * | 2010-05-05 | 2010-09-08 | 天津大学 | 欠采样下的正弦信号频率的高效测量方法及实施装置 |
CN103207319A (zh) * | 2013-03-12 | 2013-07-17 | 广东电网公司电力科学研究院 | 数字化变电站电力信号非同步采样条件下的谐波测量方法 |
Non-Patent Citations (3)
Title |
---|
Interharmonics:basic concepts and techniques for their detection and measurement;Chun Li et al.;《Electric Power Systems Research》;20031231(第66期);第39-48页 * |
一种新的电力***谐波间谐波两步检测法;惠锦 等;《电力***保护与控制》;20091201;第37卷(第23期);第28-34页 * |
电力***谐波分析的多层DFT插值校正法;刘亚梅 等;《中国电机工程学报》;20120905;第32卷(第25期);第182-188页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103941090A (zh) | 2014-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103941090B (zh) | 基于谱线能量插值的谐波测量方法 | |
CN103454497B (zh) | 基于改进加窗离散傅立叶变换的相位差测量方法 | |
CN106483374B (zh) | 一种基于Nuttall双窗全相位FFT的谐波间谐波检测方法 | |
CN103869162B (zh) | 一种基于时域准同步的动态信号相量测量方法 | |
CN107247182B (zh) | 一种基于量测相量数据的间谐波分量还原方法 | |
CN104122443B (zh) | Iec框架下的邻近谐波间谐波分离测量方法 | |
CN110837001B (zh) | 一种电力***中谐波和间谐波的分析方法与装置 | |
CN109856455B (zh) | 一种实复转换式衰减信号参数估计方法 | |
CN103308766A (zh) | 一种基于凯撒自卷积窗双谱线插值fft谐波分析方法及其装置 | |
CN103018555B (zh) | 一种高精度的电力参数软件同步采样方法 | |
CN107543962A (zh) | 主导间谐波频谱分布的计算方法 | |
CN106771591A (zh) | 一种复杂电力谐波的参数估计方法 | |
CN106324340A (zh) | 一种同步相量和频率测量动态性能的方法 | |
CN109655665A (zh) | 基于布莱克曼窗的全相位傅里叶谐波分析方法 | |
CN106501602A (zh) | 一种基于滑窗频谱分离的基波参数测量方法 | |
CN103245830B (zh) | 一种结合ar谱估计与非线性优化的间谐波检测方法 | |
Jin et al. | A novel power harmonic analysis method based on Nuttall-Kaiser combination window double spectrum interpolated FFT algorithm | |
CN109541304A (zh) | 基于六项最小旁瓣窗插值的电网高次弱幅值谐波检测方法 | |
CN112557751B (zh) | 一种基于dft迭代法的谐波参数估计方法 | |
Zhang et al. | Frequency shifting and filtering algorithm for power system harmonic estimation | |
CN105372492B (zh) | 基于三条dft复数谱线的信号频率测量方法 | |
CN110320400B (zh) | 准同步采样和改进能量算子的电压闪变包络参数提取方法 | |
CN105467209A (zh) | 一种新的金属氧化物避雷器泄漏电流分析方法 | |
CN112180161A (zh) | 一种非同步高采样率采样条件下谐波间谐波群测量方法 | |
CN112946374B (zh) | 基于卷积窗函数的三相不平衡度检测方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |