有什么好的库可以解决C ++中的非线性方程组?

问题描述 投票:9回答:8

在我正在编码的C ++应用程序中,我需要求解一个非线性方程组(N个方程,N个未知数)。

我正在解决的系统相当小(最多10个方程/未知数),因此性能不会成为真正的问题。我在网上搜索了一些非线性解算器库,我无法找到一些看起来很容易使用的东西(得到了NOX和C/C++ Minpack,但两者似乎对我的需要来说太过分了)。

为此目的,易于使用的库的任何想法和想法?

c++ math libraries equation-solving
8个回答
2
投票

有两种选择,你可以使用包含非线性求解器的日with包,用C语言编写。我发现的唯一问题是你需要给它很好的初步估计。第二种选择是使用NLEQ或NLEQ2,我认为它是优越的(在FORTRAN中写作但很容易链接到C语言。但是我刚才找到了一些问题。有一个很好的网站,列出了可能的选项at:http://plato.asu.edu/sub/zero.html


3
投票

有一点应该清楚:非线性方程解决方案并不容易。它与求解线性方程式不同。您无法始终保证获得解决方案。您选择的初始条件和增量策略会对您获得的解决方案产生深远的影响。

话虽如此,我不能推荐一个特定的库,但你应该在其选择菜单中寻找包含Newton-Raphson迭代的线性代数包。


1
投票

Numerical Recipes有一个例行程序可以帮你完成工作。


1
投票

这取决于方程的非线性。如果它们具有一些“好”的特性......最明显的是正半确定矩阵或凸性,可能有专门的算法可用。我使用IBM / ILOG CPLEX来满足大多数线性编程需求。提供了可以引入C ++应用程序的库。虽然我没有使用它们的二次编程模块,但它确实是高马力线性和(良好行为)非线性编程的最先进技术。



0
投票

你看过COIN-OR吗?如果您将问题提交给OR-Exchange,可能会有所帮助。


0
投票

它无论如何都不是免费的,但Solver会在这里工作。


0
投票

微软Z3 https://github.com/Z3Prover/z3/blob/master/examples/c%2B%2B/example.cpp

还要考虑omnn :: math:https://github.com/ohhmm/openmind/blob/master/omnn/math/test/08_System.cpp

让我们说方程系统是这样的:

(x-a1)^2 + (y-b1)^2 = c1
(x-a2)^2 + (y-b2)^2 = c2

然后你有几个选择:

Valuable a1, a2, b1, b2; // init with values

System sys;
Variable x,y;
sys << (x-a1)^2 + (y-b1)^2 - c1; // addin an equation as an equality to 0
sys << (x-a2)^2 + (y-b2)^2 - c2;

for(auto& solution : sys.Solve(x))
        std::cout << solution;

另一种方法是制作单个方程(见why):

((x-a1)^ 2 +(y-b1)^ 2 - c1)^ 2 +((x-a2)^ 2 +(y-b2)^ 2 - c2)^ 2 = 0

Variable x,y;
Valuable a1, a2, b1, b2; // init with values
auto eq = ((x-a1)^2 + (y-b1)^2 - c1)^2 + ((x-a2)^2 + (y-b2)^2 - c2)^2;
eq.SetView(Valuable::View::Equation);  // optional: equation optimizations
// get y function:
auto fn = eq(y);

// show
std::cout << fn << std::endl;

// evaluate
auto evaluate = fn;
evaluate.eval(x, 10);
evaluate.optimize(); // calculate
// show calculated value at x=10:
std::cout << evaluate << std::endl;
© www.soinside.com 2019 - 2024. All rights reserved.