数值分析第24讲(ExCht1-9).ppt

,第24讲 第1-9章 习题课,一、基本内容及基本要求,第一章、绪论 了解数值分析的研究对象与特点。 了解误差来源与分类,会求有效数字; 会简单误差估计。 了解误差的定性分析及避免误差危害。,第1-3章 习题课 (绪论、插值、逼近),第二章、插值法 了解插值的概念。 掌握拉格朗日(Lagrange)插值法及其余项公式。 了解均差的概念及基本性质，掌握牛顿插值法。 了解差分的概念，会牛顿前插公式、后插公式。 会埃尔米特(Hermite)插值及其余项公式。 知道高次插值的病态性质,会分段线性插值和分段埃尔米特插值及其误差和收敛性。 会三次样条插值,知道其误差和收敛性。,第三章、函数逼近与曲线拟合 了解函数逼近的基本概念,了解范数和内积空间。 了解正交多项式的概念,了解切比雪夫多项式和勒让德多项式以及它们的性质,知道其他常用正交多项式。 理解最佳一致逼近的概念和切比雪夫定理,掌握最佳一次一致逼近多项式的求法。 理解最佳平方逼近的概念,掌握最佳平方逼近多项式的求法,了解用正交多项式做最佳平方逼近的方法。 了解曲线拟合的最小二乘法并会计算,了解用正交多项式做最小二乘拟合。 了解最小二乘三角逼近与快速傅里叶变换*。,二、练习,P19, 5,9.,P59, 6,8.,7、P59, 4.,P115,1,4(2),6,8,13,15,17(1),19，按基本方法即可， 22(选作).,一、数值积分与数值微分,第4-5章 习题课 (数值积分和数值微分,解线性方程组的直接法),基本内容及基本要求,了解数值求积的基本思想、代数精度的概念、插值型求积公式及其代数精度、求积公式的收敛性和稳定性。 掌握牛顿-柯特斯公式及其性质和余项。 掌握复化梯形公式和复化辛普森公式及其余项。 了解龙贝格(Romberg)求积算法，知道外推法。 会高斯求积公式,了解高斯-勒让德求积公式和高斯-切比雪夫求积公式。 了解几种常用的数值微分方法。,二、练习,三、解线性方程组的直接方法,基本内容及基本要求 了解求解方程组的两类方法,了解矩阵基础知识。 掌握高斯消去法,会矩阵的三角分解。 掌握高斯列主元素消去法,了解高斯-若当消去法。 掌握直接三角分解法,了解平方根法,会追赶法,了解有关结论。 了解向量和矩阵的几种范数。 了解矩阵和方程组的性态,会求其条件数。 会初等反射阵和平面旋转阵,了解QR分解,了解用正交约化法解超定方程组。,分别用顺序Gauss消去法和直接三角分解法(杜利脱尔分解)求解线性方程组,四、练习,2. 设A为n阶对称正定阵，试证: (1) A的对角元素aii0; (2) 设L为非奇异阵，则LALT是对称正定阵； (3) 经顺序Gauss消去法A化为,求证A2为对称正定.,证明: (1)由正定二次型理论,aii=eiAei0. (或因所有主子式0),(2) 因(LALT)T= LALT , 故LALT是对称的; 又因对于任意x0,则有y=LTx0,从而 xTLALTx=(LTx)TA(LTx)=yTAy 0, 故LALT是对称正定阵.,(3) 经顺序Gauss消去法A化为,A2是对称的,因为,A2是正定的,这是因为经顺序Gauss消去法A的各阶顺序主子式的值不变, a11A2的k阶顺序主子式=A的k+1阶顺序主子式0, 且a110, 于是得出A2的各阶顺序主子式0.,3. 用带行交换的杜利脱尔分解计算线性代数方程组,AX=b，其中,4. 用追赶法求解三对角方程组,一、解线性方程组的迭代法,第6-8章 习题课 (线性方程组迭代解法,解非线性方程,矩阵特征值),基本内容及基本要求,了解迭代法及其收敛性的概念。 掌握雅可比(Jacobi)迭代法、高斯-赛德尔(Gauss-Seidel)迭代法和超松弛(SOR)迭代法。 3. 了解一阶定常迭代法的基本定理，掌握特殊方程组迭代法的收敛条件。 4. 知道分块迭代法。,雅可比迭代法计算公式：对k=0,1,高斯塞德尔迭代法计算公式：对k=0,1,SOR迭代法的计算公式:对k=0,1,定理7 若矩阵A按行(或列)严格对角占优,或按行(或列)弱对角占优不可约；则Jacobi迭代、Gauss-Seidel迭代都收敛。,定理9 对于线性方程组Ax=b，若A为对称正定矩阵，则当02时，SOR迭代收敛.,定理10 对于线性代数方程组Ax=b， 若A按行(或列)严格对角占优，或按行(或列)弱对角占优不可约；则当0w1时，SOR迭代收敛。,二、非线性方程求根,基本内容及基本要求,了解求根问题和二分法。 了解不动点迭代法,及不动点存在性和迭代收 敛性； 了解收敛阶的概念和有关结论。 3. 了解加速迭代收敛的埃特金方法和斯蒂芬森方法。 4. 掌握牛顿法及其收敛性、了解简化牛顿法和牛顿法 下山法,了解重根情形。 5. 掌握弦截法，了解抛物线法。 6. 了解非线性方程组的迭代解法。,三、矩阵特征值问题计算,了解特征值和特征向量的概念和性质, 了解圆盘定理、Schur定理和Rayleigh商。 2. 掌握乘幂法，了解其加速收敛技术,会反幂法。 3. 了解豪斯霍尔德方法。 4. 了解QR方法。,基本内容及基本要求,反幂法计算公式:,v1=(10, 8,1)T, mu1=10 u1=(1.0000, 0.8000, 0.1000)T v2=(7.2000,5.4000, -0.8000)T mu2=7.2000 u2=(1.0000, 0.7500, -0.1111)T v3=(6.5000, 4.7500, -1.2222)T mu3=6.5000 u3= (1.0000, 0.7308, -0.1880)T,第九章 常微分方程初值问题数值解法 关键词：欧拉法、后退欧拉法、梯形法、显式法、隐式法、（2、3、4阶）龙格库塔法、单步法、线性多步法、ADAMS法、预测-校正法、相容性、收敛性、稳定性（判别法）、收敛阶、局部截断误差、全局截断误差、刚性方程,本章讨论形如,的初值问题的数值解法，其中f（x,y）是已知连续函数，,为给定的初值，假设该初值问题的解存在唯一。,单步法 欧拉（Euler）法及其改进形式 向前欧拉公式（欧拉折线法或欧拉显格式）,其中h为步长，该方法的局部截断误差阶为O（h2）， 是一阶方法。 2.向后欧拉公式（后退欧拉公式）,3. 梯形公式（欧拉公式的改进）,这是二阶隐格式方法。实用中常按下述爹带进行求解：,该方法是二阶方法；如果关于k仅迭代一步，并换,则称该方法为欧拉预估校正方法。,2.显式龙格-库塔方法（Runge-Kutta）方法,一般形式为：,其中Ci为待定权因子，满足,m为所使用的f值的个数，a1=0，ai（i1），bij均