%--------------------------------------------------------------------------
%?This?function?takes?a?DxN?matrix?of?N?data?points?in?a?D-dimensional?
%?space?and?returns?a?NxN?coefficient?matrix?of?the?sparse?representation?
%?of?each?data?point?in?terms?of?the?rest?of?the?points
%?Y:?DxN?data?matrix
%?affine:?if?true?then?enforce?the?affine?constraint
%?thr1:?stopping?threshold?for?the?coefficient?error?||Z-C||
%?thr2:?stopping?threshold?for?the?linear?system?error?||Y-YZ||
%?maxIter:?maximum?number?of?iterations?of?ADMM
%?C2:?NxN?sparse?coefficient?matrix
%--------------------------------------------------------------------------
%?Copyright?@?Ehsan?Elhamifar?2012
%--------------------------------------------------------------------------

function?C2?=?admmLasso_mat_func（YaffinealphathrmaxIter）

if?（nargin?????%?default?subspaces?are?linear
????affine?=?false;?
end
if?（nargin?????%?default?regularizarion?parameters
????alpha?=?800;
end
if?（nargin?????%?default?coefficient?error?threshold?to?stop?ADMM
????%?default?linear?system?error?threshold?to?stop?ADMM
????thr?=?2*10^-4;?
end
if?（nargin?????%?default?maximum?number?of?iterations?of?ADMM
????maxIter?=?200;?
end

if?（length（alpha）?==?1）
????alpha1?=?alpha（1）;
????alpha2?=?alpha（1）;
elseif?（length（alpha）?==?2）
????alpha1?=?alpha（1）;
????alpha2?=?alpha（2）;
end

if?（length（thr）?==?1）
????thr1?=?thr（1）;
????thr2?=?thr（1）;
elseif?（length（thr）?==?2）
????thr1?=?thr（1）;
????thr2?=?thr（2）;
end

N?=?size（Y2）;

%?setting?penalty?parameters?for?the?ADMM
mu1?=?alpha1?*?1/computeLambda_mat（Y）;
mu2?=?alpha2?*?1;

if?（~affine）
????%?initialization
????A?=?inv（mu1*（Y‘*Y）+mu2*eye（N））;
????C1?=?zeros（NN）;
????Lambda2?=?zeros（NN）;
????err1?=?10*thr1;?err2?=?10*thr2;
????i?=?1;
????%?ADMM?iterations
????while?（?err1（i）?>?thr1?&&?i?????????%?updating?Z
????????Z?=?A?*?（mu1*（Y‘*Y）+mu2*（C1-Lambda2/mu2））;
????????Z?=?Z?-?diag（diag（Z））;
????????%?updating?C
????????C2?=?max（0（abs（Z+Lambda2/mu2）?-?1/mu2*ones（N）））?.*?sign（Z+Lambda2/mu2）;
????????C2?=?C2?-?diag（diag（C2））;
????????%?updating?Lagrange?multipliers
????????Lambda2?=?Lambda2?+?mu2?*?（Z?-?C2）;
????????%?computing?errors
????????err1（i+1）?=?errorCoef（ZC2）;
????????err2（i+1）?=?errorLinSys（YZ）;
????????%
????????C1?=?C2;
????????i?=?i?+?1;
????end
????fprintf（‘err1:?%2.4f?err2:?%2.4f?iter:?%3.0f?\n‘err1（end）err2（end）i）;
else
????%?initialization
????A?=?inv（mu1*（Y‘*Y）+mu2*eye（N）+mu2*ones（NN））;
????C1?=?zeros（NN）;
????Lambda2?=?zeros（NN）;
????lambda3?=?zeros（1N）;
????err1?=?10*thr1;?err2?=?10*thr2;?err3?=?10*thr1;
????i?=?1;
????%?ADMM?iterations
????while?（?（err1（i）?>?thr1?||?err3（i）?>?thr1）?&&?i?????????%?updating?Z
????????Z?=?A?*?（mu1*（Y‘*Y）+mu2*（C1-Lambda2/mu2）+mu2*ones（N1）*（ones（1N）-lambda3/mu2））;
????????Z?=?Z?-?diag（diag（Z））;
????????%?updating?C
????????C2?=?max（0（abs（Z+Lambda2/mu2）?-?1/mu2*ones（N）））?.*?sign（Z+Lambda2/mu2）;
????????C2?=?C2?-?diag（diag（C2））;
????????%?updating?Lagrange?multi

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-08-31?18:51??SSC_ADMM_v1.1\
?????文件????????6148??2012-10-03?01:46??SSC_ADMM_v1.1\.DS_Store
?????文件?????????969??2012-10-02?01:49??SSC_ADMM_v1.1\BuildAdjacency.m
?????文件?????????733??2012-10-04?23:58??SSC_ADMM_v1.1\DataProjection.m
?????文件?????????747??2012-10-02?01:43??SSC_ADMM_v1.1\Misclassification.m
?????文件????????1189??2012-10-04?23:59??SSC_ADMM_v1.1\SSC.m
?????文件?????????255??2018-08-31?18:52??SSC_ADMM_v1.1\SSC_Faces.mat
?????文件????????1218??2012-10-02?01:40??SSC_ADMM_v1.1\SpectralClustering.m
?????文件?????8018017??2012-06-29?16:19??SSC_ADMM_v1.1\YaleBCrop025.mat
?????文件????????3427??2012-10-05?00:01??SSC_ADMM_v1.1\admmLasso_mat_func.m
?????文件????????3980??2012-10-05?00:03??SSC_ADMM_v1.1\admmOutlier_mat_func.m
?????文件?????????678??2012-10-02?01:43??SSC_ADMM_v1.1\computeLambda_mat.m
?????文件?????????548??2012-10-02?01:43??SSC_ADMM_v1.1\errorCoef.m
?????文件?????????825??2012-10-02?01:43??SSC_ADMM_v1.1\errorLinSys.m
?????文件?????????580??2012-10-03?12:30??SSC_ADMM_v1.1\matrixNormalize.m
?????文件????????1099??2012-10-02?01:41??SSC_ADMM_v1.1\missclassGroups.m
?????文件???????17712??2012-10-05?00:26??SSC_ADMM_v1.1\readme.pdf
?????文件????????1558??2012-10-05?00:25??SSC_ADMM_v1.1\readme.rtf
?????文件????????1362??2012-10-05?00:08??SSC_ADMM_v1.1\run_SSC_Faces.m
?????文件????????2863??2012-10-05?00:28??SSC_ADMM_v1.1\run_SSC_MS.m
?????文件?????????680??2012-10-03?01:04??SSC_ADMM_v1.1\thrC.m