专利转让平台_买专利_卖专利_中国高校专利技术交易-买卖发明专利上知查网

全部分类
全部分类
巴特沃斯数字滤波器及利用其实现频率自适应的滤波方法

巴特沃斯数字滤波器及利用其实现频率自适应的滤波方法

IPC分类号 : H03H17/00,

申请号
CN201510375214.1
可选规格
  • 专利类型: 发明专利
  • 法律状态: 有权
  • 申请日: 2015-06-30
  • 公开号: 105048995A
  • 公开日: 2015-11-11
  • 主分类号: H03H17/00
  • 专利权人: 西安理工大学

专利摘要

本发明公开的巴特沃斯数字滤波器,由相互串联的两个二阶带通滤波器组成,本发明公开的利用巴特沃斯数字滤波器实现频率自适应的滤波方法,其采用两个双二阶巴特沃斯带通滤波器,电网电压信号通过第一个双二阶巴特沃斯带通滤波器得到电网电压基波信号,再通过三点计算方法和均值滤波得到该信号的实时频率,该实时频率作为第二个双二阶巴特沃斯带通滤波器的中心频率,使该滤波器能够实时跟踪电网频率。本发明的巴特沃斯数字滤波器及利用其实现频率自适应的滤波方法能够实时跟踪电网的相位和频率,并能有效消除因电网频率变化而导致的相位和频率检测误差,从而使滤波效果最佳,且最终输出的电网电压基波信号幅值无衰减、相位无偏差。

权利要求

1.巴特沃斯数字滤波器,其特征在于,由相互串联的两个二阶带通滤波器组成,所述两个二阶带通滤波器的传递函数分别为:

H1(s)=(K1ω0C/Q)ss2+(0/E)s+D2ω02---(2),

H2(s)=(K2ω0C/Q)ss2+(ω0/DE)s+ω02/D2---(3),

在公式(2)和(3)中,

E=1BC+4Q2+(C+4Q2)2-(2BQ)22---(4),

D=12[BEQ+(BEQ)2-4]---(5),

其中,Q为品质因数,ω0为滤波器的角频率,B和C为低通系数,K1和K2为对应两个二阶带通滤波器的增益。

2.如权利要求1所述的巴特沃斯数字滤波器,其特征在于,所述Q=5,K1=K2=1,B=1.4142,C=1,中心频率f0=50Hz,则所述两个二阶带通滤波器的传递函数为:

H1(s)=62.82ss2+47.57s+113720---(7),

H2(S)=62.82ss2+41.29s+85660---(8).

3.利用如权利要求1或2所述的巴特沃斯数字滤波器实现频率自适应的滤波方法,其特征在于,包括如下步骤:

步骤一,将电网电压Uin经过零阶保持器离散化处理后依次通过前置的巴特沃斯数字滤波器中串联的两个二阶带通滤波器,得到电网电压的基波信号Uf0=50Hz,此信号幅值有衰减,相位有偏差,并通过采样得到Uf0=50Hz最近的三个采样点的瞬时值,记为U1,U2,U3

步骤二,利用U1,U2,U3通过三点计算方法得到所述电网电压基波信号Uf0=50Hz的频率f,幅值U和相位角α,三点计算方法的公式如下:

cosωT=U1+U32U2---(9),

tanα=U3sinωTU3cosωT-U2---(10),

U=U22+U32-2U1U3cosωT2sinωT---(11),

其中,ω是电网电压基波信号的角频率,T为三点计算方法的采样时间;

步骤三,将所述电网电压基波信号Uf0=50Hz的频率f经过均值滤波后得到所述电网电压基波信号Uf0=50Hz的实时频率fav

步骤四,将所述电网电压基波信号Uf0=50Hz的实时频率fav作为后置的巴特沃斯数字滤波器的中心频率,将所述步骤一中的电网电压Uin经过零阶保持器离散化处理后的信号和将所述电网电压基波信号Uf0=50Hz的实时频率fav同时通过后置的巴特沃斯数字滤波器中串联的两个二阶带通滤波器,得到幅值无衰减,相位无偏差的电网电压基波信号,实现对电网电压Uin频率自适应的滤波处理。

4.如权利要求3所述的利用巴特沃斯数字滤波器实现频率自适应的滤波方法,其特征在于,所述前置的巴特沃斯数字滤波器和后置的巴特沃斯数字滤波器的品质因数均为5,所述前置的巴特沃斯数字滤波器和后置的巴特沃斯数字滤波器中二阶带通滤波器的函数Z变换公式为:

