南京工业大学 计算机在材料中应用 2.方程求根.pptVIP

2 方程求根; ?线性的(一次解) ?单个方程? ?多项式(n个解) 代数方程? ?非线性? ? ?超越的(解的数目不定) ? ?线性(一组解) ? 方程组? ?非线性(多组解) 实际应用中只需求得满足一定精度的近似根即可。;方程求根分三步进行;(2)根的分布(有根区间);通常有两种做法来确定隔离区间： 作y=f(x)的草图，看f(x)在x轴的交点位置来定区间[a, b] 逐步搜索，在连续区间[a, b]内，选取适当的 x1,x2?(a，b)，若f(x1)*f(x2)0，则[x1,x2]内有根。;例2.2 求f(x)= xex-2 =0的有根区间 ( P18) 解：变形ex =2/x, 则分别做f(x)= ex，g(x)=2/x的草图， f(x)与g(x)的交点横坐标即是根的大致位置[0.5,1]内有唯一实根。; 根的精确化 已知根的近似值后，将其精确化，满足精度要求。 如例2.1中若取根x*=0.5*(0+2)=1，不精确 但若将[0,2]进一步缩小，则f(1.5)0，故[1.5,2]内有根，取x*=12(1.5+2)=1.75，这样精度就提高了。;2.2 二分法; ;例2.3 用二分法求方程f(x)= x3-x2-2x+1=0在区间[0,1]内的1个实根，要求有3位有效数字 (P19) 二分法的优点：计算简单，方法可靠，只要求f(x)连续。 缺点：不能求偶数重根，也不能求复根，收敛速度与以1/2为比值的等比级数相同，不算太快。 因此，一般在求方程近似根时，不大单独使用，常用来为其它方法求方程近似根时提供好的初值。 对于方程求根，最常用方法是迭代法。;2.2.2 结构化流程图(N－S图);2.2.3 二分法的N－S图及程序设计;2.3 迭代法;考察方程：f(x)=x2-2x-3=0 解：∵ x2-2x-3=x2-2x+1-4=(x-1)2-22=0 ∴ (x-1)2 = 22 ∴ x-1=?2，即x1=3，x2=-1。 用迭代法求解： (1)变形f(x)成x=(2x+3)1/2 即?(x)= (2x+3)1/2 取x0=4，则xk+1=(2xk+3)1/2 (k=0，1，2……) x1=(2*4+3)1/2 =3.316， x2=(2*x1+3)1/2 = 3.104 x3=(2*x2+3)1/2 =3.034， x4=(2*x3+3)1/2 = 3.011， x5=(2*x4+3)1/2 = 3.004 当k?∞， xk ?3 (2)另变形成x=1/2*(x2-3) 即?(x)= 1/2 *(x2-3) 则xk= 1/2 *( xk 2-3) 取x0=4 ，x1= 1/2( x0 2-3)=6.5 ，x2= 1/2( x1 2-3)=19.625 x3= 1/2( x2 2-3)=191.070 当 k?∞， xk ? ∞ 两种方法中方法(1)是收敛的，而方法(2)发散。;2.3.1.2 收敛性; 通常迭代法的误差可以用前后两次迭代之值差的绝对值来估计误差，从而控制迭代次数。 即当 ?xk+1-xk?＜?时，迭代停止。 则取xk+1作为近似根。 这是一种事后误差估计法。;2.3.2 迭代法的几何解释;P1(x1，x2);例2.4 用迭代法求f(x)=x(x+1)2-1=0，在[0,1]内的一个实根。(补充) 解：变形x=1/(x+1)2 ，即?(x)=1/(x+1)2 ?‘(x)=-2/(x+1)3 。 ∵ ?‘(0)=-2 在[0,1]内??‘(x)?≮1，要缩小区间，用逐步搜索法 从x0=0开始，以h=0.2为步长 f(0.2)0，f(0.4)0, f(0.6)0 故[0.4,0.6]有根存在, 且??‘(x)?=? -2/(x+1)3?1 ∴ xk+1 =?(xk)=1/(xk+1)2 在[0.4,0.6]收敛 取x0=0.4，计算，当K=19时，有四位有效数字，取x*≈0.4656 则? x-x*?≤1/2*10-4 ? 逐步搜索以缩小有根区间