1.本站不保证该用户上传的文档完整性，不预览、不比对内容而直接下载产生的反悔问题本站不予受理。

2.该文档所得收入(下载+内容+预览三)归上传者、原创者。

3.登录后可充值，立即自动返金币，充值渠道很便利

数值分析绪论数值分析绪论一、数值分析的特点数值分析实际上就是介绍在计算机上解决数学问题的数值计算方法及其理论。这门课程又称为数值计算方法求数学问题的数值解称为数值问题一、数值分析的特点一、数值分析的特点计算机硬件软件功能算术与逻辑运算核心算法数值算法非数值算法计算机硬件的特点是快软件就是利用计算机高速的简单运算去实现各种复杂的功能。一、数值分析的特点一、数值分析的特点在建立了数学模型之后并不能立刻用计算机直接求解还必须寻找用计算机计算这些数学模型的数值方法即将数学模型中的连续变量离散化转化成一系列相应的算法步骤编制出正确的计算程序再上机计算得出满意的数值结果。一、数值分析的特点一、数值分析的特点数值分析这门课具有以下几个特点：()数值分析是一门与计算机应用密切结合的实用性很强的学科()数值分析这门课程即要讨论连续变量问题又要讨论离散变量问题关心的是数值结果()数值分析这门课程已成为近代数学的一个重要分支专门研究数学问题的数值解法。思维方法是归纳法核心问题是“误差”计算数学,计算物理学,计算力学,计算化学,计算生物学,计算地质学,计算经济学,等等科学计算的核心内容是以现代化的计算机及数学软件为工具以数学模型为基础进行模拟研究。现代科学的三个组成部分:科学理论,科学实验,科学计算一、数值分析的特点二、数值分析的内容二、数值分析的内容误差分析线性代数方程组的数值解法非线性方程组的数值解法数值积分与数值微分连续函数的最佳逼近插值与拟合常微分方程的数值解法最小二乘问题的数值解法代数特征值问题矩阵分解学习重点:数值方法的计算机实现（计算机实习）构造数值方法的原理（支撑理论）评价数值方法的好坏（研究数值方法的性态、可靠性、效率）迭代法,以直代曲,化整为零,外推法本课程的基本目的是使学生通过学习和实验初步建立并理解数值计算特别是科学与工程计算的基本概念为进一步深入的学习打下坚实的基础。要掌握高级编程语言:FORTRAN,CMatlabMatlab几个显著特点用Matlab处理矩阵容易用Matlab绘图轻松用Matlab编程简洁Matlab具有丰富的工具箱内容多任务重难度大！三、基本要求三、基本要求考试评分：平时作业程序占总成绩的期末考试占总成绩的闭卷考试。作业要求：每周有课外练习两周交一次作业一学期完成个综合程序课题设计。第一章数值计算与误差分析第一章数值计算与误差分析第三节常用数学软件工具第一节数值算法第二节数值计算的误差分析第一节数值算法第一节数值算法算法：从给定的已知量出发经过有限次四则运算及规定的运算顺序最后求出未知量的数值解这样构成的完整计算步骤称为算法。评价算法的两个主要标准：速度和精度一个面向计算机,计算复杂性好,又有可靠理论分析的算法就是一个好算法计算复杂性包括时间复杂性和空间复杂性时间复杂性即计算量：一个算法所需四则运算总次数一个算法所需的乘除运算总次数单位是flop空间复杂性即存储量第一节数值算法(输入x,输出y）存储量=需乘法次加法次存储单元个。需乘法次加法次存储单元个。一般地计算n次多项式的值算法、秦九韶算法(又称为Horner算法)需乘法n次加法n次存储单元n个。有递推公式算法、需乘法n次加法n次存储单元n个。算法(输入a(i)(i=,…,n),x输出y）注意其原理为算法（秦九韶算法）(输入a(i)(i=,…,n),x输出y）注意aaa…anaaa…an…amamam…amnbbb…bsbbb…bs…bnbnbn…bns=cijmsAB的计算量为N=(mns)flopAB例矩阵乘积AB的计算量分析例：求解n元线性方程组axax…anxn=b()ananx…annxn=bn由线性方程组的克莱姆(Cramer)规则可知,如果方程组()的系数矩阵A的行列式(一般记为D=|A|)不等于零,那末,这个方程组有唯一解,而且它们可以表示为xi=DiD(i=,…,n)这里,Di是指D中第i列元素用右端（b,…bn）代替构成的行列式。Ax=bA可逆克莱姆算法步骤N=(n)n!nflopn=,N=flop第二节数值计算的误差分析第二节数值计算的误差分析用计算机解决科学计算问题时需要经历以下几个环节：实际问题的精确解与用计算机计算出来的数值结果之间就有差异这种差异在数学上称为误差。数值结果是指在选择某种数值方法之后编制程序正确输入初始数据正确的情形下所获得的数值结果。一、误差的来源一、误差的来源、数学模型数学模型是通过科学实验或者观察分析一系列数据后用数学作为工具近似地描述客观事物的一种数学表达式。在数学模型中往往包含了若干参量如物体比重、阻力系数、热交换系数等这些物理参数通常由实验仪器测得根据仪器的精密程度物理参数的确定也会产生一定的误差。一、误差的来源一、误差的来源、数学模型一、误差的来源一、误差的来源、四种误差模型误差数学模型与实际问题之间出现的这种误差称为模型误差。在例中就是模型误差。、四种误差二、截断误差分析二、截断误差分析例：(截断误差)解：利用展开式的前三项取n=截断误差二、截断误差分析二、截断误差分析三、舍入误差分析三、舍入误差分析例：舍入误差设在一台虚构的位数字的计算机上计算舍入误差为例：考虑Matlab简单程序formatlongx=y=*xz=y舍入误差对计算结果影响很大例：考虑Matlab程序x=::y=x^*x^*x^*x^*x^*x^*xplot(x,y)、计算机数系、计算机数系在实数系中每一个实数可以有无穷位不同的实数代表数轴上不同的点在计算机数系中每一个数只有有限位只有部分有理数能被计算机数系中的数精确表示。浮点数：==这种允许小数点位置浮动的表示法称为数的浮点形式。实数x的十进制浮点形式为x=aaak…c,尾数阶码ai{,,,…,},cZ基数a,()称为x的规格化的浮点形式()x的k位十进制机器数fl(x)可用两种方法定义：（）截断式fl(x)=aaakcx=aa…akak…c()四舍五入式x的k位规格化十进制机器数y=aaakc,y=fl(x)ai{,,,…,},aLcU,k是机器数的字长L、U是常数。一般数制情况：k位规格化机器数y=aaakc,=,,,,ai{,,,…,},LcU,aF(,kL,U)表示以上数集全体加数它是计算机中使用的有限离散数集(机器数系）。F(,k,L,U)中的数称为机器数。F(,,,)y=aaaac例在机器数系F(,,,)中表示fl(п)若浮点数的阶码不在L,U内则出现上溢或下溢。采用截断式fl(п)=采用四舍五入式fl(п)=但是例如在位机器数系F(,,,)中输入出现下溢输入出现上溢。、绝对误差和相对误差定义：设数a是精确值x是a的一个近似值记、有效数字、有效数字例已知准确值a=…是一个无限不循环小数求截取不同位数后的近似值和误差界。解：例:设x=,x=,x=是其近似值,问它们分别有几位有效数字位位位例:考虑x=,a=,e=,er=x=,a=,e=,er=x=,a=,e=,er=当某个量的准确值很小或很大时相对误差比绝对误差更能反映准确数与近似数的差异。一个近似值的准确程度不仅与绝对误差的大小有关而且与准确值本身的大小有关。例：设计算机数系为F(,t,L,U),将实数x=aa…atat…c,(a),用四舍五入法表为机器数fl(x)求其有效数字、绝对误差限、相对误差限。解：相对误差为绝对误差为：机器数的相对误差与x无关只与字长t有关。定义eps=t为机器精度。fl(x)有t位有效数字。注：（）在数值计算中尽可能多地保留近似数的有效数字有效数字越多相对误差越小计算结果越精确。（）对IEEE国际通用标准双精度数系机器精度Eps=e最小实数Realmin=e最大实数Realmin=e、计算机中数的计算特点:、计算机中数的计算特点:()加减法先对阶,后运算,再舍入fl(x)fl(y)=fl(fl(x)fl(y))例:在F(,,,)的计算机上计算解:==(对阶靠高阶)===()乘除法先运算,再舍入fl(x)fl(y)=fl(fl(x)fl(y))fl(x)fl(y)=fl(fl(x)fl(y))()不在计算机数系中的数做四舍五入处理注：浮点数加法不满足结合律四、数值运算中的误差估计、数值运算的绝对误差和相对误差四、数值运算中的误差估计、数值运算的绝对误差和相对误差四、数值运算中的误差估计、和、差、积、商的误差估计式四、数值运算中的误差估计、和、差、积、商的误差估计式例：设x和y的相对误差为则xy的相对误差约为四、数值运算中的误差估计数值计算中值得注意的问题四、数值运算中的误差估计数值计算中值得注意的问题()防止相近的两数相减例:当x较大时,计算例计算解：当x很小时分子出现相近数相减将以上算式变形()防止大数吃小数当两个绝对值相差很大的数进行加法或减法运算时,绝对值小的数有可能被绝对值大的数"吃掉"从而引起计算结果很不可靠在上式中重新排序计算上式====例:在F(,,,)中计算上式==()防止接近零的数做除数分母接近零的数会产生溢出错误,因而产生大的误差,此时可以用数学公式化简后再做四、数值运算中的误差估计数值计算中值得注意的问题()注意计算步骤的简化,减小运算次数简化计算步骤是提高程序执行速度的关键它不仅可以节省时间还能减少舍入误差。例：设A、B、C、D分别是、、、的矩阵试按不同的算法求矩阵乘积E=ABCD解：由矩阵乘法的结合律可有如下算法E=((AB)C)D计算量N=flopE=A(B(CD))计算量N=flopE=(A(BC))D计算量N=flop五、算法的数值稳定性五、算法的数值稳定性误差的定性分析法：即研究算法的数值稳定性。误差的定量分析法：定义：一个算法如果输入数据有扰动（即有误差）而计算过程中舍入误差不增长则称此算法是数值稳定的否则称此算法为不稳定的。此算法是不稳定的无条件稳定与条件稳定对任何输入数据都是稳定的算法称为无条件稳定。对某些数据稳定而对另一些数据不稳定的算法称为条件稳定。例在F(,,,)数系中求解二次方程：解法按求根公式解得x=,x=解法解得x=x=精确解为x=x=b>>ac六、病态问题与条件数六、病态问题与条件数对数学问题而言如果输入数据有微小扰动引起输出数据（即数学问题的解）有很大扰动则称数学问题是病态问题否则称为良态问题。条件数推导计算函数值的条件数第三节数学软件工具第三节数学软件工具一、几种常用的数学软件目前流行的数学软件主要有以下几种：符号运算软件：Mathematica,Maple矩阵处理软件：Matlab统计处理软件：SAS,Spss,Origin数学CAD软件：MathCAD、符号运算软件：Mathematica,Maple、符号运算软件：Mathematica,Maple一提起计算机求解人们就会想到数值计算。但随着符号运算及计算机代数理论的日益成熟人们又利用计算机来进行符号计算。年代人们开发了几个较为成功的通用的计算机代数系统如Reduce系统和MACSYMA系统。但由于速度慢开发不方便因此没有流行起来。在年推出的Mathematica软件才使得符号运算系统的使用普及起来。在这方面代表性的软件就是Mathematica和Maple。虽然在Matlab里也有符号计算系统但它采用的是Maple内核。、符号运算软件：Mathematica,Maple、符号运算软件：Mathematica,MapleMathematica:美国Wolfram公司出品。年推出第一版版本刚推出时受到极大好评NewYork时代周刊这样评论：“theimportanceoftheprogramcannotbeoverlooked”BusinessWeek把Mathematica列为当年十大最重要的新产品之一。Mathematica在技术界也被作为智力和应用的革命性产品受到极大欢迎。人们常说Mathematica的出现标志着现代技术计算的开始。目前版本为。、符号运算软件：Mathematica,Maple、符号运算软件：Mathematica,MapleMaple：加拿大Mathsoft公司出品。年加拿大Waterloo大学开始研究符号运算系统。取这个名字是为了表明这个软件是加拿大生产的。目前最新的版本为MapleRelease。Maple功能基本和Mathematica一致。由于这两个软件功能相似因此统一介绍这两个软件的功能和特点。功能：数值计算：初等计算、线性代数计算、数值方法等符号运算：多项式因式分解与展开、微分、积分、级数、数列等。绘图功能：方便、强大、直观的二、三维绘图。特点：简单易用的交互式操作方式:两种软件都使用英文单词来命名函数非常直观、易用。对每一个命令都能给出显示结果非常方便。)功能强大:强大的数值计算功能和威力无比的符号运算功能。)结构化的程序设计语言)Tex输出格式:Tex是科技文章的一种格式科技界应用十分广泛。、矩阵处理软件：Matlab、矩阵处理软件：Matlab简介：MATLAB的含义是矩阵实验室,是MatrixLaboratory的缩写。它的前身是LINPACK（解线性方程）和EISPACK（解特征值问题）的FORTRAN子程序库。由于它把矩阵当成一个对象因此编写程序更加直观、方便。年正式推出,最新版本为VReleaseMATLAB具有非常强大和直观的计算功能并且由于其有非常好的扩展性能现在已经成为世界上应用最广泛的工程计算软件之一特点：()强大的数值运算功能在MATLAB环境中有超过种数学、统计、科学及工程方面的函数可使用函数的命名表示自然使得问题和解答像数学公式一般简单明了让用户可全力发挥在解题方面而非浪费在电脑操作上。()数据分析和可视化功能、文字处理功能MATLAB可以绘制二、三维图形与Mathematic和Maple相比它还能处理光照模型制作出高品质的图形。功能十分强大。MATLABNotebook为用户提供了强大的文字处理功能并允许WORD访问MATLAB的数值计算和可视化结果制作科学性或工程性图文并茂的文章()高级、简单、高效的程序环境做为一种解释型的程序语言MATLAB允许使用者在短时间内写完程序,所花的时间约为用FORTRAN或C的几分之一而且不需要编译(compile)及连接(link)即能执行同时包含了更多及更容易使用的内建功能。()开放及可延伸的架构MATLAB允许使用者接触它的大多数的数学源代码检查运算法更改现有函数甚至加入自己的函数使MATLAB成为使用者所需要的环境。()丰富的工具箱MATLAB的工具箱融合了套装前软体的优点与一个灵活的开放但容易操作之环境这些工具箱提供了使用者在特别应用领域所需的许多函数。现有工具箱有：符号运算（利用MapleV的计算核心执行）、图像处理、统计分析、信号处理、通信、线性矩阵不等式、偏微分方程、高阶谱分析、财政金融、神经网络、模拟分析、控制系统、实时控制、小波分析、最优化、模糊逻辑、μ分析及合成等多种。、数学CAD软件：MathCAD、数学CAD软件：MathCADMathcad是由MathSoft公司推出的一种交互式数值系统。在输入一个数学公式、方程组、矩阵之后计算机能直接给出结果而无须去考虑中间计算过程。最令人激动的是在加入软件包自带的Maple插件后能直接支持符号运算。你可以在计算机上输入数学公式、符号和等式等很容易地算出代数、积分、三角以及很多科技领域中的复杂表达式的值并可显示数学表格和图形,通过对图形结果的分析使我们对问题的理解更加形象。特点：公式编辑功能强大、直观性强。功能：公式编辑：MathCAD的最吸引人的地方在于编辑公式非常方便、实用几乎没有不能表示的数学公式。矩阵功能定义的矩阵或矢量矩阵操作：矩阵转置、求逆矩阵、求矩阵的行列式的值。Mathcad的用户主要针对具备应用数学知识但并不要求具有较多的计算机知识的用户如工程研究人员、学生等。最新版本为MATHCAD数学计算功能解方程：一般方程、线性方程组、不定方程、常微分方程、偏微分方程等数理统计与数据处理：统计函数、统计分布函数、插值预测、曲线拟合（回归）等。积分变换：Fourier变换、Laplace变换、Z变换、小波变换等。符号处理：借用的是Maple内核。图形功能：包括绘制二维和三维图形。第三版习题P（）,,,第三版习题P（）,,,思考题P,,第二版习题P,(),思考题P()(),