SVD-Krylov based techniques for structure-preserving reduced order modelling of second-order systems

Md. Motlubar Rahman; Mahtab Uddin; M. Monir Uddin; L. S. Andallah; Md. Motlubar Rahman; Mahtab Uddin; M. Monir Uddin; L. S. Andallah

doi:10.3934/mmc.2021006

Mathematical Modelling and Control

2021, Volume 1, Issue 2: 79-89. doi: 10.3934/mmc.2021006

Previous Article Next Article

Research article

SVD-Krylov based techniques for structure-preserving reduced order modelling of second-order systems

1.
Department of Mathematics, Jahangirnagar University, Savar, Dhaka-1342, Bangladesh
2.
Institute of Natural Sciences, United International University, Dhaka-1212, Bangladesh
3.
Department of Mathematics and Physics, North South University, Dhaka-1229, Bangladesh
Correction on: Mathematical Modelling and Control 1: 134-135

Received: 21 March 2021 Accepted: 16 June 2021 Published: 22 June 2021

We introduce an efficient structure-preserving model-order reduction technique for the large-scale second-order linear dynamical systems by imposing two-sided projection matrices. The projectors are formed based on the features of the singular value decomposition (SVD) and Krylov-based model-order reduction methods. The left projector is constructed by utilizing the concept of the observability Gramian of the systems and the right one is made by following the notion of the interpolation-based technique iterative rational Krylov algorithm (IRKA). It is well-known that the proficient model-order reduction technique IRKA cannot ensure system stability, and the Gramian based methods are computationally expensive. Another issue is preserving the second-order structure in the reduced-order model. The structure-preserving model-order reduction provides a more exact approximation to the original model with maintaining some significant physical properties. In terms of these perspectives, the proposed method can perform better by preserving the second-order structure and stability of the system with minimized $\mathcal{H}_2$ -norm. Several model examples are presented that illustrated the capability and accuracy of the introducing technique.

Keywords:

Citation: Md. Motlubar Rahman, Mahtab Uddin, M. Monir Uddin, L. S. Andallah. SVD-Krylov based techniques for structure-preserving reduced order modelling of second-order systems[J]. Mathematical Modelling and Control, 2021, 1(2): 79-89. doi: 10.3934/mmc.2021006

Related Papers:

[1]	Shanshan Yu, Jiang Liu, Xiaojie Lin . Multiple positive periodic solutions of a Gause-type predator-prey model with Allee effect and functional responses. AIMS Mathematics, 2020, 5(6): 6135-6148. doi: 10.3934/math.2020394
[2]	Heping Jiang . Complex dynamics induced by harvesting rate and delay in a diffusive Leslie-Gower predator-prey model. AIMS Mathematics, 2023, 8(9): 20718-20730. doi: 10.3934/math.20231056
[3]	Chuangliang Qin, Jinji Du, Yuanxian Hui . Dynamical behavior of a stochastic predator-prey model with Holling-type III functional response and infectious predator. AIMS Mathematics, 2022, 7(5): 7403-7418. doi: 10.3934/math.2022413
[4]	Binfeng Xie, Na Zhang . Influence of fear effect on a Holling type III prey-predator system with the prey refuge. AIMS Mathematics, 2022, 7(2): 1811-1830. doi: 10.3934/math.2022104
[5]	Manal Alqhtani, J.F. Gómez-Aguilar, Khaled M. Saad, Zulqurnain Sabir, Eduardo Pérez-Careta . A scale conjugate neural network learning process for the nonlinear malaria disease model. AIMS Mathematics, 2023, 8(9): 21106-21122. doi: 10.3934/math.20231075
[6]	Chuanfu Chai, Yuanfu Shao, Yaping Wang . Analysis of a Holling-type IV stochastic prey-predator system with anti-predatory behavior and Lévy noise. AIMS Mathematics, 2023, 8(9): 21033-21054. doi: 10.3934/math.20231071
[7]	Huahai Qiu, Li Wan, Zhigang Zhou, Qunjiao Zhang, Qinghua Zhou . Global exponential periodicity of nonlinear neural networks with multiple time-varying delays. AIMS Mathematics, 2023, 8(5): 12472-12485. doi: 10.3934/math.2023626
[8]	Chenxuan Nie, Dan Jin, Ruizhi Yang . Hopf bifurcation analysis in a delayed diffusive predator-prey system with nonlocal competition and generalist predator. AIMS Mathematics, 2022, 7(7): 13344-13360. doi: 10.3934/math.2022737
[9]	Xuejie Ma, Songhua Wang . A hybrid approach to conjugate gradient algorithms for nonlinear systems of equations with applications in signal restoration. AIMS Mathematics, 2024, 9(12): 36167-36190. doi: 10.3934/math.20241717
[10]	Qinghua Zhou, Li Wan, Hongbo Fu, Qunjiao Zhang . Exponential stability of stochastic Hopfield neural network with mixed multiple delays. AIMS Mathematics, 2021, 6(4): 4142-4155. doi: 10.3934/math.2021245

