位置: IT常识 - 正文

最优化理论笔记及期末复习(《数值最优化》——高立)(最优化理论pdf)

编辑:rootadmin
最优化理论笔记及期末复习(《数值最优化》——高立) 目录一、预备知识二、无约束最优化方法的基本结构三、凸集和凸函数四、负梯度方法和Newton型方法五、共轭梯度法六、约束最优化问题的最优性理论七、罚函数方法八、期末复习8.1 知识点复习8.2 习题复习8.3 大实验代码8.3.1实验内容8.3.2实验目的8.3.3算法描述8.3.4程序中的参数设置、终止准则、关键技术(语句)等说明8.3.5实验代码8.3.5.1 目标函数8.3.5.2 计算梯度8.3.5.3 Armijo准则更新步长8.3.5.4最速下降法8.3.5.5 BFGS法8.3.5.6 FR共轭梯度法8.3.5.7 主程序九、总结一、预备知识

推荐整理分享最优化理论笔记及期末复习(《数值最优化》——高立)(最优化理论pdf),希望有所帮助,仅作参考,欢迎阅读内容。

文章相关热门搜索词:最优化理论笔记PDFCSDN,最优化理论课后答案,最优化理论课后答案,最优化理论pdf,最优化理论与方法视频,最优化理论笔记怎么写,最优化理论pdf,最优化理论pdf,内容如对您有帮助,希望把文章链接给更多的朋友!

二、无约束最优化方法的基本结构

三、凸集和凸函数

四、负梯度方法和Newton型方法

五、共轭梯度法

六、约束最优化问题的最优性理论最优化理论笔记及期末复习(《数值最优化》——高立)(最优化理论pdf)

七、罚函数方法

八、期末复习8.1 知识点复习

8.2 习题复习

8.3 大实验代码8.3.1实验内容

利用Matlab编程,实现采用简单Armijo非精确线搜索求步长的三种方法:负梯度法、BFGS法及FR共轭梯度法,并求解如下无约束优化问题: minf(x)=10(x13−x2)2+(x1−1)2min f(x) =10(x_1^3-x_2)^2+(x_1-1)^2minf(x)=10(x13​−x2​)2+(x1​−1)2 通过实验过程进一步理解三种方法的原理和步骤,并对实验结果进行分析比较。

8.3.2实验目的

掌握无约束最优化算法的基本架构,并能熟练使用Matlab软件实现一些基本实用的算法并进行数值试验分析。

8.3.3算法描述

8.3.4程序中的参数设置、终止准则、关键技术(语句)等说明

