模擬建模分析(2015)課件

1、模擬建模模擬建模模擬的概念及作用 現(xiàn)實系統(tǒng)的數(shù)學或邏輯模型可能十分復雜，現(xiàn)實系統(tǒng)的數(shù)學或邏輯模型可能十分復雜，例如大多數(shù)例如大多數(shù)具有隨機因素具有隨機因素的復雜系統(tǒng)，其中的復雜系統(tǒng)，其中的一些隨機性因素很難用準確的數(shù)學公式表的一些隨機性因素很難用準確的數(shù)學公式表述，從而也無法對整個系統(tǒng)采用解析法求解。述，從而也無法對整個系統(tǒng)采用解析法求解。模擬是處理這類實際問題的有力工具。模擬是處理這類實際問題的有力工具。應用領域：應用領域：運輸系統(tǒng)模擬運輸系統(tǒng)模擬摩天大樓安全疏散系統(tǒng)模擬摩天大樓安全疏散系統(tǒng)模擬國民經(jīng)濟發(fā)展模擬國民經(jīng)濟發(fā)展模擬人口增長系統(tǒng)模擬人口增長系統(tǒng)模擬供水系統(tǒng)模擬供水系統(tǒng)模擬蒙特卡洛(

2、MonteCarlo)法，或稱統(tǒng)計試驗法、計算機隨機模擬方法 蒙特卡洛法蒙特卡洛法是一種基于“隨機數(shù)”，采用統(tǒng)計抽樣方法，近似求解數(shù)學問題或物理問題的過程。 當所求問題的解是某個事件的概率，或者是某個隨機變量的數(shù)學期望，或者是與概率、數(shù)學期望有關的量時，通過某種試驗的方法，得出該事件發(fā)生的頻率，或者該隨機變量若干個具體觀察值的算術平均值，通過它得到問題的解。這就是蒙特卡羅方法的基本思想。 計算機模擬試驗過程，就是將試驗過程（如投針，射擊）化為數(shù)學問題，在計算機上實現(xiàn)。蒙特卡羅方法常以一個“概率模型”為基礎，按照它所描述的過程，使用由已知分布抽樣的方法，得到部分試驗結果的觀察值，求得問題的近似解

3、。 從這個意義上講，蒙特卡羅方法可以部分代替物理實驗，甚至可以得到物理實驗難以得到的結果。用蒙特卡羅方法解決實際問題，可以直接從實際問題本直接從實際問題本身出發(fā)，而不從方程或數(shù)學表達式出發(fā)。身出發(fā)，而不從方程或數(shù)學表達式出發(fā)。它有直觀、形象的特點。 對于那些需要計算多個方案的問題，使用蒙特卡羅方法有時不需要像常規(guī)方法那樣逐個計算，而可以同時計算所有的方案，其全部計算量幾乎與計算一個方案的計算量相當。 另外，使用蒙特卡羅方法還可以同時得到若干個所求量，而不像常規(guī)方法那樣，需要逐一計算所求量。 計算機模擬計算機模擬：在已經(jīng)建立的數(shù)學、邏輯模型的基礎之上，通過：在已經(jīng)建立的數(shù)學、邏輯模型的基礎之上，

4、通過計算機試驗，對一個系統(tǒng)按照一定的決策原則或作業(yè)規(guī)則，由計算機試驗，對一個系統(tǒng)按照一定的決策原則或作業(yè)規(guī)則，由一個狀態(tài)變換為另一個狀態(tài)的行為進行描述和分析。一個狀態(tài)變換為另一個狀態(tài)的行為進行描述和分析。系統(tǒng)狀態(tài)隨時間而變化的動態(tài)寫照系統(tǒng)狀態(tài)隨時間而變化的動態(tài)寫照模擬的作用 對于很難用解析方法加以處理的問題，對于很難用解析方法加以處理的問題，模擬是一種有效的技術；模擬是一種有效的技術； 對建模過程中的假設進行鑒定，對理論對建模過程中的假設進行鑒定，對理論研究的結論加以檢驗；研究的結論加以檢驗； 對不同的實現(xiàn)方案進行多次模擬，按照對不同的實現(xiàn)方案進行多次模擬，按照既定的目標函數(shù)對不同方案進行比較

