試驗(yàn)四運(yùn)籌學(xué)

1、實(shí)驗(yàn)四網(wǎng)絡(luò)模型的建模與求解一、實(shí)驗(yàn)?zāi)康模涸砗蛻?yīng)掌握最短路問題、最大流問題的建模及求解方法。理解統(tǒng)籌方法的基本概念、 用，熟練運(yùn)用R軟件求解最短路問題、最大流問題以及計(jì)劃網(wǎng)絡(luò)圖相關(guān)問題。 二、實(shí)驗(yàn)內(nèi)容：(1)教材P269習(xí)題2,用R求解最短路問題；(2)用R求解教材P271習(xí)題9的最大流問題；(3)繪制教材P297習(xí)題4的計(jì)劃網(wǎng)絡(luò)圖，并用 R求解習(xí)題5。三、實(shí)驗(yàn)要求：(1)根據(jù)現(xiàn)實(shí)問題建立網(wǎng)絡(luò)模型；(2)繪制計(jì)劃網(wǎng)絡(luò)圖，計(jì)算網(wǎng)絡(luò)圖中的時(shí)間參數(shù)，確定關(guān)鍵路線；(3)運(yùn)用R軟件求解最短路問題、最大流問題以及計(jì)劃網(wǎng)絡(luò)圖相關(guān)問題；(4)對軟件運(yùn)行結(jié)果進(jìn)行分析。四、實(shí)驗(yàn)報(bào)告要求：實(shí)驗(yàn)過程描述(包括分析過

2、程、分析結(jié)果及其解釋、實(shí)驗(yàn)過程遇到的問題及體會)(1)教材P269習(xí)題2,用R求解最短路問題； edge_v library(igraph) edge_v edge_w-c(4,18,12,16,6,2,8,7,5,6)g-graph(edge_v, directed=F)g g#求出最短路 distances(g),1 ,2 ,3 ,4 ,5 ,6 ,7 TOC o 1-5 h z 1,0416182225272,4012141821233,1612026911422186801156,25219711067, 27 23 11 13 5 6 0 shortest_

3、path all_shortest_paths(g,1,to=7)$res$res1+ 5/7 vertices, from 5421d73:1 2 3 5 7$nrgeo1 1 1 1 1 1 1#畫出原圖形，用紅色畫出最短路圖形 edge_color edge_coloras.vector(shortest_path$epath1) layout(matrix(c(1,2),ncol=2, byrow = TRUE) plot(g,edge.label=E(g)$label,edge.label.color=blue) plot(g,edge.label=E(g)$label,edge.l

4、abel.color=blue,edge.color=edge_color)解：路程最短為27分鐘(2)用R求解教材P271習(xí)題9的最大流問題；library(igraph) edge_v-c(1,2,1,3,1,4,2,4,2,5,3,4,3,6,4,5,4,6,5,6)g-graph(edge_v, directed=F)edge_c-c(6,10,6,4,8,6,5,5,6,12)E(g)$capacity-edgecg_mflow-graph.maxflow(g,1,6)#畫出最大流量圖。e_w library(igraph) source(c:/users/Administrator

5、/Desktop/myfunctions.R) edge_v edge_label edge_time-c(2,4,5,4,3,2,4)g-graph(edge_v, directed=T)res k_edges-rep(否”,length(edge_label) ; k_edgesres$EF=res$LF A A-Aorder(A$工序,decreasing=F),;A 最早EF最晚LF關(guān)鍵工序1工序ES LSa 0 426否2b 0 044是3c 6 71112否4d 6 61010是5e 4 7710否6f 11 121314否7g 10 101414是8虛工序4466是e_w,res

6、$LF, ), sep=)edge.color-rep(gray,length(E(g)edge.colorres$EF=res$LF library(igraph)source(c:/users/Administrator/Desktop/myfunctions.R) edge_v-c(1,3,1,2,3,4,3,5,2,5,4,6,5,6,2,3) edge_label-c(a,bTc,d,e,fVg, edge_time-c(2,4,5,4,3,2,4) g-graph(edge_v, directed=T) res k_edges-rep(否”,length(edge_label) ;

7、 k_edgesres$EF=res$LF A-data.frame(工序=edge_label,ES=res$ES,LS=res$LS,=res$EF,最晚 LF=res$LF, 關(guān)鍵工序=k_edges)最早EF工ES；LS最早EF最晚LF關(guān)鍵工序1a0426否2b0044是3c671112否4d661010是5e47710否6f 111121314否7g10101414是8 J晅T4466是工序,decreasing=F),;Ae_w A-Aorder(A$,res$LF, ), sep=)edge.color-rep(gray,length(E(g) edge.colorres$EF=res$LF plot(g, edge.label=e w, ed