Matlab供應(yīng)與選址問題(附詳細(xì)編程).ppt

案例研究 供應(yīng)與選址 某公司有6個建筑工地要開工 每個工地的位置 用平面坐標(biāo)系a b表示 距離單位 千米 及水泥日用量d 噸 由下表給出 目前有兩個臨時料場位于A 5 1 B 2 7 日儲量各有20噸 假設(shè)從料場到工地之間均有直線道路相連 1 試制定每天的供應(yīng)計劃 即從A B兩料場分別向各工地運送多少噸水泥 使總的噸千米數(shù)最小 2 為了進(jìn)一步減少噸千米數(shù) 打算舍棄兩個臨時料場 改建兩個新的 日儲量各為20噸 問應(yīng)建在何處 節(jié)省的噸千米數(shù)有多大 一 建立模型 記工地的位置為 ai bi 水泥日用量為di i 1 6 料場位置為 xj yj 日儲量為ej j 1 2 從料場j向工地i的運送量為Xij 當(dāng)用臨時料場時決策變量為 Xij 當(dāng)不用臨時料場時決策變量為 Xij xj yj 二 使用臨時料場的情形 使用兩個臨時料場A 5 1 B 2 7 求從料場j向工地i的運送量為Xij 在各工地用量必須滿足和各料場運送量不超過日儲量的條件下 使總的噸千米數(shù)最小 這是線性規(guī)劃問題 線性規(guī)劃模型為 設(shè)X11 X1 X21 X2 X31 X3 X41 X4 X51 X5 X61 X6X21 X7 X22 X8 X32 X9 X42 X10 X52 X11 X62 X12編寫程序gying1 m 其詳細(xì)的程序為 c1 sqrt 5 1 25 2 1 1 25 2 c2 sqrt 5 8 75 2 1 0 75 2 c3 sqrt 5 0 5 2 1 4 75 2 c4 sqrt 5 5 75 2 1 5 2 c5 sqrt 5 3 2 1 6 5 2 c6 sqrt 5 7 25 2 1 7 25 2 c7 sqrt 2 1 25 2 7 1 25 2 c8 sqrt 2 8 75 2 7 0 75 2 c9 sqrt 2 0 5 2 7 4 75 2 c10 sqrt 2 5 75 2 7 5 2 c11 sqrt 2 3 2 7 6 5 2 c12 sqrt 2 7 25 2 7 7 25 2 c c1 c2 c3 c4 c5 c6 c7 c8 c9 c10 c11 c12 A 111111000000 000000111111 B 20 20 Aeq 100000100000 010000010000 001000001000 000100000100 000010000010 000001000001 Beq 3 5 4 7 6 11 vlb 0 0 0 0 0 0 0 0 0 0 0 0 x f linprog c A B Aeq Beq vlb 計算結(jié)果為 x 3 00005 00000 00007 00000 00001 00000 00000 00004 00000 00006 000010 0000 fval 135 2815 三 改建兩個新料場的情形 改建兩個新料場 要同時確定料場的位置 xj yj 和運送量Xij 在同樣條件下使總噸千米數(shù)最小 這是非線性規(guī)劃問題 非線性規(guī)劃模型為 設(shè)X11 X1 X21 X2 X31 X3 X41 X4 X51 X5 X61 X6X21 X7 X22 X8 X32 X9 X42 X10 X52 X11 X62 X12x1 X13 y1 X14 x2 X15 y2 X16 1 先編寫M文件liaoch m定義目標(biāo)函數(shù) functiony liaoch x y x 1 sqrt x 13 1 25 2 x 14 1 25 2 x 2 sqrt x 13 8 75 2 x 14 0 75 2 x 3 sqrt x 13 0 5 2 x 14 4 75 2 x 4 sqrt x 13 5 75 2 x 14 5 2 x 5 sqrt x 13 3 2 x 14 6 5 2 x 6 sqrt x 13 7 25 2 x 14 7 25 2 x 7 sqrt x 15 1 25 2 x 16 1 25 2 x 8 sqrt x 15 8 75 2 x 16 0 75 2 x 9 sqrt x 15 0 5 2 x 16 4 75 2 x 10 sqrt x 15 5 75 2 x 16 5 2 x 11 sqrt x 15 3 2 x 16 6 5 2 x 12 sqrt x 15 7 25 2 x 16 7 25 2 2 取初值為線性規(guī)劃的計算結(jié)果及臨時料場的坐標(biāo) x0 35070100406105127 編寫主程序gying2 m 3 計算結(jié)果為 x 3 00005 00004 00007 00001 0000000005 000011 00005 69594 92857 25007 7500 fval 89 8835 注意初始點的選取 某廠向用戶提供發(fā)動機 合同規(guī)定 第一 二 三季度末分別交貨40臺 60臺 80臺 每季度的生產(chǎn)費用為 元 其中x是該季生產(chǎn)的臺數(shù) 若交貨后有剩余 可用于下季度交貨 但需支付存儲費 每臺每季度c元 已知工廠每季度最大生產(chǎn)能力為100臺 第一季度開始時無存貨 設(shè)a 50 b 0 2 c 4 問工廠應(yīng)如何安排生產(chǎn)計