1 files changed, 57 insertions, 0 deletions
diff --git a/l1dantzig_example.m b/l1dantzig_example.m
new file mode 100644
index 0000000..0de90e1
--- /dev/null
+++ b/l1dantzig_example.m
@@ -0,0 +1,57 @@
+% l1dantzig_example.m
+%
+% Test out l1dantzig code (l1 minimization with bounded residual correlation).
+%
+% Written by: Justin Romberg, Caltech
+% Email: jrom@acm.caltech.edu
+% Created: October 2005
+%
+
+% put optimization code in path if not already there
+path(path, './Optimization');
+
+% signal length
+N = 512;
+
+% number of spikes to put down
+T = 20;
+
+% number of observations to make
+K = 120;
+
+% random +/- 1 signal
+x = zeros(N,1);
+q = randperm(N);
+x(q(1:T)) = sign(randn(T,1));
+
+% measurement matrix = random projection
+disp('Creating measurment matrix...');
+A = randn(K,N);
+A = orth(A')';
+disp('Done.');
+	
+% noisy observations
+sigma = 0.005;
+e = sigma*randn(K,1);
+y = A*x + e;
+
+% initial guess = min energy
+x0 = A'*y;
+
+% Dantzig selection
+epsilon = 3e-3;
+tic
+xp = l1dantzig_pd(x0, A, [], y, epsilon, 5e-2);
+toc
+
+
+% large scale
+% Afun = @(z) A*z;
+% Atfun = @(z) A'*z;
+% tic
+% xp = l1dantzig_pd(x0, Afun, Atfun, y, epsilon, 5e-2, 50, 1e-8, 500);
+% toc
+
+
+
+