方存洋， 施勝丹

(國電南瑞科技股份有限公司，江蘇 南京 211106)

基于SVG控制應用的數(shù)字濾波器實用設計

方存洋， 施勝丹

(國電南瑞科技股份有限公司，江蘇 南京 211106)

IIR型數(shù)字濾波器憑借其優(yōu)良的衰減特性在利用軟件處理信號的場合應用廣泛。針對無功發(fā)生器(svg)控制保護裝置應用程序開發(fā)過程中所使用的IIR數(shù)字濾波器，利用MATLAB中的FDAtool工具獲得階數(shù)和傳遞函數(shù)系數(shù)，并通過微調零極點位置獲得較為理想的幅頻曲線，并輸出系數(shù)寫入控制程序的函數(shù)中；分析采樣周期對濾波器零極點位置的影響，在不同采樣周期下，數(shù)字濾波器零極點位置發(fā)生改變，產(chǎn)生新的傳遞函數(shù)系數(shù)。經(jīng)過在實際裝置中的測試，獲得了良好的濾波效果。

數(shù)字濾波器；采樣周期；FDAtool；零極點；系數(shù)微調

0 引 言

數(shù)字濾波器所具有的運用靈活方便，幅頻、相頻特性可以分別給出等優(yōu)點，使得利用程序實現(xiàn)的數(shù)字濾波器成為信號處理的主要手段之一[1]。常用的數(shù)字濾波器有無限沖擊響應數(shù)字濾波器(IIR)和有限沖擊響應數(shù)字濾波器(FIR)。IIR輸出為過去輸出與現(xiàn)時輸入的加權和，計算量遠小于FIR，且可借助模擬濾波器進行設計，有大量圖表可查，方便簡單，但其幅頻特性受采樣周期的影響很大，存在混迭效應；FIR相位線性，但要取得好的衰減特性，一般階次要高[2-4]。在實際產(chǎn)品開發(fā)過程中，對計算量、響應時間等有較高的要求，為此一般采用IIR型數(shù)字濾波器。在應用IIR數(shù)字濾波器時，不同的采樣周期對應不同的傳遞函數(shù)系數(shù)，且由于其輸出與過去輸出有關，對系數(shù)要求非常高，一旦出現(xiàn)偏差(或舍入誤差不合適)，將影響濾波器效果，甚至造成系統(tǒng)不穩(wěn)定。本文分析了采樣周期變化對IIR數(shù)字濾波器傳遞函數(shù)中系數(shù)的影響，借助于MATLAB仿真軟件中的FDAtool工具計算系數(shù)，并在SVG控制保護裝置的實際程序中應用，通過錄波獲得濾波效果的偏差，以微調零極點的位置，最終獲得此控制系統(tǒng)下理想的濾波性能。

1 IIR數(shù)字濾波器原理

采用直接Ⅱ型網(wǎng)絡結構(如圖1所示)的N階IIR數(shù)字濾波器，其傳遞函數(shù)的標準形式為[5-6]：

(1)

上式具有N個零點zi和N個極點pi。將式(1)寫成如下形式：

(2)

若|pi|<1，當z→∞時，H(z)→0，系統(tǒng)穩(wěn)定，此時極點位于單位圓內；若|pi|>1，當z→∞時，H(z)→∞，系統(tǒng)不穩(wěn)定，此時極點位于單位圓外。

圖1 IIR直接Ⅱ型網(wǎng)絡結構圖

2 采樣周期對系數(shù)的影響分析

設計IIR低通濾波器，需要知道以下目標參數(shù)：采樣周期，通帶邊緣(截止)頻率fp，阻帶邊緣(起始)頻率，通帶內波動，阻帶內最小衰減。設計一個IIR數(shù)字濾波器的步驟如圖2所示。

圖2 IIR數(shù)字濾波器設計步驟示意圖

其中采樣周期對IIR數(shù)字濾波器性能的影響主要體現(xiàn)在模擬濾波器離散化過程中，模擬濾波器的離散化主要有沖擊不變變換法和雙線性變換法[7]。

模擬濾波器傳遞函數(shù)的標準形式為：

(3)

其中ξi是H(s)在pi處的留數(shù)。

對式(3)進行拉普拉斯逆變換，可得：

(4)

對h(t)進行抽樣，可得差分方程：

(5)

對式(5)進行z變換即可得數(shù)字濾波器的傳遞函數(shù)：

(6)

