import numpy as np
np.set_printoptions(precision=10)
mul=np.array([[ 2.79358336e-03, -4.33801272e-03,  1.49405606e-01, -1.49494674e-01],
[-4.33801272e-03,  6.73628032e-03, -2.32004399e-01,  2.32142708e-01],
[ 1.49405606e-01, -2.32004399e-01,  7.99047014e+00, -7.99523365e+00],
[-1.49494674e-01,  2.32142708e-01, -7.99523365e+00 , 8.00000000e+00]])
length = np.linalg.matrix_rank(mul) ##求秩
print(length)
#print(np.linalg.inv(mul))
ceshi = np.round(np.linalg.inv(mul),decimals=9)
print(ceshi)
ni = np.array([[1.15688790e+09, 8.31687162e+08, 6.45472327e+08, 6.42572567e+08],
[8.31687176e+08, 2.31447075e+09, 1.72043424e+09, 1.66778988e+09],
[6.45472390e+08, 1.72043428e+09, 5.47332438e+09, 5.43220196e+09],
[6.42572630e+08, 1.66778992e+09, 5.43220196e+09, 5.39257748e+09]])
# print(mul*ni)
# print(ni*mul)
a=np.array([[1,2],
            [2,3]])
b=np.array([[3,6],
            [4,1]])
print(a*b)
# print(np.dot(a,b))
print(np.dot(mul,ni))
# print(np.dot(mul,np.linalg.inv(mul)))
# #print(np.dot(a,np.linalg.inv(a)))
# print(np.dot(mul,np.linalg.inv(mul)))
--
FROM 120.242.253.*