正交化矩阵numpy

时间:2012-12-18 19:26:28

标签: python numpy orthogonal

我尝试将2d-numpy array正交化并失败。我使用this方法并将其翻译为此代码:

def sym(w):
    return w.dot((w.T.dot(w))**(-.5))

In [1]: a
Out[2]: 
array([[ 1.1,  0.1],
       [ 0.1,  1.1]])
In [3]: a = sym(a)

In [4]: a
Out[5]: 
array([[ 1.20909392,  2.43574363],
       [ 2.43574363,  1.20909392]])

In [6]: a.dot(a.T)
Out[7]: 
array([[ 7.39475513,  5.89008563],
       [ 5.89008563,  7.39475513]])

a.dot(a.T)应该输出身份。

0 个答案:

没有答案