matlab非线性方程求解要点.pdfVIP

非线性方程的解法 1 引 言 数学物理中的许多问题归结为解函数方程的问题，即， f (x) 0 (1.1) 这里， f (x) 可以是代数多项式，也可以是超越函数。若有数 x* 为方程 f (x ) 0 的根，或称函 数 f (x) 的零点。 设函数 f (x) 在 [ a, b] 内连续，且 f (a ) f (b) 0 。根据连续函数的性质知道， 方程 f ( x) 0 在 区间 [ a ,b] 内至少有一个实根；我们又知道，方程 f (x ) 0 的根，除了极少简单方程的根可以 用解析式表达外，一般方程的根很难用一个式子表达。即使能表示成解析式的，往往也很复 杂，不便计算。所以，具体求根时，一般先寻求根的某一个初始近似值，然后再将初始近似 值逐步加工成满足精度要求为止。 如何寻求根的初始值呢？简单述之，为了明确起见，不妨设 f (x) 在区间 [ a,b] 内有一个实 的单根，且 f ( a) 0, f (b) 0 。我们从左端出点 x0 a 出发，按某一预定的步长 h 一步一步地 向右跨，每跨一步进行一次根的 “搜索 ”，即检查每一步的起点 xk 和 xk 1 (即， xk h )的函数值 是否同号。若有： f (x ) * f (x h) 0 (1.2) k k 那么所求的根必在 (x ,x h) 内，这时可取 x 或 x h 作为根的初始近似值。 这种方法通常称 k k k k 为 “定步长搜索法 ”。另外，还是图解法、近似方程法和解析法。 2 迭代法 2.1 迭代法的一般概念 迭代法是数值计算中一类典型方法，不仅用于方程求根，而且用于方程组求解，矩阵求 特征值等方面。迭代法的基本思想是一种逐次逼近的方法。首先取一个精糙的近似值，然后 1 用同一个递推公式，反复校正这个初值，直到满足预先给定的精度要求为止。 对于迭代法，一般需要讨论的基本问题是：迭代法的构造、迭代序列的收敛性天收敛速 度以及误差估计。这里，主要看看解方程迭代式的构造。 对方程 (1.1) ，在区间 [ a, b] 内，可改写成为： x (x ) (2.1) 取 [ , ] x0 a b ，用递推公式： x (x ) , k 0,1,2, (2.2) k 1 k 可得到序列： x , x ,x , x , { x } (2.3) 0 1 2 k k k 0 ~ ~ 当 k 时，序列 { x } 有极限 x ，且 (x) 在 x 附近连续，则在式 (2.2)两边极限，得， k k 0 ~ ~ x ( x ) ~ 即， x 为方程 (2.1)的根。由于方式 (1.1)和方程 (2.1)等价，