中国传媒大学《MATLAB基础与应用》课件-第3-4章误差理论.pptVIP

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 的根。 逐次逼近： 将隐式方程归结为显式计算 4.2 简单迭代法 fixpt.m x y y = x x y y = x x y y = x x y y = x x* x* x* x* y=g(x) y=g(x) y=g(x) y=g(x) x0 p0 x1 p1 ? x0 p0 x1 p1 ? x0 p0 x1 p1 ? x0 p0 x1 p1 ? 原理：将非线性方程线性化 —— Taylor 展开 取 x0 ? x*，将 f (x)在 x0 做一阶Taylor展开: ，? 在 x0 和 x 之间. 将 (x* ? x0)2 看成高阶小量，则有： 线性 /* linear */ x y x* x0 只要 f ?C1，每一步迭代都有f ’( xk ) ? 0， 而且 ，则 x*就是 f 的根。 切线法 4.3 Newton法 newton.m 牛顿下山法 ——Newton’s Method 局部微调： 原理：若由 xk 得到的 xk+1 不能使 | f | 减小，则在 xk 和 xk+1 之间找一个更好的点 ，使得 。 xk xk+1 注：? = 1 时就是Newton’s Method 公式。 当 ? = 1 代入效果不好时，将 ? 减半计算。 newton_down.m 割线法 Newton’s Method 一步要计算 f 和 f ’，相当于2个函数值,比较费时。现用差商(f的值)近似f ’，可少算一个函数值。 x0 x1 切线 割线 切线斜率　?　割线斜率 需要2个初值 x0 和 x1。 收敛比Newton’s Method 慢，且对初值要求同样高。 secant.m 4.4 抛物线法 抛物线法是过曲线上的三点作一条 抛物线，用抛物线与x 轴的一个交点来 作为f (x)=0与x 轴交点。抛物线方法亦称 为Muller方法。抛物线法的迭代公式为： 其中， 。 parabola.m 4.5 非线性方程组的求解 非线性方程组可以看作非线性方程的推广，而非线性方程就是非线性方程组的特例。非线性方程组的一般数学描述为： 为叙述方便，记 这样上述方程组即可写为： 对于方程组的求解仍可以用牛顿法求解。 newtong.m 非线性方程的MATLAB函数求解 1、fzero()函数 MATLAB优化工具箱提供的fzero()函数是专门用于求解单变量非线性方程根的函数，该函数的调用格式为： [x,fval,exitflag,output]=fzero(fun,x0,options,p1,p2,…) 其中，fun表示函数表达式，x0是初始值，可以是标量或长度为2的向量，options是设置的过程参数，它主要包括Display和TolX两个选项，options选项可以通过函数optimset来设定，p1,p2,…是函数表达式中附加的参数，x是返回的根，fval是根x处的目标函数的值，exitflag表明解存在的情况，正数表明解存在，负数表示解不存在(遇到复数、NaN或者无穷大等)。Output包含计算过程中的信息，它是一个结构体，output.algorithm是所选用的算法，output.funcCount是函数赋值次数，output.iterations是迭代次数。 2、fsolve()函数 MATLAB最优化工具箱提供的fsolve()函数是专门用来求解多元方程的实数根的函数。它的调用格式如下： [x,fval,exitflag,output,jacobian]=fsolve(fun,x0,options,p1,p2,…) 其中，大多数参数的意义同fzero()函数，返回的jacobian是函数在x处的Jacobi矩阵， 第3章 误差理论 3.1 误差的来源 3.2 误差的基本概念 3.3 有效数字 3.4 误差的传播 3.5 数值计算中应注意的问题 3.6 MATLAB语言的数值计算精度 3.7 实例解析 本章目标：了解误差的来源、表示法、积累与传播等 中国传媒大学《MATLAB