function?out?=?deconvtv（g?H?mu?opts）
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%?out?=?deconvtvl1（g?H?mu?opts）
%?deconvolves?image?g?by?solving?the?following?TV?minimization?problem
%
%?min??mu??||?Hf?-?g?||_1?+?||f||_TV
%?min?mu/2?||?Hf?-?g?||^2?+?||f||_TV
%
%?where?||f||_TV?=?sum_{xyt}?sqrt（?a||Dxf||^2?+?b||Dyf||^2?+?c||Dtf||^2）
%?Dxf?=?f（x+1y?t）?-?f（xyt）
%?Dyf?=?f（xy+1?t）?-?f（xyt）
%?Dtf?=?f（xy?t+1）?-?f（xyt）
%
%?Input:??????g??????-?the?observed?image?can?be?gray?scale?color?or?images
%?????????????H??????-?point?spread?function
%????????????mu??????-?regularization?parameter
%?????opts.method????-?either?‘l1‘?or?{‘l2‘}
%?????opts.rho_r?????-?initial?penalty?parameter?for?||u-Df||???{2}
%?????opts.rho_o?????-?initial?penalty?parameter?for?||Hf-g-r||?{50}
%?????opts.beta??????-?regularization?parameter?[a?b?c]?for?weighted?TV?norm?{[1?1?0]}
%?????opts.gamma?????-?update?constant?for?rho_r?{2}
%?????opts.max_itr???-?maximum?iteration?{20}
%?????opts.alpha?????-?constant?that?determines?constraint?violation?{0.7}
%?????opts.tol???????-?tolerance?level?on?relative?change?{1e-3}
%?????opts.print?????-?print?screen?option?{false}
%?????opts.f?????????-?initial??f?{g}
%?????opts.y1????????-?initial?y1?{0}
%?????opts.y2????????-?initial?y2?{0}
%?????opts.y3????????-?initial?y3?{0}
%?????opts.z?????????-?initial??z?{0}
%?????**?default?values?of?opts?are?given?in?{?}.
%
%?Output:?out.f??????-?output?images
%?????????out.itr????-?total?number?of?iterations?elapsed
%?????????out.relchg?-?final?relative?change
%?????????out.Df1????-?Dxf?f?is?the?output?images
%?????????out.Df2????-?Dyf?f?is?the?output?images
%?????????out.Df3????-?Dtf?f?is?the?output?images
%?????????out.y1?????-?Lagrange?multiplier?for?Df1
%?????????out.y2?????-?Lagrange?multiplier?for?Df2
%?????????out.y3?????-?Lagrange?multiplier?for?Df3
%?????????out.rho_r??-?final?penalty?parameter
%?%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
path（pathgenpath（pwd））;

if?nargin<3
????error（‘not?enough?inputs?try?again?\n‘）;
elseif?nargin==3
????opts?=?[];
end

if?~isnumeric（mu）
????error（‘mu?must?be?a?numeric?value!?\n‘）;
end

[rowscolsframes]?=?size（g）;
memory_condition?=?memory;
max_array_memory?=?memory_condition.MaxPossibleArrayBytes/16;
if?rows*cols*frames>0.1*max_array_memory
????fprintf（‘Warning:?possible?memory?issue?\n‘）;
????reply?=?input（‘Do?you?want?to?continue??[y/n]:?‘?‘s‘）;
????if?isequal（reply?‘n‘）
????????out.f?=?0;
????????return
????end
end

if?~isfield（opts‘method‘）
????method?=?‘l2‘;
else
????method?=?opts.method;
end

switch?method
????case?‘l2‘
????????out?=?deconvtvl2（gHmuopts）;
????case?‘l1‘
????????out?=?deconvtvl1（gHmuopts）;
????otherwise
????????error（‘unknown?method?\n‘）;
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????6489??2018-04-26?16:26??deconvtvl2.m

?????文件????????917??2018-04-25?20:59??Example_image_deblur.m

?????文件????????962??2018-04-25?20:55??Example_image_denoise.m

?????文件????????412??2018-04-25?22:02??psnr.m

?????文件??????76309??2011-01-21?06:52??data\building.jpg

?????文件??????12708??2011-01-21?07:02??data\data0001.jpg

?????文件??????12855??2011-01-21?07:02??data\data0002.jpg

?????文件??????12599??2011-01-21?07:02??data\data0003.jpg

?????文件??????12720??2011-01-21?07:02??data\data0004.jpg

?????文件??????12508??2011-01-21?07:02??data\data0005.jpg

?????文件??????13025??2011-01-21?07:02??data\data0006.jpg

?????文件??????12927??2011-01-21?07:02??data\data0007.jpg

?????文件??????12748??2011-01-21?07:02??data\data0008.jpg

?????文件??????12901??2011-01-21?07:02??data\data0009.jpg

?????文件??????12601??2011-01-21?07:02??data\data0010.jpg

?????文件?????387443??2018-04-18?09:52??data\lena_color_512.png

?????文件?????153785??2016-11-29?00:34??data\lena_gray_512.png

?????文件???????8806??2011-01-21?06:53??data\wind.jpg

?????文件?????105816??2018-04-25?20:20??output\guassdenoise.jpg

?????文件??????45866??2018-04-25?21:00??output\output-001.jpg

?????文件???????2814??2018-04-25?21:42??deconvtv.m

?????文件???????6707??2018-04-26?16:26??deconvtvl1.m

?????目錄??????????0??2018-05-20?22:29??data

?????目錄??????????0??2018-05-20?22:29??output

-----------?---------??----------?-----??----

???????????????923918????????????????????24