1 引入粒子聚合度的改進(jìn)粒子群優(yōu)化算法
1．1 粒子群優(yōu)化算法(PSO)
    PSO算法是美國Kennedy和Eberhart受鳥群覓食行為的啟發(fā)，于1995年提出的。該算法的思想是通過種群粒子間的合作與競爭，產(chǎn)生群體智能指導(dǎo)優(yōu)化搜索。PSO算法可用式(1)表示。

   
式中：vidk是粒子i在第k次迭代中第d維速度；xidk是粒子i在第k次迭代中第d維的位置；ω是慣性權(quán)值系數(shù)；pbestidk，是粒子i在第k次迭代中第d維個體極值點的位置(即個體最優(yōu))；gbestdk是整個種群在第k次迭代中第d維全局極值點的位置(即全局最優(yōu))；r1，r2是[0，1]之間的隨機數(shù)；c1，c2是加速系數(shù)，或稱學(xué)習(xí)因子。
1．2 帶粒子聚合度的改進(jìn)粒子群優(yōu)化算法
    由式(1)可知，如果粒子的當(dāng)前位置在gbest，此時個體極值點與全局極值點為同一點，即pbest與gbest相同。這時粒子速度若等于零，則種群的粒子將會出現(xiàn)進(jìn)化停滯，算法只能收斂到種群目前尋找到的最優(yōu)解gbest。假如這時gbest對應(yīng)的只是一個局部最優(yōu)解，那么算法就出現(xiàn)了早熟收斂現(xiàn)象。
    針對PSO算法存在早熟和局部收斂的問題，在基本PSO的基礎(chǔ)上，加入粒子聚合度n和一個線性遞減的慣性權(quán)值系數(shù)ω，對PSO算法進(jìn)行改進(jìn)。
    聚合度n是用來反映粒子群聚集程度的一個系數(shù)。當(dāng)粒子群出現(xiàn)高度聚集，進(jìn)化停滯時，n隨迭代次數(shù)遞增；當(dāng)n大于一個閾值λ(此閾值根據(jù)具體情況選擇)時，對粒子進(jìn)行變異，使變異粒子跳離當(dāng)前位置，進(jìn)入其他區(qū)域。在其后的搜索中，算法有新的個體極值pbest和全局極值gbest，從而跳出局部收斂。多次循環(huán)迭代后，就能找到全局最優(yōu)。
    改進(jìn)的算法可用式(2)和式(3)表示：

   
式(2)中rand是[0，1]間的隨機數(shù)：
式中：max Xd和min Xd分別是粒子在d維空間上的最大值和最小值。
    慣性權(quán)值系數(shù)叫決定控制算法的收斂特性，當(dāng)ω較大時，全局搜索能力強；當(dāng)ω較小時，局部搜索能力強。文獻(xiàn)[6]通過大量實驗證明，如果ω隨算法迭代的進(jìn)行而線性減小，將顯著改善算法的收斂性能。在此，?。?/p>

   
式中：(ωmax為最大慣性權(quán)值系數(shù)；ωmin為最小慣性權(quán)值系數(shù)；k為迭代次數(shù)；ksum為迭代總數(shù)。
    

2 用IMPSO設(shè)計FIR數(shù)字濾波器
2．1 FIR數(shù)字濾波器分析
    N階FIR數(shù)字濾波器的單位抽樣響應(yīng)為k(0)，k(1)，…，k(N-1)，其傳遞函數(shù)可表示為：

   
    取z=ejω，可得到數(shù)字濾波器的頻率響應(yīng)為：

   
    如果設(shè)計FIR數(shù)字濾波器的理想頻率響應(yīng)為Hd(ejω)，則設(shè)計濾波器與理想FIR濾波器的誤差e可通過對兩濾波器的幅度在一定量的離散點上的誤差平方和來表示，即取M個離散點時：

    由式(7)容易知得，誤差e是濾波器N個系數(shù)h(n)(n=0，1，…，N-1)的函數(shù)。對FIR濾波器的設(shè)計，就要選取合適的濾波器系數(shù)h(n)，使誤差e最小。顯然，h(n)的選取是一個組合優(yōu)化問題，因此可通過優(yōu)化算法求解濾波器系數(shù)h(n)，實現(xiàn)FIR設(shè)計。
2．2 適應(yīng)度函數(shù)
    IMPSO通過適應(yīng)度函數(shù)來確定粒子當(dāng)前位置的優(yōu)劣，因此選式(7)作為優(yōu)化設(shè)計FIR數(shù)字濾波器的適應(yīng)度函數(shù)。即：

   
    顯然，F(xiàn)ithess函數(shù)值越小，則對應(yīng)濾波器的幅度均方誤差就越小，該粒子就對應(yīng)更佳的濾波器系數(shù)。算法結(jié)束后，適應(yīng)度最小的粒子所代表的參數(shù)值就是濾波器的最優(yōu)系數(shù)。
2．3 算法編碼及流程
    為了用IMPSO算法求解h(n)，應(yīng)對優(yōu)化參數(shù)h(0)，h(1)，…，h(N-1)進(jìn)行適當(dāng)?shù)木幋a，以形成IMPSO算法中的粒子。算法用實數(shù)來表示各參數(shù)，h(0)，h(1)，…，h(N-1)分別表示N個粒子當(dāng)前的位置；vh(0)，vh(1)，…，vh(N-1)分別表示當(dāng)前粒子的速度；pbest(0)，pbest(1)，…，pbest(N-1)表示各粒子的個體最優(yōu)，gbest表示全體的最優(yōu)解。算法流程如圖1所示。

