import numpy as np
from function import schmidt
from sympy import GramSchmidt, Matrix
def schmidt(ma):
MA = [Matrix(col) for col in ma.T]
gram = GramSchmidt(MA,orthonormal=False) # orthonormal=True表示需要做归一化处理
return np.asarray(gram)
a = np.array([[1,0,0],[2,1,4],[0,1,3]])
print('a: \n {}'.format(a))
b = schmidt(a)
print("b: \n {}".format(b))
使用Python实现施密特正交
于 2023-08-15 18:59:27 首次发布