L利用Matlab进行数值分析().pptVIP

* * Lecture 5 : Numerical methods using Matlab(1) Roots of equation Previous Lectures: Basic matlab operation Matlab graphic Matlab programming Next lectures: Numerical methods using matlab Roots of equation Numerical methods using matlab: 求解方程的根。 数据插值和拟合 数值积分和微分(下载源码就到源码网:） Roots of equation 1. 二分法 2. 一般迭代法 3. 牛顿迭代法 本课程只解决区间 [a,b]之间只有一个根。 Roots of equation 1.二分法 定理：若函数 在区间[a, b]上连续, 且 , 则 在此区间内至少有一个根。 Matlab应用适合于[a,b]之间只有一个根。 Roots of equation 1.二分法 基本思想： 首先取区间[a,b]的中点c=(a+b)/2，保留有根的半个区间[a,c]或[c,b]；再取新的区间的中点，保留有根的半个区间，依次类推，直到区间长度减小到给定的精度error。此时，该区间内任意一点可以作为方程的根的近似值。 Roots of equation 1.二分法 保留有根区间算法 f(c)f(b)0 ，则保留[a,c] f(c)f(b)0 ，则保留[c,b] Roots of equation 1.二分法 算法设计 Input a, b, error Do repeat: c=(a+b)/2 if f(c)f(b)0 then a=c otherwise b=c Until (b-a)error End Roots of equation f=inline(x^2-2); a=1;b=2; error=10^(-2); while abs(b-a)error c=(a+b)/2; if f(c)==0 break elseif f(c)*f(b)0 a=c; else b=c; end end c bisec.m 取中点c，If f(c)=0则c就是根，无需后面的算法了。 求方程x2-2=0的正根 为什么选择a=1, b=2? Roots of equation f=inline(x^2-2); ezplot(f) 画图可以确定根的大致位置，从而确定a,b的值 Roots of equation 2. 一般迭代法 定义：不动点 函数g(x)在点xz邻域有定义且连续，如果xz=g(xz)，则点xz称为函数g(x)的固定点（fixed point) 显然，从几何上来看，点xz是函数y=x和y=g(x)的交点。 Roots of equation 2. 一般迭代法 定义：不动点迭代 Xn+1=g(xn)，n=0,1,2, ··· Roots of equation 2. 一般迭代法 不动点定理：函数g(x)在区间[a,b]连续，且满足： Roots of equation 2. 一般迭代法 求解f(x)=0的根的一般迭代法的基本思想： 将方程f(x)=0等价地转换为x=g(x)，转换要遵循不动点定理g(x)的满足条件。然后运用固定点迭代算法，在满足一定精度的条件下，迭代求出xz。 因为f(x)=0和x=g(x)等价，所以xz就是f(x)=0根。 Roots of equation 2. 一般迭代法 Roots of equation 2. 一般迭代法 关键在于由f(x)=0构建出等价的x=g(x) 原则上说，存在多种构建x=g(x)的方法，但迭代收敛的速度会不一样。 Roots of equation 算法设计： Begin input x, error Do repeat x=g(x) Until satisfy the condition End Roots of equation 2. 一般迭代法 求方程f(x)=x2-2=0的正根 根位于区间[1,1.5] Roots of equation 2. 一般迭代法 迭代收敛 Roots of equation 2. 一般迭代法 请完成matlab编程。 Roots of equation 2. 一般迭代法 g(x)的构建仍然是个问题，是否存在程式化的构建