摘 要：提出一種在不增加分數階微分濾波器復雜度的前提下，能有效提高分數階微分濾波器性能的方法。該方法利用幾種基于典型微分算子的分數階微分濾波器之間的互補性，通過相互內插結合的方式，用于提高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算子表示為：