tv_driver.m

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

%  tvdriver.m
%

%  Initialization and setting of default stopping tolerances
%  is done in tv_param.m.

DEBUG=0; % 1 for Diagnostic output.

  newt_params=[newt_maxiter;newt_steptol;newt_gradtol;newt_tab_flag;newt_fig];
  cg_params = [cg_maxiter; cg_steptol; cg_residtol; cg_tab_flag; cg_fig];
%  Solve TV unconstrained min problem.
if DEBUG==1,
fprintf('newt_maxiter=%d\nnewt_steptol=%d\nnewt_gradtol=%d\nnewt_tab_flag=%d\nnewt_fig=%d\n\n\n',newt_params)
fprintf('cg_maxiter=%d\ncg_steptol=%d\ncg_residtol=%d\ncg_tab_flag=%d\ncg_fig=%d\n\n\n',cg_params)
end


  if pcgflag > 4
    disp(' *** Error.  Illegal pcg option. ***');
  elseif pcgflag == 0         %  CG without preconditioning.
    %  Leave cg_params vector as is.
  elseif pcgflag == 1         %  Operator Splitting PCG.
    cg_params = [cg_params; gam];
  else                        %  MultiLevel PCG
    p = log2(nx);
    n0 = 2^p0;
    cg_params = [cg_params; pcgflag-1; p0; Q_iter; nu];
  end
if DEBUG == 1,
	alpha,beta,
	pcgflag,
	cg_params
end  
  u_tv = pls_uncnst_tv(uh,bh,k_hat,alpha,beta,newt_params,pcgflag,cg_params);