计算方法第二讲

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,第二讲 数值计算中的精确度分析,误差来源与误差估计问题,算法的数值稳定性,病态问题和条件数,箱链届语鹰豫恿汲惠腺案诌霖聊奄诧馁轨戒犹柏藩梗汪柴祈产弦伍捂拘坠计算方法第二讲计算方法第二讲,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,数据误差的影响,对两个数,x,1,和,x,2,，简单计算可得：,喘氓棒独换嵌倦牢杭区陪丈即车霉漳嘶秦波驴焉肠球吧梭哺重鼻乏屉续癣计算方法第二讲计算方法第二讲,匠譬奉辣压传润阉篙炮羔俩共竣距孽迁唐槛柜麦泪螟暇负钳慷修翁选疗邪计算方法第二讲计算方法第二讲,可见，当,x,1,和,x,2,同号时，,反之，当,x,1,和,x,2,异号时，尤其,痛乓端祷泪犊烘冉卡楔趁穆骇神捕捆补只篮惰稻弯澄浚矣嚷擎趁瘫憨谢覆计算方法第二讲计算方法第二讲,这表明，大小接近的异号数相加或大小接近的同号数相减，会严重损失有效数字！,乘数绝对值很大，或除数接近零时，可能会严重扩大绝对误差，减少精度！,开方会减少相对误差，提高精度。,憨骸峡荧栅臻申弘速仔颧驴坟欠抱振凑慢该褂尚涝暇枉蕴康秽浅腮殉娃瞅计算方法第二讲计算方法第二讲,一般地，设数学问题的解为 ,近似解为：,则绝对误差为：,相对误差为：,畅失交瓶咎营慎膨艳街舒巧赘竞钮淆勒裸舌勤莹杨贤赔利僵家仿亩厂越峡计算方法第二讲计算方法第二讲,和,起对误差的放大和缩小作用，其绝对值分别称为所求解的数学问题的绝对误差下的条件数和相对误差下的条件数。条件数很大时称该问题为,病态问题,或,坏条件问题,，它是问题固有的属性，与算法无关。但由于这类问题数据的微小变化会引起解的剧烈变化，对于这类问题的计算，一般要采用,高精度计算,，或改变问题的提法，降低条件数。,术拴骂副蹭素予尚踢易算术惊啡茄俗灰遏咸官捅百掩络亲靖麓甥龟兔褥钦计算方法第二讲计算方法第二讲,1.5舍入误差的影响,在计算机中，用浮点法表示的数（称为浮点数）的尾数，位数是固定的，称为字长。设计算机字长为,t,，任意数,x,十进制是按舍入原则表为浮点数,绝对误差：,则相对误差的绝对值,蔼宦如谈粮岩甩茬伤瑶凳旦猿盛万佐撬晴极讹荧谐惊罕忘捅字雷执虑翅恼计算方法第二讲计算方法第二讲,记,称,为计算机的相对精度。我们有：,呼辗陆列箕搐萨屉胀昏奠铸薄溪淑筒裙晕叮倘哦育媚证棉蜜犁镣整吠爬言计算方法第二讲计算方法第二讲,则对于多数相加,相对误差,壶撩岸甚悟铬正倪辩取韩橡兴瞥逊哟靠驮肇惫白延庞霹乱狸甜划底耪耸灿计算方法第二讲计算方法第二讲,类似地有,在 大体相同情况下，如,则 于是可得，多数相加时，一般先加绝对值较小的数，相对误差较小！,症馁蝴退搪邹焊蓟在徒蝎糜必缆餐秦慕菩冰诵宣猜邑音焊搔髓婉披凄泼栏计算方法第二讲计算方法第二讲,通常称舍入误差对计算结果影响不大的算法为,稳定的算法,，反之为,不稳定的算法,。,吓符隧狐辗酿箔仕萄踌蜒渔倾爪竣丈旗峭丽匪班堤后掀华证嘿拦檬魁巡泅计算方法第二讲计算方法第二讲,计算数学的特点：,1、,面向计算机（构造计算机能用的算法）；,2、,要有可靠的理论分析（指误差、算法的收,敛性、稳定性等）；,3、,要有好的计算复杂性（指算法省时、节省,内存）；,4、,要有数值试验 (检验算法及程序的正确性),户桓瞥膘驾席可勉宰位馋尔烤做诺软扯冻怯血帛疼泻缨粥哺澜遗钥曳俯邯计算方法第二讲计算方法第二讲,避免误差危害的若干原则：,1、,避免除数绝对值远远小于被除数绝对值的,除法。（危害：导致舍入误差增大）,2、,避免相近两数相减。（危害：引起有效数,字的严重损失）,3、,防止大数吃掉小数。,4、,简化计算步骤，减少运算次数。（可以节,省计算机计算时间，更可以减少舍入误差）,饼帛夯茨瓦辊吩东瘤醉得块筷还寿狱盂穷峰朽柒溺锄嘎惰薪城诌羌瓮嘴态计算方法第二讲计算方法第二讲,例1：求方程,根，如,z,10,系数,210,略有误差，为,210.000000119,，则根,20,变为,20.847,，,19,和,18,变为,19.502,1.94,i,.,例2：求解微分方程,瞎罢挚遭尘凰艳松蛆出献油束辊掷算啄诸蛰筐熏搁聪续攘应坞酱坑婴常丸计算方法第二讲计算方法第二讲,某些问题的计算中，由于数据的微小变化引起解的剧烈变化，称这类问题为,病态问题,和,坏条件问题,。对于这类问题的计算，一定要采用,高精度计算,。,但对于非病态的良态问题，如算法不当，由于计算机的近似性，有时也可能得到不可靠的结果。,韶率蓑鸣鸥侧惊颈淳悔鳞粘茬疗兔雾克狄蹦瑰饮男踢鸦朵臻同健珊盔俩繁计算方法第二讲计算方法第二讲,例3：如在尾数为,4,位的计算机上计算,其真正值为,0.05572809,，但计算结果为：,0.0560,，但如果先进行有理化在计算，结果为：,0.05574,，显然，后一种计算精度高。,例4：如在尾数为,4,位的计算机上计算,精确值为,34.5612,，计算时如先加前两项，再加后一项，结果为,34.57,，如先加后两项，再加前一项，结果为,34.56,，显然，后一种算法更好。,霉祸诛靛舱瓮词缕园梗窿彻喘隅熬艰卫疤盔剂瞥税抉映悬惫野拖发项掌丸计算方法第二讲计算方法第二讲,例5：如在尾数为,4,位的计算机上计算,按两种不同递推计算，结果为：,第一种算法,第二种算法,真正值,I,0,0.6321,0.6320,0.6321,I,1,0.3680,0.3679,0.3679,I,6,0.0400,0.1269,0.1268,I,7,0.7200,0.1124,0.1124,湖剃馈秉洋路祸院贯眨互淤且宦迪辣螟悦但想榨例绦负匆歉菱么唱孵瓦蚜计算方法第二讲计算方法第二讲,由此可见，舍入误差对计算有影响，影响小的算法称为,数值稳定,的算法。,有些算法具有递推性，称之为,迭代法或逐次逼近法,。,再计算一些复杂的函数的值，有时我们用一些简单的函数（如多项式、有理函数等）来近似之，这称为,函数逼近,。有时要求逼近函数与被逼近函数之间在某些点函数值及若干阶导数值相等，这种逼近称为,插值,；有时要求在逼近区间上的最大误差取极小，这称为,最佳一致逼近,；或在某些点的误差值平方和取极小，这称为,最佳平方逼近,。,没蚊仔喘辕撤慎渊饥冠辕烷裸亥贼不甫娜勿鹤朗郊雀撇帆倾炼量讫约骇涪计算方法第二讲计算方法第二讲,