#système 2

import numpy as np
A = np.array([[4,-1,1],[4,-8,1],[-2,1,5]])
B=np.linalg.matrix_power(A,2)

B=np.linalg.matrix_power(A,3)
#transposition
print(np.transpose(A))

#multiplication
print(np.dot(A,A))

#déterminant
print(np.linalg.det(A)) # -154
#inversion
print(np.linalg.inv(A))

Y=np.array([[7],[-21],[15]])
print(A)
print(Y)
A_inv = np.linalg.inv(A)
print(A_inv)
Solution_x=np.linalg.solve(A,Y)
print(Solution_x)
#valeurs propres et vecteurs propres
print(np.linalg.eigh(A))
np.tril(A, k=0)
np.triu(A, k=0)
[[ 4 4 -2]
[-1 -8 1]
[ 1 1 5]]
[[ 10 5 8]
[-18 61 1]
[-14 -1 24]]

-154.00000000000006
[[ 0.26623377 -0.03896104 -0.04545455]
[ 0.14285714 -0.14285714 -0. ]
[ 0.07792208 0.01298701 0.18181818]]
[[ 4 -1 1]
[ 4 -8 1]
[-2 1 5]]
[[ 7]
[-21]
[ 15]]
[[ 0.26623377 -0.03896104 -0.04545455]
[ 0.14285714 -0.14285714 -0. ]
[ 0.07792208 0.01298701 0.18181818]]
[[2.]
[4.]
[3.]]
[[2.]
[4.]
[3.]]

#valeurs propres et vecteurs propres print(np.linalg.eigh(A))

#système 3
import numpy as np
A = np.array([[1,3,3],[2,2,5],[3,2,6]])
B=np.linalg.matrix_power(A,2)
B=np.linalg.matrix_power(A,3)
#transposition
print(np.transpose(A))
#multiplication
print(np.dot(A,A))
Y=np.array([[-2],[7],[12]])
print(A)
print(Y)
#déterminant
print(np.linalg.det(A)) # 5
#inversion
print(np.linalg.inv(A))
A_inv = np.linalg.inv(A)
print(A_inv)
Solution_x=np.linalg.solve(A,Y)
print(Solution_x)
#valeurs propres et vecteurs propres
print(np.linalg.eigh(A))
np.tril(A, k=0)
np.triu(A, k=0)
[[1 2 3]
[3 2 2]
[3 5 6]]
[[16 15 36]
[21 20 46]
[25 25 55]]
[[1 3 3]
[2 2 5]
[3 2 6]]
[[-2]
[ 7]
[12]]
5.000000000000002
[[ 0.4 -2.4 1.8]
[ 0.6 -0.6 0.2]
[-0.4 1.4 -0.8]]
[[ 0.4 -2.4 1.8]
[ 0.6 -0.6 0.2]
[-0.4 1.4 -0.8]]
[[ 4.]
[-3.]
[ 1.]]
(array([-0.8661408 , 1.35679908, 8.50934171]), array([[ 0.86736726,
0.25247328, 0.4288721 ],
[-0.42774337, 0.81868546, 0.38313147],
[-0.25438089, -0.51576289, 0.81809474]]))
array([[1, 3, 3],
[0, 2, 5],
[0, 0, 6]])

#système 4
import numpy as np
A = np.array([[-1,3,-4],[1,-5,4],[2,1,8]])
#déterminant
print(np.linalg.det(A)) # 0, la matrice n'est pas inversible
#inversion
print(np.linalg.inv(A))
#transposition
print(np.transpose(A))

#multiplication
print(np.dot(A,A))

#déterminant
print(np.linalg.det(A)) # 0

#inversion
print(np.linalg.inv(A))

Y=np.array([[2],[-1],[5]])
print(A)
print(Y)
A_inv = np.linalg.inv(A)
print(A_inv)
Solution_x=np.linalg.solve(A,Y)
print(Solution_x)

0.0

LinAlgError Traceback (most recent call

last)
<ipython-input-5-9944941bdbec> in <module>()
7
8 #inversion
----> 9 print(np.linalg.inv(A))
10
11 Y=np.array([[2],[-1],[5]])

< array_function internals> in inv(*args, **kwargs)

/usr/local/lib/python3.7/dist-packages/numpy/linalg/linalg.py in
inv(a)
543 signature = 'D->D' if isComplexType(t) else 'd->d'
544 extobj =
get_linalg_error_extobj(_raise_linalgerror_singular)
--> 545 ainv = _umath_linalg.inv(a, signature=signature,
extobj=extobj)
546 return wrap(ainv.astype(result_t, copy=False))
547

/usr/local/lib/python3.7/dist-packages/numpy/linalg/linalg.py in
_raise_linalgerror_singular(err, flag)
86
87 def _raise_linalgerror_singular(err, flag):
---> 88 raise LinAlgError("Singular matrix")
89
90 def _raise_linalgerror_nonposdef(err, flag):

LinAlgError: Singular matrix

#valeurs propres et vecteurs propres print(np.linalg.eigh(A)) np.tril(A, k=0) np.triu(A, k=0)

#système(S3)
Le système admet aucune solution car ce n'est pas une matrice carrée
import numpy as np
A = np.array([[4,2,-2],[3,-1,1],[1,1,1],[1,-1,1]])
#le système admet aucune solution car ce n'est pas une matrice carrée
B=np.linalg.matrix_power(A,2)
B=np.linalg.matrix_power(A,3)
#transposition
print(np.transpose(A))
#multiplication
print(np.dot(A,A))
Y=np.array([[0],[3],[1],[-2]])
print(A)
print(Y)
#déterminant
print(np.linalg.det(A))
#inversion
print(np.linalg.inv(A))
A_inv = np.linalg.inv(A)
print(A_inv)
Solution_x=np.linalg.solve(A,Y)
print(Solution_x)
#valeurs propres et vecteurs propres
print(np.linalg.eigh(A))
np.tril(A, k=0)
np.triu(A, k=0)

LinAlgError Traceback (most recent call

last)
<ipython-input-15-8f6254f1dc8f> in <module>()
5 #le système admet aucune solution car ce n'est pas une matrice
carrée
6
----> 7 B=np.linalg.matrix_power(A,2)
8
9 B=np.linalg.matrix_power(A,3)
< array_function internals> in matrix_power(*args, **kwargs)
/usr/local/lib/python3.7/dist-packages/numpy/linalg/linalg.py in
matrix_power(a, n)
618 a = asanyarray(a)
619 _assert_stacked_2d(a)
--> 620 _assert_stacked_square(a)
621
622 try:
/usr/local/lib/python3.7/dist-packages/numpy/linalg/linalg.py in
_assert_stacked_square(*arrays)
201 m, n = a.shape[-2:]
202 if m != n:
--> 203 raise LinAlgError('Last 2 dimensions of the array
must be square')
204
205 def _assert_finite(*arrays):
LinAlgError: Last 2 dimensions of the array must be squar