河海大学生数值试题.ppt

数值分析 编者 李庆扬等 课件 董祖引 Ch 1 绪论 1.1 数值分析研究对象与特点 科学技术领域的三大环节： （1）实验； （2）科学计算； （3）理论。 所谓的科学计算是指使用计算机进行数值计算的工 作，它可以部分地替代科学实验。 科学计算的主要过程： 实际问题 数学模型 数值计算方法 程序设计 上机计算求出结果 其中数值计算方法是数值分析研究的对象。 主要包括： （1）函数的数值逼近（包括插值法）； （2）数值微分和数值积分； （3）非线性方程（组）数值解； （4）数值线性代数（如线性方程组数值解、矩阵特征值特征向量的计算）； （5）（偏）微分方程数值解。 数值分析的主要特点： （1）可行性； （2）时效性； （3）可靠性； （4）实验性。 1.2 数值计算的误差 误差的来源与分类： （1）模型误差； （2）观察误差； （3）截断误差（也称方法误差）；例如： 其截断误差： （4）舍入误差（包括由十进制转换为二进制引起 的误差）。例如： 其舍入误差： 绝对误差： 并记： 相对误差： 实际计算中通常用 相对误差限： 它可表示为： 例1 按四舍五入原则下列各数： 187.9325， 0 8.000033， 2.7182818， 2/3 具有5位有效数字的近似值分别为： 187.93， 0.037856， 8.0000， 2.7183， 0.66667 至少具有n位有效数字。 （证明见黑板） （见黑板） 一般的，利用Taylor展式，有 （见黑板） 1.3 误差定性分析与避免误差危害 一般的，误差的定量分析是困难的，我们考虑下 列三各方面的定性分析。 （1）病态问题与条件数。 若输入数据有微小扰动（即误差），引起输出数 据（即问题的解）相对误差很大，则称为病态问题。 病态程度可用相对误差比值来描述，如： 它表示相对误差大约放大n 倍。 (2)算法的数值稳定性。 一个算法如果输入数据有误差，而在计算过程中 舍入误差不增长，则称此算法是数值稳定的。 比如，page9例5的A算法是数值不稳定的，B算法是稳定的。 (3)避免误差危害的若干原则。 1.避免小数作除数。 2. 避免两相近数相减。 例4 利用中心差商公式计算 （详见黑板） 3. 防止大数“吃”小数。 4. 简化计算步骤，减少运算次数。 例5 计算多项式 的值。 直接计算共需要做 次乘法和 n 次加法。 若用秦九韶算法 只要n次乘法和n 次加法即可。 Ch2 插值法 2.1 引言 上的值 插值区间。 常用的有拉格朗日插值、牛顿插值、埃尔米特插 值、埃特金插值、三次样条插值等。 2.2 拉格朗日插值 定义1 若n次多项式 在n+1个节点 上满足 事实上，拉格朗日插值基函数 令 则满足： 特别的，当n=1， 称线性插值。 当n=2， 称抛物插值二次插值。 若引入记号 则拉格朗日插值多项式可写成如下的紧凑形式： x 0 1 3 2 0 -4 （见黑板） 证明 只要证明唯一性。 则 这与代数基本定理（ n 次多项式至多有 n 个零点）矛 盾。 特别的，若令 则得 （证明见黑板） 例2 已给 sin0.32=0.314567，sin0.34=0.333487， sin0.36=0.352274，试用抛物插值计算 sin0.3367，并估 计截断误差。 （见黑板） （详见黑板） 最后，关于插值多项式的数值计算的稳定性，我 们有结论： （1）线性插值数值计算是稳定的； （2）高次（非线性）插值数值计算是不稳定的。 （详见黑板） 2.3 均差与牛顿插值公式 牛顿插值法的思想是将插值多项式表写成 可由插值条件确定。 它们正好是下述定义的均差。 定义2 称 利用归纳法可以证明： 均差具有如下基本性质： （1）（对称性）均差与节点的排序无关，即 （2） （证明略） 具体计算时可列均差表（详见黑板）。 将后式代入前式，即得 其中 显然满足插值条件，称为牛顿（均差）插值多项式。 插值余项 注意：根据插值多项式的存在唯一性定