quasifun.m

来自「Hansen最新（2008）修订后的关于反演正则化方法的电子书与相应Matlab」· M 代码 · 共 11 行

function Q = quasifun(lambda,s,xi,dsvd)

% Auxiliary routine for quasiopt.  PCH, IMM, 12/29/97.

if (nargin==3)
   f = (s.^2)./(s.^2 + lambda^2);
else
   f = s./(s + lambda);
end

Q = norm((1 - f).*f.*xi);