/ani/mrses

To get this branch, use:
bzr branch http://suren.me/webbzr/ani/mrses

« back to all changes in this revision

Viewing changes to test.m

  • Committer: Suren A. Chilingaryan
  • Date: 2010-04-28 04:30:08 UTC
  • Revision ID: csa@dside.dyndns.org-20100428043008-vd9z0nso9axezvlp
Initial import

Show diffs side-by-side

added added

removed removed

Lines of Context:
test.m
 
1
A=double((unifrnd(0,1,50,2048)));
 
2
B=double((unifrnd(0,1,40,2048)));
 
3
A(:,5) = A(:,5) + 9; A(:,7) = A(:,7) + 9; A(:,9) = A(:,9) + 9; A(:,11) = A(:,11) + 9; A(:,13) = A(:,13) + 9;
 
4
%matrix division would fail due 0 determinant if sum of covariance matrixes is constant for selected properties
 
5
%A(:,5) = 9; A(:,7) = 9; A(:,9) = 9; A(:,11) = 9; A(:,13) = 9;
 
6
%B(:,5) = 0; B(:,7) = 0; B(:,9) = 0; B(:,11) = 0; B(:,13) = 0;
 
7
%mrses(A,B,5,1,60);
 
8
%mrses(A,B,5,10,60);
 
9
%mrses(A,B,5,500,1000)
 
10
%res=mrses(A,B,5,1,1000);
 
11
 
 
12
%res=mrses(A,B,5,10,500);
 
13
%res=mrses(A,B,5,500,1000);
 
14
%res=mrses_orig(A,B,5,500,1000);
 
15
%res=mrses_octave(A,B,5,500,1000);
 
16
 
 
17
%{
 
18
res=mrses_hw_distance(A,B,5,500,1000);
 
19
 
 
20
fprintf('\n')
 
21
disp(' top properties / percentage')
 
22
fprintf('%8u ', int32(res(1:10,2))')
 
23
fprintf('\n')
 
24
fprintf('%7.2f%% ', 100*res(1:10,1)')
 
25
fprintf('\n\n')
 
26
%}
 
27
 
 
28
res=mrses_mtx(A,B,5,500,1000);
 
29
(res+1)'