最小平方差预估应用於主动前馈式噪音消法实做研究.docVIP

簡介 本文以一維長管為受控制系統，控制標的為由長管一端經由長管傳到另一端出口的噪音。在理論方面，我們以格子濾波器演算法則及主動噪音消除理論為基礎，在實際訊號處理及控制訊號演算方面則以ADSP2101定點數位訊號處理器為核心，其它量測及控制硬體包括兩個麥克風作為訊號感測器，一顆低音喇叭作為音器，前後級訊號放大器，以及AD/DA介面電路作為數位類比轉換之用。在實際測試之前，我們以質塊彈簧系統的電腦模擬系統來替代實際受控系統，進行初級可行性測試，如事先預期，電腦模擬總會得到不錯的結果。對實際一維長管系統的實際控制，我們除了上述的電硬體外，　另加一組喇叭及放大器作為噪音源，實測時噪音源的電路與控制電路獨立，經由實際測試所得結果顯示以格子濾波器演算法作主動噪音消除，的確有效果，但是並不如控制電腦模擬系統來得大。 數學模型 2.1 實際系統硬體配置 本實驗以汽車排氣管設計作為考量，在長管末端加裝主動式消音器，並在控制喇叭中間穿洞，使得長管可直通排氣，如此可以消音，亦不會像交錯放置的消音器，產生妨礙氣流通路、增加氣流阻力的被動消音之缺點。實驗裝置如圖一所示，假設有一單一噪音源從左至右通過長管，並由數個單一指向性的麥克風所收音。聲音波前依序經過M1、M2，在經過一段空氣衰減及時間延遲後，由M3收集音波。假設M1、M2所收集到的聲音訊號並不會被M3訊號所 影響，如此則M1、M2訊號間即存在FIR模型的關係。此外，為避免當作為破壞性干涉音源的喇叭所發出控制訊號，經由長管反向傳播影響到M1所收集到的訊號，而使得M1、M2訊號間不能符合單純的FIR模型。故噪音源喇叭及控制的喇叭擺置時向著同一方向，同時所有的麥克風均使用單一指向式，並將指向對著噪音來源方向。 圖一 實際系統示意圖 2.2 質量彈簧模擬系統 在測試控制實際長管系統之前，我們先將我們所提出的控制方法試用於一質量、彈簧系統的電腦模擬系統，作為初步的測試。此一模擬系的的示意如及參數如圖二所示。 圖二 質量彈簧模擬系統 系統設定: D Position Input 位置輸入 ζ Damping Coefficient 阻尼係數 K Spring Constant 彈簧彈性係數 ω Block Position 各Block 的位置 M Block Mass 各Block的質量 u Control Force 控制端之力 y Output 位置輸出 將上述的系統以State Space Model來表示可得: State space model : 主動式消音器之設計分析 3.1 主動式雜訊消除控制迴圈 雜訊消除演算法控制迴圈如下: 圖三 演算法控制迴圈 其中y為噪音源，轉移函式H為長管的動態特性，為未知，轉移函式G為致動器的動態特性，將透過事先的系統鑑別得知，而u為前饋控制命令，e為系統的輸出，亦為修正適應性濾波器參數的依據。 此演算法中，依據Track Error e與M1所感測到的訊號送入適應性系統以計算控制訊號。目的是希望將系統輸出e達到最小。 在此處，適應性控制器的修正法則為FIR格子濾波器。由以上演算法方塊器的表示中，若與是兩資料序列，n階自適應FIR格子濾波器鑑別出一個n階線性FIR濾波器，使得以下的式子得到最小: 其中‖?‖2表示均方誤差。 此外線性FIR濾波器的一般表示式如下: 值得注意的是，在此處由於格子濾波器在整理FIR模型系統參數時，十分耗時，所以我們直接使用格子濾波器的反射係數(Reflection Coefficients)做更新疊代，而不用FIR濾波器的係數ai。 雜訊消除演算法中各參數的關係如下: 根據圖三可得到以下兩 個關係式: 輸出誤差方程式(Output Error Equation) 前饋控制法則(Feedforward Control Law) 其中，F(z)為自適應系統控制器之FIR濾波器，而HLP為一低通濾波器，由於所要消音的雜訊以低頻訊號為主，故加此低通濾波器可使訊號不會因為環境或是電路所產生的高頻雜訊受到干擾。 由於系統特性的轉移函數未知，所以運用格子濾波器時的兩輸入訊號必須以下式決定: 由上面的式子，我們可以得到: 因此，當格子濾波器收斂而使得平均平方誤差最小時，系統輸出即為最小。 as F converges. 3.2預視窗系統鑑別格子濾波器 格子濾波器是根據最小平方問題，將所需處理的資料，在空間上做連續的垂直