Next Article in Journal
Optimal Control Systems Using Evolutionary Algorithm-Control Input Range Estimation
Previous Article in Journal
The Stability of Systems with Command Saturation, Command Delay, and State Delay
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Phase Preserving Balanced Truncation for Order Reduction of Positive Real Systems

by
Zeinab Salehi
1,
Paknoosh Karimaghaee
1,*,
Shabnam Salehi
2 and
Mohammad-Hassan Khooban
3
1
Department of Electrical Engineering, Shiraz University, Shiraz 51154-71348, Iran
2
Department of Automation and Instrumentation, Petroleum University of Technology, Tehran 14539-63341, Iran
3
Department of Electrical and Computer Engineering, Aarhus University, 8200 Aarhus, Denmark
*
Author to whom correspondence should be addressed.
Automation 2022, 3(1), 84-94; https://doi.org/10.3390/automation3010004
Submission received: 25 December 2021 / Revised: 16 January 2022 / Accepted: 17 January 2022 / Published: 19 January 2022

Abstract

:
This paper presents a new passivity-preserving order reduction method for linear time-invariant passive systems, which are also called positive real (PR) systems, with the aid of the balanced truncation (BT) method. The proposed method stems from the conic positive real balanced truncation (CPRBT) method, which is a modification of the BT method for PR systems. CPRBT presents an algorithm in which the reduced models are obtained from some Riccati equations in which the phase angle of the transfer function has been taken into consideration. Although CPRBT is a powerful algorithm for obtaining accurate PR reduced-order models, it cannot guarantee that the phase diagram of the reduced model remains inside the same interval as that of the original full-order system. We aim to address such a problem by modifying CPRBT in the way that the phase angle of the reduced transfer function always remains inside the conic and homolographic phase interval of the original system. This is proven through some matrix manipulations, which has added mathematical value to the paper. Finally, in order to assess the efficacy of the proposed method, two numerical examples are simulated.

1. Introduction

Model order reduction (MOR) plays a fundamental role in the analysis, control, and simulation of real-world systems. As most of these systems are naturally high order, a lot of computational effort is needed to analyze them. MOR tries to address this problem by simplifying the formulation of such high-order systems, creating reduced-order models with the desired accuracy. The input–output behavior of these reduced-order models is expected to be as close as possible to that of the original high-order systems. More importantly, MOR algorithms should guarantee that basic features of the original system, such as stability and passivity, remain unchanged in the reduced model.
Order reduction has been an interesting research topic since the 1980s. Therefore, a lot of order reduction methods were published during this period, some of which are balanced truncation (BT) [1], proper orthogonal decomposition (POD) [2], Hankel–Norm reduction [3], H model reduction [4], Padé Via Lanczos (PVL) [5], asymptotic waveform evaluation (AWE) [6,7], and PRIMA [8]. While these methods are designed to reduce the order of systems in all frequencies, some other methods are specifically developed to reduce the order of systems over a special frequency interval, such as Enns’ [9] and Gawronski’s [10] frequency-weighted methods. Several other model reduction procedures have been developed based on this concept including [11,12,13]. In addition, instead of reducing the order of systems in an infinite time horizon, some works consider reducing the order of systems over limited time intervals [14]. In order to find more publications in the area of model reduction, readers are referred to [15,16,17,18].
Among order reduction schemes, BT is so popular for reducing the order of linear time-invariant (LTI) systems, and it has a great reputation for preserving stability and offering an error bound. This method, however, cannot preserve some other features of the system such as passivity, bounded realness, and the negative imaginary property. In order to meet these requirements, some modifications to the BT method have been developed over the years such as bounded real balanced truncation (BRBT) [19], negative imaginary balanced truncation (NIBT) [20], positive real balanced truncation (PRBT) [21], mixed positive-bounded balanced truncation (MPBBT) [22], and conic positive real balanced truncation (CPRBT) [23].
Passive systems are a well-known class of dynamical systems that have wide applications in circuits, systems, and control theory [24,25,26]. The main feature of these systems is that they cannot produce energy, and they just store or dissipate energy [27]. Therefore, it is very important to maintain this property in the reduced models. In the case of LTI systems, passivity is also known as positive realness, and such systems are called positive real (PR) systems [27]. As mentioned before, the BT method cannot preserve passivity, and this is why some other modifications to this method have been developed which can preserve passivity, such as PRBT [21], mixed Gramian balanced truncation with error bounds [28], and CPRBT [23]. Some other passivity preserving order reduction methods are [22,29,30].
In this paper, we focus on the order reduction of PR systems whose transfer function lies inside a conic sector with an inner angle 2 θ , as indicated in Figure 1, based on balancing methods. The phase angle of the transfer function gives us valuable information about the capacitive and inductive properties of that system. Using such information contributes to achieving more accurate reduced models. The CPRBT method is a passivity-preserving balancing-based method which has used the phase angle of the transfer function in order reduction for the first time [23]. This method, however, cannot guarantee that the transfer function of the reduced model remains in the conic sector with the inner angle 2 θ displayed in Figure 1. In order to overcome this issue, this paper presents a novel order reduction scheme for PR systems that not only preserves passivity but also guarantees that the reduced model’s phase angle lies inside ( θ , θ ) as shown in Figure 1. The proposed method, named “Phase Preserving Balanced Truncation”, is inspired by the so-called CPRBT method and uses the same Riccati equations as CPRBT, but the choice of Gramians for balancing and order reduction is different. With the aid of these Gramians, a new realization is defined for PR systems, in which by selecting the states corresponding to the largest singular values, two parallel reduced models can be obtained. Separating these two parallel reduced models, we reach the final reduced model. Subsequently, it is proved via a theorem that not only the transfer function of the reduced model is passive, but it also remains in the conic sector with the inner angle 2 θ . In order to assess the efficacy of the proposed method, two numerical examples are also provided, which confirm the mathematical results.
The rest of this paper is arranged in the following sequence. The two order reduction methods, BT and CPRBT, are studied in Section 2. Then, a new order reduction method named “Phase Preserving Balanced Truncation” is introduced in Section 3. Subsequently, two numerical examples are included to assess the efficacy of the proposed method in Section 4. Finally, this paper ends with the conclusion in Section 5.
Notations: The fields of real and complex numbers are shown by symbols R and C , respectively. The set of all n × m real matrices is represented by R n × m . Furthermore, I n denotes n × n identity matrices. The transpose of matrix A is shown by A T . Moreover, the inverse of A is represented by A 1 , and the inverse of A T is indicated by A T . Moreover, A > 0 ( A 0 ) means that matrix A is positive (semi) definite.

