当条件数很大并且矩阵是对称的和不定的时候，有没有ojAlgo求解器？

文章/答案/技术大牛

发布

1回答

java、solver、linear-equation、ojalgo、indefinite

我使用ojAlgo来解线性方程组。在一种情况下，我得到了一个RecoverableCondition异常。可能是因为矩阵是病态的，所以条件数大约是1e15。我使用ojAlgo来解决它，如下面的代码所示。它通常是有效的，但在这种情况下不行。 有没有其他的解算器可以用于对称不定(病态)矩阵？目前失败的大小是18x18，但之后可能需

浏览 23提问于2021-10-14得票数 1

回答已采纳

2回答

matlab如何求解大型、对称和稀疏线性系统

matlab、linear-algebra、solver、sparse-matrix

也就是说，当我对一个非常大的，对称的，稀疏的A做A\b时，matlab使用什么算法？

浏览 6提问于2011-11-11得票数 2

回答已采纳

1回答

对于3x3只对称的正定线性系统，Cholesky比Householder还要快吗？

linear-algebra、eigen、solver、symmetric

我试图解一个线性系统Ax=b，其中A是3x3对称正定的。线性系统有许多求解者(C++，通过本征)。我个人更喜欢：HouseholderQr().solve()和llt().solve()，ldlt().solve()。我知道当n很大的时候，基于Cholesky分解的求解器</

浏览 0提问于2013-12-20得票数 1

回答已采纳

1回答

C++本征线性系统解决，数值问题？

c++、matlab、eigen

