Balanced truncation model order reduction in limited time intervals for large systems

被引:0
作者
Patrick Kürschner
机构
[1] Max Planck Institute for Dynamics of Complex Technical Systems,Computational Methods in Systems and Control Theory
来源
Advances in Computational Mathematics | 2018年 / 44卷
关键词
Lyapunov equation; Rational Krylov subspaces; Model order reduction; Balanced truncation; Matrix exponential; 15A16; 15A18; 15A24; 65F60; 93A15; 93C;
D O I
暂无
中图分类号
学科分类号
摘要
In this article we investigate model order reduction of large-scale systems using time-limited balanced truncation, which restricts the well known balanced truncation framework to prescribed finite time intervals. The main emphasis is on the efficient numerical realization of this model reduction approach in case of large system dimensions. We discuss numerical methods to deal with the resulting matrix exponential functions and Lyapunov equations which are solved for low-rank approximations. Our main tool for this purpose are rational Krylov subspace methods. We also discuss the eigenvalue decay and numerical rank of the solutions of the Lyapunov equations. These results, and also numerical experiments, will show that depending on the final time horizon, the numerical rank of the Lyapunov solutions in time-limited balanced truncation can be smaller compared to standard balanced truncation. In numerical experiments we test the approaches for computing low-rank factors of the involved Lyapunov solutions and illustrate that time-limited balanced truncation can generate reduced order models having a higher accuracy in the considered time region.
引用
收藏
页码:1821 / 1844
页数:23
相关论文
共 73 条
[1]  
Al-Mohy AH(2011)Computing the action of the matrix exponential, with an application to exponential integrators SIAM J. Sci. Comput. 33 488-511
[2]  
Higham NJ(2002)On the decay rate of Hankel singular values and related issues Syst. Cont. Lett. 46 323-342
[3]  
Antoulas AC(2015)Fast singular value decay for Lyapunov solutions with nonnormal coefficients SIAM J. Matrix Anal. Appl. 36 656-668
[4]  
Sorensen DC(1972)Solution of the matrix equation AX + XB = C: Algorithm 432 Comm. ACM 15 820-826
[5]  
Zhou Y(2017)Model reduction for systems with inhomogeneous initial conditions Syst. Control Lett. 99 99-106
[6]  
Baker J(2009)Error estimates and evaluation of matrix functions via the faber transform SIAM J. Numer. Anal. 47 3849-3883
[7]  
Embree M(2004)Solving large-scale control problems IEEE Control Syst. Mag. 14 44-59
[8]  
Sabino J(2016)Frequency-limited balanced truncation with low-rank approximations SIAM J. Sci. Comput. 38 A471-A499
[9]  
Bartels RH(2013)Numerical solution of large and sparse continuous time algebraic matrix Riccati and Lyapunov equations: a state of the art survey GAMM Mitteilungen 36 32-52
[10]  
Stewart GW(2015)Generalized rational krylov decompositions with an application to rational approximation SIAM J. Matrix Anal. Appl. 36 894-916