% signal matrix 

%   (           L+      L-      l+      l-      ls  )                  
%   (SP, f1                                         )
%   (AP, f2                                         )
%   (SP, f1xf2                                      )
%   (AP, f1xf2                                      )
%   (PO, f1xf2                                      )
% s.kat is necessary


! perl mkat s

load s_run1.out
load s_run2.out
load s_run3.out
load s_run4.out
load s_run5.out

matrix_r=zeros(5);
for m = 1 : 5
matrix_r(m,1)=s_run1(1,1+m);
matrix_r(m,2)=s_run2(1,1+m);
matrix_r(m,3)=s_run3(1,1+m);
matrix_r(m,4)=s_run4(1,1+m);
matrix_r(m,5)=s_run5(1,1+m);
end    

matrix_n=zeros(5);
for m = 1 : 5
matrix_n(1,m)=matrix_r(1,m)/matrix_r(1,1);
matrix_n(2,m)=matrix_r(2,m)/matrix_r(2,2);
matrix_n(3,m)=matrix_r(3,m)/matrix_r(3,3);
matrix_n(4,m)=matrix_r(4,m)/matrix_r(4,4);
matrix_n(5,m)=matrix_r(5,m)/matrix_r(5,5);
end    

dlmwrite('matrix_raw.txt',matrix_r, '\t')
dlmwrite('matrix_norm.txt',matrix_n, '\t')