Z[bs1s+bs0as2s2+as1s+as0·1-e-Tsss]=bz1z+bz0az2z2-az1z+az0---(12),

公式(12)中,系数bz1、bz0、az2、az1、az0为:

bz1=aec1Ts+bec2Ts+c---(13),

bz0=aec2Ts+bec1Ts+ce(c1+c2)Ts---(14),

az2=1(15),

az1=ec1Ts+ec2Ts---(16),

az0=e(c1+c2)T2---(17),

公式(13)-(17)中的系数a、b、c、c1、c2为:

a=bs0+bs1cs1cs1(cs1-cs2)---(18),

b=-bs0-bs1cs2cs2(cs1-cs2)---(19),

c=bs0cs1cs2---(20),

c1=-as1+as12-4as0as22as2---(21),

c2=-as1-as12-4as0as22as2---(22),

其中,Ts为采样周期。

说明书

技术领域

本发明属于电能质量监测与控制技术领域,具体涉及一种巴特沃斯数字滤波器,还涉及一种利用巴特沃斯数字滤波器实现频率自适应的滤波方法。

背景技术

我们正处于数字化时代,数字信号处理技术受到了人们的关注,其理论及算法随着计算机技术和微电子技术的发展得到了飞速地发展,并被广泛应用于语音和图象处理、数字通信、谱分析、模式识别和自动控制等领域。在数字信号处理中,数字滤波技术占有重要的地位,几乎出现在所有的数字信号处理系统中。相对于模拟滤波器来说,数字滤波器稳定、没有漂移、精度高、参数容易控制等,这些优势决定了数字滤波器的应用越来越广泛。

随着电力电子技术的飞速发展,越来越多的电力电子装置被广泛应用到各个领域。近年来,配电网中整流器、变频调速装置、电弧炉等非线性负荷不断增加,这些负荷的非线性、冲击性和不平衡的用电特性,使电网中频率波动、暂态冲击、无功功率、谐波污染及三相不平衡问题日趋严重。这些问题的出现导致装置的频率和相位检测存在稳态误差,从而使得装置的性能变差。因此,频率自适应滤波器设计是近年来的一个热点研究方向。

发明内容

本发明的目的在于提供一种巴特沃斯数字滤波器,其在中心频率处增益为0dB,相移几乎为零,具有十分理想的幅频特性和相频特性。

本发的目的还在于提供一种利用该巴特沃斯数字滤波器实现频率自适应的滤波方法,可以在频率变化、谐波未知的畸变电网中,准确快速的检测电网电压基波信号的频率和相位,可应用于电压锁相环和特定次谐波检测中。

本发明所采用的一种技术方案是:巴特沃斯数字滤波器,由相互串联的两个二阶带通滤波器组成,两个二阶带通滤波器的传递函数分别为:

H1(s)=(K1ω0C/Q)ss2+(0/E)s+D2ω02---(2),]]>

H2(s)=(K2ω0C/Q)ss2+(ω0/DE)s+ω02/D2---(3),]]>

公式(2)和(3)中,

E=1BC+4Q2+(C+4Q2)2-(2BQ)22---(4),]]>

D=12[BEQ+(BEQ)2-4]---(5),]]>

其中,Q为品质因数,ω0为滤波器的角频率,B和C为低通系数,K1和K2为对应两个二阶带通滤波器的增益。

本发明的特点还在于,

Q=5,K1=K2=1,B=1.4142,C=1,中心频率f0=50Hz,则两个二阶带通滤波器的传递函数为:

H1(s)=62.82ss2+47.57s+113720---(7),]]>

H2(s)=62.82ss2+41.29s+85660---(8)]]>

本发明所采用的另一种技术方案为:利用巴特沃斯数字滤波器实现频率自适应的滤波方法,包括如下步骤:

步骤一,将电网电压Uin经过零阶保持器离散化处理后依次通过前置的巴特沃斯数字滤波器中串联的两个二阶带通滤波器,得到电网电压的基波信号Uf0=50Hz,此信号幅值有衰减,相位有偏差,并通过采样得到Uf0=50Hz最近的三个采样点的瞬时值,记为U1,U2,U3

步骤二,利用U1,U2,U3通过三点计算方法得到所述电网电压基波信号Uf0=50Hz的频率f,幅值U和相位角α,三点计算方法的公式如下:

cosωT=U1+U32U2---(9),]]>

tanα=U3sinωTU3cosωT-U2---(10),]]>

U=U22+U32-2U1U3cosωT2sinωT---(11),]]>

