Mercurial > hg > pycsalgos
comparison tests/RecomTSTgentest.m @ 4:4393ad5bffc1
Implemented the RecommendedTST algorithm
Test file written, but test not yet working.
author | nikcleju |
---|---|
date | Mon, 24 Oct 2011 23:39:53 +0000 |
parents | |
children | 4a4e5204ecf5 |
comparison
equal
deleted
inserted
replaced
3:537f7798e186 | 4:4393ad5bffc1 |
---|---|
1 % Run BP and save parameters and solutions as reference test data | |
2 % to check if other algorithms are correct | |
3 | |
4 numA = 10; | |
5 numY = 100; | |
6 | |
7 sizesA{1} = [50 100]; | |
8 sizesA{2} = [20 25]; | |
9 sizesA{3} = [10 120]; | |
10 sizesA{4} = [15 100]; | |
11 sizesA{5} = [70 100]; | |
12 sizesA{6} = [80 100]; | |
13 sizesA{7} = [90 100]; | |
14 sizesA{8} = [99 100]; | |
15 sizesA{9} = [100 100]; | |
16 sizesA{10} = [250 400]; | |
17 | |
18 for i = 1:numA | |
19 sz = sizesA{i}; | |
20 cellA{i} = randn(sz); | |
21 cellY{i} = randn(sz(1), numY); | |
22 for j = 1:numY | |
23 cellTol{i}(j) = rand / 5; % restrict from 0 to 20% if measurements | |
24 %cellX0{i}(:,j) = cellA{i} \ cellY{i}(:,j); | |
25 cellX0{i}(:,j) = zeros(size(cellA{i},2),1); | |
26 end | |
27 end | |
28 %load BPtestdata | |
29 | |
30 for i = 1:numA | |
31 for j = 1:numY | |
32 cellXr{i}(:,j) = RecommendedTST(cellA{i}, cellY{i}(:,j), 300, cellTol{i}(j), cellX0{i}(:,j)); | |
33 end | |
34 end | |
35 | |
36 save RecomTSTtestdata |