5、，既定的目標函數(shù)對不同方案進行比較，從中選擇最優(yōu)方案。從中選擇最優(yōu)方案。模擬的分類 通常，通常，模擬時間是模擬的主要自變量模擬時間是模擬的主要自變量。其他變。其他變量作為因變量來看待。量作為因變量來看待。 推進模擬時間的基本方法：推進模擬時間的基本方法： 下次事件法下次事件法：將模擬時間由一個事件發(fā)生的時：將模擬時間由一個事件發(fā)生的時間點推進到緊接著的下一次事件發(fā)生的時間點。間點推進到緊接著的下一次事件發(fā)生的時間點。-系統(tǒng)的狀態(tài)僅在系統(tǒng)的狀態(tài)僅在 固定時間步長法固定時間步長法：模擬時間每次均以相等的固：模擬時間每次均以相等的固定步長向前推進，每到達一個新的模擬時間點定步長向前推進，每到達一個新

6、的模擬時間點需檢查相應時間段內(nèi)是否發(fā)生了事件。需根據(jù)需檢查相應時間段內(nèi)是否發(fā)生了事件。需根據(jù)實際問題合理設置模擬時間發(fā)生改變的步長。實際問題合理設置模擬時間發(fā)生改變的步長。系統(tǒng)狀態(tài)隨時間而變化的動態(tài)寫照系統(tǒng)狀態(tài)隨時間而變化的動態(tài)寫照模擬的分類 根據(jù)模擬過程中因變量的變化情況進行分類：根據(jù)模擬過程中因變量的變化情況進行分類： 第一類：離散型模擬第一類：離散型模擬：因變量在與事件時間有：因變量在與事件時間有關的具體模擬時間點呈離散性變化。大多數(shù)系關的具體模擬時間點呈離散性變化。大多數(shù)系統(tǒng)（如排隊服務系統(tǒng)）可采用離散型模擬。統(tǒng)（如排隊服務系統(tǒng)）可采用離散型模擬。 排隊系統(tǒng)通常采用離散型模擬模型。排隊

7、系統(tǒng)通常采用離散型模擬模型。tQQ0ENDt只有新顧客到來或有顧客完成服務后離去，函數(shù)值才發(fā)生變化任意時刻的排隊的顧客數(shù)其中，發(fā)生系統(tǒng)狀態(tài)變化的事件有兩個：其中，發(fā)生系統(tǒng)狀態(tài)變化的事件有兩個：一是有顧客到達；二是服務員完成服務。一是有顧客到達；二是服務員完成服務。將最近發(fā)生上述兩種事件之一的時刻設置為下次事件發(fā)將最近發(fā)生上述兩種事件之一的時刻設置為下次事件發(fā)生點，就可將服務過程描述為下圖所示的模擬模型。生點，就可將服務過程描述為下圖所示的模擬模型。這個問題中有兩類狀態(tài)變量：在等待的顧客的人數(shù)（離散的非負整數(shù)）；服務員Ai是否正在工作（是或否）；對任意一個模擬，首先要作的是: 找出能完全描述任意

8、時刻的系統(tǒng)的狀態(tài)變量的集合： 從而得到平均隊伍的長度，最長的隊伍，顧客等待的平均時間以及兩個理發(fā)員的忙閑程度等，給出能從時刻 t 的狀態(tài)變量算出時刻 t+1 的新的狀態(tài)變量的程序。模擬的分類 第二類：連續(xù)型模擬第二類：連續(xù)型模擬。因變量隨時間的改變呈連續(xù)性。因變量隨時間的改變呈連續(xù)性變化。在大多數(shù)計算機模擬過程中，按固定的步長推變化。在大多數(shù)計算機模擬過程中，按固定的步長推進模擬時間。進模擬時間。 通常需建立一系列的由系統(tǒng)狀態(tài)變量組成的狀態(tài)方程通常需建立一系列的由系統(tǒng)狀態(tài)變量組成的狀態(tài)方程組，以描述狀態(tài)變量與模擬時間的關系。組，以描述狀態(tài)變量與模擬時間的關系。 第三類：混合型模擬第三類：混合型

