Python第三方库scipy.optimize中的fsolve和root这两个函数可用于非线性方程组的求解。但其更适用于结构形式稍显简单的方程组求解。对于多变量结构复杂的非线性方程组的求解操作起来稍显困难。因此,自己动手编写求解非线性方程组的程序。
本文采用了Newton-Raphson求解方法,其中利用了一阶差商代替偏导数的求解得到雅可比矩阵。
算法流程:
首先取初值 ,然后做以下迭代:
(1) 计算
Python第三方库scipy.optimize中的fsolve和root这两个函数可用于非线性方程组的求解。但其更适用于结构形式稍显简单的方程组求解。对于多变量结构复杂的非线性方程组的求解操作起来稍显困难。因此,自己动手编写求解非线性方程组的程序。
本文采用了Newton-Raphson求解方法,其中利用了一阶差商代替偏导数的求解得到雅可比矩阵。
首先取初值 ,然后做以下迭代:
(1) 计算