实验名称非线性方程求根方法111汇编.doc

探索实验1 非线性方程求根方法 实验目的 理解非线性方程求根的过程和算法，学习用计算机求根的一些科学计算方法和简单的编程技术。 二、概念与结论 1. 零点定理：设f(x)在区间[a,b]连续,且f(a)f(b)0, 则至少存在一点x ((a,b(,使得f(x )=0。 2. 非线性方程: 函数方程f（x）= 0 中若f(x)不是x的线性函数, 则称f（x）= 0 为非线性方程, 特别若f(x)是n次多项式，则称f（x）= 0为n次多项式方程或代数方程；若f(x)是超越函数，则称f（x）= 0为超越方程。 3. 方程((x)=0 的根： 使方程((x)=0成立的x值称为方程((x)=0的根或函数((x)的零点，特别地，如果函数((x)可分解为 ((x)=（x((）mg(x) 且g(( )(0, 则称(是((x)的m重零点或((x)=0的m重根。 当m=1时，称(是((x)的单根 或单零点。 4. 简单迭代法收敛定理： 定理1：假定迭代函数((x) 满足下列条件： 1. 对任意 x∈[a,b]时,有((x) ∈[a,b] 2. 存在正常数 L1，使对任意x1,x2 ∈[a,b] 有 | ((x1) - ((x2 )|(L| x1 - x2 | 则((x)在[a,b]内有唯一的不动点x*，且对于任意初值 x0 ∈[a,b] 由迭代公式 xk+1= ((xk)产生的数列{xk}均收敛于方程根x*。 如果将条件2改为| (((x) |(L ， x ∈[a,b]，也有同样的收敛结果。 定理2：设x*是迭代函数((x)的不动点，m为正整数，且(（m）(x)在x*的邻域N(x*)内连续,并有如下关系 (（k）(x*) = 0 ,k=0,1,…m-1,(（m）(x*) (0 则由 xk+1= ((xk)产生的数列{xk}在邻域N(x*)内是m阶收敛的，且有极限 k(( 时 （ x* - xk+1 )/ （ x* - xk ) m ( (（m）(x*) /m! 5. 简单迭代法误差定理： 在简单迭代法收敛定理的条件下，有如下误差估计式 1) | x* - xk | ( | xk - xk-1 | L /(1-L) 2) | x* - xk | ( | x1 - x0 | Lk /(1-L) 6．Newton迭代法收敛定理： 假定函数f (x)在 x∈[a,b]时满足下列条件： 1. f(a)f(b)0 2.f((x ) (0, x∈[a,b] 3. f(( (x ) 存在且不变号, x∈[a,b] 则对于任意初值 x0 ∈[a,b] ，只要f(x0)f(( (x0 )(0,那么由Newton 迭代公式产生的数列{xk}一定收敛于方程根x*。 三、程序中Mathematica语句解释： f[x_]=Input[“键入函数f(x)=”] 从键盘上输入函数f(x) 变量= Input[ ] 从键盘上给变量赋值 3.N[精确数x] 或 精确数x //N　 将精确数x 转化成近似实数 4.N[精确数x, 正整数n] 将精确数x 转化成具有n位有效数字的近似实数 5.Print[表达式1,表达式2, …, 表达式n ] 在屏幕某一行上依次输出表达式1,表达式2,… 表达式n的值 6. D[f[x],x] 求函数f(x)对x的偏导数f((x)。 7. 如果条件成立,则执行对应的语句1，并将语句执行结果作为If语句的值，如果条件不成立，不执行语句1。8.If [条件, 语句1, 语句2] 根据条件的成立与否确定执行哪一个语句，具体执行为:条件成立时，执行语句1，否则，执行语句2，并将语句执行结果作为If语句的值。9. While[test,body] 当test为True时,计算body,重复对test的判断和body的计算,直到test不为True时终止 。test为条件, body为循环体,通常由body控制test值的变化。test不为True,则循环体不做任何工作. 四、方法与程序 非线性方程 f(x)=0求根的方法有区间法和迭代法两大类，二分法、弦位法就是区间法；简单迭代法和Newton迭代法及其变形是迭代法。这里只给出二分法、简单迭代法和Newton迭代法构造过程及程序。 1．二分法 二分法是最简单的求根方法，它是利用连续函数的零点定理，将含根区间逐此减半缩小，取区间的中点构造收敛点列｛xk｝来逼近