Решите задачу без якобиана

Функция помощника myfun заданный в конце этого примера реализует векторно оцененную целевую функцию без информации о производной. Решить минимизацию начиная с точки x0.

x0 = [0.3,0.4]; % Starting guess
[x,resnorm,res,eflag,output] = lsqnonlin(@myfun,x0); % Invoke optimizer

Local minimum possible.
lsqnonlin stopped because the size of the current step is less than
the value of the step size tolerance.

Исследуйте решение и количество вычислений функции.

disp(x)

    0.2578    0.2578

disp(resnorm)

  124.3622

disp(output.funcCount)

    72

Решите задачу включая якобиан

Целевая функция достаточно проста, чтобы можно было вычислить ее якобиан. Следуя определению в якобианах Вектора Functions, якобианская функция представляет матрицу

Вот, $$F_k(x)$$ является k-м компонентом целевой функции. Этот пример имеет

так

Функция помощника myfun2 заданный в конце этого примера реализует целевую функцию с якобианом. Установите опции, чтобы решатель использовал якобиан.

opts = optimoptions(@lsqnonlin,'SpecifyObjectiveGradient',true);

Запустите решатель.

lb = []; % No bounds
ub = [];
[x2,resnorm2,res2,eflag2,output2] = lsqnonlin(@myfun2,x0,lb,ub,opts);

Local minimum possible.
lsqnonlin stopped because the size of the current step is less than
the value of the step size tolerance.

Решение аналогично предыдущему решению.

disp(x2)

    0.2578    0.2578

disp(resnorm2)

  124.3622

Преимущество использования якобиана в том, что решатель принимает намного меньше вычислений функции.

disp(output2.funcCount)

    24

Вспомогательные функции

Этот код создает myfun вспомогательная функция.

function F = myfun(x)
k = 1:10;
F = 2 + 2*k-exp(k*x(1))-exp(k*x(2));
end

Этот код создает myfun2 вспомогательная функция.

function [F,J] = myfun2(x)
k = 1:10;
F = 2 + 2*k-exp(k*x(1))-exp(k*x(2));
if nargout > 1
    J = zeros(10,2);
    J(k,1) = -k.*exp(k*x(1));
    J(k,2) = -k.*exp(k*x(2));
end
end