/ani/mrses

To get this branch, use:
bzr branch http://darksoft.org/webbzr/ani/mrses
1 by Suren A. Chilingaryan
Initial import
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)'