×

基于椭球方法的多级马尔可夫链算法的单范数最小化迭代重组。 (英语) Zbl 1308.65015号

摘要:最近,研究表明,通过添加基于迭代重组的外部迭代,可以加快计算稀疏不可约马尔可夫链平稳分布的一类多重网格方法的收敛速度。加速是通过使用二次规划方法选择具有概率约束的先前精细级迭代的线性组合来实现的,以最小化残差的两范数。在本文中,我们研究了最小化残差的一个范数的替代方法。这就产生了一个非线性凸规划,必须在每个加速步骤中进行求解。为了解决这个最小化问题,我们建议对非线性凸规划使用深截椭球方法。本文的主要目的是研究是否可以通过这种方式获得在执行时间和鲁棒性方面具有竞争力的迭代重组方法。我们推导了单范数目标函数和约束函数的次梯度公式,并说明了如何构造一个初始椭球体,该椭球体保证包含精确解,并给出了其存在的条件。我们还研究了使用椭球方法最小化两范数。数值试验表明,单范数和双范数加速程序在多重网格循环数方面也有类似的减少。测试还表明,在运行时间方面,单范数椭球加速度与双范数二次规划加速度具有竞争力,并且具有更好的鲁棒性。

MSC公司:

65立方厘米 马尔可夫链的数值分析或方法
60J22型 马尔可夫链中的计算方法
90C25型 凸面编程
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Bause F.,Kritzinger P.:随机Petri网。施普林格,德国(1996)·Zbl 0872.60068号
[2] Bazaraa M.S.、Sherali H.D.、Shetty C.M.:《非线性规划:理论与算法》,第三版。新泽西州威利(2006)·Zbl 1140.90040号
[3] 伯曼A.,普莱蒙斯R.J.:数学科学中的非负矩阵。宾夕法尼亚州费城SIAM(1987)·兹伯利048.41516
[4] Bertsimas D.,Tsitisklis J.N.:线性优化导论。Athena Scientific,马萨诸塞州贝尔蒙特(1997)
[5] Bland R.G.、Goldfarb D.、Todd M.J.:椭球体方法:调查。操作。第29(6)号决议,1039–1091(1981)·Zbl 0474.90056号 ·doi:10.1287/opre.29.6.1039
[6] Brandt A.,Mikulinsky V.:关于多重网格算法中的重组迭代和小岛问题。SIAM J.科学。计算。16, 20–28 (1995) ·Zbl 0826.65022号 ·doi:10.1137/0916002
[7] Briggs W.L.、Henson V.E.、McCormick S.F.:《多重网格教程》,第二版。宾夕法尼亚州费城SIAM(2000)·Zbl 0958.65128号
[8] Buchholz P.:结构化马尔可夫链的多级解。SIAM J.矩阵分析。申请。22(2), 342–357 (2000) ·Zbl 0967.60075号 ·doi:10.1137/S0895479898342419
[9] Cao W.L.,Stewart W.J.:几乎不耦合马尔可夫链的迭代聚合/分解技术。雅克32(3),702–719(1985)·Zbl 0628.65145号 ·doi:10.1145/3828.214137
[10] Chatelin F.,Miranker W.L.:通过逐次近似方法的聚合加速。线性代数应用。43, 17–47 (1982) ·Zbl 0485.65023号 ·doi:10.1016/0024-3795(82)90242-7
[11] Sterck H.、Manteuffel T.、McCormick S.F.、Miller K.、Pearson J.、Ruge J.、Sanders G.:马尔可夫链的平滑聚合多重网格。SIAM J.科学。计算。32, 40–61 (2010) ·Zbl 1209.65011号 ·doi:10.1137/080719157
[12] Sterck H.、Manteuffel T.、McCormick S.F.、Miller K.、Ruge J.、Sanders G.:马尔可夫链的代数多重网格。SIAM J.科学。计算。32, 544–562 (2010) ·兹比尔1210.65016 ·doi:10.1137/090753589
[13] Sterck H.、Manteuffel T.、McCormick S.F.、Nguyen Q.、Ruge J.:马尔可夫链的多级自适应聚合,及其在网络排名中的应用。SIAM J.科学。计算。30, 2235–2262 (2008) ·Zbl 1173.65028号 ·doi:10.1137/070685142
[14] De Sterck H.,Manteufel T.,Miller K.,Sanders G.:马尔可夫链稳态解的自适应代数多级方法的顶级加速。高级计算。数学。35, 375–403 (2010) ·Zbl 1230.65013号 ·数字对象标识代码:10.1007/s10444-010-9168-x
[15] Sterck H.,Miller K.,Sanders G.,Winlaw M.:马尔可夫链的递归加速多级聚合。SIAM J.科学。计算。32(3), 1652–1671 (2010) ·Zbl 1213.65018号 ·数字对象标识代码:10.1137/090770114
[16] Dziuban S.T.、Ecker J.G.、Kupferschmid M.:在椭球算法中使用深度切割进行非线性规划。数学。程序。螺柱25、93–107(1985)·兹伯利0582.90084 ·doi:10.1007/BFb0121078
[17] Ecker J.G.,Kupferschmid M.:非线性规划的椭球算法。数学。程序。27, 83–106 (1983) ·Zbl 0526.90074号 ·doi:10.1007/BF02591966
[18] Frenk J.B.G.,Gromicho J.,Zhang S.:凸规划的深切椭球算法:理论与应用。数学。程序。63, 83–108 (1994) ·Zbl 0791.90041号 ·doi:10.1007/BF01582060
[19] Goldfarb D.,Todd M.J.:线性规划椭球算法的修改和实现。数学。程序。23, 1–19 (1982) ·Zbl 0477.90038号 ·doi:10.1007/BF01583776
[20] Grassmann W.,Taksar M.,Heyman D.:马尔可夫链的再生分析和稳态分布。操作。第33(5)号决议,1107-1116(1985)·兹比尔0576.60083 ·doi:10.1287/opre.3.5.1107
[21] Haviv M.:计算马尔可夫链平稳分布的聚合/分解方法。SIAM J.数字。分析。24(4), 952–966 (1987) ·Zbl 0637.65147号 ·doi:10.1137/0724062
[22] Horn R.A.、Johnson C.R.:矩阵分析。剑桥大学出版社,纽约州纽约市(1985)·Zbl 0576.15001号
[23] Horton,G.,Leutenegger,S.T.:稳态马尔可夫链的多级求解算法。1994年ACM SIGMETRICS计算机系统测量和建模会议记录,第191-200页(1994)
[24] Iudin D.B.,Nemirovskii A.S.:凸极值问题的信息复杂性和有效解决方法。Matekon Transl.公司。俄罗斯、东欧、数学。经济。13, 3–25 (1976)
[25] Khachiyan L.G.:线性规划中的多项式算法。苏联。数学。多克拉迪20,191–194(1976年)
[26] Koury J.R.,McAllister D.F.,Stewart W.J.:计算几乎完全可分解马尔可夫链平稳分布的迭代方法。SIAM J.阿尔及利亚。光盘。方法。5(2), 164–186 (1984) ·Zbl 0605.60064号 ·doi:10.1137/0605019
[27] Krieger U.R.:用代数多重网格技术数值求解大型有限马尔可夫链。收录:Stewart,W.(编辑)《马尔可夫链的数值解》,第403-424页。多德雷赫特·克鲁沃(1995)·Zbl 0860.90061号
[28] Krieger U.R.:关于马尔可夫链的两级多重网格求解方法。线性代数应用。223–224, 415–438 (1995) ·Zbl 0831.65149号 ·doi:10.1016/0024-3795(95)00166-O
[29] Leutenger,S.T.,Horton,G.:关于求解几乎完全可分解马尔可夫链的多级算法的效用。技术代表94-44,ICASE(1994)·Zbl 0862.60052号
[30] Lüthi H.-J.:关于用椭球方法求解变分不等式。数学。操作。第10(3)号决议,515–522(1985)·Zbl 0586.49016号 ·doi:10.1287/门10.3.515
[31] Mandel J.,Sekerka B.:迭代聚合方法的局部收敛证明。线性代数应用。51, 163–172 (1983) ·Zbl 0494.65014号 ·doi:10.1016/0024-3795(83)90157-X
[32] Marek I.,Mayer P.:计算随机矩阵平稳概率向量的迭代聚合/分解方法的收敛性分析。数字。线性代数应用。5, 253–274 (1998) ·Zbl 0937.65047号 ·doi:10.1002/(SICI)1099-1506(199807/08)5:4<253::AID-NLA124>3.0.CO;2-B型
[33] Marek I.,Mayer P.:计算随机矩阵平稳概率向量的几类迭代聚集/分解方法的收敛理论。线性代数应用。363, 177–200 (2003) ·Zbl 1018.65048号 ·doi:10.1016/S0024-3795(02)00333-6
[34] Molloy M.K.:IEEE传输。计算。C-31、913–917(1982)·doi:10.1109/TC.1982.1676110
[35] Philippe B.、Saad Y.、Stewart W.J.:马尔可夫链建模中的数值方法。操作。第40(6)号决议,1156–1179(1992)·兹比尔0764.65095 ·doi:10.1287/opre.40.6.1156
[36] Rockafellar R.T.:凸分析。新泽西州普林斯顿大学出版社(1970)·Zbl 0193.18401号
[37] Shor N.Z.:凸规划问题中带空间扩展的截割方法。控制论13,94–96(1977)
[38] Simon H.A.,Ando A.:动态系统中变量的集合。《计量经济学》29、111–138(1961)·Zbl 0121.5103号 ·doi:10.2307/1909285
[39] Stewart W.J.:马尔可夫链数值解简介。普林斯顿大学出版社,新泽西州普林斯顿(1994)·Zbl 0821.65099号
[40] Takahashi,Y.:马尔可夫链平稳分布数值计算的集总方法。东京理工大学信息科学系技术代表B-18(1975年)
[41] Treister E.,Yavneh I.:马尔可夫链的动态自适应平滑聚合。SISC 33、2927–2949(2011)·Zbl 1232.65167号 ·数字对象标识代码:10.1137/100799034
[42] Treister E.,Yavneh I.:随机矩阵特征问题的平方和拉伸多重网格。数字。线性代数。17, 229–251 (2010) ·Zbl 1240.65124号
[43] Trottenberg U.,Oosterlee C.W.,Schüller A.:多重网格。爱思唯尔学术出版社,加利福尼亚州圣地亚哥(2001)·Zbl 0976.65106号
[44] Washio T.,Oosterlee C.W.:非线性多重网格方案的Krylov子空间加速。电子。事务处理。数字。分析。6, 271–290 (1997) ·Zbl 0903.65096号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。