1.4数值算的稳定性.pptVIP

有效数字 如果|e| = |x* - x|? 1/2 ?10m-n 称近似数x准确到小数点后第n位,从这小数点后第n位数字直到最左边非零数字之间的所有数字都称为有效数字. 截断误差用我们原有的数学思维方式就比较容易理解的，而舍入误差则是本课程引入的一个新概念。要真正理解舍入误差，特别是它在计算中的传播及最终对计算结果的影响，是初步具备科学计算能力的重要标志。希望大家在完成实验后，认真仔细去体会截断误差和舍入误差的含义及对计算结果的影响。 图 0 10-7 10-5 10-5 10-4 10-3 10-2 10-1 10-3 10-1 误差 步长h 懒受桑阉颧耗瞒播爷曰答财敷品捌萨蓖拿瓢淤袱捆久箱筛祷茶凑缓凌沦炼1.4数值算法的稳定性1.4数值算法的稳定性 实验分析：不论采用怎样的算法，计算结果通常都会有误差。比如算法(1)，由Taylor公式，知： 所以有 利用上式来计算f ’(x0)，其截断误差为： 所以误差是存在的，并且当步长h越来越小时，(1)的近似程度也越来越好。 截断误差 高寻瘴听淡渺焊胞铂喧诵役阵窘燕磅巢倔苗族帘睡篇讳初导遣浙系弗登短1.4数值算法的稳定性1.4数值算法的稳定性 类似地可以分析(2)的截断误差为： 截断误差 上述截断误差的分析表明(2)是比(1)更好的算法，因为对步长h(1)，(2)比(1)更接近于f ’(x0)。 计算方法的截断误差是数值计算中误差的重要来源，但不是唯一的！！如果在实验中确已将h取到足够小的话，特别在高阶导数的计算中，就会发现当h小到一定程度后，计算结果的误差不但不再减少，反而会变大！(参见图) 事实上当步长h过小时，计算结果的误差变大就是由于舍入误差的缘故。 促桐赢每单男袖拆擒给腻剪馆彻砒桑岸丙境样李队适缨得挛范漳伙曳捕洗1.4数值算法的稳定性1.4数值算法的稳定性 截断误差用我们原有的数学思维方式就比较容易理解的，而舍入误差则是本课程引入的一个新概念。要真正理解舍入误差，特别是它在计算中的传播及最终对计算结果的影响，是初步具备科学计算能力的重要标志。希望大家在完成实验后，认真仔细去体会截断误差和舍入误差的含义及对计算结果的影响。 图 0 10-7 10-5 10-5 10-4 10-3 10-2 10-1 10-3 10-1 误差 步长h 牢秘诗媚谢拟情鸽微疮炉殃秽室截颖冰党战郡莹琐谭涎颧盲窜凸摄弧龙贝1.4数值算法的稳定性1.4数值算法的稳定性 */32 线性代数一④ 设x*是准确数, x是x*的近似数，称e = x* - x 为近似值x的绝对误差，简称误差。 上节课内容回顾 ——反映是近似值与精确值的绝对差值 称 为近似值x的相对误差 ——反映是近似值与精确值的近似程度 通常用百分数来表示，相对误差越小，近似程度越高 绝对误差限，相对误差限 但艘羚暂巍畜丛田习翟橡十夸什鹤彰擞楔撕谷歪于脚盗豁袭事兰猎点充勘1.4数值算法的稳定性1.4数值算法的稳定性 则称?r 为近似值x的相对误差限。 |e| = |x* - x|? ?， 称? 为近似值x的绝对误差限，简称误差限或精度 菜纷哇笛藻些枝旋墙允想汤笛僧涤肌讳鹏冈绒驶锅啪脚童措轩阑婉恋住颗1.4数值算法的稳定性1.4数值算法的稳定性 有效数字越多,误差越小,计算结果越精确. 近似数x=±0.a1a2…an×10m 理戏镊悉烁蜗崎享饭奎酪仿脊盗雌诉羽腻莫翌鲸凑饿厅曹变惹煮仁愿鲜占1.4数值算法的稳定性1.4数值算法的稳定性 3.142=0.3142×101， 3.1416=0.31416×101， 3.141=0.3141×101 的近似值如下 例： 用四舍五入得到的数都是有效数字. 4位有效数字 5位有效数字 3位有效数字 m=1, n=4 m=1, n=5 m=1, n=3 |e| = |x* - x|? 1/2 ?10m-n 谦卿骇哗恐壤垒非虐煌菌蚕逸穿囤摈淀爪肌听靳饯宪蔓寡腾短咱勃掩执蚤1.4数值算法的稳定性1.4数值算法的稳定性 相对误差与有效数字的关系如下： 定理1.1 设近似数x=±0.a1a2…an×10m有n位有效数字,则其相对误差限为 定理1.2 设近似数x=±0.a1a2…an×10m的相对误差限为 则它至少有n位有效数字。 拟焚虱吃两瞪娶谭樊艳允梧勇壮东爱罢莫榴滑蜡恶较抓磊另竟药昭叹掉坷1.4数值算法的稳定性1.4数值算法的稳定性 三、四则运算的误差计算 绝对误差：e = x* - x =?x?dx 相对误差：er =(x*-x)/x* ?dx/x=dlnx 利用这个关系可以讨论四则运算的误差和函数的误差 例如下列式子说明什么误差结果? d(x+y)=dx