计算方法130098.pptVIP

* 《计算方法》第一讲 课程内容简介 误差概念与有效数字 递推算法的数值稳定性问题 MATLAB应用 数值计算对计算机发展的影响 世界上第一台电子计算机ENIAC研制 目的是为了解决大型数值计算问题 程序设计中的数学原理和背景 算法构造、算法分析、方法的近似性 计算方法课主要内容 数值逼近、数值代数、微分方程数值解 冯·诺依曼 人类的计算能力不断受到挑战 大规模科学和工程计算、数值模拟 计算n次多项式秦九韶算法 算法II 置Sn = an, 计算 Sk－1 = a k －1+ xSk (k = n，n－1，…，1 ) 算法I 置S 0 = a0 ; Sk = Sk－1 + ak xk ( k = 1，2，…，n ) 输入 x; a0,a1,…,an S ?an k从n到1循环 S?ak－1+x×S 输出数据S ;结束 误差概念 例1.5 立方倍积问题.作一个立方体,使其体积为已知立方体的二倍. V=1 V=2 V=1: h=1 V=2 : h= ≈1.25992104989487(15位) 例1.3 Fibonacci数列递推公式: Fn = Fn – 1 + Fn – 2 (F0 = 1，F1= 1) F10=89.00000000000003? 模型误差: 建立数学模型时所引起的误差； 误差来源： 舍入误差：计算机表示的数的位数有限，通常用四舍五入的办法取近似值，由此引起的误差. 截断误差：用数值法求解数学模型时,常用有限过程代替无限过程,得到的近似解和模型准确解之间的误差； 观测误差：测量工具的限制或在数据的获取时随机因素所引起的物理量的误差； 机器数系只是有限、离散的数集 实数系则是无限、 稠密、 连续的数集. 在科学计算中实数通常被表示成如下形式 ±0.a1a2…an ×10 m 其中, a1≠0.且 a1, a2, …, an 都是 0 ～ 9 中的任一整数, 0.a1a2…an 称为尾数. 用于确定小数点位置，指数 m 称为阶码. 例如,计算机用64个位表示一实数.按惯例10个位表示阶码,52个位表示尾数,另两个位表示正负号.此时,10进制浮点数具15位,阶码不超过308 假设某一数据的准确值为x*,其近似值为x,则称E(x)= x - x* 为 x 的绝对误差. 而称 为 x 的相对误差 绝对误差与相对误差 如果存在一个适当小的正数ε ，使得 则称ε为绝对误差限。 则称εr为相对误差限 如果存在一个适当小的正数εr ，使得 例 取 π 的有限位数如下： 取一位：x1 = 3，误差限不超过0.5; 取三位：x2 = 3.14，误差限不超过0.005 取五位：x3=3.1415, 误差限不超过0.00005 若近似值 x 的绝对误差限是某一位上的半个单位，该位到 x 的第一位非零数字一共有 n 位，则称近似值 x 有 n 位有效数字. 一个带有误差的数 绝对误差限满足： 相对误差限满足： 例9. 递推计算定积分 ( n =1, …, 20 )的值. 解: I0= 1 – e – 1 ≈0.63212055882856 In = 1 – n In – 1 ( n = 1,2,······,20) ∴ 0.36787944117144 0.26424111765712 0.20727664702865 0.17089341188538 0.14553294057308 0.12680235656152 0.11238350406936 0.10093196744509 0.09161229299417 0.08387707005829 0.07735222935878 0.07177324769464 0.06694777996972 0.06273108042387 0.05903379364190 0.05545930172957 0.05719187059731 -0.02945367075154 1.55961974427919 -30.19239488558378 S0=0.63212055882856, Sn=1 – nSn-1 S1 S2 S3 S4 S5 S6 S7 S8 S9 S10 S11 S12 S13 S14 S15 S16 S17 S18 S19 S20 使用递推公式 In=1-nIn-1,取准确度较高的初值(误差限小于10-15),计算过程中误差不断增大。 当n=20时，有S20=-30.19239488558378 使用逆向递推公式 算法在执行过程中误差不断