對比式(3)和式(6)，發(fā)現(xiàn)與模擬濾波器極點pi對應的數(shù)字濾波器極點為epiTs。由此可見，在IIR數(shù)字濾波器設計過程中，傳遞函數(shù)(濾波器性能)受采樣頻率影響。

在實際應用程序中，數(shù)字濾波器的實現(xiàn)依靠其傳遞函數(shù)實現(xiàn)，在已知階數(shù)得到傳遞函數(shù)的標準形式后，獲得一組傳遞函數(shù)系數(shù)即可完成一個數(shù)字濾波器的設計。由式(2)可知，傳遞函數(shù)分母系數(shù)ai=epiTs，一旦采樣周期發(fā)生變化，傳遞函數(shù)的分母系數(shù)也要跟著改變。

令原來的采樣周期為Ts1，新的采樣周期為Ts2，對應的系數(shù)為ai1=epiTs1，則pi=(logai)/Ts1，濾波特性不變，則模擬濾波器的極點位置不變。

(7)

3 系數(shù)計算

3.1 系數(shù)的基本確定步驟

在實際應用中，設計數(shù)字濾波器的具體步驟如圖3所示。

圖3 實際應用中數(shù)字濾波器的設計過程

由圖3 可見，在FDAtool工具中輸入目標參數(shù)，獲得數(shù)字濾波器階數(shù)和傳遞函數(shù)系數(shù)后，寫入實際應用程序中獲得的濾波波形會與理想情況下有偏差。這是由于實際程序存在有限字長效應，運算過程中將對系數(shù)進行舍入處理，從而引起極點的偏移，這種情況很容易帶來穩(wěn)定性問題，甚至造成寄生振蕩。為此，需要在獲得系數(shù)后，結合程序運行帶來的實際效果進行微調，使得系數(shù)舍入誤差在穩(wěn)定范圍內。

在實際數(shù)字濾波器設計過程中，必須進行傳遞函數(shù)系數(shù)的微調。本文采用零極點位置微調法，具體為：

1)利用FDAtool工具進行系數(shù)計算，同時查看此刻零極點位置，幅頻、相頻特性曲線(FDAtool界面如圖4所示)；

2)與理想的幅頻、相頻特性曲線進行比較(主要查看截止頻率)，如不滿足要求，可通過移動零極點位置進行修正，直至獲得希望的幅頻、相頻特性曲線(性能曲線通過圖4所示右上角圖標查看)；

3)查看濾波器傳遞函數(shù)的各系數(shù)ai、bi，寫入實際程序后在裝置中運行，錄波查看數(shù)字濾波器輸出波形，分析諧波以驗證此濾波器性能。若濾波效果不理想，返回第2)步，直至獲得理想濾波效果。

圖4 FDAtool工具菜單界面

3.2 僅采樣周期變化時系數(shù)的確定方法

在SVG開發(fā)過程中，發(fā)生僅采樣周期改變，而濾波性能要求不變的情況，可以通過兩種方式獲得新的傳遞函數(shù)系數(shù)：

1)在第一步通過目標參數(shù)獲得原始傳遞函數(shù)系數(shù)時，即輸入新的采樣周期，而后按照3.1節(jié)所述步驟進行濾波器設計；

2)根據(jù)式(7)獲得傳遞函數(shù)系數(shù)，再利用零極點位置微調獲得較為理想的濾波器。

4 零極點微調方法

在系數(shù)計算部分已經(jīng)提到，在根據(jù)目標參數(shù)獲得濾波器階數(shù)和傳遞函數(shù)原始系數(shù)后，需要通過微調零極點位置獲得相對理想的濾波效果。FDAtool調整零極點位置的界面如圖5所示。

圖5 零極點位置調整界面

零極點的個數(shù)與濾波器階數(shù)有關，此處為系數(shù)微調，為此一般零極點的個數(shù)保持不變，僅通過移動零極點位置完成調整。如圖5所示，界面上半部分顯示幅頻特性曲線，界面下半部分為零極點調整界面，零極點位置發(fā)生改變后，可在上半界面的特性曲線中查看截止頻率等信息。

零極點位置微調過程中可參考如下規(guī)律：

1)極點必須在單位圓內部；

2)單位圓附近的零點會在濾波器幅頻特性的相應頻率處產(chǎn)生陷落，零點離單位圓越近，陷落越深；

3)單位圓附近的極點會在濾波器幅頻特性的相應頻率處產(chǎn)生凸峰，極點離單位圓越近，凸峰越高。

5 實驗驗證

