数学建模线性方程组的数值解法.pptVIP

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 线性方程组数值解法的MATLAB实现 若A为可逆方阵，输出原方程的解x 若A为n?m矩阵（nm），且ATA可逆，输出原方程的最小二乘解x 线性方程组数值解法的MATLAB实现 [x,y,p]=lu(A) 若A可逆，输出x为单位下三角阵L， y为上三角阵U，p为一交换阵P，使PA=LU. u =chol(A) 对正定对称矩阵A的Cholesky分解，输出u为上三角阵U，使A=UTU 2. 矩阵LU分解 [x,y]=lu(A) 若A可逆且顺序主子式不为零， 输出x为单位下三角阵L，y为上三角阵U，使A=LU； 若A可逆，x为一交换阵与单位下三角阵之积. 例. 解 A=[10 3 1;2 -10 3;1 3 10], b=[14 -5 14], x=A\b, [L1,U1]=lu(A); L1,U1, A1=L1*U1, [L2,U2,P]=lu(A); L2,U2,P, A2=L2*U2, A3=inv(P)*A2 并对系数矩阵作LU分解 shiyan51 若第1个方程改为 3x2+x3=14 结果如何 当n很大时Hilbert矩阵呈病态 线性方程组数值解法的MATLAB实现 3. 范数 条