function?[ZE]?=?alm_lrr_l1（XAlambda）
%?This?routine?solves?the?following?nuclear-norm?optimization?problem
%?min?|Z|_*+lambda*|E|_1
%?s.t.?X?=?AZ+E
%?inputs:
%????????X?--?D*N?data?matrix?D?is?the?data?dimension?and?N?is?the?number
%?????????????of?data?vectors.
%????????A?--?D*M?matrix?of?a?dictionary?M?is?the?size?of?the?dictionary
?
tol?=?1e-8;
maxIter?=?1e6;
[d?n]?=?size（X）;
m?=?size（A2）;
rho?=?1.1;
max_mu?=?1e10;
mu?=?1e-6;
atx?=?A‘*X;
inv_a?=?inv（A‘*A+eye（m））;
%%?Initializing?optimization?variables
%?intialize
J?=?zeros（mn）;
Z?=?zeros（mn）;
E?=?sparse（dn）;
?
Y1?=?zeros（dn）;
Y2?=?zeros（mn）;
%%?Start?main?loop
iter?=?0;
disp（[‘initialrank=‘?num2str（rank（Z））]）;
while?iter????iter?=?iter?+?1;
????%update?J
????temp?=?Z?+?Y2/mu;
????