gradient.m

来自「共轭梯度法」· M 代码 · 共 6 行

function [g_aim]=gradient(x,y)
ddx=0.00000001;
ddy=0.00000001;
g_aim=zeros(2,1);
g_aim(1,1)=(forward(x+ddx,y)-forward(x,y))/ddx;
g_aim(2,1)=(forward(x,y+ddy)-forward(x,y))/ddy;