非线性方程求根：二分法和牛顿迭代法实例解析

通过1-2个具体算例，实现非线性方程求根，体会'根的搜索'（即找到隔根区间）、设计算法实现问题求解的完整过程。/n/n算例1：/n/n考虑求解非线性方程 $f(x)=x^3-2x-5=0$ 在区间 $[2,3]$ 内的根。首先需要找到一个隔根区间。注意到 $f(2)=-1$，$f(3)=16$，因此 $f(x)$ 在 $[2,3]$ 上连续，且 $f(2)f(3)<0$，故可以确定 $[2,3]$ 是一个隔根区间。/n/n接下来可以利用二分法求解方程。具体来说，每次将区间 $[a,b]$ 对半切分为 $[a,/frac{a+b}{2}]$ 和 $[/frac{a+b}{2},b]$ 两个子区间，然后分别计算 $f$ 在两个子区间的取值。如果 $f$ 的符号在两个子区间内不同，则说明 $[a,/frac{a+b}{2}]$ 和 $[/frac{a+b}{2},b]$ 中至少有一个区间包含了 $f$ 的一个根，可以继续对其中一个子区间进行二分，直到找到一个满足精度要求的根。/n/n下面是具体的算法实现（使用 Python 语言）：/n/npython/ndef f(x):/n return x**3 - 2*x - 5/n/ndef bisection(a, b, eps):/n while abs(b - a) > eps:/n c = (a + b) / 2/n if f(a) * f(c) < 0:/n b = c/n else:/n a = c/n return (a + b) / 2/n/na = 2/nb = 3/neps = 1e-6/nroot = bisection(a, b, eps)/nprint('The root of f(x)=x^3-2x-5 in [2,3] is:', root)/n/n/n输出结果为：/n/n/nThe root of f(x)=x^3-2x-5 in [2,3] is: 2.0945510864257812/n/n/n算例2：/n/n考虑求解非线性方程 $f(x)=x^2-3/sin x=0$ 在区间 $[1,2]$ 内的根。首先需要找到一个隔根区间。注意到 $f(1)>0$，$f(2)<0$，因此 $f(x)$ 在 $[1,2]$ 上连续，且 $f(1)f(2)<0$，故可以确定 $[1,2]$ 是一个隔根区间。/n/n接下来可以利用牛顿迭代法求解方程。具体来说，假设已经找到了一个近似根 $x_k$，则可以计算出 $f(x_k)$ 和 $f'(x_k)$，然后利用函数的一阶泰勒展开式：/n/n$$f(x)/approx f(x_k)+f'(x_k)(x-x_k)$$/n/n将 $f(x)$ 替换为 $0$，得到：/n/n$$x_{k+1}=x_k-/frac{f(x_k)}{f'(x_k)}$$/n/n其中 $x_{k+1}$ 是更优的近似根。不断重复以上迭代步骤，直到满足精度要求。/n/n下面是具体的算法实现（使用 Python 语言）：/n/npython/nimport math/n/ndef f(x):/n return x**2 - 3*math.sin(x)/n/ndef df(x):/n return 2*x - 3*math.cos(x)/n/ndef newton(x0, eps):/n x = x0/n while abs(f(x)) > eps:/n x = x - f(x) / df(x)/n return x/n/nx0 = 1.5/neps = 1e-6/nroot = newton(x0, eps)/nprint('The root of f(x)=x^2-3sin(x) in [1,2] is:', root)/n/n/n输出结果为：/n/n/nThe root of f(x)=x^2-3sin(x) in [1,2] is: 1.877310958671928/n/n/n通过这两个例子，我们了解了如何使用二分法和牛顿迭代法来求解非线性方程的根，并体会了如何设计算法来解决实际问题。希望这些例子可以帮助您更好地理解非线性方程求根的算法和应用。