交通系統(tǒng)分析課程設(shè)計(jì)

1、<p><b>　　目 錄</b></p><p><b>　　1 線性規(guī)劃1</b></p><p>　　1.1 模型及分析1</p><p>　　1.2 Matlab求解方法2</p><p>　　1.3 Lingo求解方法3</p><p><b

2、>　　2 運(yùn)輸規(guī)劃5</b></p><p>　　2.1 模型及分析5</p><p>　　2.2 Lingo求解方法6</p><p><b>　　3 整數(shù)規(guī)劃8</b></p><p>　　3.1 模型及分析9</p><p>　　3.2 LINGGO求解方法9&

3、lt;/p><p>　　4 圖與網(wǎng)絡(luò)分析11</p><p>　　4.1 模型及分析11</p><p>　　4.2 Matlab求解方法11</p><p><b>　　5 預(yù)測(cè)分析14</b></p><p>　　5.1 貨運(yùn)量預(yù)測(cè)14</p><p>　　5.1

4、.1 模型及分析14</p><p>　　5.1.2 R軟件求解方法14</p><p>　　5.1.3 Excel求解方法15</p><p>　　5.2 綜合客運(yùn)量預(yù)測(cè)17</p><p>　　5.2.1模型及分析17</p><p>　　5.2.2用Excel里的模型求解17</p>&

5、lt;p><b>　　6參考文獻(xiàn)19</b></p><p><b>　　1 線性規(guī)劃</b></p><p>　　某地段的地面剖面圖如圖1所示（折線ABCD），擬在AD之間修建一條公路。修筑公路除一般的建造費(fèi)用外，由于填挖土方不平衡而需要增加的額外費(fèi)用為元/m3 ，其中為填挖不平衡土方量（公路填挖寬度為10m）；由于縱坡而引起汽車(chē)額外的

6、油料費(fèi)用（設(shè)計(jì)年限內(nèi)的總費(fèi)用）為元/m，其中i為縱坡度。問(wèn)如何設(shè)計(jì)縱坡才能使這些附加的費(fèi)用為最少？</p><p>　　要求最大縱坡不大于10%，并且。因坡度不大，公路長(zhǎng)度可按水平距離計(jì)算，即。</p><p>　　圖1 某路段的地面線高程</p><p><b>　　1.1 模型及分析</b></p><p>　　原問(wèn)

7、題可用如下的數(shù)學(xué)模型來(lái)表達(dá)：</p><p>　　當(dāng)時(shí)，則目標(biāo)函數(shù)為：</p><p>　　這時(shí)，需增加一個(gè)附加約束條件： </p><p><b>　　所以數(shù)學(xué)模型為：</b></p><p>　　該問(wèn)題為線形規(guī)劃問(wèn)題，為求得最優(yōu)解，可用MATLAB和LINGO求解。</p><p>　　1.2

8、 Matlab求解方法</p><p>　　將上述列出的數(shù)學(xué)模型轉(zhuǎn)成標(biāo)準(zhǔn)模型，如下所示。</p><p>　　用命令：[x，fval]= =linprog（z，A，b，A1，b1，LB，UB）在MATLAB中求解。編寫(xiě)M文件如下：（如圖2所示）</p><p>　　z=[30000,18000];</p><p>　　A=[1,0;0,-1;

9、1,-1; -1,0;0,1;-1,1;-1,-1];</p><p>　　b=[90;-10;40;-50;50;0;-120];</p><p>　　A1=[]; b1=[];</p><p>　　LB=[0;0]; UB=[];</p><p>　　[x,fval]=linprog(z,A,b,A1,b1,LB,UB)</p>

10、;<p>　　圖2 MATLAB求解結(jié)果</p><p>　　由于MATLAB軟件不能代入計(jì)算常數(shù)項(xiàng)，所以用3000000-2880000=120000（元），得到最優(yōu)解為：</p><p><b>　　，， </b></p><p>　　1.3 Lingo求解方法</p><p>　　在模型窗口中輸入如下

11、代碼：</p><p>　　min=-2880000+30000*x1+18000*x2;</p><p><b>　　x1<=90;</b></p><p><b>　　x2>=10;</b></p><p>　　x1-x2<=40;</p><p>&l

12、t;b>　　x1>=50;</b></p><p><b>　　x2<=50;</b></p><p><b>　　x1-x2>=0;</b></p><p>　　x1+x2>=120;</p><p><b>　　x1>=0;</b&

13、gt;</p><p><b>　　x2>=0;</b></p><p>　　輸入過(guò)程和計(jì)算結(jié)果見(jiàn)下圖3和圖4.</p><p>　　圖3 LINGO輸入過(guò)程</p><p>　　圖4 LINGO計(jì)算結(jié)果</p><p><b>　　2 運(yùn)輸規(guī)劃</b></p&g

