NA2方程组.pptVIP

根据矩阵乘法及相等的定义,有 得公式 u1j=a1j j=1,2,…,n li1=ai1 / u11 i=2,3,…,n 2.3 Matrix Factorization 2.3 Matrix Factorization 由LU=A及对L和U的要求可以得到分解的计算公式．根据下式(Crout分解)： l11 l21 l22 l31 l32 l33 … … … ln1 ln2 … lnn-1 lnn 1 u12 u13 … u1n 1 u23 … u2n … 1 u(n-1)n 1 = … … … … a ann L U n1 an3 … a11 a12 a13 … a1n a21 a22 a23 … a2n a31 a32 a33 … a3n an2 A … … … 2.3 Matrix Factorization 在计算机程序中常常用这种方法解线性代数方程组。 它的优点是存储量很省。L和U中的三角零元素都不 必存储，就是U的对角元素也因为都是1没有必要再 记录在程序中，这样只用一个n阶方阵就可以把L和 U贮存起来。即：下三角（包括对角元）存储L各元 素 而上三角存储U的元素。 再考察公式S会发现A中任一元素aij只在计算lij（j=i） 和uij（ji)中用到一次以后就不再出现了，因而完全 可以利用原始数组A的单元，一个个逐次贮存L或U中 的相应元素，即： a11 a12 a13 … a1n u11 u12 u13 … u1n a21 a22 a23 … a2n l21 u22 u23 … u2n a31 a32 a33 … a3n l31 l32 u33 … u3n an1 an2 an3 … ann ln1 ln2 ln3 … unn … … … ... … … … … ... ... (1) (3) (5) (2n-1) (2) (4) (6) (2n) * 引言 为便于以下讨论,把涉及到的有关名词及问题的引出暂介绍如下: 如果未知量的个数为 n ,而且关于这些未知量x1,x2, … ,xn 的幂次都是一次的(线性的),那末, n 个方程 a11x1+a12x2+ … +a1nxn=b1 ┆ ┆ ┆ (1) an1x1+an2x2+ … +annxn=bn 构成一个含 n 个未知量的线性方程组,称为 n 阶线性方程组 其中,系数a11,…,a1n,a21, …,a2n, …,an1, …,ann 是给定的常数;b1, …,bn 也是给定的常数,通常称为常数项,或称为方程组的右端. 方程组(1)也常用矩阵的形式表示,写为 Ax=b 第二章 解线性方程组的直接法 /* Direct Method for Solving Linear Systems */ 其中,A是由系数按次序排列构成的一个n阶矩阵, 称为方程组的系数矩阵,x和b都是n维向量,称为方程 组的右端向量 . 使方程组(1)中每一个方程都成立的一组数x