Vježbe i ponavljanje

from IPython.lib.display import YouTubeVideo
vid = YouTubeVideo("JDX8-i0Zgps ")
display(vid)

Pogledajmo sljedeći primjer $$A=\begin{bmatrix} 121734 & 169217& 176624& 166662\\ 169217& 235222& 245505& 231653\\ 176624 & 245505& 256423& 242029\\ 166662& 231653& 242029& 228474\end{bmatrix}\quad b=\begin{bmatrix}634237\\ 881597\\ 920581\\ 868818\end{bmatrix}$$

import numpy as np
A = np.array([[121734, 169217, 176624, 166662], [169217, 235222, 245505, 231653], [176624, 245505, 256423, 242029],
              [166662, 231653, 242029, 228474]])
b = np.array([634237,881597, 920581, 868818])
x = np.linalg.solve(A,b)
print("x=",x)

x= [1.00170481 0.99897034 0.99957186 1.00025394]

Ako uzmemo da je desna strana sada b=\begin{bmatrix}634237\\ 881597\\ 920580\\ 868818\end{bmatrix}

c = np.array([634237,881597, 920580, 868818])
x = np.linalg.solve(A,c)
print("x=",x)

x= [16316026.0332614  -9854426.48680124 -4097527.25568108  2430330.02081108]

eig=np.linalg.eigvalsh(A)
print("Eigenvalues=", eig)

Eigenvalues= [1.06252778e-08 5.05278814e-01 2.21551535e+02 8.41630943e+05]

Matrica je loše uvjetovana!

print(np.linalg.cond(A,2))

79288124894603.9

Još malo o linearnim sustavima i primjeri

vid = YouTubeVideo("DMLummq67g8")
display(vid)

import numpy as np
A = np.array([[-4000, 2000, 2000], [2000, 0.78125, 0], [2000, 0, 0]])
b = np.array([400,1.3816, 1.9273])
x = np.linalg.solve(A,b)
print("x=",x)
np.matmul(A,x)

x= [ 0.00096365 -0.698496    0.9004233 ]

array([400.    ,   1.3816,   1.9273])

D = np.array([[0.0005, 0, 0], [0, 1, 0], [0, 0, 1]])

B=np.matmul(D,A)
G=np.matmul(B,D)
print("G=", G)
print("Uvjetovanost matrice A je", np.linalg.cond(A,2))
print("Uvjetovanost matrice G je",np.linalg.cond(G,2))

G= [[-0.001    1.       1.     ]
 [ 1.       0.78125  0.     ]
 [ 1.       0.       0.     ]]
Uvjetovanost matrice A je 13990.621644372157
Uvjetovanost matrice G je 4.656945664193629