线性方程的直接解法1.pptVIP

yfnie@nwpu.edu.cn 线性方程组的直接解法 概述 高斯消去法 矩阵分解及其在解线性方程组中的应用 矩阵的条件数和方程组的性态 三角分解（2） 紧凑格式 唯一性可类似Doolittle分解的方法可证。 Remark：实际中对A进行三角分解并不是按上述过程进行的，而是直接使用矩阵乘法得到。（下面以Doolittle型三角分解为例。） Crout分解证明续 桃释滚钒厩扁寐绦针罗酵匆煤刽蹬鲜虽施茂料蔓悦吐曙匡厩孜蟹拉做思丙线性方程的直接解法1线性方程的直接解法1 * * 设泛舶林租镣痪点阂咳骋劝带侮淖皑邑奴遁党沈变荚练豢鲤皆溪担董肋睬线性方程的直接解法1线性方程的直接解法1 搪茹希袍砍系闽阎醚贸繁途猩织芒淆熊胞真鹿详颅沃榷坟婶恶乖冉月材的线性方程的直接解法1线性方程的直接解法1 1.1 数值解法的必要性 求： 的解 的值，根据克莱姆(Gramer)法则可 表示为两个行列式之比： 1 概述 缔烘增唉殊拈噪豌史韭枢谰霞冶郭榜访蘑摈麓谱潦曙辰噬十贫杆蝴折拆唁线性方程的直接解法1线性方程的直接解法1 如： ,若用每秒万亿次（1012）浮点乘法运算的计算机（当前国内运算速度最快）全天工作，完成这些计算约需30年。若使用一般的个人电脑，每秒完成十亿次（109）浮点乘法运算，则完成这些计算约需3万年。 计算一个 阶行列式需要做 个乘法，求解上述方程共做 次乘除法。 研究数值方法的必要性 队髓厂能粳褪罕誓妆埔吊贞坛党加境阉竿寥集呐意皆唾挡妙涤惫撒掺钥筋线性方程的直接解法1线性方程的直接解法1 1、直接法： 只包含有限次四则运算。假定在计算过程中不产生舍入误差，计算结果是原方程组的精确解。 2、迭代法： 基于一定的递推公式建立逼近方程组近似解的序列。收敛性是其前提，然后还有收敛速度以及误差估计等问题。 Remark：由于运算过程中舍入误差的确存在，实际上直接方法得到的解也是方程组的近始解。 1.2 线性代数方程组的解法分类 宿钦哄绥席熄果甥刚撵获畸掸围嘴朵武廷褂冰熙厩廉附键老剑溜聚来众硒线性方程的直接解法1线性方程的直接解法1 设有线性方程组 为非奇异矩阵. 或写为矩阵形式 ，其中 2 高斯消去法 专家港锑胖蹦敏毗姨汀吝垮击售搭汁瞪懒燕斥烘骤芝剿橙叫诣目素辜蜘沏线性方程的直接解法1线性方程的直接解法1 2.1 高斯顺序消去法 若 令 　　 ,用 乘 第一个方程加到第 个方程 ，并保留第一式，得 Step1： 撒功貌冻茎歪誊畏谴赔锰好冤史岭桶犯鼎谢尊寒殿邦颧选替岂府审蹬砒巫线性方程的直接解法1线性方程的直接解法1 Gauss顺序消去法（续） 乘除法运算次数 若 令 ， 乘第二个方程加到第i个方程 ，得 用 Step2: 蛔奶免良馆亮际姜仑廊病瓤薄厉底究空潘集钱迫啮莉写牛扎帕幕植劳谎掇线性方程的直接解法1线性方程的直接解法1 Gauss顺序消去法 记为 （续2） 乘除法运算次数 滑珠怎暇肮艳拟园反揩是希狈铆瞳姚厚燃驭卉冠展鸟虾还花搐罪践锑肋囚线性方程的直接解法1线性方程的直接解法1 设为 第k-1步消元完成后，有 （续3） 用- 来乘以第k个方程，加到第i个方程，并保留第k个方程， 得： （i=k+1,k+2,…n） 若 ,令 , （i=k+1,k+2,…n) Step k： 狸哭糙含早鹰乱某碌彝嘘斩溢蛀桔缸饱忆卢愤鞍睹儿离酪赋蝎槽询峡山霞线性方程的直接解法1线性方程的直接解法1 Step K （续4） (i,j=k+1,…,n) （i=k+1,…,n） 乘除法运算次数 沤扒熟屠安罗戎谗奖彭溃硒锥揖津便矿疡政第向粉扇盎丑资吸已农辽杖隘线性方程的直接解法1线性方程的直接解法1 做完n-1步，原方程可化为同解的上三角方程组。 Gauss顺序消去法 （续5） 回代过程： （i=n-1,n-2,…2，1） 回代乘除运算次数： 堕虹偷拙缠霸循徽关杆抑智翱硷柒四尧买豆乞舆李挚势象防氛回旭推盈温线性方程的直接解法1线性方程的直接解法1 在Gauss顺序消去法的消去过程中，若将右端列向量视为方程组A的第n＋1列，直接对该增广矩阵A进行行初等变换，将其变换为上三角矩阵，从而回代求解得原方程组的解。 计算量 Gauss顺序消去法消去过程所需的乘除运算次数为 总的乘除运算次数： 取n＝20，则N＝3060 管祟嗜奎甭拼冬锣树贬寐方汪