一种幅度非衰减均衡相位的分数阶椭圆滤波器设计方法技术

本发明专利技术提供一种幅度非衰减均衡相位的分数阶椭圆滤波器设计方法，具体包括：首先算法初始化，设计了一种幅度为衰减的全通滤波器，使得椭圆滤波器的阶数由只能取整数转变为任意正实数，然后引入幅度未衰减逆向滤波器，估计误差信号幅度未衰减，而相位呈现负相位，有效补偿分数阶模拟滤波器的相位；因此，本发明专利技术提出的幅度非衰减均衡相位的分数阶椭圆滤波器设计方法，滤波器的阶数调节变得更加灵活，由于逆向滤波器的引入，有效补偿了滤波器的相位延迟，提高了椭圆滤波器的线性度，同时滤波器的幅度未衰减，因此椭圆滤波器陡峭的过渡带特性并未改变，滤波器的综合性能得以提高，具有较好的工程实用性。有较好的工程实用性。有较好的工程实用性。

【技术实现步骤摘要】
一种幅度非衰减均衡相位的分数阶椭圆滤波器设计方法


[0001]本专利技术涉及滤波器设计领域，特别是指一种幅度非衰减均衡相位的分数阶椭圆滤波器设计方法。

技术介绍

[0002]在现代智能系统设计中，椭圆滤波器得到广泛应用。相比于其他滤波器，诸如切比雪夫I型和II型滤波器具有阻带和通带平滑特性，以及巴特沃斯滤波器通带和阻带同时具有平滑特性，椭圆滤波器牺牲了通带和阻带的平滑线和线性度以期获得过渡带陡峭的特性。理论研究表明，在滤波器其他性能条件相同时，椭圆滤波器具有最优的通带和阻带等纹波，以及最陡峭的过渡带。究其原因，它的滤波器传递函数同时存在零点和极点。大量的零极点位置相互制约加快了过渡带内信号衰减速率，但同时增大了设计难度，降低了设计效率。
[0003]自然界存在分数阶现象，分数阶运算在很多问题的处理上具有优于整数阶运算的特点。分数阶傅里叶变换在系统分析和现代信号处理领域应用越来越广泛，已经成为近期国内外学者研究热点。然而，目前研究和设计分数阶滤波器才刚起步。整数阶滤波器虽然能够滤除大部分杂散信号及噪声信号，但是通常会使其丢失大部分边缘低频有用信号。同时，由于滤波器的加入，使系统中的信号延迟增大，造成相位失真。因此，通过滤波器的信号的信噪比虽然增大，但是引入的失真也会极大降低滤波器输出信号的准确度。为了不产生信号相位的失真，通常要求滤波器具有零相位。零相位滤波器具有截止频率明确、计算量小的优点。当滤波器是零相位滤波器时，虽然滤波器的输入与输出没有相位差，但是存在幅度差，也即是滤波器的输入与输出的幅度不相等。/>[0004]由于传统椭圆滤波器都是直接或间接利用整数阶滤波器构建相应滤波器模型，较难平衡有效低频信号与抑制噪声信号。分数阶运算的幅度特性与相位特性均由分数阶次调节，将其应用到椭圆滤波器设计可以抑制低频噪声，对高频杂散信号有效去除的同时保留部分甚高频信号进行非线性保留。但是，由于分数阶微分的引入，导致椭圆滤波器的相位特性较差，输出时延较大。

技术实现思路

