一数值计算方法

《一数值计算方法》由会员分享，可在线阅读，更多相关《一数值计算方法（59页珍藏版）》请在金锄头文库上搜索。

1、引言,用计算机解决的问题分为数值问题：以科学计算为主非数值问题：以排序、检索等位事务处理为主计算机的能力有限，其基本运算部件只能做简单的四则运算和简单的逻辑的运算，怎样用计算机解决实际问题？,数值问题的解决步骤,数值分析是研究各种数学问题的数值解法及其理论的一门学科。涉及 数值计算方法（可用于编程的数值计算公式） 可靠性分析 复杂性衡量与评价,数值计算以及计算机模拟（包括当前流行的虚拟现实的方法），已经是在工程技术研究和经济、社会科学中广泛应用的方法，带来巨大的经济效益 天气预报与亿次计算机 波音777的无纸设计与有限元 CT、核磁共振 计算流体力学与爆炸工程 能源问题与大型计算,计算作为工程

2、技术研究方法,计算方法课程主要讨论如何构造求数学模型近似解的算法，讨论算法的数学原理、误差和复杂性，配合程序设计进行计算试验并分析试验结果。与纯数学的理论方法不同，用数值计算方法所求出的结果一般不是解的精确值或者准确的解析表达式，而是所求真解的某些近似值或近似曲线。,例如方程 x2=2sinx，在区间(1,2)内有唯一根, 但找不出求根的解析式, 只能用数值计算方法求其近似解。 有些数学问题虽有理论上的准确的公式解, 但不一定实用, 例如行列式解法的Cramer法则原则上可用来求解线性方程组,用这种方法解一个n元方程组,要算n+1个阶行列式的值,总共需要n!(n-1)(n+1)次乘法,当n=2

3、0时,其乘除法运算次数约需1021次方,即使用每秒千亿次的计算机也得需要上百年,而用高斯（Guass）消去法约需2660次乘除法运算,并且愈大,相差就愈大。可见研究和选择好的算法是非常重要的。,“精确”与“近似”,算法(数值算法):是指有步骤地完成解数值问 题的过程。数值算法的特点 目的性，条件和结论、输入和输出数据均要有明 确的规定与要求。 确定性，精确地给出每一步的操作(不一定都是运 算)定义, 不容许有歧义。 可执行性，算法中的每个操作都是可执行的 有穷性，在有限步内能够结束解题过程计算机上的算法，按面向求解问题的不同， 分为数值算法和非数值算法。,21世纪人才要求,科学素质：拓宽对21

4、世纪科学的了解； 加深对数学思想的理解； 培养用数学思考世界的习惯；数学能力：数学知识的运用能力； 对专业问题建立数学求解方法与 实际计算能力； 应用问题中数学创造性能力;计算知识：常用算法的数学理论； 在“误差、存贮、速度”之下的实 际计算方法； 对结果的数值分析方法;,1. 认识建立算法和对每个算法进行理论分析是基本任务，主动适应“公式多”的特点；2. 注重各章建立算法的问题的提法，搞清问题的基本提法，逐步深入；3. 理解每个算法建立的数学背景，数学原理和基本线索，对最基本的算法要非常熟悉；4. 认真进行数值计算的训练，学习各章算法完全是为用于实际计算，必须真会算。,如何进行学习？, 记好

5、课堂笔记, 保证课堂纪律, 按时完成作业, 按时上课，不迟到早退,几点要求,如何掌握每类数值方法背景：数学模型到数值问题的转化算法介绍：基本思想；算法描述算法评价：估计算法精度；分析误差的积累和传播算法改进：提高精度；减少内存占用算法比较：分析算法的优缺点,本课程对计算机相关专业的基本要求一本院校工科专业基本都要学习本课程掌握数值方法的基本原理掌握常用的科学与工程计算的基本方法能用所学方法在计算机上算出正确结果（编程）,本章内容1误差的来源及分类 2误差的度量 3误差的传播 4减少运算误差的原则,第一章计算方法与误差,小结,要求掌握的内容,第一章计算方法与误差,概念 包括有效数字、绝对误差、绝

6、对误差限、相对误差、相 对误差限等误差截断误差、舍入误差的详细内容，误差种类等分析运算误差的方法和减少运算误差的若干原则,第一章计算方法与误差,1.1 误差的来源及分类 早在中学我们就接触过误差的概念，如在做热力学实验中，从温度计上读出的温度是23.4度，就不是一个精确的值，而是含有误差的近似值。事实上，误差在我们的日常生活中无处不在，无处不有。如量体裁衣，量与裁的结果都不是精确无误的，都含有误差。,在用数值方法解题过程中可能产生的误差归纳起来有如下几类：1. 模型误差2. 观测误差3. 截断误差4. 舍入误差,第一章计算方法与误差,用数学方法解决一个具体的实际问题，首先要建立数学模型，这就要

