西北工業(yè)大學數(shù)學建模課程論文（個人習作，文庫共享，僅供參考）Skyone企業(yè)事故管理能力形成機制探討摘要企業(yè)事故管理能力是指企業(yè)預防事故發(fā)生和事故發(fā)生后如何處理的能力。這種能力的形成需要企業(yè)外部壓力與激勵機制、內(nèi)部的動力轉(zhuǎn)換機制和時間來共同完成。問題一研究企業(yè)內(nèi)部機制與外部機制之間響應的數(shù)學模型。首先要將內(nèi)部機制與外部機制進行量化，本文用事故經(jīng)費的投入數(shù)額來表示企業(yè)內(nèi)部轉(zhuǎn)換機制的大小，通過層次分析法將外部機制的5個因素進行量化。然后通過函數(shù)擬合分別找出內(nèi)部機制與時間以及外部機制與時間的關系。最終以時間為中間變量，建立了企業(yè)內(nèi)部機制與外部機制變化的函數(shù)關系：，并用數(shù)據(jù)進行了檢驗。問題二要建立外部機制、時間與事故管理能力的關系。首先對事故管理能力進行量化，用各年死亡率變化的絕對值之和來表示該年的事故管理能力。然后通過微分方程建立事故管理能力與時間的關系，對數(shù)據(jù)進行擬合得到k=0.3899，m=0.9552。并利用問題一中外部因素與時間的關系：，通過函數(shù)復合計算得到事故管理能力與外部機制的關系，最終利用控制變量法，控制外部五個因素依次變化，得到影響事故管理能力的主要因素：強制性因素，行政壓力與激勵，經(jīng)濟激勵與壓力。問題三需要在問題二的基礎上考慮內(nèi)部動力轉(zhuǎn)換機制對企業(yè)事故管理能力的影響，問題一和二已經(jīng)分別計算出了企業(yè)內(nèi)部轉(zhuǎn)化機制、外部機制因素向量以及企業(yè)事故管理能力與時間的關系，由此得到了四個因素之間的關系函數(shù)模型，見公式（15）。最后發(fā)現(xiàn)在最初的30多年間，內(nèi)部轉(zhuǎn)換機制的影響起主要作用，到了后期事故管理能力受制于綜合因素，趨于穩(wěn)定值8.5137。問題四針對如何提高企業(yè)事故管理能力，降低事故死亡率給出了符合實際的幾條建議。關鍵詞：內(nèi)部轉(zhuǎn)換機制企業(yè)事故管理能力函數(shù)復合阻滯增長微分方程一、問題重述企業(yè)事故管理能力是指企業(yè)預防事故發(fā)生和事故發(fā)生后如何處理的能力。這種能力的形成需要幾個方面的機制來加以保證：一是企業(yè)外部的壓力機制和激勵機制，壓力機制主要包括法律法規(guī)的1強制性壓力、行政性壓力、經(jīng)濟性處罰、社會輿論壓力、市場負面反應的壓力等，激勵機制主要包括經(jīng)濟性激勵、社會輿論的褒獎、市場的積極反應、行政性激勵、社會支持等。這種壓力機制的作用機理是：外部激勵與壓力因素作用于企業(yè)管理團隊尤其是高管團隊，使他們產(chǎn)生一定的意識，從而采取提高事故管理能力的行動，逐步提高企業(yè)事故管理的能力。這種能力的形成表現(xiàn)就是在一定時期中單位生產(chǎn)量的事故數(shù)及傷亡人數(shù)的變動。二是企業(yè)內(nèi)部的動力轉(zhuǎn)換機制，即企業(yè)內(nèi)部如何將外部的壓力轉(zhuǎn)換為企業(yè)提高事故管理能力的積極性問題，這個動力機制的主要對象就是企業(yè)的高管團隊。只有高管團隊受到個人以外的刺激產(chǎn)生一定的意識，并認識到事故管理能力的重要性以及是否發(fā)生事故對管理團隊和管理者個人的獎勵和懲處之后，他們才有動力去驅(qū)動企業(yè)事故管理能力的形成與提高。動力轉(zhuǎn)換機制應該主要包括對高管團隊的經(jīng)濟激勵與壓力、行政激勵與壓力、經(jīng)理人市場激勵與壓力（就是說在經(jīng)理人市場上，雇主會以其領導的企業(yè)是否發(fā)生過事故來判斷經(jīng)理人的管理能力強弱，從而決定對高管人員的聘用與否或者薪酬高低）、法律壓力和社會壓力等。判別轉(zhuǎn)換機制是否啟動并發(fā)揮作用的主要指標是企業(yè)事故管理的經(jīng)費投入數(shù)量，以及是否采用了公認、規(guī)范的事故管理體系。通過轉(zhuǎn)換機制使得企業(yè)內(nèi)部部門和人員提高自身事故管理的意識與技能，實現(xiàn)事故預防而不發(fā)生事故、或者即使事故發(fā)生了也能夠?qū)⑹鹿屎蟮奶幚淼玫缴鐣蛦T工等多方面的認可，避免因事故而引發(fā)更大損失的目標。除上述機制外，企業(yè)事故管理能力的形成需要考慮的另外一個重要因素就是時間，能力形成機制的啟動與管理能力的形成與提高在時間上具有連續(xù)性與滯后性。請針對以下問題進行研究：1.請建立一個描述動力轉(zhuǎn)換機制對外部壓力與激勵機制響應的數(shù)學模型，并通過數(shù)值試驗檢驗模型的有效性。2.請在外部機制、時間與事故管理能力之間建立適當?shù)臄?shù)學模型，考察外部機制、時間對事故管理能力的影響，找出影響事故管理能力形成的主要因素。3.請綜合考慮企業(yè)外部機制、轉(zhuǎn)換機制、時間等影響企業(yè)事故管理能力的因素，建立適當?shù)臄?shù)學模型，找出影響事故管理能力形成的主要因素。4.要使企業(yè)事故管理能力持續(xù)增強、減少事故發(fā)生的數(shù)量、避免人員傷亡，給出可行性的政策建議。二、問題分析企業(yè)事故管理能力是指企業(yè)預防事故發(fā)生以及事故后如何處理的能力。這種管理能力的形成主要受到兩個機制的影響，一是企業(yè)外部的激勵機制與動力轉(zhuǎn)換機制；二是企業(yè)內(nèi)部的動力轉(zhuǎn)換機制。本題主要探討外部機制與內(nèi)部機制響應的數(shù)學模型，外部機制隨時間對企業(yè)事故管理能力的影響，以及內(nèi)外機制共同作用機制對事故管理能力形成的影響。問題一：由于企業(yè)的壓力機制和激勵機制是成對出現(xiàn)的，因此將題目中的10個外部因素歸為5個。企業(yè)的內(nèi)部動力轉(zhuǎn)換機制（以下簡稱內(nèi)部機制）的作用用各年份投入的經(jīng)費數(shù)額和單位產(chǎn)量死亡率來描述。首先建立與內(nèi)部機制與時間的關系，再建立5個因素與時間的關系，然后利用復合函數(shù)的方法，以時間作為中間變量，復合得到了內(nèi)部機制與5個外部因素之間的關系，即是內(nèi)部的動力轉(zhuǎn)換機制對外部的壓力與激勵機制的響應模型。問題二：問題一中已經(jīng)建立了時間與外部機制的關系，因此只要建立時間與事故管理能力的關系，就可以得到外部機制、時間共同對事故管理能力的作用。關于事故管理能力的描述，考慮到能力形成機制的啟動與管理能力的形成與提高在時間上具有的連續(xù)性與滯后性，本文采用了相鄰兩年的死亡率的變動來描述事故管理能力的形成，對進行進一步分析后，建立了微分模型公式，然后基于微分模型公式對時間與事故管理能力的數(shù)值進行擬合。問題三：基于問題一二得到的結(jié)果，建立了時間、企業(yè)內(nèi)部轉(zhuǎn)換機制、企業(yè)外部機制（包括五個因素）以及事故管理能力四個變量兩兩之間的任意關系（見公式15），最后通過公式（15）分析了事故管理能力與企業(yè)內(nèi)部轉(zhuǎn)換機制的關系。三、模型假設與符號說明3.1模型假設1、假設企業(yè)外部機制的作用隨時間不斷增強；2、假設同一因素的激勵和壓力同時存在；3、假設死亡人數(shù)或傷亡數(shù)不變化，則事故管理能力就不會變化。4、假設事故管理能力隨時間變化為單調(diào)不減；5、內(nèi)部動力轉(zhuǎn)換機制的大小等于各年份投入的經(jīng)費數(shù)額6、假設外部的壓力機制和激勵機制都會促進企業(yè)高管產(chǎn)生意識，提高事故管理能力；3.2符號說明符號意義單位11企業(yè)投入的事故管理費千萬(元)時間年產(chǎn)量千萬單位千萬單位產(chǎn)量死亡率1相鄰兩年死亡率的變動1四、模型分析與求解4.1數(shù)據(jù)預處理對附件中給的數(shù)據(jù)進行預處理，根據(jù)，即，求得企業(yè)從1949年到2008年各年的產(chǎn)量數(shù)據(jù)，單位千萬。根據(jù)，即，得到企業(yè)歷年投入的事故管理經(jīng)費的數(shù)據(jù)。通過上面數(shù)據(jù)，做出的散點圖，如圖1。圖1：事故管理經(jīng)費隨時間變化散點圖4.2問題一4.2.2模型建立根據(jù)上面對問題一的分析，以及觀察圖4-1發(fā)現(xiàn)，事故管理經(jīng)費隨時間逐年增長，擬用指數(shù)模型和阻滯模型進行擬合。指數(shù)模型建立函數(shù)關系(1)為了找到事故管理經(jīng)費與5個外部激勵或壓力因素變化的關系，基于假設1（假設企業(yè)外部機制的作用隨時間不斷增強），在此為了便于計算，進一步假設外部機制與時間成正比，建立了時間與5個外部因素的函數(shù)關系：（2）通過MATLAB進行數(shù)據(jù)擬合后，將公式(1)(2)進行復合，可以得到事故管理經(jīng)費與外部因素之間的關系（3）即是內(nèi)部動力轉(zhuǎn)換機制對外部壓力因素與激勵因素的響應模型。阻滯型擬合考慮阻滯型擬合，建立經(jīng)費與時間之間的函數(shù)關系：（4）運用與指數(shù)型模型相同的方法，建立時間與5個外部因素的量化關系。將公式(2)(4)進行復合，可以得到阻滯模型下事故管理經(jīng)費與外部因素之間的關系（5）即是阻滯模型下內(nèi)部動力轉(zhuǎn)換機制對外部壓力因素與激勵因素的響應模型。4.2.3模型求解指數(shù)模型根據(jù)表格中的數(shù)據(jù)，對與時間按照進行擬合，得到函數(shù)擬合的圖形如圖2所示。均方根誤差為：，并且分別得到兩個參數(shù)的值如下：將兩個參數(shù)代入公式（1）得到內(nèi)部轉(zhuǎn)換機制與時間的指數(shù)增長函數(shù)如下：（6）圖2：內(nèi)部轉(zhuǎn)換機制與時間的指數(shù)擬合阻滯型擬合對內(nèi)部動力轉(zhuǎn)換機制與時間按照進行擬合，函數(shù)擬合的圖形如圖3所示，擬合的均方根誤差為：。分別得到三個參數(shù)的值如下：將三個參數(shù)代入公式（4），得到內(nèi)部動力轉(zhuǎn)換機制與時間的阻滯型函數(shù)關系:(7)圖3：內(nèi)部轉(zhuǎn)換機制與時間的阻滯型擬合通過對比兩種模型的擬合參數(shù)均方根誤差（指數(shù)擬合與阻滯型擬合），發(fā)現(xiàn)阻滯型模型擬合效果更好，這也與實際的情況相符。因此采用阻滯型擬合得到的公式（7）作為最終結(jié)果。.時間與外部因素的關系模型根據(jù)公式（2），，外部的5個激勵與壓力因素，分別對應社會，行政，經(jīng)濟，輿論與市場反應。由于這5個因素并無可用數(shù)據(jù)，且即使在實際中也很難量化，因此本文采用如下處理方法。。由于各個因素通過時間作用于動力轉(zhuǎn)換機制，然而實際中反映這些因素的數(shù)據(jù)幾乎沒有，在建立成對比較矩陣時，我們通過對題中給出的各個因素（包括法律法規(guī)、經(jīng)濟處罰與社會輿論等）進行了詳細的分析，進一步選用了廣泛使用且精確度較高的1—9尺度來表示他們之間的差別。下面介紹用層次分析法來確定系數(shù)的取值。構(gòu)造成對比較因素的矩陣用MALAB計算得到特征向量、特征根，從而得到五個變量的權(quán)重如下表：表1：外部因素的各個權(quán)重因素強制性行政性經(jīng)濟性社會輿論市場性權(quán)值根據(jù)公式，最終得到內(nèi)部動力轉(zhuǎn)換機制對外部壓力與激勵機制的響應模型為:即：（8）4.2.3數(shù)值檢驗由于內(nèi)部轉(zhuǎn)換機制與外部因素的關系是一個多元函數(shù)，不同于一元函數(shù)的關系可以通過求反函數(shù)的方法直接找到結(jié)果，而且在此外部因素的大小值不唯一。所以為了找到各個年份一個因素向量的大小，通過對因素向量隨機選取區(qū)間之間的數(shù)值，代入公式：,得到內(nèi)部轉(zhuǎn)換機制的計算值，通過與實際值比較，得到了誤差最小的一組因素向量的數(shù)值大?。ú晃ㄒ唬?，給出了一組數(shù)值見附錄一，圖4表示與的函數(shù)關系。圖4：內(nèi)部對外部響應的數(shù)值檢驗結(jié)果4.3問題二4.3.1理論分析企業(yè)高管事故管理能力隨動力轉(zhuǎn)換機制的啟動而產(chǎn)生，隨工作動力與積極性在時間上的持續(xù)而提高。不過企業(yè)事故管理能力不可能無限量增加，會有一定的定上限值，一般認為能力的變化率隨時間的增加會逐漸減小，最終趨近0。通過上述分析，指數(shù)函數(shù)符合這一變化規(guī)律，故可建立微分模型如下： 求解上述方程，可建立企業(yè)事故管理能力對動力機制在時間上的積分模型： （8）而且在由此分析，企業(yè)事故管理能力對動力機制的積分模型有其合理性。下面通過函數(shù)擬合事故管理能力于時間的關系來確定參數(shù)k和m。4.3.2模型建立事故管理能力模型根據(jù)以上對問題二的分析，用死亡率的變動來描述事故管理能力的形成?？紤]到它的滯后性以及連續(xù)性，以及假設（3），第年的事故管理能力為上一年的事故管理能力加上該年事故管理能力的變化，即或表示為(9)其中，事故管理能力與時間模型由于公式（8）的函數(shù)很難用軟件進行直接的擬合，所以為了完整擬合模型的函數(shù)，進一步建立如下函數(shù)關系，可以證明，公式（8）可以化簡為公式（10），所以理論上公式（8）與公式（10）是等價的。(10)在問題一中已經(jīng)建立了時間與外部因素之間的關系，見公式（2），因此通過復合可以得到事故管理能力與外部因素的關系。(11)主要影響因素模型為了確定主要因素，采用控制變量法，使得一個因素不動，另一個因素變化，觀察隨時間變化的趨勢，取變化速率快的作為主要影響因素。4.3.3模型求解求解函數(shù)模型根據(jù)上述模型得到的公式（10），利用MATLAB對事故管理能力與時間這兩個變量進行數(shù)據(jù)擬合，得到擬合關系圖4-4，擬合均方根誤差為：，得到三個參數(shù)如下：代入公式（10），得到事故管理能力隨時間變化的函數(shù)關系，（12）.圖5：事故管理能力隨時間變化的關系擬合主要因素將所求的結(jié)果代入公式（11），得到外部因素對事故管理能力的影響函數(shù)關系如下：（13）為了尋找影響事故管理費投入的主要因素，對于，，，，這五個影響因素，采用控制變量法，根據(jù)上述公式，分別令一個因素由0到1變化，保持其他因素不變（設之為0），得到事故管理能力關于各個因素的控制變量之后的函數(shù)關系公式（14），如圖6所示。（14）圖6：外部5個因素與事故管理能力的關系分析圖6，容易發(fā)現(xiàn)事故管理能力隨外界社會支持因素變化的速度最快，是最主要因素。另外行政因素和經(jīng)濟因素也是主要因素。因此可認為主要因素為社會支持（或強制因素）、行政因素和經(jīng)濟因素。4.4問題三綜合問題一與問題二的結(jié)果，我們得到了企業(yè)內(nèi)部轉(zhuǎn)化機制與時間的函數(shù)關系，企業(yè)企業(yè)外部機制與時間的函數(shù)關系，以及時間與企業(yè)事故管理能力之間的關系。以時間作為中間變量，通過函數(shù)復合可以得到時間、企業(yè)內(nèi)部轉(zhuǎn)換機制、企業(yè)外部機制（包括五個因素）以及事故管理能力四個變量兩兩之間的任意關系。四者之間的關系如下表達式：（15）對于外部機制、時間與事故管理能力之間的關系，第二問已經(jīng)做了討論，下面重點研究討論事故管理能力與內(nèi)部動力轉(zhuǎn)換機制之間的關系。即討論下面的隱函數(shù)。由于很難直接求解出與的關系，所以本文以時間為中間變量，將函數(shù)表達式的數(shù)值解離散量化后，對數(shù)值做圖得到如下圖7（16）圖7：事故管理能力與內(nèi)部轉(zhuǎn)換機制的關系圖通過對圖7的分析，發(fā)現(xiàn)當內(nèi)部轉(zhuǎn)換機制大概處于在0到之間時（此時的時間正好為1949年到1980年前后），增長速度非?？?，即增加單位的內(nèi)部轉(zhuǎn)換機制的大小，事故管理能力增長效果明顯，而到了后期即1980年之后，事故管理能力趨于穩(wěn)定值。由公式（15）可知，事故管理能力受到了各種綜合因素的影響，包括時間、企業(yè)內(nèi)部轉(zhuǎn)換機制、企業(yè)外部機制（包括五個因素），所以不能僅僅的找出某一因素的決定性影響大小。經(jīng)過對公式（10）的計算得到公式（8）的兩個參數(shù)值：，進一步計算，即最終事故管理能力趨于8.5137。4.5問題四通過以上分析，我們發(fā)現(xiàn)，要想提高企業(yè)的事故管理能力，首先要加大法律法規(guī)的強制性壓力，即加重對嚴重事故領導人的刑事處罰，或者更多給與事故管理知識體系支持，更多的給與社會支持與幫助；其次是行政性激勵，更多的激勵企業(yè)高管，進行企業(yè)領導行政職務的提升；最后可以進行經(jīng)濟性的激勵與壓力。1.加強法律法規(guī)的建設，進一步細化法律條文，增加法律的指導能力，為企業(yè)的生產(chǎn)與政府的管理奠定明確的法律基礎，可以逐步出臺針對企業(yè)高管人員對于企業(yè)事故的相關約束性法律，促使企業(yè)從決策的角度就開始重視企業(yè)事故管理能力的建設，為了從源頭上實現(xiàn)企業(yè)事故的逐步下降奠定基礎。2.應該注意到實踐性影響因素的重要性，企業(yè)應該重視時間對于企業(yè)是管理能力投入的連續(xù)性與滯后性，現(xiàn)在足額的事故管理能力投入是確保未來生產(chǎn)活動正常進行和企業(yè)正常運轉(zhuǎn)的必須條件，不能為了一時的企業(yè)利益的增長而忽視企業(yè)事故管理能力建設上應有的投入而為企業(yè)的發(fā)展埋下隱患。3.國家應該在輿論導向上增加人民大眾對于企業(yè)事故的關注度，從而形成正確而廣泛的群眾性基礎，為全社會形成安全生產(chǎn)的社會氣氛創(chuàng)造條件，為企業(yè)重視和加強企業(yè)事故管理能力創(chuàng)造積極地社會輿論。五、模型優(yōu)缺點優(yōu)點：1.將難以量化的企業(yè)內(nèi)部轉(zhuǎn)換機制、企業(yè)外部機制因素（包括五個因素）以及事故管理能力，通過對實際問題的合理分析，簡化進行合理的數(shù)值量化，使得論文變得簡潔明了，具有較強的說服力。2.在對企業(yè)外部機制因素進行量化時，考慮到五個因素強制性激勵與壓力，行政性激勵與壓力，經(jīng)濟性激勵與壓力，以及輿論與市場的激勵與壓力。這五個因素只給出了相關的法律，難以量化比較，故采取了層次分析法量化。3.在所有因素之中，時間是一個最基本的因素，將所有變量均和時間產(chǎn)生聯(lián)系，對問題的整個描述簡化到了最終的公式（15），這一公式很好的描述了時間、企業(yè)內(nèi)部轉(zhuǎn)換機制、企業(yè)外部機制因素（包括五個因素）以及事故管理能力四者之間的關系。缺點1.在對問題進行數(shù)值量化的過程中，雖然簡化了問題，但是同樣也忽略了一些其他的因素。2.在使用層次分析法的過程中，主觀性比較大。六、參考文獻[1]肖華勇，《實用數(shù)學建模與軟件應用》，西安：西北工業(yè)大學出版社，2008[2]孫蓬，《MATLAB基礎教程》，北京：清華大學出版社，2011.七、附錄附錄一一組因素向量的大?。ú晃ㄒ唬┮蛩貜娭菩孕姓越?jīng)濟性社會輿論市場性19490000019500.00210.00360.00530.00920.009219510.00350.00340.00670.01810.065319520.04430.07330.08340.00880.046119530.08440.09050.09760.06150.088419540.14540.10550.11170.18080.152019550.34620.34470.30850.34160.383519560.42890.41390.41830.47570.498919570.43220.47350.49190.45290.441219580.36950.32320.31340.33710.348819590.29770.28780.26710.20620.218219600.28250.27830.28710.24630.256119610.27610.27060.29730.25700.239719620.29870.29270.29320.29760.281019630.30960.33500.30420.38790.375419640.29740.22140.20890.26310.214419650.29620.25850.28080.25470.255019660.32900.33150.35600.34550.385719670.40840.40810.44010.46870.440919680.45670.44800.40110.48000.448719690.48580.42000.43470.41860.462319700.49690.43400.41260.42030.436219710.47550.47270.42450.42920.424019720.49880.49690.49600.48130.456019730.50760.50760.50830.50050.529819740.53420.54540.56910.58160.515519750.59000.58580.56160.50700.507119760.58560.58730.59560.55180.565419770.59720.59520.50700.56160.523019780.59780.56990.54460.51500.578519790.60210.61040.60700.60200.603619800.61970.61430.60150.60960.692619810.61190.69810.69360.64180.698519820.69270.62180.68280.69600.613419830.70000.67700.62380.69510.656919840.69810.69840.69140.67280.658819850.70010.71360.72080.72590.705519860.72100.73170.74600.74990.719119870.72990.74770.70460.79370.785019880.72460.72450.77990.76830.745319890.74350.74810.73330.77420.738519900.78420.71280.70140.74320.745719910.79970.74120.72180.76460.749119920.79930.79520.72110.79910.794619931.00001.00001.00001.00001.000019940.79380.78240.79650.77800.790719950.78190.74650.77850.73250.752819960.70560.73950.73970.77950.786119970.70270.71170.72320.70390.771719980.73320.74060.78680.76390.748419990.80620.80890.80980.83370.808920010.86590.89830.86010.83640.884920030.90050.92790.90470.94070.954220040.91180.90530.98720.97100.914420050.93580.97780.98850.96430.930220060.99950.99240.91300.95590.942320070.98610.99580.95620.98950.9792200811111附錄二程序附錄問題一阻滯型擬合源程序x=[12161.25160951990524933.752613031372.536862.54138549023.75101250138296.25148953.75104107.582331.2581401.2580463.758692594301.2577137.582346.2599731.25132746.25147112.5153926.25156363.75154938.75180840181293.75206505231697.5238327.5232548.75233111.25249870267948.75295961.25327105335265348033.75367451.25395306.25404737.5406605417956.25431763.75461073.75484567.5514239.8287496544.1176458217.1315390424.5283371032.4232419352.1595524625650336.9272733568.6465791519.7461871999.0201956060.6061967005.0761];n=60;xx=x(1:n);t=[ones(n,1),(1:n)'];y=log(xx(1:n));[b,bint,r,rint,stats]=regress(y,t);RR=stats(1);F=stats(2);prob=stats(3);x0=exp(b(1));r=b(2);py=x0*exp(r*t(:,2));err=xx-py;rmse=sqrt(sum(err.^2)/n)plot(1:n,xx,'*',1:n,py);問題一的數(shù)值實驗程序clcsymsty=3.0506*10^4*exp(0.0586*t);n=1000;X=xlsread('經(jīng)費.xls');x=zeros([1,5]);v=10^6*ones([60,1]);Y=ones([60,5]);fori=1:nforj=0:0.1:0.9fori=1:5x(i)=0.1*rand+j;endtime=60*(0.5540*x(1)+0.1903*x(2)+0.1563*x(3)+0.0603*x(4)+0.0391*x(5));m=subs(y,time);[k,l]=min(abs(X-m));if(abs(v(l)-X(l))>abs(m-X(l)))v(l)=m;fork=1:5Y(l,k)=x(k);endendendendYplot(1949:2008,v,'+');holdonplot(1949:2008,X)問題二的擬合程序y=[00.344671.036031.7283272.0545622.0603442.3091142.4195642.5053143.3012763.7259043.8325963.8576464.1431364.4804964