CN203012015U - 一种改进傅里叶变换的谐波电能计量*** - Google Patents
一种改进傅里叶变换的谐波电能计量*** Download PDFInfo
- Publication number
- CN203012015U CN203012015U CN 201320032423 CN201320032423U CN203012015U CN 203012015 U CN203012015 U CN 203012015U CN 201320032423 CN201320032423 CN 201320032423 CN 201320032423 U CN201320032423 U CN 201320032423U CN 203012015 U CN203012015 U CN 203012015U
- Authority
- CN
- China
- Prior art keywords
- module
- harmonic
- energy metering
- metering system
- data processing
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Measurement Of Resistance Or Impedance (AREA)
Abstract
一种改进傅里叶变换的谐波电能计量***,由三相电网数据采集模块、数据处理模块、数据管理模块以及***模块组成。其中,三相电网数据采集模块通过并行接口与数据处理模块相连,数据处理模块与数据管理模块通过异步串口相连,数据管理模块与***模块直接相连。采用本实用新型能够很好克服基波频率波动以及间谐波的影响,提高谐波参数计算的准确度。
Description
技术领域
本实用新型涉及一种信号处理技术领域,具体是一种改进傅里叶变换的谐波电能计量***。
背景技术
随着智能电网建设的全面铺开,各种分布式发电和储能设备并网运行,改善电网供电可靠性的同时也使得谐波污染问题更加严重。实际测量中,由于电网频率偏差导致谐波频率的不确定性,电压波动以及间谐波等的影响都会降低谐波分析的准确度。非稳态下的高准确度谐波检测与分析能够为谐波双向潮流计算、谐波电能计量、谐波抑制与补偿等提供科学依据,对提高能源利用率、维护绿色电力环境、确保电网安全运行具有重要意义。
电力谐波分析已经发展了多种基于数字信号处理的时域、频域方法。其中,快速傅里叶变换(快速傅里叶变换)因其计算简单且易于嵌入式***实现而得到了广泛应用。但现有测量仪器多采用恒定采样率,加上间谐波、次谐波等的存在,非稳态下的同步采样难以严格实现,导致的频谱泄漏和栅栏效应影响谐波分析的准确度。窗函数加权可减小频谱泄漏的影响,结合相应的改进算法可在一定程度上提高谐波分析准确度,但在分析非稳态下的谐波信号时其计算误差仍然较大。
发明内容
本实用新型目的在于克服谐波分析在非稳态下导致频谱泄漏和栅栏效应的影响,提出一种改进傅里叶变换的谐波电能计量***,能有效克服基波频率波动以及间谐波的影响,谐波参数计算准确度明显提高,适合于非稳态条件下的谐波分析,从而提高信号谐波分析的准确度和实用性,为进一步进行的信号参数识别与谐波电能计量提供可靠依据。
本实用新型通过以下技术方案实现上述目的:一种改进傅里叶变换的谐波电能计量***,包括三相电网数据采集模块、数据处理模块、数据管理模块以及***模块。其中,三相电网数据采集电路通过并行接口与数据处理电路相连,数据处理电路与数据管理单元通过异步串口相连、数据管理单元与***模块相连。
所述三相电网数据采集模块包括电阻分压网络、TA与电阻网络以及采样芯片ADS8364,用来完成三相电网电压、电流信号的同步采样以及A/D转换。
所述数据处理模块采用信号处理器ADSP-BF533,并且采用基于余弦偶次幂窗改进傅里叶变换的谐波分析方法,对电压、电流波形进行频谱分析,分离出基波及各次谐波分量,完成电参量测量、谐波分析及谐波电能计量。
所述数据管理单元采用单片机M30624FGPFP,完成***的显示、存储、通信以及功能选择等功能。
本实用新型的原理是:
典型余弦函数窗的离散时域表达式为
wcos(n)=sinα(nπ/N),n=0,1,2,...,N-1(1)
式中,整数α的取值决定窗函数的形式。当参数α满足α=2p(p为自然数)时,式(1)形式的余弦函数窗可表示为
式中,n=0,1,2,....,N-1;Jm满足约束条件:
则称此时的窗函数为余弦偶次幂窗,由式(2)可得其离散傅里叶变换(DFT)表达式为
式中,WR(.)为矩形窗的DFT表达式
一般情况下N>>1,可得EOCW窗的DFT简化表达式为
窗函数的主瓣宽度反应其频率分辨率,主瓣越宽,频率分辨率越低;窗函数的频谱泄漏抑制能力可通过其旁瓣峰值电平及旁瓣衰减速率来描述,峰值电平越小且旁瓣衰减速率越大,其频谱泄漏抑制能力越强。
典型EOCW窗的主瓣宽度BW、旁瓣峰值电平Asl以及旁瓣衰减速率Dsl与参数p有如下近似关系:
BW=4(p+1)π/N(7)
Asl≈-16(p+1)(8)
Dsl≈6(2p+1)(9)
含谐波与间谐波分量的多频率信号x(t)经采样率为fs的数据采集***后得到的离散序列为
为简单起见,忽略间谐波的影响,且不失一般性,以第k次谐波为例。信号x(n)经余弦偶次幂窗Wcos(n)加权截短后得到N点长
式中,λ∈[1,K],fk=λkf1,f1为基波频率。
假设以下两个条件得到满足:
1)10f1<fk<fs/2-10f1
2)|λk-λq|>BW,k,q∈[1,K],k≥2,k≠q
则负频点旁瓣的影响(式11中第2项)以及其余各次谐波对第k次谐波的泄漏影响(式11中第3项)均可忽略,此时式(11)可变为
谐波信号参数可以通过Xcos(λ)来表征,同步采样时,第k次谐波对应频点λk处的一条单一谱线;非同步采样下,能量泄漏到整个频带,由于栅栏效应,第k次谐波对应的峰值点偏离离散频点。此时,λk为非整数,假设
λk=lk+δk,δk∈[0,1](13)
式中,lk、δk分别对应实部和小数部分。
假设峰值点附近对应的幅值最大和次最大谱线分别为第lk和lk+1条谱线,其对应的幅值分别为yk1=|Xcos(lk)|、yk2=|Xcos(lk+1)|。
定义μk=δk-0.5、εk=(yk1-yk2)/(yk1+yk2),经过变量代换可得:
由于|-μk±0.5|≤1,且N一般较大,由式(6)可得:
利用最小二乘曲线拟合,可得式(15)反函数μk=f1(εk)的逼近式μk=F(εk),求出参数μk,则第k次谐波的频率修正式为
fk=λkfs/N=(lk+μk+0.5)fs/N(16)
对第lk和lk+1条谱线的幅值进行加权平均,即可得到第k次谐波的幅值修正式
当N较大时,由式(15)可得式(17)的简化表达式为Ak=N-1(yk1+yk2)v(μk),进一步采用最小二乘拟合逼近,可得其逼近式:
Ak=N-1(yk1+yk2)R(μk)(18)
结合式(12)、(15)可得到第k次谐波的相位修正式为
附图说明
图1是本实用新型所述改进傅里叶变换的谐波电能计量***结构方框图。
图2是本实用新型所述改进傅里叶变换的谐波电能计量***的控制方法流程图。
具体实施方案
以下通过实施例对本实用新型的技术方案进一步说明。
如图1所示,本实用新型所述的改进傅里叶变换的谐波电能计量***,包括三相电网数据采集模块、数据处理模块、数据管理模块以及***模块四个部分。其中,三相电网数据采集模块包括电阻分压网络、TA与电阻网络以及ADS8364芯片,实现对电网三相电压、电流信号的实时同步采样;采样数据通过并行接口送往数据处理模块中;数据处理模块采用ADSP-BF533芯片,基于余弦偶次幂窗改进傅里叶变换谐波分析方法对电压、电流波形进行频谱分析,分离出基波及各次谐波分量,完成电参量测量、谐波分析及谐波电能计量,处理结果通过异步串口发送到数据管理模块,由M30624FGPFP芯片完成数据统计、存储、通信以及功能选择等,最后通过***模块完成数据的显示输出。
如图2所示,本实用新型所述的改进傅里叶变换谐波电能计量***的控制方法流程如下:首先对含有谐波与间谐波分量的多频率信号按采样率fs进行数据采集,得到其离散序列;然后基于余弦偶次幂窗改进傅里叶变换谐波计量方法对离散序列进行加窗处理,得到谐波信号的表征公式;再运用最小二乘法对峰值谱线进行拟合,推导出信号基波与各次谐波的频率、幅值和相位计算修正式;最后通过谐波幅值、频率和初相角修正式分析计算出谐波的幅值、频率、初相角。
Claims (4)
1.一种改进傅里叶变换的谐波电能计量***,其特征在于,该***由三相电网数据采集模块、数据处理模块、数据管理模块以及***模块组成,其中,三相电网数据采集模块通过并行接口与数据处理模块相连,数据处理模块与数据管理模块通过异步串口相连、数据管理模块与***模块相连。
2.如权利要求1所述的改进傅里叶变换的谐波电能计量***,其特征在于,所述三相电网数据采集模块包括电阻分压网络、TA与电阻网络和采样芯片ADS8364。
3.根据权利要求1所述的改进傅里叶变换的谐波电能计量***,其特征在于,所述数据处理模块采用信号处理器ADSP-BF533。
4.如权利要求1所述的改进傅里叶变换的谐波电能计量***,其特征在于,所述数据管理模块采用单片机M30624FGPFP。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201320032423 CN203012015U (zh) | 2013-01-22 | 2013-01-22 | 一种改进傅里叶变换的谐波电能计量*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201320032423 CN203012015U (zh) | 2013-01-22 | 2013-01-22 | 一种改进傅里叶变换的谐波电能计量*** |
Publications (1)
Publication Number | Publication Date |
---|---|
CN203012015U true CN203012015U (zh) | 2013-06-19 |
Family
ID=48603597
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201320032423 Expired - Fee Related CN203012015U (zh) | 2013-01-22 | 2013-01-22 | 一种改进傅里叶变换的谐波电能计量*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN203012015U (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103105532A (zh) * | 2013-01-22 | 2013-05-15 | 广西电网公司电力科学研究院 | 一种改进傅里叶变换的谐波电能计量***及其控制方法 |
CN110392827A (zh) * | 2017-03-14 | 2019-10-29 | 赫姆霍兹-森德拉姆德雷斯顿-罗森多夫研究中心 | 用于表征测量对象的电阻的设备 |
-
2013
- 2013-01-22 CN CN 201320032423 patent/CN203012015U/zh not_active Expired - Fee Related
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103105532A (zh) * | 2013-01-22 | 2013-05-15 | 广西电网公司电力科学研究院 | 一种改进傅里叶变换的谐波电能计量***及其控制方法 |
CN110392827A (zh) * | 2017-03-14 | 2019-10-29 | 赫姆霍兹-森德拉姆德雷斯顿-罗森多夫研究中心 | 用于表征测量对象的电阻的设备 |
US11156574B2 (en) | 2017-03-14 | 2021-10-26 | Helmholtz-Zentrum Dresden-Rossendorf E.V. | Apparatus for characterizing the electrical resistance of a measurement object |
CN110392827B (zh) * | 2017-03-14 | 2022-05-13 | 赫姆霍兹-森德拉姆德雷斯顿-罗森多夫研究中心 | 用于表征测量对象的电阻的设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102435844B (zh) | 一种频率无关的正弦信号相量计算方法 | |
CN102890190B (zh) | 基于多小波的随机非线性负荷有功电能计量方法 | |
CN103995181B (zh) | 数字化变电站电能质量谐波分析方法 | |
CN103257271A (zh) | 一种基于stm32f107vct6的微电网谐波与间谐波检测装置及检测方法 | |
CN102323477A (zh) | 一种电网间谐波测量方法及测量仪 | |
CN102331526B (zh) | 汉宁窗函数连续频谱内插电力谐波参数获取方法 | |
CN103308766A (zh) | 一种基于凯撒自卷积窗双谱线插值fft谐波分析方法及其装置 | |
CN103454497A (zh) | 基于改进加窗离散傅立叶变换的相位差测量方法 | |
CN203133168U (zh) | 电力谐波检测仪 | |
CN102288807A (zh) | 一种测量电网电压闪变的方法 | |
CN203287435U (zh) | 一种基于stm32f107vct6的微电网谐波与间谐波检测装置 | |
CN103575984A (zh) | 基于凯塞窗双谱线插值fft的谐波分析方法 | |
CN106841778A (zh) | 基于pmu实现的次同步和超同步谐波参数的处理方法 | |
CN103926462B (zh) | 一种电力***谐波快速分析方法及运行装置 | |
CN105004939A (zh) | 一种复合电能质量扰动信号定量分析方法 | |
CN103198184A (zh) | 一种电力***低频振荡特征类噪声辨识方法 | |
CN107064744A (zh) | 一种谐波源定位方法 | |
CN103018555A (zh) | 一种高精度的电力参数软件同步采样方法 | |
CN103105529A (zh) | 一种基于参数分析的谐波电能计量***及其控制方法 | |
CN105445541A (zh) | 一种任意频率下自适应功率计算方法 | |
CN101701985A (zh) | 定频变点电网谐波检测方法及其测量仪 | |
CN106018960B (zh) | 一种基于压缩传感的同步相量测量方法 | |
CN105137181A (zh) | 基于Nuttall-Kaiser组合窗双谱线插值的谐波分析方法 | |
CN103091545A (zh) | 一种频率无关的正弦信号相量半波计算方法 | |
CN103412209B (zh) | 一种调谐滤波器的失谐度检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130619 Termination date: 20210122 |
|
CF01 | Termination of patent right due to non-payment of annual fee |