万中数值分析综合复习课modified摘要.ppt

数值分析 综合复习;第一章 基本概念;1、误差的分类（绝对误差，相对误差）;2、有效数字;例1-2 下列近似数是通过四舍五入的方法得到的，试判定它们各有几位有效数字：;例1-3 已知 e =2.718281828……, 试判断下面两个近似数各有几位有效数字？; 定理1 设;如果;4、误差的传播;例1-4设有三个近似数;5、向量和矩阵范数;;例1-6 求矩阵;简述题： 1. 叙述在数值运算中，减少误差的原则是什么？ 解：减少误差的原则： 1）要避免除数绝对值远远小于被除数的绝对值 2）要避免相近数相减； 3）求和中要防止大数吃掉小数； 4）简化计算步骤，减少运算次数。;第二章 插值与最优逼近;Lagrange插值多项式;y 0 x;y 0 x;例2-1 给定数据表;例2-2 要制作三角函数sin x的值表.要求表值有四位小数。若用线性插值，则为保证截断误差不超过表值的舍入误差，其最大允许的步长是多少？;Newton插值多项式：;由Newton公式的递推式得到：;解：;例2-5;推广：等距节点(h),二次插值的误差界是;例2-6 ： 设f(x)＝x4，试利用拉格朗日插值余项定理写出以-1, 0, 1, 2为插值节点的三次插值多项式． ;例2-7：证明由下列插值条件 所确定的拉格朗日插值多项式是一个二次多项式． 该例说明了什么问题?; 容易验证 因而6个点(xi, yi)，i＝0 1,…,5均在二次曲线p(x)＝x2-1 上． 换句话说，满足所给插值条件的拉格朗日插值多项式为 p(x)＝x2-1.; 分析: 这是一个非标准插值问题,我们可以按各种思路去做．可按两种方法去做:一种是先求牛顿或拉格朗目型插值,再通过待定系数法求Pn(x)；另一种是先求埃尔米特插值,再通过待定系数法确定Pn(x) ．下面给出三种做法．;解法二 先作满足埃尔米特插值多项式 H3(x)．;易知;例2-10;例2-11;例2-12;例2-13;例2-14： 已知函数y=f(x) 的如下数据,试求其在区间[0,3]上的三次样条插值函数S(x)。 ;已知;这样便求得;3、最优平方逼近;考虑特殊情形-(1)用多项式｛1,x,x2,…,xn｝作n次最佳平方多项式p*(x)逼近步骤/方法(权函数为１时，［a,b］=[0,1]);(2)用正交多项式作最佳平方逼近;其中，;例2-16;最小二乘逼近步骤：;平方误差有与（2.4.15）相同形式的表达式。; (2)多项式的拟合; 例 2.19 用多项式拟合表2-4中的离散数据。;;;正交多项式拟合;例2-18： 对如下数据作形如 y = a eb x 的拟合曲线;对 z=A+Bx 作线性拟合曲线，取;得正则方程组;第三章数值积分;一、数值积分;引理：n阶Newton-Cotes公式的代数精度至少是n.;结论：当n 为奇数时，n阶Newton-Cotes公式的代数精度为n； 　　　当n 为偶数时，n阶Newton-Cotes公式的代数精度为n+1。;梯形公式;复化梯形公式;龙贝格(Romberg)求积算法; 稳定性定理 若求积公式(3.1.1)中系数Ak＞0 (k＝0，1，…，n)，则此求积公式是稳定的．;Gauss型求积公式;常用的Gauss求积公式;定理 3.4.2 设 　,则Guass公式(3.4.1)的余项是;解：;例3-2 试构造形如? f(x)dx ? A0f(0)+ A1f(h)+ A2f(2h) 的数值求积公式,使其代数精度尽可能高,并指出其代数精度的阶数.;而当f(x)=x3时, 公式的左边=81h4 /4, 右边=18h4, 公式的左边?右边,说明此公式对 f(x)=x3不能准确成立. 因此,公式只具有2次代数精度.;例3-3 ;例3-4 ;例3-5;例3-6;的近似值，要求误差;如此继续算得：;则三点Gauss-Legendre求积公式为：;例如，若要计算;例3-10;例3-11;第4-6章 综合复习课1、线性方程组数值解法2、非线性方程求根;一、解线性方程组的直接方法;解线性方程组的直接方法 一般的线性方程组解法： 列（全）主元素Gauss消元法 LU分解（直接三角分解法） 特殊的线性方程组解法