2. Preliminaries

This section describes two balancing-based order reduction methods: the BT method and the CPRBT method. In the first place, we assume that an asymptotically stable, minimal, and LTI system is represented by
x ˙ ( t ) = A x ( t ) + B u ( t ) y ( t ) = C x ( t ) + D u ( t ) ,
where A R n × n is the state matrix, B R n × m is the input matrix, C R p × n is the output matrix, and D R p × m is the feed-forward matrix. In addition, x ( t ) R n is the state vector, u ( t ) R m is the input vector, and y ( t ) R p is the output vector of the system. The transfer function of this system is described by
H ( s ) = C ( s I A ) 1 B + D .

2.1. Balanced Truncation (BT)

The main focus of the BT method is on finding a realization in which the controllability Gramian P and the observability Gramian Q are diagonal and equal in value [1]. For asymptotically stable and minimal systems, these two positive definite (PD) Gramians are the unique solutions of
A P + P A T + B B T = 0 ,
A T Q + Q A + C T C = 0 .
The BT method finds a specific similarity transformation T R n × n such that
T T Q T = T 1 P T T = d i a g ( λ 1 , λ 2 , , λ n ) ,
where λ 1 λ 2 λ n > 0 . λ i is known as the Hankel singular value of the system for i = 1 , , n . Transforming the original system into the new coordination obtains
H ( s ) A ^ B ^ C ^ D = A ^ 11 A ^ 12 B ^ 1 A ^ 21 A ^ 22 B ^ 2 C ^ 1 C ^ 2 D .
To reach the rth order reduced model H r ( s ) , BT selects the first r states which are associated with the r largest Hankel singular values; that is,
H r ( s ) A ^ 11 B ^ 1 C ^ 1 D .
The reduced model H r ( s ) , obtained by the BT method, is asymptotically stable and minimal [1,31]. Although this method preserves asymptotic stability, it is incapable of preserving passivity.

2.2. Conic Positive Real Balanced Truncation (CPRBT)