Abstract

1. Introduction

The predator-prey interactions perform the major evolution that is used to consider the mechanisms effects in terms of the population's collaboration based on the ecological civilizations. Many scientists have constructed the mathematical systems using the inspirations of the population of the predator into the prey together with the ecology form of the density conditions. Numerous predator-prey models represent the predator's growth that is directly related to the density of the prey together and vice versa. The growth of the prey is directly related to the predator density with the current time. It is predictable that the growth of the predator depends upon the density of the prey in the present and the recent past ^[1,2]. Therefore, the biological influences are more realistic to perform the fading impacts of the memory. Some relevant recent studies based on the past prey density using the predator-prey interactions are discussed in ^[3,4,5,6].

The Allee effects defined by Allee provide the positive impacts based on the correlation between the precise fitness and population size ^{[7,8,9,10,11,12]}. Most of the relations based on the predator-prey are assumed to the growth of the logistic prey. The Allee impacts using the dynamics of the prey signify the population problems, genetic drift, including insufficient alimentation at low densities, inbreeding depression, mate restriction, predator evading of resistance ^{[13,14,15,16,17]}. Therefore, the Allee effects are provided as a serious factor based on the biological controller and commonly improves the loss prospect of the individual ^[18]. These impacts can be applied to indicate the weak and strong ratio based on the growth. The strong Allee's impacts ^[19,20,21] provide the growth (positive) ratio per capita, whereas the weak Allee effects indicate the growth rate (negative) per capita. The occurrence of Allee effects has already been provided in numerous biological individuals, like as marine invertebrates (gastropod) ^[22], mammals (suricates) ^[23], and insects (Glanville Fritillary butterfly) ^[24].

Few of the transition states cannot be prompt and the predator-prey interactions based on the time delay systems have achieved huge significance over the last few years ^[25,26,27]. The natural population dynamics indicates the reliability that is associated with the species response. On the constancy of the population, there are various factors of the time delay ^[28,29]. Subsequently, the time delay factor is applied in diverse ordinary form of the differential models that can be provided more accurate form. Few of the models based on the time delay factors indicate the destabilization of the model using the accompaniment of the Hopf bifurcation with the strong oscillation dynamics ^[30,31]. There are a few investigations in the literature that have been presented based on the Allee effects, which makes the system stable with the delay factors ^[32]. The biological models, cooperation and competition are associated with the population dynamics along with the factor of time delay. The competition of the population usually rises based on the food, whereas the collaboration in the population normally indicates to assume the prey or to avoid from the predator.

The aim of such investigations is to present the numerical performances through the artificial neural networks (ANNs) along with the scaled conjugate gradient (SCG), ANNs-SCG for solving the predator-prey delay differential system of Holling form-III. The predator-prey dynamics along with the delay factors in the prey competition/cooperation using the Allee effects with fading memory has never been presented before through the ANNs-SCG technique. The stochastic based computational investigations have been presented in this study to solve various complex and complicated types of dynamical systems. To mention some of them are HIV dynamical models ^[33], a food chain system of models ^[34,35], coronavirus dynamical systems ^[36], differential form of the eye surgery model ^[37,38], singular thermal explosion theory model ^[39], singular differential models ^[40], and nonlinear smoking differential system ^[41]. Few novel features of the proposed stochastic study are indicated as:

● The design of the predator-prey delay differential system with Holling form-III based on two delay terms, cooperation/competition dynamics and the prey Allee effects is presented.

● The numerical performances of the predator-prey delay differential system with Holling form-III are provided by using the stochastic procedures of the ANNs-SCG solver.

● Three different variations based on the predator-prey delay differential system with Holling form-III are provided have been performed numerically through the ANNs-SCG solver.

