2.12 求下列矩陣的特徵值和特徵向量的數值解和符號解 (-1 1 0 -4 3 0 1 0 2)的3*3矩陣

點選檢視程式碼

import numpy as np  
from scipy.linalg import eig  
  
# 定義矩陣  
A = np.array([[-1, 1, 0],  
              [-4, 3, 0],  
              [1, 0, 2]])  
  
# 計算特徵值和特徵向量  
eigenvalues, eigenvectors = eig(A)  
  
# 列印特徵值  
print("特徵值:")  
print(eigenvalues)  
  
# 列印特徵向量  
print("特徵向量:")  
for i in range(eigenvectors.shape[1]):  
    print(f"特徵值 {eigenvalues[i]:.2f} 對應的特徵向量:")  
    print(eigenvectors[:, i].real)  # 取實部，因為有時特徵向量會有複數部分，但在這個特定例子中它們是實數

print("學號：2023310143004")