CPRBT is a balancing-based order reduction method for PR systems which preserve passivity. This method is the first one ever which has used the phase angle of the transfer function in order reduction to reach more accurate reduced models.
Definition 1
([32]). The transfer function H ( s ) in (2) is called positive real if
m = p , D T + D 0 and H T ( j ω ) + H ( j ω ) 0 , for ω R .
Lemma 1
([23]). The transfer function of the positive real single-input single-output (SISO) system H ( s ) lies in the conic sector ( θ , θ ) , as shown in Figure 1, if and only if there exists F a = F a T > 0 R 2 n × 2 n such that
A a T F a + F a A a + F a B a C a T sin ( θ ) cos ( θ ) cos ( θ ) sin ( θ ) × N a 1 F a B a C a T sin ( θ ) cos ( θ ) cos ( θ ) sin ( θ ) T = 0 ,
where
A a B a C a D a : = A 0 B 0 0 A 0 B C 0 D 0 0 C 0 D , N a : = D a T sin ( θ ) cos ( θ ) cos ( θ ) sin ( θ ) + sin ( θ ) cos ( θ ) cos ( θ ) sin ( θ ) D a ,
and A a R 2 n × 2 n , B a R 2 n × 2 , C a R 2 × 2 n , and D a R 2 × 2 .
The dual of (7) is given by
A a E a + E a A a T + E a C a T B a sin ( θ ) cos ( θ ) cos ( θ ) sin ( θ ) × N b 1 E a C a T B a sin ( θ ) cos ( θ ) cos ( θ ) sin ( θ ) T = 0 ,
where E a = E a T > 0 R 2 n × 2 n and
N b : = D a sin ( θ ) cos ( θ ) cos ( θ ) sin ( θ ) + sin ( θ ) cos ( θ ) cos ( θ ) sin ( θ ) D a T .
We assume that N a and N b in (8) and (10) are nonsingular matrices. Equations (7) and (9) have two extremal solutions, which means any solutions F a and E a of (7) and (9), respectively, satisfy 0 < F a m i n F a F a m a x and 0 < E a m i n E a E a m a x  [23]. For the rest of this paper by F a and E a , we mean F a m i n and E a m i n , respectively.
CPRBT partitions F a and E a into four blocks as
F a = F W W T F ¯ , E a = E V V T E ¯ ,
where F , F ¯ , W R n × n and E , E ¯ , V R n × n . Following this, CPRBT balances the two Gramians F and E , which are in the upper-left corner of F a and E a , respectively, by finding a similarity transformation T R n × n such that
T T F T = T 1 E T T = Δ = d i a g Δ 1 , Δ 2 ,
where Δ 1 = d i a g ( δ 1 , , δ r ) , Δ 2 = d i a g ( δ r + 1 , , δ n ) , and δ 1 δ 2 δ n > 0 . Applying the similarity transformation T to the original system, one obtains
H ( s ) A ^ B ^ C ^ D = A ^ 11 A ^ 12 B ^ 1 A ^ 21 A ^ 22 B ^ 2 C ^ 1 C ^ 2 D .
Finally, CPRBT offers the reduced model H r ( s ) by choosing
H r ( s ) A ^ 11 B ^ 1 C ^ 1 D .
It is proved that the reduced model H r ( s ) obtained from CPRBT is positive real (passive) [23]. However, this method cannot guarantee that the phase angle of the reduced model remains inside ( θ , θ ) .
Remark 1
([23]). CPRBT obtains the same results for positive real multiple-input multiple-output (MIMO) systems. For more information please refer to Remark 4 in [23].

3. Phase Preserving Balanced Truncation

