当前位置:高等教育资讯网  >  中国高校课件下载中心  >  大学文库  >  浏览文档

中国科技大学:《数值计算方法》课程教学资源(PPT课件讲稿)第4章 非线性方程求根

资源类别:文库,文档格式:PPT,文档页数:20,文件大小:492KB,团购合买
非线性科学是当今科学发展的一个重要研究方向,而非线性方程的求根也成了一 个不可缺的内容。但是,非线性方程的求根非常复杂。 通常非线性方程的根的情况非常复杂:
点击下载完整版文档(PPT)

中图苔技术大荸数学系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 第4章非线性方程求根 非线性科学是当今科学发展的一个重要研究方向,而非线性方程的求根也成了 个不可缺的内容。但是,非线性方程的求根非常复杂。 通常非线性方程的根的情况非常复杂: a=1无解 sin( x)=y x+a 个解 无穷组解 4 2 a 0两个解 a=-1四个解

数 学 系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 第4章 非线性方程求根 非线性科学是当今科学发展的一个重要研究方向,而非线性方程的求根也成了一 个不可缺的内容。但是,非线性方程的求根非常复杂。 通常非线性方程的根的情况非常复杂:      = = 2 1 ) 2 sin( y x y  无穷组解        = − = = =     = + = + 1 0 4 1 1 2 2 a a a a x y a y x a 无解 一个解 两个解 四个解

) 中图学技术大荸学系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 所以,只在某个区域内可能解存在唯一,而且经常很简单的形式得不到精确解: COS(a)=0 因此,通常我们用迭代法解非线性方程 看迭代法之前,先看看一种简单直观的方法 原理:f(a)f(b)<0→3x,st,f(x)=0

数 学 系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 所以,只在某个区域内可能解存在唯一,而且经常很简单的形式得不到精确解: 因此,通常我们用迭代法解非线性方程 看迭代法之前,先看看一种简单直观的方法 原理: f (a) f (b)  0  x,s.t., f (x) = 0 e −cos( x) = 0 x 

) 中图学技术大荸学系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 41对分法 什么时候停止? b 血2b x+-xk<a1或f(x)<a2

数 学 系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 4.1对分法 a b x1 x2 a b 什么时候停止? 1 1 x x ε k+ − k  2 或 f (x)  ε x*

中图苔技术大荸数学系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 算法 While(a-b>eps) X=(a+b)/2 f(×) 若(f(x)<eps)x为解 若f(x)*f(b)<0修正区间为[x,b] 若f(a)*(x)<0修正区间为[ax End while 每次缩小一倍的区间,收敛速度为Ⅰ2,较慢,且只能求一个根,使用条件限制较大 不能保证x的精度

数 学 系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS While(|a-b|>eps) x=(a+b)/2 f(x) 若(|f(x)|<eps) x为解 若f(x)*f(b)<0 修正区间为[x,b] 若f(a)*f(x)<0 修正区间为[a,x] End while 每次缩小一倍的区间,收敛速度为1/2,较慢,且只能求一个根,使用条件限制较大 算法  2 x* x 不能保证 x 的精度

) 中图学技术大荸学系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 42迭代法 等价变换 f∫(x)=0 x=g(x) ∫(x)的根 g(x)的不动点 从一个初值x出发,计算x1=g(x),x2=g(x1 思xk1=g(x),若{x}收敛,即存在x使得 路mx=x,且g连续,则由imx,=问知Xx) k→0 gx*),即x是g的不动点,也就是f的根

数 学 系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 4.2 迭代法 f (x) = 0 x = g (x) 等价变换 f (x) 的根 g (x) 的不动点 思 路 从一个初值 x0 出发,计算 x1 = g(x0 ), x2 = g(x1 ), …, xk+1 = g(xk ), … 若 收敛,即存在 x* 使得 ,且 g 连续,则由 可知 x* = g(x* ),即x* 是 g 的不动点,也就是f 的根。    k k=0 x lim x x * k k = → ( ) k k k k x g x → + → lim 1 = lim

) 中图学技术大荸学系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 迭代法的基本步骤如下 、给出方程的局部等价形式f(x)=0x=(x 2、取合适的初值,产生迭代序列x2x+1=(x1) 3、求极限 mxn易知,该值为方程的根 n→)+ ∧ 定收敛吗?

