计算方法_复习课.pptxVIP

计算方法ComputationalMethodSchoolofAutomation.HuazhongUniversityofScienceandTechnology华中科技大学自动化学院蔡超caichao@hust.edu.cn第0章　绪论误差的定义误差限有效数字第0章典型例题例1．已知数x=2.718281828...,取近似值x*=2.7182,那麽x具有几位有效数字解:故有四位有效数点评：考查有效数字的概念。解：点评:此题考查相对误差的传播。例3．sin1有2位有效数字的近似值0.84的相对误差限是.解法1：根据有效数字与相对误差限的关系解法2：相对误差限的概念点评:此题考查相对误差与·有效数字关系第二种按定义得到的结果更好第1章　插值方法1问题的提出2拉格朗日插值公式3插值余项4埃特金插值方法5牛顿插值公式6埃尔米特插值7分段插值法8样条函数9曲线拟合的最小二乘法知识结构图典型例题提示：利用插值的拉格朗日余项说明当被插值函数为x的k次方时，误差为零。例：求下列超定方程组的最小二乘解：第2章　数值积分1机械求积2牛顿-柯特斯公式3龙贝格算法4高斯求积公式5数值微分知识结构图解以为插值节点，做的插值多项式解得从所求方程组看出，该公式至少具有2次代数精确度。具有3次代数精确度。解因为两点Gauss型求积公式的代数精确度为3．例用复化Simpson公式计算积分的近似值，并估计误差。（n=5,共11个节点）解：n=5，h=(1-0)/n=0.2，节点列为则复化Simpson公式为截断误差估计：30.9456909kT2k00.920735510.939793320946145909460833094608310.9460831第3章常微分方程的差分方法1欧拉方法2改进的欧拉方法3龙格-库塔方法4亚当姆斯方法5收敛性与稳定性6方程组与高阶方程的情形7边值问题知识结构图例1用欧拉法求初值问题当h=0.02时在区间[0,0.10]上的数值解。具体计算结果如下表：nxnyny(xn)n=y(xn)-yn001.00001.0000010.020.98200.98250.000520.040.96500.96600.000530.060.94890.95030.001440.080.93360.93540.001850.100.91920.9230.0021例2．取h=0.1,用改进欧拉法预报－校正公式求初值问题在x=0.1,0.2处的近似值.计算过程保留3位小数.解欧拉预报－校正公式为h=0.1,x0=0，y0=1，x1=0.1，于是有h=0.1,x1=0.1，y1=1.227，x2=0.2，于是有第4章　方程求根的迭代法非线性方程数值解法基本概念　（单根、重根、收敛阶）求根方法二分法及其收敛性简单迭代法简单迭代法的加速Newton迭代法Newton迭代法的变形重根Newton迭代法Newton下山法割线法迭代格式收敛性定理误差估计迭代格式收敛性定理误差估计知识结构图由于x0时，f’(x)0，且f(x)0，根据定理知：kxk00.8810.8846920.8846830.88468解：由可得由于证明该迭代公式平方收敛。根据定理该迭代公式平方收敛。第6章　直接法解线性方程组1消去法2追赶法3平方根法4误差分析第三章典型例题结 束数值运算的误差估计相对误差与有效数字间的关系定理设x*的近似数x具有标准形式：①若x具有n位有效数字，则相对误差§4复化求积公式一、复化梯形公式二、复化辛普森公式牛顿插值公式数值分析简明教程例1：构造三次埃尔米特插值多项式，要求:条件函数函数值导数值x0x1x0x10(x)10001(x)01000(x)00101(x)0001解：采用基函数法，构造三次埃米特插值多项式：由可将它写成数值分析简明教程条件函数函数值导数值x0x1x0x10(x)10001(x)01000(x)00101(x)0001即可得满足条件的三次埃尔米特插值多项式为一、样条插值的概念注：三次样条与分段埃尔米特插值的根本区别在于S(x)自身光滑，不需要知道f的导数值（除了在2个端点可能需要）；而埃尔米特插值依赖于f在所有插值点的导数值。若假设S(xi)=mi,i=0,1,…,n,利用分段Hermite插值多项式,当x[xi-1,xi]时,有其中hi=xi-xi-1.为了确定S(x),只需确定mi,i=0,1,…,n.可利用S(xi-0)=S(xi+0)来求出mi.当x[xi-1,xi]时,二、样条插值函数的建立——三转角插值法于是有由连续性条件S(xi-0)=S(xi+0)可得两侧同除以3(i[xi-1,xi]+i[xi,xi+1])=