MATLAB课件第十二章非线性极值.doc

第十二章 非线性极值 非线性规划与线性规划的区别是：其目标函数与约束条件至少有一处是非线性的。非线性规划问题可分为无约束问题和有约束问题。 1．无约束非线性规划问题： 1.1有界单变量优化 函数fminbnd 调用格式：[x,val]=fminbnd(f,x1,x2),其中f是用来求极值的函数，可以是函数名，也可以是函数表达式，意思是求函数f在区间[x1,x2]上的极小值（不是最小值）。 【例1】 求函数的极值 为了能更方便的找出极值点，先用plot函数画出该函数的曲线图，输入如下命令： x=-2:0.1:2;f=(x.^2-1).^3+1;plot(x,f) 图1 从图中可以看出函数有极小值.输入如下命令： f1='(x.^2-1).^3+1'; [x,val]=fminbnd(f1,-2,2) 回车后可得：x = 4.4409e-016；val =0，即可知极小值为f(0)=0 【例2】 求解f(x)=x^2-2x-1的极值 解答：先作图： x=-1:0.1:3;f=x.^2-2*x-1;plot(x,f) 图2 从图中可以看出函数有极小值.输入如下命令： f1=' x.^2-2*x-1'; [x,val]=fminbnd(f1,-1,3) 回车后可得：x =1；val =-2，即可知极小值为f(1)=-2 2．求解多元无约束最优化问题 函数fminunc 调用格式： x=fminunc(fun,x0) x=fminunc(fun,x0,options) x=fminunc(fun,x0,options,P1,P2) [x,fval]=fminunc(…) [x,fval, exitflag]=fminunc(…) [x,fval, exitflag,output]=fminunc(…) [x,fval, exitflag,output,grad]=fminunc(…) [x,fval, exitflag,output,grad,hessian]=fminunc(…) 说明：fun为需最小化的目标函数，x0为给定的搜索的初始点。options指定优化参数。 返回的x为最优解向量；fval为x处的目标函数值；exitflag描述函数的输出条件；output返回优化信息；grad返回目标函数在x处的梯度。Hessian返回在x处目标函数的Hessian矩阵信息。 【例3】 ： 求 通过绘图确定一个初始点： [x,y]=meshgrid(-10:0.5:10);z= 8*x-4*y +x.^2+3*y.^2; surf(x,y,z) 图3 选初始点：x0=(0,0) x0=[0,0]; 编辑ff1.m文件 function f=ff1(x) f=8*x(1)-4*x(2) +x(1)^2+3*x(2)^2; 程序调用: x0=[0,0]; [x,fval,exitflag]=fminunc(@ff1,x0) ? 结果：x = -4.0000 0.6667 fval = -17.3333 exitflag = 1 【例4】： 解：通过绘图确定一个初始点： [x,y]=meshgrid(-10:0.5:10); z= 4*x.^2+5*x.*y +2*y.^2; surf(x,y,z) 图4 选初始点：x0=(0,0) x0=[0,0]; 编辑ff2.m文件 function f=ff2(x) f=4*x(1).^2+5*x(1).*x(2) +2*x(2).^2; 程序调用: x0=[0,0]; [x,fval,exitflag]=fminunc(@ff2,x0) 结果： x = 0 0 fval =0 exitflag = 1 【例5】：计算的极值 解：1)求极小值 编写目标函数: function f=li5_16(x) f=x(1)^3-x(2)^3+3*x(1)^2+3*x(2)^2-9*x(1); 主程序窗口调用： x0=[0 0]; [x,fval, exitflag]=fminunc(@li5_16 ,x0) 计算结果： x = 1.0000 -0.0000 fval = -5 exitflag = 1 即f(1,0)=-5为极小值。 2）求极大值：可以变成求目标函数负的极小值，然后反号即可求出极大值。 编写目标函数: function f=li5_161(x) f=-x(1)^3+x(2)^3-3*x(1)^2-3*x(2)^2+9*x(1); 主程序窗口调用： x0=[-1 1]; [x,fval, exitflag]=fminunc(@li5_161 ,x0),maxf=-fval