9、模擬。因變量隨時間的推移而作連續(xù)。因變量隨時間的推移而作連續(xù)性的變化并具有離散性的突變，如庫存控制系統(tǒng)。性的變化并具有離散性的突變，如庫存控制系統(tǒng)。模擬的方式 終態(tài)模擬終態(tài)模擬：在規(guī)定的時間：在規(guī)定的時間T內(nèi)進行模擬運行，內(nèi)進行模擬運行，時間達到時間達到T時，模擬終止。其性能指標明顯取時，模擬終止。其性能指標明顯取決于系統(tǒng)的初始狀態(tài)。決于系統(tǒng)的初始狀態(tài)。 穩(wěn)態(tài)模擬穩(wěn)態(tài)模擬：隨著模擬時間的推移，系統(tǒng)的性能：隨著模擬時間的推移，系統(tǒng)的性能逐漸趨于平穩(wěn)。其目的是研究非終態(tài)系統(tǒng)長期逐漸趨于平穩(wěn)。其目的是研究非終態(tài)系統(tǒng)長期運行條件下的穩(wěn)態(tài)性能，模擬時間的長短取決運行條件下的穩(wěn)態(tài)性能，模擬時間的長短取決于

10、能否獲得系統(tǒng)性能的優(yōu)良估計（可由模擬輸于能否獲得系統(tǒng)性能的優(yōu)良估計（可由模擬輸出的精度確定）。出的精度確定）。模擬的一般步驟模擬的一般步驟 明確問題，建立模型。明確問題，建立模型。明確模擬目的，確定模擬輸出結果的目標函數(shù)分析各狀態(tài)變量之間關系，建立系統(tǒng)模型 收集和整理數(shù)據(jù)資料。特別是隨機性資料。收集和整理數(shù)據(jù)資料。特別是隨機性資料。分析收集的隨機數(shù)據(jù)，確定系統(tǒng)中隨機性因素的概率分布特征，以此為依據(jù)產(chǎn)生抽樣數(shù)據(jù) 編制程序，模擬運行。編制程序，模擬運行。編程、設定初始狀態(tài)，模擬運行時間、隨機樣本量、模擬運行次數(shù) 分析模擬輸出結果：分析模擬輸出結果：模擬結果的統(tǒng)計特性（模擬結果的統(tǒng)計特性（樣本均值、

11、方差、置信區(qū)間等）靈敏性分析（靈敏性分析（輸入?yún)?shù)做微小變化是否導致輸出結果巨大改變）選擇最優(yōu)方案。選擇最優(yōu)方案。系統(tǒng)用其他方法是是否用仿真否構造模型數(shù)據(jù)準備編程程序運行分析仿真結果結果是否充分是停止否模型是否有效修改模型策略或重復運行是否修改模型系統(tǒng)模擬注意事項：系統(tǒng)模擬注意事項：l一次模擬結果毫無意義！一次模擬結果毫無意義！l模擬是試驗性的模擬是試驗性的, ,是思維結果的驗證。是思維結果的驗證。l必須進行足夠多次的模擬，并對結果進行統(tǒng)計分必須進行足夠多次的模擬，并對結果進行統(tǒng)計分析。析。 隨機數(shù)的產(chǎn)生隨機數(shù)的產(chǎn)生 數(shù)學方法常見產(chǎn)生均勻分布隨機數(shù)的幾種方法的有平數(shù)學方法常見產(chǎn)生均勻分布隨機數(shù)