In this section, a new order reduction algorithm is presented for PR systems which not only preserves passivity but also guarantees that the phase angle of the reduced model remains inside ( θ , θ ) , as shown in Figure 1. This means that the upper bound of the absolute value of the reduced model’s phase angle would not be bigger than that of the original system’s phase angle.
For the purpose of simplification, we consider SISO systems; however, the results are applicable to MIMO systems as well. Consider Riccati Equations (7) and (9). Unlike CPRBT which looks for a similarity transformation T R n × n which balances F R n × n and E R n × n of (11), the proposed method searches for a similarity transformation T a R 2 n × 2 n which balances F a R 2 n × 2 n and E a R 2 n × 2 n of (7) and (9) such that
T a T F a T a = T a 1 E a T a T = Ψ = d i a g Ψ 1 , Ψ 2 ,
where Ψ 1 = d i a g ( ψ 1 , , ψ 2 r ) , Ψ 2 = d i a g ( ψ 2 r + 1 , , ψ 2 n ) , and ψ 1 ψ 2 ψ 2 n > 0 . We call ψ i as the phase preserving singular value of the system for i = 1 , , 2 n .
Remark 2.
Note that the above statement is valid because when we apply the similarity transformation T a to the system ( A a , B a , C a , D a ) , the Gramians F a and E a change based on the rules F a T a T F a T a and E a T a 1 E a T a T  [23]. Therefore, E a F a T a 1 E a F a T a , which means that the eigenvalues of E a F a remain invariant under similarity transformation [23]. As a consequence, there exists a similarity transformation T a which balances F a and E a . It is clear that such a similarity transformation can be obtained in the same way as the BT method.
Applying the similarity transformation T a R 2 n × 2 n to the parallel system H ( s ) I 2 , one obtains
H ( s ) I 2 = H ( s ) 0 0 H ( s ) T a 1 A a T a T a 1 B a C a T a D a = A ^ a B ^ a C ^ a D a .
Then, partitioning A ^ a , B ^ a , C ^ a based on the 2 r largest ψ i obtains
H ( s ) I 2 A ^ a B ^ a C ^ a D a = A ^ a 11 A ^ a 12 B ^ a 1 A ^ a 21 A ^ a 22 B ^ a 2 C ^ a 1 C ^ a 2 D a .
Now, we select the parallel reduced model H r ( s ) I 2 of order 2 r by selecting the states corresponding to the 2 r largest ψ i as
H r ( s ) 0 0 H r ( s ) A ^ a 11 B ^ a 1 C ^ a 1 D a .
Since H r ( s ) I 2 represents two isolated parallel systems, it is so straightforward to extract the reduced model H r ( s ) . For the sake of clarification, the algorithmic steps of the proposed method are summarized in Algorithm 1.
Algorithm 1: Phase preserving balanced truncation
1:
Consider the positive real and SISO system H ( s ) whose phase angle is inside ( θ , θ ) .
2:
Compute F a R 2 n × 2 n and E a R 2 n × 2 n by solving (7) and (9), respectively.
3:
Obtain the balancing similarity transformation T a R 2 n × 2 n such that
T a T F a T a = T a 1 E a T a T = Ψ = d i a g Ψ 1 , Ψ 2 ,
where Ψ 1 = d i a g ( ψ 1 , , ψ 2 r ) , Ψ 2 = d i a g ( ψ 2 r + 1 , , ψ 2 n ) , and ψ 1 ψ 2 ψ 2 n > 0 .
4:
Find the state space form of H ( s ) I 2 in the new coordination as
H ( s ) I 2 = H ( s ) 0 0 H ( s ) T a 1 A a T a T a 1 B a C a T a D a = A ^ a B ^ a C ^ a D a .
5:
Partition A ^ a , B ^ a , C ^ a based on the 2 r largest ψ i as
H ( s ) I 2 A ^ a B ^ a C ^ a D a = A ^ a 11 A ^ a 12 B ^ a 1 A ^ a 21 A ^ a 22 B ^ a 2 C ^ a 1 C ^ a 2 D a .
6:
Select H r ( s ) I 2 as
H r ( s ) 0 0 H r ( s ) A ^ a 11 B ^ a 1 C ^ a 1 D a .
7:
As H r ( s ) I 2 represents two separate systems, extract the rth order H r ( s ) from (19) by eliminating the redundant states.
Theorem 1.
The phase diagram of the reduced model H r ( s ) obtained from the proposed algorithm always remains inside ( θ , θ ) .
Proof. 
Assume that the system H ( s ) I 2 is in the balanced realization (17). The Riccati equation (7) in the new coordination would be
A ^ a T Ψ + Ψ A ^ a + Ψ B ^ a C ^ a T sin ( θ ) cos ( θ ) cos ( θ ) sin ( θ ) × N a 1 Ψ B ^ a C ^ a T sin ( θ ) cos ( θ ) cos ( θ ) sin ( θ ) T = 0 .
Considering Ψ = d i a g ( Ψ 1 , Ψ 2 ) , we can extract the ( 1 , 1 ) block of (20) as
A ^ a 11 T Ψ 1 + Ψ 1 A ^ a 11 + Ψ 1 B ^ a 1 C ^ a 1 T sin ( θ ) cos ( θ ) cos ( θ ) sin ( θ ) × N a 1 Ψ 1 B ^ a 1 C ^ a 1 T sin ( θ ) cos ( θ ) cos ( θ ) sin ( θ ) T = 0 .
Since ( A ^ a 11 , B ^ a 1 , C ^ a 1 , D a ) represents two separate parallel systems, there exists a similarity transformation T ¯ a R 2 r × 2 r which can extract these two systems as
A ¯ a B ¯ a C ¯ a D a : = T ¯ a 1 A ^ a 11 T ¯ a T ¯ a 1 B ^ a 1 C ^ a 1 T ¯ a D a = A r 0 B r 0 0 A r 0 B r C r 0 D 0 0 C r 0 D .
Applying T ¯ a to the system ( A ^ a 11 , B ^ a 1 , C ^ a 1 , D a ) also yields
A ¯ a T Ψ ¯ + Ψ ¯ A ¯ a + Ψ ¯ B ¯ a C ¯ a T sin ( θ ) cos ( θ ) cos ( θ ) sin ( θ ) × N a 1 Ψ ¯ B ¯ a C ¯ a T sin ( θ ) cos ( θ ) cos ( θ ) sin ( θ ) T = 0 ,
where Ψ ¯ = T ¯ a T Ψ 1 T ¯ a . In the meantime, because of the fact that Ψ 1 > 0 , we yield Ψ ¯ > 0 . As a consequence, based on the Riccati equation (23), we can say Lemma 1 holds for the reduced-order model H r ( s ) with the state-space form
H r ( s ) A r B r C r D .
This completes the proof. □
Corollary 1.
The reduced model H r ( s ) obtained from the proposed algorithm preserves passivity.
Remark 3.
The same results can be held for positive real MIMO systems although we have only studied SISO systems for simplification. To this end, we have to use the associated Riccati equations for MIMO systems which are defined in [23] instead of (7) and (9) and then use the same procedure as the proposed algorithm to obtain the reduced model.

