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

案例研究 供應(yīng)與選址 某公司有6個(gè)建筑工地要開工 每個(gè)工地的位置 用平面坐標(biāo)系a b表示 距離單位 千米 及水泥日用量d 噸 由下表給出 目前有兩個(gè)臨時(shí)料場(chǎng)位于A 5 1 B 2 7 日儲(chǔ)量各有20噸 假設(shè)從料場(chǎng)到工地之間均有直線道路相連 1 試制定每天的供應(yīng)計(jì)劃 即從A B兩料場(chǎng)分別向各工地運(yùn)送多少噸水泥 使總的噸千米數(shù)最小 2 為了進(jìn)一步減少噸千米數(shù) 打算舍棄兩個(gè)臨時(shí)料場(chǎng) 改建兩個(gè)新的 日儲(chǔ)量各為20噸 問應(yīng)建在何處 節(jié)省的噸千米數(shù)有多大 一 建立模型 記工地的位置為 ai bi 水泥日用量為di i 1 6 料場(chǎng)位置為 xj yj 日儲(chǔ)量為ej j 1 2 從料場(chǎng)j向工地i的運(yùn)送量為Xij 當(dāng)用臨時(shí)料場(chǎng)時(shí)決策變量為 Xij 當(dāng)不用臨時(shí)料場(chǎng)時(shí)決策變量為 Xij xj yj 二 使用臨時(shí)料場(chǎng)的情形 使用兩個(gè)臨時(shí)料場(chǎng)A 5 1 B 2 7 求從料場(chǎng)j向工地i的運(yùn)送量為Xij 在各工地用量必須滿足和各料場(chǎng)運(yùn)送量不超過日儲(chǔ)量的條件下 使總的噸千米數(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ì)的程序?yàn)?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 計(jì)算結(jié)果為 x 3 00005 00000 00007 00000 00001 00000 00000 00004 00000 00006 000010 0000 fval 135 2815 三 改建兩個(gè)新料場(chǎng)的情形 改建兩個(gè)新料場(chǎng) 要同時(shí)確定料場(chǎng)的位置 xj yj 和運(yùn)送量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ī)劃的計(jì)算結(jié)果及臨時(shí)料場(chǎng)的坐標(biāo) x0 35070100406105127 編寫主程序gying2 m 3 計(jì)算結(jié)果為 x 3 00005 00004 00007 00001 0000000005 000011 00005 69594 92857 25007 7500 fval 89 8835 注意初始點(diǎn)的選取 某廠向用戶提供發(fā)動(dòng)機(jī) 合同規(guī)定 第一 二 三季度末分別交貨40臺(tái) 60臺(tái) 80臺(tái) 每季度的生產(chǎn)費(fèi)用為 元 其中x是該季生產(chǎn)的臺(tái)數(shù) 若交貨后有剩余 可用于下季度交貨 但需支付存儲(chǔ)費(fèi) 每臺(tái)每季度c元 已知工廠每季度最大生產(chǎn)能力為100臺(tái) 第一季度開始時(shí)無(wú)存貨 設(shè)a 50 b 0 2 c 4 問工廠應(yīng)如何安排生產(chǎn)計(jì)