12、的幾種方法的有平方取中法、倍積取中法、乘同余法、二階與三階線性方取中法、倍積取中法、乘同余法、二階與三階線性同余法同余法 由均勻分布產(chǎn)生各種分布的隨機數(shù)、反函數(shù)法、取舍由均勻分布產(chǎn)生各種分布的隨機數(shù)、反函數(shù)法、取舍法、法、Box-Muller方法和極方法方法和極方法。 。 隨機數(shù)的產(chǎn)生隨機數(shù)的產(chǎn)生 :二項分布隨機數(shù)據(jù)二項分布隨機數(shù)據(jù) 正態(tài)分布的隨機數(shù)據(jù)的產(chǎn)生正態(tài)分布的隨機數(shù)據(jù)的產(chǎn)生 randperm(n)：產(chǎn)生一個1到n的隨機順序。randint(m,n,1 N)：生成mn的在1到N之間的隨機整數(shù)矩陣，rand(m,n)：生成0到1之間的mn的隨機數(shù)矩陣 randn是均值為0方差為1的正態(tài)分布

13、 常見分布的隨機數(shù)產(chǎn)生常見分布的隨機數(shù)產(chǎn)生 例例在我方某前沿防守地域，敵人以一個炮排（含兩在我方某前沿防守地域，敵人以一個炮排（含兩門火炮）為單位對我方進行干擾和破壞為躲避我方門火炮）為單位對我方進行干擾和破壞為躲避我方打擊，敵方對其陣地進行了偽裝并經(jīng)常變換射擊地打擊，敵方對其陣地進行了偽裝并經(jīng)常變換射擊地點點 經(jīng)過長期觀察發(fā)現(xiàn)，我方指揮所對敵方目標的指經(jīng)過長期觀察發(fā)現(xiàn)，我方指揮所對敵方目標的指示有示有5050是準確的，而我方火力單位，在指示正確是準確的，而我方火力單位，在指示正確時，有時，有1/31/3的射擊效果能毀傷敵人一門火炮，有的射擊效果能毀傷敵人一門火炮，有1/61/6的射擊效果能全

14、部消滅敵人的射擊效果能全部消滅敵人 現(xiàn)在希望能用某種方式把我方將要對敵人實施現(xiàn)在希望能用某種方式把我方將要對敵人實施的的2020次打擊結果顯現(xiàn)出來，確定有效射擊的比率及次打擊結果顯現(xiàn)出來，確定有效射擊的比率及毀傷敵方火炮的平均值。毀傷敵方火炮的平均值。分析：分析：這是一個概率問題，可以通過理論計算得到這是一個概率問題，可以通過理論計算得到相應的概率和期望值相應的概率和期望值. .但這樣只能給出作戰(zhàn)行動的但這樣只能給出作戰(zhàn)行動的最終靜態(tài)結果，而顯示不出作戰(zhàn)行動的動態(tài)過程最終靜態(tài)結果，而顯示不出作戰(zhàn)行動的動態(tài)過程. . 為了能顯示我方為了能顯示我方2020次射擊的過程，現(xiàn)采用模擬的方式。次射擊的過

15、程，現(xiàn)采用模擬的方式。 需要模擬出以下兩件事：需要模擬出以下兩件事： 1. 1. 問題分析問題分析2 2 當指示正確時，我方火力單位的射擊結果情況當指示正確時，我方火力單位的射擊結果情況1 1 觀察所對目標的指示正確與否觀察所對目標的指示正確與否模擬試驗有兩種結果，每一種結果出現(xiàn)的概率都是模擬試驗有兩種結果，每一種結果出現(xiàn)的概率都是1/21/2 因此，因此，可用投擲一枚硬幣的方式予以確定可用投擲一枚硬幣的方式予以確定，當硬幣出現(xiàn)正面時為，當硬幣出現(xiàn)正面時為指示正確，反之為不正確指示正確，反之為不正確 模擬試驗有三種結果：毀傷一門火炮的可能性為模擬試驗有三種結果：毀傷一門火炮的可能性為1/3(1

16、/3(即即2/6)2/6)，毀傷兩門的可能性為毀傷兩門的可能性為1/61/6，沒能毀傷敵火炮的可能性為，沒能毀傷敵火炮的可能性為1/2(1/2(即即3/6)3/6) 這時這時可用投擲骰子的方法來確定可用投擲骰子的方法來確定：如果出現(xiàn)的是、三個點：則認為沒能擊中敵人；如果出現(xiàn)的是、三個點：則認為沒能擊中敵人；如果出現(xiàn)的是、點：則認為毀傷敵人一門火炮；如果出現(xiàn)的是、點：則認為毀傷敵人一門火炮；若出現(xiàn)的是點：則認為毀傷敵人兩門火炮若出現(xiàn)的是點：則認為毀傷敵人兩門火炮2. 2. 符號假設符號假設i i：要模擬的打擊次數(shù)；要模擬的打擊次數(shù)； k k1 1：沒擊中敵人火炮的射擊總數(shù)；沒擊中敵人火炮的射擊總