4. Illustrative Examples

This section provides two illustrative examples to investigate the efficacy of the proposed method.

4.1. Example 1

In this example, a 10th order PR system is examined whose state-space form is
A = b l o c k d i a g k = 1 , , 5 0 1 π k 2 2 μ π k , π k = k 2 , k = 1 , , 5 , μ = 0.01 , C = B T = [ 0 0.9877 0 0.309 0 0.891 0 0.5878 0 0.7071 ] , D = 0.2 .
We want to reach the reduced model of order r = 4 utilizing two methods: CPRBT and the proposed method (phase preserving balanced truncation). As the phase angle of the transfer function is inside ( 83 , 83 ) , we chose θ = 83 in both methods. The state-space form of the 4th order reduced model obtained by the proposed method is given by
A r B r C r D = 0.0423 2.6976 0.0036 0.0151 1.0438 2.6027 0.0545 1.3232 0.0217 1.2909 0.0399 1.3668 0.0104 1.1544 0.2793 0.0389 0.0508 1.1624 0.0122 1.1533 0.5282 0.6426 0.1335 0.5749 0.2 .
The bode diagrams of the reduced-order models, along with the full-order system, are depicted in Figure 2.
Comparing the phase angle of the reduced models, we observe that the phase angle of the proposed reduced model remains inside ( 83 , 83 ) , while the phase angle of the reduced model obtained by CPRBT exceeds this bound and comes inside ( 83.145 , 83.145 ) . This clearly shows the efficacy of the proposed method. The two reduced models approximately have the same levels of accuracy in terms of the input–output behavior of the system, and their magnitude in the bode plot almost matches each other. Figure 2 completely confirms the results.

4.2. Example 2

