第3讲 线性方程组直接求解.pptVIP

* 3.7 LU分解法 对方阵A进行LU分解A=LU的通项公式为： 三、对方阵进行LU分解的过程 = － ，（j=i,i+1,…,n） =( － )/ ，（i=j+1,j+2,…,n） 在用通项公式对方阵A进行LU分解时，一个较好的计算次序是： ① U第1行， ② L第1列， ③ U第2行， ④ L第2列， ⑤ U第3行， ⑥ L第3列， …… * 3.7 LU分解法 ⑴ 对系数矩阵A进行LU分解，得到L和U。 四、LU分解法求解线性方程组的过程 方程组Ax=b?(LU)x=b?L(Ux)=b。 ⑵ 令向量y=Ux，代入上式得：Ly=b。 回代，求解Ly=b，得到y。 通项公式为：yi=bi－ yk，i=1,2,……,n。 ⑶ 回代，求解Ux=y，得到解向量x。 通项公式为：xi=(yi－ xk)/uii，i=n,n-1,……,1。 * 3.7 LU分解法 五、LU分解法的算法。 输入原方程组的阶数n。 输入原方程组的系数矩阵a[n][n]，右端向量x[n]。 for(k=0;k=n-2;k++) for(i=k+1;i=n-1;i++)，循环1次求出L的1个元素。 s=0;