17、數(shù)； k k2 2：擊中敵人一門火炮的射擊總數(shù)；擊中敵人一門火炮的射擊總數(shù)；k k3 3：擊中敵人兩門火炮的射擊總數(shù)擊中敵人兩門火炮的射擊總數(shù)E E：有效射擊比率；有效射擊比率； E E1 1：2020次射擊平均每次毀傷敵人的火炮數(shù)次射擊平均每次毀傷敵人的火炮數(shù)3. 3. 模擬框圖模擬框圖初始化:i=0,k1=0,k2=0,k3=0i=i+1骰子點數(shù)?k1=k1+1k2=k2+1k3=k3+1k1=k1+1i20?E=(k2+k3)/20 E1=0*k1/20+1*k2/20+2*k3/20停止硬幣正面?YNNY1，2，34，56 理論計算和模擬結果的比較 分類 項目 無效射擊 有效射擊 平均

18、值 模 擬 74 26 34 理 論 R2=unidrnd(6) 企業(yè)生產(chǎn)的庫存系統(tǒng)企業(yè)生產(chǎn)的庫存系統(tǒng) 在銷售部門、工廠等領域中都存在庫存問題，庫存太多造成在銷售部門、工廠等領域中都存在庫存問題，庫存太多造成浪費以及資金積壓，庫存少了不能滿足需求也造成損失，部門的浪費以及資金積壓，庫存少了不能滿足需求也造成損失，部門的工作人員需要決定何時進貨，進多少，使得所花費的平均費用最工作人員需要決定何時進貨，進多少，使得所花費的平均費用最少，而收益最大，這就是庫存問題少，而收益最大，這就是庫存問題 某企業(yè)生產(chǎn)易變質的產(chǎn)品當天生產(chǎn)的產(chǎn)品必須當天售出，某企業(yè)生產(chǎn)易變質的產(chǎn)品當天生產(chǎn)的產(chǎn)品必須當天售出，否則就

