本帖最后由 zhangxueqing 于 2017-3-26 16:36 编辑
谢谢!但是如果不放循环里,那么图形画不出来?
t_start= tic;
%Global constants and defaults
QUIET = 0;
MAX_ITER = 1000;
ABSTOL = 1e-5;
RELTOL = 1e-3;
%Data preprocessing
[m, n] = size(A);
fprintf('(1) GADMM1 for lasso problem \n');
% save a matrix-vector multiply
Atb = A'*b;
%ADMM solver
x = zeros(n,1);
z = zeros(n,1);
u = zeros(n,1);
I=eye(n);
% cache the factorization
[L U] = factor(A, rho1);
if ~QUIET
fprintf('%3s\t%10s\t%10s\t%10s\t%10s\t%10s\n', 'iter', ...
'r norm', 'eps pri', 's norm', 'eps dual', 'objective');
end
for k = 1:MAX_ITER
% x-update
q = Atb &