【正文】
^4。options(6)=0。x=fminu(39。,x0,options), 輸出結(jié)果 x = 執(zhí)行程序 BFGS立方插值 x0=[3,1,0,1]。options(7)=1。f439。options(6)=1。x=fminu(39。,x0,options), x = 輸出結(jié)果 執(zhí)行程序 DFP立方插值 x0=[3,1,0,1]。options(7)=1。f439。或 leastsq(‘fun’,x0)。 c=[,]。 執(zhí)行程序 x0=[10,]。ct39。 f 和約束條件 g 。 f 和約束條件 g 的梯度。 例題 1 f(x1,x2)=100(x2- x12)2+(1- x1)2 x12+x12≤,x1+x2 ≥ 0 編寫程序 function [f,g]=f2(x) f=100*(x(2)x(1)^2)^2+(1x(1))^2。g(2)=x(1)x(2)。opt(1)=1。f239。opt(1)=1。f239。grad39。 dg=[2*x(1),1。 輸出結(jié)果 梯度 最優(yōu)解 最優(yōu)值 迭代次數(shù) 不給出 ( ,) 69 給出 (,) 36 輸出結(jié)果 x = f = 編寫程序 x0=[1,1]。[x,opt]=constr(39。,x0,opt)。 g(1)=+x(1)*x(2)x(1)x(2)。 例題 3 求解下列條件極值問(wèn)題 1221 2 1 2 21 2 1 21212m in ( ) ( 4 2 4 2 1 ). . 1. 5 0 10 0 0xf x e x x x x xsu b to x x x xxxxx? ? ? ? ?? ? ? ?????編寫程序 x0=[1,1]。vlb=[0。vub=[]。f539。 x,f=opt(8),n=opt(10) 輸出結(jié)果 x = 0 f = n = 10 用戶提供導(dǎo)數(shù)時(shí) ,則問(wèn)題會(huì)更有效更準(zhǔn)確 . 例題 4 求解下列極值 10..)12424()(mi n21212122122211????????????????????????????????????????????xxxxxxtos u bxxxxxexf x編寫程序 ?x0=[1,1]。 [x,opt]=constr(39。,x0,opt,[],[],39。)。 dg=[x(2)1,x(2)。 x = f = n = 11 輸出結(jié)果 等式約束問(wèn)題 例 5: 求下列極值 100..)12424()(m i n2121212122122211????????????????????????????????????xxxxxxxxtos u bxxxxxexfx編寫程序 function [f,g]=f6(x) f=exp(x(1))*(4*x(1)^2+2*x(2)^2+4*x(1)*x(2)+2*x(2)+1)。g(2)=+x(1)*x(2)x(1)x(2)。 執(zhí)行程序 x0=[1,1]。opt(13)=1。f639。 x,f=opt(8),n=opt(10) 輸出結(jié)果 x = f = n = 13