《计算方法电子教案第一章序言误差.ppt》由会员分享，可在线阅读，更多相关《计算方法电子教案第一章序言误差.ppt（30页珍藏版）》请在三一办公上搜索。
1、2023/10/22,1,计算机科学与工程学院,计算方法电子教案,2023/10/22,2,同学们好！,2023/10/22,3,Computing Methods,计算方法,学习和了解科学计算的桥梁,2023/10/22,4,计算方法 能够做什么？,序言,2023/10/22,5,计算机解决实际问题的步骤建立数学模型选择数值方法编写程序上机计算,2023/10/22,6,在计算机上是否根据数学公式编程就能得到正确结果?,研究例子:求解线性方程组其准确解为x1=x2=x3=1,如把方程组的系数舍入成两位有效数字,它的解为x1=-6.222.x2=38.25 x3=-33.65.,2023/10
2、/22,7,计算方法（数值分析）研究的对象,研究数值方法的设计、分析和有关理论基础与软件实现。,计算方法又称：计算数学、数值方法、数值分析等。计算方法的分枝有最优化方法、计算几何、计算概率统计等,2023/10/22,8,课程性质,计算方法它属于数学的范畴，在很多大学的数学系有一个计算数学专业（信息与计算科学），它二年的专业课（数值代数、数值逼近、微分方程的数值解）就浓缩为计算方法这门课。计算方法是一门抽象的理论课，同时又是一门实用的技术课，它是数学与计算机联系起来的桥梁。,2023/10/22,9,课程内容,连续系统的离散化 离散性方程的数值求解,2023/10/22,10,授课目的,使学生
3、掌握用计算机解决 数学问题的一般思想方法和技术。计算方法是怎样产生的呢？我们就得从经典数学谈起了。,2023/10/22,11,经典数学,经典数学在实际应用中存在着一些难点，下面我们举几个例子。方程求根 代数方程，超越方程 我们知道代数方程一、二次都很好求解，但当次数增加时就困难了，实际上，有人论证过，当n 4时无统一的求根方式，这不涉及到求方程的技巧问题。超越方程就更是如此了。,2023/10/22,12,经典数学,B.方程组求解（Axb）我们在线性代数中学过比较著名的 Cramer法则，当det（A）!=0时，方程组有且仅有一个解。形式非常简单，但计算量非常大。这样用计算机去求解时，也将花
4、费成千上万年（几十阶的方程组）更何况上百阶的方程组在生产应用中也比比皆是。,2023/10/22,13,经典数学,用Laplace展开计算行列式，共需（n1）n！=(n+1)!以上的乘法。对于一个20阶的方程组，就需要21!5.1110E19以上的乘法。所以，在每秒作30亿次乘法的计算机上，用Cramer法则解20阶的方程组，需要的时间大约为540年。,2023/10/22,14,经典数学,C.求定积分 难点：1.原函数很难求出 2.原函数不能用有限的初等函数来表示。,2023/10/22,15,经典数学,这样的话，经典数学直接应用于工程上很难办，不能解决实际问题。人们寻找一种更简便、近似的算
5、法，这样就慢慢发展形成了计算方法这门课。其实，从计算机诞生开始，它与计算机就机密地联系在一起。,2023/10/22,16,什么是算法和计算量？,算法：从给定的已知量出发，经过有限次四则运算及规定的运算顺序，最后求出未知量的数值解，这样构成的完整计算步骤称为算法。计算量：一个算法所需的乘除运算总次数，单位是flop(Floating Octal Point 浮点八进制).计算量是衡量一个算法好坏的重要标准。,2023/10/22,17,1、模型误差2、观测误差3、截断误差4、舍入误差,第一章 误 差1.1 误差的来源,2023/10/22,18,十进制规格化浮点数 x=0.a1 a2.at10
6、J ai0,1,2,9,a10，LJU一般情况：x=0.a1 a2.atJ,=2,8,10,16,ai0,1,2,-1,LJUF(,t,L,U)表示以上数集全体加数0，它是计算机中使用有限离散集。,计算机中的数系,阶,尾数,2023/10/22,19,计算机中数的计算特点:1.加法先对阶,后运算,再舍入2.乘法先运算,再舍入3.不在计算机数系中的数做四舍五入处理,例如:在四位浮点十进制数的计算机上计算1+104 解:1+104=0.1000 101+0.1000 105=0.00001 105+0.1000 105(对阶计算)=0.10001 105=0.1000 105=104,2023/1
7、0/22,20,绝对误差：e*=x*-x,x*是近似数,x是准确数 绝对误差限*：|e*|=|x*-x|*常表示为x=x*或x*-x x*+相对误差：e*r=(x*-x)/x*,x*是近似数,x是准确数相对误差限*r：|e*/x*|=|x*-x|/|x*|*r相对误差比绝对误差更能反映准确数与近似数的差异例:考虑 1.x*=10,x=11 e*=-1 e*r=-0.1 2.x*=1000,x=1001 e*=-1 e*r=-0.001,误差定义,2023/10/22,21,如果|e*|=|x*-x|0.5 10-k 称近似数x*准确到小数点后第k位,从这小数点后第k位数字直到最左边非零数字之间
8、的所有数字都称为有效数字.,有效数字,用四舍五入得到的数都是有效数字有效数字越多,误差越小,计算结果越精确,2023/10/22,22,例如:设x1=1.73,x2=1.7321,x3=1.7320是其近似值,问它们分别有几位有效数字?,2023/10/22,23,四则运算的误差,绝对误差：e=x*-x=xdx 相对误差：er=(x*-x)/x*dx/x*=dlnx利用这个关系可以讨论四则运算的误差和函数的误差例如下列式子说明什么误差结果?d(x+y)=dx+dy dln(xy)=dlnx+dlny dln(xn)=ndlnx,2023/10/22,24,数值计算中值得注意的问题,一、防止相近
9、的两数相减(会耗失许多有效数字,可以用数学公式化简后再做).例1:各有五位有效数字的23.034与22.993相减.23.034-22.993=0.041 0.041只有两位有效数字,有效数字的耗失,说明准确度减小,因此,在计算时需要加工计算公式,以免这种情况发生.例2:当x较大时,计算,2023/10/22,25,二、防止大数吃小数.当两个绝对值相差很大的数进行加法或减法运算时,绝对值小的数有可能被绝对值大的数吃掉从而引起计算结果很不可靠.例:求一元二次方程x2-(109+1)x+109=0 的实数根.采用因式分解法,很容易得到两个根为x1=109,x2=1.如采用字长为8位的单精度计算机来
10、计算,求得根为x1109,x20.(怎样计算可得较好的结果?)两者结果不同,因为计算机计算时做加减法要“对阶”,“对阶”的结果使大数吃掉了小数.产生了误差.为了避免由于上述原因引起的计算结果严重失真,可以根据一些具体情况,存在需要把某些算式改写成另一种等价的形式.,三、防止接近零的数做除数分母接近零的数会产生溢出错误,因而产生大的误差,此时可以用数学公式化简后再做.,.,2023/10/22,26,四、注意计算步骤的简化,减小运算次数.,简化计算步骤是提高程序执行速度的关键，它不仅可以节省时间，还能减少舍入误差。例1：计算9255的值，若逐个相乘要用254次乘法，但若写成 9255=9 92
11、94 98 916 932 964 9128只需做14次乘法运算即可。例2：设A、B、C、D分别是1020、2050、501、1100的矩阵，试按不同的算法求矩阵乘积E=ABCD.解：由矩阵乘法的结合律，可有如下算法1.E=(AB)C)D.计算量N=11500flop2.E=A(B(CD).计算量N=125000flop3.E=(A(BC)D.计算量N=2200flop,2023/10/22,27,矩阵乘积AB的计算量分析,a11 a12 a13 a1na21 a22 a23 a2n.am1 am2 amm-1 amn,b11 b12 b13 b1sb21 b22 b23 b2s.bn1 bn
12、2 bnn-1 bns,=cijms,因为 cij=aik bkj 计算量为n所以上面A mn B ns的计算量为N=m n s,2023/10/22,28,学习方法,1.注意掌握各种方法的基本原理2.注意各种方法的构造手法3.重视各种方法的误差分析4.做一定量的习题5.注意与实际问题相联系,2023/10/22,29,参考书,1.计算方法引论，徐萃薇等编,高教出版社,2007.42.计算方法与算法语言上册，张德荣等编,高教出版 社3.计算方法典型题分析解集封建湖等编，西北工业大学出版社,2023/10/22,30,考试方法,1.开卷考试占60%2.平时作业及课堂回答问题占30%3.试验占10%,
