mult_hbar.m

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

  function Hbar_u = mult_hbar(umat,khat,alpha,L,Active)

%  Hbar_u = mult_hbar(umat,khat,alpha,Active)
%
%  Form product Hbar*u, where
%      Hbar = alpha*diag(Active) + diag(Inactive)*H*diag(Inactive),
%      H = K'*K + alpha*L,
%      Active is an array of 0's and 1's. 1's indicate the constraint
%         u(x,y) = 0 is active.
%      Inactive is an array of 0's and 1's. 1's indicate inactive constraints.
%      K'*K is SPD and block Toeplitz with Toeplitz blocks.

  [nx,ny] = size(umat);
  u_Inactive = (1-Active) .* umat;
  Ku_Inactive = integral_op(u_Inactive,khat);
  KstarKu_Inactive = integral_op(Ku_Inactive,conj(khat));
  Lu_Inactive = reshape( L*u_Inactive(:), nx, ny );
  Hbar_u = alpha * Active.*umat + ...
      (1-Active).*(KstarKu_Inactive + alpha*Lu_Inactive);