maxrouting.m

来自「这是一个超宽带的源码」· M 代码 · 共 13 行

function [val,der,H] = maxrouting(x,m,nf)
eps = 1e-6;
f = x(1:nf);
val = -sum(log(eps+f));

if nargout > 1
	der = [-1./(eps+f);zeros(m,1)];
	if nargout > 2
		H = diag(1./((eps+f).^2));
		H = [H,zeros(nf,m);zeros(m,nf),zeros(m,m)];
	end
end