梯度算法.doc

function y=fs2steep(f,e,a,b) %返回的是点坐标的2个分量 2 % fs2steep函数 最速下降法 3 % x=fs2steep(f,e,a,b)为输入函数 f为函数 e为允许误差 (a,b)为初始点; 4 % fsx TJPU 2008.6.15 5 x1=a;x2=b; 6 Q=fs2hesse(f,x1,x2); 7 x0=[x1 x2]; 8 fx1=diff(f,x1); %对x1求偏导数 9 fx2=diff(f,x2); %对x2求偏导数 10 g=[fx1 fx2]; %梯度 11 g1=subs(g); %把符号变量转为数值 12 d=-g1;%d为搜索方向 13 while (abs(norm(g1))=e) %norm(g1)为g1的2范数，即sqrt(x1^2+x2^2)，因为梯度其各分量=0，所以其梯度幅值=0 14 t=(-d)*d/((-d)*Q*d);%求搜索步长，方法是？ 15 x0=x0-t*g1; %搜索到的点 16 v=x0; 17 a=[1 0]*x0; 18 b=[0 1]*x0; 19 x1=a; 20 x2=b; 21 Q=fs2hesse(f,x1,x2); 22 x0=[x1 x2]; 23 fx1=diff(f,x1); %对x1求偏导数 24 fx2=diff(f,x2); %对x2求偏导数 25 g=[fx1 fx2]; %梯度 26 g1=subs(g); 27 d=-g1; 28 end; 29 y=v; 30 31 function x=fs2hesse(f,a,b) 32 % fs2hesse函数 求函数的hesse矩阵； 33 % 本程序仅是简单的求二次函数的hesse矩阵！； 34 % x=fs2hesse(f)为输入函数 f为二次函数 x1,x2为自变量； 35 % fsx TJPU 2008.6.15 36 x1=a;x2=b; 37 fx=diff(f,x1); %求f对x1偏导数 38 fy=diff(f,x2); %求f对x2偏导数 39 fxx=diff(fx,x1); %求二阶偏导数 对x1再对x1 40 fxy=diff(fx,x2); %求二阶偏导数 对x1再对x2 41 fyx=diff(fy,x1); %求二阶偏导数 对x2再对x1 42 fyy=diff(fy,x2); %求二阶偏导数 对x2再对x2 43 fxx=subs(fxx); %将符号变量转化为数值 44 fxy=subs(fxy); 45 fyx=subs(fyx); 46 fyy=subs(fyy); 47 x=[fxx,fxy;fyx,fyy]; %求hesse矩阵