1 % signal matrix 
   2 
   3 %   (           L+      L-      l+      l-      ls  )                  
   4 %   (SP, f1                                         )
   5 %   (AP, f2                                         )
   6 %   (SP, f1xf2                                      )
   7 %   (AP, f1xf2                                      )
   8 %   (PO, f1xf2                                      )
   9 % s.kat is necessary
  10 
  11 
  12 ! perl mkat s
  13 
  14 load s_run1.out
  15 load s_run2.out
  16 load s_run3.out
  17 load s_run4.out
  18 load s_run5.out
  19 
  20 matrix_r=zeros(5);
  21 for m = 1 : 5
  22 matrix_r(m,1)=s_run1(1,1+m);
  23 matrix_r(m,2)=s_run2(1,1+m);
  24 matrix_r(m,3)=s_run3(1,1+m);
  25 matrix_r(m,4)=s_run4(1,1+m);
  26 matrix_r(m,5)=s_run5(1,1+m);
  27 end    
  28 
  29 matrix_n=zeros(5);
  30 for m = 1 : 5
  31 matrix_n(1,m)=matrix_r(1,m)/matrix_r(1,1);
  32 matrix_n(2,m)=matrix_r(2,m)/matrix_r(2,2);
  33 matrix_n(3,m)=matrix_r(3,m)/matrix_r(3,3);
  34 matrix_n(4,m)=matrix_r(4,m)/matrix_r(4,4);
  35 matrix_n(5,m)=matrix_r(5,m)/matrix_r(5,5);
  36 end    
  37 
  38 dlmwrite('matrix_raw.txt',matrix_r, '\t')
  39 dlmwrite('matrix_norm.txt',matrix_n, '\t')