CN108362940B - 一种半周波畸变信号的动态谐波频率提取方法 - Google Patents
一种半周波畸变信号的动态谐波频率提取方法 Download PDFInfo
- Publication number
- CN108362940B CN108362940B CN201810206902.9A CN201810206902A CN108362940B CN 108362940 B CN108362940 B CN 108362940B CN 201810206902 A CN201810206902 A CN 201810206902A CN 108362940 B CN108362940 B CN 108362940B
- Authority
- CN
- China
- Prior art keywords
- signal
- harmonic
- cycle
- input signal
- sampling
- 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
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/02—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
- Compression, Expansion, Code Conversion, And Decoders (AREA)
Abstract
本发明公开了一种半周波畸变信号的动态谐波频率提取方法,它包括:选择采样频率和采样长度,对输入信号进行等间隔采样,得到信号离散采样序列,并由系统自动搜寻信号最大值与最小值;设置m个水平阈值Yi;利用输入信号过水平阈值与过零点的时间差判断离散采样序列的半周波内是否含有谐波;若输入信号存在谐波,则将信号进行周期延拓,延拓生成为不少于2个的完整周期信号;对延拓后生成的信号进行谐波分析,提取出输入信号相应的正半或负半周波所含的谐波频率、幅值、和相位;解决了实际电网中的谐波往往具有持续时间短、发生随机性强等非平稳特点;传统的谐波分析一般采用傅立叶变换无法准确区分较短时间内的频谱成分等技术问题。
Description
技术领域
本发明属于电力技术领域,尤其涉及一种半周波畸变信号的动态谐波频率提取方法。
背景技术
随着电力电子技术的日益发展,非线性负荷的大量使用导致电力系统中电压电流波形发生畸变,谐波问题变得尤为突出。在现有的谐波分析方法中,通常谐波都被认为是稳态信号,而在实际情况下,由于非线性负载的连接,导致谐波信号往往呈动态出现。在IR46电能表国际建议标准中,其对谐波影响下的电能计量误差分析提出了新要求,在谐波对电能计量误差的影响研究中,均将谐波视为稳态信号,而实际电网中的谐波往往具有持续时间短、发生随机性强等非平稳特点。传统的谐波分析一般采用傅立叶变换,但傅立叶变换无法准确区分较短时间内的频谱成分。
发明内容
本发明要解决的技术问题:提供一种半周波畸变信号的动态谐波频率提取方法,以解决现有技术在谐波对电能计量误差的影响研究中,均将谐波视为稳态信号,而实际电网中的谐波往往具有持续时间短、发生随机性强等非平稳特点。传统的谐波分析一般采用傅立叶变换,但傅立叶变换无法准确区分较短时间内的频谱成分等技术问题。
本发明技术方案:
一种半周波畸变信号的动态谐波频率提取方法,它包括:
步骤一:选择采样频率fs和采样长度N,对输入信号进行等间隔采样,得到信号离散采样序列x(tn),n=0,1,2,…,N-1,并由系统自动搜寻信号最大值max与最小值min;
步骤二:设置m个水平阈值Yi,m为大于等于2的整数,i=1,2,…,m,阈值取值区间为0.8max<Yi<max,或min<Yi<0.8min;
步骤三:利用输入信号过水平阈值与过零点的时间差判断离散采样序列的半周波内是否含有谐波;
步骤四:若输入信号的正半或负半周波存在谐波,将相应的正半或负半周波信号进行周期延拓,延拓生成为不少于2个的完整周期信号;
步骤五:对延拓后生成的信号进行谐波分析,提取出输入信号相应的正半或负半周波所含的谐波频率、幅值、和相位。
步骤三中,判断离散采样序列的半周波内是否含有谐波的方法为:根据阈值的取值区间,确定设置的m个水平阈值Yi与输入信号曲线相交且处于输入信号曲线升降序相同的两个交点,如果某个阈值的2个交点均与信号采样点重合,则直接计算时间差值△ti,如果某个阈值的2个交点不与信号采样点重合,则将该阈值与采样信号上述交点附近的信号采样点代入Lagrange插值公式,计算该阈值相应的时间差值△ti;求出输入信号过零点时处于输入信号曲线升降序相同的两点之间的时间差,并记为△t0,令σi=|△t0–△ti|,v=|σi/σj|式中i=1,2,…,m;j=1,2,…,m;i≠j;若v>10或v<0.1,则判断输入信号的正半或负半周波存在谐波,否则判断输入信号的正半或负半周波不存在谐波。
步骤四中仅对步骤三中所判断的含有谐波的正半或负半周波信号进行周期延拓。
本发明有益效果:
本发明对在半周波出现的谐波信号进行快速判别,然后仅将含谐波的半周波信号周期延拓为整周期信号,再采用FFT算法对其进行谐波分析,便可提取出较短时间内出现的谐波信号的频率、幅值与相位参数;本发明具有简单易懂、谐波检测速度更快、更准确等特点,采用输入信号过水平阈值与过零点的时间差的差值比较法检测谐波,然后对相应含有谐波的半周波信号进行周期延拓,最后对延拓生成的信号进行谐波分析,完成半周波畸变信号的动态谐波频率、幅值与相位参数的提取,该方法检测谐波速度快,计算简单,易于微机实现等优点。解决了现有技术在谐波对电能计量误差的影响研究中,均将谐波视为稳态信号,而实际电网中的谐波往往具有持续时间短、发生随机性强等非平稳特点。传统的谐波分析一般采用傅立叶变换,但傅立叶变换无法准确区分较短时间内的频谱成分等技术问题。
附图说明
图1为本发明流程示意图。
具体实施方式
本发明实现半周波畸变信号的动态谐波频率快速提取方法的流程图如图1所示,其具体步骤包括:
第一步,选取合适的采样频率fs和采样点数N对输入信号进行等间隔取样,并由系统搜寻采样信号的最大值max与最小值min;
第二步,设置3个水平阈值,阈值取值区间为0.8max<Yi<max,或min<Yi<0.8min;
第三步,根据阈值取值区间,选取三个水平阈值Y1、Y2、Y3,将水平阈值Y1、Y2、Y3分别与输入信号曲线相交,找出处于信号曲线升降序相同的两点,如果某阈值的两个交点均与采样点重合,则可以直接计算出时间差,如果阈值的两个交点与信号采样点均不重合,则将该阈值与两交点附近的两个采样点代入Lagrange插值公式,计算相应的时间差,将计算出的时间差分别记为△t1、△t2、△t3。计算出输入信号曲线过零点时处于信号升降序相同的两个零点之间的时间差,记为△t0,检测采样信号半周波是否含有谐波的依据为:
令σ1=|△t0–△t1|,σ2=|△t0–△t2|,σ3=|△t0–△t3|,v1=|σ1/σ2|,v2=|σ1/σ3|,v3=|σ2/σ3|。
若v1、v2、v3中存在vi>10或vi<0.1,则判断输入信号正半或负半周波含有谐波;
第四步,将输入信号的正半或负半周波延拓生成8个整周期信号;
第五步,对延拓生成的整周期信号采用加窗插值FFT改进算法进行谐波分析,计算出谐波频率、幅值与相位。
至此,完成对半周波畸变信号的动态谐波频率的快速准确提取。
Claims (3)
1.一种半周波畸变信号的动态谐波频率提取方法,它包括:
步骤一:选择采样频率fs和采样长度N,对输入信号进行等间隔采样,得到信号离散采样序列x(tn),n=0,1,2,…,N-1,并由系统自动搜寻信号最大值max与最小值min;
步骤二:设置m个水平阈值Yi,m为大于等于2的整数,i=1,2,…,m,阈值取值区间为0.8max<Yi<max,或min<Yi<0.8min;
步骤三:利用输入信号过水平阈值与过零点的时间差判断离散采样序列的半周波内是否含有谐波;
步骤四:若输入信号的正半或负半周波存在谐波,将相应的正半或负半周波信号进行周期延拓,延拓生成为不少于2个的完整周期信号;
步骤五:对延拓后生成的信号进行谐波分析,提取出输入信号相应的正半或负半周波所含的谐波频率、幅值和相位。
2.根据权利要求1所述的一种半周波畸变信号的动态谐波频率提取方法,其特征在于:步骤三中,判断离散采样序列的半周波内是否含有谐波的方法为:根据阈值的取值区间,确定设置的m个水平阈值Yi与输入信号曲线相交且处于输入信号曲线升降序相同的两个交点,如果某个阈值的2个交点均与信号采样点重合,则直接计算时间差值△ti,如果某个阈值的2个交点不与信号采样点重合,则将该阈值与采样信号上述交点附近的信号采样点代入Lagrange插值公式,计算该阈值相应的时间差值△ti;求出输入信号过零点时处于输入信号曲线升降序相同的两点之间的时间差,并记为△t0,令σi=|△t0–△ti|,v=|σi/σj|式中i=1,2,…,m;j=1,2,…,m;i≠j;若v>10或v<0.1,则判断输入信号的正半或负半周波存在谐波,否则判断输入信号的正半或负半周波不存在谐波。
3.根据权利要求1所述的一种半周波畸变信号的动态谐波频率提取方法,其特征在于:步骤四中仅对步骤三中所判断的含有谐波的正半或负半周波信号进行周期延拓。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810206902.9A CN108362940B (zh) | 2018-03-14 | 2018-03-14 | 一种半周波畸变信号的动态谐波频率提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810206902.9A CN108362940B (zh) | 2018-03-14 | 2018-03-14 | 一种半周波畸变信号的动态谐波频率提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108362940A CN108362940A (zh) | 2018-08-03 |
CN108362940B true CN108362940B (zh) | 2020-07-14 |
Family
ID=63000457
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810206902.9A Active CN108362940B (zh) | 2018-03-14 | 2018-03-14 | 一种半周波畸变信号的动态谐波频率提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108362940B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109283383B (zh) * | 2018-10-19 | 2020-08-28 | 深圳市计量质量检测研究院 | 矩形调制电压的电压波动值测量方法与装置、存储介质 |
CN109613324B (zh) * | 2018-11-07 | 2021-10-08 | 全球能源互联网研究院有限公司 | 一种谐波放大的检测方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101334431A (zh) * | 2007-12-20 | 2008-12-31 | 复旦大学 | 电网谐波的频谱插值校正分析方法 |
CN102095929A (zh) * | 2010-12-17 | 2011-06-15 | 南京邮电大学 | 一种快速测量交流电信号频率的方法 |
CN103592513A (zh) * | 2013-11-20 | 2014-02-19 | 北京富根智能电表有限公司 | 电力信号谐波分析方法和装置 |
CN106018958A (zh) * | 2016-05-11 | 2016-10-12 | 西安西热节能技术有限公司 | 孤网变频系统中二次侧电流电压频率跟踪方法 |
-
2018
- 2018-03-14 CN CN201810206902.9A patent/CN108362940B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101334431A (zh) * | 2007-12-20 | 2008-12-31 | 复旦大学 | 电网谐波的频谱插值校正分析方法 |
CN102095929A (zh) * | 2010-12-17 | 2011-06-15 | 南京邮电大学 | 一种快速测量交流电信号频率的方法 |
CN103592513A (zh) * | 2013-11-20 | 2014-02-19 | 北京富根智能电表有限公司 | 电力信号谐波分析方法和装置 |
CN106018958A (zh) * | 2016-05-11 | 2016-10-12 | 西安西热节能技术有限公司 | 孤网变频系统中二次侧电流电压频率跟踪方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108362940A (zh) | 2018-08-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103257271A (zh) | 一种基于stm32f107vct6的微电网谐波与间谐波检测装置及检测方法 | |
CN107543962B (zh) | 主导间谐波频谱分布的计算方法 | |
CN103487652A (zh) | 一种频率自适应实时分次谐波检测方法 | |
CN102353500B (zh) | 一种用于动平衡测量的不平衡信号提取方法 | |
CN106501602B (zh) | 一种基于滑窗频谱分离的基波参数测量方法 | |
CN103018555B (zh) | 一种高精度的电力参数软件同步采样方法 | |
CN103941088A (zh) | 一种基于三相信号的电力系统频率快速测量方法 | |
CN203287435U (zh) | 一种基于stm32f107vct6的微电网谐波与间谐波检测装置 | |
CN108362940B (zh) | 一种半周波畸变信号的动态谐波频率提取方法 | |
WO2023071672A1 (zh) | 一种基于非线性鲁棒估计的电网信息快速精准感知方法 | |
CN108333426A (zh) | 基于傅氏算法的电力系统频率测量方法 | |
CN110687399A (zh) | 一种配电网故障指示器波形故障开始时刻判断方法 | |
CN108414838B (zh) | 一种逆变器并联系统线路阻抗测量方法 | |
CN103543331B (zh) | 一种计算电信号谐波和间谐波的方法 | |
CN203104409U (zh) | 一种市电精确检测跟踪数字锁相的系统 | |
CN107703358B (zh) | 一种基于改进二阶广义积分器的锁相算法 | |
CN109581045A (zh) | 一种满足iec标准框架的间谐波功率计量方法 | |
CN108090270A (zh) | 一种基于形态学滤波和盲源分离的暂态振荡参数识别方法 | |
CN107271772A (zh) | 一种高精度且抗噪声干扰的电网频率快速检测方法 | |
CN103560509B (zh) | 一种基于小波分析的电压下陷检测装置及该装置的控制方法 | |
CN110007166B (zh) | 一种变频器能效的快速测量方法 | |
CN108445295B (zh) | 一种基于Goertzel算法的准被动电网阻抗识别系统及方法 | |
CN108445288B (zh) | 一种适合于半周波畸变情况的电能计量方法 | |
CN113030570A (zh) | 一种谐波电能检测方法和自适应滤波器 | |
Men et al. | A new time domain filtering method for calculating the RMS value of vibration signals |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
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 |