数 学 系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 迭代法的基本步骤如下: 1、给出方程的局部等价形式 f (x) = 0  x =(x) 2、取合适的初值,产生迭代序列 , ( ) 0 i 1 i x x = x + 3、求极限 n n x x →+ * = lim 易知,该值为方程的根 一定收敛吗?

) 中图学技术大荸学系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS y= y y=x y=glx) y=g(r)

数 学 系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS x y y = x x* y=g(x) x0 p0 x1 p1 ✓ x y y = x x* y=g(x) x0 p0 x1 p1 

中图苔技术大荸数学系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 定理 qp(x),x∈[a,b]若满足 1、a≤φ(x)≤b,x∈[a,b 2、0(x)可导,且存在正数L<1,使得对任意的x,有(x)≤L 则有 1、存在唯一的点x*,x*=(x*) 2、Vx∈[a,b]迭代收敛,且有误差估计

数 学 系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS (x), x[a,b] 若满足: 1、 a (x)  b, x[a,b] 2、 (x) 可导,且存在正数L<1,使得对任意的x,有 '(x)  L 则有: 1、存在唯一的点 x*, x* =(x*) 2、 x a,b  0  迭代收敛,且有误差估计 1 0 1 * x x L L x x k k − − −  定理

中图苔技术大荸数学系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS 证明 ①存在唯一性 做辅助函数v(x)=x-q(x),则有v(a)≤0,v(b)≥0 所以,存在点x*,S,t2y(x*)=0→x*=((x*) 若x*=(D(x**),则有: x*一x*料={(x*)一9(x*=9(5x*x*)≤Lx*一x*料 又,D<1→→x*=x** Vx∈[a,b]则 k+1-x*=p(xk)0(x*)=(5)xk-x*) k+-x8≤Dx-x米≤…≤+ Do-x 所以,任意的初值都收敛

数 学 系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS ①存在唯一性 ② 做辅助函数  (x) = x −(x) ,则有  (a)  0, (b)  0 所以,存在点 x*,s.t., (x*) = 0  x* =(x*) 若 x** =(x**) ,则有: x*−x** = (x*)−(x**) = '()(x*−x**)  L x*−x** 又, L 1  x*= x** [ , ] x0  a b 则 * ( ) ( *) '( )( *) 1 x x x x x x k+ − = k − =  k − * * 0 * 1 1 x x L x x L x x k k −  k −   − + +  所以,任意的初值都收敛 证明:

) 中图斜学技术大荸学系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS ③误差估计 1-x|=9(x)-0(x-1)≤儿x-xk1≤…≤Dx1 +…+x P x k+p k+p-1 k+1 ≤(Z+P++…+Z (1-LP k 1-L 中p的任意性,令p->+ x 证毕

数 学 系 University of Science and Technology of China DEPARTMENT OF MATHEMATICS ③误差估计 1 1 1 1 0 x x (x ) (x ) L x x L x x k k+ − k =  k − k−  k − k−  − k p k k p k p k k  x − x  x − x + + x − x + + + −1  +1 ( ) 1 0 1 L L x x k p k  + + − + +  ( ) 1 0 1 0 1 1 1 x x L L x x L L L k p k − − −  − − = 由p的任意性,令 1 0 1 * x x L L x x k k − − −  p →+ 证毕

点击下载完整版文档(PPT)VIP每日下载上限内不扣除下载券和下载次数;
按次数下载不扣除下载券;
24小时内重复下载只扣除一次;
顺序:VIP每日次数-->可用次数-->下载券;
共20页,试读已结束,阅读完整版请下载
相关文档

关于我们|帮助中心|下载说明|相关软件|意见反馈|联系我们

Copyright © 2008-现在 cucdc.com 高等教育资讯网 版权所有