factorization_luinc.m

来自「matlab基础程序 便于初学者学习」· M 代码 · 共 40 行

% factorization_luinc.m
% 矩阵的LU不完全分解示例
load west0479
A=west0479;
% 绘制稀疏矩阵图形
subplot(1,2,1);
spy(A);
title('A');
% LU分解结果图形
LU=lu(A);
subplot(1,2,2);
spy(LU);
title('LU');
% 0级LU分解示例
figure;
[L,U,P] = luinc(A,'0');
D=(L*U).*spones(P*A)-P*A;
subplot(2,2,1);
spy(L);
title('L:luinc(A,0)');
subplot(2,2,2);
spy(U);
title('U:luinc(A,0)');
subplot(2,2,3);
spy(P*A);
title('P*A:luinc(A,0)');
subplot(2,2,4);
spy(L*U);
title('L*U:luinc(A,0)');
% 丢失容限的比较
figure;
[L1,U1,P1]=luinc(A,1e-8);
[L2,U2,P2]=luinc(A,1e-2);
subplot(1,2,1);
spy(L1*U1);
title('luinc(A,1e-8)');
subplot(1,2,2);
spy(L2*U2);
title('luinc(A,1e-2)');