%--------------------------------------------------------------------------
%?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

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

?????文件???????6148??2012-10-03?01:46??SSC_ADMM_v1.1\SSC_ADMM_v1.1\.DS_Store

?????文件???????3427??2012-10-05?00:01??SSC_ADMM_v1.1\SSC_ADMM_v1.1\admmLasso_mat_func.m

?????文件???????3980??2012-10-05?00:03??SSC_ADMM_v1.1\SSC_ADMM_v1.1\admmOutlier_mat_func.m

?????文件????????969??2012-10-02?01:49??SSC_ADMM_v1.1\SSC_ADMM_v1.1\BuildAdjacency.m

?????文件????????678??2012-10-02?01:43??SSC_ADMM_v1.1\SSC_ADMM_v1.1\computeLambda_mat.m

?????文件????????733??2012-10-04?23:58??SSC_ADMM_v1.1\SSC_ADMM_v1.1\DataProjection.m

?????文件????????548??2012-10-02?01:43??SSC_ADMM_v1.1\SSC_ADMM_v1.1\errorCoef.m

?????文件????????825??2012-10-02?01:43??SSC_ADMM_v1.1\SSC_ADMM_v1.1\errorLinSys.m

?????文件????????580??2012-10-03?12:30??SSC_ADMM_v1.1\SSC_ADMM_v1.1\matrixNormalize.m

?????文件????????747??2012-10-02?01:43??SSC_ADMM_v1.1\SSC_ADMM_v1.1\Misclassification.m

?????文件???????1099??2012-10-02?01:41??SSC_ADMM_v1.1\SSC_ADMM_v1.1\missclassGroups.m

?????文件??????17712??2012-10-05?00:26??SSC_ADMM_v1.1\SSC_ADMM_v1.1\readme.pdf

?????文件???????1558??2012-10-05?00:25??SSC_ADMM_v1.1\SSC_ADMM_v1.1\readme.rtf

?????文件???????1362??2012-10-05?00:08??SSC_ADMM_v1.1\SSC_ADMM_v1.1\run_SSC_Faces.m

?????文件???????2863??2012-10-05?00:28??SSC_ADMM_v1.1\SSC_ADMM_v1.1\run_SSC_MS.m

?????文件???????1218??2012-10-02?01:40??SSC_ADMM_v1.1\SSC_ADMM_v1.1\SpectralClustering.m

?????文件???????1189??2012-10-04?23:59??SSC_ADMM_v1.1\SSC_ADMM_v1.1\SSC.m

?????文件????????680??2012-10-03?01:04??SSC_ADMM_v1.1\SSC_ADMM_v1.1\thrC.m

?????文件????8018017??2012-06-29?16:19??SSC_ADMM_v1.1\SSC_ADMM_v1.1\YaleBCrop025.mat

?????文件?????????82??2012-10-03?01:46??SSC_ADMM_v1.1\__MACOSX\SSC_ADMM_v1.1\._.DS_Store

?????文件?????????82??2012-10-05?00:01??SSC_ADMM_v1.1\__MACOSX\SSC_ADMM_v1.1\._admmLasso_mat_func.m

?????文件?????????82??2012-10-05?00:03??SSC_ADMM_v1.1\__MACOSX\SSC_ADMM_v1.1\._admmOutlier_mat_func.m

?????文件????????717??2012-10-02?01:49??SSC_ADMM_v1.1\__MACOSX\SSC_ADMM_v1.1\._BuildAdjacency.m

?????文件????????720??2012-10-02?01:43??SSC_ADMM_v1.1\__MACOSX\SSC_ADMM_v1.1\._computeLambda_mat.m

?????文件????????369??2012-10-04?23:58??SSC_ADMM_v1.1\__MACOSX\SSC_ADMM_v1.1\._DataProjection.m

?????文件????????712??2012-10-02?01:43??SSC_ADMM_v1.1\__MACOSX\SSC_ADMM_v1.1\._errorCoef.m

?????文件????????714??2012-10-02?01:43??SSC_ADMM_v1.1\__MACOSX\SSC_ADMM_v1.1\._errorLinSys.m

?????文件????????718??2012-10-03?12:30??SSC_ADMM_v1.1\__MACOSX\SSC_ADMM_v1.1\._matrixNormalize.m

?????文件????????720??2012-10-02?01:43??SSC_ADMM_v1.1\__MACOSX\SSC_ADMM_v1.1\._Misclassification.m

?????文件????????718??2012-10-02?01:41??SSC_ADMM_v1.1\__MACOSX\SSC_ADMM_v1.1\._missclassGroups.m

............此處省略15個文件信息