计算物理与lecture1104 .ppt

例如.求解方程组 Doolittle分解 若矩阵A有分解：A=LU，其中L为单位下三角阵，U为上三角阵，则称该分解为Doolittle分解，可以证明，当A的各阶顺序主子式均不为零时，Doolittle分解可以实现并且唯一。 A的各阶顺序主子式均不为零，即 Doolittle分解 Doolittle分解 Doolittle分解 Doolittle分解 例题 例题 ú ú ú ú ú ú ? ù ê ê ê ê ê ê ? é = - - - n n n n n b a c b a c b a c b A 1 1 1 2 2 2 1 1 O O O ú ú ú ú ú ú ? ù ê ê ê ê ê ê ? é - - - n n n n n b a c b a c b a c b 1 1 1 2 2 2 1 1 O O O 三对角方程组的追赶法 例 用追赶法解三对角线方程组 u1=1/2, game1=0 game2=-1/(2-0)=-1/2 u2=(0+1/2)/(2-1/2)=1/3 game3=-1/(2-1/2)=-2/3 u3=1/3/(2-2/3)=1/4 Game4=-1/(2-2/3)=3/4 u4=(1+1/4)/(2+3/4)=5/11 u3=1/4-15/44=-1/11 u2=1/3-2/33=3/11 u1=1/2+3/22=7/11 Finally, an important property of hermitian and symmetric matrices is that they have real eigenvalues. Real orthogonal matrices is 矩阵特征值和特征向量的数值解法 幂法 逆幂法 古典Jacobi法 Jacobi法的改进 QR算法 QR算法是目前求一般矩阵全部特征值和特征向量行之有效的一种方法，它适合于对称矩阵，也适合于非对称矩阵。QR算法最早在1961年由J.G.Francis提出的，后来经一系列数学家进行深入讨论,并作出了做有成效的改进与发展。 QR算法 矩阵的QR分解 定理: 设矩阵A非奇异，则一定存在正交矩阵Q，上三角矩阵R，使 且当要求R的主对角元素均为正数时，则上述分解式是唯一的。 QR算法是对A进行一系列的正交相似变换，达到求出矩阵A的全部特征值和相应的特征向量。 将矩阵A进行LU分解 Y1=12 Y2=7 Y3=-2 x3=2, x2=1, x1=3 将矩阵A进行LU分解 Ax=w Ax=BCx=w By=w y=Cx This equation can be calculated in two steps For our four-dimentional example this takes the form By=w y=Cx forward substitution Back substitution ú ú ú ? ù ê ê ê ? é - = ú ú ú ? ù ê ê ê ? é ú ú ú ? ù ê ê ê ? é - - - - - 30 19 10 6 3 6 19 13 4 6 5 2 3 2 1 x x x . LU . 1 分解求解方程组 试用 例 所以方程组的解为 。 T x ) 1 , 2 , 3 ( = 解得： x x x 3 , 2 , 1 1 2 3 = = = x x x 4 1 10 4 7 3 6 5 2 3 2 1 ú ú ú ? ù ê ê ê ? é - = ú ú ú ? ù ê ê ê ? é ú ú ú ? ù ê ê ê ? é - - 解 y Ux ) 2 ( = 5.4 The determinant and Inverse of a matrix 矩阵的行列式和逆矩阵 定义矩阵A的行列式 . 1 The determinant of a matrix = The basic definition of the determinant of A is 1 The determinant of a matrix The basic definition of the determinant of A is L= 1 ... ... ... ... 1 1 1 1 2 1 32 31 21 l l l l l l nn n n ú ú ú ú ú ú ? ù ê ê ê ê ê ê ? é - O u u u u u u nn n n ... ... 2 22 1 12 11 ú ú ú ú ? ù ê ê ê ê ? é M O U= A=L?U For the sake of simplicity, let us look at a (4x4