bfgs_u2_f_celu_grad.asv

Control optimisation. It is example of use BFGS algorithm to control satellite form Earth atmosphere

% gradient funkcji celu wzgledem u2
% u2 == x
function q = bfgs_u2_f_celu_grad(x,main_x0,main_h0,main_tau)
h = 1e-5;
for i = 1:length(x)
    e = zeros(length(x),1);
    e(i) = 1;
    xeh = x+e*h;
    F1 = bfgs_f_celu2(main_tau',main_x0,main_h0,xeh);
    F2 = bfgs_f_celu(main_tau',main_x0,main_h0,x);
    Q_grad(i) = (F1-F2)/h;
end
q = Q_grad';

end % koniec funkcji