数值分析第1,2章new

《数值分析第1,2章new》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、1数值分析NumericalAnalysis主讲张学莹zhangxy@hhu.edu.cn2教材(TextBook)李庆扬、王能超、易大义，数值分析（5版）,北京：清华大学出版社,2008参考书目(Reference)易大义，陈道琦，数值分析引论,浙江大学出版社,2000蔡大用，白峰彬，高等数值分析，清华大学大学出版社，1997J.Stoer&R.Bulirsch，IntroductiontoNumericalAnalysis(SecondEdition)，NewYork:Springer,1993其他（略）3课程简介数值计算方法是计算数学的一个分支,又称数值

2、分析或计算方法,它是研究用计算机求解各种数学问题的数值方法及其理论的一门学科,是程序设计和对数值结果进行分析的依据和基础。应用计算机解决科学技术和工程问题的步骤：(1)提出实际问题(2)建立数学模型(3)选用数值计算方法(4)编程上机计算得出数据结果。4提问：数值分析是做什么用的？5由于计算机和科学计算的发展，求解各种数学问题的数值方法也愈来愈应用于科学计算的各个领域，新的交叉学科不断涌现，如计算力学、计算物理、计算生物学、计算经济学等，统称为科学计算。科学计算、理论研究、试验研究6计算数学研究适用于计算机编程的算法基础性、应用性、边缘性面向计算机的算法：串行、并行好

3、的算法的标准：可靠的理论性、计算复杂性好（时间复杂性、空间复杂性）数值分析的特点：面向计算机、有可靠的理论分析、好的计算复杂性、数值试验数学软件软件包LAPACK（EISPACK,LINPACK)、IMSL、NGA(1000多个）、MATLAB7一、误差的来源1.模型误差:在建立数学模型过程中,不可能将所有因素均考虑,必然要进行必要的简化,这就带来了与实际问题的误差。2.观测误差:测量已知参数时,数据带来的误差。3.截断误差:在设计算法时,近似处理带来的误差。第二节数值计算的误差8函数用泰勒多项式近似代替时，有误差其中在与之间。这种误差就是截断误差。例如：94.舍

4、入误差:计算机的字长是有限的,每一步运算均需四舍五入,由此产出的误差。例如：用3.14159近似代替,产生的误差就是舍入误差。10二、误差的基本概念1.误差和误差限设是准确值x的一个近似值,称为近似值的绝对误差,简称误差.又简记.误差是无法计算的(因为准确值x不知道),但可以估计出它的一个上界。即,称是近似值的绝对误差限,简称误差限.误差是有量纲的，可正可负。112.相对误差和相对误差限称为近似值的相对误差。简记为。实际计算中,由于准确值x总是未知的,且由于12相对误差是无量纲的,也可正可负,它的绝对值的上界称为该近似值的相对误差限,记作简记为是的平方项级,故当较小时

5、,常取13三、有效数字如果近似值x*的误差限是其某一位的半个单位,该位到x*的第一位非零数字共有n位,我们称x*有n位有效数字。定义：误差限不超过该位的半个单位自左向右看,第一位非零数字14所以x*=3.14作为π的近似值,有3位有效数字；而取x*=3.1416时,所以x*=3.1416作为π的近似值,有5位有效数字。例π=3.1415926535,取x*=3.14时，15下面给出有效数字的另一等价定义用表示x的近似值，并将表示成若其误差限,则称具有n位有效数字,这里m是整数,a1,a2,,ak为0~9中的一个数字,且a10.定义：自左向右看,第一位非零数字16

6、例=3.1415926535,取=3.14时，即m-n=-2,m=1,n=3,所以=3.14作为近似值时,就有3位有效数字。17四、相对误差限与有效数字的关系证明故定理1有n位有效数字，。则其相对误差限为设近似值18则它至少有n位有效数字。定理2设近似值的相对误差限不大于故至少有n位有效数字。已知证明19例解由于所以由定理1有即得要使的近似值的相对误差限小于0.1%,要取几位有效数字。故只要对的近似数取4位有效数字,其相对误差因此,可取就可小于0.1%,20五、数值运算的误差估计（一）、算术运算的误差所以和或差的误差限是误差限之和。以上的结论适用于任意多个近似数的和

7、或差。同理可得：乘、除运算的误差,以两数为例写出1.21它是对数函数的微分。设u=xy,则lnu=lnx+lny,因而dlnu=dlnx+dlny这就是说,乘积的相对误差是各乘数的相对误差之和,相对误差限是各乘数的相对误差限之和。2.的相对误差是22即商的相对误差是被除数与除数的相对误差之差,但相对误差限是各乘数的相对误差限之和.由此可得:任意多次连乘、连除所得结果的相对误差限等于各乘数和除数的相对误差限之和。同样,若u=x/y,则lnu=lnx–lny,因此dlnu=dlnx–dlny即23例1解所以从而得到u的相对误差限等于乘数x、y和除数z、ω