19、會變質該產(chǎn)品單位成本為否則就會變質該產(chǎn)品單位成本為25 元，單位產(chǎn)品售價為元，單位產(chǎn)品售價為5 元企業(yè)為避免存貨過多而造成損失，擬從以下兩種庫存方案中元企業(yè)為避免存貨過多而造成損失，擬從以下兩種庫存方案中選出一個較優(yōu)的方案選出一個較優(yōu)的方案 方案甲：按前一天的銷售量作為當天的貨存量；方案甲：按前一天的銷售量作為當天的貨存量； 方案乙：按前二天的平均銷售量作為當天的貨存量；方案乙：按前二天的平均銷售量作為當天的貨存量； 假定市場對該產(chǎn)品的每天需求量是一個隨機變量，但從以往假定市場對該產(chǎn)品的每天需求量是一個隨機變量，但從以往的統(tǒng)計分析得知它服從正態(tài)分布的統(tǒng)計分析得知它服從正態(tài)分布 2(135,22

20、.4 )N用蒙特卡洛法解非線性規(guī)劃問題用蒙特卡洛法解非線性規(guī)劃問題對于非非線線性性規(guī)規(guī)劃劃問問題題： min f(X) XnE s.t. 0)(Xgi i=1,2, m jjjbxa j=1,2, n基本假設基本假設 試驗點的第試驗點的第j j個分量個分量x xj j服從服從 a aj j ,b ,bj j 內(nèi)的均勻分布內(nèi)的均勻分布符號假設 P: 試驗點總數(shù)；MAXP：最大試驗點總數(shù)；K: 可行點總數(shù)； MAXK：最大可行點數(shù);X*：迭代產(chǎn)生的最優(yōu)點； Q：迭代產(chǎn)生的最小值 f(X*)，其初始值為計算機所能表示的最大數(shù)求解過程求解過程 先產(chǎn)生一個隨機數(shù)作為初始試驗點先產(chǎn)生一個隨機數(shù)作為初始試驗

21、點, ,以后則將上一個試驗以后則將上一個試驗點的第點的第j j個分量隨機產(chǎn)生，其它分量不變而產(chǎn)生一新的試驗個分量隨機產(chǎn)生，其它分量不變而產(chǎn)生一新的試驗點這樣，每產(chǎn)生一個新試驗點只需一個新的隨機數(shù)分點這樣，每產(chǎn)生一個新試驗點只需一個新的隨機數(shù)分量當量當KMAXKKMAXK或或PMAXPPMAXP時停止迭代時停止迭代框框 圖圖初始化:給定MAXK,MAXP;k=0,p=0,Q:大整數(shù)xj=aj+R(bj-aj) j=1,2,nj=0j=j+1,p=p+1PMAXP?YNxj=aj+R(bj-aj)gi(X)0?i=1,2nYNjMAXK?YN輸出X,Q,停止YN 例例 max 2121222138

22、2xxxxxxz s.t 10321 xx 01x 02x 在在MatlabMatlab軟件包中編程，共需三個文件軟件包中編程，共需三個文件: :randlp.mrandlp.m, , mylp.m,mylp.m, lpconst.m.lpconst.m.主程序為主程序為randlp.m.randlp.m.% mylp.mfunction z=mylp(x) %目標函數(shù)目標函數(shù)z=2*x(1)2+x(2)2-x(1)*x(2)-8*x(1)-3*x(2); %轉化為求最小值問題轉化為求最小值問題% lpconst.m% lpconst.m function lpc=lpconst(x) %約束

23、條件 if 3*x(1)+x(2) -10= -0.5 %約束條件的誤差為5 .0 lpc=1; else lpc=0; end % randlp.m function sol,r1,r2=randlp(a,b,n) %隨機模擬解非線性規(guī)劃隨機模擬解非線性規(guī)劃debug=1;a=0; %試驗點下界試驗點下界b=10; %試驗點上界試驗點上界n=1000; %試驗點個數(shù)試驗點個數(shù)r1=unifrnd(a,b,n,1);%n1階的階的a,b均勻分布隨機數(shù)矩陣均勻分布隨機數(shù)矩陣r2=unifrnd(a,b,n,1);sol=r1(1) r2(1);z0=inf;for i=1:n x1=r1(i);

24、 x2=r2(i); lpc=lpconst(x1 x2)； if lpc=1 z=mylp(x1 x2)； if zz0 z0=z； sol=x1 x2； end endend 在計算機上產(chǎn)生隨機數(shù)最實用、最常見的方法是數(shù)學方法，即用如下遞推公式： 產(chǎn)生隨機數(shù)序列。對于給定的初始值1,2，k，確定n+k，=1，2，。, 2 , 1),(11nTknnnkn 用數(shù)學方法產(chǎn)生的隨機數(shù)，存在兩個問題：遞推公式和初始值遞推公式和初始值1,2，k確定后，整個隨機數(shù)序確定后，整個隨機數(shù)序列便被唯一確定。列便被唯一確定。不滿足隨機數(shù)相互獨立的要求。由于隨機數(shù)序列是由遞推公式確定的，而在計算機上所能表示的0

25、，1上的數(shù)又是有限的，因此，這種方法產(chǎn)生的隨機數(shù)序列就不可能不出現(xiàn)無限重復。a) 隨機數(shù)序列便出現(xiàn)了周期性的循環(huán)現(xiàn)象出現(xiàn)了周期性的循環(huán)現(xiàn)象。 由于這兩個問題的存在，常稱用數(shù)學方法產(chǎn)生的隨機數(shù)為偽隨機數(shù)。對于以上存在的兩個問題，作如下具體分析。 關于第一個問題，不能從本質上加以改變，但只要遞推公式選得比較好，隨機數(shù)間的相互獨立性是可相互獨立性是可以近似滿足以近似滿足的。至于第二個問題，則不是本質的。因為用蒙特卡羅方法解任何具體問題時，所使用的隨機數(shù)的個數(shù)總是有限的，只要所用隨機數(shù)的個數(shù)不超過隨機數(shù)的個數(shù)不超過偽隨機數(shù)序列出現(xiàn)循環(huán)現(xiàn)象時的長度偽隨機數(shù)序列出現(xiàn)循環(huán)現(xiàn)象時的長度就可以了。 用數(shù)學方法產(chǎn)

