收藏 分享(赏)

数值分析V4提要.doc

上传人:hwpkd79526 文档编号:7187983 上传时间:2019-05-09 格式:DOC 页数:12 大小:374KB
下载 相关 举报
数值分析V4提要.doc_第1页
第1页 / 共12页
数值分析V4提要.doc_第2页
第2页 / 共12页
数值分析V4提要.doc_第3页
第3页 / 共12页
数值分析V4提要.doc_第4页
第4页 / 共12页
数值分析V4提要.doc_第5页
第5页 / 共12页
点击查看更多>>
资源描述

1、1数值分析提要第一章: 绪论 关键词:误差、绝对误差、相对误差、有效数字、稳定性1 科学计算的重要性:当今世界,科学与工程计算是高科技,它已经成为衡量一个国家的实力的指标之一;它在科学、工程、社会经济活动中的作用日益重要,成为继理论研究、科学实验之后,三大科学技术研究手段之一,是不可缺少的重要技术;2误差的来源:模型误差、观测误差、截断误差(方法误差) 、舍入误差(计算误差) ;误差分类:绝对误差、相对误差;强近似、弱近似、误差限、相对误差限;3 精度的定义:有效数字 N(精确到第 N 位) ;4 衡量算法的指标:稳定性、计算精度、计算速度;5 四则运算的误差估计;有效数字与相对误差的关系;6

