symrcm
稀疏的反向Cuthill-McKee订购
语法
r = symrcm (S)
描述
r = symrcm (S)返回的对称反向Cuthill-McKee顺序年代。这是一个排列r这样年代(r, r)往往有其非零元素接近对角线。这是一个很好的预订陆或柯列斯基矩阵的分解来自长,瘦的问题。对称和非对称的订购工作年代。
对于一个真实的,对称稀疏矩阵,年代的特征值年代(r, r)是相同的吗年代,但eig (S (r, r))可能花费更少的时间来计算eig (S)。
例子
反向Cuthill-McKee订购
该声明
B =巴基;
使用一个功能演示工具箱生成邻接图截二十面体。这就是人们所熟知的是足球,穹顶建筑巴克明斯特·富勒(因此得名巴基),或者,最近,要想作为一个分子。有60个顶点。编号其中一半的顶点已被命令从一个半球,由国防部五角大楼;然后反射到另一个半球并将两部分组合在一起。
用这个编号,矩阵没有特别窄的带宽,作为第一个间谍阴谋显示:
图();次要情节(1、2、1),间谍(B)、标题(“B”)
获得反向Cuthill-McKee排序:
p = symrcm (B);R = B (p, p);
的间谍图显示了一个更窄的带宽。
次要情节(1、2、2),间谍(R)、标题(“B (p, p)”)
这个例子是仍在页面的引用symamd。
带宽也可以计算:
(i, j) =找到(B);bw = max (i j) + 1;
的带宽B和R分别是35岁和12岁。
算法
算法首先发现pseudoperipheral矩阵图的顶点。然后生成一个顶点水平结构通过广度优先搜索和订单减少pseudoperipheral顶点的距离。实现紧密地基于SPARSPAK实现所描述的乔治和刘。
引用
[1]乔治,艾伦和约瑟夫·刘计算机解决大型稀疏正定系统,普伦蒂斯·霍尔出版社,1981年。
[2]吉尔伯特,约翰·R。,克里夫硅藻土,罗伯特•施赖伯“稀疏矩阵在MATLAB:设计和实现”,暹罗在矩阵分析》杂志上,1992年。稍微扩展版本也可作为一个技术报告从施乐公司帕洛阿尔托研究中心。
之前介绍过的R2006a
这个主题有帮助吗?
