Implementing Eigenvectors & Eigenvalues in Python
Stryg for at vise menuen
Computing Eigenvalues and Eigenvectors
12345678910111213import numpy as np from numpy.linalg import eig # Define matrix A (square matrix) A = np.array([[2, 1], [1, 2]]) # Solve for eigenvalues and eigenvectors eigenvalues, eigenvectors = eig(A) # Print eigenvalues and eigenvectors print(f'Eigenvalues:\n{eigenvalues}') print(f'Eigenvectors:\n{eigenvectors}')
eig() from the numpy library computes the solutions to the equation:
eigenvalues: a list of scalars λ that scale eigenvectors;eigenvectors: columns representing v (directions that don't change under transformation).
Validating Each Pair (Key Step)
1234567891011121314151617import numpy as np from numpy.linalg import eig # Define matrix A (square matrix) A = np.array([[2, 1], [1, 2]]) # Solve for eigenvalues and eigenvectors eigenvalues, eigenvectors = eig(A) # Verify that A @ v = λ * v for each eigenpair for i in range(len(eigenvalues)): print(f'Pair {i + 1}:') λ = eigenvalues[i] v = eigenvectors[:, i].reshape(-1, 1) print(f'A * v:\n{A @ v}') print(f'lambda * v:\n{λ * v}')
This checks if:
Av=λvThe two sides should match closely, which confirms correctness. This is how we validate theoretical properties numerically.
Var alt klart?
Tak for dine kommentarer!
Sektion 1. Kapitel 39
Spørg AI
Spørg AI
Spørg om hvad som helst eller prøv et af de foreslåede spørgsmål for at starte vores chat
Sektion 1. Kapitel 39