Here, we provide an example which illustrates the effectiveness of the proposed method regarding accuracy in the input–output behavior of the system. Consider the 6th order PR system with the state-space form
A = 0 0 0 100 100 0 0 0 0 0 100 0 0 0 110 0 0 100 100 0 0 10 0 100 100 100 0 0 10 0 0 0 100 100 0 0 , B = 0 0 100 0 0 0 T , C = 0 0 1 0 0 0 , D = 1 .
The phase angle of the above system is inside ( 18.1 , 18.1 ) . Therefore, we choose θ = 18.1 for the proposed method and CPRBT. We also select r = 4 to reach the 4th order reduced models. The reduced model obtained by the proposed method is given by
A r B r C r D = 100.34 108.79 14.35 23.06 13.22 111.56 3.42 59.26 12.07 0.65 62.33 64.26 7.20 32.65 3.37 50.91 13.71 48.68 14.16 4.67 6.40 0.18 1.52 2.50 1 .
The bode diagrams of the reduced models, along with the full-order system, are illustrated in Figure 3. In addition, the frequency response of the error system H ( s ) H r ( s ) , obtained from CPRBT and the proposed method, is depicted in Figure 4 which enables us to compare the results precisely.
As shown in Figure 3 and Figure 4, the proposed method creates a more accurate reduced model than CPRBT, although both methods provide reduced models whose phase angles remain inside ( 18.1 , 18.1 ) . The H norm of the error system H ( s ) H r ( s ) is 0.3714 for the CPRBT reduced model and 0.3685 for the proposed reduced model, which confirms that the proposed method has offered a better reduced-model than CPRBT.

5. Conclusions

This paper proposes a new passivity-preserving order reduction method for PR systems whose phase angle is inside the interval ( θ , θ ) . The proposed method is a modification to the CPRBT method. What makes this method different from CPRBT is that it proves the phase angle of the reduced models remains inside the same interval as that of the original system, in contrast to CPRBT. In addition, in some cases, the proposed method can produce more accurate reduced models than CPRBT. In order to assess this method, two numerical examples are simulated. Example 1 perfectly verifies that the phase angle of the reduced model obtained by the proposed method remains inside ( θ , θ ) , while the reduced model obtained by CPRBT exceeds this interval. In addition, Example 2 shows that the proposed method can even find reduced models of higher accuracy than CPRBT. These two points show the effectiveness of the proposed method. Although the provided examples are simple, the proposed method can be applied to high-dimensional systems as well. As future work, it is of interest to find a computational H error bound for the proposed method. It is also suggested to extend this method to discrete-time, nonlinear, and descriptor systems. Furthermore, the robustness of the proposed method can be the topic of future research.

Author Contributions

