Line | |
---|
1 | % Reads the file |
---|
2 | K = textread('sectormap', '%f', 'headerlines', 9); |
---|
3 | M = K(3:260)'; % reads K + R + T of RF-Multipole |
---|
4 | |
---|
5 | % matrix J |
---|
6 | J = [ 0 1 0 0 0 0 ; -1 0 0 0 0 0 ; 0 0 0 1 0 0 ; 0 0 -1 0 0 0 ; 0 ... |
---|
7 | 0 0 0 0 1 ; 0 0 0 0 -1 0 ] |
---|
8 | |
---|
9 | % extract the kick |
---|
10 | K = M(1:6); |
---|
11 | |
---|
12 | % extract the R matrix |
---|
13 | R = reshape(M(6+(1:36)), 6, 6) |
---|
14 | |
---|
15 | % extract the T matrix |
---|
16 | T = {}; |
---|
17 | I = zeros(1, 216); |
---|
18 | index=1; |
---|
19 | for k=0:5 |
---|
20 | for j=0:5 |
---|
21 | for i=0:5 |
---|
22 | I(index++) = 6 + 36 + 1 + (i*36) + (j*6) + k; |
---|
23 | end |
---|
24 | end |
---|
25 | end |
---|
26 | |
---|
27 | for k=0:5 |
---|
28 | T{k+1} = zeros(6,6); |
---|
29 | for i=0:5 |
---|
30 | for j=0:5 |
---|
31 | T{k+1}(i+1,j+1) = M(I((i*36) + (j*6) + k + 1)); |
---|
32 | end |
---|
33 | end |
---|
34 | end |
---|
35 | |
---|
36 | disp('Test 1: transpose(R) * J * R == J') |
---|
37 | |
---|
38 | R' * J * R == J |
---|
39 | |
---|
40 | disp(['Test 2: transpose(Tk) * J * R + transpose(R) * J * Tk == 0 ' ... |
---|
41 | 'for k=1..6']); |
---|
42 | |
---|
43 | for k=1:6 |
---|
44 | k |
---|
45 | T{k}' * J * R + R' * J * T{k} |
---|
46 | end |
---|
47 | |
---|
48 | keyboard |
---|
Note: See
TracBrowser
for help on using the repository browser.