grad_ex6_1.m

来自「good code for matlab by mili , i than yo」· M 代码 · 共 19 行

function Return = Grad_Ex6_1(functname,x)
%
% numerical computation of gradient
% this allows automatic gradient computation
% 
%
% first forward finite difference
% hstep = 0.001; - programmed in
%
hstep = 0.001;
n = length(x);
f = feval(functname,x);

for i = 1:n
   xs(i) = x;
   xs(i) = xs(i) + hstep;
   gradx(i)= (feval(functname,xs) -f)/hstep;
end
Return = gradx;