施密特正交化

Python3

sympy.matrices

import numpy as np
from sympy.matrices import Matrix,GramSchmidt
 
A = np.array([[1,0,0],[0,2,1],[0,1,3]])
# 将numpy 矩阵转为sympy的列向量集合
MA = [Matrix(col) for col in A.T]
# 施密特正交化
gram = GramSchmidt(MA)

print(gram)