其中,ω是电网电压基波信号的角频率,T为三点计算方法的采样时间;

步骤三,将电网电压基波信号Uf0=50Hz的频率f经过均值滤波后得到电网电压基波信号Uf0=50Hz的实时频率fav

步骤四,将电网电压基波信号Uf0=50Hz的实时频率fav作为后置的巴特沃斯数字滤波器的中心频率,将步骤一中的电网电压Uin经过零阶保持器离散化处理后的信号和将电网电压基波信号Uf0=50Hz的实时频率fav同时通过后置的巴特沃斯数字滤波器中串联的两个二阶带通滤波器,得到幅值无衰减,相位无偏差的电网电压基波信号,实现对电网电压Uin频率自适应的滤波处理。

本发明的特点还在于,

前置的巴特沃斯数字滤波器和后置的巴特沃斯数字滤波器的品质因数均为5,前置的巴特沃斯数字滤波器和后置的巴特沃斯数字滤波器中二阶带通滤波器的函数Z变换公式为:

Z[bs1s+bs0as2s2+as1s+as0·1-e-TssS]=bz1z+bz0az2z2-az1z+az0---(12),]]>

公式(12)中,系数bz1、bz0、az2、az1、az0为:

bz1=aec1Ts+bec2Ts+c---(13),]]>

bz0=aec2Ts+bec1Ts+ce(c1+c2)Ts---(14),]]>

az2=1(15),

az1=ec1Ts+ec2Ts---(16),]]>

az0=e(c1+c2)Ts---(17),]]>

公式(13)-(17)中的系数a、b、c、c1、c2为:

a=bs0+bs1cs1cs1(cs1-cs2)---(18),]]>

b=-bs0-bs1cs2cs2(cs1-cs2)---(19),]]>

c=bs0cs1cs2---(20),]]>

c1=-as1+as12-4as0as22as2---(21),]]>

c2=-as1-as12-4as0as22as2---(22),]]>

其中,Ts为采样周期。

本发明的有益效果是:本发明的巴特沃斯数字滤波器及利用其实现频率自适应的滤波方法针对电网频率波动导致相位和频率检测存在误差从而滤波效果变差的缺点,采用两个双二阶的巴特沃斯带通滤波器,电网电压信号通过第一个双二阶巴特沃斯带通滤波器得到幅值有衰减,相位有偏差的电网电压基波信号,再通过三点计算方法和均值滤波得到该信号的实时频率,该实时频率作为第二个双二阶巴特沃斯带通滤波器的中心频率,得到幅值无衰减,相位无偏差的电网电压基波信号。使该滤波器能够实时跟踪电网的相位和频率,并能有效消除因电网频率变化而导致的相位和频率检测误差,从而使滤波效果最佳,且最终输出的电网电压基波信号幅值无衰减、相位无偏差,尤其在频率变化、谐波未知的畸变电网中,能够实时精确检测电网电压基波信号的频率和相位,达到频率自适应的目的。并具有以下有益效果:

1、参数的整定简洁方便,特别是在频率变化、谐波未知的畸变电网中,可准确快速的检测电网电压基波信号的频率和相位;

2、具有稳态精度高,动态响应较快的优点;

3、操作灵活,可应用于电压锁相环和特定次谐波检测;

4、可适用于单相或三相电网系统。

附图说明

图1是本发明的利用巴特沃斯数字滤波器实现频率自适应的滤波方法的结构框图;

图2是二阶带通滤波器在不同品质因数时的幅频和相频特性曲线图;

图3是本发明的巴特沃斯数字滤波器的幅频和相频特性曲线图;

图4是电压锁相环的结构框图;

图5是特定次谐波检测框图;

图6a是本发明的利用巴特沃斯数字滤波器实现频率自适应的滤波方法中电网频率由50Hz跃变到45Hz的仿真波形图;

图6b是本发明的利用巴特沃斯数字滤波器实现频率自适应的滤波方法中电网频率由50Hz跃变到55Hz的仿真波形图;

图7是本发明的利用巴特沃斯数字滤波器实现频率自适应的滤波方法中电网电压畸变仿真波形图;

图8a是本发明的利用巴特沃斯数字滤波器实现频率自适应的滤波方法中电网频率由50Hz跃变到45Hz的实验波形图;

图8b是本发明的利用巴特沃斯数字滤波器实现频率自适应的滤波方法中电网频率由50Hz跃变到55Hz的实验波形图;

图9a是本发明的利用巴特沃斯数字滤波器实现频率自适应的滤波方法中电网频率由50Hz跃变到45Hz且电压畸变的实验波形图;

