python怎么找一个矩阵_Python(NumPy,SciPy),找到一个矩阵的零空间

我试图找到给定矩阵的零空间(Ax = 0的解空间)。我发现了两个例子，但我似乎无法上班。此外，我不明白他们在做什么，所以我无法调试。我希望有人能够跟我走过。

文档页(numpy.linalg.svd和numpy.compress)对我来说是不透明的。我学会了通过创建矩阵C = [A | 0]，找到减少的行梯队形式和逐行求解变量来做到这一点。我似乎没有遵循这些例子中的做法。

感谢任何和所有的帮助！

A = matrix([

[2,3,5],

[-4,2,3]

])

import scipy

from scipy import linalg, matrixr

def null(A, eps=1e-15):

u, s, vh = scipy.linalg.svd(A)

null_mask = (s <= eps)

null_space = scipy.compress(null_mask, vh, axis=0)

return scipy.transpose(null_space)

当我尝试它，我得到一个空矩阵：

Python 2.6.6 (r266:84292, Sep 15 2010, 16:22:56)

[GCC 4.4.5] on linux2

Type "help", "copyright", "credits" or "license" for more information.

>>> import scipy

>>> from scipy import linalg, matrix

>>> def null(A, eps=1e-15):

... u, s, vh = scipy.linalg.svd(A)

... null_mask = (s <= eps)

... null_space = scipy.compress(null_mask, vh, axis=0)

... return scipy.transpose(null_space)

...

>>> A = matrix([

... [2,3,5],

... [-4,2,3]

... ])

>>>

>>> null(A)

array([], shape=(3, 0), dtype=float64)

>>>