habe ein kleines Problem:
ich will in meinem Programm ein lineares Gleichungssystem (3x3) lösen, das sollte ja eigentlich kein großes Problem sein, soweit ich das sehe steht dazu ja auch numpy.linalg.solve bereit. Nur spuckt das mir leider falsche Ergebnisse heraus und ich kann mir nicht so wirklich erklären warum, zumal Beispiele mit ganzzahligen Komponenten richtig gerechnet werden.
Hier mal ein Minimalbeispiel:
Code: Alles auswählen
import numpy as np
from numpy import *
J_res = array([[ 0. ,-0.45464871, -0.45464871],
[-0.46242567, 0.24564775, -0.30389665],
[ 0.83722241, -0.15772861, 0.83722241]])
res = array([[ 1.29192658],
[ 0.83722241],
[ 0.46242567]])
x = np.linalg.solve(J_res, res)
print x
print (np.dot(J_res, x) == res).all()
[[ -3.02783731e+08]
[ -2.54783725e+08]
[ 2.54783723e+08]]
False
Grüße
Takuan