南京航空航天大学计算方法期末考试.ppt

《南京航空航天大学计算方法期末考试.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第二章非线性方程的数值解法常用方法1二分法2一般迭代法3牛顿迭代法4弦截法根的隔离；误差估计；迭代收敛阶2一般迭代法(1)迭代法(1)把（1）等价变换为如下形式(2)建立迭代格式(3)适当选取初始值x0，递推计算出所需的解。定理2.2（非局部收敛定理）如果在上连续可微且以下条件满足：命题2.2若在区间内，则对任何，迭代格式不收敛。推论设x*=g(x*)，若g(x)在x*附近连续可微且，则迭代格式xk+1=g(xk)在x*附近局部收敛。(2)迭代法的收敛性简单地代之以(3)迭代法的误差估计3牛顿迭代法其迭代函数为牛顿迭代法4弦截法弦截法第三章线性代数方程组的数值解法解线性方程组的消

2、去法解线性方程组的矩阵分解法3解线性方程组的迭代法给定一个线性方程组求解向量x。(1)高斯消去法1.解线性方程组的消去法1）消元过程：对k=1,2,…,n依次计算2)回代过程：这一无回代的消去法称为高斯-若当（Jordan）消去法(2)高斯-若当（Jordan）消去法高斯-若当（Jordan）消去法一般公式：定理3.1如果的各阶顺序主子式均不为零，即有即消去法可行。推论若系数矩阵严格对角占优，即有(3)选主元素的消去法主元素的选取通常采用两种方法：一种是全主元消去法；另一种是列主元消去法。2解线性方程组的矩阵分解法一、非对称矩阵的三角分解法矩阵分解法的基本思想是：可逆下三角矩阵可

3、逆上三角矩阵对于给定的线性方程组（1）分解——解两个三角形方程组。矩阵的Crout分解的计算公式(3-12)注：3.3.3对称正定矩阵的三角分解定义3.1若n阶方矩阵A具有性质且对任何n维向量成立，则称A为对称正定矩阵。定理3.4若A为对称正定矩阵，则(1)A的k阶顺序主子式(2)有且仅有一个单位下三角矩阵L和对角矩阵D使得（3-16）这称为矩阵的乔里斯基（Cholesky）分解。(3)有且仅有一个下三角矩阵，使（3-17）这称为分解矩阵的平方根法。3解线性方程组的迭代法迭代法思想:(1)Ax=b(3-1)(2)建立迭代格式这称为一阶定常迭代格式，M称为迭代矩阵。约化便得从而可建

4、立迭代格式对（3-23）以分量表示即(1)、Jacob迭代法雅可比（Jacobi）迭代则雅可比迭代格式（3-24）可用矩阵表示为MJfJ用矩阵表示为对雅可比迭代格式修改得高斯-塞德尔（Gauss-Seidel）迭代fG-SMG-S(2)Gauss-Seidel迭代法例3.10分别用雅可比迭代法和高斯-塞德尔迭代法求解线性方程组解相应的迭代公式为雅可比迭代高斯-塞德尔迭代令取四位小数迭代计算由雅可比迭代得由高斯-塞德尔迭代得定理3.5若一阶定常迭代格式（3-26）的迭代矩阵满足条件则该迭代格式对任何初始向量均收敛。则该迭代格式对任何初始向量均收敛。定理3.6若一阶定常迭代格式（3-

5、26）的迭代矩阵满足条件迭代法的收敛性推论如果线性代数方程组Ax=b的系数矩阵A为严格对角占优矩阵，即则相应的雅可比迭代法与高斯-塞德尔迭代法对任何初始向量均收敛。定理3.8一阶定常迭代格式对任何初始向量均收敛的充分必要条件为其迭代矩阵的谱半径小于1，即这里为M的特征值第四章函数的插值与拟合法1插值多项式的构造2最小二乘法定义4.1设y=f(x)在区间[a,b]上连续，在[a,b]内n+1个互不相同的点上取值.求一代数多项式P(x)，使得则称P(x)为f(x)的插值函数1插值多项式定理4.1在n+1个互异点上满足插值条件(4-1)的次数不超过n次的插值多项式存在且惟一。两种插值多

6、项式形式(1)拉格朗日插值多项式下列列表函数的多项式Ln(x)xx0x1---xi-1xixi+1---xnyy0y1---yi-1yiyi+1---yn线性插值(n=1)，抛物插值(n=2)(2)牛顿均差插值多项式Ln(x)和N(x)插值多项式的余项例：已知列表函数,并计算f(0.5)的计算值。解：x-1012y111-5x-1012y111-5(1)由数据表构造均差表kxkf(xk)一阶均差二阶均差三阶均差0-1100-11010-3211-632-5又解：2数据的多项式最小二乘拟合xx0x1---xi-1xixi+1---xnyy0y1---yi-1yiyi+1---yn已

7、知一组数据：------这个多项式称为这组数据的最小二乘拟合多项式求最小二乘多项式的步骤例4.5试对以下数据进行多项式拟合xi12345678yi1.13.88.715.624.637.449.664.2解