8.3.5实验代码8.3.5.1 目标函数%%计算函数值function f=func(X)f=10.*(X(1).3-X(2)).2+(X(1)-1).2;end8.3.5.2 计算梯度%计算梯度值function g=grd(X)%计算梯度表达式% syms x1 x2;% f=10*(x13-x2)2+(x1-1)2;% diff(f,x1)% diff(f,x2)% ans = 2*x1 - 60*x12*(- x13 + x2) - 2% ans = - 20*x13 + 20*x2g=[2*X(1) - 60*X(1).2*(- X(1).3 + X(2)) - 2;- 20*X(1).3 + 20*X(2)];end8.3.5.3 Armijo准则更新步长function x=armijo(func,xk,dk,gk)m=0;max_m=1000;rho=0.001;alpha=1;belta=0.618;gd=gk'*dk;fk=feval(func,xk);%初始化条件while m<max_m x=xk+alpha*dk;%试探点 f=feval(func,x);%试探点的函数值 if f<=fk+alpha*rho*gd%终止条件 break; end alpha=alpha*belta;%修改alpha的值 m=m+1;end8.3.5.4最速下降法function [x1 fval1 k1]=fd(x0,func,gfunc,eps,kmax)k1 = 0;x1 = x0;%设置初始条件while k1 < kmax g = feval(gfunc,x1);%计算梯度,x改变时更新梯度 if norm(g)<eps%迭代终止条件 break; end d=-g;%更新方向 x1=armijo(func,x1,d,g);%采用Armijo搜索计算当前点x,最终找到近似最优解 k1=k1+1;endfval1=feval(func,x1);%计算目标函数值8.3.5.5 BFGS法function [x2,fval2,k2]=bfgs(x0,func,grd,H0,eps,kmax)k2=0;H=H0;x2=x0;g=feval(grd,x2);%设置初始条件while k2<kmax if norm(g)<eps%终止条件 break; end d=-H*g;%更新方向 x_=x2;%原来的x x2=armijo(func,x2,d,g);%更新后的x g_=g;%原来的g g=feval(grd, x2);%更新后的梯度 s=x2-x_; y=g-g_; if s'*y>0 v=y'*s; H=H+(1+(y'*H*y)/v)*(s*s')/v-(s*y'*H+H*y*s')/v; %采用BFGS方法更新H end k2=k2+1;endfval2=feval(func,x2);%计算目标函数值8.3.5.6 FR共轭梯度法function [x3,fval3,k3]=FR(x0,func,gfunc,eps,kmax)n=9;k3=0;x3=x0;%设置初始条件while k3<kmax g=feval(gfunc,x3);%更新g m=g'*g;%更新后的g*g if norm(g)<eps%终止条件 break; end if mod(k3,n)==0%n步重新开始策略 d=-g; else belta=m/q;%belta的计算 d=-g+belta*d;%更新d的值 if g'*d>=0 d=-g; end end x3=armijo(func,x3,d,g);%采用Armijo搜索计算当前点,最终找到近似最优解 q=g'*g;%更新前的g*g k3=k3+1;endfval3=feval(func,x3);%计算目标函数值8.3.5.7 主程序clear;clcx0=unifrnd(-5,5,2,1);%产生满足[-5, 5]均匀分布的初始点%x0=[3.4913;-1.0777];%[-5,5]均匀分布产生的初始点...x0=[0.2753;-0.1224];x0=[0.1232;1.1167];x0=[-1.1955;0.6782];x0=[-3.7301;4.1338];x0=[1.3236;-4.0246];...x0=[2.9221;4.3399];x0=[4.5949;1.7874];x0=[1.5574;2.5774];x0=[-4.6429;2.4313];x0=[3.4913;-1.0777]eps=1.e-8;%设置精度1.e-4,1.e-5;1.e-6;1.e-7;1.e-8;kmax=100000;%设置迭代上限H0=eye(2);%H初始为一个2×2的单位矩阵%%采用Armijo搜索的负梯度法程序tic[x1,fval1,k1]=fd(x0,'func','grd',eps,kmax);t1=toc;%%采用Armijo搜索的BFGS法程序tic[x2,fval2,k2]=bfgs(x0,'func','grd',H0,eps,kmax)t2=toc;%%采用Armijo搜索的FR共轭梯度法程序tic[x3,fval3,k3]=FR(x0,'func','grd',eps,kmax);t3=toc;SSE1=sqrt(sum((x1-[1;1]).2,1));%负梯度法下近似解与精确解的2范数下的误差SSE2=sqrt(sum((x2-[1;1]).2,1));%BFGS法下近似解与精确解的2范数下的误差SSE3=sqrt(sum((x3-[1;1]).2,1));%FR共轭梯度法下近似解与精确解的2范数下的误差A=[SSE1 fval1 k1 t1;SSE2 fval2 k2 t2;SSE3 fval3 k3 t3]'%分别记录【误差,函数值,迭代次数,运行时间】

九、总结

本篇文章详细的讲解最优化理论的一些常见方法,有了这些基础的最优化知识,方便我们以后深入学习最优化理论以及人工智能方面的知识。

本文链接地址:https://www.jiuchutong.com/zhishi/299929.html 转载请保留说明!

上一篇:vue3 销毁组件方法(vue destroyed销毁组件)

下一篇:【相机标定】相机内参(什么叫相机标定)

  • clashx使用教程(clash 使用)

    clashx使用教程(clash 使用)

  • 小度如何连接网络(小度如何连接网络机顶盒)

    小度如何连接网络(小度如何连接网络机顶盒)

  • 全民k歌怎么加好友(全民k歌怎么加视频到作品)

    全民k歌怎么加好友(全民k歌怎么加视频到作品)

  • 华为p40视频聊天可以开美颜吗(华为P40视频聊天怎么关闭美颜)

    华为p40视频聊天可以开美颜吗(华为P40视频聊天怎么关闭美颜)

  • 华为nova5EMUI9.1如何隐藏微信电话悬浮窗(华为nova5如何)

    华为nova5EMUI9.1如何隐藏微信电话悬浮窗(华为nova5如何)

  • 微信朋友圈怎么艾特人(微信朋友圈怎么批量删除)

    微信朋友圈怎么艾特人(微信朋友圈怎么批量删除)

  • 快手为什么一个关注都没有还显示2个(快手为什么一个人有十多个号怎么做到的)

    快手为什么一个关注都没有还显示2个(快手为什么一个人有十多个号怎么做到的)

  • 抖音橱窗不显示是怎么回事(为什么抖音橱窗不显示)

    抖音橱窗不显示是怎么回事(为什么抖音橱窗不显示)

  • 淘宝上买了东西以后,店铺不存在了,怎么办(淘宝以前买的东西在哪里找)

    淘宝上买了东西以后,店铺不存在了,怎么办(淘宝以前买的东西在哪里找)

  • 抖音直播分成是多少

    抖音直播分成是多少

  • spnal00华为什么型号(华为手机spnal00什么型号)

    spnal00华为什么型号(华为手机spnal00什么型号)

  • 1000毫安的充电宝能充手机几次(1000毫安的充电宝多少钱)

    1000毫安的充电宝能充手机几次(1000毫安的充电宝多少钱)

  • p20无线充电器能不能用(p20无线充电器能快充吗)

    p20无线充电器能不能用(p20无线充电器能快充吗)

  • 微信换绑手机号需要原手机号的验证码吗(微信换绑手机号收不到验证码)

    微信换绑手机号需要原手机号的验证码吗(微信换绑手机号收不到验证码)

  • vivoy66能改存储位置吗(vivoy67更改存储位置)

    vivoy66能改存储位置吗(vivoy67更改存储位置)

  • 天猫扣分多久清零(天猫扣48分一年多了还在待处理违规里面)

    天猫扣分多久清零(天猫扣48分一年多了还在待处理违规里面)

  • 手机卡显示hd怎么回事(手机卡显示hd怎么调成5g)

    手机卡显示hd怎么回事(手机卡显示hd怎么调成5g)

  • 免拼成功啥意思(免拼成功什么意思)

    免拼成功啥意思(免拼成功什么意思)

  • 微信支出明细怎么查(微信支出明细怎么删除)

    微信支出明细怎么查(微信支出明细怎么删除)

  • xr是什么基带(iphone xr是什么基带芯片)

    xr是什么基带(iphone xr是什么基带芯片)

  • ipad如何显示wifi密码(iPad如何显示电池健康)

    ipad如何显示wifi密码(iPad如何显示电池健康)

  • 微信显示自己手机号码(微信显示自己手机登录的,但是自己没有登录过)

    微信显示自己手机号码(微信显示自己手机登录的,但是自己没有登录过)

  • 为什么网页总是弹广告 DNS被运营商劫持?
	
(为什么网页总是跳掉)

    为什么网页总是弹广告 DNS被运营商劫持? (为什么网页总是跳掉)

  • chown命令  改变文件或目录的用户和用户组(chown命令和chmod)

    chown命令 改变文件或目录的用户和用户组(chown命令和chmod)

  • 财务软件增值税
  • 什么是离岛免税如何购买有哪些特殊规定
  • 组成计税价格是销售额吗
  • 食堂辅助帐怎么记帐
  • 三免三减半如何申报企业所得税
  • 应付利息增加会计分录
  • 营业执照备案需要什么资料
  • 学校食堂支出
  • 开发票时怎么操作收款人 复核人 开票人
  • 社保公司部分交多少比例
  • 预付款什么时候确认收入
  • 从外面买回来的菜怎么消毒
  • 银行贷款利息是百分比还是千分比
  • 公司开给个人的发票有什么用
  • 财产租赁合同印花税计税依据及税率
  • 在建工程预转固后又有工程款
  • 资产相关的递延收益
  • 给研发人员福利怎么做账
  • 个体工商户交企业所得
  • 白酒消费税最低计税价格
  • harmonyos2.0缺点
  • 结转法定盈余公积会计分录怎么写
  • 应交税金年终怎么计算
  • 付款后收到发票怎么写摘要
  • 商品房的维修金该怎么交
  • 企业注销时无力缴纳社保
  • 累计计税折旧如何调整
  • 仙客来养殖方法和注意事项视频
  • 城镇土地使用税纳税义务发生时间
  • webssh github
  • 零售价是含税价还是不含税价
  • 暂估原材料的数量怎么算
  • 委托代理出口账务处理
  • 企业所得税的征税对象包括( )
  • 外管证预缴税款怎么做分录
  • 印花税缴纳方式不包括哪些
  • 调整以前年度亏损
  • 法人股东分红要缴纳多少税
  • 企业一直零申报不经营
  • 发票未报销可以退税吗
  • mongodb开启远程连接
  • 福利费要分部门吗
  • 什么是,利息
  • 防伪税控开票系统年费
  • 金融资产经营资产
  • 会计分录记录于什么中
  • 银行贷款是应付账款吗
  • 专项应付款和政府补助的区别
  • 进口贴息资金
  • 研发废料收入会计科目
  • 怎么把应付账款调成应收
  • 期间损益科目如何结转
  • 固定资产的核算包括
  • 记账凭证模板
  • 企业建账的流程是怎样的
  • 总账建账的原则包括
  • win10怎么关闭securboot
  • cygwin下载教程
  • win1020h2版好不好
  • windowsxp桌面没有我的电脑
  • macos 网速慢
  • fsck.mode
  • Linux系统中矢量图ai格式怎么打开?
  • 深入了解linux内核
  • shell脚本tr
  • js组件是什么
  • python 解析
  • jQuery插件库
  • javascript基础题目
  • 2016年最火的歌曲排行榜前十名
  • 使用jquery实现的项目
  • jq 判断元素是否存在
  • 江苏省地方税务局网上办税
  • 安家费购房补贴区别
  • 北京国家税务局待遇
  • 浙江 国家地理
  • 关于返回
  • 河北电子税务局社保缴费流程
  • 印花税票在哪里可以查到
  • 免责声明:网站部分图片文字素材来源于网络,如有侵权,请及时告知,我们会第一时间删除,谢谢! 邮箱:opceo@qq.com

    鄂ICP备2023003026号

    网站地图: 企业信息 工商信息 财税知识 网络常识 编程技术

    友情链接: 武汉网站建设