[0005]本专利技术的主要目的在于克服现有技术中的上述缺陷，提出一种幅度非衰减均衡相位的分数阶椭圆滤波器线性设计方法，该方法设计一种幅度非衰减的全通滤波器，校正分数阶椭圆滤波器的相位延迟，同时滤波器的幅度未衰减，因此椭圆滤波器陡峭的过渡带特性并未改变，滤波器的综合性能得以提高，具有较好的工程实用性。
[0006]本专利技术采用如下技术方案：
[0007]一种幅度非衰减均衡相位的分数阶椭圆滤波器设计方法,其特征在于，包括如下步骤
[0008]S1：算法初始化，设定分数阶椭圆滤波器的阶数为L(L>0)，最大通带纹波衰减为
Ω
p
dB，最小阻带纹波衰减为Ω
a
dB，选择系数k，设置初始迭代次数ξ＝1，椭圆滤波器最大迭代次数为ξ
max
，椭圆滤波器的设计参数初始值为搜索步长分别为
△
a
i
，
△
b
i
，
△
c
i
，设置最大迭代次数分别为ξ
1m
，ξ
2m
和ξ
3m
，迭代初始数值分别为ξ1＝0，ξ2＝0和ξ3＝0，幅度非衰减逆向滤波器的时间周期为N＝N0，采样周期为T＝T0；
[0009]S2：其中
[0010]S3：执行令如果则更新当前ξ1＝ξ1+1，否则
[0011]S4：执行如果则更新当前ξ2＝ξ2+1，否则
[0012]S5：执行如果则更新当前ξ3＝ξ3+1，否则
[0013]S6：对r个零极点重复执行步骤S3-S5，直至迭代次数同时满足ξ1＝ξ
1m
，ξ2＝ξ
2m
和ξ3＝ξ
3m
，输出分数阶椭圆滤波器的最优a
i
，b
i
和c
i
，i＝1,2,3,
…
,r。
[0014]具体地，步骤S2中具体包括：
[0015]定义为：
[0016][0017]具体地，步骤S2中f(a
i
,b
i
,c
i
)具体为：
[0018][0019]其中T为采样周期，z＝e
jΩT
为复数变量，j为虚数单位，Ω为频率，
[0020][0021][
·
]代表取整。
[0022]由上述对本专利技术的描述可知，与现有技术相比，本专利技术具有如下有益效果：
[0023]本专利技术提出一种幅度非衰减分数阶椭圆滤波器设计方法，使得椭圆滤波器的阶数由只能取整数转变为任意正实数。因此，滤波器的阶数调节变得更加灵活，同时，由于逆向滤波器的引入，有效补偿了滤波器的相位延迟，提高了椭圆滤波器的线性度。
附图说明
[0024]图1幅度非衰竭分数阶椭圆滤波器示意图。
[0025]以下结合附图和具体实施例对本专利技术作进一步详述。
具体实施方式
[0026]本专利技术提出一种幅度非衰减分数阶椭圆滤波器设计方法，使得椭圆滤波器的阶数由只能取整数转变为任意正实数。因此，滤波器的阶数调节变得更加灵活，同时，由于逆向滤波器的引入，有效补偿了滤波器的相位延迟，提高了椭圆滤波器的线性度。
[0027]L(L>0)阶分数阶椭圆低通滤波器的传递函数模的平方为：
[0028][0029]其中，R
L
(Ω)为改进的雅可比函数，(Ω)为改进的雅可比函数，ε为椭圆滤波器的纹波系数。
[0030]已知椭圆低通滤波器的最大通带纹波衰减为Ω
p
dB，最小阻带纹波衰减为Ω
a
dB，选择系数k，那么分数阶椭圆滤波器的设计目标为设计a
i
，b
i
和c
i
(i＝1,2,3,
…
,r)使得
[0031][0032]其中，[
·
]代表取整，t＝l+2l5+15l9+150l
13
，，
[0033]联立(1)和(2)，并对分数阶椭圆滤波器做双线性变换，将代入，得到
[0034][0035]其中T为采样周期，z＝e
jΩT
为复数变量(j为虚数单位，Ω为频率)。对(3)化简，得到
[0036][0037]椭圆滤波器的相位为
[0038][0039]对于幅度非衰减逆向滤波器：
[0040]假定分数阶模拟滤波器的传递函数为H(e
jΩT
)＝f(a
i
,b
i
,c
i
)，其中f(a
i
,b
i
,c
i
)为(4)所示，估计误差信号为m
e
(n)，n∈[0,N]。那么它首先经过一次时间周期为N的逆向变换后本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种幅度非衰减均衡相位的分数阶椭圆滤波器设计方法,其特征在于，包括如下步骤S1：算法初始化，设定分数阶椭圆滤波器的阶数为L(L>0)，最大通带纹波衰减为Ω
p
dB，最小阻带纹波衰减为Ω
a
dB，选择系数k，设置初始迭代次数ξ＝1，椭圆滤波器最大迭代次数为ξ
max
，椭圆滤波器的设计参数初始值为搜索步长分别为
△
a
i
，
△
b
i
，
△
c
i
，设置最大迭代次数分别为ξ
1m
，ξ
2m
和ξ
3m
，迭代初始数值分别为ξ1＝0，ξ2＝0和ξ3＝0，幅度非衰减逆向滤波器的时间周期为N＝N0，采样周期为T＝T0；S2：其中S3：执行令如果则更新当前ξ1＝ξ1+1，否则S4：执行如果则更新当前ξ2＝...

【专利技术属性】
技术研发人员：傅文渊，周静敏，
申请(专利权)人：华侨大学，
类型：发明
国别省市：