第九讲--求代数方程的近似根(解).

《第九讲--求代数方程的近似根(解).》由会员分享，可在线阅读，更多相关《第九讲--求代数方程的近似根(解).（38页珍藏版）》请在文档大全上搜索。

1、求代数方程的近似根(解)数学实验q 问题背景和实验目的u 解方程（代数方程）是最常见的数学问题之一，也是众多应用领域中不可避免的问题之一。u 目前还没有一般的解析方法来求解非线性方程，但如果在任意给定的精度下，能够解出方程的近似解，则可以认为求解问题已基本解决，至少可以满足实际需要。u 本实验主要介绍一些有效的求解方程的数值方法：对分法，迭代法 和 牛顿法。同时要求大家学会如何利用Matlab 来求方程的近似解。相关概念相关概念0( )f x u 如果如果 f(x) 是一次多项式，称上面的方程为是一次多项式，称上面的方程为线性方线性方程程；否则称之为；否则称之为非线性方程非线性方程。q 线性方

2、程线性方程 与与 非线性方程非线性方程q 基本思想基本思想对分法对分法将有根区间进行对分，判断出解在某个分段内，然后再将有根区间进行对分，判断出解在某个分段内，然后再对该段对分，依次类推，直到满足给定的精度为止。对该段对分，依次类推，直到满足给定的精度为止。q 适用范围适用范围求有根区间内的求有根区间内的 单根单根 或或 奇重实根奇重实根。q 数学原理：数学原理：介值定理介值定理设设 f(x) 在在 a, b 上连续，且上连续，且 f(a) f(b)0，则由介值定，则由介值定理可得，在理可得，在 (a, b) 内至少存在一点内至少存在一点 使得使得 f( )=0。q 具体步骤具体步骤对分法对分

3、法设方程在区间设方程在区间 a,b 内连续，且内连续，且 f(a)f(b)0，给定，给定精度要求精度要求 ，若有，若有 |f(x)| ，则则 x 就是我们所需要就是我们所需要的的 f(x) 在区间在区间 (a,b) 内的内的 近似根近似根。；，计算令)( 2/ )( ) 1 (00 xfbax；输出结果停止计算，的近似根，就是我们所要，则若000 | )(| )2(xxxxf；否则令，令若bbxaxbaaxfaf1010110, ;, 0)()( )3(；输出结果，则停止计算，若令11111 | )(|, 2/ )( )4(xxxfbax；否则令，令若1212121211, ;, 0)()(

4、bbxaxbaaxfaf. .q 收敛性分析收敛性分析对分法收敛性对分法收敛性=11111 11|()()()22 22kkkkkkxbababa 设方程的根为设方程的根为 x* (ak , bk ) ，又，又 ，所以，所以2kkkabx 0(k )对分法总是收敛的对分法总是收敛的u 但对分法的收敛速度但对分法的收敛速度较慢较慢u 通常用来试探实根的通常用来试探实根的分布区间分布区间， 或给出根的一个较为或给出根的一个较为粗糙的近似粗糙的近似。根据上面的算法，我们可以得到一个每次缩小一半的根据上面的算法，我们可以得到一个每次缩小一半的区间序列区间序列 ak , bk ，在，在 (ak , bk

5、 ) 中含有方程的根。中含有方程的根。迭代法迭代法q 基本思想基本思想u 构造构造 f (x) = 0 的一个等价方程：的一个等价方程： ( )xx u 从某个近似根从某个近似根 x0 出发，计算出发，计算得到一个迭代序列得到一个迭代序列 0kkx 1()kkxx k = 0, 1, 2, . . (x) 的不动点的不动点f (x) = 0 x = (x)等价变换等价变换f (x) 的零点的零点u 若若 收敛，即收敛，即 ，假设，假设 (x) 连续，则连续，则q 收敛性分析收敛性分析迭代法的收敛性迭代法的收敛性 1limlim ()limkkkkkkxxx lim*kkxx *x( *)x k

6、x*( *)xx ( *)0f x 即即注：若得到的点列发散，则迭代法失效！注：若得到的点列发散，则迭代法失效！q 定义：定义：迭代法收敛性判断迭代法收敛性判断q 定理定理 2：如果定理如果定理 1 的条件成立，则有如下估计的条件成立，则有如下估计10|* |1kkqxxxxq 11|* |1kkkxxxxq 如果存在如果存在 x* 的的某个某个 邻域邻域 =(x*- , x* + ), 使使得对得对 x0 开始的迭代开始的迭代 xk+1 = (xk) 都收敛都收敛, 则称该迭代法在则称该迭代法在 x* 附近附近局部收敛局部收敛。q 定理定理 1：设设 x* = (x*)，的的某个某个 邻域邻

7、域 内连续，且对内连续，且对 x 都有都有 | (x)| q 1, 则对则对 x0 ，由由迭迭代代 xk+1 = (xk) 得到的点列都收敛。得到的点列都收敛。迭代法收敛性判断迭代法收敛性判断10|* |1kkqxxxxq q 定理定理 3：已知方程已知方程 x = (x)，且且(1) 对对 x a, b，有有 (x) a, b；(1) 对对 x a, b，有有| (x)| q p=2,-1,0,3; q=2,1; k=conv(p,q); q 多项式除法运算：多项式除法运算：k,r = deconv(p,q)其中其中 k 返回的是多项式返回的是多项式 p 除以除以 q 的商的商，r 是余式是

8、余式。k,r=deconv(p,q)p=conv(q,k)+r多项式的多项式的求求导导q polyderk=polyder(p) : 多项式多项式 p 的导数；的导数；k=polyder(p,q): p*q 的导数；的导数；k,d=polyder(p,q): p/q 的导数，的导数，k 是分子，是分子，d 是分母是分母 k1=polyder(2,-1,0,3); k2=polyder(2,-1,0,3,2,1); k2,d=polyder(2,-1,0,3,2,1);例：已知例：已知 ， ， 求求32)(23xxxp12)( xxq)/( ,)( , qpqpp多项式的多项式的值值q 计算计算

9、多项式多项式在给定点的值在给定点的值u 代数多项式代数多项式求值求值y = polyval(p,x): 计算多项式计算多项式 p 在在 x 点的值点的值注：若注：若 x 是向量或矩阵，则采用是向量或矩阵，则采用数组运算数组运算 (点运算点运算)！ p=2,-1,0,3; x=2; y=polyval(p,x) x=-1, 2;-2,1; y=polyval(p,x)例：已知例：已知 ，分别取，分别取 x=2 和一个和一个 2 2 矩阵，矩阵， 求求 p(x) 在在 x 处的值处的值32)(23xxxp多项式的多项式的值值u 矩阵多项式矩阵多项式求值求值Y=polyvalm(p,X)l 采用的是