logo

您所在位置网站首页 > 海量文档  > 专业论文 > 论文指导设计

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

本文档一共被下载: ,您可全文免费在线阅读后下载本文档。

  • 支付并下载
  • 收藏该文档
  • 百度一下本文档
  • 修改文档简介
全屏预览

下载提示

1.本站不保证该用户上传的文档完整性,不预览、不比对内容而直接下载产生的反悔问题本站不予受理。
2.该文档所得收入(下载+内容+预览三)归上传者、原创者。
3.登录后可充值,立即自动返金币,充值渠道很便利
第十二章 非线性极值 非线性规划与线性规划的区别是:其目标函数与约束条件至少有一处是非线性的。非线性规划问题可分为无约束问题和有约束问题。 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

发表评论

请自觉遵守互联网相关的政策法规,严禁发布色情、暴力、反动的言论。
用户名: 验证码: 点击我更换图片

“原创力文档”前称为“文档投稿赚钱网”,本站为“文档C2C交易模式”,即用户上传的文档直接卖给(下载)用户,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有【成交的100%(原创)】。原创力文档是网络服务平台方,若您的权利被侵害,侵权客服QQ:3005833200 电话:19940600175 欢迎举报,上传者QQ群:784321556