我在Eigen中的线性求解计算中得到了意想不到的无效结果。我有一个向量x和一个矩阵P。在MATLAB表示法中，我这样做：(类似于x'*inv(P)*x，但没有直接反转) 这是二次形式，应该产生一个积极的结果。矩阵P是正定的，并且没有病态。在MATLAB中，结果是肯定的，尽管很大。在我现在看到的例子中，

浏览 11提问于2018-08-23得票数 1

回答已采纳

5回答

大矩阵反演方法

c++、algorithm、math、linear-algebra、matrix-inverse

嗨，我一直在研究矩阵求逆(线性代数)，我想用C++模板编程来实现这个算法，我发现有很多种方法，比如:Gauss消去法或LU分解法，我找到了函数LU_factorize (c++ boost library我想知道是否还有其他方法，从程序员或数学家的角度来看，哪种方法更好(优点/缺点)？如果没有其他更快的方法，那么boost库中是否已经有了(矩阵)反演函数?

浏览 8提问于2010-07-28得票数 4

回答已采纳

3回答

Boost::uBLAS对Eigen

c++、boost、eigen

仔细观察这两种情况并不能使我对两者有更清楚的区别： uBLAS为稠密、单位和稀疏向量、稠密、恒等式、三角形、带状、对称、hermitian和稀疏矩阵提供模板化的C++该库涉及向量和矩阵上常用的基本线性代数运算:约简(如不同范数)、向量和矩阵的加减法以及用标量乘积、内积和外积乘法、矩阵向量和<e

浏览 2提问于2016-05-14得票数 21

1回答

解线性系统的最佳方法(就结果的准确性而言)是什么？

c++11、linear-algebra

解线性系统的最佳方法(就结果的准确性而言)是什么？它是否使用迭代方法(如共轭梯度...)或分解方法(如SVD或LU)？提前谢谢你。

浏览 11提问于2018-08-02得票数 0

2回答

Mathematica a追加矩阵

wolfram-mathematica、append

我必须在很大的参数空间范围内求解一个非线性方程组。我使用的是FindRoot，它对初始起点很敏感，所以我必须手动完成，通过反复试验和绘图来完成，而不是将方程放在循环或表格中。所以我想要做的是创建一个数据库或一个矩阵，它有固定的列数，但是可变的行数，这样我就可以在我求解它们的时候，不断地向它追加新的结果。FindRoot[{f1(

浏览 0提问于2012-02-03得票数 2

1回答

或者-工具谷歌如何告诉求解器在达到某个结果时停止

python、or-tools

我有一个简单的3d分配问题。我为成本矩阵随机生成0到100之间的整数。我注意到，当问题的大小足够大时，比如40个任务。(40X40X40成本矩阵).The总成本为0。这是有道理的。随着矩阵的增长，我假设有更多的方法可以达到预期的结果，并且我不希望求解时间呈指数级增长。但id确实如此。我假设这是因为算法不会满足于零结果，直到它绝对确定它是最优的

浏览 25提问于2020-10-12得票数 1

回答已采纳

2回答

为什么这个多项式方程条件不好？

matlab、math、polynomial-math、polynomial-approximations

我有1x1024矩阵。所以我想估计一个多项式方程。Y= acquired data.A 1024 element vectorpolyfit(x,y,5)我有很多这样的数字。在原始曲线和多拟合曲线之间是否有匹配的选择？

浏览 0提问于2015-08-24得票数 1

回答已采纳

2回答

求解一个大型(n>1000000)线性方程组

matlab、large-data、inverse、finite-element-analysis

在以下问题上，我已达到极限：我有两个HPC集群可供我处理，一个有5.7TB的RAM，另一个有1.5TB (但更快<em

浏览 5提问于2020-02-13得票数 2

7回答

ublas与矩阵模板库(MTL4)

c++、math、linear-algebra

我正在用c++写一个双曲型偏微分方程的软件。几乎所有的符号都是向量和矩阵符号。最重要的是，我需要线性代数求解器。是的，向量和矩阵的大小可以有很大的变化(从1000到只能通过分布式内存计算才能解决的大小，例如。集群或类似的体系结构)。如果我生活在乌托邦，我已经有了线性求解器，它对集群、GPU和多核具有很好

浏览 0提问于2009-07-01得票数 14

回答已采纳

2回答

针对奇异矩阵求解()崩溃

c++、armadillo、singular、matrix-inverse、numerical-stability

本文利用armadillo及其解函数，计算了抛物线拟合的三个系数。0.0221 0.0002 0.0000 |+--------------------------------++------------+| -0.0065 |+------------+ 显然求解()并不能解决它，甚至Matlab也警告说有没有任何解决办法，或更强大的&#

浏览 4提问于2013-01-23得票数 3

6回答

Ax=b线性代数系统的C++内存有效解

c++、boost、linear-algebra、lapack、umfpack

我正在使用Boost UBlas的数值库绑定来解决一个简单的线性系统。下面的代码运行良好，只是它仅限于处理相对较小的'm‘的矩阵A(m X m)。实际上，我有一个更大的矩阵，维度为m= 10^6 (最大可达10^7)。是否有现有的C++方法来解决有效使用内存的Ax=b。

浏览 3提问于2009-08-07得票数 9

回答已采纳

2回答

我想用python创建一个nxn对称矩阵。假设是n=9，那么我想要的内容如下： array[[0,1,0,0,0,1,1,0,1],[1,0,1,1,1,0,0,0,0],[0,1,0,1,1,0,0,0,0]….].我知道如何做到这一点，首先在python中创建一个in零矩阵(np.zeros((9，9)，然后使用循环填充它的1和0。但我觉得这不是一种巨蟒式的方式。所以我正在寻找一种优化的方法，如果矩阵很大

浏览 1提问于2015-03-19得票数 0

3回答

求逆协方差矩阵的马氏距离

c、algorithm、matrix、linear-algebra

我正在写一个函数来计算两个向量之间的马氏距离。我知道这是使用方程a'*C^-1*b来实现的，其中a和b是向量，C是协方差矩阵。我的问题是，有没有一种有效的方法可以在不使用高斯-乔丹消元法的情况下找到矩阵的逆矩阵，或者没有办法绕过这一点？我正在寻找一种自己做这件事的方法，而不是使用任何预定义的函数。我知

浏览 2提问于2012-08-03得票数 10

回答已采纳

1回答

末端执行器定位的逆运动学？

opengl、inverse-kinematics、skeletal-animation

我试图实现逆运动学求解，但这一次，即使在末端执行器的方向。当末端执行器只需要这个位置时，我就成功了。我了解到，在这种情况下，你可以构造这样的雅可比矩阵，其中w_i是全局空间中的i_th旋转轴，p_i是从i_th轴到目标位置的向量。问题是，当我必须在下面的方程中计算x_dot时。当x_dot只有需要考虑的位置，并且</em

浏览 2提问于2018-10-28得票数 3

回答已采纳

6回答

只有在找不到文件时才重定向请求？

apache、.htaccess、mod-rewrite

我希望有一种方法可以用mod_rewrite和Apache来实现，但是也许还有另外一种方法可以考虑。例如，假设客户机不需要自己的index.html。请记住，我不想在任何时候将客户端重定向到web根，如果客户端目录没有指定自己的</e

浏览 4提问于2009-02-27得票数 21

回答已采纳

1回答

绘制像鼠标光标一样始终位于顶部的移动图标，并适用于全屏应用程序

c++、windows、winapi、mouse

我需要绘制一个图标，移动大约与鼠标光标同步，并始终在所有窗口的顶部。操作系统: Windows 7 我有一个解决方案，通过在顶部最透明的窗口中绘制我的图标来在某种程度上工作。这种解决方案有一些主要的缺点，因为最上面的窗口会干扰其他最上面的窗口，并且一些全屏应用程序不能正常工作。例如，开始菜单和任务栏将覆盖我的窗口，如果我不定期将其设置为最顶部。对于一些全屏应用程序，更新图标窗口位置的性能下降很大，并且</em

浏览 0提问于2012-08-23得票数 1

3回答

Mat.inv()生成opencv中的所有零

c++、opencv

format(temp1, "numpy" ) << std::endl;转置正确工作，矩阵乘法也是正确的因此，Mat temp1的大小为960x960。然而，当我做temp2 =temp1.inv()时，我会在temp2中接收所有的零。我的意思是零是所有960x960的细胞。我不

浏览 4提问于2013-02-24得票数 7

回答已采纳

点击加载更多