资源描述：
《第2章非线形方程及其非线性方程组解法ppt课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、第二章方程(组)的迭代解法§1引言§2迭代解法§3迭代公式的改进§4联立方程组的迭代解法§5联立方程组的延拓解法§6联立方程组的牛顿解法1求f(x)=0的根§1引言1.1涉及到的概念f(x)既可以是代数多项式,也可以是超越函数f(x)=a0xn+a1xn-1+…+an-1x+an(a0≠0)如三角函数,指数函数的复合函数等方程的根:满足f(x)=0的x重根和单根:如果f(x)=(x-)mg(x)且g()≠0,则称为f(x)=0的m重根.m=1称为单根,m>1称为重根.2§1引言1.2本章重点介绍求方程实根的迭代解法(适用于求解代数方程和超越方程)代
2、数方程:根的个数与其最高次数相同,有成熟的圈定根的方法超越方程:可能有一个,几个根或者无解,无固定的圈定根的方法3§2迭代解法1本节重点（关键问题）根的初值的确定方法；迭代法的求解过程迭代法的收敛性迭代序列的误差估计4§2迭代解法2.1根的初值确定方法求方程根的几何意义:求曲线y=f(x)与x轴交点的横坐标。求根的具体步骤为:确定根的初值x0将x0进一步精确到所需要的精度52.1根的初值确定方法定理2.1设f(x)为区间[a,b]上的单值连续函数如果:f(a)·f(b)<0(2.1)则:[a,b]中至少有一个实根。如果:f(x)在[a,b]上还是单调地递
3、增或递减，则:仅有一个实根(有单根的条件)ba62.1根的初值确定方法2.1.1画图法画出y=f(x)的略图,从而看出曲线与x轴交点的大致位置。也可将f(x)=0分解为1(x)=2(x)的形式,1(x)与2(x)两曲线交点的横坐标所在的子区间即为含根区间.1.具体步骤72.1根的初值确定方法2.1.1画图法2.实例已知:f(x)=x㏒x–1=0;求根的初值范围(1)可以改写为:㏒x=1/x(2)画出对数曲线y=㏒x,与双曲线y=1/x,它们交点的横坐标位于区间[2,3]内解:82.1根的初值确定方法2.1.1画图法023yx192.1根的初值确定
4、方法2.1.2扫描法1.原理对于给定的f(x)，设有根区间为[A,B],从x0=A出发，以步长h=(B-A)/n(n是正整数)，在[A,B]内取定节点：xi=x0＋ih(i=0，1，2，…，n)，从左至右检查f(xi)的符号，如发现xi与端点xi-1的函数值异号，则得到一个缩小的有根子区间［xi-1,xi］。关键是选取步长hh太小，资源耗费大；h过大，可能遗漏根。102.1根的初值确定方法2.1.3对分(二分)法1.具体步骤若f(r)·f(a)>0，取a=r；否则取b=r设[xk-1,xk]为含根子区间，初值对于根的误差要求为ε,令a=xk-1,b=xk
5、,计算出f(a),f(b)后,进行如下:取[a,b]的中点r=(a+b)/2,计算f(r)若b-a>ε,转向Begin;否则结束Begin:112.1根的初值确定方法2.1.3对分(二分)法abx1x2abx*122.1根的初值确定方法2.1.3对分(二分)法13误差分析：第1步产生的有误差第k步产生的xk有误差对于给定的精度,可估计二分法所需的步数k：①简单;②对f(x)要求不高(只要连续即可).①无法求复根及偶重根②收敛慢注：用二分法求根，最好先给出f(x)草图以确定根的大概位置。或用搜索程序，将[a,b]分为若干小区间，对每一个满足f(ak)·f
6、(bk)<0的区间调用二分法程序，可找出区间[a,b]内的多个根，且不必要求f(a)·f(b)<0。14f(x)=0x=φ(x)等价变换思路:2.2迭代法的求解过程f(x)的根φ(x)的不动点由公式f(x)=0出发将其分解为等价形式x=φ(x)2.2.1建立迭代公式迭代函数例如：f(x)=x3+2x2-4=0可以分解为:x=x+f(x)=x3+2x2+x-4;x=2(1/(2+x))1/2x=x-(x3+2x2-4)/(3x2-4x)152.2迭代法的求解过程2.2.2迭代解法(简单迭代法)由初值x0出发,按迭代函数xn+1=φ(xn)(n=0,1,2,
7、…)进行计算迭代公式x0,x1,x2,…,xn,称为迭代序列；迭代序列的值相应地称为根的0次,1次,2次,…,n次近似值；序列的计算过程称为迭代过程；如果序列x0,x1,x2,…收敛于,即则为方程的根.证明：162.2迭代法的求解过程例2.2：用迭代法求方程f(x)=x3-x-1=0在x=1.5附近的根，要求根的近似值稳定至小数点后5位.解:(1)将方程改写为x=(1+x)1/3(2)按上式建立迭代公式xn+1=(1+xn)1/3x0=1.5(3)取x0=1.5逐次迭代得:x1=1.35721,x2=1.33086,x3=1.32588,x4=1.324
8、94,x5=1.32476,x6=1.32473,x7=1.32472,x8=1
