分享
分享赚钱 收藏 举报 版权申诉 / 2

类型MATLAB最速下降法(梯度法)程序.doc

  • 上传人:kpmy5893
  • 文档编号:7654049
  • 上传时间:2019-05-23
  • 格式:DOC
  • 页数:2
  • 大小:27.50KB
  • 配套讲稿:

    如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。

    特殊限制:

    部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。

    关 键  词:
    MATLAB最速下降法(梯度法)程序.doc
    资源描述:

    1、 查看文章 matlab最速下降法2010-08-18 17:13function x=fsxsteep(f,e,a,b)% fsxsteep函数 最速下降法% x=fsxsteep(f,e,a,b)为输入函数 f 为函数 e 为允许误差 (a,b)为初始点;% fsx TJPU 2008.6.15x1=a;x2=b;Q=fsxhesse(f,x1,x2);x0=x1 x2;fx1=diff(f,x1); %对 x1求偏导数fx2=diff(f,x2); %对 x2求偏导数g=fx1 fx2; %梯度g1=subs(g); %把符号变量转为数值d=-g1;while (abs(norm(g1)

    2、=e)t=(-d)*d/(-d)*Q*d);t=(-d)*d/(-d)*Q*d); %求搜索方向x0=x0-t*g1; %搜索到的点v=x0;a=1 0*x0;b=0 1*x0;x1=a;x2=b;g1=subs(g);d=-g1;end;x=v;function x=fsxhesse(f,a,b)% fsxhesse函数 求函数的 hesse矩阵;% 本程序仅是简单的求二次函数的 hesse矩阵!;% x=fsxhesse(f)为输入函数 f 为二次函数 x1,x2 为自变量;% fsx TJPU 2008.6.15x1=a;x2=b;fx=diff(f,x1); %求 f对 x1偏导数fy

    3、=diff(f,x2); %求 f对 x2偏导数fxx=diff(fx,x1); %求二阶偏导数 对 x1再对 x1fxy=diff(fx,x2); %求二阶偏导数 对 x1再对 x2fyx=diff(fy,x1); %求二阶偏导数 对 x2再对 x1fyy=diff(fy,x2); %求二阶偏导数 对 x2再对 x2fxx=subs(fxx); %将符号变量转化为数值fxy=subs(fxy);fyx=subs(fyx);fyy=subs(fyy);x=fxx,fxy;fyx,fyy; %求 hesse矩阵syms x1 x2;X=x1,x2;fx=X(1)2+2*X(2)2;z=fsxsteep(fx,0.001,1,1)

    展开阅读全文
    提示  道客多多所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    关于本文
    本文标题:MATLAB最速下降法(梯度法)程序.doc
    链接地址:https://www.docduoduo.com/p-7654049.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    道客多多用户QQ群:832276834  微博官方号:道客多多官方   知乎号:道客多多

    Copyright© 2025 道客多多 docduoduo.com 网站版权所有世界地图

    经营许可证编号:粤ICP备2021046453号    营业执照商标

    1.png 2.png 3.png 4.png 5.png 6.png 7.png 8.png 9.png 10.png



    收起
    展开