选用算法时应遵循原则.pptxVIP

计算方法当x很大时可作相应的变换 则用右端来代替左端。 第1页/共18页第一页，共19页。 计算方法（2）防止大数“吃掉”小数例 求二次方程x2-105x+1=0的根 解：按二次方程求根公式 x1=(105+(1010-4)1/2)/2 x2=(105-(1010-4)1/2)/2 在8位浮点数计算得 x1=(105+105 )/2=105 (正确）, x2=(105-105 )/2=0 (错误） 第2页/共18页第二页，共19页。 计算方法产生错误的原因① 出现大数 吃掉小数4的情况② 分子部分出现两个相近数相减而丧失有效数位常称为灾难性的抵消第3页/共18页第三页，共19页。 计算方法（3）绝对值太小的数不宜做除数当分母为两个相近数相减时,会丧失有效数字这里分子的误差被扩大104倍,再如若将分母变为0.0011,即分母只有0.0001的变化时,计算结果却有了很大变化 第4页/共18页第四页，共19页。 计算方法 解: 分子分母分别计算后相除(取9位小数) A=0.0005×0.0143×0.0012=00.0012 =0.000000009(有舍入) B=0.0003×0.0125×0.0135=00.0135 =0.000000051(有舍入) D=A/B=0.17647 真值为0，所以D只准确到小数后一位例： 计算第5页/共18页第五页，共19页。 计算方法 算法2。分成三组因子。每组只取六位小数计算 a=0.0005/0.0003=1.666667(有舍入) b=0.0143/0.0125=1.144000 c=0.0012/0.0135=0.088889 (有舍入) D=a×b×c=1. 666667× 1.144000× 0.088889 =0.169482，准确到小数后5位。bca第6页/共18页第六页，共19页。 计算方法（4）简化计算步骤，减少运算次数 x255=xx2x4x8x16x32x64x128 原先要做254次乘法现只需14次即可 又如计算多项式 p(x)=anxn? an-1xn-1 ? … ? a1x ? a0 的值 若直接计算akxk,再逐项相加，一共要做 n+(n-1)+…+2+1=n(n+1)/2次乘法和n次加法 第7页/共18页第七页，共19页。 计算方法如果将前n项提出x，则有 p(x)=（anxn-1? an-1xn-2 ? … ? a1 ）x ? a0 =((anxn-2?an-1xn-3?…? a2)x?a1）x ? a0 =(…(anx ? an-1)x?…?a2)x ? a1)x ? a0写成递推公式 第8页/共18页第八页，共19页。 计算方法于是 ,这种多项式求值的算法称为秦九韶算法,只做n次乘法和n次加法,程序实现简单 第9页/共18页第九页，共19页。 计算方法 （5）选用数值稳定性好的公式，以控制舍入误差的传播 对于一个数学问题的求解往往有多种数值方法在选择数值方法时，要注意所用的数值方法不应将计算过程中难以避免的误差放大的较快，造成计算结果完全失真。第10页/共18页第十页，共19页。 计算方法例13 计算积分 并估计误差。解 容易得到递推公式第11页/共18页第十一页，共19页。 计算方法即 为 第12页/共18页第十二页，共19页。 计算方法则准确的理论递推式 实际运算的递推式 两式相减有 第13页/共18页第十三页，共19页。 计算方法 这就是说,若 与 的误差为 = - ,即 ，则误差的递推规律为 于是 第14页/共18页第十四页，共19页。 计算方法计算 时的误差被扩大了 倍,显然算法是数值不稳定的。 如果将递推公式变换一种形式第15页/共18页第十五页，共19页。 计算方法准确的理论递推式实际运算的递推式从而有 即第16页/共18页第十六页，共19页。 计算方法于是有则这个算法的误差传递规律为 即每计算一步的误差的绝对值是上一步的十分之一，误差的传播逐步缩小，得到很好的控制，这个算法是数值稳定的 。第17页/共18页第十七页，共19页。 感谢您的观看！第18页/共18页第十八页，共19页。 内容总结计算方法。第1页/共18页。第2页/共18页。② 分子部分出现两