% factorization_svds.m
% 稀疏矩陣的奇異值分解操作示例
load west0479
s=svd(full(west0479));
disp('最大的四個奇異值：')
sl=svds(west0479,4)
disp('最小的六個奇異值：')
ss=svds(west0479,6,0)
% 圖形顯示結果
plot(s1,'ro')
hold on
plot(ss,'g*')
hold on
plot(s,'bd')
set(gca,'Xtick',[0:1:5]);