7、对实际问题进行抽象、简化，因而数学模型本身总含有误差，这种误差叫做模型误差数学模型是指那些利用数学语言模拟现实而建立起来的有关量的描述数学模型的准确解与实际问题的真解不同,1. 模型误差,在数学模型中通常包含各种各样的参变量，如温度、长度、电压等，这些参数往往是通过观测得到的，因此也带来了误差，这种误差叫观测误差数学模型中的参数和原始数据，是由观测和试验得到的由于测量工具的精度、观测方法或客观条件的限制,使数据含有测量误差,这类误差叫做观测误差或数据误差根据实际情况可以得到误差上下界数值方法中需要了解观测误差,以便选择合理的数值方法与之适应,2. 观测误差,精确公式用近似公式代替时,所产生的误

8、差叫截断误差 例如, 函数f(x)用泰勒(Taylor)多项式,3. 截断误差,（介于0与x之间）,近似代替，则数值方法的截断误差是,截断误差的大小直接影响计算结果的精度和计算 工作量，是数值计算中必须考虑的一类误差,在数值计算中只能对有限位字长的数值进行运算需要对参数、中间结果、最终结果作有限位字长的处理工作，这种处理工作称作舍入处理用有限位数字代替精确数，这种误差叫做舍入误差，是数值计算中必须考虑的一类误差,4. 舍入误差,第一章计算方法与误差,例如在计算时用3.14159近似代替，产生的误差R= -3.14159=0.0000026就是舍入误差。 上述种种误差都会影响计算结果的准确性，因

9、此需要了解与研究误差，在数值计算中将着重研究截断误差、舍入误差，并对它们的传播与积累作出分析。,1.3 误差的度量,1.3.1 绝对误差和绝对误差限 定义1.1 设精确值x的近似值 x* ，称差 e(x*) =x-x* 近似值x*的绝对误差，简称误差。 e(x*)又记为e* 当e*0时，x*称为弱近似值，当e*0时，x*称为强近似值|e*|越小， x*的精度越高 由于精确值一般是未知的,因而e* 不能求出来, 但可以根据测量误差或计算情况设法估计出它的取值范围，即误差绝对值的一个上界或称误差限。,1.3 误差的度量,定义1.2 设存在一个正数，使,则称为近似值的绝对误差限，简称误差限或精度。

10、实际应用中经常使用这个量来衡量误差限, 这就是说, 如果近似数 的误差限为 , 则表明准确值 x 必落在 上, 常采用下面的写法,来表示近似值的精度或准确值x所在的范围。,1.3 误差的度量,例1 设x =3.1415926 近似值x* =3.14,它的绝 对误差是 0.001 592 6，有, x-x*=0.0015926 0.002=0.210-2例2 又近似值x* =3.1416，它的绝对误差是 0.0000074，有 x-x*=0.0000074 0.000008=0.810-5例3 而近似值x* =3.1415，它的绝对误差是 0.0000926，有 x-x*=0.0000926 0

11、.0001=0.110-3可见，绝对误差限*不是唯一的，但*越小越好,1.3.2 相对误差和相对误差限,只用绝对误差还不能说明数的近似程度,例如甲打字每100个错一个,乙打字每1000个错一个,他们的误差都是错一个,但显然乙要准确些,这就启发我们除了要看绝对误差外,还必须顾及量的本身。定义1.3 绝对误差与精确值x的比值,称为相对误差。 简记为,1.3.2 相对误差和相对误差限,相对误差越小,精度就越高,实际计算时,x通常是不知道的,因此可用下列公式计算相对误差,定义1.4 设存在一个正数 ，使,则称 为近似值 的相对误差限。 简记为,1.3.2 相对误差和相对误差限,例4. 甲打字每100个

12、错一个，乙打字每1000个 错一个，求其相对误差解： 根椐定义:甲打字时的相对误差 乙打字时的相对误差,1.3.3 有效数字,定义1.5 设x的近似值,其中 是0到9之间的任一个数,但p是正整数, m是整数,若,则称 为x的具有n位有效数字的近似值， 准确到第n位， 是 的有效数字。,规格化,1.3.3 有效数字,例5. 3.142作为的近似值时有几位有效数字解： m = 1 m n =1n =-3 所以 n =4，具有4位有效数字,例6. 当取3.141作为的近似值时 -3.141=0.3141592101 -0.3141101 0.0000592 101 0.0005=1/2 10-2 m

13、-n=1-n=-2 所以n=3具有3位有效数字推论 如果近似数x*误差限是某一位的半个单位, 由该位到x*的第一位非零数字一共有n位, x*就有n位有效数字,也就是说准确到该位.,再如3.1416作为的近似值时 -3.1416 = 0.3141592101-0.31416101 0.00000074 101 0.00000740.00005 0.5 10-4 m-n=1-n=-4 所以 n=5x*= 3.1416有5位有效数字,关于有效数字说明 用四舍五入取准确值的前n位x*作为近似值,则 x*必有n位有效数字。如3.142作为 的近似值 有4位有效数字，而3.141为3位有效数字 有效数字相同的两个近似数，绝对误差不一定 相同。例如，设x1*=12345,设x2*=12.345,两者 均有5位有效数字但绝对误差不一样 x- x1* =x- 12345 0.5= 1/2 100 x- x2* =x- 12.3450.0005=1/210-3 把任何数乘以10p(p=0,1,)不影响有效位数 准确值具有无穷多位有效数字,如三角形面积 S=1/2ah=0.5ah 因为0.5是真值,没有误差 *=0,因此n,准确值具有无穷位有效数字,