【正文】 有一個(gè)目標(biāo)函數(shù)，是單目標(biāo)最優(yōu)化方法。但是，在許多實(shí)際工程問(wèn)題中，往往希望多個(gè)指標(biāo)都達(dá)到最優(yōu)值，所以它有多個(gè)目標(biāo)函數(shù)。這種問(wèn)題稱(chēng)為多目標(biāo)最優(yōu)化問(wèn)題。 多目標(biāo)最優(yōu)化問(wèn)題的數(shù)學(xué)模型為 )(m i n xFnRx ? 1,2,10)( mixGi??? mmmixG i ,2,10)( 11 ????? uixxx ?? 其中)( xF為目標(biāo)函數(shù)向量。 此優(yōu)化問(wèn)題在 Ma t l a b 中主要由函數(shù) f g oal at t a i n 來(lái)實(shí)現(xiàn)。此問(wèn)題在控制系統(tǒng)中有廣泛的應(yīng)用。 函數(shù) fgoalattain min v . fi(X)wiv ≤goali i=1,2,…,t AX≤b （線性不等式約束） AeqX=beq （線性等式約束） C(X)≤0 （非線性不等式約束條件） Ceq(X)=0 （非線性等式約束） Lb ≤X ≤Ub （邊界約束條件） 一、多目標(biāo)優(yōu)化問(wèn)題數(shù)學(xué)模型 標(biāo)量變量 各分目標(biāo)函數(shù) 分目標(biāo)函數(shù)的權(quán)重 各分目標(biāo)函數(shù)的目標(biāo)值 二、優(yōu)化函數(shù)使用格式 [x,fval,exitflag,output, grad,hessian]= fgoalattain(fun,x0,goal,w,A,b,Aeq,beq,Lb,Ub,’ Nlc’ ,options,P1,P2…) 設(shè)置優(yōu)化選項(xiàng)參數(shù) 各分目標(biāo)權(quán)重 各分目標(biāo)期望值 目標(biāo)函數(shù)在最優(yōu)解的海色矩陣 返回目標(biāo)函數(shù)在最優(yōu)解的梯度 優(yōu)化算法信息的一個(gè)數(shù)據(jù)結(jié)構(gòu) 返回算法的終止標(biāo)志 返回目標(biāo)函數(shù)的最優(yōu)值 返回目標(biāo)函數(shù)的最優(yōu)解 附加參數(shù) 非線性約束條件的函數(shù)名 設(shè)計(jì)變量的下界和上界 線性等式約束的常數(shù)向量 線性等式約束的系數(shù)矩陣 線性不等式約束的常數(shù)向量 線性不等式約束的系數(shù)矩陣 無(wú)定義時(shí)以空矩陣 符號(hào)“ [ ]”代替 函數(shù) fgoalattain 初始點(diǎn) 目標(biāo)函數(shù)文件名 ? 格式： x = fgoalattain(fun,x0,goal,weight) ? x = fgoalattain(fun,x0,goal,weight,A,b) ? x = fgoalattain(fun,x0,goal,weight,A,b,Aeq,beq) ? x = fgoalattain(fun,x0,goal,weight,A,b,Aeq,beq,lb,ub) ? x = fgoalattain(fun,x0,goal,weight,A,b,Aeq,beq,lb,ub,nonlcon) ? x = fgoalattain(fun,x0,goal,weight,A,b,Aeq,beq, lb,ub,nonlcon,options) ? x = fgoalattain(fun,x0,goal,weight,A,b,Aeq,beq, lb,ub,nonlcon,options,P1,P2,...) x = fgoalattain(fun,x0,goal,weight,A,b,Aeq,beq,lb,ub,nonlcon) ? 說(shuō)明： ? fgoalattain 求解多目標(biāo)達(dá)到問(wèn)題。 ? x = fgoalattain(fun,x0,goal,weight) 試圖通過(guò)變化 x來(lái)使目標(biāo)函數(shù) fun達(dá)到 goal指定的目標(biāo)。初值為 x0， weight參數(shù)指定權(quán)重。 ? x = fgoalattain(fun,x0,goal,weight,A,b) 求解目標(biāo)達(dá)到問(wèn)題，約束條件為線性不等式 A*x = b。 ? x = fgoalattain(fun,x0,goal,weight,A,b,Aeq,beq) 求解目標(biāo)達(dá)到問(wèn)題，除提供上面的線性不等式外，還提供線性等式 Aeq*x = beq。當(dāng)沒(méi)有不等式存在時(shí)，設(shè)置 A=[ ]和 b=[ ]。 ? x = fgoalattain(fun,x0,goal,weight,A,b,Aeq,beq,lb,ub) 為設(shè)計(jì)變量 x定義下界 lb和上界 ub集合，這樣始終有 lb = x = ub。 ? x = fgoalattain(fun,x0,goal,weight,A,b,Aeq,beq,lb,ub,nonlcon) 將目標(biāo)達(dá)到問(wèn)題歸結(jié)為 nonlcon參數(shù)定義的非線性不等式 c (x)或非線性等式 ceq (x)。 fgoalattain優(yōu)化的約束條件為 c(x) = 0 和 ceq(x) = 0。若不存在邊界，設(shè)置 lb=[ ] 和（或） ub=[ ]。 ? x = fgoalattain(fun,x0,goal,weight,A,b,Aeq,beq,lb,ub,nonlcon,options) 用 options中設(shè)置的優(yōu)化參數(shù)進(jìn)行最小化。 ? x = fgoalattain(fun,x0,goal,weight,A,b,Aeq,beq,lb,ub,nonlcon,options,P1,P2,...) 將問(wèn)題參數(shù) P1， P2等直接傳遞給函數(shù) fun和 nonlcon。若不需要參數(shù) A， b， Aeq，beq， lb， ub， nonlcon和 options，將它們?cè)O(shè)置為空矩陣。 例 10 15 某化工廠擬生產(chǎn)兩種新產(chǎn)品 A 和 B ，其生產(chǎn)設(shè)備費(fèi)用分別為： A ， 2 萬(wàn)元 / 噸； B ， 5 萬(wàn)元 / 噸。這兩種產(chǎn)品均將造成 環(huán)境污染，設(shè)由公害所造成的損失可折算為： A ， 4 萬(wàn)元 / 噸； B ， 1 萬(wàn)元 / 噸。由于條件限制，工廠生產(chǎn)產(chǎn)品 A 和 B 的最大生產(chǎn)能力各為每月 5 噸和 6噸，而市場(chǎng)需要這兩種產(chǎn)品的總量每月不少于 7 噸。試問(wèn)工廠如何安排生產(chǎn)計(jì)劃，在滿(mǎn)足市場(chǎng)需要的前提下，使設(shè)備投資和公害損失均達(dá)到最小。該工廠決策認(rèn)為，這兩個(gè)目標(biāo)中環(huán)境污染應(yīng)優(yōu)先考慮，設(shè)備投資的目標(biāo)值為 20 萬(wàn)元，公害損失的目標(biāo)為 12 萬(wàn)元。 解： 設(shè)工廠每月生產(chǎn)產(chǎn)品 A 為 x1噸， B 為 x2噸，設(shè)備投資費(fèi)為 f1( x) ，公害損失費(fèi)為 f2( x ) ，則這個(gè)問(wèn)題可表達(dá)為多目標(biāo)優(yōu)化問(wèn)題： m i n f1 ( x) = 2x1+5x2 m i n f2 ( x) = 4x1+x2 x1≤ 5 x2≤ 6 x1+x2≥ 7 x1， x2≥ 0 ? 需要編寫(xiě)目標(biāo)函數(shù)的 M文件 ，返回目標(biāo)計(jì)算值： ? function f=Ex1015(x) ? f(1)=2*x(1)+5*x(2)。 ? f(2)=4*x(1)+x(2)。 ? 給定目標(biāo)，權(quán)重按目標(biāo)比例確定，給出初值，在 Matlab中實(shí)現(xiàn)為： ? goal=[20 12]。 ? weight=[20 12]。 ? x0=[2 5]。 ? A=[1 0。0 1。 1 1]。 ? b=[5 6 7]。 ? lb=zeros(2,1)。 ? [x,fval,attainfactor,exitflag]=fgoalattain(Ex1015,x0,goal,weight,A,b,[],[],lb,[]) ? x = ? ? fval = ? ? attainfactor = ? ? exitflag = ? 1 ? 故工廠每月生產(chǎn)產(chǎn)品 A為 ， B為 。設(shè)備投資費(fèi)和公害損失費(fèi)的目標(biāo)值分別為 和 。 ? 例 1016 某工廠因生產(chǎn)需要欲采購(gòu)一種原材料，市場(chǎng)上的這種原料有兩個(gè)等級(jí)，甲級(jí)單價(jià) 2元 /kg，乙級(jí)單價(jià) 1元 /kg。要求所花總費(fèi)用不超過(guò) 200元，購(gòu)得原料總量不少于 100kg，其中甲級(jí)原料不少于 50kg，問(wèn)如何確定最好的采購(gòu)方案。 ? 解： 設(shè) x1和 x2分別為采購(gòu)甲級(jí)和乙級(jí)原料的數(shù)量（ kg），要求采購(gòu)總費(fèi)用盡量少，采購(gòu)總量盡量多，采購(gòu)甲級(jí)原料盡量多。由題意建立下面的數(shù)學(xué)模型： ? min z1 = 2x1+x2 ? max z2 = x1+x2 ? max z3 = x1 ? 2x1+x2≤200 ? x1+x2≥100 ? x1≥50 ? x1， x2≥0 ? 需要編寫(xiě)目標(biāo)函數(shù)的 M文件 ，返回目標(biāo)計(jì)算值： ? function f=Ex1016(x) ? f(1)=2*x(1)+x(2)。 ? f(2)=x(1)x(2)。 ? f(3)=x(1)。 ? 給定目標(biāo)，權(quán)重按目標(biāo)比例確定，給出初值，在 Matlab中實(shí)現(xiàn)： ? goal=[200 100 50]。 ? weight=[200 100 50]。 ? x0=[50 50]。 ? A=[2 1。1 1。1 0]。 ? b=[200 100 50]。 ? lb=zeros(2,1)。 ? [x,fval,attainfactor,exitflag]=fgoalattain(Ex1016,x0,goal,weight,A,b,[],[],lb,[]) ? x = ? 50 50 ? fval = ? 150 100 50 ? 所以，最好的采購(gòu)方案是采購(gòu)甲級(jí)原料和乙級(jí)原料各 50kg。此時(shí)采購(gòu)總費(fèi)用為 150元，總重量為 100kg，甲級(jí)原料總重量為 50kg。 The end! 謝謝觀看 /歡迎下載 BY FAITH I MEAN A VISION OF GOOD ONE CHERISHES AND THE ENTHUSIASM THAT PUSHES ONE TO SEEK ITS FULFILLMENT REGARDLESS OF OBSTACLES. BY FAITH I BY FAITH