14、t;<p>　　假設(shè)某平衡物資問(wèn)題有三個(gè)產(chǎn)地(i=1，2，3)和四個(gè)銷(xiāo)地(j=1,2,3，,4)，始點(diǎn)需要運(yùn)出的物資量為、終點(diǎn)需要此物資的總量為;及各產(chǎn)銷(xiāo)點(diǎn)之間的運(yùn)輸費(fèi)用單價(jià)如表2所示，出行總量。試求系統(tǒng)運(yùn)輸費(fèi)用最小的運(yùn)輸費(fèi)用方案(i=1，2，3，4)。</p><p>　　表1 各OD點(diǎn)間出行時(shí)耗表 </p><p&g

15、t;<b>　　2.1 模型及分析</b></p><p>　　在平衡物資運(yùn)輸?shù)难芯恐?，?jīng)常遇到這樣的分配問(wèn)題。設(shè)，，…，為物資產(chǎn)地，相應(yīng)地，，…,相應(yīng)的物資運(yùn)出量。，，…，為物資銷(xiāo)地，，,…，為需要此物資的總量?？偟倪\(yùn)輸量為N。那么，設(shè)從產(chǎn)地到銷(xiāo)地的運(yùn)輸量為，運(yùn)輸費(fèi)用為，則總的運(yùn)輸費(fèi)用為: ?，F(xiàn)在的問(wèn)題是如何分配運(yùn)量使得總的運(yùn)輸費(fèi)用為最少。即找出，滿足 </p>

16、<p>　　(i=1,2,…，m； j=1,2,…,n)</p><p>　?。╥=1,2，…，m）</p><p>　?。╦=1,2，…，n）</p><p><b>　　且使最小。</b></p><p>　　2.2 Lingo求解方法</p><p><b>　　(1

17、)程序</b></p><p><b>　　sets:</b></p><p>　　row/1,2,3/:a;</p><p>　　arrange/1,2,3,4/:b;</p><p>　　link(row,arrange): c,x;</p><p><b>　　ends

18、ets</b></p><p><b>　　data:</b></p><p><b>　　a=5,4,6;</b></p><p>　　b=4,2,3,6;</p><p>　　c=6,22,5,6,</p><p><b>　　3,10,4,8,&l

19、t;/b></p><p><b>　　1,8,2,1;</b></p><p><b>　　enddata</b></p><p>　　[OBJ]min=@sum(link(i,j): c(i,j)*x(i,j));</p><p>　　@for(row(i): @sum(arrange(j

20、):x(i,j))=a(i););</p><p>　　@for(arrange(j): @sum(row(i):x(i,j))=b(j););</p><p>　　@for(link(i,j):x(i,j)>=0;);</p><p><b>　　End</b></p><p>　　在模型窗口中輸入上述代碼，然后

21、點(diǎn)擊工具條上的solve按鈕即可。如圖5所示：</p><p>　　圖5 運(yùn)輸規(guī)劃模型LINGO程序圖</p><p><b>　　(2)計(jì)算結(jié)果</b></p><p>　　由上述過(guò)程解得該系統(tǒng)最小總運(yùn)輸費(fèi)為59，如圖6所示。</p><p>　　圖6 運(yùn)輸規(guī)劃模型LINGO總運(yùn)輸費(fèi)用圖</p><

22、p>　　由圖7可看出最優(yōu)系統(tǒng)相應(yīng)的分配情況是:從到的出行量為2，從到的出行量是3；從到的出行量是2，從到為2；從到的出行量為6，其余始點(diǎn)到終點(diǎn)的出行量均為0。</p><p>　　圖7 運(yùn)輸規(guī)劃模型交通分配圖</p><p><b>　　3 整數(shù)規(guī)劃</b></p><p>　　現(xiàn)用集裝箱托運(yùn)甲、乙兩種貨物，每箱的體積、質(zhì)量、可獲得利潤(rùn)及

23、托運(yùn)所受限制見(jiàn)表2。問(wèn)兩種貨物各托運(yùn)都少箱可獲利最大？</p><p>　　表2 每箱貨物的體積、質(zhì)量、可獲利潤(rùn)及托運(yùn)所受限制</p><p><b>　　3.1 模型及分析</b></p><p>　　設(shè)、分別為甲、乙兩種貨物的托運(yùn)箱數(shù)，則此問(wèn)題的線性規(guī)劃數(shù)學(xué)模型為：</p><p>　　3.2 LINGGO求解方法&

24、lt;/p><p><b>　　(1)程序</b></p><p><b>　　sets:</b></p><p>　　num_i/1,2/:b;</p><p>　　num_j/1,2/:x,c;</p><p>　　link(num_i,num_j):a;</p>

25、<p><b>　　endsets</b></p><p><b>　　data:</b></p><p><b>　　b=20,8;</b></p><p><b>　　c=12,9;</b></p><p>　　a=4,5,2,1;<

