我想用Python检查矩阵是正定的还是正半定的.
我怎样才能做到这一点? SciPy中是否有针对该模块或其他模块的专用功能?
解决方法:
我假设你已经知道你的矩阵是对称的.
对正定性的一个很好的测试(实际上是标准的!)是试图计算其Cholesky分解.如果你的矩阵是正定的,它就会成功.
这是最直接的方式,因为它需要O(n ^ 3)运算(具有小常数),并且您需要至少n个矩阵向量乘法来“直接”测试.
标签:python,scipy,matrix,math,linear-algebra
来源: https://codeday.me/bug/20191005/1854681.html