python牛顿法解非线性方程组_利用python求非线性方程

最近在做的东西中有一件任务，相当于一个函数已知y来求x，网上找了各种办法最终得以实现。在此说明方法，并记录一些坑。

要求的函数比如：*log(x) - log(1-x) + 2.2 * (1 -2x) *

最好用的方法，利用Scipy.optimize中的fsolve函数。

在该方法中，我们可以调用scipy.optimize.fsolve来求解非线性方程（组），具体方法如下：

from scipy.optimize import fsolve

import numpy as np

# 按格式要求定义我们需要求的函数

def f(x):

return np.log(x) -np.log(1-x) + 2.2*(1-2x)

# 调用fsolve函数

sol_fsolve = fsolve(f, [0.1, 0.9]) # 第一个参数为我们需要求解的方程，第二个参数为方程解的估计值

print(sol_fsolve)

[0.17071517 0.82928483]

# 输入两个解意味着根据你的估计值来进行梯度下降等算法找到的方程解。解的值域为（0，1），因此我估计为[0.1, 0.9]从而让函数从两边开始梯度下降，找到左右两个解

# 其实该方程有三个解，只是我不需要中间值的解，因此从左右两端估计

手动实现牛顿迭代法

牛顿迭代法是求非线性方程常用方法之一，具体原理如下：

随后，附上python实现代码：

from sympy import *

x = symbols('x')

f = log(x