TP-optimisation-numerique/Hess_f_C14.m
2023-06-10 20:58:09 +02:00

33 lines
986 B
Matlab
Executable file

function [H_f, res, J_res] = Hess_f_C14(beta, donnees, residu, J_residu)
%
% Paramètres en entrés
% --------------------
% beta : vecteur des paramètres
% real(p)
% donnees : Données
% real(n,2)
% residu : fonction qui code les résidus
% res_beta = residus(beta)
% J_residu : fonction qui code la matrice jacobienne
% J_res_beta = J_residu(beta);
%
% Paramètres en sortie
% --------------------
% H_f : Matrice hessienne
% real(p,p)
% res : vecteur des résidus
% real(n)
% J_res : Matrice jacobienne des résiduis
% real(n,p)
%
res = residu(beta);
J_res = J_residu(beta);
% cf formule du cours
S = [ sum(0*donnees(:,1)) , sum(donnees(:,1).*exp(-beta(2)*donnees(:, 1)));
sum(donnees(:, 1).*exp(-beta(2)*donnees(:, 1))), sum(-beta(1)*(donnees(:, 1).^2).*exp(-beta(2)*donnees(:, 1))) ];
H_f = S + J_res.'*J_res;
end