1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
|
% l1decode_example.m
%
% Test out l1decode code.
%
% Written by: Justin Romberg, Caltech
% Email: jrom@acm.caltech.edu
% Created: October 2005
%
path(path, './Optimization');
% source length
N = 256;
% codeword length
M = 4*N;
% number of perturbations
T = round(.2*M);
% coding matrix
G = randn(M,N);
% source word
x = randn(N,1);
% code word
c = G*x;
% channel: perturb T randomly chosen entries
q = randperm(M);
y = c;
y(q(1:T)) = randn(T,1);
% recover
x0 = inv(G'*G)*G'*y;
xp = l1decode_pd(x0, G, [], y, 1e-4, 30);
% large scale
% gfun = @(z) G*z;
% gtfun = @(z) G'*z;
% xp = l1decode_pd(x0, gfun, gtfun, y, 1e-3, 25, 1e-8, 200);
|
