一种分数阶微分IIR滤波器的设计方法和改进

摘 要：提出一种在不增加分数阶微分滤波器复杂度的前提下，能有效提高分数阶微分滤波器性能的方法。该方法利用几种基于典型微分算子的分数阶微分滤波器之间的互补性，通过相互内插结合的方式，用于提高IIR分数阶数字滤波器的性能。改进后的分数阶微分滤波器频率响应更接近理想分数阶微分滤波器，表明所提方法的有效性。
关键词：分数阶微积分；数字微分器；IIR滤波器；微分算子；连续分数扩充

0 引 言
分数阶微积分是一个既古老又现代的话题。早在整数阶微积分产生的时候分数阶微积分就产生了，该问题曾被许多数学家，如Leibniz(1695)，Euler(1738)，Liouville(1850)，Hardy和Littlewood(1925)等涉及和探究过。虽然分数阶微积分的研究难度很大，但近三百年在众多科学家的不懈努力下，分数阶微积分作为纯数学分支已经发展渐成体系，但其物理意义不明确，阻碍了分数维微积分的应用，目前在工程技术界中没有得到广泛应用。从Mandelbrot提出分形学说，将Rie―mann―Liouville分数阶微积分用以分析和研究分形媒介中的布朗运动以来，分数阶微积分才在许多学科，特别是在化学、电磁学、控制学、材料科学和力学中引起广泛关注并尝试着应用。随信息科学的变革和迅猛发展，分数阶运算在很多问题的处理过程中所拥有整数阶运算无可比拟的优点正逐渐显露出来。
目前分数阶滤波器已经在分数阶控制器、信号处理、图像压缩和处理等领域得到成功应用。分数阶数字分数阶微分滤波器的设计和改进，正成为分数阶微积分研究领域的一个热点。数字微分滤波器的设计方法通常可以归为2类：第一种是线性相位F1R滤波器方法；另一种是IIR滤波器法。考虑到滤波器设计复杂度因素，FIR微分滤波器阶数会受到限制，影响了其频率响应对理想频率响应的逼近效果，因此这里考虑使用IIR分数阶微分滤波器来实现分数阶运算。
IIR分数阶数字微分滤波器设计的重点是实现分数阶算子的离散化，即是找到一个函数Gv(z)，使其频率响应无限逼近理想分数阶数字微分器的频率响应Hv(ω)=(jω)v。基本步骤可以归纳为：首先，找到频率响应接近理想一阶微分的算子；然后基于所选用的微分算子，推导出分数阶微分滤波器传输函数；最后通过各种展开方法把传输函数的分数阶形式转化为整数阶滤波器形式。完成分数阶展开的常用方法有幂级数展开(PSE)和连续分数扩充(CFE)，其中连续分数扩充方法对函数的逼近更好，收敛更快。
首先对Rectangular算子、Tustin算子、Simpson算子这几种典型微分算子通过连续分数扩充，得到相应的0．5阶微分滤波器频率响应。通过分析这几种算子的频率响应表明，基于这几种典型算子的分数阶微分滤波器各有优缺点和具有互补性，将这几种典型算子进行结合可得到更接近理想分数阶微分算子频率响应的算子。

1 典型IIR分数阶微分滤波器
1．1 基于Simpson算子的IIR分数阶数字微分滤波器
Simpson微分算子表示为：

GvSn(z)中v表示微分阶数；n表示滤波器阶数。图1是基于Simpson算子的O．5阶微分滤波器的频率响应曲线图。

在此使用连续分数扩充(CFE)方法完成对上式的展开，这里简要介绍分数阶算子实现过程中使用到的CFE方法。对于任何一个函数D(z)，可以用下面连续分数的形式来表示：

式中，系数ai，bi是关于变量z的有理函数或常数。只需要通过截断操作，就能得到有限阶逼近函数。下面列出T=0．001 s时，使用连续分数扩展(CFE)完成上式的展开，得到0．5阶微分的Simpson分数阶微分滤波器传递函数GvSn(z)：

通过对比和分析，从误差和计算复杂度两个方面均衡考虑分数阶微分滤波器阶数的选为5阶比较合适。因此这里滤波器的阶数都选为5阶。
1．2 基于Rectangular算子的IIR分数阶数字微分滤波器
．Rectangular算子表示为：

基于Rectangular算子的分数阶微分器传输函数可以写为：

这里使用连续分数扩充(CFE)法将展开上式，实现对函数的有限阶逼近。下面列出T=0．001 s时，O．5阶微分Rectangular分数阶微分滤波器传递函数GvRn(z)：

GvRn(z)中v表示微分阶数；n表示滤波器阶数。
1．3 基于Tustin算子的IIR分数阶数字微分滤波器
Tustin算子表示为：

基于Tustin算子的分数阶微分器传输函数可以写为：

使用连续分数扩充(CFE)方法将上式展开，完成对函数的有限阶逼近。下面列出了T=0．001 s时，0．5阶微分Tustin分数阶微分滤波器传递函数GvTn(z)：

GvTn中v表示微分阶数；n表示滤波器阶数。
图2是基于典型Rectangular算子、Tustin算子和simpson算子的0．5阶微分滤波器的频率特性曲线，所实现的滤波器阶数都是5阶。从图2中可以看出3种滤波器在低频区域，幅度曲线还能与理想幅度一致，但随着频率增加，特别是在高频区域，误差迅速增大。