mikas blog

Die Gerschgorin-Kreise finde ich sehr elegant und via python.codefetch.com bin ich auf den Geschmack gekommen, die numerische und lineare Algebra mittels Python ein wenig nachzuvollziehen. Nur ein paar ganz simple Code-Snippets (unter Verwendung der interaktiven Python-Shell ipython) zur Anregung:

% ipython

from numarray import *
import numarray.linear_algebra as la

# Eigenwerte, Eigenvektoren, Inverse + Determinante:
a=array(((2,0,-1),(0,1,0),(-1,0,2)))
print a
print la.eigenvalues(a)
print la.eigenvectors(a)
print la.inverse(a)
print la.determinant(a)

# Pseudo-Inverse:
m=array(((0,1),(1,-1),(-1,-1)))
print m
print la.generalized_inverse(m,rcond=1e-10)

quit()

Oder alternativ auch via Scipy:

% ipython -p scipy

from scipy.linalg import *

a = array([[2,0,-1],[0,1,0],[-1,0,2]])

# Eigenwerte + Eigenverkten:
eig(a)

# Inverse:
inv(a)

# Determinante:
det(a)

quit()

This entry was posted by mika on Sunday, January 27th, 2008 at 11:47 and is filed under Computer. You can follow any responses to this entry through the RSS 2.0 feed. Both comments and pings are currently closed.