数值分析第版第一章课件模板.ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《数值分析第版第一章课件模板.ppt》由用户(晟晟文业)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数值 分析 第一章 课件 模板
- 资源描述:
-
1、第第1章章数值分析与科学计算引论数值分析与科学计算引论 1.1 数值分析的对象、作用与特点 1.2 数值计算的误差 1.3 误差定性分析与避免误差危害 1.4 数值计算中算法设计的技术 1.5 数学软件(略)1 数值分析的定义:数值分析的定义:数值分析的主要内容:数值分析的主要内容:本课程主要内容包括插值与数据逼近、数值微分与数值积分、线性方程组的数值求解、非线性方程与方程组求解、特征值计算、常微分方程数值解等.数值分析也称计算数学计算数学,是数学科学的一个分支,主要研究用计算机求解各种数学问题的数值计算方法及其理论与软件实现.1.1 1.1 数值分析研究对象与特点数值分析研究对象与特点2 数
2、值分析既有纯数学高度抽象性与严密科学性的特点,数值分析是一门内容丰富,研究方法深刻,有自身理论体系的课程.虽然数值分析也是以数学问题为研究对象,但它不像纯数学那样只研究数学本身的理论,而是把理论与计算紧密结合,着重研究数学问题的数值方法及其理论.又有应用数学的广泛性与实际试验的高度技术性的特点,是一门与计算机使用密切结合的实用性很强的数学课程.3 三、要有好的计算复杂性,时间复杂性好是指节省时间,空间复杂性好是指节省存储量,这也是建立算法要研究的问题,它关系到算法能否在计算机上实现.数值分析的特点:数值分析的特点:一、面向计算机,能根据计算机特点提供切实可行的有效算法.二、有可靠的理论分析,能
3、任意逼近并达到精度要求,对近似算法要保证收敛性和数值稳定性,还要对误差进行分析.4 四、要有数值实验,即任何一个算法除了从理论上要满足上述三点外,还要通过数值试验证明是行之有效的.51.2 1.2 数值计算的误差数值计算的误差 1.2.1 1.2.1 误差来源与分类误差来源与分类 用计算机解决科学计算问题的过程如下:首先要建立数学模型数学模型,它是对被描述的实际问题进行抽象、简化而得到的,因而是近似的.数学模型与实际问题之间出现的误差称为模型误差模型误差.实际问题数学模型6 以上两种误差不在“数值分析”的讨论范围.实际问题数学模型 在数学模型中往往还有一些根据观测得到的物理量,如温度、长度、电
4、压等等,这些参量显然也包含误差.这种由观测产生的误差称为观测误差观测误差.数值分析只研究用数值方法求解数学模型产生的误差.当数学模型不能得到精确解时,通常要用数值方法求它的近似解.7 近似解与精确解之间的误差称为截断误差截断误差或方法误差方法误差.实际问题数学模型上机计算求出结果数值计算方法8 例如,用泰勒(Taylor)多项式 nnnxnfxfxffxP!)0(!2)0(!1)0()0()()(2,)!1()()()()(1)(nnnnxnfxPxfxR近似代替可微函数 ,)(xf则数值方法的截断误差是.0之间与在x 有了计算公式后,在用计算机做数值计算时,还要受计算机字长的限制,原始数据在
5、计算机上表示会产生误差,计算过程又可能产生新的误差,这种误差称为舍入误差舍入误差.9产生的误差用 近似代替 ,14159.3就是舍入误差.此外由原始数据或机器中的十进制数转化为二进制数产生的初始误差对数值计算也将造成影响.14159.3 R例如,0000026.0 分析初始数据的误差通常也归结为舍入误差.研究计算结果的误差是否满足精度要求就是误差估计问题.这里主要讨论算法的截断误差与舍入误差,而截断误差将结合具体算法讨论.10 若能根据测量工具或计算情况估计出误差绝对值的一个上界,即 1.2.2 1.2.2 误差与有效数字误差与有效数字 设 为准确值,x为 的一个近似值,*xxxxe*通常准确
6、值 是未知的,x因此误差 也是未知的.*e为近似值的绝对误差绝对误差,定义定义1 1称简称误差误差.*xxe 则 叫做近似值的误差限误差限,*它总是正数.11 例如,用毫米刻度的米尺测量一长度 ,读出和该长度接近的刻度 ,x*x是 的近似值,*xx它的误差限是 ,mm5.0于是 mm.5.0*xx如读出的长度为 ,mm765则有 .5.0765 x 虽然从这个不等式不能知道准确的 是多少,但可知x,5.7655.764 x结果说明 在区间 内.x5.765,5.76412 对于一般情形 ,*xx即*,*xxx也可以表示为.*xx 需要注意注意的是误差限的大小并不能完全表示近似值的好坏.13 例
7、如,有两个量 ,110 x,51000 y;1,10*xx则.5,1000*yy虽然 比 大 4 倍,*y*x但1000/5*/*yy%5.0比 10/1*/*xx%10要小得多,这说明 近似 的程度比 近似 的程度好.*yy*xx 所以除考虑误差的大小外,还应考虑准确值 本身的大小.x14 实际计算中,由于真值 总是未知的,x*xxxxeer 把近似值的误差 与准确值 的比值*exxxxxe*称为近似值 的相对误差相对误差,*x记作 .*re作为 的相对误差,*x条件是 较小,*xeer通常取此时利用,*xxe知15 相对误差也可正可负,它的绝对值上界叫做相对误差限相对误差限,xxxxexe
8、xe*)*(*是 的平方项级,*re记作 ,*r*)*(*)(2exxe*)/*(1*)/*(2xexe故可忽略不计.*xr即16%,10*xx%,5.0*yy上例中 与 的相对误差限分别为xy可见 近似 的程度比 近似 的程度好.*yy*xx根据定义,17 当准确值 位数比较多时,常常按四舍五入的原则得到 的前几位近似值 ,xx*x14159265.3 x 取3位,14.3*3x 取5位,1416.3*5x它们的误差都不超过末位数字的半个单位,,102114.32例如,002.0*3,000008.0*5即.10211416.3418 若近似值 的误差限是某一位的半个单位,*x该位到 的第一
9、位非零数字共有 位,就说 有 位有效有效数字数字.*xn*xn 表示为),1010(10*)1(121nnmaaax(2.1)其中 是0到9中的一个数字,为整数,),1(niaima,01.1021*1nmxx(2.2)定义定义2 2且19如取 作为 的近似值,14.3*x取 ,3.1416*x按这个定义,*x就有3位有效数字,*x就有5位有效数字.20 按定义,187.93,0.037856,8.0000,2.7183.的5位有效数字近似数是8.0000,而不是8,000033.8x 例例1 1 按四舍五入原则写出下列各数具有5位有效数字的近似数:187.9325,0.03785551,8.
10、000033,2.7182818.上述各数具有5位有效数字的近似数分别是因为8只有1位有效数字.注意:21如果以 m/s2 为单位,重力常数g,,m/s80.92g若以km/s2为单位,它们都具有3位有效数字,2km/s00980.0g,102180.92g按(2.1)的表示方法,,3,0nm,102100980.05g这里.3,3nm 它们虽然写法不同,但都具有3位有效数字.例例2 2因为按第一种写法按第二种写法),1010(10*)1(121nnmaaax(2.1)22 至于绝对误差限,由于单位不同所以结果也不同,m/s102122*1但相对误差都是 80.9/005.0*r 注意相对误差
11、与相对误差限是无量纲的,而绝对误差与误差限是有量纲的.例2说明有效位数与小数点后有多少位数无关.,m/s102125*2.00980.0/00005.023 从(2.2)可得到具有 位有效数字的近似数 ,其绝对误差限为 n*x,1021*1nm在 相同的情况下,越大则 越小,故有效位数越多,绝对误差限越小.mn110 nm.1021*1nmxx(2.2)24若 的相对误差限 ,*x)1(1*10)1(21nra设近似数 表示为*x)1.2(),1010(10*)1(121llmaaax其中 是0到9中的一个数字,),1(liai;1021)1(1*nra反之,则 至少具有 位有效数字.*xn
12、若 具有 位有效数字,n*x定理定理1 1ma,01为整数.则其相对误差限为25由(2.1)可得,10)1(1ma当 有 位有效数字时*xn*xxxr 反之,由*rxxx,105.01nm证明证明*101xammnma10105.011;102111na11110)1(2110)1(nmaa)1.2(),1010(10*)1(121llmaaax26知 至少有 位有效数字.*xn 定理说明,有效位数越多,相对误差限越小.27由于,4.420知 ,41a故只要取 ,4n3*10125.0r即只要对 的近似值取4位有效数字,其相对误差限就小于0.1%.20此时由开方表得 .472.420 设取 位
13、有效数字,n.1021)1(1*nra例例3 3 要使 的近似值的相对误差限小于0.1%,需取 20几位有效数字?由定理1就有%,1.010328 1.2.3 1.2.3 数值运算的误差估计数值运算的误差估计 两个近似数 与 ,其误差限分别为 及 ,*1x*2x)(*1x)(*2x);()()(*2*1*2*1xxxx);()()(*1*2*2*1*2*1xxxxxx).0()()()/(*22*2*1*2*2*1*2*1xxxxxxxx它们进行加、减、乘、除运算得到的误差限分别为29 设 是一元函数,的近似值为 ,以 近似 ,其误差界记作 ,)(xfx*x*)(xf)(xf*)(xf 一般情
14、况下,当自变量有误差时函数值也产生误差,之间介于*,*)(2)(*)*)(*)()(2xxxxfxxxfxfxf 取绝对值得*).(2)(*)(*)(*)()(2xfxxfxfxf 其误差限可利用函数的泰勒展开式进行估计.利用泰勒展开30*).(*)(*)(xxfxf 假定 与 的比值不太大,可忽略 *)(xf*)(xf *)(x的高阶项,于是可得计算函数的误差限31 当 为多元函数,如计算 时.f),(1nxxfA 的近似值为 ,nxx,1*1,nxx则 的近似值为A),(*1nxxfA于是由泰勒展开,函数值 的误差 为*A*)(Ae如果AAAe*)(nkkkknxxxxxf1*1)(),(
15、,1*nkkkexf),(),(1*1nnxxfxxf32于是误差限;)(*)(1*nkkkxxfA(2.3)而 的相对误差限为*A*)(*Arr.*)(1*nkkkAxxf(2.4)*)(AA33 已测得某场地长 的值为 ,lm110*l宽 的值d为 ,m80*d已知 .m1.0*,m2.0*ddll试求面积 的绝对误差限与相对误差限.lds 因,ldsdlslds),()(*)(*ddsllss;)(*)(1*nkkkxxfA知例例4 4解解由34其中,m80*dls而,m2.0*)(l于是绝对误差限);m(27)1.0(110)2.0(80*)(2s相对误差限*)(*)(sssrm,11
16、0*lds,m1.0*)(d*)(dls%.31.0880027351.3 1.3 误差定性分析与避免误差危害误差定性分析与避免误差危害 一个工程或科学计算问题往往要运算千万次,由于每步运算都有误差,如果每步都做误差分析是不可能的,也不科学.因为误差积累有正有负,绝对值有大有小,都按最坏情况估计误差限得到的结果比实际误差大得多,这种保守的误差估计不反映实际误差积累.36 考虑到误差分布的随机性,有人用概率统计方法,将数据和运算中的舍入误差视为适合某种分布的随机变量,20世纪60年代以后对舍入误差分析提出了一些新方法,较重要的有威尔金森(Wilkinson)的向后误差分析法向后误差分析法和穆尔(
17、Moore)的区间分析法区间分析法两种.然后确定计算结果的误差分布,这样得到的误差估计更接近实际,这种方法称为概率分析法概率分析法.37 1.3.1 1.3.1 算法的数值稳定性算法的数值稳定性 用一个算法进行计算,如果初始数据误差在计算中传播使计算结果的误差增长很快,这个算法就是数值不稳定数值不稳定的.计算 并估计误差.),1,0(ee101ndxxIxnn 由分部积分可得计算 的递推公式 nI),2,1(11nnIInn若计算出 ,0I代入(3.2),可逐次求出 的值.,21II(3.1).e1ee11010dxIx例例5 538 而要算出 就要先计算 .0I1e,!)1(!2)1()1(
18、1e21kk并取 ,7k则得 ,3679.0e13679.0e17R计算过程中小数点后第5位的数字按四舍五入原则舍入.若用泰勒多项式展开部分和 用4位小数计算,截断误差!81.1041439 当初值取为 时,用(3.2)递推006321.0II).,2,1(1;6321.0)A(10nInIInn计算结果见表1-1的 列.nI 用 近似 产生的误差 就是初值误差,0I0I000IIE它对后面计算结果是有影响的.),2,1(11nnIInn(3.2).e1ee11010dxIx计算公式为 7.55290.170440.728080.207430.216070.264220.112060.3679
展开阅读全文