讲：ch直接三角分解法.pptVIP

通过下列运算可看出其置换的情况 a1 = {{1, 0, 0}, {-1, 1, 0}, {-3, 14/3, 1}}; MatrixForm[%] a2 = {{-1, 4, 3}, {0, 3, 6}, {0, 0, -14}}; MatrixForm[%] a1.a2 MatrixForm[%] 三、直接三角分解法 提示 LU分解补充练习 Clear[P1,P2,P3,L,U,A,B] A={{4,2,1,5},{8,7,2,10},{4,8,3,6},{12,6,11,20}}; MatrixForm[%]; P1={{1,0,0,0},{-A[[2,1]]/A[[1,1]],1,0,0},{-A[[3,1]]/A[[1,1]],0,1,0}, {- A[[4,1]]/A[[1,1]],0,0,1}}; MatrixForm[%]; A1=P1.A; MatrixForm[%]; P2={{1,0,0,0},{0,1,0,0},{0,-A1[[3,2]]/A1[[2,2]],1,0},{0,-A1[[4,2]]/A1[[2,2]],0,1}}; MatrixForm[%]; A2=P2.A1; MatrixForm[%]; P3={{1,0,0,0},{0,1,0,0},{0,0,1,0},{0,0,-A2[[4,3]]/A2[[3,3]],1}}; L