2回答
在不使用numpy.linalg的情况下在Python中实现Choleky分解。*(是的，这是一个赋值)我遇到了一个关于使用浮点数的问题。non-positive definite")我事先测试了矩阵的对称性；n是维数，L是下三角Cholesky使用随机浮点n×n矩阵乘以它们的转置(为了获得正定矩阵)，两个ValueErrors都被提出，即w/out提高ValueErrors。L
浏览 34提问于2019-12-01得票数 3
回答已采纳
我正在将一个脚本从Matlab转换到Python，在Python中我有一个矩阵VarCov，并且我得到了它的Cholesky分解。有时，由于浮点近似，一个应该是正数(PD)的矩阵不是，我必须在对角线上加上一个小数字。VarCov + 0.0001 * eye(size(VarCov,1));end
Matlab很方便，因为它可以在进行Cholesky分解时返回矩阵是否为PD。这样做时，它不会引发错误。似乎在P
浏览 9提问于2018-03-05得票数 1
回答已采纳
1回答
我想知道每个表单的浮点操作的数量。维基百科引用了的一篇论文
n*(n+1)*(n+2)/6.
这就是我认为报纸得到n^3/6的地方。这是2*n个<em
浏览 1提问于2014-04-16得票数 3
回答已采纳
如果矩阵不是对称的正(半)定的，而Hessians是对称的，Cholesky分解就失败了。因此，我可以在我的矩阵上使用numpy.linalg.cholesky，这似乎是最有效的检查方法之一。见。现在，我不知道如何在没有错误的情况下执行Cholesky分解。例如
if np.linalg.cholesky(H) ...???
浏览 0提问于2019-07-07得票数 0
回答已采纳
6回答
高斯消去
这样的系统将被解决数百万次，因此算法在维数(n=9)方面应该是相当快的。只是提出一个错误，还是试图返回一些明智的答案？
提高误差是可以的。
浏览 4提问于2012-11-13得票数 7
回答已采纳
我在numpy.linalg.cholesky中看到了cholesky分解，但找不到LDU分解。有谁能推荐一个函数来使用吗？
浏览 0提问于2017-08-02得票数 10
回答已采纳
5回答
经过几次计算，我得到了一个浮点数矩阵，其中包含一些非常小的负项(例如-2.6e-9)，它们接近于0。我知道这个矩阵应该是正定的，即所有的正特征值，然而，我怀疑这些大约为0但负的项使得矩阵不是正定的，因此我不能做我想做的进一步的计算(Cholesky分解，对于那些想知道的人)。，我的问题是这个：我是否可以强制这些极小的条目
我想我误解了浮点数！
浏览 23提问于2021-12-06得票数 7
回答已采纳
2回答
在我的matlab代码中，我必须处理给定矩阵的Cholesky分解。我通常调用chol(A,'lower')来生成下三角因子。还需要注意的是，chol被称为4000时间，因为我需要4000不同配置的cholesky因子。
浏览 4提问于2013-02-26得票数 6
回答已采纳
1回答
我读到过在Numpy/Scipy中矩阵的Cholesky分解只有在它是正定的情况下才有效。实际上，下面的方法不起作用，因为矩阵是半正定的 np.linalg.cholesky([[1, 0], [0, 0]])array我
浏览 216提问于2019-03-11得票数 1
我的矩阵在做Cholesky分解时仍然是正定的，numpy给出的错误如下
浏览 49提问于2020-11-14得票数 0
我正在尝试做一个就地的Cholesky分解，但要么这个特性实际上没有在scipy中实现，要么就是有一些我不理解的东西。我在这里发帖，以防是后者。numpy.random.seed(0)V = numpy.dot(X.transpose(),X)scipy.linalg.cholesky然而，当我运行这段代码时，我的V和R是相同的(cholesky没有覆盖R)。我是误解了overwrite_a的用途，还是犯了其他错
浏览 0提问于2013-01-19得票数 4
回答已采纳
下一个可重复的代码生成183名正常人的394个观测值，并试图将它们与Cholesky分解联系起来：d <- 211
}v <- runif(183,0.6,0.8)生成cholesky矩阵相关的法线
浏览 3提问于2020-01-31得票数 0
9回答
我需要知道矩阵是不是。我的矩阵是numpy矩阵。我希望在numpy库中找到任何相关的方法，但没有成功。我很感谢你的帮助。
浏览 7提问于2013-04-29得票数 71
回答已采纳
我想用本征来计算稀疏矩阵的cholesky分解。然而，结果是不正确的，我找不到理由。我怎样才能得到正确的答案？decomposition cholesky.analyzePattern(XX);
std::cout << Eigen::MatrixXd(XX) &l
浏览 0提问于2018-11-09得票数 1
回答已采纳
1回答
我正在开发一些代码，可以从浮点或固定点从HW获取其数据。目前我们把它作为浮点。我认为使用浮点数会导致更大的舍入误差。
浏览 3提问于2016-11-24得票数 2
回答已采纳
我试图使用Julia中的isposdef()作为一种先验检验矩阵是否可以被cholesky分解分解的方法。
看起来isposdef并不总是有效的。我是不是用错了？
浏览 4提问于2014-06-19得票数 9
回答已采纳
然而，通常当我想使用Cholesky分解时，我这样做是为了解决一个线性系统，其中直接矩阵求逆可能是不稳定的。
因此，我正在寻找一种类似于在中实现的方法。
浏览 1提问于2015-11-14得票数 7
1回答
目前，我正在对不同的S使用我自己写的LDL分解函数，反演L，并计算变换，因为我正在处理少量不同的S和大n的情况。我使用Fortran，我宁愿使用一些LAPACK函数来分解因子，而不是我自己的LDL分解，我不能保证它是稳定的等等，但是我不确定哪种方法会是好的。
浏览 2提问于2012-12-13得票数 0
1回答
我试图建立一个超过183个变量的随机矩阵关联来计算Cholesky分解，并将183个随机法线关联起来。for (i in 1:183) { if (i < j) { }}我得到以下错误：
误差在chol.default(corr)中: 14阶的前导小
浏览 2提问于2020-01-31得票数 2
回答已采纳
云服务器
ICP备案
实时音视频
对象存储
即时通信 IM
腾讯云开发者
Copyright © 2013 - 2026 Tencent Cloud. All Rights Reserved. 腾讯云 版权所有
深圳市腾讯计算机系统有限公司 ICP备案/许可证号：粤B2-20090059 
粤公网安备44030502008569号
腾讯云计算（北京）有限责任公司 京ICP证150476号 | 京ICP备11018762号