Conceptualization, Z.S., P.K., S.S., and M.-H.K.; methodology, Z.S., P.K., S.S., and M.-H.K.; software, Z.S.; validation, Z.S. and P.K.; formal analysis, P.K. and M.-H.K.; investigation, Z.S. and S.S.; resources, Z.S., P.K., S.S., and M.-H.K.; writing—original draft preparation, Z.S. and S.S.; writing—review and editing, P.K. and M.-H.K.; visualization, Z.S. and P.K.; supervision, P.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Moore, B. Principal component analysis in linear systems: Controllability, observability, and model reduction. IEEE Trans. Autom. Control 1981, 26, 17–32. [Google Scholar] [CrossRef]
  2. Sirovich, L. Turbulence and the dynamics of coherent structures. I. Coherent structures. Q. Appl. Math. 1987, 45, 561–571. [Google Scholar] [CrossRef] [Green Version]
  3. Glover, K. All optimal Hankel-norm approximations of linear multivariable systems and their L error bounds. Int. J. Control 1984, 39, 1115–1193. [Google Scholar] [CrossRef]
  4. Kavranoǧlu, D.; Bettayeb, M. Characterization of the solution to the optimal H model reduction problem. Syst. Control Lett. 1993, 20, 99–107. [Google Scholar] [CrossRef]
  5. Feldmann, P.; Freund, R.W. Efficient linear circuit analysis by padé approximation via the Lanczos process. IEEE Trans. Comput. Aided Des. Integr. Circuits Syst. 1995, 14, 639–649. [Google Scholar] [CrossRef]
  6. Pillage, L.T.; Rohrer, R.A. Asymptotic waveform evaluation for timing analysis. IEEE Trans. Comput. Aided Des. Integr. Circuits Syst. 1990, 9, 352–366. [Google Scholar] [CrossRef] [Green Version]
  7. Chiprout, E.; Nakhla, M.S. Asymptotic Waveform Evaluation and Moment Matching for Interconnect Analysis; Kluwer: Norwell, MT, USA, 1994. [Google Scholar]
  8. Odabasioglu, A.; Celik, M.; Pileggi, L.T. PRIMA: Passive reduced-order interconnect macromodeling algorithm. IEEE Trans. Comput. Aided Des. Integr. Circuits Syst. 1998, 17, 645–654. [Google Scholar] [CrossRef] [Green Version]
  9. Enns, D.F. Model reduction with balanced realizations: An error bound and a frequency weighted generalization. In Proceedings of the 23rd IEEE Conference on Decision and Control, Las Vegas Hilton, LV, USA, 12–14 December 1984; pp. 127–132. [Google Scholar]
  10. Gawronski, W.; Juang, J.-N. Model reduction in limited time and frequency intervals. Int. J. Syst. Sci. 1990, 21, 349–376. [Google Scholar] [CrossRef]
  11. Wang, G.; Sreeram, V.; Liu, W. A new frequency-weighted balanced truncation method and an error bound. IEEE Trans. Autom. Control 1999, 44, 1734–1737. [Google Scholar] [CrossRef] [Green Version]
  12. Kumar, D.; Sreeram, V. Factorization-based frequency-weighted optimal Hankel-norm model reduction. Asian J. Control 2020, 22, 2106–2118. [Google Scholar] [CrossRef]
  13. Toor, H.I.; Imran, M.; Ghafoor, A.; Kumar, D.; Sreeram, V.; Rauf, A. Frequency limited model reduction techniques for discrete-time systems. IEEE Trans. Circuits Syst. II Express Briefs 2019, 67, 345–349. [Google Scholar] [CrossRef]
  14. Haider, S.; Ghafoor, A.; Imran, M.; Malik, F.M. Time-limited Gramians-based model order reduction for second-order form systems. Trans. Inst. Meas. Control 2019, 41, 2310–2318. [Google Scholar] [CrossRef]
  15. Liang, W.; Chen, H.-B.; He, G.; Chen, J. Model order reduction based on dynamic relative gain array for mimo systems. IEEE Trans. Circuits Syst. II Express Briefs 2020, 67, 2507–2511. [Google Scholar] [CrossRef]
  16. Alsmadi, O.; Al-Smadi, A.; Gharaibeh, E. Firefly artificial intelligence technique for model order reduction with substructure preservation. Trans. Inst. Meas. Control 2019, 10, 2875–2885. [Google Scholar] [CrossRef]
  17. Lu, K.; Jin, Y.; Huang, P.; Zhang, F.; Zhang, H.; Fu, C.; Chen, Y. The applications of POD method in dual rotor-bearing systems with coupling misalignment. Mech. Syst. Signal Process. 2021, 150, 107236. [Google Scholar] [CrossRef]
  18. Bui-Thanh, T.; Willcox, K.; Ghattas, O. Model reduction for large-scale systems with high-dimensional parametric input space. SIAM J. Sci. Comput. 2008, 30, 3270–3288. [Google Scholar] [CrossRef]
  19. Opdenacker, P.C.; Jonckheere, E.A. A contraction mapping preserving balanced reduction scheme and its infinity norm error bounds. IEEE Trans. Circuits Syst. 1988, 35, 184–189. [Google Scholar] [CrossRef]
  20. Tu, K.; Du, X.; Fan, P. Negative imaginary balancing for mode reduction of LTI negative-imaginary systems. In Proceedings of the 26th Chinese Control and Decision Conference, Changsha, China, 31 May–2 June 2014; pp. 4234–4239. [Google Scholar]
  21. Phillips, J.R.; Daniel, L.; Silveira, L.M. Guaranteed passive balancing transformations for model order reduction. IEEE Trans. Comput. Aided Des. Integr. Circuits Syst. 2003, 22, 1027–1041. [Google Scholar] [CrossRef]
  22. Salehi, Z.; Karimaghaee, P.; Khooban, M.-H. Mixed positive-bounded balanced truncation. IEEE Trans. Circuits Syst. Ii: Express Briefs 2021, 68, 2488–2492. [Google Scholar] [CrossRef]
  23. Salehi, Z.; Karimaghaee, P.; Khooban, M.-H. A new passivity preserving model order reduction method: Conic positive real balanced truncation method. IEEE Trans. Syst. Man Cybern. Syst. 2021. [Google Scholar] [CrossRef]
  24. Huang, C.H.; Ioannou, P.A.; Maroulas, J.; Safonov, M.G. Design of strictly positive real systems using constant output feedback. IEEE Trans. Autom. Control 1999, 44, 569–573. [Google Scholar] [CrossRef]
  25. Liu, M.; Lam, J.; Zhu, B.; Kwok, K.W. On positive realness, negative imaginariness, and H control of state-space symmetric systems. Automatica 2019, 101, 190–196. [Google Scholar] [CrossRef]
  26. Misgeld, B.J.; Hewing, L.; Liu, L.; Leonhardt, S. Closed-loop positive real optimal control of variable stiffness actuators. Control Eng. Pract. 2019, 82, 142–150. [Google Scholar] [CrossRef]
  27. Brogliato, B.; Lozano, R.; Maschke, B.; Egeland, O. Positive Real Systems. In Dissipative Systems Analysis and Control; Springer: Cham, Switzerland, 2020. [Google Scholar]
  28. Salehi, Z.; Karimaghaee, P.; Khooban, M.-H. Model order reduction of positive real systems based on mixed gramian balanced truncation with error bounds. Circuits Syst. Signal Process. 2021, 40, 5309–5327. [Google Scholar] [CrossRef]
  29. Cheng, X.; Scherpen, J.M.; Besselink, B. Balanced truncation of networked linear passive systems. Automatica 2019, 104, 17–25. [Google Scholar] [CrossRef] [Green Version]
  30. Zulfiqar, U.; Imran, M.; Ghafoor, A.; Liaqat, M. Time/frequency-limited positive-real truncated balanced realizations. IMA J. Math. Control Inf. 2020, 37, 64–81. [Google Scholar] [CrossRef]
  31. Pernebo, L.; Silverman, L. Model reduction via balanced state space representations. IEEE Trans. Autom. Control 1982, 27, 382–387. [Google Scholar] [CrossRef]
  32. Anderson, B.D.; Vongpanitlerd, S. Network Analysis and Synthesis: A Modern Systems Theory Approach; Dover: New York, NY, USA, 2013. [Google Scholar]