26、;/p><p><b>　　enddata</b></p><p>　　[OBJ]max=@sum(num_j(j):c(j)*x(j));</p><p>　　@for(num_i(i):@sum(num_j(j):a(i,j)*x(j))<=b(i););</p><p>　　@for(num_j(j):x(j)&

27、gt;=0;);</p><p>　　@for(num_j(j):@gin(x(j)););</p><p><b>　　End</b></p><p>　　在模型窗口中輸入上述代碼，然后點(diǎn)擊工具條上的solve按鈕即可。如圖8所示:</p><p>　　圖8 整數(shù)規(guī)劃模型LINGO程序圖</p><

28、p><b>　?。?）計(jì)算結(jié)果</b></p><p>　　從圖9中可得，甲貨物的托運(yùn)箱數(shù)為4，乙貨物的托運(yùn)箱數(shù)0，則可獲得的最大利潤(rùn)為48百元，如圖9所示。</p><p>　　圖9 整數(shù)規(guī)劃模型LINGO總獲利圖</p><p><b>　　4 圖與網(wǎng)絡(luò)分析 </b></p><p>　　

29、某地區(qū)七個(gè)城鎮(zhèn)間的公路交通網(wǎng)如圖10所示，試用標(biāo)號(hào)法計(jì)算從A城到G城的最短路線。圖中弧旁數(shù)據(jù)為公路長(zhǎng)度。</p><p>　　圖10 某地區(qū)七個(gè)城鎮(zhèn)間的公路交通網(wǎng)</p><p><b>　　4.1 模型及分析</b></p><p>　　最短路問(wèn)題可借助于距離矩陣求解，先構(gòu)造一個(gè)距離矩陣D：</p><p>　　D中的

30、元素定義如下：</p><p>　　依題意可得，距離矩陣為：</p><p>　　4.2 Matlab求解方法</p><p><b>　?。?）程序</b></p><p>　　新建M-file，在窗口中輸入以下代碼：如圖11所示，然后保存文件至默認(rèn)文件夾。</p><p><b>　

31、　圖11 代碼輸入圖</b></p><p><b>　?。?）計(jì)算結(jié)果</b></p><p>　　最后在Command Window窗口輸入以下代碼，如圖12所示。</p><p>　　>>a=[0,6,3,inf,inf,inf,inf;inf,0,inf,8,7,inf,inf;inf,2,0,8,inf,12,

32、inf;inf,inf,inf,0,2,4,inf;inf,inf,inf,inf,0,inf,3;inf,inf,inf,inf,3,0,6;inf,inf,inf,inf,inf,inf,0];</p><p>　　>> [long,path]=floyd(a,1,7)</p><p><b>　　圖12代碼輸入窗口</b></p>&l

33、t;p>　　則自動(dòng)彈出結(jié)果，結(jié)果顯示：A到G的最短路長(zhǎng)為15，最短路線為：，如圖13所示。 </p><p><b>　　圖13 計(jì)算結(jié)果</b></p><p><b>　　5 預(yù)測(cè)分析</b></p><p><b>　　5.1 貨運(yùn)量預(yù)測(cè)</b></p><p>

34、　　某市社會(huì)總產(chǎn)值與貨運(yùn)量之間有線性相關(guān)關(guān)系，見(jiàn)表3.試建立數(shù)學(xué)模型，并預(yù)測(cè)當(dāng)該市社會(huì)總產(chǎn)值達(dá)60億元時(shí)，該市的貨運(yùn)量是多少。</p><p>　　表3 某社會(huì)總產(chǎn)值與貨運(yùn)量之間關(guān)系 </p><p>　　5.1.1 模型及分析</p><p>　　根據(jù)題意可知，社會(huì)總產(chǎn)值（億元）與貨運(yùn)量（千萬(wàn)t）存在相關(guān)關(guān)系，可以用

35、二元線性回歸方程進(jìn)行分析?？山⒎匠倘缦拢?lt;/p><p>　　式中：X——社會(huì)總產(chǎn)值；</p><p><b>　　Y——貨運(yùn)量；</b></p><p>　　a,b——參數(shù)。</p><p>　　可用R軟件和Excel計(jì)算回歸方程中的系數(shù)，求解過(guò)程如下。</p><p>　　5.1

36、.2 R軟件求解方法</p><p>　　（1）在R軟件中編寫(xiě)如下圖14中程序(第一塊紅色部分)，運(yùn)行計(jì)算，得到a=-53.4341，b=1.9385。</p><p>　?。?）再在R軟件中編寫(xiě)如下代碼，進(jìn)行求解當(dāng)社會(huì)總產(chǎn)值x=60時(shí)的貨運(yùn)量。</p><p><b>　　編程：</b></p><p>　　new&l

