《数值分析 绪论.pptx》由会员分享，可在线阅读，更多相关《数值分析 绪论.pptx（62页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、数值分析 能够做什么？11 IntroductionIntroduction第1页/共62页 研究使用计算机求解各种科学与工程计算问题的数值方法（近似方法），对求得的解的精度进行评估，以及如何在计算机上实现求解等。数值分析课程中所讲述的各种数值方法在科学与工程计算、信息科学、管理科学、生命科学等交叉学科中有着广泛的应用第2页/共62页应用问题举例应用问题举例第3页/共62页今有上禾三秉，中禾二秉，下禾一秉，实三十九斗；上禾二秉，中禾三秉，下禾一秉，实三十四斗；上禾一秉，中禾二秉，下禾三秉，实二十六斗。问上、中、下禾实一秉各几何？答曰：上禾一秉九斗四分斗之一。中禾一秉四斗四分斗之一。下禾一秉二斗
2、四分斗之三。-九章算术1、一个两千年前的例子第4页/共62页第5页/共62页2 2、天体力学中的、天体力学中的KeplerKepler方程方程x是行星运动的轨道，它是时间t 的函数.第6页/共62页全球定位系统：全球定位系统：在地球的任何在地球的任何一个位置，至一个位置，至少可以同时收少可以同时收到到4 4颗以上卫星颗以上卫星发射的信号发射的信号 3、全球定位系统（全球定位系统（Global Positioning Global Positioning System,GPS)System,GPS)第7页/共62页 表示地球上一个接收点R的当前位置，卫星Si的位置为 ，则得到下列非线性方程组第8
3、页/共62页记为其中，第9页/共62页4 4、已经测得在某处海洋不同深度处的水温如下：、已经测得在某处海洋不同深度处的水温如下：深度（深度（M M）466 741 950 1422 1634466 741 950 1422 1634水温（水温（o oC C）7.04 4.28 3.40 2.54 2.137.04 4.28 3.40 2.54 2.13根据这些数据，希望合理地估计出其它深度（如根据这些数据，希望合理地估计出其它深度（如500500米，米，600600米，米，10001000米米）处的水温）处的水温第10页/共62页5 5、用比较简单的函数代替复杂的函数、用比较简单的函数代替复杂
4、的函数误差为最小，即距离为最小（在不同的度量意义下）第11页/共62页6 6、人口预测、人口预测 下面给出的是中国下面给出的是中国19001900年到年到20002000年的人口数，年的人口数，我们的目标是预测未来我们的目标是预测未来的人口数（数据量较大的人口数（数据量较大时）时）19505519619606620719708299219809870519901143332000126743第12页/共62页第13页/共62页7 7、铝制波纹瓦的长度问题、铝制波纹瓦的长度问题 建筑上用的一种铝制波纹瓦是用一种机器将一块平整的铝板压制而成的.假若要求波纹瓦长4英尺,每个波纹的高度(从中心线)为1
5、英寸,且每个波纹以近似2英寸为一个周期.求制做一块波纹瓦所需铝板的长度L.第14页/共62页 这个问题就是要求由函数这个问题就是要求由函数f f f f(x x x x)=)=)=)=sin x sin x sin x sin x 给定的给定的曲线从曲线从x x=0=0到到x x=48=48英寸间的弧长英寸间的弧长L.L.由微积分学我们知道由微积分学我们知道,所求的弧长可表示为所求的弧长可表示为:上述积分称为第二类椭圆积分,它不能用普通方法来计算.第15页/共62页数值计算方法的意义、内容与方法数值计算方法的意义、内容与方法软件的核心就是算法。软件的核心就是算法。20 20 世纪最伟大的科学技
6、术发明世纪最伟大的科学技术发明-计算机计算机 计算机是对人脑的模拟，它强化了人的思维智能；计算机是对人脑的模拟，它强化了人的思维智能；计算机的发展和应用，已不仅仅是一种科学技术计算机的发展和应用，已不仅仅是一种科学技术现象，而且成了一种政治、军事、经济和社会现象；现象，而且成了一种政治、军事、经济和社会现象；没有软件的支持，超级计算机只是一堆废铁而已；没有软件的支持，超级计算机只是一堆废铁而已；算法犹如乐谱，算法犹如乐谱，软件犹如软件犹如CDCD盘片，盘片，而硬件如同而硬件如同CDCD唱机。唱机。第16页/共62页理论研究科学实验科学计算计算数学 现代科学研究的三大支柱第17页/共62页212
7、1世纪信息社会的两个主要特征：世纪信息社会的两个主要特征：“计算机无处不在计算机无处不在”“数学无处不在数学无处不在”2121世纪信息社会对科技人才的要求：世纪信息社会对科技人才的要求：-会用数学解决实际问题会用数学解决实际问题-会用计算机进行科学计算会用计算机进行科学计算第18页/共62页建立数学模型选取计算方法编写上机程序计算得出结果科学计算解题过程第19页/共62页一、计算数学的产生和早期发展计算数学是数学的一个古老的分支，虽然数学不仅仅是计算，但推动数学产生和发展的最直接原因还是计算问题计算问题。二、二十世纪计算数学的发展数值代数 最优化计算 数值逼近 计算几何 概率统计计算 蒙特卡罗
8、方法 微分方程的数值解法 微分方程的反演问题 第20页/共62页数值代数：方程求根、线性方程组求解、特征值和特征向量的计算、非线性方程组的求解；数值逼近：插值与函数逼近、数值微分 和积分、最小二乘法；微分方程数值解：常微分方程数值解；偏微分方程数值解：差分法 有限元法 有限体积法三、数值计算的主要内容第21页/共62页&教材数值分析(21世纪数学系列教材)李庆扬等（华中科技大学出版社）&参考书目 数值计算方法 徐涛 编著（吉林科学技术出版社）应用数值方法 使用MATLAB和C语言 Robert J.Schilling&Sandra L.Harris （机械工业出版社）数值分析基础教程 李庆扬
9、编（高等教育出版社）现代数值分析 李庆扬、易大义、王能超 编著 （高等教育出版社）数值分析与科学计算 Jeffery J.Leader 著，张威，刘志军，李艳红等译，（清华大学出版社)第22页/共62页2 算 法一、算法的概念 描述算法可以有不同的方式。例如，可以用日常语言和数学语言加以叙述，也可以借助形式语言（算法语言）给出精确的说明，也可以用框图直观地显示算法的全貌。定义：由基本运算及运算顺序的规定所构成的完整的 解题步骤，称为算法算法。第23页/共62页例：求解二元一次联立方程组用行列式解法：首先判别 （1）如果 ，则令计算机计算 输出计算的结果x1，x2。（2）如果D=0，则或是无解，
10、或有无穷多组解。是否为零，存在两种可能：第24页/共62页令通过求解过程，可以总结出算法步骤如下：S2 计算S3 如果则输出原方程无解或有无穷多组解的信息;否则S1 输入S4 输出计算的结果第25页/共62页输入 D=a11a22-a12a21D=0开始输出 x1,x2 结 束 No输出无解信息Yes第26页/共62页二、算法优劣的判别 计算量的大小 存贮量 逻辑结构例：用行列式解法求解线性方程组:n阶方程组，要计算n+1个n阶行列式的值，总共需要做n!(n-1)(n+1)次乘法运算。n=20 需要运算多少次？n=100?第27页/共62页一、误差的来源与分类 从实际问题中抽象出数学模型 模型
11、误差例:质量为m的物体，在重力作用下,自由下落，其下落距离s 与时间t 的关系是：其中 g 为重力加速度。3 误 差第28页/共62页 通过测量得到模型中参数的值 观测误差 求近似解 方法误差(截断误差）例如，当函数 用Taylor多项式 近似代替时，数值方法的截断误差是 与0之间。在第29页/共62页机器字长有限 舍入误差 用计算机、计算器和笔算，都只能用有限位=3.1415926 小数来代替无穷小数或用位数较少的小数来代替位数较多的有限小数，如：第30页/共62页四舍五入后在数值计算方法中，主要研究截断误差截断误差和舍入误差舍入误差（包括初始数据的误差）对计算结果的影响！第31页/共62页
12、二、误差的概念1、绝对误差与绝对误差限例:若用以厘米为最小刻度的尺去量桌子的长，大约为1.45米，求1.45米的绝对误差。1.45米的绝对误差=？不知道！是近似值 的绝对误差绝对误差,简称为误差误差。定义：设 是准确值，为 的一个近似值，称 第32页/共62页但实际问题往往可以估计出但实际问题往往可以估计出 不超过某个正数不超过某个正数 ，即即 则称则称 为绝对误差限，有了绝为绝对误差限，有了绝对误差限对误差限就可以知道就可以知道 的范围为的范围为即即 落在落在 内。内。在应用上，常常采用下列写法来刻划在应用上，常常采用下列写法来刻划 的精度。的精度。第33页/共62页2、相对误差与相对误差限
13、定义：设 是准确值，是近似值，是近似值的误差，通常取为近似值 的相对误差相对误差，记作 ，称 一般情况下是不知道一般情况下是不知道 的，怎么办？的，怎么办？第34页/共62页事实上，当 较小时是 的二次方项级,故可忽略不计.相应地，若正数满足 则称 为 的相对误差限。第35页/共62页3、有效数字定义：如果则说 近似表示 准确到小数后第 位，并从这由上述定义第 位起直到最左边的非零数字之间的一切数字都称为有效数字有效数字，并把有效数字的位数称为有效位数有效位数。第36页/共62页定义:也即，若有 位有效数字。则称则称其中，是1到9中的一个数字；是0到9中一个数字；为整数，且 若近似值 的误差限
14、是某一位的半个单位,该位到 的左边第一位非零数字共有 位,就说 有 位有效数字。第37页/共62页取 作 的近似值，就有三位有效数字；取 作 的近似值，就有五位有效数字。例如：注：注：若一近似数是由原真值经四舍五入得到，若一近似数是由原真值经四舍五入得到，则必为有效数则必为有效数.第38页/共62页4、误差限与有效数字的关系 则 至少具有 位有效数字。Th1.1：对于用 式表示的近似数 ，若 具有 位有效数字，则其相对误差限为反之，若 的相对误差限为第39页/共62页Th1.2：设反之,若 的相对误差的绝对值大于 ，其中 为整数,为正整数,。有 位有效数字。则 至多若 至多有 位有效数字,即
15、是有效数字,而 不是有效数字,则 的相对误差的绝对值必大于 ;第40页/共62页证明：不是有效数字 反之,若 则 不是有效数字，即 至多有 位有效数字.第41页/共62页4 数值运算的误差估计一一、四则运算四则运算两个近似数两个近似数 与与 ,其误差限分别为其误差限分别为 及及 ,它们进行加减乘除运算得到的误差限分别为它们进行加减乘除运算得到的误差限分别为第42页/共62页二二、函数误差估计函数误差估计当自变量有误差时，计算函数值也会产生误差，其误差限可利用函数的Taylor展开式进行估计。设 是一元函数，的近似值为 ,以 近似 ，其误差限记作 ，可用Taylor展开 介于介于之间之间.取绝对
16、值得取绝对值得第43页/共62页假定 与 的比值不太大,可忽略 的高阶项,于是可得计算函数的误差限为 当 为多元函数时计算 ,如果的近似值为 ,则 的近似为于是函数值 的误差 由Taylor展开,得：得：第44页/共62页于是误差限为于是误差限为而而 的相对误差限为的相对误差限为(1.3.1)(1.3.1)(1.3.2)(1.3.2)第45页/共62页例例:已测得某场地长已测得某场地长 的值为的值为 ,宽宽 的值为的值为 ,已知已知 ,.试求试求面积面积 的绝对误差限与相对误差限的绝对误差限与相对误差限.解解:因因 其中其中由式由式(1.3.1)(1.3.1)得得第46页/共62页而而于是绝对
17、误差限为于是绝对误差限为相对误差限为相对误差限为第47页/共62页5 数值计算中应该注意的一些原则1.要使用数值稳定的算法例：求的值的值.解：由于初值第48页/共62页递推公式按公式就可以逐步算出注意此公式精确成立What happened?!不稳定的算法！这就是误差传播所引起的危害这就是误差传播所引起的危害 !第49页/共62页由题设中的递推公式可看出，的误差扩大了5倍后传给 ，因而初值 的误差对以后各步这就造成 的计算结果严重失真。计算结果的影响，随着 的增大愈来愈严重。要怎么做才能解决这个问要怎么做才能解决这个问题呢题呢?第50页/共62页可求得I9 0.017，按改写后的公式可逐次求得
18、不妨设I9 I10，于是由将公式变为第51页/共62页I8 0.019 I7 0.021I6 0.024 I8 0.028I4 0.034 I3 0.043I2 0.058 I1 0.088I0 0.182 稳定的算法！在我们今后的讨论中，误差将不可回避，算法的稳定性会是一个非常重要的话题。第52页/共62页2.要避免两个相近的数相减在数值计算中，两个相近的数作减法时有效数字会损失。例:求的值。当x=1000，y 的准确值为0.01580(1)、直接相减第53页/共62页类似地(2)将原式改写为则 y=0.01581 第54页/共62页3.尽量避免绝对值太小的数作分母例：如分母变为0.0011
19、，也即分母只有0.0001的变化时结果相差这么结果相差这么大大!第55页/共62页4.避免大数吃小数精确解为 算法1：利用求根公式例：例：用单精度计算用单精度计算 的根。的根。第56页/共62页在计算机内，109存为0.11010，1存为0.1101。做加法时，两加数的指数先向大指数对齐，再将浮点部分相加。即1 的指数部分须变为1010，则：1=0.0000000001 1010，取单精度时就成为：109+1=0.100000001010+0.00000000 1010=0.10000000 1010第57页/共62页算法2：先解出再利用注：注：求和时从小到大相加，可使和的误差减小。例：按从小
20、到大、以及从大到小的顺序分别计算1+2+3+40+109第58页/共62页5.简化计算步骤，避免误差积累。一般来说，计算机处理下列运算的速度为例：多项式求值：给定的x 求下列n 次多项式的值。解：1.用一般算法，即直接求和法；2.逐项求和法；3.秦九韶方法；第59页/共62页算法的递推性计算机上使用的算法常采用递推化的形式，递推化的基本思想是把一个复杂的计算过程归结为简单过程的多次重复。这种重复在程序上表现为循环。递推化的优点是简化结构和节省计算量。第60页/共62页例：用秦九韶方法求多项式解：Ka5-KvK00.008330.00833v0=a510.041670.04v1=v0 x+a420.166670.15867v2=v1x+a330.50.46827v3=v2x+a2410.90635v4=v3x+a1510.81873v5=v4x+a0第61页/共62页感谢您的观看！第62页/共62页