1、非線性優(yōu)化問題函數(shù) fmincons. t AX-bAq X = beg G(X)MO Cq(X0(線性不等式約束)(統(tǒng)性錚式約束)非稅性不等式約束)(非線性等式約束)，VLB三X VUB 有算約束使用格式：x,fval,exitflag,output,lambda,grad,hessian = fmincon(fun,x0,A,b,Aeq,beq,lb,ub,nonlcon,options)x0為開始運算的初始值x為最優(yōu)解時x的值fval為z的最優(yōu)值exitflag判斷算法停止的原因grad梯度這個函數(shù)的基本形式為x = fmincon(fufiaxOaAab,Aeq.beq.lbbonlc

2、on； options)其中firn為你要求最小值的函數(shù)，可以單寫一個文件設(shè)置函數(shù)，如以上給的例子中。L如果版中有N個變量，如X”,或者是什么的，自己排個順序，在Wi中統(tǒng)一都是用x(l)ax(2) ._x(n)表示的。表示初始的猜測值，大小要與變量數(shù)目相同3. Ab為線性不等約束，A*x= b, A應(yīng)為n*n階矩陣，學(xué)過線性代數(shù)應(yīng)不難寫出斯口卜Aeqbcq為線性相等約束，Aeq*x = bcq。Aeqbeq同上可求Ibub為變量的上下邊界，正負(fù)無窮用-In琲口hi臻示,Ibub應(yīng)為X階數(shù)組nonIcon為非線性約束，可分為兩部分，非技性不等約束6非線性相等約束，ceq可按下面的例子設(shè)置func

3、tion c:ce = nonlconi)c = -x( 5xQ)24;oe = :% no nonlinear equality constraints7,最后是options,可以用OFTMEE底數(shù)設(shè)置，見上例具體可見OPTIMSE底數(shù)的幫助文件。題求下列函數(shù)的域小值例1.二瀆優(yōu)化*min f - T - 2r; 十；& + ；工s t 2xi+3k20m文件function y=funn(x)y=-x(1)-2*x(2)+0.5*(x(1).八2+x(2).八2)主程序A=2 3;1 4;b=6;5;vlb=0;0;x0=0;0;x,fval,exitflag,output,lambda

4、,grad,hession=fmincon(funn,x0,A,b,，口，vlb)結(jié)果如下Active inequalities (to within options.TolCon = 1e-006):lower upper ineqlin ineqnonlin2x =0.76471.0588fval =-2.0294exitflag =output =iterations: 1funcCount: 6Issteplength: 1stepsize: 1.3061algorithm: medium-scale: SQP, Quasi-Newton, line-search firstorder

5、opt: 1.7969e-008constrviolation: 0message: 1x788 charlambda =lower: 2x1 doubleupper: 2x1 doubleeqlin: 0 x1 doubleeqnonlin: 0 x1 doubleineqlin: 2x1 doubleineqnonlin: 0 x1 doublegrad =-0.2353-0.9412hession =1001假如約束條件是非線性約束條件，即2a;6，應(yīng)該怎樣解？首先增加一個M文件 function y=nonline(x) y=2*x(1)+3*x(2).八2-6主程序A=2 3;1 4

6、;b=6;5;vlb=0;0;x0=0;0;x,fval,exitflag,output,lambda,grad,hession=fmincon(funn,x0,A,b,vlb,nonline)x =0.76471.0588fval =-2.0294exitflag =1output =iterations: 1funcCount: 6lssteplength: 1stepsize: 1.3061algorithm: medium-scale: SQP, Quasi-Newton, line-search firstorderopt: 1.7969e-008constrviolation: 0message: 1x788 charlambda =lower: 2x1 double upper: 2x1 double eqlin: 0 x1 double