安徽工業(yè)大學管理科學與工程學院 系統(tǒng)工程實驗指導書系統(tǒng)工程實驗報告 姓 名： 班 級： 學 號： 指導老師： 2013年 12 月 6 日 實驗一 簡單庫存模型一、 實驗目的1、 熟悉STELLA軟件的基本操作2、 加深對系統(tǒng)動力學主要要素和基本思想的理解3、 學會利用STELLA軟件建立一階反饋系統(tǒng)模型、仿真運行及結(jié)果分析二、 實驗要求1、簡單庫存模型各變量及其因果關(guān)系圖如下圖：庫存量庫存差額訂貨量+()R1DI+期望庫存Y2、各變量之間的關(guān)系可用如下方程表示：L IK=IJ+DT*R1JKN I=1000R R1KL=DK/ZA DK=Y-IKC Z=5C Y=60003、要求利用STELLA建立上述庫存模型的流圖，仿真計算并分析結(jié)果三、實驗步驟1、確定水準變量、速率變量、輔助變量、常量及水準變量初值；2、熟悉STELLA軟件操作指導，建立模型的四個基本構(gòu)造塊為：棧（stock）、流（flow）、轉(zhuǎn)換器(converter) 、連接器（connector），設(shè)置仿真參數(shù)（采用默認值）；2、根據(jù)因果關(guān)系圖連接流；3、確定水準方程、速率方程、輔助方程、賦初值方程和常量方程；4、建立模型仿真結(jié)果分析所需的數(shù)據(jù)模塊；5、仿真及結(jié)果分析四、實驗結(jié)果1.簡易庫存模型：2.簡易庫存模型輸出特性示意圖：(1) 表格：.（2）曲線：3.量化分析模型及仿真計算公式：I(t) = I(t - dt) + (R1) * dtINIT I = 1000INFLOWS:R1 = D/ZD = Y-IY = 6000Z = 5五、實驗分析：實驗一通過簡易的流圖模型以及dynamo方程導出了如上所述的簡單庫存系統(tǒng)輸出特性示意圖以及簡單庫存系統(tǒng)仿真計算結(jié)果表。實驗二 二階庫存系統(tǒng)一、 實驗目的1、加深對系統(tǒng)動力學主要要素和基本思想的理解2、學會利用STELLA軟件建立二階反饋系統(tǒng)模型、仿真運行及結(jié)果分析3、分析二階系統(tǒng)的特點二、 實驗要求出庫量R31、二階庫存系統(tǒng)模型各變量及其因果關(guān)系圖如下圖：庫存量入庫量途中存貨量訂貨量庫存差額I（）+（）GR2 +R1+ 2、各變量之間的關(guān)系可用如下方程表示：L GK=GJ+DT*(R1KL-R2JK)L IK=IJ+DTR2JKR R1KL=D/ZR R2KL=G/WR R3KL=200A D=Y-IKC Y=6000C W=10,Z=5C I=1000 ，G=80003、要求利用STELLA建立上述庫存模型的流圖，仿真計算并分析結(jié)果三、實驗步驟1、確定水準變量、速率變量、輔助變量、常量及水準變量初值；2、建立模型的四個基本構(gòu)造塊為：棧（stock）、流（flow）、轉(zhuǎn)換器(converter) 、連接器（connector），設(shè)置仿真參數(shù)（采用默認值）；2、根據(jù)因果關(guān)系圖連接流；3、確定水準方程、速率方程、輔助方程、賦初值方程和常量方程；4、建立模型仿真結(jié)果分析所需的數(shù)據(jù)模塊；5、仿真及結(jié)果分析四、練習 改變各參數(shù)（Z、W及水準變量初始值）的值，觀察水準變量和速率變量變化情況，分析參數(shù)對其仿真結(jié)果的影響五、實驗結(jié)果：1. 二階庫存系統(tǒng)模型:2. 二階庫存系統(tǒng)模型輸出簡易模型示意圖：（1）圖表：（2）曲線：3.量化分析模型及仿真計算公式：G(t) = G(t - dt) + (R1 - R2) * dtINIT G = 8000INFLOWS:R1 = D/ZOUTFLOWS:R2 = G/WI(t) = I(t - dt) + (R2) * dtINIT I = 1000INFLOWS:R2 = G/WD = Y-IW = 10Y = 6000Z = 5六、實驗分析：實驗二是以簡單庫存系統(tǒng)為基礎(chǔ)，繪制了二階負反饋結(jié)構(gòu)模型,通過dynamo方程的量化分析模型及仿真計算,得到如上圖的二階庫存系統(tǒng)輸出特性示意圖及二階庫存系統(tǒng)仿真計算結(jié)果表。實驗三 高校發(fā)展規(guī)模問題一、 實驗目的1、學會根據(jù)具體問題建立模型的因果關(guān)系圖；2、應用STELLA軟件建立流圖統(tǒng)模型、仿真運行及結(jié)果分析；3、能夠根據(jù)數(shù)據(jù)分析結(jié)果解釋具體模型。二、 實驗要求1、 高校發(fā)展規(guī)模問題：教學型高校的在校本科生和教師（S和T）是按照一定的比例而相互增長的。已知某高校現(xiàn)有本科生10000名，且每年以SR的增幅增加，每一名教師可引起本科生人數(shù)增加的速率是1人/年。學校現(xiàn)有教師1500名，每個本科生可引起教師增加的速率（TR）是0.05人/年。用STELLA軟件分析該校未來幾年的發(fā)展規(guī)模，要求：(1) 獨立畫出因果關(guān)系圖和流圖（包括相應的DYNAMO方程）；(2) 對流圖進行仿真，根據(jù)仿真結(jié)果分析未來35年的在校本科生和教師的人數(shù)變化情況。三、 實驗步驟1、確定水準變量、速率變量、輔助變量、常量及水準變量初值；2、建立模型的四個基本構(gòu)造塊為：棧（stock）、流（flow）、轉(zhuǎn)換器(converter) 、連接器（connector），設(shè)置仿真參數(shù)（采用默認值）；2、根據(jù)問題繪制因果關(guān)系圖，并根據(jù)因果關(guān)系圖連接流；3、確定水準方程、速率方程、輔助方程、賦初值方程和常量方程；4、建立模型仿真結(jié)果分析所需的數(shù)據(jù)模塊；5、仿真及結(jié)果分析四、實驗結(jié)果：1.數(shù)據(jù)流圖：2.數(shù)據(jù)表與曲線圖：（1）表：（2）曲線：3.相關(guān)函數(shù)：實驗四 AHP建模訓練一、 實驗目的1、加深對AHP建模步驟的理解；2、學會使用Yaahp軟件二、 實驗要求及步驟1、在Yaahp中建立以下的投資方案選擇模型A2、 根據(jù)判斷矩陣標度定義在Yaahp中建立以下判斷矩陣：A1B1B2B3B1C1C2C3B1153C111/28B21/511/3C2219B31/331C31/81/91B2C1C2C3B3C1C2C3C1147C1147C21/416C21/416C31/71/61C31/71/61結(jié)果：（1）實驗圖：2.計算結(jié)果：3、 計算最優(yōu)投資方案;因為權(quán)重0.47050.40870.1208。對應備選方案分別為興建工廠，購置資產(chǎn)，存入銀行。所以，投資方案最好選擇興建工廠。4、熟悉Yaahp軟件的數(shù)據(jù)輸出功能：標度類型：e(0/5)e(8/5)最終結(jié)果-備選方案權(quán)重興建工廠0.4705購置資產(chǎn)0.4087存入銀行0.1208-1. 投資效果 判斷矩陣一致性比例：0.0000; 對總目標的權(quán)重：1.0000; lambda_max：3.0000-投資效果風險收益流轉(zhuǎn)Wi風險1.00002.22551.49180.4718收益0.44931.00000.67030.2120流轉(zhuǎn)0.67031.49181.00000.3162-2. 風險 判斷矩陣一致性比例：0.0000; 對總目標的權(quán)重：0.4718; lambda_max：3.0000-風險興建工廠購置資產(chǎn)存入銀行Wi興建工廠1.00000.81874.05520.4052購置資產(chǎn)1.22141.00004.95300.4949存入銀行0.24660.20191.00000.0999-3. 收益 判斷矩陣一致性比例：0.0171; 對總目標的權(quán)重：0.2120; lambda_max：3.0178-收益興建工廠購置資產(chǎn)存入銀行Wi興建工廠1.00001.82213.32010.5289購置資產(chǎn)0.54881.00002.71830.3317存入銀行0.30120.36791.00000.1394-4. 流轉(zhuǎn) 判斷矩陣一致性比例：0.0171; 對總目標的權(quán)重：0.3162; lambda_max：3.0178-流轉(zhuǎn)興建工廠購置資產(chǎn)存入銀行Wi興建工廠1.00001.82213.32010.5289購置資產(chǎn)0.54881.00002.71830.3317存入銀行0.30120.36791.00000.1394-附錄：判斷矩陣標度定義實驗五 應用AHP進行超市選址一、 實驗目的1、學會根據(jù)具體問題建立AHP模型；2、學會應用Yaahp軟件進行綜合分析；3、能夠根據(jù)數(shù)據(jù)分析結(jié)果解釋具體問題。二、 實驗要求1、超市選址問題描述：某超市連鎖集團打算在上海D1、青島D2、成都D3三地選擇一個城市進行投資，經(jīng)過專家分析得到影響超市選址的評價指標主要有市場吸引力B1、零售環(huán)境B2和競爭態(tài)勢B3，其中評價市場吸引力的主要有投資成本C1、物流C2和協(xié)同效應C3三個指標，評價零售環(huán)境的主要有市場增長潛力C4、進入壁壘C5和人均零售額C6三個指標，評價競爭態(tài)勢的主要有主要競爭對手C7、市場飽和度C8和網(wǎng)點密度C9三個指標，經(jīng)過專家打分可以確定各判斷矩陣如下表：A1B1B2B3B1C1C2C3B1153C1122B21/511/3C21/212B31/331C31/21/21B2C4C5C6B3C7C8C8C411/32C7147C5315C81/416C61/21/51C91/71/612、要求：(1) 在Yaahp中繪制選址問題的AHP評價模型；(2) 利用Yaahp計算各子指標層在整個評價體系中所占的權(quán)重系數(shù)；(3) 建立超市城市選址的評價體系（評價指標及權(quán)重分配、打分等級劃分）3、如果經(jīng)過專家進一步判斷，得出方案層相對各子指標的評價矩陣如下:C1D1D2D3C2D1D2D3D1127D111/28D21/216D2219D31/71/61D31/81/91C3D1D2D3C4D1D2D3D1148D1147D21/416D21/416D31/81/61D31/71/61C5D1D2D3C6D1D2D3D1147D1122D21/416D21/212D31/71/61D31/21/21C7D1D2D3C8D1D2D3D1153D1137D21/511/3D21/316D31/331D31/71/61C9D1D2D3D1153D21/511/3D31/331利用Yaahp選擇最優(yōu)城市。三、 實驗步驟1、確定超市選址問題的目標層、準則層和方案層；2、在在Yaahp中繪制選址問題的AHP評價模型；3、建立各級判斷矩陣（注意要根據(jù)具體的評價尺度）；4、針對要求2利用Yaahp計算各子指標層在整個評價體系中所占的權(quán)重系數(shù)；5、針對要求3計算方案層各方案相對于目標層的權(quán)重系數(shù)，選擇最優(yōu)方案（若一致性判斷不滿足要求，可局部自行調(diào)整相應的判斷矩陣）。四、實驗結(jié)果：1.AHP圖：2.最終結(jié)果備選方案權(quán)重上海D10.4902青島D20.3400成都D30.16981. 超市選址 判斷矩陣一致性比例：0.0000; 對總目標的權(quán)重：1.0000; lambda_max：3.0000超市選址市場吸引力零售環(huán)境B2競爭態(tài)勢B3Wi市場吸引力1.00002.22551.49180.4718零售環(huán)境B20.44931.00000.67030.2120競爭態(tài)勢B30.67031.49181.00000.31622. 市場吸引力 判斷矩陣一致性比例：0.0043; 對總目標的權(quán)重：0.4718; lambda_max：3.0044市場吸引力投資成本C1物流C2協(xié)同效應C3Wi投資成本C11.00001.22141.22140.3786物流C20.81871.00001.22140.3314協(xié)同效應C30.81870.81871.00000.29003. 零售環(huán)境B2 判斷矩陣一致性比例：0.0043; 對總目標的權(quán)重：0.2120; lambda_max：3.0044零售環(huán)境B2市場增長潛力C4進入壁壘C5人均零售額C6Wi市場增長潛力C41.00000.67031.22140.2976進入壁壘C51.49181.00002.22550.4745人均零售額C60.81870.44931.00000.22794. 競爭態(tài)勢B3 判斷矩陣一致性比例：0.0386; 對總目標的權(quán)重：0.3162; lambda_max：3.0401競爭態(tài)勢B3網(wǎng)點密度C9市場飽和度C8競爭對手C7Wi網(wǎng)點密度C91.00001.82213.32010.5217市場飽和度C80.54881.00003.32010.3497競爭對手C70.30120.30121.00000.12865. 投資成本C1 判斷矩陣一致性比例：0.0000; 對總目標的權(quán)重：0.1786; lambda_max：3.0000投資成本C1上海D1青島D2成都D3Wi上海D11.00001.22143.32010.4717青島D20.81871.00002.71830.3862成都D30.30120.36791.00000.14216. 物流C2 判斷矩陣一致性比例：0.0171; 對總目標的權(quán)重：0.1563; lambda_max：3.0178物流C2上海D1青島D2成都D3Wi上海D11.00001.22144.05520.4648青島D20.81871.00004.95300.4349成都D30.24660.20191.00000.10037. 協(xié)同效應C3 判斷矩陣一致性比例：0.0043; 對總目標的權(quán)重：0.1368; lambda_max：3.0044協(xié)同效應C3上海D1青島D2成都D3Wi上海D11.00001.82214.05520.5503青島D20.54881.00002.71830.3228成都D30.24660.36791.00000.12698. 市場增長潛力C4 判斷矩陣一致性比例：0.0171; 對總目標的權(quán)重：0.0631; lambda_max：3.0178市場增長潛力C4上海D1青島D2成都D3Wi上海D11.00001.82213.32010.5289青島D20.54881.00002.71830.3317成都D30.30120.36791.00000.13949. 進入壁壘C5 判斷矩陣一致性比例：0.0171; 對總目標的權(quán)重：0.1006; lambda_max：3.0178進入壁壘C5上海D1青島D2成都D3Wi上海D11.00001.82213.32010.5289青島D20.54881.00002.71830.3317成都D30.30120.36791.00000.139410. 人均零售額C6 判斷矩陣一致性比例：0.0043; 對總目標的權(quán)重：0.0483; lambda_max：3.0044人均零售額C6上海D1青島D2成都D3Wi上海D11.00001.22141.22140.3786青島D20.81871.00001.22140.3314成都D30.81870.81871.00000.