Research Article  Open Access
MohammadSahadet Hossain, M. Monir Uddin, "Efficient Techniques for Solving the Periodic Projected Lyapunov Equations and Model Reduction of Periodic Systems", Mathematical Problems in Engineering, vol. 2017, Article ID 4362641, 11 pages, 2017. https://doi.org/10.1155/2017/4362641
Efficient Techniques for Solving the Periodic Projected Lyapunov Equations and Model Reduction of Periodic Systems
Abstract
We have presented the efficient techniques for the solutions of largescale sparse projected periodic discretetime Lyapunov equations in lifted form. These types of problems arise in model reduction and state feedback problems of periodic descriptor systems. Two most popular techniques to solve such Lyapunov equations iteratively are the lowrank alternating direction implicit (LRADI) method and the lowrank Smith method. The main contribution of this paper is to update the LRADI method by exploiting the ideas of the adaptive shift parameters computation and the efficient handling of complex shift parameters. These approaches efficiently reduce the computational cost with respect to time and memory. We also apply these iterative Lyapunov solvers in balanced truncation model reduction of periodic discretetime descriptor systems. We illustrate numerical results to show the performance and accuracy of the proposed methods.
1. Introduction
Periodic systems and control theory have received a lot of attention in many areas of science and engineering, especially in the areas where the periodic control is deserved, such as spacecraft attitude control, control of industrial processes and mechanical systems, helicopter design, and nonlinear circuit design [1–6]. In this paper, we consider the linear periodic discretetime descriptor systems with timevarying dimensions in the formwhere , , , are periodic with a period , , , and . If all the matrices are nonsingular, then (1) can be transformed into a periodic standard system.
Model reduction means to replace the given mathematical model (1) by a much smaller model which preserves certain crucial properties of the original model, such as stability and passivity, and it is then amenable for simulation and analysis. Analysis and reducedorder modeling of system (1) may require solving largescale sparse projected periodic discretetime Lyapunov equations. Model reduction procedures of such periodic systems using direct solution approaches have been considered in [6–10]. Recently, attention has been devoted to iterative solutions of largescale sparse Lyapunov equations using the alternating directions implicit (ADI) method [11, 12], the Smith method [12–14], and Krylov subspace methods [15, 16]. Generalization of all these methods has also been considered in [17, 18].
However, lowrank iterative methods and their application in model reduction of such periodic descriptor systems have been considered in [19, 20]. The drawbacks of the ADI method proposed that there are the suboptimal shifts computations and the slow convergence of the iterative solutions to the exact solutions. In this paper, we proposed the adaptive shift parameters [21] computation instead of suboptimal parameters and explored the efficient handing of the complex shift [22] parameters. These two proposed schemes reduce the computational cost and fasten the convergence of the iterative solutions to exact solutions. Numerical results are discussed to show the efficiency of the techniques.
The paper outline is as follows. In Section 2, we briefly review the analogous timeinvariant representation of the periodic systems, known as lifted representation of periodic systems. We discuss the cyclic lifted representation of periodic descriptor system (1). We reviewed the causalnoncausal decomposition, the periodic projected Lyapunov equations [8], and their corresponding lifted forms [7]. We also review some important concepts of the discretetime descriptor system in linear timeinvariant (LTI) form and the balanced truncation (BT) model reduction in LTI case. We also discuss the basics of the generalized ADI method in this section. The main idea behind this is that the concepts in the LTI structures will help more precisely to understand their corresponding periodic interpretations. In Section 3, we discuss the efficient solution of the causal lifted Lyapunov equations using the adaptive shifts based LRADI method. We also reviewed the Smith method from [20] for noncausal lifted Lyapunov equations in Section 4. We consider a balanced truncation model reduction method for periodic descriptor systems in Section 5. Section 6 contains numerical results that illustrate the efficiencies of the described iterative methods in model reduction of periodic descriptor system. Section 7 deficits the conclusion and remarks.
2. Preliminaries
In this section, we review the analogous timeinvariant representation of the periodic systems, known as lifted representation. We also introduce some notations and briefly review the fundamental concept of balanced truncation model reduction for LTI discretetime systems. We, afterward, generalize these concepts and results in the periodic setting.
2.1. Cyclic Lifted Representation of Periodic Systems
Analysis and modeling of periodic discretesystems are often described by an analogous timeinvariant representation of the periodic systems, known as lifted representation [1, 2, 6]. Here we consider the cyclic lifted representation [4] of periodic descriptor system (1), given bywhereThe descriptor vector, system input, and output of (2) are related to those of (1) via respectively. We can define the important concepts and dynamics of the periodic systems, such as regularity and asymptotic stability, by using their corresponding lifted forms.
2.2. The PPDLEs and Their Lifted Forms
Consider periodic discretetime descriptor system (1), where the set of periodic matrix pairs is periodic stable [20, 23]. Then the periodic Gramians of discretetime periodic descriptor systems satisfy a class of projected discretetime Lyapunov equations [7, 8].
Theorem 1 (see [23]). (1) The periodic projected discretetime Lyapunov equations with special right hand sideswith , have the unique symmetric, positive semidefinite periodic solutions and , respectively.
(2) The periodic projected discretetime Lyapunov equations with special right hand sideswith , have the unique symmetric, positive semidefinite periodic solutions and , respectively.
The solutions of (5) and of (6), respectively, are called the causal and noncausal reachability Gramians of system (1). Similarly, the solutions and of (7) and (8), respectively, are called the causal and noncausal observability Gramians of (1).
Theorem 2 (see [7]). Consider a periodic descriptor system (1) and the corresponding lifted form (2), where is periodic stable. Then the matrices and solve the lifted projected discretetime algebraic Lyapunov equations (LPDALEs):respectively, where , , and are as in (3), and
The matrices and are the solutions of the lifted causal and noncausal reachability Gramians, respectively. A similar result can also be stated for the lifted causal and noncausal observability Gramians [7]. The solution techniques for such Lyapunov equations are shown in [20] using lowrank ADI methods. In the following section we mainly focus on the solution techniques of (9a) more efficiently.
2.3. Balanced Truncation for LTI DiscreteTime Systems
In order to understand the BT model reduction method for the timevarying discretetime systems, we briefly review the idea for the timeinvariant case. Consider a LTI discretetime system:in which is nonsingular, and , , and . Assume that the original system is asymptotically stable; that is, all finite eigenvalues of the system pencil lie inside the unite circle; then the controllability Gramian and the observability Gramian can be computed by solving the discretetime Lyapunov equations (see, e.g., [24]):Equivalently, to compute the controllability and observability Gramians we can solve the continuoustime Lyapunov equationswhere is the Cayleytransformed pencil [25], and , and . Note that the solutions of (12) and (13) are identical to the solutions of (14) and (15), respectively [20, 25]. It is true that the Cayley transformation can destroy the system’s sparsity pattern where it needs to be restored. But, using proper permutations and efficient computation, one can recover this difficulty.
Since both the Gramians are symmetric and positive semidefinite, they have symmetric Gramians factors
The main drawback of BT based model order reduction (MOR) techniques is the computation of the solutions of the Lyapunov equations that are used in computing the Gramian factors. During the recent decades several approaches have been developed that allow exploiting the fact that often all coefficient matrices are highly sparse and the numbers of inputs and outputs are very small compared to the number of DoFs. Then usually the Gramians and can be expressed in lowrank factored forms; that is, and are thin rectangular matrices. Therefore instead of computing the full Gramian factors it is quite reasonable to compute the thin rectangular matrices and .
2.4. LowRank ADI for LTI DiscreteTime Systems
One popular iterative method for computing the lowrank and is the generalized lowrank Cholesky factoralternating directions implicit (GLRADI) iteration [11, 26, 27]. In the following we briefly introduce the GLRADI methods to solve continuoustime algebraic Lyapunov equations (14) and (15). For model reduction of system (1) we have to solve such a kind of Lyapunov equations.
The GLRADI iterations successively compute the columns of the two Gramian factors following the simple steps [28]:
In this iteration either is applied to generates , or is used for .
A set of optimal ADI shift parameters or simply shift parameters are necessary for fast convergence of the GLRADI iterations. We will discuss the shift parameter selection criterion later in this section. In these iterations we also see that if the maximum number of iterations is greater than the number of shifts , then the shift parameters are used in a cyclic way.
In the GLRADI iterations all of the input matrices , , and are real, due to the complex shift parameters in each iteration step, the updated store complex data, which increases the overall complexity and memory requirements of the method. Moreover, in the balancing based methods using these complex Gramian factors yields complex reduced systems by performing some complex arithmetic operations. This problem is resolved in [22]. In this regard, the important assumption is that the selected ADI shift parameters should be proper.
Definition 3. The ADI shift parameters are called proper if and are complex conjugates of each other when one of them is complex.
In [22] it is shown that at the st iteration of the GLRADI iterations, can be computed bywhere . This identity states that, for , can be computed explicitly from , which releases us from solving a shifted linear system with . This idea also results in the following theorem.
Theorem 4. Let us assume a set of proper ADI shift parameters. For a pair of complex conjugate shifts , the two subsequent block iterates and of GLRADI iterations satisfy [22]
This theorem reveals that for a pair of complex conjugate shifts at any iteration step in the GLRADI iteration can be updated byA version of the GLRADI algorithm is summarized in [29, Algorithm 5], which computes lowrank real Gramian factors. Additionally, the GLRADI iterations can be stopped whenever the norm of the ADIresidualbecomes very small. But computing in Frobeniusnorm or 2norm in each iteration step is an expensive task, since in each iteration the resulting residual (21) is an matrix. Recently, in [30] the authors show that, in the th iteration, the ADIresidual in (21) can be represented as withand the in the GLRADI iterations can be expressed as [30]The authors in [30] also show that the ADIresidual factor defined in (23) can be computed by the expressionIn the case of real setting when we consider , one must computewhere is defined in (18). The rank of is at most , that is, the number of columns of . Therefore, the computation of the Frobeniusnorm or 2norm of in each iteration is extremely cheap. Applying these strategies (computation of real Gramian factors and lowrank residual based stopping techniques), the updated GLRADI will be briefly discussed in Section 3.
3. Efficient Solutions of Causal LPDALEs Using LowRank ADI Methods
The basic ADI method for linear systems has been proposed in [31] and then used in [11, 12, 17, 32–34] for solving (projected) continuoustime and discretetime Lyapunov equations. Consider the LPDALE (9a), where the matrix is singular. For such equations both the ADI and Smith iterations fail to converge [20]. This problem is circumvented by considering a generalized Cayley transformation [25]:which transforms the LPDALE (9a) to an equivalent projected continuoustime algebraic Lyapunov equation (PCALE):Here is the Cayleytransformed pencil. In this transformation the finite stable eigenvalues of are mapped to the finite stable eigenvalues of , and the infinite eigenvalue of is mapped to . More details of that transformation can be found in [20, 25]. The solutions of PCALE (28) and (9a) are unique and they are solvable for every [20].
3.1. Efficient Computation of LowRank Gramian Factors
We present the updated version of lowrank ADI methods for Lyapunov equation (28), which is also the solution of (9a). An updated algorithm for the lowrank ADI method of generalized state space systems can be found in [29]. Here the idea is generalized for the periodic descriptor system.
Applying ADI method the solution of Lyapunov equation (28) can be obtained by performing the following iteration steps, beginning with and :where , with , is a shift parameter and . Steps (29a) and (29b) can be written in a single stepwhere and . The error bound in iteration is given by (see, e.g., [35])Since , this leads toIn th iteration the residual can be computed bywhere .
We observe that the computed in (30) can be decomposed aswhere and . Since the number of inputs is very small compared to the dimension of the system, the right hand side of (9a) is a lowrank matrix. Therefore often this decomposition is possible, and can be computed by following iteration steps:with , and is maximum number of iterations. In each iteration in (35b) the number of columns of is increased by the number of columns . Therefore after few steps the number of columns of grows dramatically. In order to resolve this drawback a further optimization occurs, where the number of columns to modify is kept constant by using the clever reordering of the ADIparameters. This updated version of GLRADI [36] is given below:
Unfortunately, the GLRADI iteration does not preserve the block diagonal structure at every iteration step in Algorithm 1. But has almost block diagonal form since it approximates the solution of the LPDALE (9a) given by , where are the periodic solutions of (5). We then consider the block partitioning with , and the causal reachability Gramian can be approximated by the matrices for .