● The correctness of the ANNs-SCG computing solver is performed by using the relative performances of the achieved and the reference results.

● The dependability and consistency of the ANNs-SCG computing solver is approved through the absolute error (AE) value for solving the predator-prey delay differential system with Holling form-III.

● The STs, regression values, MSE performances, EHs and correlation measures are provided using the ANNs-SCG computing for solving the predator-prey delay differential system with Holling form-III.

The paper is organized as: Section 2 shows the mathematical formulation of the predator-prey delay differential system with Holling form-III. The designed structure is given in Section 3, whereas Section 4 shows the result simulations. Concluding remarks are presented in the final section.

2. Mathematical formulations of the model

In this section, the design of the predator-prey delay differential system based Holling form-III with two delay terms, cooperation/competition dynamics and the prey Allee effects is provided. The density of the predator consists of the past and present populations of the prey. In this model, the competition delay form shows the instability, however the delay term induces the stability based on the Hopf bifurcation. The fading state based on the memory term is considered by considering the growth rate of the predator. Hence a biological model using the growth rate of the predator species is provided along with the prey density. The functional response of the predator using the predator-prey interactions has been established in four factors based on the Holling type I–IV. The Holling form III designates the functional response of the predation positive impacts that relate to the population of prey. The predation impacts increase in the enhancement of the prey population. These forms of responses are usually considered when the predator competently achieves a substitute source with low prey density. There are various studies using the Allee effects have been provided together with Holling form III and the functional response ^[42,43,44]. The predator-prey delay differential mathematical system based Holling form-III based on two delay terms is given as ^[45]:

