数值计算与MATLAB课件.pptVIP

电子信息工程学院 电子信息工程学院 第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