• <input id="zdukh"></input>
  • <b id="zdukh"><bdo id="zdukh"></bdo></b>
      <b id="zdukh"><bdo id="zdukh"></bdo></b>
    1. <i id="zdukh"><bdo id="zdukh"></bdo></i>

      <wbr id="zdukh"><table id="zdukh"></table></wbr>

      1. <input id="zdukh"></input>
        <wbr id="zdukh"><ins id="zdukh"></ins></wbr>
        <sub id="zdukh"></sub>
        公務員期刊網 精選范文 數學建模魯棒性分析范文

        數學建模魯棒性分析精選(九篇)

        前言:一篇好文章的誕生,需要你不斷地搜集資料、整理思路,本站小編為你收集了豐富的數學建模魯棒性分析主題范文,僅供參考,歡迎閱讀并收藏。

        數學建模魯棒性分析

        第1篇:數學建模魯棒性分析范文

        關鍵字:鎖相環;數學建模;穩定性分析

        1概述

        在VLSI系統中,鎖相環常被用于生成高頻時鐘。由于系統工作頻率在不斷的升高,鎖相環的輸出頻率也必須相應的提高,此外隨著工藝尺寸的縮小,晶體管集成度在不斷增加,PLL需要工作在更嘈雜的噪聲環境中,電荷注入和時鐘饋通等帶來的紋波也增加了鎖相環的抖動,所有這些都使得鎖相環的低抖動設計面臨著更大的挑戰。

        在文獻[1]中,John G. Maneatis提出了一款2階自偏置鎖相環,該系統具有大的工作帶寬,優良的低抖動特性和很好的魯棒性能。但是不足的是該PLL仍然是二階系統。二階系統濾波器設計簡單,系統穩定性較好,但是對噪聲的抑制效果遠沒有高階的好,不能滿足更高要求的低抖動設計。

        高階鎖相環可以有效地衰減高頻分量和抑制紋波,具有較好的低抖動性能,但是高階鎖相環存在多個極點,極大的增加了系統穩定性的設計難度。在系統級對PLL進行穩定性分析以及濾波器參數選擇也變得更加復雜。本文將基于一款三階自偏置鎖相環進行分析,建立其系統級數學建模,推導出該鎖相環的性能參數,并根據此數學模型,得到系統響應的波特圖,依此分析了該系統的穩定性,并由此來確定濾波器的參數。利用數學模型可以簡潔準確的分析系統特性,并有效的指導電路設計。

        本文第2節介紹了三階自偏置鎖相環的結構,第3節建立了其的數學模型并量化分析其系統性能參數,第4節中,利用波特圖分析了該系統的穩定性能,第5節給出了系統級和電路級的仿真結果。

        2三階自偏置鎖相環的結構

        三階自偏置電荷泵鎖相環的結構如圖1所示,它由鑒頻鑒相器(PFD)、電荷泵(CP)、環路濾波器(LF)、壓控振蕩器(VCO)、偏置電路(BG)和一個N分頻器構成。與傳統的鎖相環不同的是,這里增加了一個偏置電路(BG)部件,該部件可以自己選擇最佳工作偏置點,該偏置點只與工作頻率相關[1]。

        如圖1所示,PFD檢測VCO分頻后的時鐘與參考時鐘之間的相位差,其檢測結果驅動著電荷泵充電或者放電。環路濾波器將電流轉換成控制電壓,并濾掉交流電壓分量,產生一個穩定的控制電壓給偏置電路,偏置電路產生一個對噪聲不敏感的VCO控制電壓。當鎖相環鎖定時,VCO就會產生一個N倍于參考時鐘振蕩頻率的穩定時鐘。

        在現在的VLSI設計中,各種各樣的噪聲增加了低抖動的設計難度。此外,在PFD檢測到相位差時,電荷注入和時鐘潰通都引起較大的紋波,增加了抖動。出于低抖動性能的考慮,環路濾波器(LF)用了二階濾波器代替了常用的一階濾波器,所以整個鎖相環是一個三階負反饋系統,這樣可以有效地抑制溫波,減小抖動。接下來的兩章,將詳細介紹該鎖相環的數學建模和穩定性分析。

        3數學建模

        鎖相環可以用一個連續時間負反饋系統來表示,其性能可轉移到頻域區分析。本章首先將建立頻率響應函數,根據響應函數得到詳細的性能參數。

        3.1 環路濾波器(LF)響應函數

        環路濾波器的作用是將電流轉換成電壓,同時濾除高頻噪聲,對鎖相環的性能起決定性作用?,F代鎖相環中一般采用結構簡單的無源濾波器,這里我們使用了兩階濾波器,如圖2所示。

        圖2是一個二階無源阻抗型的濾波器,是一個將電流轉換到電壓的積分電路,其傳遞函數為:

        F■(s)=■ (1)

        3.2 PLL系統響應函數

        該PLL是一個三階負反饋系統。我們設定I■為電荷泵的電流,K■為電荷泵的增益,N為分頻系數。輸出函數P■(s)與P■(s)輸入的關系可以寫為:

        P■(s)=(P■(s)-P■(s)/N)

        ×■■K■■(2)

        由此可以得到PLL系統閉環傳輸函數Hclose(s)和開環傳輸函數Hopen(s)分別是:

        H■(s)=■

        =■(3)

        H■(s)=■■(4)

        3.3 量化分析

        在傳統的鎖相環中,I■ 和R都是定量,這使得決定系統性能的兩個參數ζ(阻尼因子)和 ωN(帶寬)也是一個定值。為了減小抖動,要求帶寬越大越好,但是為了保證系統的穩定性,帶寬需要做到工作頻率的十分之一以下[1]。這就導致了穩定性和抖動性能之間的沖突,在設計的時候必須保證最低工作頻率下的穩定性,有可能使得在較高頻率工作時的抖動性能不理想。而在自偏置鎖相環中,IP和R都是變量,帶寬與工作頻率的比值是一個定值,這樣在高工作頻率時也能得到很好的抖動性能。三階以上的高階鎖相環其帶寬和阻尼因子的計算相當復雜,但是如果C2是C1的1/5以下時就可以將此三階系統近似為兩階系統來分析[2]。所以公式(1)可以近似為:

        F■(s)=R■+■(5)

        那么開環傳輸函數可以簡化為:

        P■(s)=(P■(s)-■)?■

        ?R+■?K■?■(6)

        公式6可以變形為:

        ■=N?■ (7)

        其中,

        ?灼=■?■(8)

        ω■=■(9)

        在自偏置鎖相環中, R與■成反比,■與輸出頻率成正比[1],所以在鎖定的情況下,■與參考時鐘也成正比。這樣ζ和ω■ /ωREF都是固定的值。所以這個PLL在滿足穩定性的同時,在整個工作帶寬內都可實現良好的低抖動特性。

        4系統穩定性分析

        4.1穩定性概念

        對于負反饋系統,它的閉環函數可以寫為:

        ■(s)=■ (10)

        我們注意到如果NH(s)=-1,上式的增益將趨于無窮大,電路可以將自身的噪聲無限放大直到振蕩。所以反饋系統穩定必須滿足當H(jω)=1 時有下式成立[2]:

        ∠H(jω)< -180°,當H(jω)=1時(11)

        分析系統穩定性的方法有很多,其中較為簡單實用的是波特圖的方法??梢园压?變形為:

        H■(s)=K■(12)

        其中T■=C■R, K=■, T■=■

        這個系統有三個極點,ω2=1/T2是其中一個,另外兩個都在原點。是ω1=1/T1唯一的零點。ω2>ω1必須恒定成立,以保證系統的穩定性。依據公式(3)和(4),得到系統的波特圖如圖3所示。

        可以看出,該鎖相環在原點處有兩個極點,幅頻特性呈二階衰減特性下降,相移也達到-180°,然后系統引入的零點ω1使鎖相環呈單階衰減的特性,相位朝著-90°的方向移動,即有了相位裕度,使得鎖相環能夠穩定成為可能。為了有更好的高頻衰減特性,環路濾波器引入的另一個極點ω2,使鎖相環再次進入二階衰減的特性,所以高頻噪聲得到有效衰減。

        本文為全文原貌 未安裝PDF瀏覽器用戶請先下載安裝 原版全文

        從圖中也可以看出穩定性能和抖動性能的折中關系。如果ω2遠大于ω1,我們能得到較大的相位裕度,這樣系統等穩定性會很好,但是對高頻噪聲的去除效果并不好。如果ω2離ω1較近,對高頻噪聲抑制較好,但是穩定性能較差。所以濾波器參數的設計體現了穩定性、抖動性能等參數之間的折中,通常會取相位裕度在55°左右[2]。

        為了確定合理的濾波器參數,重要的是要知道在ω1

        用s=jw帶入公式(4)中可得到

        H■(s)=■■ (13)

        相位裕度公式為:

        φ■=arctg(wC■R)=arctg■+180°(14)

        將公式(14)進行求導取最大值,相位裕度函數的最大值處于頻率ωopt,它是ω1和ω2的幾何平均數,即:

        ω■=■(15)

        我們希望截止頻率具有足夠的相位裕度,因此選取截止頻率等于ω■是最佳的選擇,在實際的高性能PLL設計中,已經證明,選擇ω2/ω1比值近似為10較為有利[2],我們既可以獲得足夠的相位裕度,保證鎖相環系統穩定,而且還不至于使三階極點離中心太遠,造成去抖效果不理想。

        5仿真結果

        基于此數學模型,在采用matlab對PLL系統級進行仿真,其VCO的控制電壓曲線的仿真波形如圖4所示,可見系統大約在0.7微秒就已經穩定下來,控制電壓沒有出現寬幅的震蕩,說明此系統具有較好的穩定性。

        根據此數學模型,在65nm工藝下實現了自偏置鎖相環的電路和版圖設計。使用hspice對電路進行仿真,控制電壓的仿真曲線如圖5所示,可以看到PLL的控制電壓在0.6微妙左右也穩定下來,工作曲線與數學模型的曲線基本吻合。

        圖6是該PLL整體版圖,其中關鍵部件只占很少的面積,大部分面積被電容和去耦電容占據。

        其輸出波形及抖動測量如圖7所示,仿真工具為hspice,仿真條件是:電壓為1V,溫度為75℃,corner為TT。

        表1中列出了三個輸出頻率下的性能參數,仿真條件同上。由于工作條件和初始狀態的不同,鎖定時間也略有差別,在整個工作頻率范圍內,抖動都能控制在15ps以內。

        此外,為了測試該電路的魯棒性能,分別對電壓和電容有偏差的情況作了仿真,仿真結果顯示:壓差在15%以內,容差在30%以內時,該電路均可正常工作。

        6總結

        本文從系統級對三階自偏置鎖相環進行了數學建模,以一種簡潔的方法分析了高階鎖相環的性能參數以及穩定性。從仿真結果來看,這種建模分析的方法可以有效的表征系統特性,并指導電路設計?;诖藬祵W模型完成了三階自偏置鎖相環的電路和版圖設計,其輸出頻率范圍為400MHz~2GHz,在工作范圍內抖動可控制在15ps以內。

        本文是從系統級對鎖相環進行分析,沒有涉及到電路的設計,且對該鎖相環自偏置的原理沒有過多的說明。這些都是值得研究和探討的地方,將在后續的工作中繼續研究。

        參考文獻

        [1] Roland E. Best, Phase-Locked Loops Design, Simulation, And Applications, McGraw-Hill, 1999

        [2] J.Maneatis, Low_Jitter and Process-Independent DLL and PLL Based on self-Biased Techniques ISSCC, 1996

        [3] John G. Maneatis and Mark A. Horowitz, Precise Delay Generation Using Coupled Oscillators, IEEE JSSC, Vol. 28. No 12 . Dec 1993

        [4] 周潤德 等譯, 數字集成電路, 電子工業出版社, 2004.

        第2篇:數學建模魯棒性分析范文

        關鍵詞:卡里巴水壩 拆除 替代 詳細分析 新的水壩系統

        中圖分類號:TV631 文獻標識碼:A 文章編號:1672-3791(2017)04(a)-0095-04

        Detailed Analysis of the Acquisition and Optimization of the New Kariba Dam System

        Yang Shuhan

        (The Dalian university of technology, Dalian Liaoning, 116024, China)

        Abstract: The Kariba Dam on the Zambezi River is one of the larger dams in Africa. Its construction was controversial, and a 2015 report by the Institute of Risk Management of South Africa included a warning that the dam is in dire need of maintenance. A number of options are available to the Zambezi River Authority (ZRA) that might address the situation. One of the options is: Removing the Kariba Dam and replacing it with a series of ten to twenty smaller dams along the Zambezi River. This paper provides the detailed analysis of this option. That is under the condition that this new system of dams have the same overall water management capabilities as the existing Kariba Dam while providing the same or greater levels of protection, giving water management options for Lake Kariba that are in place with the existing dam. The analysis includes a strategy for modulating the water flow through the new multiple dam system that provides a reasonable balance between safety and costs. In addition to addressing known or predicted normal water cycles, the strategy should explains and justifies the actions that should be taken to properly handle emergency water flow situations (i.e. flooding and/or prolonged low water conditions). This paper can provide guidance to the ZRA if they choose this option to solve the problem to some extent.

        Key Words: The Kariba Dam; Removing; Replacing; Detailed analysis; New system of dams

        在大蔚氖褂悶詡洌故障是不可避免的。而且,有必要強調的是,卡里巴水壩在南非數千人的日常生活中是非常重要的。因此許多組織都非常注意卡里巴大壩的問題。2008年2月6日,BBC報告說,大雨可能會導致水從大壩中被釋放出來,這將迫使下游的50 000居民進行疏散。上升的水平線導致了2010年3月閘門的開閘,這需要居住在洪泛區的13萬人疏散,并且洪水可能蔓延到附近地區引起了廣泛的關注[1]。

        為了方案的順利實施,我們考慮了盡可能多的方面。從新壩的選址到對周圍環境的影響,從建設成本到這些小水壩的收益,從新壩建設對當地人民的風險評估到所需動力和水量的最低要求等。下面將通過設置適當的模型并將模型與實際結合給出細節分析,并借用虛擬數據來體現我們在計算機程序中的想法。我們最終目標的關鍵點是達到合理性和可行性。

        1 詳細分析過程的概述

        首先,有幾個解決卡里巴大壩故障問題的相關要點:卡里巴湖的水容量調節能力;居民財產和生命安全;新壩系統的防洪能力;優化新壩的位置;得到的所有方案的比較;變量和一般規律;結果和事實之間的關聯程度。

        然后,基于這些想法,我們確認了我們工作的順序。

        (1)根據假設找到所需數據。

        (2)使用已有數據和我們的想象力建立模型。

        (3)驗證模型并在要點的基礎上改進缺點。

        (4)總結所有所得,給出一個或多個解決方案。

        2 假設

        在建模后,通^調整變量,找出隨著資金總數的增加,水庫蓄水能力達到最優時所對應的建壩位置、規模、數量之間隱含的變化關系。從所有可行方案中選取成本合適、對水資源調節能力較強的一組,通過查找數據模擬河流以年為周期的水流變化情況,制定合理的調水方案,在保證水庫安全的情況下盡量少放水、均勻放水,同時保證人們用水。

        3 模型描述

        3.1 一級模型

        3.1.1 理想河

        根據當地部分沿河城市的海拔數據,我們做出了一定范圍內河流海拔高度與距離源頭距離的擬合曲線,將整條河流模擬成了理想化的形式。

        河道為平滑的梯形:為了使河流寬度貼近真實數據,在保證河道寬度的最大值、最小值和平均值的情況下,通過隨機數模擬河道寬度;假定在不存在大壩干擾河流水量的情況下,河流單位面積的水流量相同,且水面與河道底部平行。

        3.1.2 理想壩

        理想壩是一個建在梯形河道上的大壩,其在正常蓄水位下便可以把水完全攔截,不必擔心水從河道中溢出。正常蓄水位以下均為水壩可調節庫容。

        3.2 二級模型

        在理想河模型上我們共獲得了45個等距的可以建造水壩的點。我們可以在此基礎上獲取每個可選壩址的信息并通過下面提及的公式與函數大致算出不同大小、不同位置的大壩的性能及成本。隨后,我們利用退火算法得到最優的解決方案,即在同一工程造價下收益最多的方案。

        3.3 三級模型

        在我們找到新壩群合適的位置和規模后,原有水量將根據水庫現有水量和規模分配給每個水庫。在電腦上模擬徑流是通過基于多年最大徑流量,最小徑流量和平均徑流量的隨機數法。為計算各壩的月流入量、月流出量和現有存儲量并給出適當的分配和運輸水的方法,我們為水庫設置了最高和最低安全水位。

        4 模型建立

        4.1 一級模型

        根據我們的理解與假設,我們首先要能夠得到壩址以解決選址問題。所以我們的初級任務是基于數據模擬河流情況。該模型建立后,探索各變量之間的聯系以得到最優解。

        根據公式計算庫容:

        V=BLH/K (1)

        式中:V為水庫總庫容;B為壩址的河道寬;L為回水長度;H為正常蓄水量;K為常數。

        根據公式計算回水長度:

        L=H/C (2)

        式中:C為河道平均坡度。

        根據公式計算平均坡度:

        C=(d1-d2)/S (3)

        式中:d1、d2為海拔高度;S為1和2間的河道長度。

        4.2 二級模型

        二級模型,見表1。

        4.3 三級模型

        (1)使用搜索到的信息限制釋放水的體積;(2)建立線性規劃的數學模型;(3)通過隨機數法生成合適的引水方案。

        5 模型的檢測和敏感性分析

        5.1 參數的探討

        一些參數是在討論和查找信息后合理創造的。通過查找大量的材料我們盡可能多的使用了已有數據,在確實無法獲得實際值的情況下才會自己設置參數。我們花費了大量的時間尋找我們建模的數據,因為我們相信更多的實際數據,意味著更準確和符合現實的預測。

        5.1.1 理想河的參數

        用矩形代替梯形河道的行為魯棒性較差并且根據自然河道的情況,很容易判斷矩形河道是不合理的、無法使用的。

        5.1.2 二級模型的參數

        正常蓄水量和成本都有著低魯棒性。這是兩個非常重要的參數,因為它們與水壩性能和總成本相關。

        5.1.3 三級模型的參數

        (1) 對干旱和洪水的標準進行調整;(2)增加隨機數的波動范圍來影響河流徑流量。被增強的水庫調節能力達到了足以抵抗洪水的程度,也就是等于擁有了優越的魯棒性。

        5.2 規律

        通過對3組不同造價的建造方案在組內進行比較和進行組間比較之后得出以下結論。

        (1)隨著造價的提升,相應的最大蓄水量前期變化明顯,后期上升較為緩慢,在造價約為2.96億時,性價比較高,且對水的調節能力有顯著的提高。

        (2)在所討論的適宜造價范圍內,造價相同時,建較多的大壩性價比較低,通過對性價比較高的方案的分析,推薦建壩數量為12~14個。

        (3)通過對每個可建壩點的建壩次數的累加,可以優化適合建壩的位置,在理想模型下,考慮到我們需要大蓄水量,因此推薦選擇河道較為寬闊的地方建壩,且在這樣的位置建立水壩的規模也往往較大。

        (4)通過數據發現,大壩如果建得過于密集,其蓄水能力將大打折扣,因此在實際選擇建壩位置時,不宜過密。

        6 總結

        6.1 建模結果

        建模結果見圖1、圖2、表2、表3。

        6.2 模型優勢

        (1)通過線性規劃得到解決方案,如果在這個問題中我們可以得到更準確的數據,那么方案將更加準確,這就意味著此算法有改進的空間。

        (2)退火算法和線性規劃的方法使解決方案更具有普遍性。類似的問題可以通過僅改變少數參數解決。

        (3)通過隨機數模擬的方法,可以使模型多變,這便能幫助我們找到一般規律,這說明此模型具有較強的解決問題的能力。

        (4)所用模型都很易建立。

        6.3 模型劣勢

        (1)盡管模型已得到改進,它仍然與現實不同。

        (2)信息的缺乏使一些數據不夠準確并會影響最終結果的準確性。

        6.4 未來工作的展望

        (1)對卡里巴水庫大壩的實地考察以獲得更多實際數據。

        (2)使用新數據優化方案。

        (3)更新不準確的數據。

        (4)在情況變得更糟糕前盡快執行最佳方案。

        7 結語

        該文通過數學建模的思想給出了建造大壩隨著資金增加的變化傾向于一般規律和調節水量的一般方法。經過討論、查閱資料合理制定了一些參數,提升預測的準確性。該文給出的方案可以在安全和成本方面為卡里巴新水壩的水流量系統提供合理的策略。解決已知或預測的正常水周期,解釋和證明了正確處理應急水流情況應當采取的行動。為日后卡里巴的修繕工程提供一定的幫助。

        參考文獻

        [1] Kariba Dam[EB/OL].https:///wiki/Kariba_Dam#cite_note-21.

        第3篇:數學建模魯棒性分析范文

        本文作者:工作單位:安徽埃夫特智能裝備有限公司

        從控制系統設計角度來說,可以采用辯證法內外因基本原理來分析影響重載機器人控制品質的因素,首先,如果系統存在動力學耦合、柔性等非線性因素,僅僅采用傳統的線性控制很難獲得良好的控制品質,底層伺服回路的控制缺陷是影響機器人控制品質的內因。第二,如果運動規劃環節處理不當,傳輸給底層運動控制回路的運動指令不合理,即存在位置不連續,速度不連續,加速度躍變等情況,對系統會產生嚴重的沖擊,即便底層伺服控制設計再優秀,同樣也會嚴重影響系統控制品質,這就是所謂的外因。下面就從內外因角度對目前在機器人運動規劃和底層伺服控制方面的相關進展進行綜述。機器人運動規劃方法運動規劃與軌跡規劃是指根據一定規則和邊界條件產生一些離散的運動指令作為機器人伺服回路的輸入指令。運動規劃的輸入是工作空間中若干預設點或其他運動學和動力學的約束條件;運動規劃的輸出為一組離散的位置、速度和加速度序列。運動規劃算法設計過程中主要需要考慮以下三個問題:(1)規劃空間的選?。和ǔG闆r下,機器人軌跡規劃是在全局操作空間內進行的,因為在全局操作空間內,對運動過程的軌跡規劃、避障及幾何約束描述更為直觀。然而在一些情況下,通過運動學逆解,運動規劃會轉換到關節空間內完成。在關節空間內進行運動規劃優點如下:a.關節空間內規劃可以避免機構運動奇異點及自由度冗余所帶來種種問題[1-4];b.機器人系統控制量是各軸電機驅動力矩,用于調節各軸驅動力矩的軸伺服算法設計通常情況也是在關節空間內的,因此更容易將兩者結合起來進行統一考慮[5,6];c.關節空間運動規劃可以避免全局操作空間運動規劃帶來的每一個指令更新周期內進行運動規劃和運動學正逆計算帶來的計算量,因為如果指令更新周期較短,將會對CPU產生較大的計算負荷。(2)基礎函數光滑性保證:至少需要位置指令C2和速度指令C1連續,從而保證加速度信號連續。不充分光滑的運動指令會由于機械系統柔性激起諧振,這點對高速重載工業機器人更為明顯。在產生諧振的同時,軌跡跟蹤誤差會大幅度增加,諧振和沖擊也會加速機器人驅動部件的磨損甚至損壞[7]。針對這一問題,相關學者引入高次多項式或以高次多項式為基礎的樣條函數進行軌跡規劃,其中Boryga利用多項式多根的特性,分別采用5次、7次和9次多項式對加速度進行規劃,表達式中僅含有一個獨立參數,通過運動約束條件,最終確定參數值,并比較了各自性能[8]。Gasparetto采用五次B樣條作為規劃基礎函數,并將整個運動過程中加速度平方的積分作為目標函數進行優化,以確保運動指令足夠光滑[9]。劉松國基于B樣條曲線,在關節空間內提出了一種考慮運動約束的運動規劃算法,將運動學約束轉化為樣條曲線控制頂點約束,可保證角度、角速度和角加速度連續,起始點和終止點角速度和角加速度可以任意配置[10]。陳偉華則在Cartesian空間內分別采用三次均勻B樣條,三次非均勻B樣條,三次非均勻有理B樣條進行運動規劃[11]。(3)運動規劃中最優化問題:目前常用的目標函數主要為運行時間、運行能耗和加速度。其中關于運行時間最優的問題,較為經典是Kang和Mckay提出的考慮系統動力學模型以及電機驅動力矩上限的時間最優運動規劃算法,然而該算法加速度不連續,因此對于機器人來說力矩指令也是不連續的,即加速度為無窮大,對于真實的電驅伺服系統來說,這是無法實現的,會對系統產生較大沖擊,大幅度降低系統的跟蹤精度,對機械本體使用壽命也會產生影響[12]。針對上述問題Constantinescu提出了解決方法,在考慮動力學特性的基礎上,增加對力矩和加速度的約束,并采用可變容差法對優化問題進行求解[13]。除了以時間為優化目標外,其他指標同樣被引入最優運動規劃模型中。Martin采用B函數,以能耗最少為優化目標,并將該問題轉化為離散參數的優化問題,針對數值病態問題,提出了具有遞推格式的計算表達式[14]。Saramago則在考慮能耗最優的同時,將執行時間作為優化目標之一,構成多目標優化函數,最終的優化結果取決于兩個目標的權重系數,且優化結果對于權重系數選擇較為敏感[15]。Korayem則在考慮機器人負載能力,關節驅動力矩上限和彈性變形基礎上,同時以在整個運行過程中的位置波動,速度波動和能耗為目標,給出了一種最優運動規劃方法[6],然而該方法在求解時,收斂域較小,收斂性較差,計算量較大。

        考慮部件柔性的機器人控制算法機器人系統剛度是影響動態性能指標重要因素。一般情況下,電氣部分的系統剛度要遠遠大于機械部分。雖然重載工業機器人相對于輕型臂來說,其部件剛度已顯著增大,但對整體質量的要求不會像輕型臂那么高,而柔性環節仍然不可忽略,原因有以下兩點:(1)在重載情況下,如果要確保機器人具有足夠的剛度,必然會增加機器人部件質量。同時要達到高速高加速度要求,對驅動元件功率就會有很高的要求,實際中往往是不可實現(受電機的功率和成本限制)。(2)即使驅動元件功率能夠達到要求,機械本體質量加大會導致等效負載與電機慣量比很大,這樣就對關節剛度有較高的要求,而機器人關節剛度是有上限的(主要由減速器剛度決定)。因此這種情況下不管是開鏈串聯機構還是閉鏈機構都會體現出明顯的關節柔性[16,17],在重載搬運機器人中十分明顯。針對柔性部件帶來的系統控制復雜性問題,傳統的線性控制將難以滿足控制要求[17-19],目前主要采用非線性控制方法,可以分成以下幾大類:(1)基于奇異攝動理論的模型降階與復合控制首先針對于柔性關節控制問題,美國伊利諾伊大學香檳分校著名控制論學者MarkW.Spong教授于1987年正式提出和建立柔性關節的模型和奇異攝動降階方法。對于柔性關節的控制策略絕大多數都是在Spong模型基礎上發展起來的。由于模型的階數高,無法直接用于控制系統設計,針對這個問題,相關學者對系統模型進行了降階。Spong首先將奇異攝動理論引入了柔性關節控制,將系統分成了慢速系統和邊界層系統[20],該方法為后續的研究奠定了基礎。Wilson等人對柔性關節降階后所得的慢速系統采用了PD控制律,將快速邊界層系統近似為二階系統,對其阻尼進行控制,使其快速穩定[21]。針對慢速系統中的未建模非線性誤差,Amjadi采用模糊控制完成了對非線性環節的學習[22]。彭濟華在對邊界層系統提供足夠阻尼的同時,將神經網絡引入慢速系統控制,有效的克服了參數未知和不確定性問題。連桿柔性會導致系統動力學方程階數較高,Siciliano和Book將奇異攝動方法引入柔性連桿動力學方程的降階,其基本思想與將奇異攝動引入柔性關節系統動力學方程一致,都將柔性變形產生的振動視為暫態的快速系統,將名義剛體運動視為準靜態的慢速系統,然后分別對兩個系統進行復合控制,并應用于單柔性連桿的控制中[23]。英國Sheffield大學A.S.Morris教授領導的課題組在柔性關節奇異攝動和復合控制方面開展了持續的研究。在2002年利用Lagrange方程和假設模態以及Spong關節模型建立柔性關節和柔性連桿的耦合模型,并對奇異攝動理論降階后的慢速和快速子系統分別采用計算力矩控制和二次型最優控制[24]。2003年在解決柔性關節機器人軌跡跟蹤控制時,針對慢速系統參數不確定問題引入RBF神經網絡代替原有的計算力矩控制[25].隨后2006年在文獻[24]所得算法和子系統模型的基礎上,針對整個系統穩定性和魯棒性要求,在邊界層采用Hinf控制,在慢速系統采用神經網絡算法,并給出了系統的穩定性分析[26]。隨著相關研究的開展,有些學者開始在奇異攝動理論與復合控制的基礎上作出相應改進。由于奇異攝動的數學復雜性和計算量問題,Spong和Ghorbel提出用積分流形代替奇異攝動[27]。針對奇異攝動模型需要關節高剛度假設,在關節柔度較大的情況下,劉業超等人提出一種剛度補償算法,拓展了奇異攝動理論的適用范圍[28]。(2)狀態反饋和自適應控制在采用奇異攝動理論進行分析時,常常要同時引入自適應控制律來完成對未知或不精確參數的處理,而采用積分流形的方式最大的缺點也在于參數的不確定性,同樣需要結合自適應控制律[29,30]。因此在考慮柔性環節的機器人高動態性能控制要求下,自適應控制律的引入具有一定的必要性。目前對于柔性關節機器人自適應控制主要思路如下:首先根據Spong模型,機器人系統階數為4,然后通過相應的降階方法獲得一個二階的剛體模型子系統,而目前的大多數柔性關節自適應控制律主要針對的便是二階的剛體子系統中參數不確定性。Spong等人提出了將自適應控制律引入柔性關節控制,其基于柔性關節動力學奇異攝動方程,對降階剛體模型采用了自適應控制律,主要采用的是經典的Slotine-Li自適應控制律[31],并通過與Cambridge大學Daniel之間互相糾正和修改,確立一套較為完善的基于奇異攝動模型的柔性關節自適應控制方法[32-34]。(3)輸入整形控制輸入整形最原始的思想來自于利用PosicastControl提出的時滯濾波器,其基本思想可以概括為在原有控制系統中引入一個前饋單元,包含一系列不同幅值和時滯的脈沖序列。將期望的系統輸入和脈沖序列進行卷積,產生一個整形的輸入來驅動系統。最原始的輸入整形方法要求系統是線性的,并且方法魯棒性較差,因此其使用受到限制。直到二十世紀九十年初由MIT的Signer博士大幅度提高該方法魯棒性,并正式將該方法命名為輸入整形法后[35],才逐漸為人們重視,并在柔性機器人和柔性結構控制方面取得了一系列不錯的控制效果[36-39]。輸入整形技術在處理柔性機器人控制時,可以統一考慮關節柔性和連桿柔性。對于柔性機器人的點對點控制問題,要求快速消除殘余振蕩,使機器人快速精確定位。

        這類問題對于輸入整形控制來說是較容易實現的,但由于機器人柔性環節較多,呈現出多個系統模態,因此必須解決多模態輸入整形問題。相關學者對多模態系統的輸入整形進行了深入研究。多模態系統的輸入整形設計方法一般有:a)級聯法:為每個模態設計相應的濾波器,然后將所有模態的時滯濾波器進行級聯,組合成一個完整的濾波器,以抑制所有模態的振蕩;b)聯立方程法:直接根據系統的靈敏度曲線建立一系列的約束方程,通過求解方程組來得到濾波器。這兩種方法對系統的兩種模態誤差均有很好的魯棒性。級聯法設計簡單,且對高模態的不敏感性比聯立方程法要好;聯立方程法比較直接,濾波器包含的脈沖個數少,減少了運行時間。對于多模態輸入整形控制Singer博士提出了一種高效的輸入整形方法,其基本思想為:首先在靈敏度曲線上選擇一些滿足殘留振蕩最大幅值的頻段,在這些特定的頻帶中分別選擇一些采樣頻率,計算其殘留振蕩;然后將各頻率段的殘留振蕩與期望振蕩值的差平方后累加求和,構成目標函數,求取保證目標函數最小的輸入整形序列。將頻率選擇轉化為優化問題,對于多模態系統,則在每個模態處分別選擇頻率采樣點和不同的阻尼系數,再按上述方法求解[40]。SungsooRhim和WayneBook在2004年針對多模態振動問題提出了一種新的時延整形濾波器,并以控制對象柔性模態為變量的函數形式給出了要消除殘余振動所需最基本條件。同時指出當濾波器項數滿足基本條件時,濾波器的時延可以任意設定,消除任何給定范圍內的任意多個柔性振動模態產生的殘余振動,為輸入整形控制器實現自適應提供了理論基礎[41],同時針對原有輸入整形所通常處理的點對點控制問題進行了有益補充,M.C.Reynolds和P.H.Meckl等人將輸入整形應用于關節空間的軌跡控制,提出了一種時間和輸入能量最優的軌跡控制方法[42]。(4)不基于模型的軟計算智能控制針對含有柔性關節機器人動力學系統的復雜性和無法精確建模,神經網絡等智能計算方法更多地被引入用于對機器人動力學模型進行近似。Ge等人利用高斯徑向函數神經網絡完成柔性關節機器人系統的反饋線性化,仿真結果表明相比于傳統的基于模型的反饋線性化控制,采用該方法系統動態跟蹤性能較好,對于參數不確定性和動力學模型的變化魯棒性較強,但是整個算法所用的神經網絡由于所需節點較多,計算量較大,并且需要全狀態反饋,狀態反饋量獲取存在一定困難[43]。孫富春等人對于只具有關節傳感器的機器人系統在輸出反饋控制的基礎上引入神經網絡,用于逼近機器人模型,克服無法精確建模的非線性環節帶來的影響,從而提高機器人系統的動態跟蹤性能[44]。A.S.Morris針對整個柔性機器人動力學模型提出了相應的模糊控制器,并用GA算法對控制器參數進行了優化,之后在模糊控制器的基礎上,綜合了神經網絡的逼近功能對剛柔耦合運動進行了補償[45]。除采用神經網絡外,模糊控制也在柔性機器人控制中得以應用。具有代表性的研究成果有V.G.Moudgal設計了一種具有參數自學習能力的柔性連桿模糊控制器,對系統進行了穩定性分析,并與常規的模糊控制策略進行了實驗比較[46]。Lin和F.L.Lewis等人在利用奇異攝動方法基礎上引入模糊控制器,對所得的快速子系統和慢速子系統分別進行模糊控制[4748]。快速子系統的模糊控制器采用最優控制方法使柔性系統的振動快速消退,慢速子系統的模糊控制器完成名義軌跡的追蹤,并對單柔性梁進行了實驗研究。Trabia和Shi提出將關節轉角和末端振動變形分別設計模糊控制器進行控制,由于對每個子系統只有一個控制目標,所以模糊規則相對簡單,最后將兩個控制器的輸出進行合成,完成復合控制,其思想與奇異攝動方法下進行復合控制類似[49]。隨后又對該算法進行改進,同樣采用分布式結構,通過對輸出變量重要性進行評估,得出關節和末端點的速度量要比位置量更為重要,因此將模糊控制器分成兩部分,分別對速度和位置進行控制,并利用NelderandMeadSimplex搜索方法對隸屬度函數進行更新[50]。采用基于軟計算的智能控制方法相對于基于模型的控制方法具有很多優勢,特別是可以與傳統控制方法相結合,完成對傳統方法無法精確建模的非線性環節進行逼近,但是目前這些方法的研究絕大部分還處于仿真階段,或在較簡單的機器人(如單自由度或兩自由度機器人)進行相關實驗研究。其應用和工程實現受限的主要原因在于計算量大,但隨著處理器計算能力的提高,這些方法還有廣泛的應用前景。

        第4篇:數學建模魯棒性分析范文

        關鍵詞:三維人臉建模;人臉檢測;特征提取

        中圖分類號:TP18 文獻標識碼:A文章編號:1009-3044(2010)17-4755-02

        Research on the Thoughts of Design of 3D Face Detection Based on Image

        ZHU Jun-jun, CHENG Tao

        (College of Computer Science, Sichuan normal University, Chengdu 610101, China)

        Abstract: Automatic face detection has long been an active research area because of it's potential for applications such as law enforcement, security access, and man-machine interaction. Most of the extensive recognition methods are based on 2D methods. So, it is not the real 3D methods. This paper focused on the research of the technique of 3D face diction based on image. It means, detect 3D face with two original photos of human's head and other faces. In this paper, I find a method witch is the real 3D method. It is based on the image itself. What we need to do is about three steps to obtain the detection. Three steps: processing, modeling, detecting.

        Key words: 3D face modeling; face detection; feature extraction

        廣義的人臉識別實際包括構建人臉識別系統的一系列相關技術,包括人臉圖像采集、人臉定位、人臉識別預處理、身份確認以及身份查找等;而狹義的人臉識別特指通過人臉進行身份確認或者身份查找的技術或系統?;谌四樚卣鬟M行身份驗證是目前比較流行的身份驗證技術,它具有、友好、 直接等優秀的特點,能夠被廣大用戶所接受。因此計算機人臉識別技術是生物特征最為活躍最有挑戰性的領域之一。它結合了認知科學、圖象處理、計算機圖形學、機器視覺和模式識別等多個研究領域,研究的成果有著廣闊的應用前景。

        1 目前人臉識別中存在的問題

        到目前為止,人臉識別研究雖然取得很大進步,形成了各種方法理論,但仍然存在著各種各校的局限,比如識別的結果容易受到面部表情、人物姿勢、自然光照、天氣狀況、人物生理特征變化等各方面綜合因素的影響。

        主要體現在:1) 人臉人為變形(如表情等)的不確定性;2) 人臉外觀的復雜性(如發型、斑點、膚色等);3) 圖像生成結果過程中的不確定性(如太陽光照的強度、光照來源方向等);4) 人臉的膚色生理復雜性(如婦女生了小孩后臉型發胖、長斑等)。

        2 三維人臉識別國內外研究狀況

        三維人臉識別一起以來有著突飛猛進的發展,自從上個世紀以來,真正標志著三維人臉的識別是1980年開始,國內的三維人臉識別的展開要稍微慢一些,還處于剛起步的狀態。在國內要實現三維人臉識別,首先要對三維數據信息的獲得,然后已經成功地運用在目前的生產實際應用之中,例如我們國內的大型醫院里面用到的采用PHILIPS公司的儀器的CT醫學成像,三維激光切割技術,結構光學方法,但是應用不是不夠成熟,技術依賴太強。而在國外,三維人臉識別技術已經有了大的進展,特別是在三維數據提取過程中有著非常大的發展,比如在對外界光照情況下如何克服不同光源影響,Pentland等提出采用多視角圖像解決姿態變化的人臉識別問題,Georghiade等提出亮度錐方法處理姿態光照變化問題,Blanz等提出3D形變模型方法。國外對于處理三維人臉識別的基本方法是采取數學方面的幾何特征對圖像自身進取處理,利用深度圖象處理技術,分析面貌曲面的曲率等幾何特征,對面貌曲面進行凹凸區域的分割、正側面輪廓邊緣的提取。

        3 三維人臉建模與識別介紹

        三維人臉建模作為計算機視覺和計算機圖形學領域中的一個主要基礎,它是達到識別的先決條件,而人臉憑借其特有的普遍性和易用性成為眾多專家學者進行三維建模的研究對象。從20世紀70年代Parke建立第一個臉部模型開始,許多研究人員致力于三維人臉建模的研究。特別是20世紀90年代以來,計算機視覺和計算機圖形學技術的迅猛發展為三維人臉建模提供了技術支持,學多學者開始嘗試建立逼真的三維人臉模型,并取得了顯著的成績。

        近年來,隨著三維感知技術和三維幾何學學科了解的逐步深入,人們進行人臉識別也開始受三維人臉建模技術的影響而將眼光轉移到關注三維以及二維和三維混合領域上。使用三維人臉模型進行人臉識別研究具有很多優勢,首先,基于3D模型的人臉識別將對姿態和光照變化具有更強的魯棒性。其次,引入3D模型,識別工作就可以在2.5D或者2D深度圖像和3D人臉模型之間展開,從而一定程度上提高識別效率和精確度。雖然事實證明3D感知設備如3D深度掃描儀輔助的3D人臉建模方法在3D人臉識別方面具有非常優秀的效果,但其受可用性和感知設備的高昂費用所限,目前單純利用3D感知設置進行人臉識別的商業應用仍然很少。

        4 三維人臉識別過程所需要的基本步驟

        綜上所述,為了達到一個自動的人臉識別的過程,如圖1所示,需要的步驟有:

        1) 人臉檢測與人臉分割。從自然界給定的場景里面檢測出人臉,然后進行定位,最后提取。

        2) 人臉的規范化。首先要校正人臉在位置、光照和角度等外界環境因素影響下的變化。

        3) 人臉表征。采用某種方法表示出數據庫中的已知人臉和檢測出的人臉,通常的方法有幾何特征、代數特征、特征臉、固定特征模板等。

        4) 人臉識別。根據人臉的表征方法,選擇適當的匹配策略將得到的人臉與數據庫中的已知人臉相比較。

        5 三維人臉識別系統的構成

        5.1 三維人臉的定位系統構成

        人臉定位方法基本上可以分成兩類,分別是對于彩色圖像的和基于灰度圖像的,基于彩色圖像的定位方法有:例如基于膚色的定位方法,該方法利用膚色與背景色彩的區別從人臉圖像中定位出人臉的大致位置,利用模板匹配從候選區域中檢測并提取出人臉;基于灰度圖像的方法:基于人臉輪廓的人臉定位方法,該方法利用對人臉的先驗知識,在snake曲線提取邊緣的算法的基礎上,研究了ASM算法,這種算法先對人臉樣本進行訓練,然后進行曲線檢測,對人臉檢測定位更具通用性。

        5.2 三維人臉的重建的方法介紹

        三維重構人臉的方法有基于三維結構光,基于3D形變模型和基于通用模型等3D形變模型是一種參數化的通用模型,本文就是在3D形變模型標準模型的基礎上重構特定人臉3D形狀。在通用模型上重建出特定的人臉的過程就是參數調整的過程。把表述人臉的特征點分成不同的優先級,利用ASM算法思想,完成人臉中不同優先級的特征點的調整,從而在3D形模型的基礎上重構出特定的人臉。

        5.3 三維人臉識別系統的方法介紹

        1) 主成份分析方法(Principal Component Analysis,PCA):通過正交變換得到從高維降維到低維子空間的變換矩陣,這些正交矩陣基又被稱作“特征臉”。該方法是一種簡單實用的基于代數特征的人臉識別方法,保留了人臉面部器官之間的拓撲關系,也保留了各器官部件的信息,得到廣泛的應用。但它對光照、視角的變化影響比較大,因此對預處理的要求比較高。

        2) 線性判別分析方法(Linear Discriminant Analysis, LDA):該方法是模式識別領域里非常有效的降維手段,它利用訓練樣本集的類別信息,定義了類內擴散矩陣和類間擴散矩陣,策略就是使類內擴散矩陣盡量小,類間擴散矩陣盡量大,以此達到降維、分類的目的。近年來,衍生出許多改進的或擴展的LDA算法,如多類問題LDA方法、Foley-Sammon最佳鑒別矢量集以及UODV算法等。

        3) Gabo小波提取特征法:該方法是模式識別領域里應用較為廣泛的特征提取方法之一,它充分利用圖像中的不同方向和頻率上的特征信息,有效提取能夠表述圖像的細節信息,以此來表征圖像,達到分類的目的。線性分析方法在將高維矢量映射為低維矢量時,其變換矩陣是線性的。然而,人臉圖像實際上非常復雜,很難用線性模型對其進行完全的刻畫和表示,因此,相應的非線性方法得到了很好的研究。常見的非線性方法有:基于核的主成份分析方法(KPCA-Kernel PCA),其主要思想是基于某選定的核函數,構造從輸人空間到特征空間的一個非線性映射。

        具體的流程圖如圖2所示。

        6 三維人臉檢測系統設計的難點與總結

        三維人臉識別中面部的關鍵特征點定位是一個具有相當挑戰性并且很復雜的問題,也是計算機視覺和圖形學領域的一個最基本問題。目前,自動三維人臉識別系統仍然存在著很多的局限和不足。主要存在的問題有:人臉多種姿態、人臉表情異同、環境光照強弱等。因此特征點的精確定位是解決這些問題的基礎。本文知識對人臉特征點定位進行了一些嘗試性的研究,但是還是有很多要待下一步的研究。人臉圖像上的關鍵特征點是臉部主要器官輪廓上的點,這些器官輪廓在人臉部區域具有明顯的邊緣特征。因此,在目標的搜索過程中加入了邊緣約束局部灰度模型,使人臉特征點收斂到邊緣性強的人臉輪廓上,并且在搜索的過程中引入了多分辨率的搜索策略,提高了人臉特征點定位速度和精度。在本文中,由于使用的IMM標準人臉庫,光照比較均勻且庫中的人臉圖像大多是正面人臉圖像,因此目前的算法對光照和人臉的姿態適應能力較弱為了解決光照對算法的影響。在下一步的工作中,首先人臉圖像的光照強度進行校正,將其校正到一個標準的光照條件下,其次要建立多姿態人臉的統計模型。在進行人臉特征點定位之前,先對人臉圖像的姿態進行判定,然后選用不同的統計模型對人臉圖像進行特征點的定位;但是光照的校正和人臉姿態的判定是比較復雜的問題,因此將算法擴展到多姿態的人臉面部的關鍵特征點的自動定位上,并使之具備良好的光照適應能力將是下一步研究的重點。IMM人臉庫中的人臉圖像都是無遮擋的,因此在后續的工作中要對部分遮擋人臉進行特征點的定位。本文中的人臉定位點方法在人臉的表情夸張的情況下的定位不是很正確,一方面因為目前所用的人臉訓練樣本圖像不夠豐富,另一方面是因為特征點處的局部紋理特征模型不夠精細。因此,后續工作主要是豐富人臉訓練樣本集,使訓練所得模型更加具有普遍性。另一方面嘗試對人臉關鍵特征點用較復雜的局部模型進行建模,以便能夠得到更精確的關鍵特征點定位結果。在基于人臉特征點定位的人臉三維重建過程中,為了減少運算量,我們使用的是比較簡單的Candide_3模型,且沒有進行紋理的映射,在后續的工作中,將采用更加復雜的通用三維模型,從而重構處真實感人臉。在人臉識別中,僅僅使用了特征點集之間的線段Hausdorff進行人臉的識別,在后續的工作中可以加入特征點處的局部Gabon特征來進行人臉的識別。

        參考文獻:

        [1] Bowyer K W, Chang K, Flynn P. A survey of approaches and challenges in 3D and multi-model 2D+3D face recognition[J].Computer Vision and Image Understanding,2006,101(1):1215.

        [2] 陳倩.多生物特征融合身份識別研究[D].浙江大學,2007.

        [3] Zhuang Lian-sheng, Wang Yong, Li Bin. A robust facial feature points location algorithm under variable illumination[C]//Proceedings of the IEEE International Conference on Mechatronics and Automation, Luoyang,2006(1):1227-1231.

        [4] 范玉華,馬建偉.ASM及其改進的人臉面部特征定位算法[J].計算機輔助設圖形學學報,2007,19(11):1411-1415.

        第5篇:數學建模魯棒性分析范文

        關鍵詞:穩健性設計;質量控制;田口質量

        中圖分類號:TH122 文獻標識碼:A 文章編號:1009-2374(2013)07-

        產品設計是決定產品的第一也是最重要的環節。產品設計帶來的質量問題如果不及時處理,會引起連鎖反應,其解決需要的時間和費用很高。在設計過程中考慮得全面、合理、仔細能夠有效地降低成本,減少質量問題發生。通過穩健性設計不僅能夠提高質量,還能使產品特性對不可控因素的敏感性降低。

        1 穩健性產品設計技術

        1.1 穩健性設計的基本原理

        產品的質量在其生命周期內會被各種因素影響,這些影響具有不確定性,會導致產品的質量特性波動。直接消除干擾因素,雖然可以解決問題但是實現難度過大、成本過高??梢员M量降低干擾因素,使質量與因素之間關聯變弱,對干擾變得不敏感,這就是穩健性設計的原理。

        1.2 穩健性設計典型方法

        穩健性來源于控制理論中的魯棒性,是指變量對因素發生微小差變的不敏感性。如何定量地度量設計的穩健性是穩健性設計的基礎??尚蟹€健性是指產品性能質量在印象因素作用下穩定在所允許的范圍內的能力;敏感穩健性是指產品性能質量在噪聲因素作用下保持穩定的能力。穩健性的指標有質量損失函數、信噪比、質量信息熵等。經過長期研究和應用,穩健性設計的技術取得了很大的進展,出現了多種穩健性設計方法。馬義中通過熵和協方差矩陣的關系,建立多元質量特性的信噪比計算公式來度量產品質量特性的整體波動,為了克服質量特性協方差不能直接反映質量特性的波動關系,利用信息熵概念度量穩健設計中多元質量特性的整體波動。比較常用的如下:

        1.2.1 田口方法。田口方法以正交試驗設計為基礎,將產品的設計分為系統設計、參數設計和容差設計三個階段,最后通過正交試驗設計確定參數值可以到達的最佳水平組合。該方法為穩健性設計提供了理論基礎,但是必須事先確定方案的大致范圍,局限性強,需要進一步研究。

        1.2.2 雙響應面法。雙響應面法可以將噪聲因素和設計變量結合,綜合考慮其對產品質量的影響。適用于噪聲因素非正太分布,求解誤差小,但是對試驗數據敏感,模型擬合較為困難。

        1.2.3 隨機模型法。隨機模型將設計變量和噪聲的隨機性作為重點考慮對象,將穩健設計引入到有隨機因素的工程中。該方法不要求變量為正太分布,實用價值很高,但是建模和計算較為困難。

        1.2.4 靈敏度法。為了保證產品對設計變量和噪聲因素的影響達到最小,將設計函數、設計變量和噪聲因素的非線性效應作為主要研究對象。該方法可以快速確定需要調節的范圍,可以結合其他穩健性設計方法,但是主要應用于質量指標問題。

        1.3 穩健性設計特點

        1.3.1 穩健性設計可以將噪聲因素分為可控與不可控因素,利用模型盡可能地消除干擾因素,提高質量。

        1.3.2 產品的平均特性可以由內表可查,穩健性由外表考察,采用先進的方法可以在盡可能少的試驗中獲取大量信息。

        1.3.3 穩健性產品設計與工業各個環節一起考慮,將其中的干擾都計算到數學模型中,所以其抗干擾的能力很強。

        1.3.4 穩健性產品設計考慮到系統的非線性特性,可以在有效控制成本的前提下,生產出較高質量的產品。

        1.3.5 該技術需要對技術人員進行專業的培訓,才能迅速轉化為成果。

        1.3.6 已經研發出響應的平臺及軟件,在計算時輔助下大大提升了工作效率。

        1.3.7 適用性強,可以應用到各個領域中,不斷地提高產品質量。

        2 基于田口質量的產品穩健性設計分析

        對于定制生成方式,需要將客戶的定制轉化各種模塊化的配合,然后通過選擇模塊來實現產品的選擇。將田口質量法與模塊化結合,可以減少產品由于質量波動造成的損失。

        2.1 定制產品的多元損失函數

        假設產品的質量特性為y,目標值為,合理偏差為。質量特性的類型可以分為望目、望大、望小三種。望目特性:希望特性值存在一個目標值,并希望實際的特性值圍繞目標波動,且波動越小越好。望小特性:希望特性和波動越小越好,如零件磨損等。望大特性:希望特性和波動越大越好,如機器效率、零件壽命等。產品質量損失函數為:

        (1)

        式中:反映了質量波動程度,k為比例常數,k的值與望目、望大、望小三種特性相關。該公式說明當質量特性在合理的偏差范圍中波動時,雖然產品合格但是仍然會給用戶產生一定程度的損失。

        2.2 定制產品的穩健性表達

        因為噪聲的影響,產品的質量具有一定的隨機性,因此穩健性必須考慮。這里根據產品的質量特性引入了信噪比的概念,通過信噪比來度量產品質量特性的穩健性。

        其中作為產品穩健性的度量指標,其值越大,產品越趨于穩健,損失就越小。因此根據不同的質量特性可以對賦予不同權重。

        2.3 定制產品的穩健選配模型

        根據數學模型,模塊實例選配的步驟如下:

        2.3.1 分析產品質量特點,判斷其類型。

        2.3.2 根據模塊實例的統計數據,計算各種質量特性的信噪比。

        2.3.3 計算質量特性權重。

        2.3.4 確定各質量特性類型,并計算產品質量特性。

        2.3.5 分析產品封閉環,計算配合偏差。

        2.3.6 最優化求解,選擇最佳方案。

        3 實際應用

        對于出頭彈簧參數設置,可以采用穩健性設計原理。已知數據:F=4N,力臂L=24mm,則力矩M=4×24=96N·mm。滿足力矩的參數不是唯一的,不同的彈簧鋼絲直徑、中徑、圈數都會對結果產生不同的作用。這里要選擇最佳信噪比、靈敏度的方案,通過分析找出穩定和調整因素。彈簧設計時質量特性的數學模型:Φ,其中d為直徑,D為中徑,n為有效圈數,Φ為旋轉角度,如表1

        所示:

        所以可以確定最佳參數為:d=1.2mm、D=7.2mm、n=2.9圈、Φ=34°。

        田口方法是質量高、可靠性高、效率高的設計方法,隨著廣泛的應用,效果會越來越好。實踐表明穩健性是產品設計和選擇的第一準則,田口方法可以很好地實現產品的穩健性。

        4 結語

        穩健性產品設計可以提高產品的設計質量,是一種快速、經濟的獲取提高產品質量的可行性解決方案。因此研究穩健性產品設計,并將其各個領域全面推廣是非常必要的。

        參考文獻

        [1] 黃自興,韓壽祖,邵丹雄,等.穩健性產品設計探討[J].石油商技,2002,(3):1-5.

        [2] 陳立周.穩健設計[M].北京:機械工業出版社,2000.

        第6篇:數學建模魯棒性分析范文

        關鍵詞:國際貿易;非線性經濟學;經濟預測;數據分析;評述

        國際貿易是世界各國經濟、政治、文化交流的重要紐帶,通過它可以擴大相互作用,促進相互的經濟合作,改善國際環境,為本國的經濟發展創造良好的外部條件。

        國際貿易是一個復雜時變非線性系統。對國際貿易問題的分析與預測是建立良性經貿環境的前提。國際貿易不僅受到一國(或地區)的經濟條件、自然條件、貿易政策等國內因素影響,還受到國際市場需求、全球經濟增長等國際因素影響,且各影響因素相互之間存在非線性關系。對進出口貿易數據分析與預測,其實質就是建立適當的數學模型并進行應用。研究表明,各個國家的外貿環境復雜多變,導致在原來數據結構之上建立的簡單線性模型失真,這給進出口貿易的分析與預測帶來了很多困難。如何建立能很好擬合該數據特征的非線性模型,是一種積極的探討,近年來日益受到國際貿易領域研究與應用者的重視。改革開放20多年來,我國外經貿事業發生了巨大變化。特別是加入世界貿易組織后,我國對外開放進入了新的階段,我國已經成為世界第三大貿易國。關于我國外貿的準確分析與預測對于促進經濟的平穩、持續增長,制定有效的外貿政策具有重要意義。本文介紹了目前用于國際貿易分析與預測的一些方法,以期供相關學者與實際工作者借鑒。

        一、國際貿易的非線性特性分析

        奴隸社會國際貿易的發展是與暴力掠奪、海上搶劫、販賣奴隸密切聯系在一起的。受自然經濟的制約,國際貿易在奴隸社會經濟中的地位并不重要。封建社會的國際貿易雖然較奴隸社會有了進一步發展,但由于占主導地位的仍然是自給自足的自然經濟,因而貿易的規模和范圍還是很有限的。國際貿易真正獲得巨大的發展,出現在資本主義生產形成和發展時期。從16世紀至第一次世界大戰,國際貿易發展迅速。全世界經歷過一次由國際貿易帶動的全球化浪潮,各國之間的經濟聯系因國際貿易的快速增加而大大加強[1]。許多國家對外交往范圍之廣、程度之深,如果以出口占GDP的比重來衡量的話,現在仍有一些國家(如日本等國)的國際化程度不能達到當時的水平。第二次世界大戰后結束初期到1973年,是國際貿易迅速發展階段。這一階段國際貿易增長速度之快在歷史上是空前的。這主要是因為:(1)戰后發達資本主義國家經濟的迅速恢復和發展;(2)國際分工和生產國際化的深入和擴大;(3)跨國公司的發展;(4)國際金融貿易組織的建立和經濟一體化;(5)關稅與貿易總協定的作用。從1973年到1985年,由于經濟危機、能源危機、貨幣制度危機、農業危機的爆發,導致國際貿易由迅速發展轉向緩慢發展,甚至停滯的階段。20世紀80年代后半期至今,是國際貿易發展速度總體回升的階段。近幾十年涌動著又一波全球化浪潮,呈現出兩個新的特征:中間產品貿易的大幅度增長(當然也是國際貿易的大幅度增長)和FDI流動的快速增加。以發達國家為例,過去20年,美國出口占GDP的比重大約翻了一番;OECD國家,即使略去集團內貿易不計,出口占GDP的比重也大致增長了一倍。與GDP相比,國際貿易出現了更快速度的非線性式增長。不管是從國際貿易的發展史,還是從各個國家的發展情況來看,國際貿易從宏觀上看,均呈非線性發展的顯著特點。

        此處著重評述一下在微觀上我國學者對非線性特征的認知與分析所做的研究與探討工作。

        黃夢橋和王濤生(2005)[2]結合使用相關維數和最大Lyapunov指數方法,采用我國進出口貿易的月度數據序列作為樣本,說明了國際貿易市場具有非線性和低維混沌特征。

        湛墾華、張永安和馮宗憲(1997)[3]應用自組織理論與方法,對國際市場演變的非線性機制作了較系統的分析,從新的視角闡述了國際市場演變的一般規律及特點。他們認為從自組織理論看,國際市場作為一個遠離平衡的、開放的非線性系統,在國際政治、地緣關系、自然資源、人口分布等構成的國際經貿市場的勢場背景下,是以一定的基核為始點,在不斷漲落的外界環境影響下,以不同的演化方式,組成有特定功能的系統。

        祝寶江(2006)[4]把國際貿易信用系統置于耗散結構狀態下進行了研究。說明在耗散結構狀態下的國際貿易信用系統不停地與經濟社會系統通過實物信用和貨幣信用在內外環境之間進行著物質、能量、信息的交換與流通。系統耗掉了一定的能量,具有微小漲落和非線性動力過程,從而使多基元、多組分、多層次的國際貿易信用非線性流動,使系統各要素之間產生協調動作和相干效應。

        二、國際貿易的非線性模型與方法

        關于國際貿易非線性復雜系統的分析預測的模型與方法目前國內外都有諸多探討,并取得了一定效果。下面按方法分類進行評述。

        (一) 神經網絡模型

        神經網絡是目前應用廣泛的非線性建模與預測方法,它具有較強的非線性映射功能,具魯棒性(robust)和容錯性,在股市預測、證券預測、GDP預測、財務報警等經濟領域內皆有應用。各個領域的研究與應用結果說明,經過恰當的結構設置與模型學習,神經網絡預測精度可優于其他傳統預測方法。

        神經網絡模型在國際貿易問題的分析與預測方面也已有一些結果。李小紅(2007)[5]將BP神經網絡應用于重慶外貿出口額的預測,建立起預測模型并進行實證預測。結果表明,所建立的預測模型具有較高的預測精度,可以作為相關部門制定出口貿易發展目標的參考依據。楊衛和平瑛(2005)[6]利用BP神經網絡建立水產品貿易模型,選取1999年至2003年中國主要出口伙伴國的GDP值、我國水產品總量、出口水產品平均價格、進口水產品平均價格,以其作為參數來考察它們與衡量出口的幾個主要指標: 出口總量、出口總額、各出口方式所占數量及總額等之間的關系。結果表明,利用BP神經網絡預測模型可以比較精確地預測水產品貿易情況,但也有些不足,如收斂不夠快,有時會受局部平坦或局部最小的影響。可考慮結合遺傳算法或其它一些全局收斂速度快的算法,來進一步提高精確度。

        近期更多的研究側重于考慮采用神經網絡和其他方法的結合,通過建立混合模型來提高預測精度。張一、徐山鷹和汪壽陽(2003)[7]分析了協整技術及在此基礎上建立的誤差校正模型在預測中的應用。結合神經元網絡的特點,他們提出了一個基于神經元網絡的非線性誤差糾正預測模型,選取了出口、外匯儲備和外國直接投資數據,建立了兩個在不同數據基礎上的出口預測模型,并對模型的結果進行比較分析,證實兩個模型反映了中國近幾年出口的變化,都是比較可信的,并對兩個模型的結果取平均值作為2003年全年中國外貿出口最終的預測結果,預測了中國2003年出口貿易的大體情況。傅曉旗和謝雯(2005)[8]運用協整分析方法和誤差校正模型及非線性BP神經元算法,采用進口、出口、外匯儲備、匯率等數據建立經濟計量模型,對2006年我國進出口進行預測。

        Lean YU、Shouyang WANG和Kin Keung LAI (2008)[9]提出了一種新穎的非線性集成學習方法,在神經網絡基礎上建立了基于神經網絡的非線性誤差糾正預測模型(EC-VAR),合理地運用計量經濟學方法、整合技術和人工智能(AI)方法,對中國的進出口貿易進行預測。

        眾所周知,目前最常用的BP神經網絡模型誤差反向傳播神經網絡中存在諸多有待解決的問題:參數較多,并且部分參數要憑經驗主觀確定;模型極難建立,需要經過反復多次訓練和實驗;模型擬合精度與泛化能力矛盾問題,也就是說,可能會對訓練樣本集“過度匹配”達到較高的擬合效果,而對于新的輸入樣本卻可能產生與目標值差別較大的輸出。各個專家都在如何改進這些問題做積極探討。向劍偉(2007)[10]在相空間重構的非線性思想基礎上,采用貝葉斯正則化方法,建立一個時滯BP神經網絡模型,提高了BP網絡的泛化能力。他將該模型應用于某電子行業進出口貿易非線性時間序列預測中,結果說明改進的模型具有較好的泛化能力,較準確地擬合了進出口貿易發展的歷史以及趨勢。

        (二) 貝葉斯方法

        貝葉斯方法源于英國學者貝葉斯于1763年在皇家學會學報上發表的論文《論機會學說中一個問題的求解》(An Essay Towards Solving a Problem in the Doctrine of Chances)。20世紀50年代后,隨著統計理論及方法的應用范圍擴大,貝葉斯理論也受到了歡迎,并迅速發展。貝葉斯基本思想是,將人們的經驗知識作為先驗信息和采集到的樣本信息一起結合到實際模型中,即利用先驗信息、樣本信息來進行預測。

        Bayes方法應用廣泛,嚴格地說并不能歸為非線性方法,但鑒于其靈活選擇先驗分布的技術,對處理復雜經濟問題有其獨特優勢,所以本文也把他作為處理國際貿易問題的工具之一來給于評述與推薦。此處主要給出幾個非線性模型的Bayes分析方法。

        Jiang Weijin和Xu Yuhui(2006)[11]在相空間重構思想下,提出了一種時滯BP神經網絡模型,并通過貝葉斯正規化方法,提高了泛化能力,效果良好。祝樹金和賴明勇(2005)[12]針對非線性時間序列預測,根據相空間重構的非線性預報思想,利用G-P算法計算飽和嵌入維,即輸入層節點數,同時結合貝葉斯正則化方法確定隱層節點數,并提高網絡的泛化能力,建立了一類時滯的BP神經網絡預報模型(TDBPNN)。該模型選擇了中國1989年1月至2003年6月進出口貿易的月度數據作為訓練樣本,進行歸一化處理,采用多步預測法,不僅能夠有效地擬合實際數據,而且可以合理地預測實際序列的發展趨勢。前文也提及到向劍偉(2006)[10]采用貝葉斯正則化方法,建立一類新的時滯神經網絡(TDBP)預測模型,選取了一電子外貿企業1989年1月至2003年6月進出口貿易月度數據作為學習樣本,運用多步預測方法預測了該企業2004年進出口貿易發展趨勢。數據結果表明精度較好,能很好地跟蹤原時間序列,具有比現有同類方法更快的響應性能,迭代次數減小,具有較強的魯棒性和泛化能力。

        (三) 灰色系統理論

        灰色系統是指元素(參數)信息不完全、結構信息不完全、關系信息不完全、運行的行為信息不完全的系統?;疑到y理論即灰色系統所做的預測的一種理論框架,它是由我國著名控制專家鄧聚龍教授于20世紀80年代首先提出的?;疑到y理論是一種新的建模方法,它可以克服線性回歸模型的不足之處,對數據量的要求比較低,并且計算簡單,建模精度高,被廣泛應用于各種分析與預測領域。

        在國際貿易方面,一些學者運用灰色系統理論,對進出口額進行分析,預測進出口貿易的發展趨勢。鄒晶和姜志新(2004)[13]建立了基開GM(1,1)的灰色預測模型,選取中國2001年10月份到2003年5月份的外貿出口額作為樣本數據進行預測計算。通過對模型進行應用,將預測值與實際的外貿出口數值進行比較,發現誤差較小,模型基本符合要求。李蘇(2008)[14]采集了我國1990年至2004年的進出口總額,對原始數據經累加生成,經過光滑性、準指數規律性檢驗后,建立GM(1,1)的灰色預測模型,預測了2010、2015、2020年我國的進出口貿易總額。綜合分析,預測結果比較符合我國進出口貿易總額的實際發展趨勢。

        (四) 支持向量機

        支持向量機(Support Vector Machines,SVM)是上世紀90年代中期提出的一種機器學習算法,由于它具有自學習、自調整模型的特點,能對各種混沌系統產生較好的預測效果,使之成為當前機器學習界的研究熱點。目前支持向量回歸模型已應用在金融時間序列預測、產品需求和銷售預測、電力負荷預測和故障診斷等方面,獲得較好效果。

        有學者運用支持向量機方法,對進出口數據進行了分析。肖智和陳婷婷(2006)[15]針對SVM以重慶市外貿出口信息為例,運用SVM方法對重慶摩托車出口進行了實證研究,建立了時序預測模型對其出口總量和發展趨勢做出了較為精確的預測。數據結果顯示,SVM方法在外貿進出口信息分析預測中的運用是可行的,而且SVM模型對樣本量小、波動性強的外貿出口時序具有較高的預測精度,尤其是對其發展趨勢的預測,模型具有較高的擬合度。

        (五) 協整分析、向量自回歸和誤差校正模型

        基于協整關系的誤差校正(VEC)模型是一個含有協整約束的向量自回歸(VAR)模型,被認為是一種精度較高的動態經濟預測模型。該模型在宏觀經濟領域中得到日益廣泛的應用,尤其是在貨幣需求、消費、投資、進出口方面取得了良好的效果。將神經網絡等非線性方法與該模型結合,建立一種非線性混合預測模型,應用于進出口貿易,能得到更精確的預測結果。

        協整分析方法的系統提出和介紹歸功于Engle和Granger(1987)。該方法的基礎思想是:如果兩個或以上的變量的值呈現非平穩性,但它們的某種線性組合卻呈現平穩性,表明變量之間存在某種長期穩定關系,即協整關系。如果非平穩的變量之間存在協整關系,必然可以建立一個誤差校正模型。

        在國際貿易中,很多學者利用了協整分析方法和誤差校正模型,對我國進出口進行預測。徐山鷹和汪壽陽(2006)[8]通過對2005年我國對外貿易形勢和2006年國內外宏觀經濟環境的分析,采用1994年1月份到2005年11月份的月度數據,建立了誤差較正模型,并使用BP神經元網絡方法進行非線性誤差較正,對2006年我國進出口進行預測。程桂云(2007)[16]選取國內生產總值(GDP)、出口額(EX)、進口額(IM)三個變量作為研究對象,選取遼寧省1979-2004年的年度數據作為樣本區間,對各變量數據取對數,然后對LNGDP、LNEX、LNIM三個變量進行平穩性檢驗和協整檢驗,建立向量誤差修正模型(VECM),最后對各變量進行GRANGER因果性檢驗,對遼寧省的對外貿易與經濟增長之間的長短期關系進行了實證分析的檢驗。菲和馬超群(2008)[17]根據1999年1月到2006年6月的月度出口貿易額數據,建立了一種基于BP神經網絡和誤差校正向量自回歸模型的非線性混合預測模型?!糐P+1〗結果顯示,該模型能夠反映經濟系統中各變量的長期均衡關系,同時非線性的協整變量能夠反映出經濟系統其他變量的短期波動對預測變量的影響,適合于經濟變量的預測。

        三、結束語

        對國際貿易非線性經濟系統的分析、研究與應用,目前已有一些成果,并在不斷地完善與發展中。國際貿易的定量分析與預測問題極其復雜,影響因素眾多,不可能建立一種完全統一的模型,應該根據不同的目的,不同的數據特征,應用多種分析方法,去選取、建立合適的模型。而且應該及時根據影響國際貿易的因素的變化來檢驗模型,不斷改進和完善模型及方法,提高分析、評估和預測的效果,為我國相關部門提供正確的決策依據。

        參考文獻

        [1]王懷民.關稅減讓、垂直專業化與國際貿易的非線性增長. 青海社會科學,2006(2):30-32.

        [2]黃夢橋,王濤生.國際貿易系統中的分數維和最大Lyapunov指數[J].經濟數學,2005,22(3):301-306.

        [3]湛墾華,張永安,馮宗憲.國際市場演變的非線性機制分析[J]. 西安交通大學學報,1997,31 (增刊1):5-9.

        [4]祝寶江.耗散結構下國際貿易信用的非線性流動[J]. 國際貿易問題,2006(2): 34-39.

        [5]李小紅.基于BP神經網絡的重慶市外貿出口預測研究[J]. 商場現代化,2007(3):214-215.

        [6]楊衛,平瑛.基于BP模型的中國水產品出口規模預測研究[J]. 中國漁業經濟,2005(4):25-29.

        [7]張一,徐山鷹,汪壽陽.一類基于神經元網絡的誤差糾正模型的應用——2003年度中國出口預測[J]. 預測,2003,22(3):21-26.

        [8]傅曉旗,謝雯,徐山鷹,汪壽陽.2006年中國進出口預測與分析[J]. 管理評論,2006,18(1): 22-25.

        [9]Lean YU,Shouyang WANG and Kin Keung Lai. Forecasting China’s Foreign Trade Volume with a Kernel-based Hybrid Econometric-AI Eensemble Learning Approach[J]. Jrl Syst Sci & Complexity ,2008 (21):1-19.

        [10]向劍偉.一種新的時滯神經網絡非線性時間序列預測方法[J]. 現代電子技術,2007(4):118-119.

        [11]Jiang Weijin and Xu Yuhui. A novel method for nonlinear time series forecasting of time-delay neural network[J]. Wuhan University Journal of Natural Sciences,2006,11(5):1357-1361.

        [12]祝樹金,賴明勇.基于貝葉斯正則化的TDBPNN模型的中國外貿預報中的應用用評估[J]. 中國管理科學,2005,13(1): 1-8.

        [13] 鄒晶,姜志新.基于GM(1.1)的灰色系統模型在我國外貿出口預測中的應用[J]. 國際貿易問題,2004(2):27-29.

        [14]李蘇.基于灰色系統模型的我國進出口貿易總額發展趨勢預測[J]. 商業研究,2008(3):113-115.

        [15]肖智,陳婷婷.基于支持向量機的外貿出口預測[J]. 科技管理研究,2006(7):231-234.

        [16]程桂云.對外貿易與經濟增長:基于遼寧省的實證分析和檢驗[J]. 國際貿易問題,2007(2):62-66.

        [17]王菲菲,馬超群.區域出口貿易非線性混合預測模型及其實證[J]. 統計與決策,2008(1):46-48.

        [18]樊重俊,張堯庭.多元自回歸模型的Bayes分析方法[J]. 工程數學學報,1991(1):143-148.

        第7篇:數學建模魯棒性分析范文

        關鍵字:數據挖掘方法

        中圖分類號: C37 文獻標識碼: A

        數據挖掘,也可以稱為數據庫中的知識發現 (Knowledge Discover Database,KDD),就是從大量的、不完全的、有噪聲的、模糊的、隨機的實際應用數據中,提取隱含在其中的、人們事先不知道的、但又是潛在有用的信息和知識的過程[16]。實際上這是一個模式提取的過程,主要基于人工智能、機器學習、統計學等技術,高度自動化的分析企業原有的數據,做出歸納行的推理,從中挖掘出潛在的模式,預測客戶的行為,幫助企業的決策者調整市場策略,減少風險,做出正確的決策。

        1 數據挖掘的分類

        數據挖掘就是發現隱藏在數據中的模式。它所發現的模式按功能模型一般可分為兩大類:描述型(descriptive)的模式和預測型(predictive)模式[21]。描述型的模式是對當前數據中存在的事實做規范的描述,它所揭示的是當前數據的一般特性;預測型模型則是以時間為關鍵參數,對于時間序列型數據,根據其歷史和當前的值去預測其未來的值。根據模式的特征,預測和描述可以通過下面的任務來完成。

        圖1 數據挖掘模型

        1) 關聯規則(Association Rules)

        關聯模式是數據項之間存在的關聯規則,是在同一事件中出現不同項之間的相關性,例如客戶在一次購買活動中所購買的不同商品之間的關聯性。在數據挖掘領域,對于關聯模式的研究開展得比較深入,人們提出了多種關聯規則挖掘算法,如Apriori,DHP,Partition,Sampling,FP-Growth等算法。這些算法能夠發現數據庫中形如“80%的客戶在一次購買活動中購買X商品的同時也購買Y商品”之類的知識。

        2) 分類分析(Classification)

        分類就是構造一個分類模型,把具有某些特征的數據項映射到某個給定的類別上。這個過程分為兩步:模型的創建和模型的使用。模型的創建是指通過對訓練數據集的學習建立分類模型;模型使用是指使用分類模型對測試數據和新的數據進行分類。其中的訓練數據集是帶有類標號的,也就是在分類之前,要劃分的類別是已經確定的,通常分類模型是以分類規則、決策樹或數學表達式的形式給出的。

        分類模式往往表現為一棵樹,從樹根開始搜索,沿著數據滿足的分支走。走到樹葉時就能確定類別。已有許多數據分類方法,如決策樹方法、統計方法及粗糙集方法等。Metha, Agrawal, Rissanen等人開始研究面向數據庫的分類方法。J. Han等人在他們開發的知識發現系統DBMiner中采用了基于概括的決策樹方法,該方法集成了面向屬性的歸納和決策歸納技術。

        3) 聚類分析(Clustering)

        聚類就是根據數據的屬性對一系列未分類數據進行類別劃分,把一組個體按照相似性分成若干個類或簇,即“物以類聚”。其目的是使類間的數據差別盡能大,類內的數據差別盡可能小,即“最小化類間的相似性,最大化類內的相似性”原則。與分類模式不同的是聚類中要劃分的類別是未知的,它是不依賴于預先定義的類和帶類標號的訓練數據集的非監督學習(unsupervised learning ),無需背景知識,其中類的數量由系統按照某種性能指標自動確定。聚類分析的方法有很多,其中包括系統聚類法、分解法、加入法、動態聚類法、模糊聚類法、運籌方法等。采用不同的聚類方法,對于相同的記錄集合可能有不同劃分結果。

        4) 回歸分析(Regression)

        回歸模式的函數定義與分類模式相似,主要差別在于分類模式采用離散預測值(例如類標號),而回歸模式則采用連續的預測值。它通過具有己知值的變量來預測其他變量的值。在最簡單的情況下,回歸采用的是類似于線性回歸的標準統計技術。但在大多數現實世界中的問題是不能用簡單的線性回歸所能預測的。如商品的銷售量、股票價格、產品合格率等,很難找到簡單有效的方法來預測,因為要完全地描述這些事件的變化需要上百個變量,而且這些變量本身往往都是非線性的。為此學術界提出了很多試圖解決這個問題方法,如邏輯回歸、決策樹和神經網絡等。

        5) 序列模式分析(Sequential)

        序列模式分析和關聯規則分析相似,它是描述基于時間或其他序列的經常發生的規律或趨勢,并對其進行建模。它能發現數據庫中形如“在某一段時間內,75%的顧客購買商品A,接著購買商品B,然后又購買商品C,即序列A-B- C出現的頻度較高”之類的知識。序列模式將關聯模式和時間序列模式結合起來,重點考慮數據之間在時間維上的關聯性。在進行序列模式挖掘時主要有以下幾個問題值得注意:首先是序列的持續時間,也就是某個時間序列的有效時間或者是用戶選擇的一個時間段:其次是時間折疊窗口,即在某一段時間內發生的事件可以被看作是同時發生的;最后是所發現的模式時間間隔。

        6) 偏差分析(Deviation)

        偏差分析是指對差異或極端特例的描述,如聚類劃分外的偏離值。在大多數據挖掘方法中都是將這些偏差信息作為噪聲而丟掉,然而在一些實際應用中,這種罕見的數據可能比正常的數據更有價值。比如網絡的入侵檢測和信用卡的欺詐檢測等。我們可在通過這些異常數據的偏差來分析其中的原因,以便對其采用相應的措施。

        2 數據挖掘的主要方法

        數據挖掘是一門新興的研究領域,其技術基礎是人工智能(Artificial Intelligence )。它借鑒了信息論、數理邏輯、進化計算、神經計算和統計學等理論和算法[22]。在此介紹幾種主流的方法。

        1) 遺傳算法

        遺傳算法是一種基于生物進化過程的組合優化方法,它是生物學和計算機相結合的產物,由美國密西根大學的D.J Holland教授和他的同事們在1975年首次提出的。根據適者生存的原則模擬自然界的生命進化機制,形成當前群體適合的規則組成新群體,以及這些規則的后代。

        基于這些思想,根據遺傳算法的最適合模型,并進一步對數據模型進行優化。

        由于遺傳算法是一種弱算法,具有高效性和靈活性的特點,在數據挖掘中也用于評估其他算法的適應度。

        遺傳算法擅長于數據聚類,通過事件的類比和空間上的類比,可以把大量繁雜的信息數據進行系統化、條理化,從而找出數據之間的內在關系,得出有用的概念和模式。再建立數據模式時,將遺傳算法與神經網絡相結合,可以更好地提高模型的適應性。因此遺傳算法廣泛應用于自動控制、機器學習、模式識別和組合優化等領域。

        2) 神經網絡

        神經網絡(neural network )是由多個神經元按照某種方式相互連接形成,靠網絡狀態對外部輸入信息的動態響應來處理信息,網絡的信息分布式存儲于連接權系數中,使網絡具有很強的容錯性和魯棒性。神經網絡的核心是結構和算法,例如Hopfield網就是以結構見長,而BP (back propagation)網是以算法見長。

        神經網絡和基于符號的傳統技術相比,具有直觀性、并行性和抗噪聲性。目前,已出現了許多網絡模型和學習算法,主要用于分類、優化、模式識別、預測和控制等領域。在數據挖掘領域,主要采用前向神經網絡提取分類規則。

        在數據挖掘領域,將神經網絡用于數據挖掘,重點要解決好以下兩個問題:一是降低訓練時間,二是挖掘結果的可理解性

        3) 統計分析方法

        統計分析方法是利用統計學、概率論的原理對數據庫中的各屬性進行統計分析,從而找出其中的關系和規律。統計分析方法是最基本的數據挖掘方法之一。常用的統計分析方法有

        判別分析法:建立一個或多個判別函數,并確定一個判別標準,然后對未知屬性的對象根據觀測值將其劃分歸為已知類別中的一類。

        因子分析法:用較少的綜合變量來表達多個觀察變量。根據相關性大小把變量分組,使得各組內的變量之間相關性較高,不同組變量的相關性較低。

        相關分析和回歸分析法:相關分析是用相關關系來度量變量間的相關程度?;貧w分析是用數學方程來表示變量間的數量關系,方法有線性回歸和非線性回歸。

        偏最小二乘回歸法:是一種新型的多元統計數據分析方法,它主要研究的是多因變量(結果)對多自變量(原因)的回歸建模,特別當各變量內部高度線性相關時,用偏最小二乘回歸法更有效。另外,偏最小二乘回歸較好地解決了樣本個數少于變量個數等問題。

        在數據挖掘中,統計分析方法適用于分類挖掘和聚類挖掘。

        4) 粗集方法

        粗集(rough set)理論的特點是不需要預先給定某些特征或屬性的數量描述,而是直接從給定問題出發,通過不可分辨關系和不可分辨類確定問題的近似域,從而找出該問題中的內在規律。粗集理論同模糊集、神經網絡、證據理論等其它理論均成為不確定性計算的一個重要分支。

        粗集理論是由波蘭華沙理工大學的Z.Pawlak教授于1982年提出的一種研究不完整、不確定知識和數據的表達、學習及歸納的理論方法。粗集理論采用了上近似集合、下近似集合和邊界來定義粗糙集。

        粗糙集合理論可以用于分類,發現不準確數據或噪聲數據內在的聯系。找出可以描述給定數據集中所有概念的屬性子集是個難題。在給定的現實世界數據中,往往有些類不能被可用的屬性區分,則可以用粗糙集合來近似地定義這些類。根據目前己有的給定問題的知識,將問題的論域進行劃分,然后對劃分后的每一個組成部分確定其對某一概念的支持度,即肯定支持此概念或不支持此概念和模糊概念。上述情況分別用3個近似集合來表示。即將知識定義為對事物的分類能力。這種能力分別由上近似集、下近似集、等價關系等概念來體現。

        5) 決策樹方法

        決策樹((decision tree)是一個類似于流程圖的樹型結構,其中樹的每一個內部節點代表對一個屬性的測試,其分支代表測試的每一個結果:樹的每一個葉子節點代表一個類別。決策樹通過一系列的規則對數據進行分類。目前許多基于規則進行歸納的商用數據挖掘系統都是采用決策樹方法。

        決策樹分類方法的優點:

        1、決策樹方法結構簡單,便于人們理解;

        2、決策樹模型效率高,對訓練及數據量大的情況較為合適;

        3、決策樹方法具有較高的分類精確度;

        4、決策樹可以清晰的顯示哪些字段比較重要。

        建立一棵決策樹可能只要對數據庫進行幾遍掃描之后就能完成,這也意味著需要計算的資源較少,而且可以很容易的處理包含很多預測變量的情況,因此決策樹模型可以建立的很快,并適用于大量的數據處理。常用的算法有CHAID,CART,Quest、C5.0和ID3算法。

        建立決策樹的過程,即樹的生長過程是不斷地把數據進行切分的過程,每次切分對應一個問題,也對應一個節點。對每個切分都要求分成的組之間的差異最大。各種決策樹算法之間的k要區別就是對這個“差異”衡量方式的區別。

        對決策樹的批評常見的是,認為其在為一個節點選擇怎樣進行分割時使用的“貪心”算法。此種算法在決定當前分割時根本不考慮此次選擇會對將來的分割產生什么樣的影響。換句話說,所有的分割都是順序完成的,一個節點完成分割之后不可能以后還有機會回頭考慮此次分割的合理性,每次分割都是依賴于它前面的分割方法,只要第一次分割有一點點不同,那么由此得到的整個決策樹就會完全不同。

        除上述方法外,還有把數據與結果轉化表達成可視化形式的可視化技術、模型方法和歸納學習等方法。

        參考文獻

        第8篇:數學建模魯棒性分析范文

        隨著各行各業對板帶材質量要求的不斷提高,軋輥偏心成為影響產品質量的不容忽視的重要因素。厚度控制過程中的軋輥偏心控制技術的開發和研究仍然是板帶材軋制所面臨的共同課題。我國對軋輥偏心控制問題的研究還不深入,本論文的工作就是試圖在這方面做些努力。

        本文的研究內容是厚度控制過程的軋輥偏心控制技術,著重探索應用重復控制抑制軋輥偏心的控制方法,從頻域和離散域兩個方面提出厚度控制系統的重復控制器的設計方案,并對方案在穩定性、穩態特性、過渡過程特性和魯棒性方面進行理論分析,同時對系統進行仿真研究。本文的主要工作如下:

        ⑴ 給出了冷軋厚度控制的數學模型和軋件硬度波動前饋補償的控制模型;對軋輥偏心進行了系統、深入的研究,給出獲得偏心信號模型的改進快速傅立葉變換的方法;

        ⑵ 針對單輸入單輸出PID厚度控制系統,首先提出了單軋輥偏心擾動重復控制頻域設計方案,在重復控制環節中引入一種補償器,有效提高了系統穩態精度。其次提出了多軋輥偏心擾動重復控制頻域設計方案,提出了軋輥偏心的并行重復控制器結構。

        ⑶ 針對多輸入多輸出厚度、張力控制系統,首先給出了系統控制對象模型,其次提出了單軋輥偏心重復控制頻域設計方案,然后擴展到多軋輥偏心控制系統,并給出了單軋輥偏心擾動和多軋輥偏心擾動時重復控制補償器的設計方法。

        ⑷ 針對重復控制對偏心擾動的基波及其諧波抑制效果較好,而對基波和諧波附近頻率擾動的抑制較差問題,提出了一種魯棒重復控制結構,這種結構對軋輥偏心擾動信號的周期不確定性有較強的魯棒性。

        ⑸ 因工程中普遍采用數字化設計,對于流量AGC、反饋AGC控制結構及流量AGC、反饋AGC、軋件硬度前饋的控制結構分別提出了單軋輥偏心、雙軋輥偏心及多軋輥偏心魯棒數字重復控制器設計方案。這些設計方案能有效地降低補償器階次。

        理論分析和仿真結果證明上述提出的系統設計方案的有效性。

        Application study on roll eccentricity control on cold rolling AGC system based on repetitive control

        Abstract

        As the tolerance requirement for the thickness of steel plate and strip products getting tighter and tighter, the roll eccentricity is becoming more and more important factor affecting the product quality. To precisely control the flat rolled products in cold rolling, an investigation on roll eccentricity becomes essential and such research is lacked in our country so far. In this thesis, attention is focused on control of gauge of steel plate and strip in the presence of periodic disturbances such as the eccentricity.

        Repetitive control system, known to be effective for periodic disturbance, seems to fit naturally with the eccentricity compensation problem. The roll eccentricity control technologies in cold rolling of flat rolled strip based on repetitive control theory are mainly studied in this dissertation. The key innovations of this paper are summarized as follows:

        ⑴ The mathematical models of steel plate and strip gauge control are presented. A control scheme of feed forward compensation for material rigidity is put forward. Through theoretical analysis, the characteristics of roll eccentricity are summarized. An modified Fast Fourier Transform algorithm of acquiring roll eccentricity signal is proposed.

        ⑵ For SISO PID gauge control system, design in the frequency domain based on repetitive controller rejecting single roll eccentricity disturbance is introduced. A compensator is included for the first time in the control scheme. In addition, a control structure of repetitive controllers resisting multi roll eccentricity disturbance is proposed. All the design schemes have been analyzed. Simulations show that proposed schemes are effective.

        ⑶ For MIMO gauge and tension control system, repetitive control for a single roll eccentricity compensation is first presented.The structure is then extended to the case of multiple roll eccentricity. The design method of compensator of repetitive controllers is introduced. Theoretical analysis and simulation results are presented to demonstrate the effectiveness of the repetitive control structure proposed.

        ⑷ Repetitive control is useful if periodic disturbances act on a control system. Perfect (asymptotic) disturbance rejection can be achieved if the period is known exactly. For those cases where the roll eccentricity period changes, a robust repetitive controller structure is proposed. It uses a robust repetitive control structure in the feedback configuration, so that small changes of period do not degrade the disturbance rejection properties. The robust repetitive controller shows good result for rejecting eccentricity.

        ⑸ The digital robust repetitive control schemes compensating single and multiple roll eccentricity are proposed. The design framework can reduce the order of compensator effectively. The gauge adopt the structure of constant volume flow AGC and feedback AGC , avoiding system instability caused by control delay of measurement of height instrument. The material rigidity feedforward is added to the gauge control structure. The theoretical analysis and simulation results on the two gauge control structure show a good performance on the rejection of disturbances such as eccentricity.

        目 錄

        摘 要 1

        Abstract 2

        1 緒論 1

        1.1 問題的提出 1 1.3 重復控制理論研究現狀 12

        1.4 本文的主要工作及各部分內容安排 14

        1.4.1 主要研究內容 14 2 軋輥偏心問題的理論分析和冷軋板板帶厚度控制模型 17

        2.1 軋輥偏心問題的理論分析 17

        2.1.1 輥身和輥徑不同軸的情況 17 2.1.3 偏心信號的采集和處理 20

        2.1.4 應用MMFFT方法的偏心控制方案 31

        2.2 帶鋼厚度控制模型 35

        2.2.1 帶鋼冷軋過程的基本方程 35

        2.2.2 厚度反饋控制模型 38

        2.2.3 前饋控制模型 41

        2.3 本章小結 44

        3 冷軋SISO板厚控制過程中軋輥偏心的重復控制 45

        3.1 單軋輥偏心擾動重復控制系統 45

        3.1.1 厚度控制系統結構及組成 45 3.1.3 重復控制環節的設計 54

        3.1.4 重復控制和魯棒PID控制混合設計 58

        3.2 雙軋輥偏心重復控制系統 61

        3.2.1 雙軋輥偏心重復控制系統的結構及仿真 61

        3.2.2 系統穩定性分析 62

        3.3 多軋輥偏心重復控制系統 65

        3.3.1 多軋輥偏心重復控制系統的結構及仿真 65

        3.3.2 系統穩定性分析 66

        3.4 本章小結 67

        4 MIMO厚度、張力控制系統的軋輥偏心重復控制 69

        4.1 厚度和張力控制系統結構和對象模型 69

        4.1.1 過程控制模型 69 4.1.3 厚度、張力及速度控制系統的解耦 75

        4.1.4 閉環控制系統仿真 78

        4.2 單軋輥偏心擾動重復控制系統 80

        4.2.1 系統結構 80

        4.2.2 系統穩定性分析 80

        4.2.3 系統品質分析 83

        4.2.4 系統魯棒性能分析 84

        4.2.5 重復控制器設計 84

        4.2.6 系統仿真 85

        4.3 多軋輥偏心擾動重復控制系統 86

        4.3.1 系統結構 86

        4.3.2 系統穩定性分析 87

        4.3.3 系統性能分析 88

        4.3.4 系統魯棒性分析 89

        4.3.5 系統設計及仿真 90

        4.4 本章小結 96

        5 周期不確定的軋輥偏心魯棒重復控制系統 97

        5.1 周期不確定軋輥偏心擾動的重復控制原理和結構 97

        5.1.1 常規重復器的結構及其對周期不確定擾動抑制分析 97 5.2 周期不確定單軋輥偏心擾動的魯棒重復控制系統 106

        5.2.1 系統穩定性分析 107

        5.2.2 系統動態性能 108

        5.2.3 系統魯棒性分析 108

        5.2.4 系統仿真 109

        5.3 多周期偏心擾動的魯棒重復控制系統 110

        5.3.1 系統的結構及穩定性 110

        5.3.2 雙軋輥偏心擾動的魯棒重復控制系統仿真 112

        5.4 本章小結 113

        6 厚度控制過程的軋輥偏心擾動數字魯棒重復控制 114

        6.1 數字重復控制器抑制擾動信號的原理 115

        6.2 單周期(基波)擾動的魯棒數字重復控制系統 116

        6.2.1 系統結構 116

        6.2.2 單周期擾動魯棒重復控制系統的穩定性 119

        6.2.3 單周期擾動魯棒重復控制系統約束條件分析 121

        6.2.4 單軋輥偏心(基波)擾動數字魯棒重復控制系統 124

        6.3 基波及二次諧波擾動魯棒數字重復控制系統 126 6.3.2 數字重復控制器的設計 129

        6.3.3 單軋輥偏心擾動(基波及二次諧波)魯棒數字重復控制系統 129

        6.4 多周期擾動魯棒數字控制系統 131

        6.4.1 系統的結構及魯棒穩定性 131

        6.4.2 雙軋輥偏心擾動魯棒數字重復控制系統 134

        6.5 本章小結 135

        7 總結與展望 136

        7.1 本文的工作總結 136

        7.2 今后研究展望 136

        參 考 文 獻 138

        在學研究成果 145

        致 謝 146

        緒論

        問題的提出 冷軋過程中,影響產品厚度精度的因素很多,但大體可分為兩大類[3~5],即軋件工藝參數的變化和軋機狀態的變化。軋件工藝參數的變化,主要包括材料的變形抗力和坯料尺寸以及張力、工藝等軋制工作條件的變化。板帶材的化學成分和組織的不均勻、焊接時的焊縫等都會造成材料變形抗力的變化,在冷軋時引起出口厚度的波動。熱軋鋼卷(來料)帶來的擾動主要有熱軋帶厚不勻,這是由于熱軋設定模型及AGC控制不良造成的,來料厚度不均勻將使實際壓下量產生波動,導致軋制壓力和彈跳的變化,進而影響產品厚度精度;熱軋卷硬度不勻(變形阻力),這是由于熱軋終軋及卷取溫度控制不良造成的。來料厚差將隨著冷軋厚度控制逐架減少。但來料硬度確具有重發性,即硬度較大或較小的該段帶鋼進入每一機架都將產生厚差。冷軋時帶鋼前后張力的變化、軋制速度的變化及摩擦系數波動等也是造成軋出厚度波動的原因。帶鋼軋制過程中的張力變化會改變變形區應力狀態,從而造成軋制壓力的波動和軋出厚度的不均。軋制速度變化主要是通過摩擦系數、軸承油膜厚度來影響軋制壓力和實際輥縫,導致軋出厚度的變化。軋機本身的擾動主要包括不同速度和壓力條件下油膜軸承的油膜厚度將不同(特別是加減速時油膜厚度的變化)、軋輥偏心、軋機各部分熱膨脹、軋輥磨損等。軋輥偏心是高頻擾動,會引起板厚周期性波動,影響產品質量。

        此外還有工藝等其它原因造成的厚差,屬于這類的有:不同軋制乳液以及不同速度條件下軋輥-軋件間軋制摩擦系數的不同(包括加減速時的摩擦系數的波動);全連續冷連軋或酸洗-冷連軋聯合機組在工藝上需要的動態變規格將產生一個楔形過渡段;酸洗焊縫或軋制焊縫通過軋機時造成的厚差。這一類屬于非正常狀態厚差,不是冷軋AGC所能解決的,是不可避免的。

        根據產生帶鋼厚度偏差的不同原因,可采取相應的厚度調節方式和措施來消除或減少它。目前,按其調節方式概括為[6,7]:

        ⑴ 調節壓下量即改變輥縫;

        ⑵ 改變帶鋼在機架前、后張力或一側的張力,即改變軋件塑性曲線的陡度;

        ⑶ 改變軋制速度;

        ⑷ 同時改變軋輥輥縫與帶鋼張力。

        在上述調節方式中,最常用的是調節壓下的厚度控制方法[8~10]。調節壓下量即調節輥縫有兩種不同方式,即:

        ① 電動桿渦輪帶動壓下螺絲轉動使工作輥之間的相對輥縫產生變化來實現帶鋼厚度控制的。由于電機、減速機的慣性很大,電機及傳動系統的啟動、制動時間長,因此,從厚度控制指令發出到軋出預定的帶鋼厚度其控制時間更長。另外,因需大的電機、減速壓下 它是通過電機、減速機、蝸機等機電設備,故軋機成本高,而且維修也不方便; 為了克服諸多因素對板帶材厚度的影響,提高產品的厚度精度,已經開發了和發展了多種厚度控制系統[15~17],如測厚儀反饋AGC、壓力AGC、流量AGC、監控AGC和前饋AGC等。傳統AGC在控制精度方面各有其獨特的特點,在軋機上得到廣泛的應用[18~20]。

        ⑴ 測厚儀反饋AGC

        測厚儀反饋AGC系統是在帶鋼從軋機軋出后,通過軋機出口測厚儀測出實際軋出厚度值,并將其與給定厚度值比較,得出厚度偏差:

        (1.1)

        再通過厚度自動控制裝置將變換為輥縫調節量的控制信號,輸出給壓下或推上機構,以消除厚度偏差。用測厚儀信號進行厚度反饋控制時,由于考慮到軋機機構的限制、測厚儀的維護以及為了防止帶鋼斷裂而損壞測厚儀,測厚儀一般裝設在離直接產生厚度變化的輥縫有一定距離的地方,這就使檢測出的厚度變化量和輥縫控制量不在同一時間發生,所以實際軋出厚度的波動不能得到及時反映。結果整個厚度控制系統的操作都有一定的時間滯后,用下式表示: 式中為滯后時間,為軋制速度,是軋輥中心線到測厚儀的距離。由于存在時間滯后,所以這種測厚儀反饋式厚度自動控制系統很難進行穩定控制。因此目前普遍采用利用彈跳方程對變形區出口厚度進行檢測,然后進行反饋控制。這將大大減少滯后,但由于彈跳方程精度不高,雖然加上油膜厚度補償等措施仍不能保證精度。這正是當前推出流量AGC的原因。安裝了激光測速儀后可精確實測前滑,因而流量方程精度大為提高,用變形區入口及變形區出口流量相等法,根據入口測厚儀及機架前后激光測厚儀可準確確定變形區出口處的實際厚度,因而提高反饋控制的精度。根據流量變形區入出口流量相等:

        (1.3)

        式中: 分為入出口帶鋼寬度; 分為入出口的速度, 分為入出口帶鋼的厚度。一般情況下,入出口寬度變化不大,因而有:

        (1.4)

        從而得到出口厚度:

        (1.5)

        ⑵ 間接測厚反饋AGC

        為了避免直接測厚儀產生的時間滯后,常采用壓力間接測厚反饋AGC系統。即借助于測量某一時刻的軋制壓力和空載輥縫,通過彈跳方程計算出此時刻的軋出厚度,亦即:

        (1.6)

        式中:為軋出厚度,為軋制壓力,為預壓靠值,為空載輥縫,為軋機剛度系數。利用此測得的厚差進行厚度自動控制就可以克服前述的傳遞時間滯后,實現穩定的反饋控制,提高產品厚度精度。然而,在計算帶鋼出口厚度的算式中,是在軋輥軸承處測出的輥縫值,軋輥偏心對實際輥縫的影響在此反映不出來,這就給控制系統帶來了誤差。假定在某一時刻,偏心對輥縫的實際影響為,那么此時的實際輥縫值為,實際造成的厚度厚度波動值應為(假設該時刻沒有其它因素使變化)

        (1.7)

        但由于此時輥縫仍為,所以由計算得出的厚度波動為:

        (1.8)

        因和符號相反,顯然。這樣,就給以作為反饋量的間接測厚AGC系統引入了誤差,造成了間接測厚AGC系統調節質量的降低甚至惡化。即當有偏心存在時,實際板厚減少了,但由于這時軋制力增大,間接測厚AGC系統反而認為板厚增加了,因此控制器就越朝著使板厚減少的方向動作,結果使得比沒有壓力的AGC系統時的板厚精度更為低劣。

        由此可見,間接測厚AGC系統克服了時間延遲,是一種實用、有效的厚度自動控制系統。但是,如前所述,間接測厚AGC系統不但不能對偏心有所抑制,而且還會由于軋輥偏心的存在而導致其控制質量的進一步變差。當產品精度要求較高或軋輥偏心較嚴重時,間接測厚AGC就不可能達到滿意的控制效果。所以,在配置有間接測厚AGC系統的軋機上,常常附加一些抑制偏心影響的措施,如設置死區、帶通濾波等。這些措施避免了軋輥偏心對間接測厚AGC系統的惡劣影響,卻不能消除軋輥偏心對軋出厚度所產生的直接不良影響。

        ⑶ 前饋AGC

        考慮到來料厚差是冷軋帶鋼產生厚差的重要原因之一,因此冷連軋機一般在第一機架前設有測厚儀,可直接量測來料厚差用于前饋控制,機架間亦設有測厚儀用于下一機架的前饋控制。前饋AGC的原理是根據來料厚度波動信號,再根據軋制速度作適當延時,在波動部分進入機架的同時調節輥縫,以消除厚度偏差。輥縫調節量為:

        (1.9)

        式中:為軋件塑性系數。

        ⑷ 張力AGC

        冷軋帶鋼,特別是后面的機架,帶鋼愈來愈硬,越來越薄,因此塑性變形越來越困難,亦即其值越來越大,因而使壓下效率越來越小。

        (1.10)

        式中:為壓下效率,當遠遠大于時,為了消除一個很小的厚差需移動一個很大的。

        采用液壓壓下后由于其動作快使這一點得到補償,但對于較硬的鋼種,軋制較薄的產品時精調AGC還是借助于張力AGC。當然張力AGC有一定的限制,當張力過大時需移動液壓壓下使張力回到極限范圍內以免拉窄甚至拉斷帶鋼。

        ⑸ 監控AGC

        機架后測厚儀雖存在大滯后但其根本優點是高精度測出成品厚度,因此一般作為監控。監控是通過對測厚儀信號的積分,以實測帶鋼厚度與設定值比較求得厚差總的趨勢(偏厚還是偏薄)。有正有負的偶然性厚差是通過積分(或累加)將相互抵消而得不到反映。如總的趨勢偏厚應對機架液壓壓下給出一個監控值,對其“系統厚差”進行糾正,使帶鋼出口厚度平均值更接近設定值。為了克服大滯后,一般調整控制回路的增益以免系統不穩定,或者放慢系統的過渡過程時間使其遠遠大于純滯后時間,為此在積分環節的增益中引入出口速度。其后果是控制效果減弱,厚度精度降低??朔鬁蟮牧硪环N辦法是加大監控控制周期,并使控制周期等于純滯后時間,亦即每次控制后,等到被控的該段帶鋼來到測厚儀下測出上一次控制效果后再對剩余厚差繼續監控,以免控制過頭。這樣做的后果亦將減弱監控的效果。為此,有些系統設計了“預測器”,通過模型預測出每一次監控效果,繼續監控時首先減去“預測”到的效果,使監控系統控制周期可以加快,并且不必擔心控制過頭而減少控制增益。

        即此偏心將使帶鋼出口厚度產生的波動,這一嚴重影響是不容忽視的。不僅如此,如前所述,軋輥偏心還會對壓力AGC系統產生不良的影響,使其調節質量惡化。所以,要想軋出高精度帶鋼,必須考慮補償軋輥偏心影響的措施。采用厚度外環和壓力內環的目的亦是為了抑制偏心的影響。軋輥偏心將明顯反映在軋制壓力信號和測厚儀信號中。對軋制力來說,實測的軋制力信號實際是由給定軋制力(其中包括來料厚度和來料硬度帶來的影響)和偏心信號綜合組成[1],考慮到這兩部分信號在控制策略上是相反的,因此在未投入偏心補償時必須通過信號處理將軋制力信號分解成兩個部分。從軋制力信號提取出的偏心信息可以用下式表示:

        (1.11)

        式中:分是幅值、頻率和初相角。頻率與轉速有關,幅值決定偏心大小,而初相角則決定于信號的初始坐標點,為此需在軋輥上設有單脈沖編碼器(多脈沖等于將軋輥轉角分成多個等分,并以其中一個坐標點作為初始坐標點)。從正弦特性可知,只有兩個幅值相等但反相,頻率相等并且初始角相同的兩個信號相加才能完全互相抵消。否則,頻率不同的正弦信號無法相加;幅值不同則無法完全消除偏心影響;初始角對不準則無法抵消,如果差還可能加劇而不是抵消。由于在實施控制時還要考慮液壓執行機構慣性問題,采用這種兩個完全相反的正弦波抵消的辦法實施起來難度較大。

        軋輥偏心控制技術的研究情況

        軋輥偏心,一般可歸納為兩種類型,一種是由輥身和輥頸不同軸度誤差所引起的偏差,另一種是由輥身橢圓度(不圓度)引起的偏差,由于軋輥偏心的干擾,輥縫偏差一般可達0.025~0.05mm。軋輥轉一周,其干擾變化一次,故軋輥偏心的干擾發生高頻周期變化,從而造成成品帶鋼厚度的波動。軋輥偏心,主要是指支撐輥偏心,因為工作輥直徑小,其偏心量只有幾個;而支撐輥直徑一般為1500mm左右,軋輥磨床加工精度能保證軋輥橢圓度約為,上下輥疊加。隨著用戶對產品質量要求日益嚴格,這種軋輥偏心的干擾越來越不能忽視。為了有效抑制偏心干擾,對系統各個部分的快速性和準確性都要求很高,任何部分的誤差和時滯都會影響補償效果,甚至可能使偏心的不良影響加劇。計算機在工業過程控制中的普遍應用和液壓壓下(推上)裝置在軋機上的應用為解決這個問題提供了硬件上的可能性。由于電動機壓下裝置慣性大,傳輸效率低(一般),對周期性高頻變化無能為力,一般只能在控制系統中設置“死區”,以避免壓下螺絲周期性頻繁動作。而液壓壓下系統慣性小,壓下速度和加速度都顯著提高(一般,同時具有設備重量輕、有過負荷保護能力等優點。對于消除由軋輥偏心所造成的這種高頻變化的周期波動,必須采用這種液壓壓下(推上)系統。 第一類解決辦法按其信號檢測和模型辨識的在線和離線方式,可分為開環控制和閉環控制。按其信號處理手段可分為簡單處理法、各種濾波器法和傅立葉級數法。早期的簡單處理方法包括用千分尺直接測定支持輥的移動或間接測定軋輥軸承座的移動,并根據這個測定值調整安裝在軋輥上的自整角機輸出的正弦波的相位和振幅,按照與支撐輥移動相反的方向實施補償。濾波器方法是一種常用的偏心信號檢測方法。各種濾波器方法都程度不同地存在些問題,不可避免地混進偏心以外的頻率成分,而又毫無辦法地漏掉了偏心信號中的諧波分量。除了濾波器以外,還有解決偏心控制問題的傅立葉分析法。這一方法一般來說要比濾波器方法的信號處理精度高,補償效果顯著。北京科技大學孫一康教授和他的博士研究生劉淑貞在20世紀90年代初以上海第三冷軋帶鋼廠的高精度四輥可逆冷軋機為試驗背景,配以必要的測量儀表和計算機系統,并利用快速傅立葉變換的偏心控制方案,利用相干時間平均方法的偏心控制方案和復合建模偏心控制方案進行大量的現場實驗,取得了滿意的實驗效果[35~37]。

        澳大利亞的E.K.Tech等提出的用于冷軋機的改進的帶鋼厚度控制器和我國原冶金部自動化院陳振宇教授等提出的冷軋機軋輥偏心自校正調節器則應屬于第二類。在消除軋輥偏心影響的同時,也抑制了其它干擾因素對帶鋼厚度均勻性的影響。Tech方案是根據軋制原理,建立一套包括支持軋輥偏心效應、軋機部件的塑性變形過程和彈性變形形變在內的控制設計模型并估計偏心信號周期。反饋控制器對軋制力、滯回、與軋機有關參數和軋制力調整機構的非線形響應進行補償。此方法在把偏心分量從厚度計法厚度誤差估計中分離出來,通過前饋方法補償偏心干擾效應的同時,也實現了準確的厚度估計,通過反饋回路完成了綜合厚度控制。這一方法在澳大利亞公司的冷軋機的初步現場實現表明,它可使軋輥偏心對軋制力和帶鋼出口厚度的影響減少30%,使總的厚度精度提高40%。但此方法要求對軋機系統各部分的機理和參數都了解得很清楚,而且對測厚儀的安裝位置等也有限制,這對有些軋機而言是難以實現的。

        國外對偏心診斷、智能和最優控制的研究較深入和富有成果,主要有:Kugi等提出基于穩定傳遞函數的因數分解逼近和最小均方算法;Aistleitner K等提出采用神經網絡進行偏心辨識的方法;Garcia等提出了采用多處理器實時偏心診斷方法和實時模糊偏心診斷方法;Fechner等提出了神經偏心濾波器,該濾波器用于在線偏心控制時對于變化的偏心周期具有較好的適應性,該方法還用到了遞歸最小二乘學習算法;Choi 等提出了偏心最優控制方法等。

        除此之外,歐美日各大公司的工程專家也提出了多種軋輥偏心的補償方法,這些方法又可以分為下面三類:

        ⑴ 被動軋輥偏心控制方法。這類方法不是試圖補償軋輥偏心對軋件厚度的影響,其主要目的是使輥縫控制系統對軋輥偏心引起的厚度干擾影響不敏感,而不需要輥縫按照輥縫偏心函數進行校正,這就排除了厚度變化增大的可能;

        ⑵ 主動軋輥偏心補償法。這類方法一般包括軋輥偏心分量檢測和隨后得出的補償信號送到輥縫調節器中以補償軋輥偏心,軋輥偏心分量是從反映主要軋制參數(如軋制力、輥縫、軋件出口厚度以及帶鋼張力等)的信號中測得的,根據檢測信號的不同處理方法,這類方法可分為下面兩種:

        ① 分析法 軋輥偏心分量是通過應用數學分析法(例如傅立葉分析法)從檢測信號中提取出來;

        ② 綜合法 軋輥偏心分量是通過復制軋輥偏心分量得到,信號復制可采用機械法和電量法;

        ⑶ 預防軋輥偏心控制法。這類方法是在軋制前創造一些條件以便能減小偏心對厚度的影響,而在軋制中不采用任何校正措施。

        國外公司典型的偏心補償方法有:

        ⑴ 死區法 死區法是一種被動偏心控制法,此法通常可消除控制信號中的周期分量;

        ⑵ 軋制力法 軋制力法是一種主動式軋輥偏心方法,把出口厚度的誤差信號轉換成附加軋制力基準信號;

        ⑶ 輥縫厚度控制法 輥縫厚度控制法是利用安裝在軋機工作輥之間的傳感器測出軋制過程中的輥縫偏差,由德國Krupp提出的輥縫控制(IGC)系統就由輥縫傳感器組成的,它們被裝在機架每側的工作輥輥頸之間,這樣,它們不會受到帶鋼的損壞;

        ⑷ 前饋控制法 已經在軋機輥縫控制中得到廣泛的應用,它包括以下三個步驟:

        ① 在上游機架的前幾機架的軋制道次中,分段測出帶鋼厚度波動;

        ② 當帶鋼每一段即將進入末尾即機架軋制輥縫中時,確定所需的厚度修正量;

        ③ 在末尾幾機架中對帶鋼每一段實施厚度修正。應用這種方法能夠補償包括軋輥偏心在內的各種因素在內的厚度偏差。一般在中間使用張力控制系統主要有兩種。第一種是通過調節上游機架的速度進行帶鋼張力控制,第二是調節下游機架的輥縫進行帶鋼張力控制。成功采用前饋控制系統控制軋輥偏心的關鍵在于軋機電機能否使速度調節器獲得適當的速度響應特性;

        ⑸ Newmann法 這種方法是由德國穆勒-紐曼公司的Newmann等人提出的,它是利用隨支承輥同時旋轉的凸輪來模擬軋輥偏心,位移傳感器測出凸輪偏心,然后發出電子信號,傳送給輥縫調節器。這種方法雖然簡單,但沒有得到廣泛應用。原因是:

        ① 在機架中安裝支承輥之前,顯然要仔細測定每一個支承輥偏心幅度和相移;

        ② 在軋輥偏心測定結束后,每個凸輪和支承輥偏心相移必須一致。由于支承輥偏心明顯非正弦變化,所以要把它和凸輪正弦變化對應起來相當困難;

        ③ 支承輥與凸輪外形的不協調性也是造成軋輥偏心不能得到補償重要原因;

        ④ 不能補償工作輥橢圓度造成的輥縫變化;

        ⑹ Alsop法 以測厚儀原理為基礎進行輥縫控制。假設帶鋼厚度發生波動,使軋制載荷產生低頻波動,而軋輥偏心使載荷產生相當高的頻率波動,載荷信號的低頻分量在任何通道都不會衰減,它將產生正反饋,正反饋大小為: 式中:為軋機縱向剛度,另一方面載荷信號的高頻分量僅能通過一個通道,就這部分來說,載荷回路中產生負反饋信號,增益大小為,這樣回路會產生信號,它被送到輥縫調節器以補償軋輥偏心;

        ⑺ Smith 法 英國戴維聯合儀器公司的Smith提出以測厚儀原理為基礎的輥縫控制系統中軋輥偏心補償法,它的缺點是使用了金屬構件類型的整流器,它會產生于控制信號的波幅差不多的噪聲信號;

        ⑻ Howard法 英國戴維聯合工程公司的Howard提出利用在軋制過程中兩個所測定的參數來測定軋輥偏心,第一個參數是安裝在軋機每側的載荷傳感器測出軋制力的波動量,第二個參數為即將進入軋機的軋件厚度波動量;

        ⑼ Shiozaki(鹽崎)、Takahashi(高橋)法 也稱為軋輥偏心傅立葉分析法(FARE),它是日本的石川島播磨公司(IHI)Shiozaki、Takahashi提出的,該方法應用了軋輥偏心量ec和軋制力變化量之間的關系:

        (1.13)

        式中:Q為軋件塑性系數,為軋機縱向剛度。因為軋輥波動量與支承輥旋轉一周周期一致,于是可得:

        (1.14)

        式中:A為偏心量幅值,為支承輥角位置與軋輥零偏心位置之間的相位角。由于軋制力波動包含有許多不同頻率的分量,對于一級諧波來說,根據簡單傅立葉級數,其變化量表達式為:

        (1.15)

        式中:B、C為常數。在支承輥旋轉一周的時間內,通過測量軋制力的變化量就可以獲得A、B、C和,按照預設定的時間間隔對測定的軋制力進行采樣,其中T是支承輥旋轉一周所需時間,為旋轉一周的采樣個數。可得:

        ,,, (1.16)

        通過FARE法測出偏心信號通過壓力控制回路可以調節輥縫,以便減小或增大偏心補償載荷,偏心補償載荷信號將持續累積到軋輥偏心載荷分量在軋制載荷信號中完全消失為止。然后,當再也測不出偏心載荷分量時,FARE輸出信號就被存儲在存貯器中。隨著軋制持續運行,FARE信號不斷存儲于存儲器中,并且持續計算;

        ⑽ Cook法 西屋電氣公司的Cook提出的方法是建立在假設軋輥偏心所起的軋制力變化為正弦變化,變化周期等于支承輥旋轉周期基礎上。假設,軋制力為

        (1.17)

        式中:為支承輥旋轉一周對應的平均軋制力,為軋制力變化的振幅,為支承輥選定零位與平均軋制力對應的支承輥位置之間的夾角,為軋輥角位置。于是得到:

        (1.18)

        式中:分為上下支承輥對應的補償信號波幅:

        , (1.19)

        式中:為軋機縱向剛度。

        ⑾ Fox法 Cook法的應用局限于雙驅動布置的電機,而檢測軋輥偏心需花費大量的時間,西屋電氣公司的Fox利用上下支承輥之間的差異產生的搖擺現象控制偏心。根據Fox法,在壓靠時將軋輥轉動但不咬入軋件時測定軋制力,假定軋輥偏心變化量呈正弦變化,此時在一個偏擺周期內,兩軋輥軋制力信號分別等于: 式中:分為上下支承輥角位置,分為偏心引起的軋制力波動幅度。

        ⑿ Ichiryu等人的方法 日本日立公司的Ichiryu等人提出提出連續測定入口帶鋼厚度和軋制力,然后使用這些測量值獲得出口厚度,根據相關函數,利用統計方法就可以測出軋輥偏心造成的干擾量,然后從控制系統中消除;

        ⒀ Hayama(葉山)方法 該方法已應用在三菱重工研制的自動軋輥偏心控制系統中,這種方法的原理是使用在線和離線方法檢測軋輥偏心,然后加權求和。離線法是在壓靠條件下利用搖擺現象測定軋輥偏心,在線法是在軋制條件下,通過使一個支承輥相連的脈沖發生器信號和所測的軋制力信號聯系起來,進行軋輥偏心檢測;

        ⒁ Yamagui(山口)法 日本日立和新日鐵公司的山口提出的軋輥偏心方法是通過出口厚度偏差采樣測得的從頭前轉期間的數據計算出軋輥偏心補償信號;

        ⒂ Weihrich和Wohld法 德國西門子公司的Weihrich和Wohld提出的軋輥偏心的方法是基于測厚儀原理,通過求和放大器用輥縫的輸出信號和載荷傳感器輸出信號來計算板帶出口厚度,而求和放大器的輸出信號也包含有軋輥偏心成分。利用輥縫值和成正比的信號就可以通過信號混合器產生軋輥偏心總的信號,同時也改變入口處板帶厚度變化成分和入口處板帶穩定成分,利用高通濾波器,從混合器輸出信號中去掉穩定成分;

        ⒃ Gerber法 伯里斯(bliss)公司的Gerber開發了一套自適應數字化偏心補償(ADEC)系統,該系統利用了聲學技術的最新成果,即具有復制信號中的任意選定交變成分的技術;

        ⒄ Ooi(大井)法 日本住友公司的Ooi利用支承輥平衡液壓缸的這些機構來控制偏心。這種方法是使帶有電動壓下結構的軋機無須進行任何顯著的的改進就可以實現系統高精度的快速效應。支承輥偏心通過傅立葉分析就可以確定出上下輥操作及驅動側位置相關的軋輥偏心成分;

        ⒅ Ginzburg法 國際軋鋼咨詢公司及聯合工程公司的Ginzburg提出兩種軋輥偏心補償方法,第一種方法是利用差拍現象,尤其是利用上下支承輥向同一方向發生偏心時軋輥偏心最小的事實;第二種方法是在軋制過程中對軋輥偏心進行連續補償。

        總之,隨著對該高質量板帶材需求的日益提高,軋輥偏心控制問題得到各國軋鋼控制界的普遍重視,各種檢測和控制方法相繼出現。國外大公司一般在這個領域獲得專利,我國在這個領域尚有差距,需要促進對軋輥偏心控制技術問題的理論分析和研究,不斷使其走向深入和完善。

        重復控制理論研究現狀 針對周期信號發生器正反饋帶來的非平凡問題在如何保證系統穩定問題,Hara等證明,如果對象是正則的且不是嚴格正則的,系統就能保證穩定[48]。為了克服這種重復控制系統不易穩定的局限性,Hara等1988年提出在重復控制環節中引用低通濾波器來濾掉高頻部分,以高頻部分犧牲一些特性來實現系統的魯棒性。因此低通濾波器的選擇對于重復控制非常重要,它的引入一方面有利于系統穩定,另一方面,卻帶來系統的穩態誤差,它反映了閉環系統特性和系統魯棒穩定性間一種折中考慮。1985年和1988年Hara等提出了基于狀態空間的設計方法。近來,魯棒優化控制和結構奇異值方法也用來設計和分析重復控制[49,50]。Peery 和 Ozbay(1993)利用無窮維優化控制原理提出了一種2步法設計優化重復控制器。他們同時提出通過優化重復控制器的濾波器進一步改善系統主要特性的方法。Guvcac(1996)對于連續時間的重復控制系統結構奇異值提出魯棒穩定和動態特性分析方法,即分別用-1和1代替系統內模的延遲部分估計結構奇異值的下確界和上確界,這樣就把原來的無窮維問題化作有窮維問題??梢岳眠@種結構估計連續時間重復控制系統的穩定性和魯棒特性。但是,得出結構奇異值的下確界比用1代替時小,上確界又比用-1代替時大。直到延遲足夠大這種估計才能得到滿意的結果。另外,這種估計還需滿足相位要求,因此這種結構不能用來綜合。

        重復控制器不斷被改進,且被數字化[51~54]。為了減小控制器離散化造成的誤差,很多研究者關注于用離散化方法直接設計重復控制器。Tomizuka等提出一種針對穩定開環對象的零相位偏差跟蹤的重復控制器(Zero Phase Error Tracking Controller),這種方法特點是濾波器的結構和對象同階并滿足時延?;谙嗤难a償器結構,Tsao和Tomizuka(1988,1994)進一步獲得使系統魯棒穩定的內模零相位低通濾波器的方法,給出了和非模型動態的關系,確定了魯棒穩定的充分條件。這種方法可以用于最小相位和非最小相位系統。Alter 和Tsao推導出基于二維模型匹配算法的重復控制算法,并它應用到線性馬達的控制過程。Kim和Tsao(1997)綜合前饋、重復和反饋控制方法,實現電液執行器的魯棒特性控制。Tsao 等把重復控制利用到凸輪機械的非圓旋轉。在極點配置方法中,Ledwich 和Bolton提出了LQ(Linear Quadratic)設計方法。Hillerstrom和Sternby(1994)提出了基于標準Bezout辨識的極點配置方法。Bamich 和 Pearson(1991)提出了采樣數據提升技術(lifting technology)并將其用于設計最優采樣數據重復控制系統。Langari 和Francis(1996)提出基于結構奇異值的采樣數據魯棒控制系統的魯棒分析方法。

        Srinivasan和Shaw提出了頻域設計方法[55,56],并提出了被稱作重構譜的頻率函數[57,58],利用它來判定重復控制系統的相對穩定性。如果在沒有重復控制環節時閉環系統穩定,則對于頻率,是系統穩定的充分條件。設計重復控制的離散時間重構譜的改進方法由Srinivasan和Shaw于1993年提出?;谥C波頻率處對象頻率響應的系統穩定改進方法在1995年由Sadegh提出。 Hanson(1996年)提出一種序貫重復控制系統。首先利用最優控制設計能增加閉環動態硬度的內環控制器,然后基于零相位偏差跟蹤控制設計外環重復控制器以保證跟蹤或抑制周期輸入。由于這是兩步設計(兩個控制器分別設計),所設計的控制器階次必然高。內環最優控制器的特性將在最大峰值2處被外環重復器降低。Guo提出利用替代基于重復控制零相位跟蹤控制中的。選擇和做為靈敏度函數進行頻率調整,以抑制磁盤驅動伺服控制的二次諧波干擾的抑制。眾所周知,基于重復控制的零相位偏差跟蹤控制需要是低通濾波器,且頻帶盡可能寬。因為的選擇必須兼顧重復控制特性和穩定魯棒性,因此靈敏度函數的頻率調整受這種因素限制。Li和Tsao成功應用魯棒重復控制于磁盤伺服控制。

        本文的主要工作及各部分內容安排

        主要研究內容

        由于重復控制對周期性信號具有很好的自學習能力,因此對周期性擾動具有很好的抑制作用。重復控制只需知道擾動信號的周期,對信號的初始狀態,如初相角和幅值等沒有要求,這樣大大簡化信號的檢測,同時降低了控制難度。重復控制的難點在于對系統穩定性要求較高。國內外將重復控制應用于軋輥偏心控制的文獻不多。圍繞研究帶鋼高精度厚度控制的目的,本文以獲得厚度精度控制為目標,重點研究厚度控制過程中應用重復控制抑制軋輥偏心擾動。本文主要做兩方面的工作。首先針對厚度控制過程中軋輥偏心補償問題的特點,將先進的的重復控制理論和自動控制理論有機結合應用到這個問題中來,提出控制方案;其次從穩態精度、穩定性和魯棒性三個方面進行理論分析,對控制方案進行計算機仿真研究。

        ⑴ 首先提出了單軋輥偏心擾動重復控制抑制的單輸入單輸出(SISO)厚度控制系統頻域設計方案,方案中為了彌補重復控制延遲環節前引入濾波器帶來的控制精度問題,提出一種補償器,給出了補償器的設計方法。厚度控制采用測厚儀測厚的反饋AGC控制方案,用Smith預估器補償被控對象滯后,補償后的廣義對象采用常規PID控制。同時還給出了一種將魯棒PID控制器和重復控制設計結合在一起的混合設計方法。其次提出了多軋輥偏心擾動重復控制補償的SISO厚度控制系統頻域設計方案。對系統的穩定性、魯棒性和系統動態品質進行了分析,同時對控制方案進行了仿真研究;

        ⑵ 針對多輸入多輸出厚度、張力控制系統,首先提出了單軋輥偏心重復控制頻域設計方案,然后擴展到多軋輥偏心控制系統。采用逆奈奎斯特方法對被控對象進行解耦。對控制方案進行了理論分析和仿真研究,證明重復控制抑制單周期和多周期偏心擾動的有效性;

        ⑶ 針對重復控制對偏心擾動的基波及其諧波抑制效果較好,而對基波和諧波附近頻率信號擾動的抑制較差,同時軋制過程中因各種原因造成軋輥偏心信號的周期可能波動或者偏心擾動信號不能準確測量或辨識情況,提出了一種魯棒重復控制結構,從理論上證明了這種魯棒重復控制較常規重復控制性能優越,對擾動信號的周期波動不敏感,具有很強的魯棒性。將這種結構用于厚度控制系統,仿真結果證明了這種結構對周期不確定軋輥偏心信號具有很強的抑制能力。

        ⑷ 因工程中普遍采用數字化設計,分別提出了單軋輥偏心、雙軋輥偏心及多軋輥偏心魯棒數字重復控制器設計方案.這種設計能有效地降低補償器階次。厚度控制采用流量AGC和反饋AGC結合的控制結構,避開因測厚儀測厚滯后造成的系統不易穩定的弊端。所有方案都進行了理論分析,同時對所提出的硬度前饋和厚度反饋的控制結構在偏心擾動和硬度擾動下進行了仿真,結果證明這些方案的有效性。

        各部分內容安排

        全文共分7個部分,每部分的具體內容安排如下:

        第一章首先闡述了冷軋板帶厚度控制方法和研究現狀,指出了抑制軋輥偏心擾動在高精度厚度控制過程中的重要性;其次,綜述了國內外軋輥偏心的研究成果及現狀。

        接著全面介紹了重復控制概念的基本內涵、應用的對象和重復控制理論的研究成果;最后給出了本文的主要研究內容。

        第二章首先全面而系統的歸納了軋輥偏心的的特點和性質,在此基礎上,給出了獲取偏心信號的改進傅立葉方法;其次給出了幾種厚度控制模型,提出了硬度波動前饋控制模型。

        第三章首先給出了單軋輥偏心擾動重復控制抑制的單輸入單輸出(SISO)厚度控制系統頻域設計方案;其次提出了多軋輥偏心擾動重復控制補償的SISO厚度控制系統頻域設計方案。對系統的穩定性、魯棒性和系統動態品質進行了分析,同時對控制方案進行了仿真研究。

        第四章針對多輸入多輸出厚度、張力控制系統,首先提出了單軋輥偏心重復控制頻域設計方案,然后擴展到多軋輥偏心控制系統。對控制方案進行了理論分析和仿真研究。

        第五章針對周期不確定軋輥偏心信號,提出了一種魯棒重復控制結構,從理論上證明了這種魯棒重復控制較常規重復控制性能優越,并對其抑制周期不確定軋輥偏心信號進行了仿真。

        第六章分別提出了單軋輥偏心、雙軋輥偏心及多軋輥偏心魯棒數字重復控制器設計方案,給出降低補償器階次的方法。對所有方案都進行了理論分析和計算機仿真。

        第七章對全文工作進行了總結,提出了下一步工作設想。

        軋輥偏心問題的理論分析和冷軋板板帶厚度控制模型

        軋輥偏心問題的理論分析

        廣義上說,軋輥和軋輥軸承形狀的不規則引起輥縫周期性變化稱為軋輥偏心。軋輥偏心會導致軋件厚度周期變化,軋輥的偏心可以歸納為兩種基本類型。一種是由輥身和輥徑的不同軸度引起的偏差所引起的;另一種是由軋輥本身所具有的橢圓度所產生的。而實際情況可能是兩者共同作用的結果。

        輥身和輥徑不同軸的情況

        圖2.1 輥身和輥徑不同軸的情況

        如圖2.1所示,為輥徑的軸心,為輥身的軸心,為輥身的半徑,X為與之間的距離。偏心運動軌跡相當于輥身表面可移動點A繞輥徑軸線轉動,即偏心波形為的軌跡。設支承輥轉動的角速度為,,在三角形中, 由余弦定理可知:

        (2.1)

        設t=0時,=0,=,。由正弦定理得:

        (2.2)

        從而有:

        (2.3)

        因而有軋輥偏心運動軌跡的參數方程為:

        (2.4)

        根據以上參數方程,得軋輥偏心波形如圖2.2所示。

        圖2.2 軋輥偏心波形

        圖2.3 輥身為橢圓時的示意圖

        軋輥具有橢圓度的情況

        如圖2.3 所示,o是軋輥的軸心,是理想輥身的半徑,a 和b 分別是實際橢圓截面的長軸和短軸。實際情況可能不是橢圓。偏心波形為橢圓周上可移動點A與理想圓周的徑向距離的軌跡,r為A到軋輥軸心線的距離。設輥身轉動的角速度為,t=0時,,則有:

        (2.5)

        又由橢圓方程 得 :

        從而

        因此有

        (2.6)

        因而得到軋輥偏心曲線方程為

        (2.7)

        得到的偏心波形類似于圖2.2。

        如果兩個輥的角速度相同,那么合成的偏心信號仍然是同頻率的周波。這是因為周期信號可以分解為一系列的正弦波之和。而兩個同頻率的正弦波之和仍是正弦波。設和為兩個角頻率為的正弦波,其中

        (2.8)

        則合成的波形為

        (2.9)

        式中:

        (2.10)

        (2.11)

        合成波形的振幅發生變化,相位發生偏移,頻率保持不變。軋輥偏心波形一般不是純粹的正弦曲線,而是包括多次諧波的復雜的周期波。它有以下特點:⑴ 周期性 軋輥每轉動一周,偏心信號重復出現一次;⑵ 頻率和幅值不是固定不變的。當軋制速度變化時,其頻率也隨之成比例變化。在軋制過程中,由于軋輥的熱膨脹和磨損,偏心信號的幅值也會發生緩慢變化;⑶ 偏心信號不僅含有多次諧波,而且還含有各種各樣的隨機干擾。

        偏心信號的采集和處理

        軋輥偏心對厚度的影響可以用出口厚度變化的頻譜分析來評估,斯太爾克利用快速傅立葉變換(FFT),從出口厚度數字化信號中分離所有周期分量,并依據所有軋輥轉速和尺寸,能夠辨別出大部分頻譜峰值,通過對頻譜選擇過濾同時結合反變換FFT技術,每個軋輥對出口厚度變化的影響都能測量出來。從上面分析中,我們知道軋輥偏心信號是包括多次諧波的高頻周期波,偏心信號的頻率與軋制速度成正比。在生產過程中,由于隨機噪聲、緩慢變化量等的存在,采集的偏心信號會出現突變、漂移等無規則變化,但總的偏心信息不會突變。軋輥更換以后,它的偏心量就基本上確定了。,并在短時間內不會突變。根據這一特點,在每次換輥以后,在正常軋制狀態下,對軋制壓力信號進行采集,從中提取偏心成分,建立偏心模型。進而對軋輥的偏心進行補償。

        將采集到的軋制力信號進行A/D轉換,然后進行去均值(去掉直流分量)和相干時間平均處理,使噪聲干擾得以減弱或消除,提高信噪比;對預處理后的信號進行快速傅立葉變換(FFT),建立軋輥偏心參數模型。在軋輥上安裝一個光碼盤,以產生兩列脈沖。一列相對軋輥某一固定點,每轉一周發出一個脈沖,此脈沖作為采樣和控制的初始定位信號;另一列是軋輥每轉一周,光碼盤發出128個脈沖數列以進行FFT,建立模型。相干時間平均方法適應于周期信號或重復信號,它將各個周期信號和噪聲信號同時疊加后加以平均,如果噪聲是隨機的,則在疊加過程中會相互抵消,而信號是有規律的,疊加平均后幅值不變。必要條件是噪聲應具有一定隨機性,而信號則具有重復性,且兩者互不相干。

        設混有噪聲的信號為,信號反映系統的某種基本特征。在相同的條件下,具有重復性。噪聲為均值為零,方差為的平穩隨機信號,且、互不相關。對第i個樣本采樣M次,然后做相干平均得:

        (2.12)

        傅立葉變換是在以時間為自變量的信號與以頻率為自變量的頻譜函數之間的變換關系。傅立葉變換可以辨別出或區分出組成任意波形的一些不同頻率的正弦波??焖俑读⑷~變換是建立在離散時間概念上的,它不單純是對離散時間付立葉變換的近似,而是從離散付立葉變換出發,有一整套自成體系的、 離散時間域中的嚴格的基本定理和數學關系。離散付立葉變換能把一個有限長度序列映射成另一個有限長度序列,因而很適合于數字計算機計算。利用離散付立葉變換的一些代數結構,可以實現高速算法,快速付立葉變換能使離散付立葉變換的計算時間成數量級的縮短。快速付立葉變換的出現使付立葉變換已不僅僅是一種理論概念,而且成為一種技術手段。

        ⑴ 離散付立葉變換[65 ,66]

        當用數字計算機對信號進行頻譜分析時,要求信號必須以離散值作為輸入,而計算機輸出所得的頻譜值,自然也是離散的。因此,必須針對各種不同形式信號的具體情況,或者在時域和頻域上同時取樣,或者在時域上取樣,或者在頻域上取樣。信號在時域上取樣導致頻域的周期函數,而在頻域上取樣導致時域的周期函數,最后將使原時間函數和頻率函數都成為周期離散的函數。

        從嚴格的數學意義上講,離散周期序列的付立葉變換是不存在的。但是,如果利用周期函數可能展開為付立葉級數的指數形式并使用沖激序列,則可以把付立葉級數逐項作積分變換,從而在形式上得到付立葉變換對。

        設為一周期連續信號,如果以抽樣間隔為的抽樣率進行抽樣,抽樣結果為,則可表示為:

        (2.13)

        設一個周期內的抽樣點數為,即到,則

        可寫成:

        于是有:

        (2.14)

        對進行抽樣等于先將它的一個周期抽樣成,然后把這一個周期進行延拓。所以有:

        (2.15)

        式中上的符號表示周期重復,它是離散時間周期沖激序列,是的一個周期內抽樣所得的數值;為抽樣序號,;為抽樣間隔;為的周期;為任意整數。

        令,并將展開成付立葉級數

        (2.16)

        式中:,的單位為,系數可表示為:

        (2.17) (2.18)

        對式(2.18)進行付立葉變換得:

        (2.19)

        定義

        (2.20)

        由于

        所以。這里是的個周期,。也就是說的周期為,在每個周期內,。于是,式(2.20)可寫成:

        (2.21)

        上式說明,周期離散時間序列經付立葉變換后在頻域中是離散頻率的周期序列,這種形式的變換也稱為離散付立葉級數變換。在數學上,離散周期序列的付立葉級數變換可簡明表示為:

        (2.22)

        (2.23)

        為了方便,令,則式(2.22)和式(2.23)可表示為:

        (2.24)

        (2.25)

        離散付立葉級數變換是周期序列,仍不便于計算機計算,但離散付立葉級數每個周期序列卻只有(一個周期內取點個數)個獨立的復值,只要知道它的一個周期的內容,其它的內容也就知道了。同時限制式(2.24)中的和式(2.25)中的都只在區間內取值,就得到了一個周期的和一個周期的之間的對應的關系:

        (2.26)

        (2.27)

        這就是有限長序的離散付立葉變換對。

        上兩式所示的離散付立葉變換對可以看成是連續函數在時域、頻域取樣所構成的變換,可以看作是連續付立葉變換的近似,是一種很有用的變換方法。然而,當數據有較長的長度時,這種變換的計算量是很大的。分析式(2.26) 和式(2.27)可知,當用直接方法計算DFT時,總運算量及總運算時間近似地比例于,這在很大時,所需的運算量及總算時間近似地比例于,這在很大時,所需的運算量非??捎^,要想用DFT方法對信號作實量處理一般是有困難的。

        ⑵ 快速付立葉變換(FFT)

        快速付立葉變換是為減少DFT計算次數的一種快速有效的算法。它使DFT的運算大為簡化,運算時間一般可縮短一至二個數量級,其突出的優點在于能夠快速高效地和比較精確地完成DFT的計算。

        FFT改善DFT運算效率的基本途徑是利用DFT中的權函數所固有的兩個特性,一個是的對稱性,即,另一個是的周期性,即。利用的對稱性,可根據正弦和余弦函數的對稱性來歸并DFT中的某些項,結果可使乘法次數約減少一半。假定是一個高復合數,可利用權系數的周期性,把點DFT進行一系列分解和組合,使整個DFT的計算過程變成一個系列迭代運算過程。因為迭代運算的計算量要比直接計算的計算量少很多,尤其是當很大時,可能成百位甚至成千倍地減少??焖俑读⑷~變換算法正是基于這一基本思想而發展起來的。權系數的周期性是導出FFT算法的一個關鍵因素,高復合性則是實現FFT算法的一個重要條件。根據不同的分解方法,可以導出多種FFT算法,如按時間抽取的FFT算法,按頻率抽取的FFT算法,的高復合性則是實現FFT算法的一個重要條件。根據不同的分解方法,可以導出多種FFT算法,如按時間抽取的FFT算法,按頻率抽取的FFT算法,為復合數的FFT算法等。時域抽點算法的迭代過程是基本在每級把輸入時間序列分解為兩個更短的子序列,頻域抽點算法的迭代過程則基于在每級把輸出頻率序列分解成兩個更短的子序列。

        以2為基時域抽點FFT算法是最基本最常用的算法,基2算法要求采樣點數為2的整數次冪。設有一個點序列,而,首先將按序號之奇偶分解為兩個點的子序列,因而得:

        (2.28)

        如采用下列變量替換:(當為偶數時),(當為奇數時),則上式可變為:

        (2.29)

        又因

        所以上式又可改寫為:

        (2.30)

        由于對于均有定義,而及只對有定義,因此,有必要就情況下對2.30作出說明。根據DFT的周期性可得:

        (2.31)

        考慮到:

        則上式可改寫為:

        (2.32)

        經整理后得:

        (2.33)

        式中:和可分別寫成序列和的點DFT。

        式(2.33)表明,一個點DFT可分解成兩個點DFT,而這兩個點DFT又可組合成為一個點DFT,效果是相同的,但是運算量卻大不相同。很明顯,如果以一次復乘和一次復加稱為一次運算,那么,計算兩個點DFT約共需運算,此外再加上按式(2.33)組合需要次運算,所以按先分解后組合的方式計算一個點DFT總共約需次運算。當較大(即)時,它的運算量比直接運算點的DFT約可減少一半。

        因為是2的冪,所以可進一步將每個點子序列按奇偶號分解為兩個點子序列,再令每兩個點子序列組合成一個點DFT……。上述分解過程還可繼續進行,直到第次分解,每個子序列都只有兩點。這樣,就把點DFT的運算轉化為級組合運算,M級組合就是M級迭代過程。每次迭代要求N/2次復乘和N次復加,M級迭代約需次復乘和次復加。每次迭代要求次復乘和點DFT的迭代運算過程是基于在每級把輸入時間序列分解成兩個更短的子序列,因此稱為時域抽點算法。圖2.4 說明了此迭代運算過程。

        圖 2.4 N點基2 FFT的M級迭代過程

        經過FFT變換結果,就可以計算出各次諧波的振幅和相角,從而建立軋輥的偏心模型,其振幅A=,相角,頻率隨軋輥速度變化而變化。

        偏心模型還必須轉換為與采集脈沖對應的離散點的模型,即將帶有三個參數的正弦波偏心模型轉換成128個脈沖對應的離散點模型。軋輥偏心控制對檢測和控制系統的準確性和快速性要求很高,定位定點采樣保證了通過數據處理獲得的偏心模型的唯一性和準確性。把正弦波的一個周期分成N段,列成表格,用步長DELTA掃過這個表,用序號作為角度參數,查表求出序列的值。假設每兩個采樣點之間的時間間隔維t,則正弦頻率為。當步長不是整數時,采用點可能落在兩表值之間,可以采用線性內插法加以修正。

        ⑶ 基2時域FFT算法的改進(MMFFT)

        針對軋輥偏心信號本身及其控制問題的特點,對傳統的基2時域FFT算法進行改進(MMFFT)。改進分兩部,第一步改進的是取消傳統FFT方法對采樣持續時間的限制,使快速付立葉變換算法適用于處理軋輥偏心波動這類周期未知或變動的周期信號,同時又能抑制FFT固有的泄漏效應。第二步改進是就偏心控制問題而言,將周期信號中各次正弦波的絕對頻率轉換為相對頻率,從而提高算法在偏心控制中應用的可靠性和實用性。

        ① 第一步改進(Modlified FFT)

        人們對DFT感興趣主要是因為它是連續付立葉變換的一個近似。近似的準確程度嚴格說來是被分析波形的一個函數,兩個變換之間的差異是因DFT需要對連續時間信號取樣和截斷而產生的。因而在應用DFT解決實際問題時,常常遇到混疊效應、柵欄效應和泄漏效應等問題。

        對一個連續信號x(t)進行數字處理時,要在計算機上進行計算,而計算機的輸入只允許是數字信號,所以必須對連續信號x(t)進行抽樣,即

        (2.34)

        式中:為對x(t)抽樣所形成的序列。T為抽樣間隔,為抽樣率,。如果抽樣率選得過高,即抽樣間隔過小,則一定的時間里抽樣點數過多,造成對計算機存貯量的需要過大和計算時間太長。但如果抽樣率過低,則在DFT運算中將在頻域出現混疊現象,形成頻譜失真,使之不能反映原理的信號。這樣將使進一步的數字處理失去依據,而且也不能從這個失真的頻譜中恢復出信號來。因此,對連續信號的抽樣率需大于奈奎斯特頻率,即抽樣率至少應等于或大于信號所含有的最高頻率的兩倍,即。

        如果x(t)是一個周期信號,它只具有離散頻譜,那么,x(t)抽樣后進行FFT運算得出的頻譜就是它的離散頻譜。但是如果x(t)是個非周期函數,它的頻譜是連續的,把x(t)的抽樣進行DFT運算得到的結果就只能是連續頻譜上的若干點。因為這就好象是從柵欄的一邊通過縫隙觀看另一邊的景象一樣,所以稱這種效應為柵欄效應。如果不附加任何特殊處理,則在兩個離散的變換線之間若有一特別大的頻譜分量,將無法檢測出來。減少柵欄效應的一個方法就是在原記錄末端填加一些零值變動時間周期內的點數,并保持記錄不變。這實質上是人為地改變了周期,從而在保持原有線連續形式不變的情況下,變更了譜線的位置。這樣,原來看不到的頻譜分量就能夠移動到可見的位置上。

        泄漏效應是由于在時域中對信號進行截斷而引起的。實際問題中,所遇到的離散時間序列x(nT)可能是非時限的,而處理這個序時時,需要將其限制為有限的N點,即將它截斷。這就相當于將序列乘以一個矩形窗口,如果對有限帶寬的周期函數抽樣后的截斷長度并不正好是其周期的整數倍,就會導致離散付立葉變換和連續付立葉變換之間出現顯著的差異。這是因為,根據頻域卷積定理,時域中的,則頻域中與進行卷積。這里,和分別是的付立葉變換,這樣將使截斷后的頻譜不同于它加窗以前的頻譜。泄漏效應的產生是由于矩形窗函數的付立葉變換中具有旁瓣亦有一定帶寬而引起的。如圖2.5所示。為了減少泄漏,應盡量尋找頻譜中窗函數,即旁瓣小、主瓣窄的窗函數?;蛘咄ㄟ^限制采樣的持續時間來抑制泄漏效應。

        圖2.5 矩形窗口的時域與頻域圖形

        第9篇:數學建模魯棒性分析范文

        摘 要 本文提出了一種標準粒子濾波器的改進算法——高斯混合采樣粒子濾波算法(gmsppf)。仿真結果表明,新算法在大幅降低計算復雜度的前提下,具有比標準粒子濾波算法(sir-ppf)更好估計性能. 關鍵詞 卡爾曼濾波;粒子濾波;序列蒙特卡洛;貝葉斯濾波;高斯混合采樣 1 引言 貝葉斯方法為動態系統的估計問題提供了一類嚴謹的解決框架。它利用已知的信息建立系統的概率密度函數可以得到對系統狀態估計的最優解。對于線性高斯的估計問題,期望的概率密度函數仍是高斯分布,它的分布特性可用均值和方差來描述。卡爾曼濾波器很好地解決了這類估計問題[1]。對于非線性系統的估計問題,最經典并得到廣泛應用的方法以擴展的卡爾曼濾波為代表,這類方法需要對模型進行線性化,同時要求期望的概率密度函數滿足高斯分布,然而在對實際系統建模時,模型往往是非線性非高斯的。此時,最優估計很難實現。 粒子(particle)濾波器——序列重要性采樣粒子濾波器,是一種適用于強非線性、無高斯約束的基于模擬的統計濾波器[2]。它利用一定數量的粒子來表示隨機變量的后驗概率分布,從而可以近似得到任意函數的數學期望,并且能應用于任意非線性隨機系統。本文介紹一種估計性能更好的粒子濾波算法——高斯混合采樣粒子濾波器(gmsppf),相比通常意義上的粒子濾波算法(sir-pf),gmsppf粒子濾波器具有更小的系統狀態估計的均方誤差和均值。 2 貝葉斯濾波問題 貝葉斯濾波用概率統計的方法從已觀察到的數據中獲得動態狀態空間(dss)模型參數。在dss模型中,包含狀態和觀測兩個方程[3][4]。其中狀態轉移方程(state equation)通常寫作 (1) 這里,是已知,且是白噪聲獨立的隨機序列,而且分布是已知的。觀測方程表達式寫為 (2) 這里:是白噪聲序列,獨立且分布已知。并且滿足。 圖1描述了dss模型中狀態轉移和似然函數的關系。假設初始時刻系統的狀態分布已知,k時刻的已知信息序列表示。 圖1 動態狀態空間模型(dssm) 這樣,貝葉斯估計的問題理解為:利用觀測到的信息yk,求解系統狀態的概率分布。若系統狀態的變化是隱馬爾柯夫過程,即當前系統的狀態信息只與上一個時刻的狀態有關,可以通過預測和更新的途徑求解。 (3) 這里: (4) 假設xk,wk是相互獨立的隨機變量,滿足 。于是,參考(1)式可以把(4)式寫為 (5) 其中,是采樣函數。當是已知時,xk可以通過確定性方程(1)得到。 依據貝葉斯準則,系統狀態估計量 (6) 其中, (7) 另外,在給定 xk,vk,分布的條件下, yk的條件概率依據測量方程(2)可以表示為如下形式 (8) 由(6)式可以看出,后驗概率密度包含3個部分。先驗概率似然函數和證據。如何獲得這三項的近似是貝葉斯濾波的核心問題。更新方程(5)中觀測值 用來對 的先驗預測值修正,從而獲得狀態 的后驗概率。方程(3)和(6)的遞歸關系構成了求解貝葉斯估計問題的兩個步驟:預測與更新。如果(1),(2)中的hk,fk是線性的,且噪聲wk,vk滿足高斯白噪聲,可以把貝葉斯估計問題簡化為卡爾曼分析解。但這類問題僅僅是實際問題中很小的一個部分。對于更多的問題,很難得到分析解。只有通過對問題的近似線性處理(擴展卡爾曼濾波)或其它途徑(蒙特卡洛方法)實現非線性、非高斯問題的解。依據后面分析問題需要,這里重點對蒙特卡洛方法積分進行說明。 3 蒙特卡洛方法 在過去的二十多年,蒙特卡洛方法得到了很大的發展。其優點就是用系列滿足條件的采樣點及其權重來表示后驗概率密度。蒙特卡洛方法采用統計抽樣和估計對數學問題進行求解。按照其用途,可以把蒙特卡洛方法分為三類[5]:蒙特卡洛抽樣、計算、優化。其中,蒙特卡洛抽樣是尋找有效的、方差很小的、用于估計的抽樣方法。蒙特卡洛計算則是設計產生滿足特定要求隨機數的隨機發生器的問題。而蒙特卡洛優化是采用蒙特卡洛思想對實際中的非凸非差分函數優化求解。對于,可以由概率空間p(x)中抽取n個樣本,用近似值作為的解。大數定理證明:收斂于,并且滿足條件。這里,是的方差。不同于確定性的數字計算,蒙特卡洛近似的一個重要特點就是估計的精度獨立于狀態空間的維數。而且,積分估計的方差與采樣點的個數成反比。顯然,蒙特卡洛近似方法的關鍵點有兩個:首先如何由一個樣本空間中抽取n個采樣點,用來表征后驗概率密度。其次就是計算。 重要性抽樣(important sampling)解決了如何借助于已知分布來對實現有效采樣的問題,由marshall 1965年提出。當數據空間十分巨大時,重要性抽樣只對其中“重要”區域進行采樣,節省了計算量。對于高維采樣空間模型,如統計物理學、貝葉斯統計量,這一點尤為重要。重要性抽樣的中心思想是選擇一個覆蓋真實分布p(x)的建議分布q(x)[8]。這樣, (9) 對q(x)作蒙特卡洛抽樣,假設粒子數目為n,有 (10) 其中,稱為重要性權重,再作歸一處理, (11) 是歸一化權重。為了減小估計的方差,選擇的建議性分布q(x)與p(x)盡可能匹配。通常,建議分布q(x)需要一個長的拖尾,這樣可以解決區間之外的干擾。確切的說,匹配的q(x)必須與p(x)f(x)成正比[9]。當q(x)與p(x)不匹配時,w(x(i))是不均勻分布的,在整個遞歸迭代的過程中,存在大量的權值極小的樣本,而這些樣本對估計的貢獻很小。事實上,權值較大的少數樣本決定蒙特卡洛采樣的估計精度。大量時間損耗在這些“無關緊要”的粒子計算上,即所謂的粒子退化現象(degeneracy problem)。目前,標準的粒子濾波器選擇先驗概率(prior)作為建議分布。 對于粒子退化現象,采樣—重要性重采樣方法給出了很好的解決途徑。其基本思想就是通過在兩次重要性采樣之間增加重采樣步驟,消除權值較小的樣本,并對權值較大的樣本復制,降低了計算的復雜度。在o(n)時間復雜度范圍內可以已排序的均勻分布序列作重采樣處理。 對重采樣(resampling)處理,新的采樣結果放在數組,具體的算法用偽碼語言寫為如下的形式: 步驟1:令這里必須注意是隨機變量的累計概率密度序列。 步驟2:初始假設,當, 產生一組序列分布。對一個固定的j,分別用逐一比較,一旦,就可以得到一組新的樣本集合。如此循環直到。需要說明的是,重采樣方法在消除粒子退化問題的 同時,也帶來了其它兩個問題:首先,降低了粒子運算并行執行的可能性;其次,由于權值較大的粒子多次被選擇,粒子的多樣性減少。這種情況尤其在小過程噪聲條件下表現更為明顯[11]。 圖2 sir-pf重要性采樣與重采樣示意圖 4 gmsppf濾波算法 如前所述,利用序列重要性采樣和重采樣的方法,粒子濾波可以有效的遞歸更新后驗概率的分布。但是,由于對粒子未加假設,大量的粒子在處理非線性、非高斯問題時出現了計算的高復雜性問題。另外,由于少數權值較大的粒子反復被選擇,粒子坍塌明顯。文獻[4]提出了在重要性采樣步驟的建議分布的生成階段“搬運”粒子到似然較高區域,可以緩解坍塌,同時提高估計的性能。但是不可避免的是對每一個粒子的后驗概率處理,使得計算的復雜性進一步加劇。鑒于此種情況,這里介紹一種新穎的高斯混合采樣粒子濾波器(gaussian mixture sigma point particle filter,gmsppf)。gmsppf算法利用有限高斯混合模型表征后驗概率分布情況,可以通過基于重要性采樣的加權的后驗粒子,借助于加權的期望最大化算法(weighted expection maximization)替換標準重采樣步驟,降低粒子坍塌效應。

        4.1 基于高斯混合近似的采樣卡爾曼濾波器 根據最優濾波理論,一個概率密度p(x)都可以寫作高斯混合模型(gaussian mixture model)。即,這里,g是高斯分量的個數,是高斯分量的權重,是以向量為均值,以p(g)為協方差矩陣的隨機向量x的高斯分布。 考慮dss狀態轉移方程和觀測方程,假設先驗概率及噪聲密度服從高斯混合模型(gmm)。這樣,預測的先驗概率密度滿足,更新后。 這里,。在此基礎之上,預測的先驗概率和后驗概率對應的均值和方差可以通過采樣卡爾曼濾波器(sigma point kf)計算。 4.2 基于觀測更新的重要性采樣(important sampling) 前已敘及重要性抽樣是一種蒙特卡洛方法,即用一組帶有權值的樣本數據來表征隨機變量的概率密度。利用dss模型的一階馬爾柯夫本質和給定狀態的觀測值依賴性,可以推導遞歸的權值更新方程,這里僅對于給定的粒子而言。在gmsppf算法中,用gmm近似來。作為建議分布。由于包含了最新的樣本數據,使得粒子聚集在高似然區域,一定程度減少了粒子坍塌效應。另外,使用預測的先驗概率平滑權值更新方程中的,這是因為gmsppf算法用gmm表示后驗概率,本次后驗同時又是下一個時間步的先驗概率,gmm模型中高斯核對后驗概率做了平滑處理?;谟^測更新步驟的重要性采樣方法中對粒子不作任何假設,對非線性、非高斯問題具有很強的魯棒性。 4.3 采用加權的em算法做重采樣和gmm還原 基于觀測更新步驟的重要性采樣輸出是一組加權的粒子,在標準的粒子濾波器中,這些粒子必須作重采樣處理丟棄小權值粒子,同時對權值較大的粒子做放大處理。通過這種處理,可以有效的防止粒子集合的方差增加太快。不幸的是,重采樣步驟只對當觀測似然微弱、大量粒子聚集極少數粒子副本情況有效。在gmsppf算法中,采用加權的期望最大(weighted expection maximization)直接得到gmm模型,實現對加權粒子的最大似然擬合,這就相當于對粒子的后驗概率做了平滑,避免了粒子坍塌問題,同時,gmm模型中的高斯核的個數減少到g,防止其呈指數級增長,降低了算法復雜度。 為了比較算法的性能,系統狀態估計的條件均值,均方誤差(error convariane)可以通過兩個方法計算,即在加權的em算法平滑之前,用下面公式

        求解,描述了系統的均值與均方誤差性能。 5 算法性能分析與結論 這里,給定系統狀態估計問題的算法評估模型 (12) 是噪聲,。另外,非平穩觀測模型 (13) ,其中,觀測噪聲服從高斯分布。如果給定含噪的系統狀態觀測值yk,采用兩種不同的算法:標準的粒子濾波算法sir-pf以及gmsppf算法對系統的狀態xk估計。每次實驗共做150次,每次的觀察樣本重新產生,sir-pf算法中粒子的個數是250個。gmsppf算法中采用兩種方案:第一種方案用5個高斯核擬合狀態后驗概率。狀態噪聲vk,觀測噪聲nk各用一個高斯核擬合。第二種方案則用3個高斯核擬合gamma(3,2)分布的拖尾狀態噪聲,這里擬合方法采用em算法。圖3、圖4描述了系統的隱狀態和觀測值及sir-pf,gmsppf算法系統狀態的估計值。 圖3 sir-pf粒子濾波器狀態估計 圖4 gmsppf粒子濾波器狀態估計 采用4.3部分的均方誤差和均值計算公式對不同算法對系統狀態估計性能作了比對。圖3、圖4曲線表明,在系統的觀測噪聲nk均方誤差很小,而過程噪聲服從具有長的拖尾 分布時,采用轉移概率作為建議分布的標準粒子濾波器性能很差。這是因為觀測方程中峰值似然函數和系統狀態急劇的跳躍變化產生的結果。盡管可以通過采樣卡爾曼(sigma-point)濾波器將粒子向似然峰值區域搬動解決這一問題,但是也使得計算量加大。gmsppf算法兩種不同方案都具有比sir-pf更好的系統狀態估計性能,均方誤差比后者數量級降低了1/103-1/104。與1個高斯核擬合過程噪聲的gmsppf算法比較,3個高斯核擬合算法性能更好,但時間復雜度同樣有所提高。 由于gmsppf算法在大幅度降低了算法的計算復雜度同時,可以獲得精確的系統估計性能。所以說,gmsppf算法為粒子濾波理論實時應用,如目標定位(單目標與多目標)、時變信道估計、圖像增強、機器故障診斷以及語音信號處理等提供了一個新的方案。 參考文獻 [1] y.c.ho and r.c.k.lee,”a bayesian approach to problems in stochastic estimation and control”ieee trans.automat.contr. vol .ac-9.pp.333-339 [2]a.doucet,n.freitas,n.gordon. sequential monte carlo methods in practice [m].springer [3]b.d.o. anderson and j.b.moore .optimal filtering . [m]prentice hall englowcod cliffo,nj. 1979 [4]n.j.gordon,d.j.salmond,a.f.m.smith,novel approach to nonlinear/non-gaussian bayesian state estimation,iee proceedings vol140,no2,april 1993 [5]muller,“monte carlo integration in general dynamic models“ contemp .math.1991,115,pp,145-163 [6]fredric gustafsson,niclas bergman,”particle filters for position,navigation and tracking “,final version for ieee transactions on signal processing special issue on monte carlo methods for statistical signal [7]arnaud doucet,simon godsill,”on sequential monte carlo sampling method for bayesian filtering “statistics and computing(2000),10,197-208,recived july 1998 and accepted august 1999 [8] jayesh h.kotecha and petar m.djurric,”gaussian particle filtering ”in proc. workshop statistical signal process.singapore,aug.2001 [9]j.s.liu&r.chen.”sequential monte carlo methods for dynamical systems”.journal of the amerian statisticalassociation,1998,volume 93.pp.1032-1044 [10] zhe chen,”bayesian filtering:from kalman filter to particle filters,and beyond” manuscript in 2003,april [11] jayesh h.kotecha and petar m.djurric,”gaussian sum particle filtering ” ieee transactions on signal processing,2003,vol.51.no.10.october [12]m.sanjeev arulampalam,simon maskell,neil gordon,and tim clapp,”a tutorial on particle filteolinear/non- gaussian bayesian tracking”,ieee transaction on signal processing,vol,50,no 2february 2002

        无码人妻一二三区久久免费_亚洲一区二区国产?变态?另类_国产精品一区免视频播放_日韩乱码人妻无码中文视频
      2. <input id="zdukh"></input>
      3. <b id="zdukh"><bdo id="zdukh"></bdo></b>
          <b id="zdukh"><bdo id="zdukh"></bdo></b>
        1. <i id="zdukh"><bdo id="zdukh"></bdo></i>

          <wbr id="zdukh"><table id="zdukh"></table></wbr>

          1. <input id="zdukh"></input>
            <wbr id="zdukh"><ins id="zdukh"></ins></wbr>
            <sub id="zdukh"></sub>
            性刺激欧美三级在线现看中文 | 日韩精品一区二区三区影院 | 中文字幕午夜乱理片 | 亚洲精品视频久久久久久 | 亚洲五月综合缴情在线 | 日韩AV大片在线一区二区 |