$\left\{ \begin{array}{l} \frac{{dP(y)}}{{dy}} = \psi P(y)\left( {1 - \frac{{{P_{{y_1}}}}}{k}} \right)({P_{{y_1}}} - h) - \frac{{b{{\left( {P(y)} \right)}^2}Q(y)}}{{d + {{\left( {P(y)} \right)}^2}}} & P(0) = {\phi _1}, \hfill \\ \frac{{dQ(y)}}{{dy}} = \frac{{sbR{{(y)}^2}Q(y)}}{{d + {{\left( {R(y)} \right)}^2}}} - \tau Q(y) - m{\left( {Q(y)} \right)^2} & Q(0) = {\phi _2}, \hfill \\ \frac{{dR(y)}}{{dy}} = \frac{1}{\delta }\left( {P(y) - R(y)} \right) & R(0) = {\phi _3}, \hfill \\ \end{array} \right.$

(1)

in the above model, ${P_{{y_1}}} = P(y - {t_1})$ and ${P_{{y_2}}} = P(y - {t_2})$ represent the time delay forms of the predator-prey delay differential system and the parameter descriptions are shown in Table 1.

Table 1. Detailed parameters using the proposed ANNs-SCG procedure for the predator-prey delay differential system of Holling form-III.

Parameter	Values
$P(y)$	Prey
$Q(y)$	Predator
$R(y)$	Secluded number of characters
$d$	Rate of development from $P(y)$ to $Q(y)$
$b$	Growth of natural population
$\psi$	Development term
$\tau$	Rate of development from $P(y)$ to $R(y)$ , cure of secluded characters
$m$	Rate of development from $P(y)$ to $R(y)$ , cure of diseased characters
${t_1}$	Delay element
$y$	Period
${t_2}$	Delay element
$\delta$	Qualities' death rate
${\phi _1}, \, {\phi _2}, {\phi _3}$	Opening conditions

| Show Table

DownLoad: CSV

3. Proposed ANNs-SCG procedure

This part of the research provides the ANNs-SCG procedure for solving the predator-prey delay differential system of Holling form-III is provided by the significant operators based on the proposed stochastic ANNs-SCG and the implemented procedures of the stochastic ANNs-SCG scheme. The simulations of three different cases of the predator-prey delay differential system of Holling form-III are provided by using the ANNs-SCG procedure. Fifteen neurons have been used to solve the predator-prey delay differential system of Holling form-III along with the selection of data chosen as 13%, 12% and 75% for testing, training, and substantiation. The structure of the input, hidden and output layers is provided in Figure 1, while the optimization procedures using the multi-layer ANNs-SCG technique are given in Figure 2.

Figure 1. Input, Hidden and output layer procedure for the predator-prey delay differential system of Holling form-III.

DownLoad: Full-Size Img PowerPoint

Figure 2. Designed ANNs-SCG procedure for solving the predator-prey delay differential system of Holling form-III.

DownLoad: Full-Size Img PowerPoint

4. Simulation of the results

This section shows the numerical performances of three different cases based on the predator-prey delay differential system of Holling form-III by using the ANNs-SCG procedure. The mathematical representation of each case is provided as:

Case 1: Consider the predator-prey delay differential system of Holling form-III with $b = 0.4$ , $d = 0.14$ , $\delta = 0.21$ , ${t_1} = 0.11$ , ${t_2} = 0.15$ , $s = 0.41$ , $m = 0.7$ , $h = 0.3,$ $\tau = 0.6$ , $k = 0.1$ , $\psi = 0.2$ , ${\phi _1} = 0.15,$ ${\phi _2} = 0.2$ and ${\phi _3} = 0.25$ .

$\left\{ \begin{array}{l} \frac{{dP(y)}}{{dy}} = 0.2P(y)\left( {1 - \frac{{P(y - 0.11)}}{{0.1}}} \right)(P(y - 0.15) - 0.3) - \frac{{0.4{{\left( {P(y)} \right)}^2}Q(y)}}{{0.14 + {{\left( {P(y)} \right)}^2}}} & P(0) = 0.15, \hfill \\ \frac{{dQ(y)}}{{dy}} = \frac{{0.16{{\left( {R(y)} \right)}^2}Q(y)}}{{0.14 + {{\left( {R(y)} \right)}^2}}} - 0.6Q(y) - 0.7{\left( {Q(y)} \right)^2} & Q(0) = 0.2, \hfill \\ \frac{{dR(y)}}{{dy}} = 4.76\left( {P(y) - R(y)} \right) & R(0) = 2.5. \hfill \\ \end{array} \right.$

(3)

Case 2: Consider the predator-prey delay differential system of Holling form-III with $b = 0.4$ , $d = 0.14$ , $\delta = 0.21$ , ${t_1} = 0.11$ , ${t_2} = 0.15$ , $s = 0.41$ , $m = 0.7$ , $h = 0.3,$ , $\tau = 0.6$ , $k = 0.1$ , $\psi = 0.2$ , ${\phi _1} = 0.2,$ ${\phi _2} = 0.25$ and ${\phi _3} = 0.25$ :

$\left\{ \begin{array}{l} \frac{{dP(y)}}{{dy}} = 0.2P(y)\left( {1 - \frac{{P(y - 0.11)}}{{0.1}}} \right)(P(y - 0.15) - 0.3) - \frac{{0.4{{\left( {P(y)} \right)}^2}Q(y)}}{{0.14 + {{\left( {P(y)} \right)}^2}}} & P(0) = 0.2, \hfill \\ \frac{{dQ(y)}}{{dy}} = \frac{{0.16{{\left( {R(y)} \right)}^2}Q(y)}}{{0.14 + {{\left( {R(y)} \right)}^2}}} - 0.6Q(y) - 0.7{\left( {Q(y)} \right)^2} & Q(0) = 0.25, \hfill \\ \frac{{dR(y)}}{{dy}} = 4.76\left( {P(y) - R(y)} \right) & R(0) = 0.3. \hfill \\ \end{array} \right.$

(4)

Case 3: Consider the predator-prey delay differential system of Holling form-III with $b = 0.4$ , $d = 0.14$ , $\delta = 0.21$ , ${t_1} = 0.11$ , ${t_2} = 0.15$ , $s = 0.41$ , $m = 0.7$ , $h = 0.3,$ $\tau = 0.6$ , $k = 0.1$ , $\psi = 0.2$ , ${\phi _1} = 0.2,$ ${\phi _2} = 0.25$ and ${\phi _3} = 0.25$ :

$\left\{ \begin{array}{l} \frac{{dP(y)}}{{dy}} = 0.2P(y)\left( {1 - \frac{{P(y - 0.11)}}{{0.1}}} \right)(P(y - 0.15) - 0.3) - \frac{{0.4{{\left( {P(y)} \right)}^2}Q(y)}}{{0.14 + {{\left( {P(y)} \right)}^2}}} & P(0) = 0.25, \hfill \\ \frac{{dQ(y)}}{{dy}} = \frac{{0.16{{\left( {R(y)} \right)}^2}Q(y)}}{{0.14 + {{\left( {R(y)} \right)}^2}}} - 0.6Q(y) - 0.7{\left( {Q(y)} \right)^2} & Q(0) = 0.3, \hfill \\ \frac{{dR(y)}}{{dy}} = 4.76\left( {P(y) - R(y)} \right) & R(0) = 0.35. \hfill \\ \end{array} \right.$

(5)

The numerical results based on the predator-prey delay differential system of Holling form-III using the stochastic ANNs-SCG procedure are presented in Figures 3–5. Figure 3 indicates the gradient measures and best assessments performances for solving the predator-prey delay differential system of Holling form-III. These plots show the STs and MSE performances of the best curves, training and verification to solve three different cases of the predator-prey delay differential system of Holling form-III. The obtained performances based on the predator-prey delay differential system have been provided using the iterations 548, 48 and 39 that are calculated as 1.0883×10^‑9, 4.312×10^‑10 and 4.669×10^‑10, respectively. The gradient values are also reported in Figure 3 for solving the predator-prey delay differential system. These gradient performances are reported as 9.9916×10^-8, 9.51×10^-8 and 9.3259×10^-8 for the predator-prey delay differential system. The graphical plots represent the convergence performances of the ANNs-SCG stochastic technique to solve the predator-prey delay differential system. Figure 4 signifies the performances of the fitting cures for solving the predator-prey delay differential system. These graphical representations indicate the result comparisons for each variation of the predator-prey delay differential system. The error plots representations using the validation, testing, and training have been provided for each variation of the predator-prey delay differential system using the ANNs-SCG solver. The EHs plots together with the regression measures are provided in Figure 4 for the predator-prey delay differential system using the ANNs-SCG solver. The values of the EHs have been reported as 8.3×10^-6, 5.3×10^-6 and 5.7×10^-6 for the respective variations of the predator-prey delay differential system. The plots based on the regression measures have been derived in Figure 5 to signify the correlation values. It is noted that the correlation performances are reported as 1 for the respective variations of the predator-prey delay differential system using the ANNs-SCG solver. The authentication, testing and training performances signify the accuracy and exactness of the ANNs-SCG technique to solve the predator-prey delay differential system using the ANNs-SCG solver. The MSE measures via training, testing, verification, complexity, epochs, and backpropagation are provided in Table 2 for the predator-prey delay differential system using the ANNs-SCG solver.

Figure 3. MSE and STs for the predator-prey delay differential system using the ANNs-SCG solver.

DownLoad: Full-Size Img PowerPoint

Figure 4. EHs and result performances for the predator-prey delay differential system using the ANNs-SCG solver.

DownLoad: Full-Size Img PowerPoint

Figure 5. Regression measures based on the predator-prey delay differential system using the ANNs-SCG solver.

DownLoad: Full-Size Img PowerPoint

Table 2. Proposed ANNs-SCG technique for the predator-prey delay differential system.

Case	MSE			Gradient	Mu	Iterations	Performance	Time
Case	Testing	Training	Endorsement	Gradient	Mu	Iterations	Performance	Time
1	3×10^-9	7.16×10^-10	1.08×10^-9	9.99×10^-8	1.00×10^-10	548	7.17×10^-10	04
2	2×10^-9	1.90×10^-10	4.312×10^-10	9.51×10^-8	1.00×10^-10	48	1.90×10^-10	02
3	2×10^-11	7.36×10^-12	4.66×10^-10	9.33×10^-8	1.00×10^-9	39	7.37×10^-12	01

| Show Table

DownLoad: CSV

The result comparison performances and the values of the AE are provided for three different cases of the predator-prey delay differential system using the ANNs-SCG solver is presented in Figures 6 and 7. The precision of the results based on the stochastic ANNs-SCG technique is provided in for solving the predator-prey delay differential system of Holling form-III. The overlapping of the results is performed for solving the predator-prey delay differential system of Holling form-III. The performances of the AE based on the ANNs-SCG computing solver for three different cases of the predator-prey delay differential system of Holling form-III in signified in . The predator-prey delay differential system of Holling form-III is classified into three dynamics, $P(y)$ , $Q(y)$ and $R(y)$ . The AE for $P(y)$ is performed as 10^-5 to 10^-6, 10^-5 to 10^-7 and 10^-6 to 10^-8 for the respective cases of the predator-prey delay differential system. The AE performances for the $Q(y)$ category is presented as 10^-4 to 10^-5, 10^-3 to 10^-5 and 10^-4 to 10^-6 for the delay differential system. The AE measures for $R(y)$ class is presented as 10^-4 to 10^-5, 10^-3 to 10^-5 and 10^-4 to 10^-6 for the respective cases of the predator-prey delay differential system.

Figure 6. Results comparisons for the predator-prey delay differential system using the ANNs-SCG solver.

DownLoad: Full-Size Img PowerPoint

Figure 7. AE values for the predator-prey delay differential system using the ANNs-SCG solve.

DownLoad: Full-Size Img PowerPoint

5. Conclusions

The aim of the current investigations is to provide the numerical computing performances on the predator-prey delay differential system of Holling form-III. The predator-prey delay form of the differential system with Holling form-III is based on the two delay terms, cooperation/competition dynamics and the prey Allee effects. The competition state helps to induce the instability, whereas the cooperation state makes the stability of the system via Hopf bifurcation. Some concluding remarks of these investigations are presented as:

● A mathematical form of the nonlinear predator-prey model along with the two delay forms of the differential system with Holling form-III is based on the two delay terms, cooperation/competition dynamics and the prey Allee effects.

● The existence of the delay terms in the nonlinear predator-prey model makes more complex. Therefore, the stochastic LVMBPNNs computing procedures is a better way to provide the numerical performances.

● Fifteen numbers of neurons have been applied to solve the predator-prey delay differential system of Holling form-III together with the data selection, which is chosen as 13%, 12% and 75% for testing, training, and substantiation.

● The exactness of the LVMBPNNs computing solver is performed through the comparison of the achieved and reference numerical results.

● The AE measures are calculated in good domains that are found as 10^-5 to 10^-8 for each dynamic of the predator-prey delay differential model in Holling form-III.

Upcoming research directions: The designed stochastic procedures can be applied to solve the fluid dynamical systems ^[46,47,48], fractional kinds of differential systems ^{[49,50,51,52,53]}, and nonlinear differential systems ^[54,55,56].

Acknowledgments

The research on "Supervised neural learning for the predator-prey delay differential system of Holling form-III'' by Khon Kaen University has received funding support from the National Science, Research and Innovation Fund.

Conflict of interest

All authors declare that there are no potential conflicts of interest.

References

[1]	V. I. Arnol'd, Mathematical methods of classical mechanics, Springer Science & Business Media, 60 (2013).
[2]	R. Riaza, Differential-algebraic systems: analytical aspects and circuit applications, World Scientific, (2008).
[3]	M. M. Uddin, Computational Methods for Approximation of Large-scale Dynamical Systems, Chapman and Hall/CRC, (2019).
[4]	M. M. Rahman, M. M. Uddin, L. Andallah, M. Uddin, Tangential interpolatory projections for a class of second-order index-1 descriptor systems and application to mechatronics, Production Engineering, 15 (2021), 9–19. doi: 10.1007/s11740-020-00995-4
[5]	E. Eich-Soellner, C. Führer, Numerical methods in multibody dynamics, Springer, 45 (1998).
[6]	R. R. Craig Jr, A. J. Kurdila, Fundamentals of structural dynamics, John Wiley & Sons, (2006).
[7]	B. Yan, S. X.-D. Tan, B. McGaughy, Second-order balanced truncation for passive-order reduction of rlck circuits, IEEE Transactions on Circuits and Systems II: Express Briefs, 55 (2008), 942–946. doi: 10.1109/TCSII.2008.925655
[8]	A. C. Antoulas, Approximation of large-scale dynamical systems, SIAM, (2005).
[9]	P. Benner, V. Mehrmann, D. C. Sorensen, Dimension reduction of large-scale systems, Springer, 45 (2005).
[10]	M. M. Uddin, Computational methods for model reduction of large-scale sparse structured descriptor systems, (2015).
[11]	B. Moore, Principal component analysis in linear systems: Controllability, observability, and model reduction, IEEE transactions on automatic control, 26 (1981), 17–32. doi: 10.1109/TAC.1981.1102568
[12]	M. S. Tombs, I. Postlethwaite, Truncated balanced realization of a stable non-minimal state-space system, International Journal of Control, 46 (1987), 1319–1330. doi: 10.1080/00207178708933971
[13]	E. Grimme, Krylov projection methods for model reduction, Ph.D. dissertation, University of Illinois at Urbana Champaign, (1997).
[14]	S. Gugercin, A. C. Antoulas, C. Beattie, H_2 model reduction for large-scale linear dynamical systems, SIAM journal on matrix analysis and applications, 30 (2008), 609–638. doi: 10.1137/060666123
[15]	S. Gugercin, An iterative svd-krylov based method for model reduction of large-scale dynamical systems, Linear Algebra and its Applications, 428 (2008), 1964–1986. doi: 10.1016/j.laa.2007.10.041
[16]	F. Tisseur, K. Meerbergen, The quadratic eigenvalue problem, SIAM review, 43 (2001), 235–286. doi: 10.1137/S0036144500381988
[17]	P. Benner, J. Saak, Efficient balancing based mor for second order systems arising in control of machine tools, Proceedings of the MathMod, (2009).
[18]	M. S. Hossain, S. G. Omar, A. Tahsin, E. H. Khan, Efficient system reduction modeling of periodic control systems with application to circuit problems, in 2017 4th International Conference on Advances in Electrical Engineering (ICAEE), IEEE, (2017), 259–264.
[19]	J.-R. Li, Model reduction of large linear systems via low rank system gramians, Ph.D. dissertation, Massachusetts Institute of Technology, (2000).
[20]	S. A. Wyatt, Issues in interpolatory model reduction: Inexact solves, second-order systems and daes, Ph.D. dissertation, Virginia Tech, (2012).
[21]	M. Khatibi, H. Zargarzadeh, M. Barzegaran, Power system dynamic model reduction by means of an iterative svd-krylov model reduction method, in 2016 IEEE Power & Energy Society Innovative Smart Grid Technologies Conference (ISGT), IEEE, (2016), 1–6.
[22]	S. Li, J. Trevelyan, Z. Wu, H. Lian, D. Wang, W. Zhang, An adaptive svd–krylov reduced order model for surrogate based structural shape optimization through isogeometric boundary element method, Computer Methods in Applied Mechanics and Engineering, 349 (2019), 312–338. doi: 10.1016/j.cma.2019.02.023
[23]	A. Lu, E. L. Wachspress, Solution of lyapunov equations by alternating direction implicit iteration, Computers & Mathematics with Applications, 21 (1991), 43–58.
[24]	P. Benner, J.-R. Li, T. Penzl, Numerical solution of large-scale lyapunov equations, riccati equations, and linear-quadratic optimal control problems, Numerical Linear Algebra with Applications, 15 (2008), 755–777. doi: 10.1002/nla.622
[25]	M. M. Rahman, M. M. Uddin, L. S. Andallah, M. Uddin, Interpolatory projection techniques for $\mathcal{H}_2$ optimal structure-preserving model order reduction of second-order systems, Advances in Science, Technology and Engineering Systems Journal, 5 (2020), 715–723. doi: 10.25046/aj050485
[26]	P. Benner, M. Köhler, J. Saak, Sparse-dense sylvester equations in $h_2$ -model order reduction, (2011).

This article has been cited by:

1.	Ashraf Adnan Thirthar, Nazmul Sk, Bapin Mondal, Manar A. Alqudah, Thabet Abdeljawad, Utilizing memory effects to enhance resilience in disease-driven prey-predator systems under the influence of global warming, 2023, 69, 1598-5865, 4617, 10.1007/s12190-023-01936-x
2.	Daniel Paul, Maranya Mayengo, Salamida Daudi, Modeling the dynamics of Diamondback Moth infestations on cabbage biomass, 2023, 42, 23529148, 101338, 10.1016/j.imu.2023.101338
3.	Beyene Chanyalew, Subrata Kumar Sahu, Endalkachew Teshome Ayele, Artificial neural networks for the food web model, 2024, 139, 2190-5444, 10.1140/epjp/s13360-024-05107-0
4.	Qi Quan, Xiangjun Dai, Jianjun Jiao, Dynamics of a Predator–Prey Model with Impulsive Diffusion and Transient/Nontransient Impulsive Harvesting, 2023, 11, 2227-7390, 3254, 10.3390/math11143254

Reader Comments

Your name:*

Email:*
© 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)