Das deutsche Python-Forum

Code: Alles auswählen

import scipy as S
import scipy.linalg as L
import random

def test(m):
	eig = L.eigvals(m)
	dim = S.shape(m)[0]
	eye = S.matrix(S.eye(dim,dim))
	for x in eig:
		eye_t = eye * x.real
		D = m - eye_t
		det = L.det(D)
		print 'Eigenwert:   ',x.real
		print 'Determinante nach Test:',det

def generate_matrix(dim, size):
	A = S.matrix(S.zeros((dim,dim)))
	for i in range(dim):
		for j in range(dim):
			A[i,j] = random.random()*size
	A *= A.transpose() #Symmetrisierung der Matrix
	print 'Groeßter Wert in der getesteten Matrix:',A.max()
	return A
	
def exe(i, size):
	for k in range(i):
		print 'Dimension der Matrix:',k+1
		A = generate_matrix(k+1,size)
		test(A)
		print'---------------'