计算方法大作业.docVIP

作业一： 线性方程组的数值解法 两类方法，第一是直接解法，得到其精确解； 第二是迭代解法，得到其近似解。 一、Gauss消去法 1．顺序Gauss消去法 记方程组为： 　　　 消元过程： 经n－1 第k步　 若 　　 回代过程： 2．Gauss—Jordan消去法 避免回代，消元时上下同时消元 3．Gauss列主元消去法 Gauss列主元消去法原理：每步消元前，选列主元，交换方程。 算法：将方程组用增广矩阵表示。 消元过程： 对k=1,2,n-1, 选主元，找 如果，则矩阵A奇异，程序结束；否则执行3。 如果，则交换第k行与第行对应的元素位置， 消元，对i=k+1, ,n,计算 对j=L+1, ,n+1,计算 回代过程： 1．若则矩阵A奇异，程序结束；否则执行。 2． 二、LU分解法 1．求解原理 线性方程组写成矩阵形式为： AX=b 若A=LU，则LUX= b， 记UX=Y 则LY= b 若L、U为特殊矩阵，则求解线性方程组变为解两个特殊线性方程组问题。 2．Doolittle分解 L为下三角矩阵, U为上三角矩阵,不一定能分解,分解也不一定唯一; 设L或U是单位三角矩阵, 若能分解,则可分解唯一. L是单位下三角矩阵,称为Doolittle分解; U是单位上三角矩阵,称为Crout分解； 定理: n阶矩阵A有唯一分解的充要条件为A的前n-1阶主子式都不为0. Doolittle分解算法： 由矩阵乘法： 得到： 算法特点：先计算U的行，再计算L的列，交替进行；存储时可用紧凑格式。 矩阵分解后，解两个三角方程组： LY= b，UX=Y 3．Crout分解 若L为下三角矩阵，U是单位上三角矩阵，则称Crout分解； 算法特点：先计算L的列，再计算U的行，交替进行。 线性方程组的迭代解法 一、Jacobi迭代公式 一般地，对线性方程组 若，则可从第i个方程中解出,得到Jacobi迭代公式： 简记为： Gauss--Seidel迭代公式 非线性方程的数值解法 二分法 1．根的存在性 定理：函数f(x)在区间[a，b]连续，且f(a).f(b)0,则方程f(x)=0在区间[a，b]有根。 方程的根存在，不一定唯一，若在区间[a，b]上有唯一根，称区间[a，b]为根隔离区间。 二分法（区间逐次分半法） 原理：通过计算根隔离区间中点，将区间分半，缩小区间，得到方程近似根数列{。 取 二、迭代法 1．迭代原理 迭代法是一种逐次逼近法，由提供的递推公式计算，逐次精确，直到满足精度要求。 方程f(x)=0变形为， 得到递推公式--------简单迭代公式 称为迭代函数 给初值计算，得到数列{，若，则称迭代收敛，否则发散。 2．迭代收敛性判定 收敛性定理：设方程的迭代函数在[a，b]满足： （1）当时，； （2）在[a，b]可导，且，， 则（1）方程在[a，b]有唯一根； （2）迭代公式收敛，即； （3）误差估计。 说明可根据迭代函数的导数判断迭代收敛性。 3．Newton 迭代法 Newton切线公式： 几何作法 迭代公式 例：利用解二次方程推导近似计算的公式。 由Newton切线公式 Newton弦截公式： Newton切线公式的缺点及改进 几何作法 迭代公式 Newton弦截公式是两步公式。 插值与拟合 1．拉格朗日插值 设已知 x0,x1,··· ,xn 及 yi = f(xi)(i = 0,1,··· ,n)，Ln(x) 为不超过 n 次的多项式，且满足插值条件 Ln(xi) = yi (i = 0,1,··· ,n)。由对 L2(x) 的构造经验，可设 , 其中，li(x)(i = 0,1,··· ,n) 均为 n 次多项式且满足 (4.9)(i,j = 0,1,··· ,n)。不难验证，这样构造出的 Ln(x) 满足插值条件。因此问题归结为求 li(x)(i = 0,1,··· ,n) 的表达式。因 xj(j 6= i) 是 n 次多项式 li(x) 的 n 个根，故可设 再由 得 故有 称为 n 阶拉格朗日插值公式，其中 li(x)(i = 0,1,··· ,n) 称为 n 阶拉格朗日插值的基函数。 2Hermite 插值 拉格朗日插值仅考虑节点的函数值约束，而一些插值问题还需要在某些节点具有插值函数与被插值函数的导数值的一致性。具有节点的导数值约束的插值称为 Hermite 插值。 3牛顿插值法 由于拉格朗日插值公式计算缺少递推关系，每次新增加节点需要重新计算，高次插值无法利用低次插值的结果。通过