2、 减少误差的一些基本原则:避免大数除小数、避免相近数相减、避免大数与小数的加减、简化运算步骤减少运算次数。一、定理 设 的近似值 有 位有效数字,x*xn,则其相对误差限 ,反之,)1010(10(2* nmaax )1(*02nra若其相对误差限满足: ,则 有 位有效数字。)1(1*(nr*xn二、误差限的运算:2*1*12*1*2*1*21|)()(|)/(|)(xx三、误差的传播 )(|)( )(,(, )()(* *1*11*xfxf xxffff iiniinn四、稳定n1第二章 代数插值关键词:拉氏插值、艾氏插值、样条插值、等距插值、分片插值、线性插值、二次插值1设 是区间 上的

3、连续函数,已知它在 上 个互异点 ,)(xfyba, ba,1nniy0若代数多项式 在点 处满足 则称 为函数 的插Pi 2,0,)(iyxPi )(xP)(xf2值多项式, 为插值节点。nix0满足上述条件的次数不超过 的多项式 是存在且惟一的。n)(xPn常用的代数插值公式有:拉格朗日型 nj jjiinj yxyxlP00)()(和牛顿型公式)()(, )(,)( 11010 102010 nnn xxxf xffx 插值余项为)()!(),)( 1)1(10 ffxR nnnn 其中 。),(),(01baiin2埃尔米特插值公式为 njjjnjjn yxyxH001 )()()(其

4、中njxlx xlxxljj jjjniijjj ,210),() )()(21)(21 22 的表达式同拉格朗日型公式中的基函数。插值余项为(lj)()!2() 12)(12 xfnxRnnn 3三次样条插值多项式,详见参考文献 1,2,5 等。第三章 函数逼近关键词:最佳一致逼近、最佳平方逼近、最小二乘逼近、数据拟合、多项式逼近、权函数、正交多项式、勒让得多项式、切比雪夫多项式1、函数逼近的概念。设 中一组线baCnixbaCxfi ,0),(,)( 为性无关的函数, 在某种范数意义下)(.,*10 xfspann与若 3距离最近,即 ,则称 为 ffxaxnii min,)()(0 *满

5、 足 )(*x在函数类 意义下的最佳逼近函数。特别当范数为 时,分别)(xf 中 在 范 数 和2为最佳平方逼近和最佳一致逼近。常用的函数类 次多项集合。 次最佳平方逼近多项式和最佳一致逼近多项式都n为 n是存在惟一的。次最佳平方逼近多项式 的系数 满足正规方程组n)(*x*ianifajji ,10,),(0*逼近误差 ,这里内积定义为2*)(fxfnn满 足 jifa0*),(。为 权 函 数,)(),(dgfgfba2、正交多项式的引入,可能简化最佳平方逼近多项式的正规方程组的求解。常用的正交多项式有:切比雪夫多项式、勒让德多项式、拉盖尔多项式以及埃尔米特多项式。3、对数据 进行曲线拟合

6、,常用最小二乘法(最佳平方逼近) 。其Niyxi ,10),(拟合曲线及误差的表达式同连续函数的情形完全一样。只需把内积理解为i iiixgfgf0)()(),(即可。第四章 数值积分和数值微分关键词:求积代数精度、机械积分公式、高斯点积分公式、牛顿-柯特斯求积公式、矩形公式、梯形公式、辛普森公式、复化求积公式、差分、向前差分、向后差分、一阶差分、二阶差分、高阶差分、外推法数值积分是求定积分 的近似算法。badxfI)(1、 牛顿-柯特斯求积分式牛顿-柯特斯求积公式是节点等距的插值型求积公式。将 等分,其节点为nba,则牛顿-柯特斯求,xfxnabhnkhax kkk )(,/)(,1,0处

7、的 函 数 值 为节 点积公式为)()()( 0)0 knkknkba fCxfAdxf 4其中 为柯特斯系数。)(nkk,CA为 求 积 系 数时,求积公式为梯形公式:1nbafbfabdxfba ),(21)(2)( 3时,求积公式为辛普森公式:2bafbffafbdxfba ),(2801)(2(4)6)( 452、 复化求积公式复化梯形公式 :nh/)( )(280)(4)(2()( 4112 fhabfxfxfaffRTdxf niniiinba 3、求积公式的代数精确度概念一个求积公式 的多项式都能精确地成立,但 bankkmxfAxf对 于 次 数0)()(对于 次多项式 不能精

8、确成立,则称该求积公式具有 次代数精确度。1m1m形如 的求积公式至少具有 次代数精确度的充要条件是,该bankkxfdxf0)()( n求积公式为插值型的。牛顿- 柯特斯公式是插值型公式。4、 高斯型求积公式以 为节点 bxxaxxb nn 101)()(, 的 零 点的 正 交 多 项 式上 带 权 作插值多项式所得到的插值型求积公式,可使该求积公式的代数精确度提高到 次。12这样的节点 称为高斯点,相应的求积公式称为高斯型求积公式,其一般表达式为ni0 banknkfRxfAdxffI0)()()(其中系数 、节点 可查表,余项kAk),(,)()!2( 21)2( badxfnfRba

9、nnn 常用的高斯型求积公式有:高斯-勒让德求积公式、高斯- 切比雪夫求积公式、高斯-拉盖尔求积公式及高斯-埃尔米特求积公式。5、 数值微分数值微分的求法,总的来说有两大类方法:泰勒展开方法和用插值函数求微分方法。5第五章 解线性代数方程组的直接法关键词:高斯消去法、LU 分解、Cholesky 分解、追赶法、范数、向量范数、1 范数、2范数、P 范数、无穷范数、矩阵范数、列范数、行范数、F 范数、谱半径、条件数、误差估计设线性方程组为 Ax=b,其中 A 为 n 阶非奇异矩阵,x,b 为 n 维列向量。1 高斯消去法、高斯-若当消去法和主元素消去法见参考文献 1。2 矩阵三角分解法。高斯消元

10、过程实现了 A 的一个三角因子分解 A=LU,其中 L 为单位下三角阵,U 为上三角阵,此分解称为 A 的 Doolittle 分解。A 有唯一 Doolittle 分解的充要条件是 A 的前 n-1 个顺序主子式均不为零。如果 A 对称正定,则有唯一分解 A= ,这里 为下三角阵,此分解称为平方根分解或TLCholesky 分解.当 A 对称正定时也可以进行如下分解 A= ,其中 L 为单位下三角阵,D 为对角阵.D解方程 Ax=b 等价于解两个三角方程组这种解方程组的方法称为直接三角分解法.yUxbL3.三对角方程组的追赶法求解见参考文献 1,4.4.向量、 矩阵的范数及方程组的性态.向量

11、与矩阵范数的定义见参考文献 1.常用的向量范数是 p-范数:, ( 正整数)nipipx11)|(,当 p 为 1,2,时分别为:1范数: niixA1|2范数: 2112)|(nii-范数 |maxiniA与上述向量范数相容的矩阵范数分别是:1范数: iijnj11|2范数: 2max2)(AT-范数 |1ijni其中 为 的最大特征值.)(maxAT6设| . | a 为非奇异矩阵 A 的某种算子范数 ,当 A 的条件数时,方程组 Ax=b 是”病态”的,此时求解比较困难.条件数1)(aaCond越大,则病态越严重,从而越难求解,如果系数矩阵的条件数不是很大,则方程是”良态”的.第六章 解

12、线性代数方程组的迭代法关键词:迭代法、雅可比迭代、高斯塞得尔迭代、 (弱、超)松弛(SOR)迭代法、收敛性判别法设有方程组 Ax=b,其中 为 n 阶非奇异阵,x,b 为 n 维列向量。ijaA)(1 简单迭代法将 Ax=b 等价变形为 x=Bx+g,其中 B 为 n 阶矩阵,g 为 n 维常向量。任取迭代格式,)0(nRxk=(0,1,2, 。 。 。 )gBxkk)()1(称为简单迭代法,B 称为迭代矩阵。特别地,若 则有所谓的雅可比迭代法),.,(0niai nijjkjiki knibxx1)()1( ,.)210)(,.21(2 高斯赛德尔迭代法把简单迭代法中的迭代矩阵 B= 写成

13、B=L+D+Unij)(其中 D=diag(b ii) ,L 和 U 分别是由 B 的对角线以下和以上的部分的元素组成的严格下、上三角阵,则高斯赛德尔迭代格式为:.)2,10(,)()1()1( kgxDxkkk其分量形式为:njikjiijkjki knib1)(1)()( ,.)210;,.(3逐次超松弛(SOR)迭代法逐次超松弛迭代法为:.)2,10()()1( ()1()( kgUxLxxkkkk其中 L 和 U 分别是雅可比迭代矩阵 的对角线以下和以上元素AadiIBi组成的严格下、上三角阵, 分 量 形 式 为 :.)(1biagi7 11)()()()1( ,.)210;,.21

14、(ij nij ikjikjikiki kniabxxaxx 称为松弛因子。4 。迭代法的收敛性判定(1)对任意初始向量 和右端项 g 简单迭代法都收敛的充要条件是)0(x .1)(B(2)若迭代矩阵 B 满足|B| 1= ,则简单迭代法和相应的高斯赛德法niijjb1|ma都收敛。(3)若迭代矩阵 B 满足|B| = ,则简单迭代法和相应的高斯赛德法niijj1|ax都收敛。(4)若系数矩阵 A 对称正定,则求解 Ax=b 的与 Jacobi 迭代法对应的高斯赛德法都收敛。(5)若系数矩阵 A 严格对角占优或 不可约弱对角占优,则求解 Ax=b 的与 Jacobi 迭代法对应的高斯赛德法都收

15、敛。(6)逐次超松弛法收敛的必要条件是 02。若系数矩阵对称正定,则 02也是充分条件。第七章 方程求根关键词:二分法、简单迭代法、不动点原理、 (方程组)牛顿迭代法、收敛阶、弦割法、抛物法、加速法、ARTEN 加速法一、内容提要本章内容主要讨论非线性方程 f(x)=0 的数值求解,其他内容未列入本书。常用的方法有:1 二分法(参阅参考文献 1。2)2 简单迭代法将 f(x)=0 变为另一种等价形式: ,选取根 x*的近似值 则递推)(x,0ba关系式 k=(0,1,2。 。 。 ) ,可产生迭代序列x k,这种迭代方法称为简单)(1kkx迭代法。简单迭代法有如下收敛性定理:设函数 在有限的区

16、间a,b 上满足(1) ,当 ;)( bxa)(,ba(2) ,1|,|)(|,1|)(| 22121 xqxbaxqx (或对则 在a,b上存在唯一的解 x*,且对 由 ,0ba8k=0,1,2。 。 。),(1kkx产生的迭代序列x k收敛于 x*,并有误差估计式|1| 0qk3 牛顿迭代法迭代公式为 ,,.21,)(1kxfxk设 f(x)在包含其零点 x*的某个闭区间 a,b内有二阶连续导数,且 则只要取0)(*xfx0 充分接近于 x*,牛顿迭代法产生的序列x k收敛于 x*,且有:)(2)(*1limfxk4 弦割法、抛物线法5 迭代收敛阶的概念设迭代过程 收敛于 的根 x*,如果

17、迭代误差 满足),(1kkx)(*xek下列关系式:(c 为常数且不为 0)epk|1li则称该迭代过程是 p 阶收敛的,p=1,01 称超线性收敛。前述的简单迭代法是线性收敛的,牛顿法是平方收敛的,双点弦割法的收敛阶为1.618,抛物线法的收敛阶是 1.839。收敛阶的判别原理:对于迭代过程 如果 在所求的根 x*的邻近连续,),(1kkx)(xp且 则称该迭代过程在点 x*邻近是 p 阶0,)(.)()( *)(*1* xxpp收敛的。第八章 矩阵的特征值与特征向量的计算关键词:圆盘定理、幂法、反幂法、原点平移法、瑞利商加速法、雅可比法、QR 分解、H 变换、HESSENNERG 矩阵、正

18、交变换法1 圆盘定理(Gerschgorin 定理)矩阵 的任何的一个特征值至少位于复平面上nijaA)(的几个圆盘nijjii izD1),.2(|,|中 的 一 个 圆 盘 上 。92、乘幂法与反幂法乘幂法是求一个实矩阵的按模最大的特征值和对应的特征向量的方法。它是一个迭代法。设 为 A 的特征值,x i 为 对应的特征向量,并且有:ii(i=1,2, 。 。 。 ,n)线性无关,则对任一非零inx及|.|21向量 ,乘幂法的计算公式为)( 0u中绝对值)()( ) 是 向 量(其 中 kkkk vvvAuvmax,.21),max(/)()(19最大的一个分量。迭代结果的各种情况讨论见参

19、考文献 1,3。特别当 ).(,)(),(/| 11)(21 kukk 时 ,反幂法和两方法的加速法见参考文献 1,2。2 实对称阵的雅可比方法雅可比方法是求一个实对称矩阵的全部特征值与特征向量的方法。起基本思想是通过正交相似变换,把实对称矩阵化为对角阵,该对角阵的对角元就是原对称矩阵的特征值。记 A0=A,施行一系列正交相似变换 的 上 半 部 分 中 选 ) , 则若 11,.2,KTkKAR1cossin1inco1 kR 行行qp其中 的选取必须满足 的 各 个 列 向 量 。 。 。变 换 矩 阵 的 特 征 向 量 是 正 交 相 似这 样 就 得 则时 , 令当 时 ,当TkT

20、kqpkqkpkpq kpqkqkpRAattc csignccaasign signa21 )()(2 2)1()1()1( )1()1()()()()( ,0.cossin,/os,| /)(ta,t.t,/2; ,2co4/.4|,/2t 4、QR 方法、Householder 变换、对称三对角阵的二分法等见参考文献 1,2。第九章 常微分方程初值问题数值解法10关键词:欧拉法、后退欧拉法、梯形法、显式法、隐式法、 (2、3、4 阶)龙格库塔法、单步法、线性多步法、ADAMS 法、预测-校正法、相容性、收敛性、稳定性(判别法) 、收敛阶、局部截断误差、全局截断误差、刚性方程本章讨论形如:

21、的初值问题的数值解法,其中 f(x,y)是已知连续函数, 为给)(,aybxfdx 定的初值,假设该初值问题的解存在唯一。(一) 单步法1 欧拉(Euler)法及其改进形式(1) 向前欧拉公式(欧拉折线法或欧拉显格式)其中 h 为步长,该方法的局部截断误01 ,.210),(ynyxhfn差阶为 O(h 2) ,是一阶方法。(2) 向后欧拉公式(后退欧拉公式)011 ,.20),(ynyxfnn(3) 梯形公式(欧拉公式的改进),.210),(),(211 nyxfyxfhnnn这也是一个隐格式,是二阶方法。实用中常按下述爹带进行求解: ,.210.;,)()(2,()1()0( nkyxfy

22、xfhyknnkn该方法是二阶方法;如果关于 k 仅迭代一步,并换 则称该方法为欧拉预,1)(nk和估校正方法。2 显式龙格-库塔方法(Runge-Kutta)方法一般形式为:其中 Ci 为待定权因子,满足miKbyhaxfKCyijjniiimin,.21)(111为所使用的 f 值的个数, a1=0,a i(i 1) ,b ij 均是待定参数,随着 m 取miC1,11不同的正整数值,便可以得到各阶显式龙格-库塔格式,特别当 m=4 时,可以得到四阶、经典的龙格-库塔格式:)3423121 4321,(),()(6KyhxfKfyxhKnnn3 单步法的收敛性与稳定性4 单步法的收敛性与稳

23、定性的概念非常重要,有关定义和结论可见参考 2,6,10。(二) 线性多步法1 一般公式其中. 0101 nknknnknkn fffhyaya ).2,10(,ii为和 n 无关的常数, 若 则称上式是 k 步方法;若),(iiyxf ,0|,0k局部截断误差为 O(h s+1) ,则称上式是 s 阶方法。2、阿达姆斯显式公式:可由数表给出,见参考文献rkknknnfy01 系 数1,2,10。特别当 r=0 时,上式是欧拉折线法:r=3 时,上式为24/93759 3211 nnn fffhy阿达姆斯显式公式的局部截断误差为这里 特别)(212)1(, nrrhryR ,1)( 101 nrnrr xds当 r=3 时, .705)()1(,3nh2 阿达姆斯隐式公式:,其中 也可查表而得,当 r=0 时上式rkknnfhy01*1*rk为欧拉后退公式;r=1 时为梯形公式;r=3 时为:其局部截断误差是:24/59111 nnnn fffhy12).(2*12)(, nrrhryR当 r=3 时 709)5()(,3线性多步法的收敛性、稳定性和预估校正法等见参考文献 10。

展开阅读全文
相关资源
猜你喜欢
相关搜索
资源标签

当前位置:首页 > 企业管理 > 管理学资料

本站链接:文库   一言   我酷   合作


客服QQ:2549714901微博号:道客多多官方知乎号:道客多多

经营许可证编号: 粤ICP备2021046453号世界地图

道客多多©版权所有2020-2025营业执照举报