数学分析6.7方程的近似解.docVIP

第六章微分中值定理及其应用

7方程的近似解

方程求解的方法主要有两种：解析法和数值法.

解析法又称为公式法，可以得到精确的解。

法国数学家伽罗瓦证明了形如：y=anxn+an-1xn-1+…+a0(an≠0)的代数方程中，当n≥5时，一般不存在求解公式.

数值解法——牛顿切线法

设f在[a,b]上二阶可导，满足：f’(x)·f”(x)≠0，f(a)·f(b)0，则

存在收敛点列{xn}，使其极限xn=ξ恰好是方程f(x)=0的解.因此，

当n充分大时，xn可作为ξ的近似值，即方程的近似解.

(1)设f’(x)0，f”(x)0，从而有f(a)0，f(b)0，并设f(ξ)=0.

∵f”(x)0，∴f在[a,b]上严格凸，∴f(x)f(a)+f’(a)(x-a),x∈[a,b]

设x0=a，则f在点a的切线与x轴的交点为x1=a-=x0-

∴f(x1)f(a)+f’(a)(x1-a)=0，如图(1).

以[x1,b]代替[a,b]重复上述步骤可得y=f(x)在点x1的切线与x轴交点为

x2=x1-，则f(x2)0，a=x0x1x2ξb.重复n次后可确定点列{xn}为

xn=xn-1-,n=1,2,….{xn}严格递增且有上界，可设xn=c.

由f和f’连续，可得：xn=(xn-1-)，即c=c-，∴f(c)=0.

由f严格单调，可知f(x)=0有唯一解：c=ξ.

又由中值定理有f(xn)=f(xn)-f(ξ)=f’(η)(xn-ξ)，xnηξ，

∴xn-ξ=，记m={|f’(x)|}，则|xn-ξ|≤.

同理可分析以下三种情形：(2)、(4)应取x0=b，得到单调递减的{xn}.

(2)若f’(x)0，f”(x)0，则有f(a)0，f(b)0；如图(2)

(3)若f’(x)0，f”(x)0，则有f(a)0，f(b)0；如图(3)

(4)若f’(x)0，f”(x)0，则有f(a)0，f(b)0.如图(4)

图(1)

图(2)

图(3)

图(4)

例：用牛顿切线法求方程x3-2x2-4x-7=0的近似解，使误差不超过0.01.

解：记f(x)=x3-2x2-4x-7，则f’(x)=3x2-4x-4=(3x+2)(x-2)；f”(x)=6x-4.

∴f’(-)=f’(2)=0；∵f”(-)=-80，f”(2)=80；

∴f有极大值f(-)0，极小值f(2)0.又f(x)=-∞，f(x)=+∞，

∴f(x)=0有且只有一个根.∵f(3)=-100，f(4)=90，

∴方程的根ξ∈(3,4).当x∈[3,4]时，f’(x)0，f”(x)0.

从点B(4,9)作切线与x由相交于x1=4-=4-≈3.68.

又m={|f’(x)|}=11，f(x1)=f(3.68)≈1.03，则

|x1-ξ|≤=0.01，不符合要求，再取点B’(x1,f(x1))作切线，得

x2=x1-=3.68-≈3.63，f(x2)≈-0.042，则

|x2-ξ|≤=0.01，∴取ξ≈3.63误差不超过0.01.

练习题

1、求方程x3-x2+2=0的实根到三位有效数字.

解：记f(x)=x3-x2+2，则f’(x)=x2-2x=x(x-2)；f”(x)=2x-2.

∴f’(0)=f’(2)=0；∵f”(0)=-20，f”(2)=20；

∴f有极大值f(0)=2，极小值f(2)=.又f(x)=-∞，f(x)=+∞，

∴f(x)=0有且只有一个根.∵f(-1)=0，f(-2)=-0，

∴方程的根ξ∈(-2,-1).当x∈[-2,-1]时，f’(x)0，f”(x)0.

从点A(-2,-)作切线与x由相交于x1=-2-=-2+=-≈-1.417.

x2=x1-≈-1.219，x3=x2-≈-1.196，x4=x3-≈-1.195，

∴取ξ≈-1.20.

2、求方程x=0.538sinx+1的根的近似值，精确到0.001.

解：记f(x)=x-0.538sinx-1，则f’(x)=1-0.538cosx0；∴f(x)在R上严格增.

又f(0)=-1，f(x)=+∞，∴方程有且只有一个根ξ.

∵f(1)=-0.538sin1=-0.415，f(2)=1-0.538sin2=0.51，∴方程的根ξ∈(1,2).

当x∈[1,2]时，f”(x)=0.538sinx0.

从点B(2,0.51)作切线与x由相交于x1=2-=2-≈1.583.

又m={|f’(x)|}=f’(1)≈0.707，f(x1)=0.583-0.538sin1.583≈0.044，

则|x1-ξ|≤==0.062，不符合要求，再取点B’(x1,f(x1))作切线，

得x2=x1-=1.583-≈1.538，f(x2)=0.538(1-sin1.538)≈0.000029，

则|x