ch3_28.m

同名书配套光盘及例程高质量MATLAB教学书籍

load west0479;                  %加载MATLAB 7.0自带的一个稀疏矩阵
S = west0479;
[L,U,P] = luinc(S,'0');             %做不完全LU分解    
figure(1);                       %新开一个图形显示窗口
subplot(2,3,1);                   %把该窗口分成2*3个字窗口，并在第1个字窗口上画图
spy(S);                         %显示稀疏矩阵S
title('S');                        %给图添上标题
subplot(2,3,2);                   %在第2个字窗口上画图
spy(L);                         %显示稀疏矩阵L
title('L');                        %给图添上标题
subplot(2,3,3);                   %在第2个字窗口上画图
spy(U);                         %显示稀疏矩阵U
title('U');                        %给图添上标题
subplot(2,3,4);                   %在第2个字窗口上画图
spy(P*S);                       %显示稀疏矩阵P*S
title('P*S');                      %给图添上标题
subplot(2,3,5);                   %在第2个字窗口上画图
spy(L*U);                       %显示稀疏矩阵L*U
title('L*U');                     %给图添上标题