Figure 1. Conic sector with the inner angle 2 θ .
Figure 1. Conic sector with the inner angle 2 θ .
Automation 03 00004 g001
Figure 2. The bode diagrams of the original system, the CPRBT reduced model, and the proposed reduced model in Example 1.
Figure 2. The bode diagrams of the original system, the CPRBT reduced model, and the proposed reduced model in Example 1.
Automation 03 00004 g002
Figure 3. The bode diagrams of the original system, the CPRBT reduced model, and the proposed reduced model in Example 2.
Figure 3. The bode diagrams of the original system, the CPRBT reduced model, and the proposed reduced model in Example 2.
Automation 03 00004 g003
Figure 4. The frequency response of the error system H ( s ) H r ( s ) obtained from CPRBT and the proposed method in Example 2.
Figure 4. The frequency response of the error system H ( s ) H r ( s ) obtained from CPRBT and the proposed method in Example 2.
Automation 03 00004 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Salehi, Z.; Karimaghaee, P.; Salehi, S.; Khooban, M.-H. Phase Preserving Balanced Truncation for Order Reduction of Positive Real Systems. Automation 2022, 3, 84-94. https://doi.org/10.3390/automation3010004

AMA Style

Salehi Z, Karimaghaee P, Salehi S, Khooban M-H. Phase Preserving Balanced Truncation for Order Reduction of Positive Real Systems. Automation. 2022; 3(1):84-94. https://doi.org/10.3390/automation3010004

Chicago/Turabian Style

Salehi, Zeinab, Paknoosh Karimaghaee, Shabnam Salehi, and Mohammad-Hassan Khooban. 2022. "Phase Preserving Balanced Truncation for Order Reduction of Positive Real Systems" Automation 3, no. 1: 84-94. https://doi.org/10.3390/automation3010004

APA Style

Salehi, Z., Karimaghaee, P., Salehi, S., & Khooban, M. -H. (2022). Phase Preserving Balanced Truncation for Order Reduction of Positive Real Systems. Automation, 3(1), 84-94. https://doi.org/10.3390/automation3010004

Article Metrics

Back to TopTop