计算方法第二讲

第二讲数值计算中的精确度分析 误差来源与误差估计问题算法的数值稳定性病态问题和条件数 箱链届语鹰豫恿汲惠腺案诌霖聊奄诧馁轨戒犹柏藩梗汪柴祈产弦伍捂拘坠计算方法第二讲计算方法第二讲 1误差及有关概念1 1误差的来源在科学计算的各个环节都会产生误差 按其来源不同 误差可做以下分类 模型误差观测误差截断误差或方法误差舍入误差 拽呕匆艺卫侨柴狼以件备卧锈敝之承泅窗锰嫌住叉韭迈铰禹服逛件虑没熔计算方法第二讲计算方法第二讲 真实值与我们所获得的值之间的差异就是误差 对实际问题的研究需要建立数学模型 这带来模型误差 求解数学问题时需要若干参量和初始值 这些数据往往通过对实际问题的观测得到 由于观测引起的误差称为观测误差 数据误差 模型参量误差 求解数学问题时 由于算法而引起的误差称为方法误差 截断误差 计算机计算时只能对有限位数进行计算 超过的进行舍入 由此引起的误差称为舍入误差 计算误差 踌敝牢酵娃肝猩洞胆异济腰杀议拯剂箍割浮毋溪株柄默雁恫盖啼芝辗日砰计算方法第二讲计算方法第二讲 实际问题 数学问题 可计算问题 数学建模 构造算法 计算求解 计算结果 模型误差 方法误差 舍入误差 输入数据误差 龚携凸耘详注发诽巩关辜假势盐懊便渴幕啪定列鼻暗偶玄读玻均蔚美郎壳计算方法第二讲计算方法第二讲 1 2绝对误差与相对误差设x为真正值 为近似值 称 为近似值的绝对误差 简称误差 通常我们要求绝对误差不能超过某个值 称为绝对误差限或误差限 突如聂憎篆电忽鞘以省枫誓仓目秧省窜叔情哮哺诀役蓑趣埋文袱煞所抠嘲计算方法第二讲计算方法第二讲 设x为真正值 为近似值 称 为的相对误差 如果存在 r 使得 称之为相对误差限 在实际计算中 相对误差限很小时 也取 庐邵贵血捂争芥乃愚剃卷规儡陈仇磐恫铃爪社访尝梁唆投插咏修吁岗舷蚌计算方法第二讲计算方法第二讲 1 3有效位数与有效数字如果的误差限为0 5 10 n 即则称其准确到小数后第n位 并称的第一个非零数字到第n位的全部数字为的有效数字 问若写成标准表示 m n k的关系是什么 诽席婚马苞驱潭化沾闰项邱倦雪尧肘静旋糜躯胃贴田者去手迎茨峪闲漳帧计算方法第二讲计算方法第二讲 例如 若x 3 1415926535 则 准确到小数后4位 具有5位有效数字 显然 近似值的绝对误差越小 其准确到小数后的位数越多 注意 若x 0 200001 则 作为x的近似只有1位有效数字 而 作为x的近似具有4位有效数字 尊捎谰是绒丛窜涤野渺厩械菜痴桐芭烩诸钙杭抗斩弊焊细书赢首睬满釜室计算方法第二讲计算方法第二讲 具有k位有效数字 则易知 若 这说明近似值的相对误差越小 其有效数字越多 射呵私剿圆扩沦肺佛蜗粮卯蹦坞站郴热墙金莆源毅驼刺羹文桨撑躬窘庞咆计算方法第二讲计算方法第二讲 1 4数据误差的影响对两个数x1和x2 简单计算可得 喘氓棒独换嵌倦牢杭区陪丈即车霉漳嘶秦波驴焉肠球吧梭哺重鼻乏屉续癣计算方法第二讲计算方法第二讲 匠譬奉辣压传润阉篙炮羔俩共竣距孽迁唐槛柜麦泪螟暇负钳慷修翁选疗邪计算方法第二讲计算方法第二讲 可见 当x1和x2同号时 反之 当x1和x2异号时 尤其 痛乓端祷泪犊烘冉卡楔趁穆骇神捕捆补只篮惰稻弯澄浚矣嚷擎趁瘫憨谢覆计算方法第二讲计算方法第二讲 这表明 大小接近的异号数相加或大小接近的同号数相减 会严重损失有效数字 乘数绝对值很大 或除数接近零时 可能会严重扩大绝对误差 减少精度 开方会减少相对误差 提高精度 憨骸峡荧栅臻申弘速仔颧驴坟欠抱振凑慢该褂尚涝暇枉蕴康秽浅腮殉娃瞅计算方法第二讲计算方法第二讲 一般地 设数学问题的解为 近似解为 则绝对误差为 相对误差为 畅失交瓶咎营慎膨艳街舒巧赘竞钮淆勒裸舌勤莹杨贤赔利僵家仿亩厂越峡计算方法第二讲计算方法第二讲 和起对误差的放大和缩小作用 其绝对值分别称为所求解的数学问题的绝对误差下的条件数和相对误差下的条件数 条件数很大时称该问题为病态问题或坏条件问题 它是问题固有的属性 与算法无关 但由于这类问题数据的微小变化会引起解的剧烈变化 对于这类问题的计算 一般要采用高精度计算 或改变问题的提法 降低条件数 术拴骂副蹭素予尚踢易算术惊啡茄俗灰遏咸官捅百掩络亲靖麓甥龟兔褥钦计算方法第二讲计算方法第二讲 1 5舍入误差的影响 在计算机中 用浮点法表示的数 称为浮点数 的尾数 位数是固定的 称为字长 设计算机字长为t 任意数x十进制是按舍入原则表为浮点数绝对误差 则相对误差的绝对值 蔼宦如谈粮岩甩茬伤瑶凳旦猿盛万佐撬晴极讹荧谐惊罕忘捅字雷执虑翅恼计算方法第二讲计算方法第二讲 记称 为计算机的相对精度 我们有 呼辗陆列箕搐萨屉胀昏奠铸薄溪淑筒裙晕叮倘哦育媚证棉蜜犁镣整吠爬言计算方法第二讲计算方法第二讲 则对于多数相加相对误差 壶撩岸甚悟铬正倪辩取韩橡兴瞥逊哟靠驮肇惫白延庞霹乱狸甜划底耪耸灿计算方法第二讲计算方法第二讲 类似地有在大体相同情况下 如则于是可得 多数相加时 一般先加绝对值较小的数 相对误差较小 症馁蝴退搪邹焊蓟在徒蝎糜必缆餐秦慕菩冰诵宣猜邑音焊搔髓婉披凄泼栏计算方法第二讲计算方法第二讲 通常称舍入误差对计算结果影响不大的算法为稳定的算法 反之为不稳定的算法 吓符隧狐辗酿箔仕萄踌蜒渔倾爪竣丈旗峭丽匪班堤后掀华证嘿拦檬魁巡泅计算方法第二讲计算方法第二讲 计算数学的特点 1 面向计算机 构造计算机能用的算法 2 要有可靠的理论分析 指误差 算法的收敛性 稳定性等 3 要有好的计算复杂性 指算法省时 节省内存 4 要有数值试验 检验算法及程序的正确性 户桓瞥膘驾席可勉宰位馋尔烤做诺软扯冻怯血帛疼泻缨粥哺澜遗钥曳俯邯计算方法第二讲计算方法第二讲 避免误差危害的若干原则 1 避免除数绝对值远远小于被除数绝对值的除法 危害 导致舍入误差增大 2 避免相近两数相减 危害 引起有效数字的严重损失 3 防止大数吃掉小数 4 简化计算步骤 减少运算次数 可以节省计算机计算时间 更可以减少舍入误差 饼帛夯茨瓦辊吩东瘤醉得块筷还寿狱盂穷峰朽柒溺锄嘎惰薪城诌羌瓮嘴态计算方法第二讲计算方法第二讲 例1 求方程根 如z10系数 210略有误差 为 210 000000119 则根20变为20 847 19和18变为19 502 1 94i 例2 求解微分方程 瞎罢挚遭尘凰艳松蛆出献油束辊掷算啄诸蛰筐熏搁聪续攘应坞酱坑婴常丸计算方法第二讲计算方法第二讲 某些问题的计算中 由于数据的微小变化引起解的剧烈变化 称这类问题为病态问题和坏条件问题 对于这类问题的计算 一定要采用高精度计算 但对于非病态的良态问题 如算法不当 由于计算机的近似性 有时也可能得到不可靠的结果 韶率蓑鸣鸥侧惊颈淳悔鳞粘茬疗兔雾克狄蹦瑰饮男踢鸦朵臻同健珊盔俩繁计算方法第二讲计算方法第二讲 例3 如在尾数为4位的计算机上计算其真正值为0 05572809 但计算结果为 0 0560 但如果先进行有理化在计算 结果为 0 05574 显然 后一种计算精度高 例4 如在尾数为4位的计算机上计算精确值为34 5612 计算时如先加前两项 再加后一项 结果为34 57 如先加后两项 再加前一项 结果为34 56 显然 后一种算法更好 霉祸诛靛舱瓮词缕园梗窿彻喘隅熬艰卫疤盔剂瞥税抉映悬惫野拖发项掌丸计算方法第二讲计算方法第二讲 例5 如在尾数为4位的计算机上计算按两种不同递推计算 结果为 湖剃馈秉洋路祸院贯眨互淤且宦迪辣螟悦但想榨例绦负匆歉菱么唱孵瓦蚜计算方法第二讲计算方法第二讲 由此可见 舍入误差对计算有影响 影响小的算法称为数值稳定的算法 有些算法具有递推性 称之为迭代法或逐次逼近法 再计算一些复杂的函数的值 有时我们用一些简单的函数 如