annotate tests/RecomTSTgentest.m @ 67:a8d96e67717e
Added the Analysis-By-Synthesis algorithms used in the papers "Analysis-based sparse reconstruction with synthesis-based solvers", "Choosing Analysis or Synthesis Recovery for Sparse Reconstruction" and "A generalization of synthesis and analysis sparsity"
author |
Nic Cleju <nikcleju@gmail.com> |
date |
Tue, 09 Jul 2013 14:21:10 +0300 |
parents |
4a4e5204ecf5 |
children |
|
rev |
line source |
nikcleju@4
|
1 % Run BP and save parameters and solutions as reference test data
|
nikcleju@4
|
2 % to check if other algorithms are correct
|
nikcleju@4
|
3
|
nikcleju@5
|
4 % numA = 10;
|
nikcleju@5
|
5 % numY = 100;
|
nikcleju@5
|
6 %
|
nikcleju@5
|
7 % sizesA{1} = [50 100];
|
nikcleju@5
|
8 % sizesA{2} = [20 25];
|
nikcleju@5
|
9 % sizesA{3} = [10 120];
|
nikcleju@5
|
10 % sizesA{4} = [15 100];
|
nikcleju@5
|
11 % sizesA{5} = [70 100];
|
nikcleju@5
|
12 % sizesA{6} = [80 100];
|
nikcleju@5
|
13 % sizesA{7} = [90 100];
|
nikcleju@5
|
14 % sizesA{8} = [99 100];
|
nikcleju@5
|
15 % sizesA{9} = [100 100];
|
nikcleju@5
|
16 % sizesA{10} = [250 400];
|
nikcleju@5
|
17 %
|
nikcleju@5
|
18 % for i = 1:numA
|
nikcleju@5
|
19 % sz = sizesA{i};
|
nikcleju@5
|
20 % cellA{i} = randn(sz);
|
nikcleju@5
|
21 % cellY{i} = randn(sz(1), numY);
|
nikcleju@5
|
22 % for j = 1:numY
|
nikcleju@5
|
23 % cellTol{i}(j) = rand / 5; % restrict from 0 to 20% if measurements
|
nikcleju@5
|
24 % %cellX0{i}(:,j) = cellA{i} \ cellY{i}(:,j);
|
nikcleju@5
|
25 % cellX0{i}(:,j) = zeros(size(cellA{i},2),1);
|
nikcleju@5
|
26 % end
|
nikcleju@5
|
27 % end
|
nikcleju@5
|
28 % %load BPtestdata
|
nikcleju@5
|
29 tic
|
nikcleju@4
|
30 for i = 1:numA
|
nikcleju@4
|
31 for j = 1:numY
|
nikcleju@4
|
32 cellXr{i}(:,j) = RecommendedTST(cellA{i}, cellY{i}(:,j), 300, cellTol{i}(j), cellX0{i}(:,j));
|
nikcleju@4
|
33 end
|
nikcleju@4
|
34 end
|
nikcleju@5
|
35 toc
|
nikcleju@4
|
36
|
nikcleju@5
|
37 % save RecomTSTtestdata |