3 仿真算例
    為了驗證所提出算法的有效性，在計算機上采用Matlab語言進(jìn)行FIR數(shù)字濾波器設(shè)計的仿真實驗。同時為了比較算法性能，還采用基本遺傳算法(BGA)和基本粒子群優(yōu)化算法(BPSO)進(jìn)行了相同的濾波器優(yōu)化設(shè)計。仿真實驗中，粒子群優(yōu)化算法的參數(shù)設(shè)置為：群體大小Size=30，參數(shù)維數(shù)Codel=30，最大慣性權(quán)重為O．9，最小為0．2，聚合度為20，最大迭代次數(shù)為200；遺傳算法參數(shù)設(shè)置為：種群30，遺傳代數(shù)200，交叉概率0．5，變異概率0．1。
    實例1 設(shè)計一個低通FIR數(shù)字濾波器，其技術(shù)指標(biāo)如下：

   
    實例2 設(shè)計一個帶通FIR數(shù)字濾波器，其技術(shù)指標(biāo)如下：

   
    圖2和圖4分別是三種算法在設(shè)計FIR低通和帶通數(shù)字濾波器的參數(shù)優(yōu)化過程圖。

    圖3和圖5則是用三種不同算法設(shè)計的FIR低通和帶通數(shù)字濾波器。從圖2和圖4中容易看出，無論是對低通還是帶通濾波器的設(shè)計，因為IMPSO對BPSO加入了聚合度檢測，能進(jìn)行智能變異，同時采用線性遞減的慣性權(quán)值系數(shù)。因此與BPSO相比，IMPSO既有最快的尋優(yōu)速度，也具有最好的適應(yīng)值，只要迭代次數(shù)設(shè)置合理，在迭代次數(shù)范圍內(nèi)，粒子總會找到全局最優(yōu)值。BPSO與BGA的尋優(yōu)速度慢，容易陷入早熟收斂，很難得到理想的最優(yōu)參數(shù)。由圖3和圖5的FIR低通與帶通數(shù)字濾波器的幅頻特性曲線可明顯看出，利用IMPSO設(shè)計的濾波器在三種算法中最接近理想的濾波器，是較好的FIR濾波器設(shè)計方法。

4 結(jié) 語
    在此通過加入聚合度，并將遺傳算法中的變異思想引入到PSO算法中，對粒子實現(xiàn)智能變異，能有效克服標(biāo)準(zhǔn)．PSO容易進(jìn)入局部收斂的缺點。同時由于該算法用到的參數(shù)少，程序?qū)崿F(xiàn)簡單，因此與BGA等其他算法相比，具有運算量少，尋優(yōu)速度快等優(yōu)點。通過實例仿真結(jié)果表明，這里提出的IMPSO算法在FIR低通與帶通數(shù)字濾波器的設(shè)計上比BPSO和BGA具有更好的收斂速度和搜索能力。在設(shè)計FIR帶阻和多通帶數(shù)字濾波器時，IMPSO也是一種有效的設(shè)計方法。