第5章_解线性方程组的直接法计算方法答辩.pptVIP

当A进行LU分解后,Ax=b就容易解了. 即Ax=b等价于: 结束 这样,Ax=b分解为解两个三角形方程组,三角形方程组是极易求解的.(5.18)即是高斯消元法的矩阵表述. * * L,U 的求法： 1.U的第一行元素 得到： * 2. L的第一列元素 得到： ……，若已计算出 U的前k-1行和L的前k-1列元素，则 * 3. U的第k行元素：ukk, uk,k+1, …, ukn 当列标j=行标k时，有 得到： (5.19) * 4. L的第k列元素：lk+1,k, lk+2,k, …, ln,k 当行标i=列标k时，有 得到： (5.20) 1.解 L y = b 2.解 U x = y 结束 * 杜利特算法实际上就是高斯消元法的另一种形式.它的计算量与高斯消元法一样.但它不是逐次对A进行变换,而是一次性地算出L和U的元素.L和U的元素算出后,不必另辟存贮单元存放,可直接存放在A的对应元素的位置,节省存贮单元,因此也称为紧凑格式法. AX=b 的解： 杜利特算法可用文字描述如下: 将A的元素划分为形如“┏”的n框, 如A的第一行和第一列元素为第一框,紧靠第一框内部的第二行和第二列元素为第二框,依此类推,ann一个元素为第n框.为便于描述我们称原A矩阵中元素为a元素,计算后每框中的每行元素为u元素(包含对角元素),每列元素为l元素(不包含对角元素). 这样我们可以在A上直接修改计算,作如下的操作: (1)第1框:u元素等于对应的a元素;l元素等于对应的a元素除以本框对角元素; (2)第2到n框:u元素等于对应的a元素减去已经算过的各框中同行同列元素的乘积;l元素除进行以上操作外,还要除以本框对角元素 结束 * 例 用Dolittle分解求解方程组 结束 * 得到的矩阵由L和U的元素拼合而成,它的上三角部分(包含对角元素)即是LU分解后的U矩阵;它的下三角部分(不包含对角元素)加上一个单位阵(即在对角元处全部写为1),就是L矩阵. 将LU分解换一个提法:要求L为一般下三角阵,U为单位上三角阵,就是克劳特分解. LT显然是单位上三角阵,记 结束 这就是克劳特分解。 有 * 实际上将A转置为AT,进行LU分解 AT=LU 则 A=UTLT 显然当A的各阶顺序主子式非零时,它是存在且唯一的. UT显然是一般下三角阵,记 5.2.2 克劳特(Crout)分解 结束 * 5.2.3 追赶法 在很多情况下,如三次样条插值,常微分方程的边值问题等都归结为求解系数矩阵为对角占优的三对角方程组Ax=f,即: 其中|i-j|1时,aij=0,且满足如下的对角占优条件: (1)|a1||b1|0,|an||cn|0 (2)|ai|≥|bi|+|ci|, bici≠0, i=2,3,…,n-1. * 其系数矩阵A进行Crout分解, 得: 比较A=LU两边的系数, 可得: 且： 结束 * 用矩阵乘法比较和Crout分解可推导总结算法步骤如下: 1.分解计算 2.解Ly=f 3.解Ux=y 结束 * 实际计算中Ax=f的阶数往往很高,应注意A的存贮技术.已知数据只用4个一维数组就可存完.即{ai},{bi},{ci},{fi}各占一个一维数组, {ui}和{vi}可存放在{bi},{ci}的位置,{yi} 和{xi}则可放在{fi}的位置,整个运算可在4个一维数组中运行.追赶法的计算量很小,只是5(n-1)次乘除法.追赶法的计算也不要选主元素. 例(见教材) 其中： 结束 * 5.2.4 对称正定矩阵的LDLT分解 若A对称正定矩阵，对A作Dolittle分解，可得： 计 算 方 法 课 件 第5章 解线性方程组的直接法 在自然科学和工程技术中,很多问题归结为解线性方程组.有的问题的数学模型中虽不直接表现为含线性方程组,但它的数值解法中将问题“离散化”或“线性化”为线性方程组.因此线性方程组的求解是数值分析课程中最基本的内容之一. 线性方程组: 结束 常记为矩阵形式 Ax=b (5.2) * 此时A是一个n×n方阵,x和b是n维列向量. 根据线性代数知识若 |A| ≠0,(5.2)的解存在且唯一. 关于线性方程组的解法一般分为两大类,一类是直接法,即经过有限次的算术运算,可以求得(5.1)的精确解(假定计算过程没有舍入误差).如线性代数课程中提到的克莱姆算法就是一种直接法.但该法对高阶方程组计算量太大,不是一种实用的算法.实用的直接法中具有代表性的算法是高斯消元法,其它算法都是它的变形和应用. 另一类是迭代法,它将(5.1)变形为某种迭代公