一种求解高阻尼PageRank问题的加权块Arnoldi算法.pdfVIP

第 3O卷第 4期 江苏师范大学学报(自然科学版) Vo1．3O，NO．4 2012年 12月 JournalofJiangsuNormalUniversity(NaturalScienceEdition) Dec．，2Ol2 一 种求解高阻尼PageRank问题的 加权块Arnoldi算法 袁 梅 (江苏师范大学 数学与统计学 院，江苏 徐州 221116) 摘要 ：提 出了一种加权块 Arnoldi方法求解 PageRank问题．为 了加快算法 的收敛速度 ，采用子空间迭代法作为加 速策略．数值实验结果表明，当阻尼因子 靠近 1时，提出的加速加权块 Arnoldi算法比现有的一些Krylov子空间 方法优越． 关键词 ：Google；PageRank；Arnoldi；块 Arnoldi；Krylov子空间 中图分类号：O241．6 文献标识码 ：A 文章编号：2095—4298(2012)04—0050—06 A weightedblockArnoldimethodforsolvingPageRankproblem withhigh dampingfactors YuanMei (SchoolofMathematics8LStatistics，JiangsuNormalUniversity，Xuzhou221116，Jiangsu，China) Abstract：Inthispaper，aweightedblockArnoldimethodforsolvingPageRankproblem withhighdampingfactors isproposed．To improvetheconvergenceperformanceofthisalgorithm ，thesubspaceiterationmethodisadopted． Numericalexperimentshowsthatthenew algorithm issuperiortOmanystate—of—。artKrylovsubspacemethodsfor PageRankproblem ，especiallywhenthedampingfactorisclosetO1． Keywords：Google；PageRank；Arnoldi；blockArnoldi；Krylovsubspace 0 引言 近年来 ，随着互联网技术的飞速发展，Google已成为 目前使用范围最广、最受欢迎的搜索引擎之一 ，其 中心算法是 PageRank，它主要用于求解 Google矩阵的主特征向量 (PageRank向量)．在互联 网中，web超 链接结构可以看作是一个有向图，马尔可夫模型可以由一个 × 阶的概率矩阵P来表示 ，它的元素 表示 在一段时间内从网页 i跳到网页J的可能性．若网页i没有出度，那么矩阵P的第 i行元素为 0，称这样的网 页为悬挂点．一种修正的方法是将 e／n加在这些行上 ，其中e 是元素全为 1的行向量．因此定义 — P+如 ／ ， 其 中如果 i是悬挂点，那么d一1；否则 d一0．但是 ，并不能保证 P是不可约的．为此 ，通过秩 1矩阵E 一ev 进行修正 ， A 一 + (1一 d)E ， 其中p是列随机矩阵，v是个性化 向量 ，a(0a1)称为阻尼 因子[1]，一般取 a一0．85．对于随机不可约的 Google矩阵A，PageRank向量是它的最大特征值 1所对应的右特征向量．实验中发现 ，当阻尼因子a一0．99 时得到的PageRank向量要 比a一0．85时的更为准确 ¨2]．因此 ，本文重点研究高阻尼因子． 求解 PageRank向量的最基本方法是幂法 ]，但是 ，当阻尼因子