MATLAB最速下降法梯度法程序.doc资源-CSDN文库

版权申诉

175 浏览量 2022-07-06 01:04:36 上传评论收藏 44KB DOC 举报

资源推荐

资源详情

资源评论

　　　

　

　　　

查看文章　　

　　

　　　

matlab 最速下降法　

2010-08-18 17:13　

function x=fsxsteep(f,e,a,b)

% fsxsteep 函数最速下降法

% x=fsxsteep(f,e,a,b)为输入函数 f 为函数 e 为允许误差 (a,b)为初始点;　　

% fsx TJPU 2008.6.15　　

x1=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))>=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.15　

x1=a;x2=b;

fx=diff(f,'x1'); %求 f 对 x1 偏导数

fy=diff(f,'x2'); %求 f 对 x2 偏导数

fxx=diff(fx,'x1'); %求二阶偏导数对 x1 再对 x1

fxy=diff(fx,'x2'); %求二阶偏导数对 x1 再对 x2

fyx=diff(fy,'x1'); %求二阶偏导数对 x2 再对 x1

fyy=diff(fy,'x2'); %求二阶偏导数对 x2 再对 x2

本内容试读结束，登录后可阅读更多

下载后可阅读完整内容，剩余2页未读，立即下载

内容反馈

版权申诉

老帽爬新坡

粉丝: 82
资源: 2万+

最新资源

资源上传下载、课程学习等过程中有任何疑问或建议，欢迎提出宝贵意见哦~我们会及时处理！点击此处反馈

feedback-tip