Open CASCADE学习|非线性方程组

非线性方程组是一组包含非线性数学表达式的方程,即方程中含有未知数的非线性项。解这类方程组通常比解线性方程组更为复杂和困难。

非线性方程组在很多领域都有应用,例如物理学、工程学、经济学等。解决非线性方程组的方法有很多种,包括数值方法和解析方法。数值方法是通过迭代或搜索来找到近似解,而解析方法则是通过对方程进行变换或展开来找到精确解。

在处理非线性方程组时,需要注意一些问题,例如初始值的选择、解的唯一性和稳定性等。同时,也需要根据具体问题的特点选择合适的求解方法。

OpenCASCADE提供了非线性方程组的类math_FunctionSet,下面给出一个具体的例子来说明其的用法。

待求解方程组:

从几何上看其解就是圆心在原点,半径为2的圆与曲线的交点:

cpp 复制代码
#include <math_FunctionSet.hxx>
#include <math_FunctionSetWithDerivatives.hxx>
#include <math_FunctionSetRoot.hxx>
​
class MyFunctionSet : public math_FunctionSetWithDerivatives
{
public:
    virtual Standard_Integer NbVariables() const
{
        return 2;
    }
​
    virtual Standard_Integer NbEquations() const
{
        return 2;
    }
​
    virtual Standard_Boolean Value(const math_Vector& X, math_Vector& F)
{
        F(1) = X(1) * X(1) + X(2) * X(2) - 4.0;
        F(2) = Pow(M_E, X(1)) + X(2) - 1.0;
        return Standard_True;
    }
​
    virtual Standard_Boolean Derivatives(const math_Vector& X, math_Matrix& D)
{
        // matrix D is Jacobi matrix.
        D(1, 1) = 2.0 * X(1);
        D(1, 2) = 2.0 * X(2);
        D(2, 1) = Pow(M_E, X(1));
        D(2, 2) = 1.0;
        return Standard_True;
    }
​
    virtual Standard_Boolean Values(const math_Vector& X, math_Vector& F, math_Matrix& D)
{
        Value(X, F);
        Derivatives(X, D);
        return Standard_True;
    }
​
private:
};
​
void test()
{
    MyFunctionSet aFunctionSet;
    math_FunctionSetRoot aSolver(aFunctionSet);
    math_Vector aStartingPoint(1, 2);
    // 1. (1.0, 1.0)
    aStartingPoint(1) = 1.0;
    aStartingPoint(2) = 1.0;
    aSolver.Perform(aFunctionSet, aStartingPoint);
    if (aSolver.IsDone())
    {
        aSolver.Dump(std::cout);
    }
​
    // 2. (1.0, -1.0)
    aStartingPoint(1) = 1.0;
    aStartingPoint(2) = -1.0;
    aSolver.Perform(aFunctionSet, aStartingPoint);
    if (aSolver.IsDone())
    {
        aSolver.Dump(std::cout);
    }
}
int main(int argc, char* argv[])
{
    test();
    return 0;
}
​

math_FunctionSetRoot Status = Done

Location value = math_Vector of Length = 2

math_Vector(1) = -1.81626

math_Vector(2) = 0.837368

Number of iterations = 14

math_FunctionSetRoot Status = Done

Location value = math_Vector of Length = 2

math_Vector(1) = 1.00417

math_Vector(2) = -1.72964

Number of iterations = 6

相关推荐
SheepMeMe几秒前
蓝桥杯2024省赛PythonB组——日期问题
python·算法·蓝桥杯
随便昵称几秒前
蓝桥杯专项复习——前缀和和差分
c++·算法·前缀和·蓝桥杯
commonbelive3 分钟前
团体程序设计天梯赛——L1-100 四项全能
c++
脑子慢且灵5 分钟前
蓝桥杯冲刺:一维前缀和
算法·leetcode·职场和发展·蓝桥杯·动态规划·一维前缀和
genispan7 分钟前
QT/C++ 多线程并发下载实践
开发语言·c++·qt
姜威鱼30 分钟前
蓝桥杯python编程每日刷题 day 21
数据结构·算法·蓝桥杯
CYRUS STUDIO38 分钟前
Unidbg Trace 反 OLLVM 控制流平坦化(fla)
android·汇编·算法·网络安全·逆向·ollvm
ゞ 正在缓冲99%…1 小时前
leetcode22.括号生成
java·算法·leetcode·回溯
小卡皮巴拉1 小时前
【力扣刷题实战】矩阵区域和
开发语言·c++·算法·leetcode·前缀和·矩阵
WG_171 小时前
第五章.图论
算法·图论