转载 matlab最速下降法 梯度下降算法 matlab

function x=fsxsteep(f,e,a,b)
% fsxsteep函数 最速下降法
[转载]matlab最速下降法 梯度下降算法 matlab
% 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
fxx=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)

原文链接:http://blog.163.com/greatboy@126/blog/static/5786079120107185130807/

  

爱华网本文地址 » http://www.413yy.cn/a/25101016/306704.html

更多阅读

qq飞车闪耀的皇族第7关怎么过 qq飞车燕子坞最快跑法

qq飞车闪耀的皇族第7关怎么过——简介主要是在多练习 , 首先不能有失误,最重要的是近道是关键。第一圈别用计时跑法,也就是视频那种。有气就W喷放了。我刚到第二圈就放气超了他们。只后不失误就行了。很抱歉,视频数据暂时无法显示

最有效减肥法从早餐开始瘦 最快最有效的减肥法

如果你想要成功减肥,早餐就一定要吃。不吃早餐不会让你瘦下来,只会毁掉你的减肥计划。下面,小编给你推荐5大最有效的早餐减肥法,助你轻松健康享瘦哦!最有效减肥法从早餐开始瘦——步骤/方法最有效减肥法从早餐开始瘦 1、不吃早饭更容易肥

声明:《转载 matlab最速下降法 梯度下降算法 matlab》为网友唯彩分享!如侵犯到您的合法权益请联系我们删除