图9b是本发明的利用巴特沃斯数字滤波器实现频率自适应的滤波方法中电网频率由50Hz跃变到55Hz且电压畸变的实验波形图;

图10a是本发明的利用巴特沃斯数字滤波器实现频率自适应的滤波方法中相位检测实验波形图;

图10b是本发明的利用巴特沃斯数字滤波器实现频率自适应的滤波方法中频率检测实验波形图。

具体实施方式

下面结合附图和具体实施方式对本发明进行详细说明。

带通滤波器采用以中心频率f0=ω0/2πHz为几何中心,且通带宽度为BW的频带信号。通带宽度BW=ωUL,其中ωL和ωU分别为低端和高端的截止频率;Q为滤波器的品质因素,Q=ω0/BW,它用来衡量滤波器的选择性。当分别选取品质因数Q为1、2、3、4、5时,二阶带通滤波器的幅频和相频特性曲线如图2所示:Q值越高,通带宽度越窄,滤波效果越好,但响应速度越慢,Q值越低,通带宽度越宽,响应速度越快,但滤波效果越差,综合考虑,本发明的巴特沃斯数字滤波器选择Q=5。

在计算机上实现巴特沃斯带通滤波器时,有限字长可能引起滤波器数值不稳,从而降低滤波器的带宽,为了避免上述问题,本发明将选择的四阶巴特沃斯数字滤波器分解成相互串联的两个二阶带通滤波器,即双二阶巴特沃斯带通滤波器,四阶巴特沃斯数字滤波器的传递函数为:

H(s)=(KCω02/Q2)s2s4+(0/Q)s3+(2+C/Q2)ω02s2+(03/Q)s+ω04---(1),]]>

其中B和C是低通系数,可查阅相关手册获得,K是电路增益。

公式(1)的传递函数可以分解为两个二阶带通滤波器的传递函数,分别为:

H1(s)=(K1ω0C/Q)ss2+(0/E)s+D2ω02---(2),]]>

H2(s)=(K2ω0C/Q)ss2+(ω0/DE)s+ω02/D2---(3),]]>

在公式(2)和(3)中,

E=1BC+4Q2+(C+4Q2)2-(2BQ)22---(4),]]>

D=12[BEQ+(BEQ)2-4]---(5),]]>

其中,Q为品质因数,ω0为滤波器的角频率,B和C为低通系数,K1和K2为对应两个二阶带通滤波器的增益,且K1K2=K。

带通滤波器,不但要能滤除常见的高次谐波,还要能滤除某些特殊负载产生的次谐波(低于基波频率的分量)。所以,品质因素Q的选取,要足够大,以保证对高次谐波和次谐波的衰减;但是如果Q太大,响应速度非常慢。

因此,本发明的巴特沃斯数字滤波器取值Q=5,K1=K2=1,查阅滤波器设计手册取B=1.4142,C=1,中心频率f0=50Hz,将以上参数分别代入公式(1)、(2)和(3)中,得到四阶巴特沃斯数字滤波器的传递函数为:

H(s)=3948s2s4+88.86s3+2.013×105s2+8.77×106s+9.741×109---(6),]]>

两个二阶带通滤波器的传递函数为:

H1(s)=<

巴特沃斯数字滤波器及利用其实现频率自适应的滤波方法专利购买费用说明

专利买卖交易资料

Q:办理专利转让的流程及所需资料

A:专利权人变更需要办理著录项目变更手续,有代理机构的,变更手续应当由代理机构办理。

1:专利变更应当使用专利局统一制作的“著录项目变更申报书”提出。

2:按规定缴纳著录项目变更手续费。

3:同时提交相关证明文件原件。

4:专利权转移的,变更后的专利权人委托新专利代理机构的,应当提交变更后的全体专利申请人签字或者盖章的委托书。

Q:专利著录项目变更费用如何缴交

A:(1)直接到国家知识产权局受理大厅收费窗口缴纳,(2)通过代办处缴纳,(3)通过邮局或者银行汇款,更多缴纳方式

Q:专利转让变更,多久能出结果

A:著录项目变更请求书递交后,一般1-2个月左右就会收到通知,国家知识产权局会下达《转让手续合格通知书》。

动态评分

0.0

没有评分数据
没有评价数据
×

打开微信,点击底部的“发现”

使用“扫一扫”即可将网页分享至朋友圈

×
复制
用户中心
我的足迹
我的收藏

您的购物车还是空的,您可以

  • 微信公众号

    微信公众号
在线留言
返回顶部