26、生的偽隨機數(shù)容易在計算機上得到，可以進行復算，而且不受計算機型號的限制。因此，這種方法雖然存在著一些問題，但仍然被廣泛地在計算機上使用，是在計算機上產(chǎn)生偽隨機數(shù)的主要方法。 一道工序用自動化車床連續(xù)加工某種零件，由于刀具損壞等原因該工序會出現(xiàn)故障，其中刀具損壞故障占95%，其它故障僅占5%。工序出現(xiàn)故障是完全隨機的，假定在生產(chǎn)任一零件時出現(xiàn)故障的機會均相同。工作人員通過檢查零件來確定工序是否出現(xiàn)故障。 現(xiàn)積累有100次刀具故障記錄，故障出現(xiàn)時該刀具完成的零件數(shù)如附表?，F(xiàn)計劃在刀具加工一定件數(shù)后定期更換新刀具。 已知生產(chǎn)工序的費用參數(shù)如下： 故障時產(chǎn)出的零件損失費用 f = 200元/件； 進行

27、檢查的費用 x = 10元/次； 發(fā)現(xiàn)故障進行調(diào)節(jié)使恢復正常的平均費用 d = 3000元/次 未發(fā)現(xiàn)故障時更換一把新刀具的費用 k = 1000元/次。1）假定工序故障時產(chǎn)出的零件均為不合格品，正常時產(chǎn)出的零件均為合格品，試對該工序設計效益最好的檢查間隔和刀具更換策略。2）如果該工序正常時產(chǎn)出的零件不全是合格品，有2%為不合格品：而工序故障時產(chǎn)出的零件有40%為合格品，60%為不合格品。工序正常而誤認為有故障停機產(chǎn)生的損失費用為1500元/次。對該工序設計效益最好的檢查間隔和刀具更換策略。1999年年A題題 自動化車床管理自動化車床管理 較復雜的非線性優(yōu)化問題1、驗證模型的正確性，采用計算機

28、模擬以零件的序列為模擬序列確定隨機事件，如主要隨機事件-刀具的損害作為，可分析證明刀具壽命服從正態(tài)分布編制程序模擬事件發(fā)生過程對不同的給定檢查間隔和換刀間隔進行模擬，得到單位零件平均工作周期內(nèi)的資金損耗比較模型結果和模擬結果2、檢驗假設如第二問中，加入了隨機因子，問題更為復雜，建模過程中不免加入一些假設使問題更簡單化。如由于非刀具故障占總故障的百分比較小，為5%，故假設其對結果影響不大，從而不考慮非刀具故障。可通過計算機模擬對此假設進行驗證。3、對多個策略進行仿真比較如對問題二通過分析提出如下幾個檢查策略（h為檢查間隔）策略1、檢查第h個零件，如為正品，則繼續(xù)生產(chǎn)；如為廢品，停止生產(chǎn)，換刀。策略2、檢查第h個零件，如為正品，則繼續(xù)生產(chǎn)；如為廢品，則檢查第h+1個，若正品，繼續(xù)生產(chǎn)，如廢品，換刀。策略3、檢查第h個零件，如為廢品，則換刀；如為正品，則檢查第h+1個，若正品，繼續(xù)生產(chǎn)，如廢品，換刀。策略4、檢查第h和h+1個零件，如兩件都為正品，則繼續(xù)生產(chǎn)；如兩件都為廢品，換刀；如一正一費，則檢查第h+2個，若正品，繼續(xù)生產(chǎn)，如廢品，換刀。至多到T個產(chǎn)品一定換刀