exm051532_1.m

来自「这是本人收集的一些程序源代码」· M 代码 · 共 16 行

%exm051532_1.m
clear all,n=200;
rand('state',1),randn('state',2)
A=sprandsym(n,0.015,0.1,1);
shg
subplot(1,2,1),spy(A,'b',10),title('Spy plot of matrix A')
subplot(1,2,2),d=symmmd(A);
spy(A(d,d),'b',10),title('Matrix A with Minimun degree ordering');
B=full(A);
% 
format short e
tic, L1=chol(B);t1=toc;
tic, L2=chol(A);t2=toc/t1;
tic, L3=chol(A(d,d));t3=toc/t1;
disp('    全元素阵      稀疏矩阵  	     最小排序阵'),disp([1,t2,t3])