[材料科学]2方程求根.pptVIP

第二章 方程求根 第一节 引言 求解f（x）=0，解x*称为方程的根，即f（x）的零点） 若f（x）=(x-x*)m g（x），g（x）?0，m?1则称x*为f（x）=0的m重根，或为f（x）的m重零点。 f（x）为超越函数时，就是超越方程，f（x）为n次多项式就是n 次代数方程。（n?5，不能用解析式表示解） ?线性的（一次解） ?单个方程? ?多项式（n个解） 代数方程? ?非线性? ? ?超越的（解的数目不定） ? ?线性（一组解） ? 方程组? ?非线性（多组解） 实际应用中只需求得满足一定精度的近似根即可。 方程求根分三步进行 1 .根的存在性 方程有没有根，有几个根？ 定理1（代数基本定理）：在复数范围内，n次代数方程至少有一个根。 f(x)=anxn+an-1xn-1+……+a1x+a0=0 n--正整数，x--复变量，a0—an为实、复常数。 定理2：n次代数方程有n个根。 2 .根的分布（有根区间） 求根的隔离区间，定一个[a, b]，使[a，b]内有且只有一个x*使f（x*）=0 定理3：设函数f（x）在[a，b]内连续，严格单调，且 f（a）*f（b）0，则在[a，b]内f（x）=0有且仅有一个实根。 通常有两种做法来确定隔离区间： 作y=f（x）的草图，看f（x）在x轴的交点位置来定区间[a，b] 逐步搜索，在连续区间[a，b]内，选取适当的x1，x2?（a，b），若f（x1）*f（x2）0，则[x1，x2]内有根。 例1：求xex-2 =0的有根区间 解：变形ex =2/x。则分别做f(x)= ex，g（x）=2/x的草图， f（x）与g（x）的交点的横坐标即是根的大致位置[0.5，1]内有唯一实根。 例2：求f（x）=x3-3x2+4x-3=0的有根区间。 解：三次多项式，作图复杂 f’（x）= 3x2-6x+4=3（x-1）2 +1 在（-∞，+∞）内f（x）=0，即f（x）为一单调递增函数。 那么f（x）=0在（-∞，+∞）内最多只有一个实根。 f（0）= 03-3*02+4*0-3=-30 f（1）= 13-3*12+4*1-3=-10 f（2）= 03-3*22+4*2-3=10 f（0）*f（2）0，故在[0，2]内有唯一实根。 根的精确化 已知根的近似值后，将其精确化，满足精度要求。 如例2中若取根x*=0.5*（0+2）=1，不精确 但若将[0，2]进一步缩小，则f（1.5）0，故[1.5，2]内有根，取x*=12（1.5+2）=1.75，这样精度就提高了。 第二节 二分法 一 方法 是一种直接法（一组非重复应用的公式），直观、简单。 设f（x）=0在[a，b]内连续，严格单调，且 f（a）*f（b）0 则f（x）=0在[a，b]内只有一个根，将区间对分，判别f（x）的符号，逐步缩小有根区间。 取xmid=0.5*（a+b） 若f（xmid）? (预先给定的精度)，则xmid即为根。 否则，若f（a）*f（xmid）0，则取a1=a，b1=xmid 若f（a）*f（xmid）0，则取a1=xmid，b1=b 此时有根区间缩小为[a1，b1]，区间长度为b1-a1=0.5*（b-a） 再对[a1，b1]进行上述过程，直到满足? 因为b1-a1=0.5*（b-a） b2-a2=(1/2)*（b1-a1）=(1/22)（b-a） …… bk-ak=(1/2k )*（b-a） f（ak）*f（bk）0 k??时， (1/2k)*（b-a）?0 取xk=(1/2)*（ak+bk）即最后一个区间的中点xk作为近似根，则误差为 ?xk-x*?? (1/2k+1)*（b-a）? (最后一个区间长度为(1/2k)*（b-a），故误差不超过长度的一半。） 例3：p17 二分法的优点：计算简单，方法可靠，只要求f（x）连续。 缺点：不能求偶数重根，也不能求复根，收敛速度与以1/2为比值的等比级数相同，不算太快。 因此，一般在求方程近似根时，不大单独使用，常用来为其它方法求方程近似根