利用python进行矩阵运算

import numpy as np

R1 = np.array([[3.68, 11.26, 33.1, 35.4, 16.55], \
               [3.91, 11.49, 34.25, 33.56, 16.78], \
               [3.68, 15.17, 33.79, 31.72, 15.63], \
               [3.68, 8.97,	31.95, 35.4, 20]])
A1 = np.array([0.254, 0.253, 0.248, 0.245])
C1 = A1 @ R1

R2 = np.array([[5.06, 11.72, 36.32, 31.03, 15.86], \
               [4.37, 12.18, 37.01, 29.89, 16.55], \
               [5.52, 17.47,	38.39, 25.52, 13.1], \
               [4.6, 14.48, 36.09, 32.64, 12.18]])
A2 = np.array([0.255, 0.245, 0.254, 0.246])
C2 = A2 @ R2

R3 = np.array([[3.22, 10.34, 29.43,	38.62, 18.39], \
               [4.37, 13.79, 35.86,	28.97, 17.01], \
               [5.52, 14.02, 35.86, 30.11, 14.48], \
               [4.83, 16.55,	38.16, 26.21, 14.25], \
               [5.98, 19.77,	34.48, 27.13, 12.64], \
               [3.91, 8.05, 35.63, 32.64, 19.77], \
               [4.14, 11.72, 36.55, 33.79, 13.79], \
               [5.29, 16.32, 40,	26.21, 12.18]])
A3 = np.array([0.124, 0.124, 0.124, 0.128, 0.123, 0.125, 0.126, 0.127])
C3 = A3 @ R3

R4 = np.array([[4.83, 17.47, 38.85,	27.59, 11.26], \
               [7.13, 14.94, 39.54, 25.75, 12.64], \
               [5.06, 15.4,	36.32, 26.44, 16.78],\
               [5.29, 14.71, 32.41,	33.1, 14.48]])
A4 = np.array([0.25, 0.246,	0.261, 0.243])
C4 = A4 @ R4

R = np.array([C1, C2, C3, C4])
A = np.array([0.244, 0.244, 0.255, 0.258])
C = A @ R

print('C = A*R = ' + str(C))