37、t;-data.frame(X=60) lm.pred<-predict(lm.sol,new,interval="prediction",level=0.95)</p><p><b>　　lm.pred</b></p><p>　　按回車(chē)鍵運(yùn)行，由圖14得出當(dāng)社會(huì)總產(chǎn)值X=60時(shí)：貨運(yùn)量Y為62.87392千萬(wàn)

38、t。</p><p>　　95%置信度的的置信區(qū)間為，且變量X與Y之間存在高度相關(guān)的關(guān)系。</p><p>　　圖14 R軟件運(yùn)行結(jié)果圖</p><p>　　5.1.3 Excel求解方法</p><p><b>　　求解過(guò)程如下：</b></p><p>　?。?）在Excel表格中輸入原始數(shù)據(jù)

39、（如圖15所示）</p><p>　?。?）依次點(diǎn)擊“工具”，“數(shù)據(jù)分析”，“回歸”，彈出如圖16所示選項(xiàng)框后進(jìn)行編輯</p><p>　?。?）點(diǎn)擊確定得到分析結(jié)果（如圖17所示）</p><p><b>　　圖15 原始數(shù)據(jù)</b></p><p><b>　　圖16 回歸命令圖</b><

40、/p><p><b>　　圖17 回歸統(tǒng)計(jì)圖</b></p><p>　　如圖17所示，有, , </p><p><b>　　故</b></p><p>　　即當(dāng)該市社會(huì)總產(chǎn)值達(dá)60億元時(shí)，該市的貨運(yùn)量是59.2957千萬(wàn)。</p><p>　　5.2 綜合客運(yùn)量預(yù)測(cè)</

41、p><p>　　某地區(qū)公路網(wǎng)規(guī)劃中需預(yù)測(cè)2010年的綜合客運(yùn)量，現(xiàn)調(diào)查收集到該地區(qū)1981--2000年綜合客運(yùn)量數(shù)據(jù)如表7-16所示，根據(jù)上述條件預(yù)測(cè)該地區(qū)2010年綜合客運(yùn)量。</p><p>　　表4 某地區(qū)歷年綜合客運(yùn)量（萬(wàn)人次/年） </p><p>　　5.2.1模型及分析</p><p>　　通過(guò)

42、對(duì)上表的數(shù)據(jù)分析發(fā)現(xiàn)，綜合客運(yùn)量的隨著時(shí)間的推移呈現(xiàn)總體增加的趨勢(shì)。所以，根據(jù)歷史統(tǒng)計(jì)資料可以以時(shí)間為自變量建立時(shí)間序列模型，對(duì)未來(lái)綜合客運(yùn)量進(jìn)行預(yù)測(cè)。</p><p>　　該模型屬于時(shí)間序列法當(dāng)中的趨勢(shì)外推法，該方法一般包括以下六個(gè)階段：</p><p><b>　?。?）數(shù)據(jù)收集；</b></p><p>　?。?）選擇預(yù)測(cè)趨勢(shì)線的函數(shù)類(lèi)型

43、；</p><p><b>　?。?）擬合曲線；</b></p><p><b>　?。?）趨勢(shì)外推；</b></p><p>　?。?）預(yù)測(cè)結(jié)果分析及說(shuō)明；</p><p>　?。?）研究預(yù)測(cè)結(jié)果在決策和規(guī)劃中的應(yīng)用。</p><p>　　5.2.2用Excel里的模型求解

44、</p><p><b>　　求解過(guò)程如下：</b></p><p>　?。?）輸入數(shù)據(jù)（如圖18所示）</p><p>　?。?）依次點(diǎn)擊“插入”，“圖表”，“散點(diǎn)圖”</p><p>　　（3）加入趨勢(shì)線（如圖18所示）</p><p>　　圖18 指數(shù)時(shí)間序列模型求解圖</p>

45、<p>　　由圖18可知，綜合客運(yùn)量隨時(shí)間變化的指數(shù)時(shí)間序列模型為： </p><p>　　式中：y——綜合客運(yùn)量；</p><p>　　x——時(shí)間序列（以1981年作為時(shí)間序列起點(diǎn)）</p><p>　　將2010年所對(duì)應(yīng)的時(shí)間序列（30）代入上述預(yù)測(cè)模型，計(jì)算得到2010年該地區(qū)的全社會(huì)綜合客運(yùn)量為：萬(wàn)人次</p><p>　

46、　所以，預(yù)測(cè)該地區(qū)2010年綜合客運(yùn)量為31216萬(wàn)人次。</p><p><b>　　6參考文獻(xiàn)</b></p><p>　　[1]王煒等. 道路交通工程系統(tǒng)分析方法. 北京：人民交通出版社，2004.</p><p>　　[2]王沫然. Matlab與科學(xué)計(jì)算. 北京：電子工業(yè)出版社，2005.</p><p>