1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51
| import numpy as np
R23 = np.array([ [0.18, 0.14, 0.18, 0.14, 0.13, 0.23], [0.15, 0.20, 0.15, 0.25, 0.10, 0.15], [0.25, 0.12, 0.13, 0.12, 0.18, 0.20], [0.16, 0.15, 0.21, 0.11, 0.20, 0.17], [0.23, 0.18, 0.17, 0.16, 0.15, 0.11], [0.19, 0.13, 0.12, 0.12, 0.11, 0.33], [0.17, 0.16, 0.15, 0.08, 0.25, 0.19]])
A23 = np.array([1/7, 1/7, 1/7, 1/7, 1/7, 1/7, 1/7])
B23 = np.dot(A23, R23)
R1 = np.array([ [0.12, 0.18, 0.17, 0.23, 0.13, 0.17], [0.15, 0.13, 0.18, 0.25, 0.12, 0.17], [0.14, 0.13, 0.16, 0.18, 0.20, 0.19], [0.12, 0.14, 0.15, 0.17, 0.19, 0.23], [0.16, 0.12, 0.13, 0.25, 0.18, 0.16]]) A1 = np.array([0.2, 0.2, 0.2, 0.2, 0.2]) B1 = np.dot(A1, R1)
R2 = np.array([ [0.13, 0.15, 0.14, 0.18, 0.16, 0.25], [0.12, 0.16, 0.13, 0.17, 0.19, 0.23], B23, [0.14, 0.13, 0.15, 0.16, 0.18, 0.24], [0.16, 0.15, 0.15, 0.17, 0.18, 0.19]]) A2 = np.array([0.2, 0.2, 0.2, 0.2, 0.2]) B2 = np.dot(A2, R2)
R3 = np.array([ [0.15, 0.14, 0.13, 0.18, 0.14, 0.26], [0.16, 0.15, 0.18, 0.14, 0.16, 0.21]]) A3 = np.array([0.5, 0.5]) B3 = np.dot(A3, R3)
R = np.array([B1, B2, B3]) A = np.array([1/3, 1/3, 1/3]) B = np.dot(A, R)
print(B)
|