【正文】 c1=sqrt((5a0).^2+(1b0).^2)。 c2=sqrt((2a0).^2+(7b0).^2)。 c=[c1,c2]。 A=[ones(1,6),zeros(1,6)。zeros(1,6),ones(1,6)]。 b=[20。20]。 Aeq=[eye(6),eye(6)]。 beq=[3 5 4 7 6 11]39。 L=zeros(1,12)。 [Z,val]=linprog(c, A, b, Aeq,beq ,L) ?單擊此處編輯母版文本樣式 ?第二級(jí) ?第三級(jí) ?第四級(jí) ?第五級(jí) 43 單擊此處編輯母版文本樣式第二級(jí)第三級(jí)第四級(jí)第五級(jí) ?單擊此處編輯母版副標(biāo)題樣式 單擊此處編輯母版文本樣式第二級(jí)第三級(jí)第四級(jí)第五級(jí)數(shù)學(xué)實(shí)驗(yàn)之 －－非線性規(guī)劃 實(shí)驗(yàn)?zāi)康? 引 例 基本概念 算法概述 軟件求解 范 例 布置實(shí)驗(yàn) 結(jié) 束 課堂延伸 最優(yōu)目標(biāo)值 f = （噸千米） 料場(chǎng) A, B運(yùn)往各工地的水泥的日運(yùn)量分別為 i 1 2 3 4 5 6 zi1(料場(chǎng) A) 3 5 0 7 0 1 zi2(料場(chǎng) B) 0 0 4 0 6 10 供應(yīng)與選址 范 例 問(wèn)題 1的求解結(jié)果 ?單擊此處編輯母版文本樣式 ?第二級(jí) ?第三級(jí) ?第四級(jí) ?第五級(jí) 44 單擊此處編輯母版文本樣式第二級(jí)第三級(jí)第四級(jí)第五級(jí) ?單擊此處編輯母版副標(biāo)題樣式 單擊此處編輯母版文本樣式第二級(jí)第三級(jí)第四級(jí)第五級(jí)數(shù)學(xué)實(shí)驗(yàn)之 －－非線性規(guī)劃 實(shí)驗(yàn)?zāi)康? 引 例 基本概念 算法概述 軟件求解 范 例 布置實(shí)驗(yàn) 結(jié) 束 課堂延伸 問(wèn)題 2的求解 要為新建料場(chǎng)選址 , 料場(chǎng)位置 (xj, yj)為未知時(shí) , 決策變量為 zij， xj, yj，模型為非線性規(guī)劃模型。 2622112161m in ( ) ( ). . , 1 , , 6, 1 , 20ij j i j ijiij ijij jiijf z x a y bs t z d iz r jz????? ? ? ??????????（工地日用量） （料場(chǎng)日儲(chǔ)量） 供應(yīng)與選址 范 例 ?單擊此處編輯母版文本樣式 ?第二級(jí) ?第三級(jí) ?第四級(jí) ?第五級(jí) 45 單擊此處編輯母版文本樣式第二級(jí)第三級(jí)第四級(jí)第五級(jí) ?單擊此處編輯母版副標(biāo)題樣式 單擊此處編輯母版文本樣式第二級(jí)第三級(jí)第四級(jí)第五級(jí)數(shù)學(xué)實(shí)驗(yàn)之 －－非線性規(guī)劃 實(shí)驗(yàn)?zāi)康? 引 例 基本概念 算法概述 軟件求解 范 例 布置實(shí)驗(yàn) 結(jié) 束 課堂延伸 目標(biāo)函數(shù)的函數(shù) M 文件 : function f=liaocmb(x) a0=[ 3 ]。 b0=[ 5 ]。 c1=sqrt((x(13)a0).^2+(x(14)b0).^2)。 c2=sqrt((x(15)a0).^2+(x(16)b0).^2)。 c=[c1,c2]。 f=c*x(1:12,1)。 供應(yīng)與選址 范 例 問(wèn)題 2的求解 ?單擊此處編輯母版文本樣式 ?第二級(jí) ?第三級(jí) ?第四級(jí) ?第五級(jí) 46 單擊此處編輯母版文本樣式第二級(jí)第三級(jí)第四級(jí)第五級(jí) ?單擊此處編輯母版副標(biāo)題樣式 單擊此處編輯母版文本樣式第二級(jí)第三級(jí)第四級(jí)第五級(jí)數(shù)學(xué)實(shí)驗(yàn)之 －－非線性規(guī)劃 實(shí)驗(yàn)?zāi)康? 引 例 基本概念 算法概述 軟件求解 范 例 布置實(shí)驗(yàn) 結(jié) 束 課堂延伸 function [c,ceq]=liaocys(x) A=[ones(1,6),zeros(1,6)。zeros(1,6),ones(1,6)]。 b=[20。20]。 Aeq=[eye(6), eye(6)]。 beq=[3 5 4 7 6 11]39。 c=A*x(1:12,1)b。 ceq=Aeq*x(1:12,1)beq。 約束條件的函數(shù) M 文件 : 供應(yīng)與選址 范 例 問(wèn)題 2的求解 ?單擊此處編輯母版文本樣式 ?第二級(jí) ?第三級(jí) ?第四級(jí) ?第五級(jí) 47 單擊此處編輯母版文本樣式第二級(jí)第三級(jí)第四級(jí)第五級(jí) ?單擊此處編輯母版副標(biāo)題樣式 單擊此處編輯母版文本樣式第二級(jí)第三級(jí)第四級(jí)第五級(jí)數(shù)學(xué)實(shí)驗(yàn)之 －－非線性規(guī)劃 實(shí)驗(yàn)?zāi)康? 引 例 基本概念 算法概述 軟件求解 范 例 布置實(shí)驗(yàn) 結(jié) 束 課堂延伸 clear L=zeros(16,1)。 x0=[zeros(1,12),5,1,2,7]39。 options=optimset(39。largescale39。,39。off39。,39。display39。,39。iter39。,39。MaxFunEval39。,2021)。 [x,val,]=fmincon(39。liaocmb39。,x0, [ ], [ ], [ ],[ ] ,L,[ ],39。liaocys39。,options) 解非線性規(guī)劃的主程序 供應(yīng)與選址 范 例 問(wèn)題 2的求解 ?單擊此處編輯母版文本樣式 ?第二級(jí) ?第三級(jí) ?第四級(jí) ?第五級(jí) 48 單擊此處編輯母版文本樣式第二級(jí)第三級(jí)第四級(jí)第五級(jí) ?單擊此處編輯母版副標(biāo)題樣式 單擊此處編輯母版文本樣式第二級(jí)第三級(jí)第四級(jí)第五級(jí)數(shù)學(xué)實(shí)驗(yàn)之 －－非線性規(guī)劃 實(shí)驗(yàn)?zāi)康? 引 例 基本概念 算法概述 軟件求解 范 例 布置實(shí)驗(yàn) 結(jié) 束 課堂延伸 i 1 2 3 4 5 6 zi1(料場(chǎng) A) 3 0 4 7 6 0 zi2(料場(chǎng) B) 0 5 0 0 0 11 計(jì)算結(jié)果： 最優(yōu)目標(biāo)值 f = （噸千米） 新料場(chǎng)位置的改變 ， 目標(biāo)值比原來(lái)減少了 。 新料場(chǎng) A, B的坐標(biāo)為 （ ， ） 和（ ， ） 。 新料場(chǎng) A, B運(yùn)往各工地的水泥的日運(yùn)量分別為 問(wèn)題 2的求解 ?單擊此處編輯母版文本樣式 ?第二級(jí) ?第三級(jí) ?第四級(jí) ?第五級(jí) 49 單擊此處編輯母版文本樣式第二級(jí)第三級(jí)第四級(jí)第五級(jí) ?單擊此處編輯母版副標(biāo)題樣式 單擊此處編輯母版文本樣式第二級(jí)第三級(jí)第四級(jí)第五級(jí)數(shù)學(xué)實(shí)驗(yàn)之 －－非線性規(guī)劃 實(shí)驗(yàn)?zāi)康? 引 例 基本概念 算法概述 軟件求解 范 例 布置實(shí)驗(yàn) 結(jié) 束 課堂延伸 優(yōu)化結(jié)果是新料場(chǎng)應(yīng)建在用量最大的工地旁邊，你預(yù)先估計(jì)到這個(gè)結(jié)果了嗎？ 供應(yīng)與選址 范 例 ?單擊此處編輯母版文本樣式 ?第二級(jí) ?第三級(jí) ?第四級(jí) ?第五級(jí) 50 單擊此處編輯母版文本樣式第二級(jí)第三級(jí)第四級(jí)第五級(jí) ?單擊此處編輯母版副標(biāo)題樣式 單擊此處編輯母版文本樣式第二級(jí)第三級(jí)第四級(jí)第五級(jí)數(shù)學(xué)實(shí)驗(yàn)之 －－非線性規(guī)劃 實(shí)驗(yàn)?zāi)康? 引 例 基本概念 算法概述 軟件求解 范 例 布置實(shí)驗(yàn) 結(jié) 束 課堂延伸 求解非線性規(guī)劃 布置實(shí)驗(yàn) 10 10111 2 3 6 104 5 6 73 7 8 9 106m in ( ) ( l n( / ) ). . 2 2 2212110 , 1 , .. ., 10j j j ijijf x x c x xs t x x x x xx x x xx x x x xxj?????? ? ? ? ?? ? ? ?? ? ? ? ?????其中 cj, j=1,… ,10是常數(shù) , 其值分別為 : {, , , , , … , , , , , } ?單擊此處編輯母版文本樣式 ?第二級(jí) ?第三級(jí) ?第四級(jí) ?第五級(jí) 51 單擊此處編輯母版文本樣式第二級(jí)第三級(jí)第四級(jí)第五級(jí) ?單擊此處編輯母版副標(biāo)題樣式 單擊此處編輯母版文本樣式第二級(jí)第三級(jí)第四級(jí)第五級(jí)數(shù)學(xué)實(shí)驗(yàn)之 －－非線性規(guī)劃 實(shí)驗(yàn)?zāi)康? 引 例 基本概念 算法概述 軟件求解 范 例 布置實(shí)驗(yàn) 結(jié) 束 課堂延伸 求解非線性規(guī)劃 5,4,3,2,1,0105。10..)e x p(m i n3231543251254321?????????????????ixixxxxxxxxtsxxxxxfiiii布置實(shí)驗(yàn) ?單擊此處編輯母版文本樣式 ?第二級(jí) ?第三級(jí) ?第四級(jí) ?第五級(jí) 52 單擊此處編輯母版文本樣式第二級(jí)第三級(jí)第四級(jí)第五級(jí) ?單擊此處編輯母版副標(biāo)題樣式 單擊此處編輯母版文本樣式第二級(jí)第三級(jí)第四級(jí)第五級(jí)數(shù)學(xué)實(shí)驗(yàn)之 －－非線性規(guī)劃 實(shí)驗(yàn)?zāi)康? 引 例 基本概念 算法概述 軟件求解 范 例 布置實(shí)驗(yàn) 結(jié) 束 課堂延伸 桃李花園服務(wù)中心選址 P158, 操練一。 設(shè) (ai, bi)(i=1~20)為第 i 棟住宅樓的坐標(biāo)； a = [ ]。 b =[ ]。 布置實(shí)驗(yàn)