贝格姆·恩塞斯·卡纳塔罗;阿尼尔·拉奥（Anil V.Rao）。;蒂莫西·戴维斯（Timothy A.Davis）。 最优控制正交配置方法中产生的Karush-Kuhn-Tucker矩阵的状态缺陷约束对图粗化方法。 （英语） Zbl 1348.49024号 计算。最佳方案。申请。 64，第3期，793-819（2016）. 摘要：描述了一种状态缺陷约束对图粗化方法，用于提高通过Legendre-Gauss-Radau正交配置法离散最优控制问题所产生的大型稀疏Karush-Kuhn-Tucker矩阵的数值分解过程中的计算效率。该方法利用了正交配置法产生的Karush-Kuhn-Tucker矩阵的特殊稀疏结构。状态缺陷约束对图粗化方法将状态的每个分量与其相应的缺陷约束进行配对，并强制配对的行在重新排序的Karush-Kuhn-Tucker矩阵中相邻。使用各种基准最优控制问题给出了聚合的状态缺陷约束配对结果，其中发现，所提出的状态缺陷约束对图粗化方法显著减少了延迟枢轴数和浮点运算数，并通过单位时间执行更多的浮点运算提高了计算效率。然后证明，当最优控制问题包含状态和控制等式路径约束时，状态缺陷约束对图粗化方法对由Legendre-Gauss-Radau配置产生的Karush-Kuhn-Tucker矩阵的效果较差，因为这样的矩阵可能有对应于缺陷和路径约束的延迟枢轴。然后使用使用最大匹配的交替图粗化方法来尝试进一步减少延迟枢轴的数量。然而，我们发现，与状态缺陷约束对图粗化方法相比，这种交替图粗化法没有提供更多的优势。 引用于2文件 MSC公司： 49平方米25 最优控制中的离散逼近 49立方米 基于非线性规划的数值方法 65千5 数值数学规划方法 90立方 非线性规划 关键词：最优控制;离散化;正交配置法;卡鲁什-库恩-塔克系统;图形粗化 软件：MA47型;MA27号机组;SOCS系统;C解析;帕迪索;AMD公司;伊波特;MUMPS公司;鼻涕;GPOPS公司;MA57型;稀疏矩阵 PDF格式BibTeX公司 XML格式引用 \textit｛B.Ş.Cannataro｝等人，计算。最佳方案。申请。64，第3号，793--819（2016；Zbl 1348.49024） 全文： 内政部 参考文献： [1] Gill，P.E.，Murray，W.，Saunders，M.A.：SNOPT：大规模约束优化的SQP算法。SIAM J.Optim公司。12(4), 979-1006 (2002) ·Zbl 1027.90111号 ·doi:10.1137/S1052623499350013 [2] Biegler，L.T.、Zavala，V.M.：使用IPOPT的大规模非线性规划：企业范围动态优化的集成框架。计算。化学。工程33（3），575-582（2009）·doi:10.1016/j.compchemeng.2008.08.006 [3] Duff，I.S.，Reid，J.K.：MA47，不定稀疏对称线性系统直接解的Fortran代码。RAL报告，第95-001页（1995年）·Zbl 1070.65534号 [4] 达夫，I.S.:MA57——稀疏对称定和不定系统解的代码。ACM事务处理。数学。柔和。30(2), 118-144 (2004) ·Zbl 1070.65525号 ·doi:10.1145/992200.992202 [5] Amestoy，P.R.，Duff，I.S.，L'Excellent，J.-Y.，Koster，J.：MUMPS:通用分布式内存稀疏矩阵解算器。斯普林格，海德堡（2001） [6] Schenk，O.，Gärtner，K.：PARDISO中的两级动态调度：提高共享内存多处理系统的可伸缩性。并行计算。28(2), 187-197 (2002) ·Zbl 0982.68195号 ·doi:10.1016/S0167-8191（01）00135-1 [7] Elnagar，G.，Kazemi，M.，Razzaghi，M.:离散最优控制问题的伪谱勒让德方法。IEEE传输。自动。控制40（10），1793-1796（1995）·Zbl 0863.49016号 ·数字对象标识代码：10.1109/9.467672 [8] Elnagar，G.，Razzaghi，M.：线性二次型最优控制问题的配置型方法。最佳方案。控制应用程序。方法18（3），227-235（1998）·Zbl 0873.49023号 ·doi:10.1002/（SICI）1099-1514（199705/06）18:3<227:：AID-OCA598>3.0.CO；2-A型 [9] Fahroo，F.，Ross，I.M.：用勒让德伪谱法估算成本。J.制导控制动态。24（2），270-277（2001）·数字对象标识代码：10.2514/2.4709 [10] Fahroo，F.，Ross，I.M.：切比雪夫伪谱法直接轨迹优化。J.制导控制动态。25（1），160-166（2002）·数字对象标识代码：10.2514/2.4862 [11] Benson，D.A.、Huntington，G.T.、Thorvaldsen，T.P.、Rao，A.V.：通过正交配置法进行直接轨迹优化和成本状态估计。J制导控制动态。29(6), 1435-1440 (2006) ·数字对象标识代码：10.2514/1.20478 [12] Huntington，G.T.、Benson，D.A.、Rao，A.V.：四面体航天器编队的最佳配置。J.宇航员。科学。55(2), 141-169 (2007) ·doi:10.1007/BF03256518 [13] Huntington，G.T.，Rao，A.V.：使用高斯伪谱方法对航天器编队进行优化重构。J制导控制Dyn。31(3), 689-698 (2008) ·数字对象标识代码：10.2514/1.31083 [14] Gong，Q.，Fahroo，F.，Ross，I.M.：最优控制中伪谱方法的谱算法。J制导控制动态。2008年5月-6月31日（3）·Zbl 1165.49034号 [15] Gong，Q.，Ross，I.M.，Kang，W.，Fahroo，F.：向量映射定理与伪谱方法收敛性之间的联系。计算。最佳方案。申请。41(3), 307-335 (2008) ·Zbl 1165.49034号 ·文件编号：10.1007/s10589-007-9102-4 [16] Rao，A.V.，Benson，D.A.，Darby，C.L.，Francolin，C.，Patterson，M.A.，Sanders，I.，Huntington，G.T.：Algorithm 902:GPOPS，一种使用高斯伪谱方法解决多相最优控制问题的matlab软件。ACM事务处理。数学。柔和。37（2），第22条，39页，2010年4月-6月·Zbl 1364.65131号 [17] Darby，C.L.，Hager，W.W.，Rao，A.V.：用于解决最优控制问题的hp自适应伪谱方法。最佳方案。控制应用程序。方法32（4），476-502（2011）·Zbl 1266.49066号 ·doi:10.1002/oca.957 [18] Darby，C.L.，Hager，W.W.，Rao，A.V.：使用可变低阶自适应伪谱方法进行直接轨迹优化。J.太空船。火箭48（3），433-445（2011）·数字对象标识代码：10.2514/1.52136 [19] Kameswaran，S.，Biegler，L.T.：使用雷达点配置直接转录最优控制问题的收敛速度。计算。最佳方案。申请。41(1), 81-126 (2008) ·Zbl 1219.49029号 ·doi:10.1007/s10589-007-9098-9 [20] Garg，D.，Patterson，M.A.，Darby，C.L.，Francolin，C.，Huntington，G.T.，Hager，W.W.，Rao，A.V.：通过雷达伪谱方法对有限时域和无限时域最优控制问题进行直接轨迹优化和代价估计。计算。最佳方案。申请。49(2), 335-358 (2011). doi:10.1007/s10589-00-09291-0·兹比尔1226.49026 ·doi:10.1007/s10589-00-09291-0 [21] Garg，D.，Patterson，M.A.，Hager，W.W.，Rao，A.V.，Benson，D.A.，Huntington，G.T.：使用伪谱方法数值求解最优控制问题的统一框架。Automatica 46（11），1843-1851（2010）。doi:10.1016/j.automatica.2010.06.048·Zbl 1219.49028号 ·doi:10.1016/j.automatica.2010.06.048 [22] Garg，D.，Hager，W.W.，Rao，A.V.：解决无限小时最优控制问题的伪谱方法。Automatica 47（4），829-837（2011）。doi:10.1016/j.automatica.2011.01.085·Zbl 1215.49040号 ·doi:10.1016/j.automatica.2011.01.085 [23] Patterson，M.A.，Rao，A.V.:GPOPS-II：一个MATLAB软件，用于使用hp自适应高斯正交配置方法和稀疏非线性规划解决多相最优控制问题。ACM事务处理。数学。柔和。41(1), 1:1-1:37 (2014) ·Zbl 1369.65201号 ·doi:10.1145/2558904 [24] Patterson，M.A.、Hager，W.W.、Rao，A.V.：用于优化控制的ph网格细化方法。最佳方案。控制应用程序。方法36（4），398-421（2015）。doi:10.1002/oca2114·Zbl 1329.49048号 [25] Francolin，C.C.，Hager，W.W.，Rao，A.V.：使用积分高斯正交配置方法的最优控制中的Costate近似。最佳方案。控制应用程序。方法36（4），381-397（2015）。doi:10.1002/oca2112·Zbl 1329.49045号 [26] Karypis，G.，Kumar，V.：一种快速、高质量的划分不规则图的多级方案。SIAM J.科学。计算。20(1), 359-392 (1998) ·Zbl 0915.68129号 ·doi:10.1137/S1064827595287997 [27] Pothen，A.、Simon，H.D.、Wang，L.、Barnard，S.T.：朝向快速实现光谱嵌套分离。参见：1992年ACM/IEEE超级计算会议记录，第42-51页。IEEE计算机学会出版社（1992）·Zbl 0436.65021号 [28] Pothen，A.，Simon，H.D.，Liou，K.-P.：用图的特征向量划分稀疏矩阵。SIAM J.矩阵分析。申请。11(3), 430-452 (1990) ·兹比尔0711.65034 ·数字对象标识代码：10.1137/0611030 [29] Miller，G.L.，Teng，S.-H.，Thurston，W.，Vavasis，S.A.：自动网格划分。柏林施普林格（1993）·兹比尔0803.68083 ·doi:10.1007/978-14613-8369-73 [30] Miller，G.L.，Teng，S.-H.，Vavasis，S.A.：图形分隔符的统一几何方法。在：诉讼。，第32届计算机科学基础年会，第538-547页。IEEE（1991） [31] Cheng，C.-K.，Wei，Y.C.：性能稳定的改进双向分区算法[VLSI]。IEEE传输。计算。辅助设计。集成。电路系统。10(12), 1502-1511 (1991) ·数字对象标识代码：10.1109/43103500 [32] Hendrickson，B.，Leland，R.W.：划分图的多级算法。SC 95、28（1995）·兹伯利0816.68093 [33] Mansour，N.，Ponnusamy，R.，Choudhary，A.，Fox，G.：物理优化方法的图形收缩：并行计算机上映射数据的质量-成本权衡。摘自：《第七届超级计算国际会议记录》，第1-10页。ACM（1993）·Zbl 1165.49034号 [34] Cong，J.，Shinnerl，J.R.：VLSICAD中的多级优化。柏林施普林格出版社（2003）·Zbl 1030.00052号 ·doi:10.1007/978-1-4757-3748-6 [35] Safro，I.，Sanders，P.，Schulz，C.：图形分区的高级粗化方案，实验算法。柏林施普林格出版社（2012）·Zbl 1347.68355号 [36] Word，D.P.、Kang，J.、Akesson，J.和Laird，C.D.：大规模非线性动态优化问题的高效并行解决方案。计算。最佳方案。申请。59, 667-688 (2014) ·Zbl 1303.90122号 ·doi:10.1007/s10589-014-9651-2 [37] Patterson，M.A.，Rao，A.V.：利用直接配置伪谱方法中的稀疏性来解决最优控制问题。J.太空船。火箭队。49(2), 364-377 (2012) ·doi:10.2514/1.A32071 [38] Nocedal，J.，Overton，M.L.：非线性约束优化的投影Hessian更新算法。SIAM J.数字。分析。22(5), 821-850 (1985) ·Zbl 0593.65043号 ·doi:10.1137/0722050 [39] 西尔维斯特，J.J.：十九。证明了每个齐次二次多项式都可以通过实正交代换化简为正负平方和的形式。伦敦。爱丁堡。都柏林菲洛斯。《科学杂志》。4（23），138-142（1852） [40] Levenberg，K.：用最小二乘法求解某些问题的方法。问：申请。数学。2, 164-168 (1944) ·Zbl 0063.03501号 [41] Wächter，A.，Biegler，L.T.：关于大规模非线性规划的点内滤波器线性搜索算法的实现。数学。程序。106(1), 25-57 (2006) ·Zbl 1134.90542号 ·doi:10.1007/s10107-004-0559-y [42] Duff，I.S.，Reid，J.K.：MA27——一组用于求解稀疏对称线性方程组的Fortran子程序。UKAEA原子能研究所（1982年）·Zbl 0863.49016号 [43] Duff，I.S.，Reid，J.K.：不定稀疏对称线性的多边解。ACM事务处理。数学。柔和。（TOMS）9（3），302-325（1983）·Zbl 0515.65022号 ·数字对象标识代码：10.1145/356044.356047 [44] Liu，J.W.：稀疏矩阵解的多波前方法：理论与实践。SIAM版本34（1），82-109（1992）·Zbl 0919.65019号 ·数字对象标识代码：10.1137/1034004 [45] Betts，J.T.：使用非线性规划进行最优控制和估计的实用方法，第19卷。费城SIAM出版社（2010）·Zbl 1189.49001号 ·数字对象标识代码：10.1137/1.9780898718577 [46] 赵勇杰：滑翔机动态飞行的最佳模式。Optim。控制应用程序。方法25（2），67-89（2004）·Zbl 1073.49033号 ·doi:10.1002/oca.739 [47] Goddard，R.H.：一种达到极端高度的方法。《自然》105，809-811（1920）·数字对象标识代码：10.1038/105809a0 [48] Ledzewicz，U.，Schättler，H.：肿瘤抗血管生成数学模型的最佳控制分析。最佳方案。控制应用程序。方法29（1），41-57（2008）·doi:10.1002/oca.814 [49] Bryson，A.E.：应用最优控制：优化、估计和控制。CRC出版社，纽约（1975年） [50] Bulirsch，R.，Nerz，E.，Pesch，H.J.，von Stryk，O.：最优控制中直接和间接方法的结合：悬挂滑翔机的最大射程。柏林施普林格（1993）·Zbl 0808.65067号 [51] Sakawa，Y.：使用最优控制的自由飞行机器人的轨迹规划。最佳方案。控制应用程序。方法20（5），235-248（1999）·doi:10.1002/（SICI）1099-1514（199909/10）20:5<235:：AID-OCA658>3.0.CO；2-I型 [52] Tinney，W.F.，Walker，J.W.：通过最优有序三角因式分解直接求解稀疏网络方程。程序。IEEE 55（11），1801-1809（1967）·doi:10.1109/PROC.1967.6011 [53] Davis，T.A.，Amestoy，P.，Duff，I.S.：近似最小度排序算法。SIAM J.矩阵分析。申请。17(4), 886-905 (1996) ·Zbl 0861.65021号 ·doi:10.1137/S0895479894278952 [54] George，A.：规则有限元网格的嵌套剖分。SIAM J.数字。分析。10(2), 345-363 (1973) ·Zbl 0259.65087号 ·doi:10.1137/0710032 [55] Gibbs，N.E.，Poole Jr，W.G.，Stockmeyer，P.K.：减少稀疏矩阵带宽和轮廓的算法。SIAM J.Numer。分析。13(2), 236-250 (1976) ·Zbl 0329.65024号 ·数字对象标识代码：10.1137/0713023 [56] Chan，W.-M.，George，A.：反向Cuthill-Makee算法的线性时间实现。位数字。数学。20(1), 8-14 (1980) ·Zbl 0436.65021号 ·doi:10.1007/BF01933580 [57] Amestoy，P.R.，Davis，T.A.，Duff，I.S.：算法837:AMD，一种近似最小度排序算法。ACM事务处理。数学。柔和。30(3), 381-388 (2004) ·Zbl 1070.65534号 ·doi:10.1145/1024074.1024081 [58] Davis，T.A.，Hu，Y.：佛罗里达大学稀疏矩阵集合。ACM事务处理。数学。柔和。（TOMS）38（1），1（2011）·Zbl 1365.65123号 [59] Davis，T.A.：稀疏线性系统的直接方法，第2卷。费城SIAM出版社（2006）·Zbl 1119.65021号 ·doi:10.137/1.9780898718881 此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。