nonneg_driver.m

采用matlab编写的数字图像恢复程序

%  nonnegdriver.m
%
  newt_params= ...
	[newt_maxiter;newt_steptol;newt_gradtol;newt_tab_flag;newt_fig;... 
	GA_param; ls_maxiter; tol0];

  cg_params = [cg_maxiter; cg_steptol; cg_residtol; cg_tab_flag; cg_fig];

  if reg_order == 0
    Lh = speye(N);
  else
    Lh = laplacian(nx,ny);
  end
  
  [uh,gradnormvec] = pls_nonneg(bh,k_hat,alpha,Lh,newt_params,cg_params);

  conv_factor = exp(mean(diff(log(gradnormvec))));
  fprintf(' Geometric mean of ||pgrad||/||old pgrad|| = %g.\n', conv_factor);