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

案例研究 供應(yīng)與選址 某公司有6個(gè)建筑工地要開(kāi)工 每個(gè)工地的位置 用平面坐標(biāo)系a b表示 距離單位 千米 及水泥日用量d 噸 由下表給出 目前有兩個(gè)臨時(shí)料場(chǎng)位于A 5 1 B 2 7 日儲(chǔ)量各有20噸 假設(shè)從料場(chǎng)到工地之間均有直線(xiàn)道路相連 1 試制定每天的供應(yīng)計(jì)劃 即從A B兩料場(chǎng)分別向各工地運(yùn)送多少?lài)嵥?使總的噸千米數(shù)最小 2 為了進(jìn)一步減少?lài)嵡讛?shù) 打算舍棄兩個(gè)臨時(shí)料場(chǎng) 改建兩個(gè)新的 日儲(chǔ)量各為20噸 問(wèn)應(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 在各工地用量必須滿(mǎn)足和各料場(chǎng)運(yùn)送量不超過(guò)日儲(chǔ)量的條件下 使總的噸千米數(shù)最小 這是線(xiàn)性規(guī)劃問(wèn)題 線(xiàn)性規(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編寫(xiě)程序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ù)最小 這是非線(xiàn)性規(guī)劃問(wèn)題 非線(xiàn)性規(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 先編寫(xiě)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 取初值為線(xiàn)性規(guī)劃的計(jì)算結(jié)果及臨時(shí)料場(chǎng)的坐標(biāo) x0 35070100406105127 編寫(xiě)主程序gying2 m 3 計(jì)算結(jié)果為 x 3 00005 00004 00007 00001 0000000005 000011 00005 69594 92857 25007 7500 fval 89 8835 注意初始點(diǎn)的選取 某廠向用戶(hù)提供發(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) 第一季度開(kāi)始時(shí)無(wú)存貨 設(shè)a 50 b 0 2 c 4 問(wèn)工廠應(yīng)如何安排生產(chǎn)計(jì)