數(shù)學建模作業(yè)5數(shù)學規(guī)劃模型-供應與選址的問題

數(shù)學建模作業(yè)5數(shù)學規(guī)劃模型----供應與選址的問題數(shù)學建模作業(yè)5數(shù)學規(guī)劃模型----供應與選址的問題數(shù)學建模作業(yè)5數(shù)學規(guī)劃模型----供應與選址的問題數(shù)學建模作業(yè)5數(shù)學規(guī)劃模型----供應與選址的問題編制僅供參考審核批準生效日期地址：電話：傳真:郵編:一、問題提出某公司有6個建筑工地要開工，每個工地的位置（用平面坐標系(a，b)表示，距離單位：km）及水泥日用量d（噸）由下表給出。目前有兩個料場位于A(5,1)，B(2,7)，日儲量各有20噸。試制定每天的供應計劃，即從A，B兩料場分別向各工地運送多少水泥，可使運輸費用（總的噸千米數(shù)）最小，并求出噸千米數(shù)。(注：先畫圖，在坐標上標出各工地位置（用藍色*標示）和料場位置（用紅色o標示）)（2）目前公司準備建立兩個新的料場，日儲量各為20噸，為使運輸費用最省，問新的料場應建在何處，并算出兩料場分別向工地運輸多少噸水泥和費用。（注：初始值取x0=[35471000005115477]’）二、問題分析對于問題（1），確定用A，B兩料場分別向各工地運送水泥，使運輸費用（總的噸千米數(shù)）最小，即要知道兩點間線段最小，料場到工地的路線是直的，而要滿足六個工地的需求，又要考慮到A、B兩個料場的供應量，即在各工地用量必須滿足和各料場運送量不超過日儲量的條件下，使總的噸千米數(shù)最小，這是線性問題。。對于問題（2），需要重新改建六個新的料場，使得在在各工地用量必須滿足和各料場運送量不超過日儲量的條件下，使總的噸千米數(shù)最小，則需要確定新的料場的具體位置，這是非線性問題。三、模型假設(shè)1、假設(shè)料場和建筑工地之間都可以由直線到達；2、運輸費用由“噸千米數(shù)”來衡量；3、兩料場的日存儲量夠向各建筑工地供應；4、運輸途中不發(fā)生意外，從料場運出的水泥總量不會超過各個料場的日存儲量。四、模型建立（顯示模型函數(shù)的構(gòu)造過程）記工地的位置為，水泥日用量為，i=1,…,6;料場位置為，日儲量為，j=1,2；料場向工地的運送量為。目標函數(shù)為：約束條件為：當用臨時料場時決策變量為：當不用臨時料場時決策變量為：，，使用臨時料場的情形：使用兩個臨時料場A(5,1)，B(2,7)．求從料場j向工地的運送量.在各工地用量必須滿足和各料場運送量不超過日儲量的條件下，使總的噸千米數(shù)最小，這是線性規(guī)劃問題。線性規(guī)劃模型為：其中，i=1,2,…,6,j=1,2,為常數(shù)設(shè)X11=X1,X21=X2,,X31=X3,X41=X4,X51=X5,,X61=X6X12=X7,X22=X8,,X32=X9,X42=X10,X52=X11,,X62=X12改建兩個新料場的情形：改建兩個新料場，要同時確定料場的位置(xj,yj)和運送量，在同樣條件下使總噸千米數(shù)最小．這是非線性規(guī)劃問題．非線性規(guī)劃模型為：設(shè)X11=X1,X21=X2,X31=X3,X41=X4,X51=X5,,X61=X6X12=X7,X22=X8,X32=X9,X42=X10,X52=X11,X62=X12x1=X13,y1=X14,x2=X15,y2=X16五、模型求解（顯示模型的求解方法、步驟及運算程序、結(jié)果）建立程序：x=[3];y=[5];x0=[5,2];y0=[1,7];plot(x,y,'*b');holdon;plot(x0,y0,'or');text,,'1¤μ?1');text,,'1¤μ?2');text,,'1¤μ?3')text,5,'1¤μ?4');text(3,,'1¤μ?5');text,,'1¤μ?6')text(5,1,'á?3?A');text(2,7,'á?3?B');使用臨時料場的情形：編寫程序如下：cleara=[3];b=[5];d=[3547611];x=[52];y=[17];e=[2020];fori=1:6forj=1:2aa(i,j)=sqrt((x(j)-a(i))^2+(y(j)-b(i))^2);endendCC=[aa(:,1);aa(:,2)]'A=[111111000000000000111111];B=[20;20];Aeq=[100000100000010000010000001000001000000100000100000010000010000001000001];beq=[d(1);d(2);d(3);d(4);d(5);d(6)];vlb=[000000000000];vub=[];x0=[123010010101];[xx,fval]=linprog(CC,A,B,Aeq,beq,vlb,vub,x0)程序截圖如下：程序的運行結(jié)果為：xx=fval=運行結(jié)果截圖如下：即由料場A、B向6個工地運料方案為：123456料場1350701料場20040610總的噸千米數(shù)為.改建兩個新料場的情形：先編寫M文件：functionf=liaoch(x)a=[3];b=[5];d=[3547611];e=[2020];f1=0;fori=1:6s(i)=sqrt((x(13)-a(i))^2+(x(14)-b(i))^2);f1=s(i)*x(i)+f1;endf2=0;fori=7:12s(i)=sqrt((x(15)-a(i-6))^2+(x(16)-b(i-6))^2);f2=s(i)*x(i)+f2;endf=f1+f2;再編寫主程序為：clearx0=[35471000005115477];A=[11111100000000000000001111110000];B=[20;20];Aeq=[100000100000000001000001000000000010000010000000000100000100000000001000001000000000010000010000];beq=[3547611]';vlb=[zeros(12,1);-inf;-inf;-inf;-inf];vub=[];[x,fval,exitflag]=fmincon('liaoch',x0,A,B,Aeq,beq,vlb,vub)程序截圖如下：程序運行結(jié)果如下：x=Columns1through8000Col