conjgrad.m

斯坦福大学的小波的工具箱 比matlab自带的工具箱更完善

function [fnew,pnew,rnew,Lp] = ConjGrad(f,p,pold,r,u,J0,Lpold);

%Lp    	= L(u,FWT_Dyad(p,J0));
Lp    	= L(u.*FWT_ATrou(p,J0));
lambda 	= (r'*p) / (p'*Lp);
fnew 	= f + lambda * p;
rnew	= r - lambda * Lp;
pnew	= Lp - (Lp'*Lp) / (p'*Lp) * p;
if ~(pold == 0),
    pnew = pnew  - (Lp'*Lpold) / (pold'*Lpold) * pold;
end

% Written by Maureen Clerc and Jerome Kalifa, 1997
% [e-mail], [e-mail]
    
    
%   
% Part of WaveLab Version 802
% Built Sunday, October 3, 1999 8:52:27 AM
% This is Copyrighted Material
% For Copying permissions see COPYING.m
% Comments? e-mail [e-mail]
%