在Mathematica中的数值计算误差.pdfVIP

第27卷 第 4期 延边大学学报 (自然科学版) Vo1．27NO．4 2001年 12月 Journa[ofYanbianUniversity(NaturalScience) Dec．2001 文章编号：1004—4353(2001)04—0257—06 在 Mathematica中的数值计算误差 张寿， 金哲， 赵之莘， 白京洙 (延边大学理工学院物理系，吉林 延吉 133002) 摘要：利用Math~mtica软件，探讨了不同表达式下的多项式、级数展开式、方程解等的运算误差 关键词 ：Mathematics；数值计算 ；误差 中图分类号 ：TB115 文献标识码：A 随着科技的发展，计算机已广泛地应用在各个领域 ，特别是在物理学研究中起着很重要 的作用，但是计算机的数值运算结果并不是绝对准确无误的．由于计算机的内存单元是由有 限的二进值 (bit)组成的，因此它不能存储无限长的数，所以在实际应用中一般采用实数形 式进行计算．这样既可以扩大数的表示范围，又可保证必要的有效位数 ，这是用整数形式所 做不到的 但一个实数化成二进制数时，会出现一些误差，例如，要用有限位数 的二进制数来 代表一个十进制数0．1是不可能的，只能近似的代表 ，所 以很多实数不能精确的表示出来． 现代科学研究不只是简单的实数代人，而是很多公式及函数反复运算，互相嵌套的结果 ，所 以很可能会出现较大的误差 ，影响结果的真实性 ．Mathematics[1．-3]是一种强大的数学计算 、 处理和分析的工具，主要应用于解决 自然科学领域中的问题．为了尽量避免或减小运算中这 些误差的出现，我们利用Mathematiea分三部分讨论了不同情况下的运算误差． 1 不同表示方式下的多项式 以(z一1)。为例，将其展开并变形为z 一3x +3x一1和 z((z一3)z+3)一1，对这 3种表达形式进行 以下误差分析 A．为了得出一个多项式的不同表达方式，首先打开 AlgebraHorner软件包 ，运用函数 Plot画出它们的二维图形 ，并根据图形分析误差的大小． AlgebraHornef Plot[ 一3 +3x一1，fx,0．99998，1．00002}，Frnme-~True，P10tS刚e—·Au自。ma6c】 h=Homer[x一3x2+3x一1] 一 I+x(3+(一3+X)X) Plot[h，{x，0．99998，1．00002t，Frame-}True，PlotStyle-~Automatic】 Plot[(x一1)，{x，0．99998，1．00002}，PlotRange--*All，Franle-~True，Pl0tstyk-． Automatic] 收稿 日期 ：2001—07—27 作者简介：张寿(1965一)．男(朝鲜旅)．吉林龙井^，延边大学理工学疏物理系副教授 延边大学学报 (自然科学版) 第 27卷 ／ ， ／ ． ／… ／／／～ ／ 围 1 函数 一3 +3 一1的囤形 围2 函数 一1十 C3+ 卜3+ )zJ的囤形 围3 函数 c 1)的图形 其 中，图1为 z3—3x +3x—l的图形 ，图2为 ((z一3)z+3)一1的图形，图 3为 ( —1)的图形 ，z的取值范围都是0．99998到 1．00002，从上面3个图形的比较中可以看 出，3条 曲线的形状基本一致 ，但图l曲线的光滑度最差 ，图3曲线的光滑度最好．也就是说， 3种表达方式在相同取值范围内产生的误差大小不同，其