SVG控制程序中，采樣周期為10 kHz，主要使用低通濾波器，階數(shù)一般為1或者2，以圖6(a)所示波形為例進行濾波，設計一個截止頻率為80 Hz左右的IIR數(shù)字濾波器。在MATLAB的FDAtool中輸入數(shù)字濾波器的目標參數(shù)，可得濾波器階數(shù)為2，該濾波器傳遞函數(shù)的系數(shù)為：b0=0.000 609 854 7，b1=0.001 219 709 4，b2=0.000 609 854 7；a0=1，a1=-1.928 942 26，a2=0.931 381 68。

將此組系數(shù)寫入程序中的二階濾波器函數(shù)，對輸出波形錄波后發(fā)現(xiàn)，高頻諧波含量明顯下降，但2次諧波增大?？梢?，此組系

數(shù)在本控制系統(tǒng)中不合適，引起了共振，為此進行零極點位置微調，最終獲得較為理想的波形，此時的系數(shù)為：b0=0.000 537 169 8，b1=0.001 074 339 6，b2=0.000 537 169 8；a0=1，a1=-1.933 380 2，a2=0.935 528 9。濾波前后的波形及諧波含量情況分別如圖6～圖9所示。

圖6 濾波前后電流波形及其諧波含量分析

圖7 濾波前諧波含量分析

圖8 與系統(tǒng)共振時諧波分布

圖9 諧波后諧波含量分布

采樣周期變?yōu)? kHz，傳遞函數(shù)的系數(shù)為：b0=0.000 660 779，b1=0.001 3215 58，b2=0.000 660 779；a0=1，a1=-1.925 984，a2=0.928 627，獲得的濾波后波形如圖10所示。

圖10 采樣周期為9 kHz時，濾波后電流波形

6 結束語

利用MATLAB軟件中的FDAtool工具可簡化IIR型數(shù)字濾波器的設計步驟，而結合實際錄波獲得的濾波效果，進行零極點位置微調則簡單直觀的抵消了有限字長帶來的混迭效應，獲得理想的濾波效果；采樣周期對傳遞函數(shù)的系數(shù)有影響，通過系數(shù)調整，可獲得相同的濾波效果。數(shù)字濾波器應用廣泛，分析采樣周

期對傳遞函數(shù)系數(shù)的影響、掌握其實用設計技術，具有深刻的意義。

[1] 王小艷.靜止無功發(fā)生器的電流檢測和控制策略研究[D].西安：西安理工大學,2010.

[2] 王群,姚為正,王兆安.高通和低通濾波器對諧波檢測電路檢測效果的影響研究[J].電工技術學報,1999,14(5):22-26.

[3] 周柯,羅安,夏向陽,等,一種改進的id-iq諧波檢測方法及數(shù)字低通濾波器的優(yōu)化設計[J].中國電機工程學報,2007,27(34):96-101.

[4] 程遠楚,徐德鴻,劉燕,等.基于瞬時無功理論的svc低通濾波器設計與應用[J].電工技術學報，2008, 23(9):138-143.

[5] 胡廣書. 數(shù)字信號處理—理論算法與實現(xiàn)[M]. 2版.北京:清華大學出版社,2003.

[6] 丁玉美,高西全.數(shù)字信號處理[M].西安:西安電子科技大學出版社,2005.

[7] 劉翔宇,楊仁剛,王文成.基于id-iq法諧波檢測中數(shù)字低通濾波器的設計及其DSP實現(xiàn)[J].電力自動化設備,2006,26(8):51-54.

Practical Design of a Digital Filter Based on SVG Control Application

Fang Cunyang, Shi Shengdan

(NARI Technology Co., Ltd., Nanjing Jiangsu 211106, China)

IIR digital filters find wide application in signal processing through software thanks to their excellent attenuation characteristics. With regards to IIR digital filters used in the development of application programs for the control protection device for Static Var Generator (SVG), FDAtool in Matlab is used to provide the order number and transfer function coefficient, a quite ideal amplitude-frequency curve is obtained through fine-tuning of the pole zero position in the tool, and the output coefficient is written into the function of the control program. This paper further analyzes the effect of sampling period upon pole zero position of the filter. For different sampling periods, the digital filter’s pole zero position differs from each other, thus generating a new transfer function coefficient. A good filtering effect is achieved as shown in the test on actual devices.

digital filter; sampling period; FDA tool; pole-zero; coefficient fine-tuning

10.3969/j.issn.1000-3886.2016.06.009

TN713+.7

A

1000-3886(2016)06-0027-04

方存洋(1968-)，男，安徽人，本科，工程師，主要研究方向：電力電子無功補償技術。

定稿日期: 2016-06-09