- 1、本文档共27页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
数值计算与MATLAB课件.ppt
电子信息工程学院 电子信息工程学院 第4章 求解非线性方程f(x)=0 电子信息工程学院 * 电子信息工程学院 * 数值计算与MATLAB 电子信息工程学院 * 第四章 求解非线性方程f(x)=0 在科学技术中,许多问题常归结为一元函数方程f(x)=0。方程按f(x)是多项式或超越函数分别称为代数方程或超越方程。 如果在区间[a,b]内只有方程f(x)=0的一个根,则称区间[a,b]为隔根区间。 描图法 逐步搜索法 电子信息工程学院 * 4.1 求解f(x)=0的MATLAB符号法 solve (s1, s2,…,sn, v1, v2,…, vn) solve (s1, s2,…,sn, v1,v2,…,vn) [z1,z2,…,zn]= solve (s1, s2,…,sn, v1, …, vn) 电子信息工程学院 * 例4-1 由方程 ax2+bx+5=0求出x和b来。 解:键入 s1=a*x^2+b*x+5; 或s1=a*x^2+b*x+5=0; 或s1=sym(a*x^2+b*x+5); 或s1=sym(a*x^2+b*x+5=0); 再键入x=solve(s1) 回车得出 若键入 b=solve(a*x^2+b*x+5,b) 或b=solve(s1, b) 回车得出 b= -(a*x^2+5)/x 电子信息工程学院 * 4.2 求方程f(x)=0数值解的基本方法: 二分法 迭代法 牛顿法 弦截法 电子信息工程学院 * 4.2.1 求实根的二分法原理 零点存在定理:设方程 f(x) = 0中的函数 f(x) 为实函数,且满足: ① 函数f(x)在[a,b]上单调、连续; ② f(a)、 f(b)异号 则(a,b)内至少存在一点ζ,使得f(ζ)=0。 电子信息工程学院 * 电子信息工程学院 * 例1 证明1-x-sinx=0在[0,1]内有根,使用二分法求误差不大于 的根需要二分多少次? clf,ezplot x-x, grid,hold,ezplot(1-x-sin(x)) 电子信息工程学院 * clear a=0;b=1; n=0; while b-a=0.00005 x=(a+b)/2; n=n+1; fa=1-a-sin(a); fx=1-x-sin(x); if fa*fx0 b=x; else a=x; end end x,n=n-1 电子信息工程学院 * x = 0.51095581054688 n= 14 电子信息工程学院 * 4.2.2 迭代法 迭代法基本原理:把方程等价地变换成形式 f(x)=x-g(x)=0 x=g(x) 若函数g(x)连续,则称上式为迭代函数。用它构造出迭代公式: xk+1=g(xk) , k=0,1,2,…… 从初始值x 0出发,便可得出迭代序列 { xk }=x0,x1,x2,…,xk,… 如果上述迭代序列收敛,且收敛于x*,则有: (g(xk)-xk+1)=(g(x*)-x*)=f(x*)=0, 可见x*便是方程 的根。 电子信息工程学院 * 电子信息工程学院 * 迭代公式收敛定理: 方程x=g(x)在(a,b)内有根x*,如果 (1)当x∈[a,b]时,g(x)∈[a,b]; (2)g(x)可导,且存在正数q1,使得对于任意x∈[a,b]都有 | |≤q1。 则: (1) 方程x=g(x)在(a, b)内有唯一的根x*; (2)迭代公式x k+1=g(x k)对(a, b)内任意初始近似根x0均收敛于x*; (3)近似根xk的误差估计公式为: | x*-x k|≤ |x1 -x0| 电子信息工程学院 * 例2 对下列方程构造收敛的迭代公式。 xex-1=0,(0,1); x3-x2-1=0, x0=1.5。 例3 编程用迭代法求1-x-sinx=0的根。 电子信息工程学院 * x1 = 0.5110 clear x0=0.5; x1=1-sin(x0); while abs(x1-x0)=0.00005 x0=x1; x1=1-sin(x0); end x1 电子信息工程学院 * 4.2.3 牛顿法(切线法) 切线法原理 :把函数f(x)在某一初始值x0点附近展开成泰勒级数,取其线性部分,近似地代替函数f
文档评论(0)