matlab中矩阵重排列,稀疏矩阵重新排序 - MATLAB & Simulink Example - MathWorks 中国

本文展示了如何使用MATLAB中的不同矩阵重排序方法,如colperm、symrcm、amd和dissect,来减少稀疏矩阵Cholesky分解的计算时间和存储需求。通过实例分析,解释了每种方法如何影响矩阵结构,特别是针对杠铃矩阵的Cholesky分解。
摘要由CSDN通过智能技术生成

可视化稀疏矩阵

spy 图可以显示矩阵中的非零元素。

下面的 spy 图显示了从杠铃矩阵的一部分得到的稀疏对称正定矩阵。此矩阵描述类似杠铃的图中的连接。

load barbellgraph.mat

S = A + speye(size(A));

pct = 100 / numel(A);

spy(S)

title('A Sparse Symmetric Matrix')

nz = nnz(S);

xlabel(sprintf('Nonzeros = %d (%.3f%%)',nz,nz*pct));

733ebe8b728d852e5e0995b1d25ba007.png

以下是杠铃图。

G = graph(S,'omitselfloops');

p = plot(G,'XData',xy(:,1),'YData',xy(:,2),'Marker','.');

axis equal

80cb22857035585682edec98a88089f6.png

计算 Cholesky 因子

计算 Cholesky 因子 L,其中 S = L*L'。请注意,L

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值