We compute lowrank Gramian factor , so that by following the steps (36a)–(36c). We summarize the resulting iteration using column compression in Algorithm 1. It should be noted that in Algorithm 1 the computation of the matrices and can be avoided by rewriting the GLRADI iteration in terms of the original matrices and [20]. We should also note that we do not compute the inverses of the matrices at step (2) of Algorithm 1 explicitly. We represent the matrices and with the original sparse matrices with and and use the sparse inversion technique in step (2) [20].
The iteration can be stopped in the th step if the residual norm, that is, , where is defined in (33), is sufficiently small. If we can compute as a lowrank factor then computing is very cheap. Here, we discuss how to compute the lowrank residualfactor from each iteration in the GLRCFADI algorithm. Following the procedure as shown in Section 2.4, iterate in (36b) can also be written asAlso,
3.2. Selection of the ADI Shift Parameters
The convergence rate of the lowrank ADI algorithms discussed above depends on a set of ADI shift parameters. The ADI shift parameters were originally introduced by Wachspress in [37] to solve Lyapunov equations using the ADI methods. The author shows that a set of optimal ADI shift parameters for system (2) can be computed by solving the socalled ADI minmax problem [37, 38]:where is the spectrum of the matrix .
For a largescale system, determining the entire spectrum of is almost impossible. Therefore, in the literature several techniques are proposed; see, for example, [26, 36, 39, 40] to solve the minmax problem (40) using a much smaller part of the spectrum. Currently, one such commonly used technique is Penzl’s heuristic approach introduced in [26], where Ritz values (see, e.g., [41]) and () reciprocal Ritz values with respect to and , respectively, are employed. A complete procedure of the heuristic approach can be found in [20]. However, for largescale descriptor systems, computing the Ritz values with respect to is a challenging task since is then not invertible.
Another ADI shift selection criterion that we emphasize here is the adaptive approach which is introduced in [21]. The performances of this adaptive shift selection approach for a largescale system are also investigated in [29, 42, 43]. There the authors show that the approach is superior to the heuristic approach, especially for the (generalized) descriptor systems. To initialize the shifts, we first compute the eigenvalues of the projected pencil where is formed by the basis of range of any () random matrix. Then select some desired number of optimal shifts by solving the ADI minmax problem like in the heuristic procedure. Then whenever all the initial shifts have been used, the pencil is projected to the span of the current (see, e.g., Algorithm 1 of [21]) and computes some optimal shifts from the eigenvalues of the projected pencil and uses them as the set of next shift parameters. This approach is repeated while the algorithm has not converged to the given tolerance.
4. Smith Method for Noncausal LPDALEs
Consider the LPDALE (9b). For nonsingular , this equation is equivalent to the LPDALE:
Note that if is the index of the system [20], then the approximate solution of the noncausal LPDALE (9b) can be obtained by the Smith iterations [14]:Therefore, the Cholesky factor of the solution of (42) and also of the LPDALE (9b) takes the formThus, the solution of (9b) can be obtained with few computations for systems of low index. The solution of the noncausal LPDALE (9b) is briefly discussed in [20].
Similar to the causal case, we consider the block partitioningsuch that the noncausal reachability Gramians of system (1) can be computed in factored form for .
Remark 5. The causal and noncausal observability Gramians and of system (1) and the corresponding the lowrank Cholesky factors and , where and , can also be determined from the corresponding LPDALEs that are dual to the LPDALEs (9a) and (9b); see [7] for details.
5. Model Reduction of Periodic Systems
For periodic descriptor system (1), we construct a reducedorder model of dimension of the formwhere , , , are periodic matrices, , and . Apart from having a much smaller state space dimension, the reducedorder model also preserves important physical properties of the original system such as regularity and stability and the small approximation error. Balanced truncation for discretetime periodic descriptor system has been considered in [7, 20]. Some important concepts related to balanced truncation of discretetime periodic descriptor systems has been presented in [8].
Consider is periodic stable, and the Cholesky factors of the causal and noncausal Gramians satisfyThen the causal and noncausal Hankel singular values and of periodic descriptor system (1) are defined asrespectively, where denotes the singular values of the corresponding product matrices. For a balanced system, one truncates the states related to the small causal Hankel singular values because they do not change system properties essentially. Unfortunately, truncation of small nonzero noncausal Hankel singular values may lead to unstable reduction [7].
Hence, we compute the singular value decompositions of the product matriceswhere , , , and are orthogonal:with , and is nonsingular for .
We then compute reducedorder system (46) aswhere the projection matrices and have the formwith and .
Let be the transfer function of lifted system (2) and be the transfer function of the corresponding reducedorder lifted system formed from the reducedorder matrices in (51). Then the norm error bound for the reducedorder system is given by where , contains the truncated causal Hankel singular values from (49).
6. Numerical Experiments
In this section, we present some numerical results. We consider an artificial periodic discretetime descriptor system which is reformed from its original model in [8, Example 1]. The reformulation of the model is described in [20]. Here all the computations are carried out using MATLAB® R2012b on an Intel CORE i5 processor with a 2.90GHz clock and 16 GB of RAM.
To show the comparisons of GLRADI (i.e., Algorithm 2 in [20]) and our UGLRADI (i.e., Algorithm 1), we compute the lowrank Gramian factors by solving the LPDALEs applying both the algorithms. We also investigate the performances of both the heuristic and adaptive shifts to implement these algorithms. For implementing the GLRADI iteration we chose 40 optimal heuristic shifts out of 50 large and 40 small magnitude Ritz values. To select 40 heuristic shifts it takes 3.98 seconds of CPU time. In case of the adaptive shifts (for implementing UGLRADI method), in each cycle 40 proper shift parameters are selected following the procedure discussed above. We require only 0.03 seconds of CPU time to compute this number of adaptive shift parameters. Note that for computing the initial shifts first we project the pencil onto the column space of an random matrix.
Figures 1 and 2, respectively, show the rate of convergence and the computational time of the GLRADI and UGLRADI methods. In both cases, we see that UGLRADI method performs better than the GLRADI method. We compute the reducedorder model of the periodic descriptor system by using the lowrank Gramian factors obtained from UGLRADI method. The reformed periodic descriptor system has , , and for . The set of periodic matrix pairs is periodic stable with and for . The original lifted system has order . The computed reducedorder model has subsystems of orders . The comparisons of the frequency responses and the error bounds of the original model and the reduced model are shown in Figure 3. Figure 3(a) shows the frequency responses of the original and reduced models match nicely. The estimated error is below the error bound which is reflected in Figure 3(b).
(a) Controllability Gramian factor
(b) Observability Gramian factor
(a) Controllability Gramian factor
(b) Observability Gramian factor
(a) Frequency responses
(b) Absolute error and error bound
7. Conclusions
In this article we have suggested the adaptive shifts based ADI method and reviewed the Smith method for computing the lowrank factors of the periodic Gramians for periodic discretetime descriptor systems. Later, these factors are used in a balancing based model reduction technique to find a reducedorder model for the periodic discretetime descriptor system. The proposed techniques save memory and the computational time to generate a reducedorder model. The capability and efficiency of the methods have been reflected in the numerical results.
Competing Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 S. Bittanti and P. Colaneri, “Invariant representations of discretetime periodic systems,” Automatica. A Journal of IFAC, the International Federation of Automatic Control, vol. 36, no. 12, pp. 1777–1793, 2000. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 M. Farhood, C. L. Beck, and G. E. Dullerud, “Model reduction of periodic systems: a lifting approach,” Automatica, vol. 41, no. 6, pp. 1085–1090, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 P. Misra, “Timeinvariant representation of discrete periodic systems,” Automatica, vol. 32, no. 2, pp. 267–272, 1996. View at: Publisher Site  Google Scholar  MathSciNet
 B. Park and E. Verriest, “Canonical forms of discrete linear periodically timevarying systems and a control application,” in Proceedings of the 28th Conference on Decision and Control, pp. 1220–1225, Tampa, Fla, USA, December 1989. View at: Google Scholar
 A. Varga, “Balancing related methods for minimal realization of periodic systems,” Systems & Control Letters, vol. 36, no. 5, pp. 339–349, 1999. View at: Publisher Site  Google Scholar  MathSciNet
 A. Varga, “A periodic systems toolbox for MATLAB,” in Proceedings of the 16th IFAC World Congress, Prague, Czech Republic, July 2005. View at: Google Scholar
 P. Benner, M.S. Hossain, and T. Stykel, “Model reduction of periodic descriptor systems using balanced truncation,” in Model Reduction for Circuit Simulation, P. Benner, M. Hinze, and J. ter Maten, Eds., vol. 74 of Lecture Notes in Electrical Engineering, pp. 193–206, Springer, Berlin, Germany, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 E. K. Chu, H.Y. Fan, and W.W. Lin, “Projected generalized discretetime periodic Lyapunov equations and balanced realization of periodic descriptor systems,” SIAM Journal on Matrix Analysis and Applications, vol. 29, no. 3, pp. 982–1006, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. Varga, “An overview of recent developments in computational methods for periodic systems,” in Proceedings of the 3rd IFAC Workshop “Periodic Control Systems” (PSYCO '07), pp. 157–162, St. Petersburg, Russia, August 2007. View at: Google Scholar
 A. Varga and P. Van Dooren, “Computing the zeros of periodic descriptor systems,” Systems & Control Letters, vol. 50, no. 5, pp. 371–381, 2003. View at: Publisher Site  Google Scholar  MathSciNet
 J.R. Li and J. White, “Low rank solution of Lyapunov equations,” SIAM Journal on Matrix Analysis and Applications, vol. 24, no. 1, pp. 260–280, 2002. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 T. Penzl, “A cyclic low rank Smith method for large sparse Lyapunov equations,” SIAM Journal on Scientific Computing, vol. 21, pp. 1401–1418, 2000. View at: Google Scholar
 S. Gugercin, D. C. Sorensen, and A. C. Antoulas, “A modified lowrank Smith method for largescale Lyapunov equations,” Numerical Algorithms, vol. 32, no. 1, pp. 27–55, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 R. Smith, “Matrix equation XA + BX = C,” SIAM Journal on Applied Mathematics, vol. 16, pp. 198–201, 1968. View at: Google Scholar
 K. Jbilou and A. J. Riquet, “Projection methods for large Lyapunov matrix equations,” Linear Algebra and its Applications, vol. 415, no. 23, pp. 344–358, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 V. Simoncini, “A new iterative method for solving largescale Lyapunov matrix equations,” SIAM Journal on Scientific Computing, vol. 29, no. 3, pp. 1268–1288, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 T. Stykel, “Lowrank iterative methods for projected generalized Lyapunov equations,” Electronic Transactions on Numerical Analysis, vol. 30, pp. 187–202, 2008. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 T. Stykel and V. Simoncini, “Krylov subspace methods for projected Lyapunov equations,” Applied Numerical Mathematics, vol. 62, no. 1, pp. 35–50, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 M.S. Hossain, Numerical methods for model reduction of timevarying descriptor systems [Ph.D. thesis], Chemnitz University Of Technology, Chemnitz, Germany, 2011, http://nbnresolving.de/urn:nbn:de:bsz:ch1qucosa74776.
 P. Benner, M.S. Hossain, and T. Stykel, “Lowrank iterative methods for periodic projected Lyapunov equations and their application in model reduction of periodic descriptor systems,” Numerical Algorithms, vol. 67, no. 3, pp. 669–690, 2014. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 P. Benner, J. Saak, and M. M. Uddin, “Balancing based model reduction for structured index2 unstable descriptor systems with application to flow control,” Numerical Algebra, Control and Optimization, vol. 6, no. 1, pp. 1–20, 2016. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 P. Benner, P. Kürschner, and J. Saak, “Efficient handling of complex shift parameters in the lowrank Cholesky factor ADI method,” Numerical Algorithms, vol. 62, no. 2, pp. 225–251, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 E.W. Chu, H.Y. Fan, and W.W. Lin, Reachability and Observability of Periodic Descriptor Systems, National Tsing Hua University, Hsinchu, Taiwan, 2005.
 M. Hinze and S. Volkwein, “Proper orthogonal decomposition surrogate models for nonlinear dynamical systems: error estimates and suboptimal control,” in Dimension Reduction of LargeScale Systems, P. Benner, V. Mehrmann, and D. Sorensen, Eds., vol. 45 of Lecture Notes in Computational Science and Engineering, pp. 261–306, Springer, Berlin, Germany, 2005. View at: Publisher Site  Google Scholar
 V. Mehrmann, “A step toward a unified treatment of continuous and discretetime control problems,” Linear Algebra and Its Applications, vol. 241–243, pp. 749–779, 1996. View at: Google Scholar
 T. Penzl, “A cyclic lowrank Smith method for large sparse Lyapunov equations,” SIAM Journal on Scientific Computing, vol. 21, no. 4, pp. 1401–1418, 2000. View at: Publisher Site  Google Scholar  MathSciNet
 D. Sorensen and Y. Zhou, “Bounds on eigenvalue decay rates and sensitivity of solutions to Lyapunov equations,” Tech. Rep. TR0207, Department of Computational and Applied Mathematics, Rice University, Houston, Tex, USA, 2002, http://www.caam.rice.edu/caam/trs/tr02.html#TR0207. View at: Google Scholar
 P. Benner, J.R. Li, and T. Penzl, “Numerical solution of largescale Lyapunov equations, Riccati equations, and linearquadratic optimal control problems,” Numerical Linear Algebra with Applications, vol. 15, no. 9, pp. 755–777, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. M. Uddin, Computational methods for model reduction of largescale sparse structured descriptor systems [Ph.D. thesis], Otto von Guericke University Magdeburg, Magdeburg, Germany, 2015.
 P. Benner, P. Kürschner, and J. Saak, “An improved numerical method for balanced truncation for symmetric secondorder systems,” Mathematical and Computer Modelling of Dynamical Systems, vol. 19, no. 6, pp. 593–615, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 D. W. Peaceman and J. Rachford, “The numerical solution of parabolic and elliptic differential equations,” Journal of the Society for Industrial and Applied Mathematics, vol. 3, pp. 28–41, 1955. View at: Google Scholar  MathSciNet
 P. Benner and H. Faßbender, “On the numerical solution of largescale sparse discretetime Riccati equations,” Advances in Computational Mathematics, vol. 35, no. 24, pp. 119–147, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 D. Calvetti, N. Levenberg, and L. Reichel, “Iterative methods for $XAXB=C$,” Journal of Computational and Applied Mathematics, vol. 86, no. 1, pp. 73–101, 1997. View at: Publisher Site  Google Scholar
 A. Lu and E. L. Wachspress, “Solution of Lyapunov equations by alternating direction implicit iteration,” Computers & Mathematics with Applications, vol. 21, no. 9, pp. 43–58, 1991. View at: Publisher Site  Google Scholar  MathSciNet
 F. D. Freitas, J. Rommes, and N. Martins, “Gramianbased reduction method applied to large sparse power system descriptor models,” IEEE Transactions on Power Systems, vol. 23, no. 3, pp. 1258–1270, 2008. View at: Publisher Site  Google Scholar
 J. Saak, Efficient numerical solution of large scale algebraic matrix equations in PDE control and model order reduction [Ph.D. dissertation], TU Chemnitz, 2009.
 E. L. Wachspress, “Iterative solution of the Lyapunov matrix equation,” Applied Mathematics Letters, vol. 1, no. 1, pp. 87–90, 1988. View at: Publisher Site  Google Scholar  MathSciNet
 E. Wachspress, The ADI Model Problem, Springer, New York, NY, USA, 2013. View at: Publisher Site  MathSciNet
 J.R. Li, Model reduction of large linear systems via low rank system gramians [Ph.D. thesis], Massachusettes Institute of Technology, 2000.
 P. Benner, H. Mena, and J. Saak, “On the parameter selection problem in the NewtonADI iteration for largescale Riccati equations,” Electronic Transactions on Numerical Analysis, vol. 29, pp. 136–149, 2008. View at: Google Scholar  MathSciNet
 G. H. Golub and C. F. Van Loan, Matrix Computations, Johns Hopkins University Press, Baltimore, Md, USA, 3rd edition, 1996. View at: MathSciNet
 P. Benner, J. Saak, and M. M. Uddin, “Structure preserving model order reduction of large sparse secondorder index1 systems and application to a mechatronics model,” Mathematical and Computer Modelling of Dynamical Systems. Methods, Tools and Applications in Engineering and Related Sciences, vol. 22, no. 6, pp. 509–523, 2016. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. M. Uddin, M. S. Hossain, and M. F. Uddin, “Rational Krylov Subspace Method (RKSM) for solving the lyapunov equations of index1 descriptor systems and application to balancing based model reduction,” in Proceedings of the 9th International Conference on Electrical Computer Engineering (ICECE '16), pp. 451–454, Dhaka, Bangladesh, December 2016. View at: Google Scholar
Copyright
Copyright © 2017 MohammadSahadet Hossain and M. Monir Uddin. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.