編號：時間：2021年x月x日書山有路勤為徑，學(xué)海無涯苦作舟頁碼：第類別倉庫名稱倉庫編號地址面積m2距CPD路程(km)需求量Q單位需求量Q/m2發(fā)貨倉庫CPDCPD上海嘉定區(qū)園汽路1000號42,00001658133.94711101上海嘉定區(qū)民豐路24號3,6006.414680.40721102上海嘉定于塘路379號15,0003.2478233.1831103上海嘉定區(qū)園工路1169號1,0003.97410.7441104上海嘉定區(qū)園國路1366號14,5003.799920.68951105昆山市淀山湖鎮(zhèn)北苑路288號25,000354750.01969106上海嘉定區(qū)泰豐路225號2,1182114965.427非發(fā)貨倉庫11001上海市嘉定區(qū)安亭鎮(zhèn)墨玉北路98號5,0005.38493116.9821002上海市嘉定區(qū)和靜東路318號9,0004.115282016.98根據(jù)具體的情況得出各倉庫同面積、需求量的關(guān)系面積m2需求量QCPD42,00016581311013,6001468110215,0004782311031,000741110414,5009992110525,00047591062,1181149610015,0008493110029,000152820各倉庫同面積、需求量的柱狀圖分析：①由柱狀圖可以得出CPD的面積以及客戶每天的需求量最大；②1105的面積明顯和每天的需求量不協(xié)調(diào)，面積偏大而需求量偏小。③兩個非發(fā)貨倉庫的需求量也偏大，從總體上看：汽車零部件的需求量明顯超過了已有倉庫提供的貨量。=4\*GB3④倉庫1005沒有充分的去利用，資源浪費慮裝卸貨物所用時間的情況下由出發(fā)地到達目的地所用的時間如下表。名稱時間/min距離/km名稱時間/min距離/km1001到110141.81105到CPD64351001到110273.99106到CPD521001到110341.91001到CPD95.31001到110462.61002到CPD124.11001到11054227.61001經(jīng)1101到CPD168.21001到910673.81001經(jīng)1102到CPD127.11002到110152.11001經(jīng)1103到CPD115.81002到110221.31001經(jīng)1104到CPD136.31002到110342.21001經(jīng)1105到CPD10662.61002到1104106.41001經(jīng)9106到CPD125.81002到11053934.81002經(jīng)1101到CPD178.51002到910642.21002經(jīng)1102到CPD74.51101到CPD126.41002經(jīng)1103到CPD115.91102到CPD53.21002經(jīng)1104到CPD1710.11103到CPD73.91002經(jīng)1105到CPD10369.81104到CPD73.71002經(jīng)9106到CPD94.2根據(jù)以上數(shù)據(jù)會出折線圖:分析：①兩個非發(fā)貨倉庫經(jīng)1105到達CPD的時間和距離是最長的；②缺貨必然會引起兩個非發(fā)貨倉庫向發(fā)貨倉庫先發(fā)貨，然后發(fā)話倉庫將貨物發(fā)向CPD，從而延長了零部件配送的時間。建立模型確定最佳存貨量模型根據(jù)需要，對模型提出以下幾點假設(shè)采用（Q，R）策略，對庫存進行連續(xù)性檢查，當(dāng)庫存降低到訂貨點水平R時，即發(fā)出一個訂貨，每次的訂貨量保持不變，都為固定值Q；需求是連續(xù)的，隨機變化的，需要根據(jù)以往的情況進行預(yù)測，需求分布函數(shù)已知；供應(yīng)過程中允許缺貨，但會產(chǎn)生缺貨成本；采購時無價格折扣，但有運輸折扣；在產(chǎn)品存儲過程中會產(chǎn)生現(xiàn)金利息成本；訂購貨物后不能瞬時到貨，需要提前一段時間進行訂購，訂貨提前期是固定不變的；參數(shù)設(shè)置D：企業(yè)對某一產(chǎn)品的年需求量。Q：企業(yè)每次訂購該產(chǎn)品的訂購批量。P：該產(chǎn)品的單位采購價格。k1：該產(chǎn)品的計劃缺貨率，因庫存過高會導(dǎo)致較高的庫存持有成本，而缺貨會產(chǎn)生缺貨成本。S：一次訂購費為S,且每次訂購費都相等。H：單位貨物在倉庫中保存一年的的儲存成本。K2：資金的年利息率。b：單位產(chǎn)品的缺貨成本。n：企業(yè)每年訂購該產(chǎn)品的次數(shù)。f(Q)：訂購量為Q時的單位運輸成本。Q1：企業(yè)的庫存年平均水平。L：訂購提前期，用以減少缺貨的發(fā)生。公式推導(dǎo)在企業(yè)年需求量為D，每次訂購量為Q的條件下，企業(yè)年訂購次數(shù)n的計算公式為：則企業(yè)年訂購成本C1為：假設(shè)企業(yè)需求速率是不變的，則企業(yè)年庫存平均水平Q1為：企業(yè)年庫存存儲成本C2為：企業(yè)年原材料成本C3為：企業(yè)年缺貨成本C4為：企業(yè)年運輸成本C5為：企業(yè)年總成本=年訂購成本+年存儲成本+年原材料成本+年缺貨成本+年運輸成本（4-8）企業(yè)年訂貨量Ｑ越大，存儲成本C2就越高，年運輸成本C5越高，而訂購成本C1越低；反之，訂貨量Q越小，儲存成本C2就越低，年運輸成本也C5越低，但訂貨成本C1越高。企業(yè)需要在增加訂貨量所節(jié)約的訂購成本和運輸成本與增加庫存量所提高的儲存成本之間進行權(quán)衡，以求得兩者的最佳組合。假定價格是固定的，原材料成本C3與批量規(guī)模Q無關(guān)，即不產(chǎn)生規(guī)模經(jīng)濟效益。年缺貨成本C4與訂購批量Q也無關(guān)。隨著批量規(guī)模的增加，年總成本先下降后上升。也就是說訂貨批量存在一個最佳值，使總成本最小，訂貨批量大于或小于這個最佳值都會使總成本增大。最佳批量規(guī)模是使總成本最小的批量規(guī)模。根據(jù)公式（4-8）對總成本求關(guān)于訂貨量Q的一階導(dǎo)數(shù)，并使其為零，就可以得到最佳訂貨批量Q*。則最佳訂貨批量Q*為總結(jié):通過對市場上每種零件需求量的分析科學(xué)的預(yù)測每個倉庫每次從兩個非發(fā)貨倉庫的最佳訂貨單量，從而減少產(chǎn)品的積壓和降低了缺貨帶來的經(jīng)濟損失，提高了經(jīng)濟效益。同時，可以有效地解決倉庫1103,1104,1105和9106均由于補貨不及時或時間太短沒辦法將零部件短駁至CPD倉庫而不能按時完成訂單。保證外庫（發(fā)貨倉庫）與總庫的發(fā)貨同步性，以及控制非發(fā)貨倉庫向發(fā)貨倉的及時補貨?；趘ensim軟件對倉庫監(jiān)控實施的仿真(沒有填寫)4.2.3倉庫貨零部件的重組安吉零部件物流存在的問題：由于市場上根據(jù)人口密度，經(jīng)濟收入，國家政策，季節(jié)性需求等要素造成了各個零部件的需求變化。而實際上安吉物流對每種零部件在每個倉庫的存放是不變的，零部件的不合理存放導(dǎo)致了發(fā)貨倉庫同非發(fā)貨倉庫之間同步性問題，進而帶來經(jīng)濟損失。（1）提出的解決方案:采取倉庫零部件儲存的種類重組（2）采取的原則（考慮的因素）：主要考慮每種零部件的需求量，體積以及零部件的價格。（3）建立的模型由于不知道具體的數(shù)值，所以將零部件的種類設(shè)為i(1～n);ai為i零件的體積；bi為i零件的體積占總零件體積的百分比；xi為i零件的需求量；yi為零件的需求量占總零件需求量的百分比；pi為i零件的價格，qi為i零件價格占總零件價格的百分比；權(quán)重zi=ai*bi+xi*yi+pi*qi；將計算的zi進行從小到大的排序就得到了每種零件重要性的一個排列。根據(jù)越重要的零件存放越靠近CPD倉庫的原則實現(xiàn)存儲的零部件重組。列出表格如下圖：占總需求百分比yii零件的體積ai(m3)占總體積的百分比bii零件的價格pi占總價格的百分比qi權(quán)重zi=ai*bi+xi*yi+pi*qiy0a0b0p0q0z0y1a1b1p1q1z1y2a2b2p2q2z2y3a3b3p3q3z3y4a4b4p4q4z4y5a5b5p5q5z5y6a6b6p6q6z6y7a7b7p7q7z7y8a8b8p8q8z8y9a9b9p9q9z9y10a10b10p10q10z10y11a11b11p11q11z11y12a12b12p12q12z12y13a13b13p13q13z13y14a14b14p14q14z14y15a15b15p15q15z15y16a16b16p16q16z16y17a17b17p17q17z17………………ynanbnpnqnzn程序算法見附件一當(dāng)n=5，ai,bixi,yi,piqi取一定值時計算結(jié)果見附件二采取的具體措施：通過計算的輸入，根據(jù)零件的實際情況，我們很明確的可以知道零件重要性的排序。對需求量比較大的零部件，盡量放在距離CPD倉庫較近的外庫，方便零件的及時供應(yīng)。對于每一個季度，我們可以通過市場調(diào)研，確定零件的參數(shù)指標，對零部件進行更為合理的分配，使經(jīng)濟達到最優(yōu)。（4）針對1105號倉庫提出如下方案：由以上分析可以得出：1105號倉庫離目標地的距離最長，市場的需求不是太大，同時1105號倉庫可以儲存的貨物非常大。實際上根據(jù)市場的需求可以看出，1105號倉庫在很大程度上沒有充分的利用。通過對1005號倉庫合理的劃分出一部分區(qū)域儲存其他種類的汽車零部件，從而減少市場其他倉庫的零件需求過大而倉庫本身的儲存不足問題。解決了倉庫1103,1104,1105和9106均由于補貨不及時或時間太短沒辦法將零部件短駁至CPD倉庫而不能按時完成訂單。4.2.4取貨線路的優(yōu)化安吉汽車零部件以往外庫向CPD倉庫供貨運輸線路以往線路圖：對已有線路的分析：=1\*GB3①1105昆山市淀山湖鎮(zhèn)北苑路288號倉庫離兩個非發(fā)貨中心的距離其他的發(fā)貨倉庫相比較是最遠的，同時1105到達CPD的距離也是所有倉庫中到達CPD中距離最遠的；=2\*GB3②非發(fā)貨倉庫1001和1002向1101,1103,1102,9106四個倉庫發(fā)貨時有相當(dāng)?shù)牡乩韮?yōu)勢，向1104發(fā)貨時就缺少了這種地理優(yōu)勢。1001,1002,1101,1102,1103,9106的地理位置比較集中，有利于進行統(tǒng)一的發(fā)貨管理。=3\*GB3③當(dāng)1104缺貨時，1104要向CPD發(fā)貨時首先由兩個非發(fā)貨倉庫向1104發(fā)貨，1104對貨物接受核對后才向CPD發(fā)貨，無形中延長了1104向CPD發(fā)貨所走的線路。=4\*GB3④如果每一個倉庫用一輛5T的車運輸，線路的重復(fù)性比較高，對小車的利用率比較低，同時增加了運輸?shù)馁M用。=5\*GB3⑤各個外庫距離不等很難滿足同步到達CPD倉庫，進而影響CPD倉庫向客戶配送的效率；線路優(yōu)化的目的總的目標降低運輸?shù)慕?jīng)濟成本和時間成本形成一套完整合理的運輸線路方案，使企業(yè)節(jié)約成本費用具體目標滿足外庫向中心倉庫送貨時的同步性，提高發(fā)單的速度減少5T小車的數(shù)量，滿足運輸要求減少小車的回程空載率運輸線路優(yōu)化模型的建立循環(huán)取貨模型的假設(shè)：假設(shè)一：CPD倉庫為虛擬的配送中心，各個外庫假設(shè)為供應(yīng)商。每個外庫向CPD的補貨數(shù)量定義為客戶的需求量，以逆向的思維去解決。假設(shè)二：運輸成本與距離成正比假設(shè)三：每一項任務(wù)只能由一輛車完成（一對一，一對多）假設(shè)四：不考慮車輛在各個倉庫之間的交接費用假設(shè)五：每個外庫僅被一輛車訪問假設(shè)六：CPD編號為0外庫編號為1到LCPD及外庫均以i循環(huán)取貨模型的建立符號的說明Yki：點i的任務(wù)由車輛k來完成。等于1時，成立。等于0時，不成立。Xijk:車輛k從i點到j(luò)點。等于1時，成立。等于0時不成立。gi：第i個倉庫向CPD倉庫的配貨量。q：小車最大的載重量，與小車本身有關(guān)Cij：從i倉庫到j(luò)倉庫的運輸費用，與距離成正比Sij：從i倉庫到j(luò)倉庫的運輸距離a：小車行駛一公里的費用K:零部件倉庫小車總的數(shù)量數(shù)學(xué)模型Cij=a×Sij模型的求解設(shè)a為一個固定值，已知每個倉庫之間的距離矩陣，運用退化算法對模型進行求解。求解三次。第一次，不考慮1001,1002的情況下，以CPD為虛擬配送中心，進行求解。第二次，不考慮1002，以1001為虛擬配送中心，進行求解。第三次，不考慮1001，以1002為虛擬配送中心，對模型進行求解。求解后的零部件配送流程圖及分析倉庫之間的距離矩陣：CPD11011102110311041105910610011002CPD03.73525.34.111016.40211023.22.8011003.92.21103.906.41105352931.331.434.1034.227.634.891063.81001.92.627.63.803.51002.26.40不考慮1001和1002，以CPD為虛擬配送中心求解程序及結(jié)果見附件五不考慮1002，以1001為虛擬配送中心進行求解程序及結(jié)果見附件六不考慮1001，以1002為虛擬配送中心進行求解程序及結(jié)果見附件七對三次求解的結(jié)果進行匯總，繪制出最基本的運輸流程圖，再次基礎(chǔ)上，根據(jù)實際情況進行求解，以滿足需求，達到最優(yōu)。CPD110411031101110511029106CPD10011104CPD91061102110311011105100110029106CPD1104110211031105110110021104110410021001CPD倉庫10021001CPD倉庫11051101091061102110311051101910611021103紅線表示CPD向其他發(fā)貨倉庫取貨的線路，黑線表示非發(fā)貨倉庫1001向發(fā)貨倉庫補貨的線路，藍線表示非發(fā)貨倉庫1002向發(fā)貨倉庫補貨的線路。采用循環(huán)取貨的方式，提高零部件的同步發(fā)貨。（4）優(yōu)化后的線路與以往線路的對比優(yōu)化后的線路節(jié)約了運輸?shù)某杀?，減少了小車的使用，也盡可能的減少了小車的空載回程率。循環(huán)的取貨方式，滿足了汽車零部件的同步配送理念，節(jié)約了不能同步配送所等待浪費的不必要的時間，使客戶的訂單可以及時的去配發(fā)送，提升了安吉物流公司的配送速度，增強了客戶對安吉物流的信任程度。4.2.5建立新倉庫（1）建設(shè)新倉庫考慮的主要問題倉庫的面積倉庫存儲的零部件中類倉庫的建設(shè)費用倉庫建設(shè)所能達到的效果倉庫的地理位置（2）建設(shè)零部件倉庫的主要目的通過建立新外庫，使外庫向CPD能及時的去配送貨物，解決CPD貨物短缺的問題，使CPD倉庫可以及時去向客戶配送貨物。緩解其他倉庫配送的壓力，均衡倉庫分布，實現(xiàn)各個倉庫之間更為緊密的聯(lián)系。貫穿同步配送的思想，增加倉庫來滿足同步配送的要求新倉庫的備選地址的選擇依據(jù)以下是安吉物流倉庫的分布圖：在進行倉庫備選地址選擇時，考慮的是倉庫的建設(shè)費用，以及新建倉庫與各個倉庫之間的聯(lián)系。在選擇時，我們主要考慮哪些距離CPD倉庫比較近的地區(qū)，方便此倉庫向CPD倉庫零部件的進行運輸，節(jié)約運輸和時間的成本。另外，我們盡量選擇以哪些和其他倉庫運輸線路在一條線上的地區(qū)，這樣可以實現(xiàn)循環(huán)取貨的理念，滿足同步配送的觀念。在此案例中，我們通過分析比較，選擇了五個比較合適建立新倉庫的地區(qū)。包括：百安公路985號，園福路358號，百安公路1558號，園區(qū)路688號，安亭鎮(zhèn)園大路185號。建立新倉庫選擇的模型模型的假設(shè)假設(shè)一：倉庫建設(shè)在上海嘉定區(qū)。假設(shè)二：新建倉庫為發(fā)貨倉庫。假設(shè)三：新建倉庫也配備一輛5T的小車。假設(shè)四：倉庫之間的運輸費用與距離成正比。假設(shè)五：倉庫的容量盡可能足夠的大，可以解決當(dāng)前缺貨的狀況。假設(shè)六：CPD對新建倉庫的需求量已知假設(shè)七：上海嘉定區(qū)的工業(yè)地價相同。假設(shè)八：新建倉庫的單位零部件的存儲費相同。模型符號的說明：n:新建倉庫備選地個數(shù)t:實際擬建倉庫個數(shù)Xi:新建倉庫到達CPD的距離Hi:從新建倉庫到CPD運輸單位距離貨物的運價Gi:從非發(fā)貨倉庫到新建倉庫的距離Ki:從非發(fā)貨倉庫到新建倉庫運輸單位距離貨物的運價Qi：從新建倉庫運到CPD的貨物數(shù)量Di:第i個倉庫存儲貨物的單位存儲費Ei:CPD按照訂單對新建倉庫的需求量Si:新建倉庫的面積Pi:備選地的地價Yi=0，表示不在第i個備選地建倉庫；1，表示在第i個備選地建倉庫數(shù)學(xué)模型說明：模型中的目標函數(shù)是使總成本最小，其中第一項表示從倉庫到CPD的運輸費用，第二項表示新建倉庫的建設(shè)費用，第三項表示倉庫的存儲成本。約束條件項，第一項表示倉庫運到CPD的貨運量要滿足CPD按照訂單對倉庫提供零部件的需求量，第二項表示在備選地中選擇t個建立倉庫，第三項表示變量的取值限制。模型求解對于此數(shù)學(xué)模型，我們引用lingo去求解=1\*GB3①Lingo所編寫的程序如下所示：已知備選地距CPD倉庫的距離、距離百安公路1588嘉定區(qū)園福路358號百安公路985號上海市嘉定區(qū)園區(qū)路688號上海嘉定區(qū)安亭鎮(zhèn)園大路185號CPD倉庫2.1公里1.2公里0.89公里1.2公里2.2公里根據(jù)距離，編寫的程序見附件三求解的結(jié)果，見附件四從結(jié)果中可知，選擇在嘉定區(qū)園福路358號最為合適4.2.6CPD倉庫的重新選擇（1）重新選擇的原因和目的新建倉庫的基礎(chǔ)上，倉庫進行了局部的重新分布，有可能導(dǎo)致CPD倉庫地址的改變。為了進一步去優(yōu)化此方案，我們必須在已有發(fā)貨倉庫的基礎(chǔ)上的重新選擇。選擇出一個最為合適的中心倉庫，以達到最優(yōu)化原則（2）選擇的方法在這里，我們只提供一種較為合理的選擇方法步驟：=1\*GB3①用Floyd算法求出各個倉庫的距離矩陣D=（dij）vv=2\*GB3②計算各倉庫做為CPD倉庫的總運力m（vij）=3\*GB3③求vk使，則vk就是CPD應(yīng)該處在的最佳位置。附頁：附件一#include<stdio.h>voidshuzhu(intn,float*p){inti;for(i=0;i<n;i++)scanf("%f",p+i);}voidshuchu(intn,float*p){inti;for(i=0;i<n;i++)printf("%f",*(p+i));}voidjisuan(float*a,float*b,float*x,float*y,float*p,float*q,float*z,intn){inti;for(i=0;i<n;i++)z[i]=a[i]*b[i]+x[i]*y[i]+p[i]*q[i];}voidpaixu(intn,float*z){for(inti=0;i<n;i++)for(intj=i+1;j<n;j++){if(*(z+i)>*(z+j))//升序{floattemp;temp=*(z+i);*(z+i)=*(z+j);*(z+j)=temp;}}}voidmain(){floata[100];floatb[100];floatx[100];floaty[100];floatp[100];floatq[100];floatz[100];intn;//數(shù)的個數(shù)printf("inputthen:\n");//輸入數(shù)組的大小scanf("%d",&n);printf("inputtheai:\n");//ai為i零件的體積shuzhu(n,&a[0]);printf("inputthebi:\n");//bi為i零件的體積占總零件體積百分比shuzhu(n,&b[0]);printf("inputthexi:\n");//xi為i零件的需求量shuzhu(n,&x[0]);printf("inputtheyi:\n");//yi為i零件的需求占總零件需求的百分比shuzhu(n,&y[0]);printf("inputthepi:\n");//pi為i零件的價格shuzhu(n,&p[0]);printf("inputtheqi:\n");//qi為i零件的價格占總零件價格百分比shuzhu(n,&q[0]);shuchu(n,&a[0]);printf("\n");shuchu(n,&b[0]);printf("\n");shuchu(n,&x[0]);printf("\n");shuchu(n,&y[0]);printf("\n");shuchu(n,&p[0]);printf("\n");shuchu(n,&q[0]);printf("\n");jisuan(&a[0],&b[0],&x[0],&y[0],&p[0],&q[0],&z[0],n);paixu(n,&z[0]);shuchu(n,&z[0]);}附件二：附件三：Lingo求解程序：model:sets:newj/1..5/:newd,newdf,news,newxul,x;!新設(shè)倉庫距離，地價，面積，需求量，0-1規(guī)劃;yuanj/1..9/:d,df,s,xul;!原有倉庫距離，地價，面積，需求量;endsetsmin=@sum(yuanj(i):xul(i))+@sum(newj(k):newxul(k)*x(k))+@sum(newj(k):news(k)*newdf(k)*x(k))+@sum(yuanj(i):s(i)*xul(i))+@sum(newj(k):news(k)*newxul(k)*x(k));@for(newj:@bin(x));@sum(newj:x)<3;@sum(newj:x)>1;!倉庫數(shù)量限制;(@sum(yuanj(i):xul(i))+@sum(newj(k):newxul(k)*x(k)))>675559;!倉庫儲存量限制;(@sum(yuanj(i):s(i))+@sum(newj(k):news(k)))>2000000;!倉庫面積限制;@for(newj(k):news(k)*newdf(k)>100000);!新建倉庫地價費用限制;!(@sum(yuanj(i):s(i)*xul(i))+@sum(newj(k):news(k)*newxul(k)*x(k)))<200000;!倉庫運費限制;data:d=03.7 3525.34.1;df=0.811.2111111;s=420003600150001000145002500021185000 9000;xul=165813 1468 47823 741 9992 475 11496 84931 152820;newd=91.22.2;newdf=0.811.211;enddataend附件四：lingo求解的結(jié)果：VariableValueReducedCostNEWD(1)2.1000000.000000NEWD(2)1.2000000.000000NEWD(3)0.89000000.000000NEWD(4)1.2000000.000000NEWD(5)2.2000000.000000NEWDF(1)0.80000000.000000NEWDF(2)1.0000000.000000NEWDF(3)1.2000000.000000NEWDF(4)1.0000000.000000NEWDF(5)1.0000000.000000NEWS(1)125000.00.000000NEWS(2)100000.00.000000NEWS(3)83333.330.000000NEWS(4)100000.00.000000NEWS(5)1474449.0.000000NEWXUL(1)0.0000000.000000NEWXUL(2)0.00000016666.65NEWXUL(3)0.1439864E+140.000000NEWXUL(4)0.0000000.000000NEWXUL(5)0.0000000.000000X(1)0.000000100000.0X(2)1.000000100000.0X(3)0.1389020E-070.000000X(4)0.000000100000.0附件五：從cpd出發(fā)的路徑程序及結(jié)果clccleard1=[70,40];%距離矩陣dd=[0 6.4 3.2 3.9 3.7 35 2 0;6.4 0 2.8 2.1 4.4 29 4.1 6.4;3.2 2.8 0 1.3 3.2 31.3 1.5 3.2;3.9 2.1 1.3 0 3.9 31.4 2.4 3.9;3.7 4.4 3.2 3.9 0 34.1 3.8 3.7;35 29 31.3 31.4 34.1 0 34.2 35;2 4.1 1.5 2.4 3.8 34.2 0 2;0 6.4 3.2 3.9 3.7 35 2 0;];S0=[];Sum=inf;forj=1:1000S=[11+randperm(6),8];%改數(shù)字（列數(shù)列數(shù)-2temp=0;fori=1:7%改數(shù)字（列數(shù)-1temp=temp+d(S(i),S(i+1));endiftemp<SumS0=S;Sum=temp;endende=0.1^30;L=20000;at=0.999;T=1;%退火過程fork=1:L%產(chǎn)生新解c=2+floor(6*rand(1,2));%改數(shù)字（列數(shù)--2c=sort(c);c1=c(1);c2=c(2);%計算代價函數(shù)值df=d(S0(c1-1),S0(c2))+d(S0(c1),S0(c2+1))-d(S0(c1-1),S0(c1))-d(S0(c2),S0(c2+1));%接受準則ifdf<0S0=[S0(1:c1-1),S0(c2:-1:c1),S0(c2+1:8)];%改數(shù)字（列數(shù)--1Sum=Sum+df;elseifexp(-df/T)>rand(1)S0=[S0(1:c1-1),S0(c2:-1:c1),S0(c2+1:8)];%改數(shù)字（列數(shù)--1Sum=Sum+df;endT=T*at;ifT<ebreak;endend%輸出路徑及路徑長度S0,Sum結(jié)果：S0=15426378Sum=73.5000附件六：從1001出發(fā)的路徑程序及結(jié)果clccleard1=[70,40];%距離矩陣dd=[0 1.8 3.9 1.9 2.6 27.6 3.8 0;1.8 0 2.8 2.1 4.4 29 4.1 1.8;3.9 2.8 0 1.3 3.2 31.3 1.5 3.9;1.9 2.1 1.3 0 3.9 31.4 2.4 1.9;2.6 4.4 3.2 3.9 0 34.1 3.8 2.6;27.6 29 31.3 31.4 34.1 0 34.2 27.6;3.8 4.1 1.5 2.4 3.8 34.2 0 3.8;0 1.8 3.9 1.9 2.6 27.6 3.8 0];S0=[];Sum=inf;forj=1:1000S=[11+randperm(8),10];%改數(shù)字79temp=0;fori=1:9temp=temp+d(S(i),S(i+1));endiftemp<SumS0=S;Sum=temp;endende=0.1^30;L=20000;at=0.999;T=1;%退火過程fork=1:L%產(chǎn)生新解c=2+floor(8*rand(1,2));%改數(shù)字7c=sort(c);c1=c(1);c2=c(2);%計算代價函數(shù)值df=d(S0(c1-1),S0(c2))+d(S0(c1),S0(c2+1))-d(S0(c1-1),S0(c1))-d(S0(c2),S0(c2+1));%接受準則ifdf<0S0=[S0(1:c1-1),S0(c2:-1:c1),S0(c2+1:10)];%改數(shù)字9Sum=Sum+df;elseifexp(-df/T)>rand(1)S0=[S0(1:c1-1),S0(c2:-1:c1),S0(c2+1:10)];%改數(shù)字9Sum=Sum+df;endT=T*at;ifT<ebreak;endend%輸出路徑及結(jié)果：S0=15734268Sum=67.9000附件七：從1002出發(fā)的路徑程序及結(jié)果clccleard1=[70,40];%距離矩陣d%距離矩陣dd=[0 2.1 1.3 2.2 6.4 34.8 2.2 0;2.1 0 2.8 2.1 4.4 29 4.1 2.1;1.3 2.8 0 1.3 3.2 31.3 1.5 1.3;2.2 2.1 1.3 0 3.9 31.4 2.4 2.2;6.4 4.4 3.2 3.9 0 34.1 3.8 6.4;34.8 29 31.3 31.4 34.1 0 34.2 34.8;2.2 4.1 1.5 2.4 3.8 34.2 0 2.2;0 2.1 1.3 2.2 6.4 34.8 2.2 0];S0=[];Sum=inf;forj=1:1000S=[11+randperm(8),10];%改數(shù)字79temp=0;fori=1:9temp=temp+d(S(i),S(i+1));endiftemp<SumS0=S;Sum=temp;endende=0.1^30;L=20000;at=0.999;T=1;%退火過程fork=1:L%產(chǎn)生新解c=2+floor(8*rand(1,2));%改數(shù)字7c=sort(c);c1=c(1);c2=c(2);%計算代價函數(shù)值df=d(S0(c1-1),S0(c2))+d(S0(c1),S0(c2+1))-d(S0(c1-1),S0(c1))-d(S0(c2),S0(c2+1));%接受準則ifdf<0S0=[S0(1:c1-1),S0(c2:-1:c1),S0(c2+1:10)];%改數(shù)字9Sum=Sum+df;elseifexp(-df/T)>rand(1)S0=[S0(1:c1-1),S0(c2:-1:c1),S0(c2+1:10)];%改數(shù)字9Sum=Sum+df;endT=T*at;ifT<ebreak;endend%輸出路徑及結(jié)果：S0=17534628Sum=73.0000附錄資料：不需要的可以自行刪除Word常用快捷鍵由于Word有定義快捷鍵的功能，在重新指定快捷鍵后，不能再使用該組合完成以前的操作。例如，按快捷鍵Ctrl+B可將選定文本改為加粗格式，如果將Ctrl+B重新指定給一個新的命令或其他命令，則不能通過按Ctrl+B為文本應(yīng)用加粗格式，除非將快捷鍵指定恢復(fù)到初始設(shè)置。下面的列表使用說明如下：（1）表中出現(xiàn)的（F）、（E）、（V）、（I）、（O）、（T）、（A）、（W）：分別代表文件、編輯、視圖、插入、格式、工具、表格和窗口菜單。（2）（切換方式）：表示重復(fù)按該鍵還原操作。（3）Num？：表示小鍵盤上的數(shù)字鍵。（4）Ctrl＋A：表示Ctrl鍵和A鍵的組合，其余類推。A.1常規(guī)快捷鍵Ctrl＋A（E）選取整篇文檔（Ctrl＋Num5）Ctrl+R可使段落右對齊Ctrl＋B（E）加粗文本（Ctrl＋Shift＋B）Ctrl＋]（O）按磅值增加所選定內(nèi)容的字號Ctrl＋C（E）復(fù)制所選內(nèi)容Ctrl＋[（O）按磅值縮小所選定內(nèi)容的字號Ctrl＋D（O）修改選定字符格式Ctrl＋=（O）將選定內(nèi)容設(shè)為下標(切換方式)Ctrl＋E（O）段落居中Ctrl＋＋（O）將選定內(nèi)容設(shè)為上標(切換方式)Ctrl＋H（E）查找并修改指定文字或格式Ctrl＋＊（V）顯示/隱藏所有非打印字符Ctrl＋I（O）傾斜所選文字（切換方式）Ctrl＋Tab（A）在單元格中插入制表符Ctrl＋M（F）調(diào)整整段縮進Ctrl+L（O）段落左對齊Ctrl＋N（F）創(chuàng)建新文檔或模板Ctrl＋Enter（O）在插入點插入一個分頁符Ctrl＋O（F）打開已有的文檔或模板Ctrl＋Up（E）將插入點上移一個段落Ctrl＋P（F）打印文檔（=Ctrl＋Shift＋F12）Ctrl＋Down（E）將插入點下移一個段落Ctrl+Q（O）刪除段落格式Ctrl＋Home（E）將插入點移到文檔開始Ctrl＋S（F）保存當(dāng)前活動文檔Ctrl＋End（E）將插入點移到文檔結(jié)尾Ctrl＋T（O）設(shè)置懸掛式縮進Alt＋/（A）快速選定整個表格Ctrl＋U（O）給所選內(nèi)容添加下劃線Shift＋F1（W）有關(guān)命令、屏幕區(qū)域的幫助信息或文字屬性Ctrl＋V（E）在插入點插入剪貼板內(nèi)容Ctrl+Alt+I打印預(yù)覽Ctrl+W關(guān)閉文檔Ctrl+K（E）插入超級鏈接Ctrl＋X（E）剪切所選內(nèi)容并將其放入剪貼板Ctrl+Shift+Space-bar創(chuàng)建不間斷空格Ctrl＋Y（E）重復(fù)上一步操作（=F4=Alt＋Enter）Ctrl+連字符創(chuàng)建不間斷連字符Ctrl＋Z（E）取消上一步操作（=Alt＋Backspace）Ctrl+Spacebar刪除字符格式A.2Ctrl與光標鍵聯(lián)合使用的快捷鍵Ctrl+←可將插入點移到上一個英文單詞或中文詞Ctrl+→可將插入點移到下一個英文單詞或中文詞Ctrl+↑可將插入點移到上一個段落Ctrl+↓可將插入點移到下一個段落A.3Ctrl與部分數(shù)字聯(lián)合使用的快捷鍵Ctrl+0可迅速在段前增加或刪除12磅的行距Ctrl+1改變段落的行距為單倍行距Ctrl+5改變段落的行距為1.5倍行距Ctrl+2改變段落的行距為雙倍行距A.4Ctrl與Shift+字母鍵聯(lián)合使用的快捷鍵Ctrl+Shift+P再通過按“↑”和“↓”鍵改變字號Ctrl+Shift+>增大所選文字的字號Ctrl+Shift+<減小所選文字的字號Ctrl+Shift+A將所有選定的字母設(shè)為大寫Ctrl+Shift+C復(fù)制選定內(nèi)容的格式Ctrl+Shift+D段落最后一行文字占滿全行Ctrl+Shift+E打開或關(guān)閉任務(wù)欄上的“修訂”標記Ctrl+Shift+F再通過按“↑”和“↓”鍵改變字體Ctrl+Shift+H應(yīng)用隱藏文字格式Ctrl+Shift+K將所選字母中的小寫字母設(shè)成大寫，但字母字體大小不變Ctrl+Shift+L應(yīng)用“列表”樣式Ctrl+Shift+M取消左側(cè)段落縮進Ctrl+Shift+N選定段落變?yōu)椤罢摹睒邮紺trl+Shift+Q將所選英文字母變?yōu)镾ymbol字體Ctrl+Shift+S再通過按“↑”和“↓”鍵改變樣式Ctrl+Shift+T減小懸掛縮進量Ctrl+Shift+V對選定內(nèi)容粘貼格式Ctrl+Shift+W只給字、詞加下劃線，不給空格加下劃線Ctrl+Shift+*顯示非打印字符Alt+Shift+D插入“日期”域Alt+Ctrl+L插入ListNum域Alt+Shift+P插入頁數(shù)域Alt+Shift+T插入時間域A.5Ctrl與Alt+字母鍵聯(lián)合使用的快捷鍵Ctrl+Alt+D在文章末尾插入尾注Ctrl+Alt+F在當(dāng)前頁插入尾注Ctrl+Alt+K啟動“自動套用格式”