www.gusucode.com > 个用于图像标定的算法matlab源码程序 > code16/lk20_common/hessian.m
function H = hessian(VI_dW_dp, N_p, w) % HESSIAN - Compute Hessian % H = HESSIAN(VI_DW_DP, N_P, W) % Iain Matthews, Simon Baker, Carnegie Mellon University, Pittsburgh % $Id: hessian.m,v 1.1.1.1 2003/08/20 03:07:36 iainm Exp $ if nargin<3 error('Not enough input arguments'); end H = zeros(N_p, N_p); for i=1:N_p h1 = VI_dW_dp(:,((i-1)*w)+1:((i-1)*w)+w); for j=1:N_p h2 = VI_dW_dp(:,((j-1)*w)+1:((j-1)*w)+w); H(j, i) = sum(sum((h1 .* h2))); end end