CN105629064B - 一种正弦波低通滤波提取方法 - Google Patents
一种正弦波低通滤波提取方法 Download PDFInfo
- Publication number
- CN105629064B CN105629064B CN201511020469.2A CN201511020469A CN105629064B CN 105629064 B CN105629064 B CN 105629064B CN 201511020469 A CN201511020469 A CN 201511020469A CN 105629064 B CN105629064 B CN 105629064B
- Authority
- CN
- China
- Prior art keywords
- low
- formula
- pass filtering
- output signal
- sine wave
- 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
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)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种正弦波低通滤波提取方法,包括如下步骤:步骤一,采样输入信号,确定采样周期t、提取目标角频率ω及滤波系数m;步骤二,利用采用得到的周期t、目标角频率ω及滤波系数m,建立离散状态方程,并根据建立的离散状态方程求得当前时刻的一对正交输出信号Yα(n)和Yβ(n);步骤三,根据求得的一对正交输出信号Yα(n)和Yβ(n)得到谐波的幅值和相位,完成谐波的提取。本发明的正弦波低通滤波提取方法,通过步骤一、步骤二和步骤三的设置,就可以有效的提取出谐波,且计算简单,计算量少,幅值衰减少。
Description
技术领域
本发明涉及一种正弦波谐波提取方法,更具体的说是涉及一种正弦波低通滤波提取方法。
背景技术
随着光伏、风电等新能源电力大规模入网,以及越来越多的非线性负载用电设备的使用,电网的谐波污染越来越严重。谐波使电能的生产、传输和利用的效率降低,使电气设备过热、产生振动和噪声,并使绝缘老化,使用寿命缩短,甚至发生故障或烧毁。谐波可引起电力系统局部并联谐振或串联谐振,使谐波含量放大,造成电容器等设备烧毁。谐波还会引起继电保护和自动装置误动作,使电能计量出现混乱。对于电力系统外部,谐波对通信设备和电子设备会产生严重干扰。所以,一方面光伏逆变器需要能够适应于高谐波污染的电网环境,另一方面需要使光伏逆变器发送到电网的谐波尽可能地减小。
为了实现上述要求,首先要做的就是实现对基波和谐波的提取。传统的低通滤波器对直流分量没有相移,滤波效果好,但是对高频交流分量相移大,幅值衰减严重。而采用傅里叶变换提取谐波的方法,计算复杂,计算量大,占用处理器内存。
发明内容
针对现有技术存在的不足,本发明的目的在于提供一种可以有效的提取谐波的正弦波低通滤波提取方法。
为实现上述目的,本发明提供了如下技术方案:一种正弦波低通滤波提取方法,包括如下步骤:
步骤一,采样输入信号,确定采样周期t、目标角频率ω及滤波系数m;
步骤二,利用采样得到的周期t、目标角频率ω及滤波系数m,建立离散状态方程,并根据建立的离散状态方程求得当前时刻的一对正交输出信号Yα(n)和Yβ(n);
步骤三,根据求得的一对正交输出信号Yα(n)和Yβ(n)得到谐波的幅值和相位,完成谐波的提取。
作为本发明的进一步改进,所述步骤二建立的离散状态方程如下:
式中:Yα(n)、Yβ(n)为一对正交输出信号,对应输入信号X(n),Yα(n-1)、Yβ(n-1)为上一时刻的输出,ω为目标角频率,t为离散时间常量,m为滤波系数,取值为0~1之间的实数。
作为本发明的进一步改进,所述步骤二中的离散状态方程采用旋转公式和低通滤波器公式相融合而成。
作为本发明的进一步改进,所述旋转公式如下:
式中:(α,β)为已知向量,(α′,β′)为已知向量(α,β)逆时针旋转θ角度后得到的向量,在融合时将旋转公式变形,改为:
变形后式中,ω为提取的目标信号角频率,t为离散时间常量,α(n)、β(n)为一定时刻的向量,α(n-1),β(n-1)为上一时刻的向量。
作为本发明的进一步改进,所述低通滤波器公式如下:
Y(n)=(1-m)×Y(n-1)+m×X(n)
式中:Y(n)为一定时刻输出信号,m为滤波系数,Y(n-1)为上一时刻输出信号,X(n)为输入信号。
作为本发明的进一步改进,所述旋转公式由下列公式推导而出:
将一个已知向量(α,β)逆时针旋转θ角度后,得到向量(α′,β′);将
代入
得到
即为旋转公式其中,为向量(α,β)的向量角度,L为向量长度。
作为本发明的进一步改进,所述低通滤波器公式由下列公式推导而出:
根据低通滤波器传递函数:
式中:X(s)为输入,Y(s)为输出,a为滤波带宽,s为时域常数,
利用上述低通滤波器传递函数,前向通道为积分:
Y(n)=Y(n-1)+aT(X(n)-Y(n-1))=(1-aT)Y(n-1)+aTX(n)
设m=aT,得低通滤波器离散公式
Y(n)=(1-m)×Y(n-1)+m×X(n)
式中:Y(n)为一定时刻输出信号,m为滤波系数,Y(n-1)为上一时刻输出信号,X(n)为输入信号,a为滤波带宽,T为离散周期。
本发明具有以下有益效果,通过步骤一的设置,就可以有效的确定采样周期t、提取目标角频率ω及滤波系数m,而通过步骤二的设置,就可以有效的利用采样得到的周期t、目标角频率ω及滤波系数m,建立离散状态方程,并根据建立的离散状态方程求得当前时刻的一对正交输出信号Yα(n)和Yβ(n),然后通过步骤三的设置,根据求得的一对正交输出信号Yα(n)和Yβ(n)得到谐波的幅值和相位,完成谐波的提取,这样就能够简单有效的提取出谐波信号了,期间并未采用傅里叶变换的方式,因而相比于现有技术中的谐波提取方式,其更加的简单方便,并且计算量小。
附图说明
图1为提取cos(ωT)得到的波形图;
图2为提取cos(3ωT)得到的波形图;
图3为提取cos(5ωT)得到的波形图。
具体实施方式
下面将结合附图所给出的实施例对本发明做进一步的详述。
本实施例的一种正弦波低通滤波提取方法,其特征在于:包括如下步骤:
步骤一,采样输入信号,确定采样周期t、目标角频率ω及滤波系数m;
步骤二,利用采样得到的周期t、目标角频率ω及滤波系数m,建立离散状态方程,并根据建立的离散状态方程求得当前时刻的一对正交输出信号Yα(n)和Yβ(n);
步骤三,根据求得的一对正交输出信号Yα(n)和Yβ(n)得到谐波的幅值和相位,完成谐波的提取,在使用本方法提取谐波的过程中,先进行步骤一,对输入信号进行相关参数的采样,确定采样周期t、提取目标角频率ω及滤波系数m,然后进行步骤二,利用确定采样周期t、提取目标角频率ω及滤波系数m建立起离散状态方程,接着根据离散状态方程求得当前时刻的一对正交输出信号Yα(n)和Yβ(n),在求得当前时刻的一对正交输出信号Yα(n)和Yβ(n)后,就可以有效的得到谐波的幅值和相位,如此便可以有效的提取出谐波了,十分的简单和方便,期间没有用到傅立叶变换,所以相比现有技术中的提取谐波的方法,计算量要少的多。
作为改进的一种具体实施方式,所述步骤二建立的离散状态方程如下:
式中:Yα(n)、Yβ(n)为一对正交输出信号,对应输入信号X(n),Yα(n-1)、Yβ(n-1)为上一时刻的输出,ω为目标角频率,t为离散时间常量,m为滤波系数,取值为0~1之间的实数,在本实施例中,确定系统采样周期为t=0.0001s,采样输入信号为X,设基波为50Hz,提取的目标频率为150Hz,则ω=2πf=300π,m=0.001,保存上一时刻输出,计算出当前时刻的输出:
Yα(n)=0.999×(cos(ωt)×Yα(n-1)-sin(ωt)×Yβ(n-1))+0.001×X(n)
Yβ(n)=0.999×(sin(ωt)×Yα(n-1)+cos(ωt)×Yβ(n-1))+0.001×Yβ(n-1)
其中Yα(n-1)为上一时刻的实部输出,Yβ(n-1)为上一时刻的虚部输出,计算出Yα(n)和Yβ(n),如此便可以有效的获得谐波的幅值和相位,完成谐波的提取了。
作为改进的一种具体实施方式,所述步骤二中的离散状态方程采用旋转公式和低通滤波器公式相融合而成,通过旋转公式和低通滤波器公式的融合便可以有效的得到离散状态方程,就可以使得整个提取过程更加的简单和方便,计算量大大的减小。
作为改进的一种具体实施方式,所述旋转公式如下:
式中:(α,β)为已知向量,(α′,β′)为已知向量(α,β)逆时针旋转θ角度后得到的向量,在融合时将旋转公式变形,改为:
变形后式中,ω为目标角频率,t为离散时间常量,α(n)、β(n)为一定时刻的向量,α(n-1),β(n-1)为上一时刻的向量,通过将旋转公式变形,就使得旋转公式能够更好的与低通滤波器公式融合,完成离散状态方程的建立了,有效的降低了计算量。
作为改进的一种具体实施方式,所述低通滤波器公式如下:
Y(n)=(1-m)×Y(n-1)+m×X(n)
式中:Y(n)为一定时刻输出信号,m为滤波系数,Y(n-1)为上一时刻输出信号,X(n)为输入信号,将低通滤波器公式设置成上式,就可以有效的利用到之前步骤一所采样到的滤波系数m了,并且公式整体简单,计算量少。
作为改进的一种具体实施方式,所述旋转公式由下列公式推导而出:
将一个已知向量(α,β)逆时针旋转θ角度后,得到向量(α′,β′);将
代入
得到
即为旋转公式其中,为向量(α,β)的向量角度,L为向量长度,如此便可以简单有效的推导出旋转公式,有效的降低了离散状态方程建立的计算量。
作为改进的一种具体实施方式,所述低通滤波器公式由下列公式推导而出:
根据低通滤波器传递函数:
式中:X(s)为输入,Y(s)为输出,a为滤波带宽,s为时域常数,
利用上述低通滤波器传递函数,前向通道为积分:
Y(n)=Y(n-1)+aT(X(n)-Y(n-1))=(1-aT)Y(n-1)+aTX(n)
设m=aT,得低通滤波器离散公式
Y(n)=(1-m)×Y(n-1)+m×X(n)
式中:Y(n)为一定时刻输出信号,m为滤波系数,Y(n-1)为上一时刻输出信号,X(n)为输入信号,a为滤波带宽,T为离散周期,如此便可以简单有效的推导出低通滤波公式,有效的降低了离散状态方程建立的计算量。
本实施例中主要例举了对信号X=6+3×cos(ωT)+2×cos(3ωT+1.8)+1×cos(5ωT)进行滤波,ω=100π所提取cos(ωT)、cos(3ωT)以及cos(5ωT)所得到的波形,分别对应如图1、图2和图3所示,可以明显看出本发明能够十分有效地提取所需要的交流分量,具有计算简单,占用内存少的优点。
综上所述,本发明采用步骤一的设置可以有效的采样到提取谐波所需要的周期t、目标角频率ω及滤波系数m,然后通过步骤二和步骤三的设置,就可以有效的建立离散状态方程,然后利用离散状态方程求得的一对正交输出信号Yα(n)和Yβ(n)得到谐波的幅值和相位,完成谐波的提取。
以上所述仅是本发明的优选实施方式,本发明的保护范围并不仅局限于上述实施例,凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理前提下的若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (6)
1.一种正弦波低通滤波提取方法,其特征在于:包括如下步骤:
步骤一,采样输入信号X(n),确定采样周期t、目标角频率ω及低通滤波系数m;
步骤二,利用采样得到的周期t、目标角频率ω及低通滤波系数m,建立离散状态方程,并根据建立的离散状态方程求得当前时刻的一对正交输出信号Yα(n)和Yβ(n);
步骤三,根据求得的一对正交输出信号Yα(n)和Yβ(n)得到谐波的幅值和相位,完成谐波的提取,所述步骤二建立的离散状态方程如下:
式中:Yα(n)、Yβ(n)为一对正交输出信号,对应输入信号X(n),Yα(n-1)、Yβ(n-1)为上一时刻的输出,ω为目标角频率,t为离散时间常量,m为低通滤波系数,取值为0~1之间的实数。
2.根据权利要求1所述的正弦波低通滤波提取方法,其特征在于:所述步骤二中的离散状态方程采用旋转公式和低通滤波器公式相融合而成。
3.根据权利要求2所述的正弦波低通滤波提取方法,其特征在于:所述旋转公式如下:
式中:(α,β)为已知向量,(α′,β′)为已知向量(α,β)逆时针旋转θ角度后得到的向量,在融合时将旋转公式变形,改为:
变形后式中,ω为目标角频率,t为离散时间常量,α(n)、β(n)为一定时刻的向量,α(n-1),β(n-1)为上一时刻的向量。
4.根据权利要求2或3所述的正弦波低通滤波提取方法,其特征在于:所述低通滤波器公式如下:
Y(n)=(1-m)×Y(n-1)+m×X(n)
式中:Y(n)为一定时刻输出信号,m为低通滤波系数,Y(n-1)为上一时刻输出信号,X(n)为输入信号。
5.根据权利要求3所述的正弦波低通滤波提取方法,其特征在于:所述旋转公式由下列公式推导而出:
将一个已知向量(α,β)逆时针旋转θ角度后,得到向量(α′,β′);将
代入
得到
即为旋转公式其中,
为向量(α,β)的向量角度,L为向量长度。
6.根据权利要求4所述的正弦波低通滤波提取方法,其特征在于:所述低通滤波器公式由下列公式推导而出:
根据低通滤波器传递函数:
式中:X(s)为输入,Y(s)为输出,a为滤波带宽,s为时域常数,
利用上述低通滤波器传递函数,前向通道为积分:
Y(n)=Y(n-1)+aT(X(n)-Y(n-1))=(1-aT)Y(n-1)+aTX(n)
设m=aT,得低通滤波器离散公式
Y(n)=(1一m)×Y(n-1)+m×X(n)
式中:Y(n)为一定时刻输出信号,m为低通滤波系数,Y(n-1)为上一时刻输出信号,X(n)为输入信号,a为滤波带宽,T为离散周期。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511020469.2A CN105629064B (zh) | 2015-12-30 | 2015-12-30 | 一种正弦波低通滤波提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511020469.2A CN105629064B (zh) | 2015-12-30 | 2015-12-30 | 一种正弦波低通滤波提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105629064A CN105629064A (zh) | 2016-06-01 |
CN105629064B true CN105629064B (zh) | 2018-10-02 |
Family
ID=56044187
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201511020469.2A Active CN105629064B (zh) | 2015-12-30 | 2015-12-30 | 一种正弦波低通滤波提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105629064B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109030938B (zh) * | 2017-06-08 | 2021-05-11 | 许继集团有限公司 | 一种基于正弦滤波的抗谐波测频方法和装置 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0337573A (ja) * | 1989-07-04 | 1991-02-18 | Fukuoka Denki Keiki Kk | ディジタル高調波測定器 |
CN103219745B (zh) * | 2013-04-19 | 2014-12-03 | 浙江埃菲生能源科技有限公司 | 一种基于正交正弦波提取器的并网逆变控制算法 |
CN103618310B (zh) * | 2013-12-05 | 2016-09-21 | 国家电网公司 | 一种大容量统一电能质量控制器及其控制方法 |
CN103825302B (zh) * | 2014-03-12 | 2016-02-17 | 浙江埃菲生能源科技有限公司 | 一种改进的离散正交信号发生器 |
CN104319818B (zh) * | 2014-11-18 | 2016-06-08 | 浙江埃菲生能源科技有限公司 | 一种基于电网电压正序提取器的并网与独立双模式切换方法 |
CN104833851B (zh) * | 2015-04-24 | 2017-11-17 | 合肥工业大学 | 基于分布式相关卡尔曼滤波的电力系统谐波估计方法 |
-
2015
- 2015-12-30 CN CN201511020469.2A patent/CN105629064B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105629064A (zh) | 2016-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106560718B (zh) | 一种接地故障电流检测保护方法及装置 | |
Zhang et al. | Realization of the transient-based boundary protection for HVDC transmission lines based on high frequency energy criteria | |
Gao et al. | Waveform-based method for fast and accurate identification of subsynchronous resonance events | |
CN105449718B (zh) | 基于改进型串联信号延迟对消算法的并网同步锁相方法 | |
WO2019109969A1 (zh) | 一种电力系统宽频带多振荡模式分量的辨识方法 | |
CN101900751B (zh) | 能量回馈单元中实时电压矢量的检测方法 | |
HE et al. | A new method of designing double-tuned filter | |
CN103197130A (zh) | 一种电网电压过零点检测方法及装置 | |
CN104833853A (zh) | 一种频率自适应的滑窗dft谐波检测方法 | |
CN105629064B (zh) | 一种正弦波低通滤波提取方法 | |
CN107991531A (zh) | 基于变分模态分解和去除噪声的电网谐波检测方法 | |
CN105137179A (zh) | 一种基于级联延迟信号消除法的指定次谐波检测方法 | |
CN108490254A (zh) | 一种基于多算法融合的宽带多频信号测量方法 | |
CN106645948A (zh) | 一种基于小波变换的自适应电网信号基波检测方法 | |
CN110460104A (zh) | 一种光伏逆变器等效正、负序故障阻抗推导方法 | |
CN105656044A (zh) | 一种基于快速傅立叶变换的抑制电流谐波的方法 | |
CN102759659B (zh) | 一种电力系统中电气信号谐波瞬时值的提取方法 | |
CN106372369B (zh) | 用于逆变器输出直流分量、谐波分量检测的瞬时计算方法 | |
CN105891632B (zh) | 一种基于二阶低通陷波器的同步信号检测方法 | |
CN108120863A (zh) | 一种新型fbd谐波电流检测方法 | |
CN108957130B (zh) | 一种基于瞬时无功功率理论的电力系统谐波检测方法 | |
CN114301175A (zh) | 基于注入信号的配电台区户变关系识别的方法及装置 | |
CN102621373A (zh) | 交流电源均流输入系统的信号采集装置及处理方法 | |
CN106771597B (zh) | 一种省略Park变换的无锁相环式SAI谐波电流检测方法 | |
Kim et al. | A novel islanding detection method using goertzel algorithm in grid-connected system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |