数值分析笔记期末复习.doc
《数值分析笔记期末复习.doc》由会员分享,可在线阅读,更多相关《数值分析笔记期末复习.doc(23页珍藏版)》请在咨信网上搜索。
1、第一章引论1、数值分析研究对象:数值分析是计算数学的一个主要部分,计算数学是数学科学的一个分支,它研究用计算机求解各种数学问题的数值计算方法及其理论与软件实现。2、数值分析特点:面向计算机,要根据计算机特点设计切实可行的有效算法有可靠的理论分析,能任意逼近并达到精度要求,对近似计算要保证收敛性和数值稳定性要有好的计算复杂性,时间复杂性好是指节省时间,空间复杂性好是指节省存贮量,这也是建立算法要研究的问题。要有数值试验,即任何一个算法除了从理论上要满足上述三点外,还要通过数值试验证明是行之有效的。3、数值分析实质:是以数学问题为研究对象,不像纯数学那样只研究数学本身的理论,而是把理论与计算紧密结
2、合,着重研究数学问题的数值方法及理论。4、用计算机解决科学计算问题通常经历以下过程实际问题-数学模型(应用数学)-数值计算方法-程序设计-上机计算结果(计算数学)5、误差来源及分类1.模型误差从实际问题中抽象出数学模型 2.观测误差通过测量得到模型中参数的值 (通常根据测量工具的精度,可以知道 这类误差的上限值。)3.截断误差当数学模型得不到精确解时,要用数值计算方法求它的近似解,由此产 生的误差称为(截断误差)或(方法误差)4.舍入误差由于计算机字长有限,原始数据的输入及浮点数运算过程中都有可能产 生误差,这样产生的误差称为舍入误差6、五个关于误差的概念1.绝对误差2.绝对误差限3.相对误差
3、4.相对误差限(1)定义:设某一量的准确值为x,近似值为x*,则x*与x之差叫做近似值x*的绝对误差(简称误差),记为(2)性质:(1)绝对误差e(x*) 可正可负 (2) |e(x*) |的大小标志着x*的精确度 (3) 绝对误差e(x*) 未知(3)判断:绝对误差是误差的绝对值?(错)(1)定义:若指定一个适当小的正数,使则称为近似值 x* 的绝对误差限。(有时用表示近似值x*的精度或准确值的所在范围。)(2)性质:(1)在实际问题中,绝对误差一般是有量纲的,绝对误差限也是有量纲的。(2)绝对误差限是正的,有无穷多个【则比大的任意正数均是绝对误差限】(1)定义:绝对误差与准确值之比称为x*
4、的相对误差。(2)性质:(1)相对误差是个无量纲量。值小者精度高。(2)由于准确值x未知,故实际问题中,当 | 较小时,常取(1)定义:若指定一个适当小的正数 ,使则称为近似值 x*的相对误差限。(2)性质:当|较小时,可用下式计算5.有效数字(1)定义:若近似值x*的绝对误差限是某一位的半个单位,该位到x*的第一位非零数字一共有n位,则称近似值x*有n位有效数字,或说x*精确到该位。注意:近似值后面的零不能随便省去!(2)例题:取x1*= 3作为的近似值,则:一个有效数字 取 x2* =3.14 作为的近似值,则:三个有效数字 取 x3* =3.1416作为的近似值,则:五个有效数字 它们的
5、误差都不超过末位数字的半个单位。(3)性质:(1)有效数字越多,则绝对误差越小 (2)有效数字越多,则相对误差越小 有效数字的位数可刻画近似数的精确度!6、一元函数的误差估计问题:设y=f(x),x的近似值为x*,则y的近似值 y*的误差如何计算?故相应的误差限计算如下7、二元函数的误差估计问题:设y=f(x1, x2), x1, x2的近似值为x1*, x2* ,则y的误差如何计算? 故绝对误差限为8、多元函数的误差估计9、加减乘除运算的误差估计加法减法乘法除法绝对误差绝对误差限相对误差相对误差限10、算法的数值稳定性概念及运算(1)定义:初始数据的误差或计算中的舍入误差在计算过程中的传播,
6、因算法不同而异。一个算法,如果计算结果受误差的影响小,就称该算法具有较好的数值稳定性11、设计算法的五个原则(一) 要避免相近两数相减(二) 要防止大数“吃掉”小数,注意保护重要数据求和时从小到大相加,可使和的误差减小。若干数相加,采用绝对值较小者先加的算法,结果的相对误差限较小(三) 注意简化计算步骤,减少运算次数,避免误差积累(秦九韶)(四) 要避免绝对值小的数作除数(五) 设法控制误差的传播许多算法具有递推性。递推法运算过程较规律,但多次递推必然导致误差的积累。 第二章 逼近问题1,函数逼近1、插值问题: 求一条曲线严格通过数据点2、曲线拟合问题: 求一条曲线在一定意义下靠近数据点 2,
7、插值问题1、定义:求一个简单函数(x)作为f(x)的近似表达式,以满足我们称这样的问题为插值问题; 并称(x)为 f (x)的插值函数; f (x)为被插函数, x0 , x1, x2, , xn是插值节(基)点;是插值原则.3,插值多项式1、定义:求一个次数不超过n的多项式使满足插值原则(条件)称Pn(x)为 f (x)的n次插值多项式2、定理:在n+1个互异节点处满足插值原则且次数不超过n的多项式Pn(x)存在并且唯一。注:若不将多项式次数限制为 n ,则插值多项式不唯一。也是一个插值多项式,其中可以是任意多项式。4,插值问题拉格朗日差值牛顿插值二次插值基函数一阶差商k阶差商零阶差商1.差
8、商与节点的排列次序无关,称为差商的对称性2.高阶差商可由低阶差商反复作一阶差商得到,计算具有递推性3.若f(x)在a, b上存在n阶导数,则为了使得|n+1(x)|尽可能小一些,插值基点的选取原则是:使x尽可能位于区间Ix的中部,这里Ix是包含x以及所用基点的最小闭区间。1.计算量省,便于程序设计2.具有承袭性的插值公式,便于理论分析埃尔米特差值插值条件中除函数值插值条件外,还有导数值插值条件,即已知:2n+2个条件求:一个次数不超过2n+1的多项式H2n+1(x)解法1:基函数法解法2:承袭法分段低次插值原因:当插值基点无限加密时,Pn(x)也只能在很小范围内收敛,这一现象称为龙格(Rung
9、e)现象,它表明通过增加基点来提高逼近程度是不宜的。定义:设在a,b上给出插值条件:求一个折线插值函数Ih(x)满足xix0x1xnf(xi)f0f1fn1Ih(x)是a,b上的连续函数2Ih(xk)=fk,k = 0,1,n3Ih(x)在每个小区间xk,xk+1上是线性函数则称Ih(x)为分段线性插值函数数学表达: 性质:1分段线性插值多项式是分段函数;2可以预见,但n充分大时,Ih(x)能很好逼近f(x)。3Ih(x)有一个缺点:在插值点处有尖点,即一阶导数不连续,不够光滑。解决办法:三次埃尔米特插值三次样条插值两种构造方法5,最小二乘法1、 定义:已知:一组实验数据(xi,yi)(i=0
10、,1,m),且观测数据有误差求:自变量x与因变量y之间的函数关系y=F(x) ,不要求y=F(x)经过所有点,而只要求在给定点上误差按某种标准最小。2、 度量标准:(1)使残差的最大绝对值为最小 (2)使残差的绝对值之和为最小 (3)使残差的平方和为最小 3、最小二乘法多项式拟合已知:一组数据(xi,yi)(i = 0,1,m) 求:在次数不超过n的多项式中找一个函数,使误差平方和最小,即这里:解: 故: 解得:4、最小二乘法非多项式拟合参数线性 已知:一组数据(xi,yi)(i = 0,1,m)求:在函数类中找一个函数 ,使误差平方和最小,即这里:已知:一组数据(xi,yi),且每个点对应权
11、因子wi 0, (i=1,2,m).求:在函数类中找一个函数 ,使误差平方和最小,即这里:最小二乘法非多项式拟合参数非线性 第三章 定积分1,求解定积分问题方法:(求曲边梯形面积)旧:(1)牛顿莱布尼兹公式 【需要寻求原函数的困难】【已知点离散】新:(2)机械求积公式 【多项式机械求积公式】*【解决原函数的困难】【中矩形公式】*【解决原函数的困难】【梯形公式】*【解决原函数的困难】【插值型求积公式】*【解决离散问题】2,代数精度 (1)目的:数值求积方法是近似方法,为了保证精度,我们自然希望公式能对“尽可能多”的函数准确成立,这就提出了所谓代数精度的概念。(2)定义:若某个求积公式对于次数m的
12、多项式均能够准确成立,但对于m+1次多项式就不一定准确,则称该求积公式有m次代数精度。若某个求积公式对于1, x, xm 均能够准确成立,但对于xm+1就不准确成立,则称该求积公式有m次代数精度。(3)定理:当n为偶数时,牛顿柯特斯公式至少有n+1次代数精度。注:在实际应用时,出于对计算复杂性和计算速度的考虑,我们常常使用低阶偶数求积公式,代替高一阶的奇数求积公式。3,插值求积公式(1)定理:具有n+1个求积节点的机械求积公式至少有n次代数精度的充分必要条件是,它是插值型的。试总结证明机械求积公式是插值型求积公式的方法。(2)求积公式的余项若求积公式的代数精度为m,则余项形如其中K是不依赖于f
13、(x)的待定参数。3,牛顿柯特斯求积公式梯形,辛普森,柯特斯1、 定义【牛顿柯特斯】梯形公式辛甫生(Simpson)公式柯特斯(Cotes)公式一阶【2次代数精度】令f(x)=x2二阶【3次代数精度】令f(x)=x4四阶【5次代数精度】4,复化求积公式1、 定义:为了提高精度通常可把积分区间分成若干子区间,再在每个子区间上用低阶求积公式。这种方法称为复化求积法。 复化求积法就是先用低阶的牛顿柯特斯公式求得每个子区间xk, xk+1上的积分 Ik,然后再求和,用 作为所求积分I的近似值。即复化梯形公式复化辛甫生公式复化柯特斯公式 复化辛甫生公式精度优于复化梯形公式5,高斯求积公式1、 定义:机械
14、求积公式含有2n+2个待定参数:若适当选择这些参数使求积公式具有2n+1次代数精度,则这类公式称为高斯公式。中矩形公式是2、 高斯点定义:高斯公式的求积节点称为高斯点。3、 求一点的高斯公式 设一点高斯公式为则其代数精度应为4、 求二点的高斯公式再设两点高斯公式为代数精度应为5、 高斯点的性质:定理:对于插值型求积公式(4.1),其节点是高斯点的充要条件是以这些点为零点的多项式与任意次数不超过n的多项式P(x)均正交,即【证明看习题】6、 高斯勒让得公式 特别地,取a, b=-1, 1,其上高斯公式为:下面求对应的高斯点。 对于任意求积区间a, b如何求?作变换可以化到区间-1,1上,这时7、
15、 带权的高斯公式1、定义:求积公式若该公式具有2n+1次代数精度,则称这类公式为带权的高斯公式。上述(x)0是权函数。2、定理:是高斯点的充要条件是是区间a, b上关于(x)的正交多项式。3、特别:若a, b = -1,1,权函数是所建立的高斯公式为【切比雪夫高斯公式】【xk是切比雪夫多项式的零点】第四章 解线性方程组1、分类【一般形式】【矩阵形式】 系数矩阵为低阶稠密矩阵(阶数大约不超过150)【直接法】(经过有限步算术运算,可求得方程组的精确解的方法,一般用于解系数矩阵为低阶稠密矩阵)系数矩阵为大型稀疏矩阵(阶数高且零元素较多)【迭代法】(用某种极限过程去逐步逼近精确解的方法,一般用于解系
16、数矩阵为大型稀疏矩阵)2、范数向量范数矩阵范数类型在Rn上的向量x =(x1, xn)TRn在Rnn上的矩阵A=(aij),定义设对任意向量 xRn,按一定的规则有一实数与之对应,记为x,若x满足则称x为向量的范数设对任意矩阵 ARnn,按一定的规则有一实数与之对应,记为A,若A满足则称A为矩阵的范数 称为范数或最大范数【绝对值最大的元素】 称为1范数 称为2范数 称为p范数【所有值绝对值的P次方求和,再开P次方】 称为范数或行范数【所有元素绝对值之和最大的一行】 称为1范数或列范数【所有元素绝对值之和最大的一列】 称为2范数【ATA的最大特征值开平方】 称为Frobennius范数【所有元素
17、平方和开平方】性质定义:如果Rn中有两个范数 |x|s 与 |x|t ,存在常数m, M0,使对任意n维向量x有则称这两个范数等价.性质:对两种等价范数而言,某向量序列在其中一种范数意义下收敛时,则在另一种范数意义下也收敛。定理:Rn上的任意两个范数等价。【注:今后研究向量序列的收敛性时,可在任何一种范数意义下研究。】3、“病态”方程组1、定义:当一个方程组,由于系数矩阵 A 或右端常数项 b 的微小变化,引起方程组Ax=b解的巨大变化,则称此方程组为“病态”方程组,矩阵A称为“病态”矩阵,否则称方程组为“良态”方程 组,A为“良态”矩阵.2、如何划分“病态”的程度:条件数:设A为非奇异阵,称
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数值 分析 笔记 期末 复习
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【天****】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【天****】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。