eigenvalue_vector_sparse.m

来自「matlab从入门到精通第3章书后源码对初学者很有帮助」· M 代码 · 共 26 行

% eigenvalue_vector_sparse.m
% 稀疏矩阵的特征值和特征向量
% A=delsq(numgrid('C',30))得到一对称正定阵，维度是632维
% 特征值分布在区间(0,8)之间，特征值4出现了18次
A=delsq(numgrid('C',30)); 
d=eig(full(A)); 
[dum,ind] = sort(abs(d)); 
dlm=eigs(A); 
dsm=eigs(A,6,'sm');
dsmt=sort(dsm);
% 画出6个最大特征值
subplot(2,1,1) 
plot(dlm,'k+') 
hold on 
plot(d(ind(end:-1:end-5)),'ks') 
hold off 
legend('eigs(A)','eig(full(A))',3) 
set(gca,'XLim',[0.5 6.5])
% 画出6个最小特征值
subplot(2,1,2) 
plot(dsmt,'k+') 
hold on 
plot(d(ind(1:6)),'ks') 
hold off 
legend('eigs(A,6,''sm'')','eig(full(A))',2) 
set(gca,'XLim',[0.5 6.5])