MATLAB实用教程 教学课件 作者 张磊 郭莲英 丛滨 02.pptVIP

* 考虑线性方程组AX=B和矩阵A的LU分解，线性方程组可以改写成L*U*X=B，由于左除算符\可以快速处理三角矩阵，因此可以快速解出： X=U\(L\B) 矩阵的行列式和逆也可以利用LU分解来计算，如 det(A)=det(L)*det(U) inv(A)=inv(U)*inv(L) 对于稀疏矩阵，在MATLAB中提供了函数luinc()来做不完全LU分解，其具体用法如下： [L U]= luinc(X,DROPTOL)，其中X、L和U的含义与函数lu()中的变量相同，DROPTOL为不完全LU分解的丢失容限。当DROPTOL设为0时，退化为完全LU分解。 [L,U] = luinc(X,0)，0级不完全LU分解。 [L,U,P] = luinc(X,0)，0级不完全LU分解。 * 3．QR分解 QR分解就是将m×n的矩阵A分解为m×n的矩阵Q和n×n的上三角矩阵R的乘积，且Q*Q=I，即A=Q*R。 在MATLAB中QR分解是由函数qr()来实现，其具体用法如下： [Q,R] = qr(A) 满足A=Q*R。 R = qr(A)， 返回上三角矩阵R。 4．奇异值分解 奇异值分解就是将m?n的矩阵A分解为A=U*S*V’，其中U为m?m的酉矩阵，V为n?n的酉矩阵，S为m?n的矩阵，并可如下表示： ，其中 ， ， 在MATLAB中奇异值分解是由函数svd()来实现，其具