Mathematical Problems in Engineering

Volume 2014 (2014), Article ID 765858, 11 pages

http://dx.doi.org/10.1155/2014/765858

Research Article

Complex-Coefficient Frequency Domain Stability Analysis Method for a Class of Cross-Coupled Antisymmetrical Systems and Its Extension in MSR Systems

Yuan Ren1,2 and Jiancheng Fang²

¹Department of Space Equipment, Equipment Academy, Beijing 101416, China

²School of Instrumentation Science & Opto-electronics Engineering, Beihang University, Beijing 100191, China

Received 25 July 2013; Revised 13 December 2013; Accepted 16 December 2013; Published 20 February 2014

Academic Editor: Asier Ibeas

Copyright © 2014 Yuan Ren and Jiancheng Fang. 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.

Abstract.

This paper develops a complex-coefficient frequency domain stability analysis method for a class of cross-coupled two-dimensional antisymmetrical systems, which can greatly simplify the stability analysis of the multiple-input multiple-output (MIMO) system. Through variable reconstruction, the multiple-input multiple-output (MIMO) system is converted into a single-input single-output (SISO) system with complex coefficients. The pole locations law of the closed-loop system after the variable reconstruction has been revealed, and the controllability as well as observability of the controlled plants before and after the variable reconstruction has been studied too, and then the classical Nyquist stability criterion is extended to the complex-coefficient frequency domain. Combined with the rigid magnetically suspended rotor (MSR) system with heavy gyroscopic effects, corresponding stability criterion has been further developed. Compared with the existing methods, the developed criterion for the rigid MSR system not only accurately predicts the absolute stability of the different whirling modes, but also directly demonstrates their relative stability, which greatly simplifies the analysis, design, and debugging of the control system.

1. Introduction

The cross-coupling effect exists widely in many systems in practice, for example, the magnetically suspended rotor (MSR) systems [1–4], spacecraft attitude control systems [4, 5], power electronic systems [6–9], optical systems [10, 11], and so on [12].

Commonly, the cross-coupled systems can be classified into two kinds: the cross-coupled symmetrical and antisymmetrical systems [13, 14]. As for the controlled plant with antisymmetrical cross-coupling characteristics, antisymmetrical cross-decoupling controller is often needed to realize the decoupling control of the system. This paper focuses on the stability analysis of the cross-coupled two-dimensional antisymmetrical system. For the sake of simplicity and clarity, the whole control system including the cross-coupled antisymmetrical controlled plant and its antisymmetrical cross-decoupling controller is called the cross-coupled antisymmetrical system in this paper.

Cross-coupled antisymmetrical system is a class of special multi-input multi-output (MIMO) systems. As for the stability analysis of the MIMO system, the most traditional method resorts to resolve the characteristic roots [15]. However, as for the high-order cross-coupled antisymmetrical system, it is very tedious and even impossible to obtain its analytical expressions. Although this issue may be alleviated by employing the Routh-Hurwitz criterion [16–18], this method cannot present the whole information about the system stability. Additionally, rebuilding a proper Lyapunov function is not easy since there is no rule to follow [19], especially for a high-order MIMO system. Therefore, when employed to judge the stability of the cross-coupled antisymmetrical system, these traditional methods are relatively complex. Moreover, the resolved relative stability information is abstract since they lack specific physical meanings. As a result, these methods have some difficulty in combining the stability analysis with the controller design.

Accordingly, it is quite interesting and significative to judge directly and simply the stability of the cross-coupled antisymmetrical system, which is also a basic theoretical issue to be urgently addressed in engineering applications. An ideal stability criterion must be combined tightly with the design process of the controller, which can effectively guide the system design and further reveal the influence of the variation of the system parameters on the system stability [20]. Obviously, the above methods can not satisfy this requirement, and it will be more difficult when considering simultaneously the variations of the multiple system parameters.

The Nyquist stability criterion of the single-variable system has many advantages [21–24]. It not only can judge the closed-loop stability, but also can directly reveal the relationships between the relative stability of the system and its gain. Although this method has been extended to the multivariable case [25], such as characteristic locus method [26–31] and the inverse Nyquist array (INA) method [32, 33], these extended methods are quite complex since they depend on some special computer-aided-design (CAD) softwares. Especially when diagonal dominance does not hold, the stability of the MIMO system cannot be judged by the INA method.

Among all the cross-coupled antisymmetrical systems, the MSR system is a class of typical cases in which stability analysis is quite important since the MSR system is an open-loop unstable system with high rotor speed [34–36]. Especially for the MSR system with significant gyroscopic effects, its whirling modes stability analysis is more challenging. However, the existing stability analysis methods [37–40] can only judge the stability of the whole whirling modes and cannot present the stability information of the different whirling modes, which is the key concern of the controller design and debugging. To solve this issue, the whirling modes stability criterion [41] based on complex coefficient frequency characteristics has been proposed, which can judge the stability of the different whirling modes. This method is of great value to the development of stability analysis of control theory. However, this method is only suited to the special MSR system, and it does not suit a class of cross-coupled antisymmetrical systems. How to extend the method to a class of cross-coupled systems is a great interest in engineering.

To simplify the absolute stability judgment of the cross-coupled antisymmetrical system and reveal directly its relative stability, a complex-coefficient frequency domain stability analysis method for a class of cross-coupled antisymmetrical systems is presented in this paper.

2. Complex-Coefficient Frequency Domain Stability Analysis Method

2.1. Variable Reconstruction of the Cross-Coupled Antisymmetric System

Figure 1 shows the structural scheme of the cross-coupled two-dimensional antisymmetrical system, where , , . Without loss of generality, the standard type of the cross-coupled antisymmetrical controlled plant can be written as

where , , and .

Figure 1: Structural scheme of the cross-coupled two-dimensional antisymmetrical system.

According to Figure 1, we assume that the state vector of the controlled plant is , where . The variable reconstruction is performed as follows: , , , . Then, after the variable reconstruction, the original MIMO system with real coefficients is converted into a SISO control system with complex coefficients, which is shown in Figure 2.

Figure 2: Structural scheme of the SISO control system with complex coefficients.

Based on the information of the original controlled plant, the system realization of the SISO controlled plant can be written as

And the open-loop transfer function of the SISO system is given by

2.2. Pole Locations Law before and after Variable Reconstruction

Theorem 1. Suppose that the linear differential equations have the form

where is the differential operator, , and . If all the solutions for the above satisfy the following differential equations:

where , then the eigenvalues of the and satisfy the following relationships:

(1)If , then ;(2)If and , then .

Proof. Assume that the basic solution set of is , according to the theory of the linear differential equations, has one of the following three forms: , , and , where is a constant real vector and .

At the same time, can also be given by

As for the complex-coefficient differential equations , its basic solution set can be described as , where has a form of , where is the constant complex vector and .

Define the following function space:

From the proposition, it can be obtained that . According to the property of the general solution for the ordinary differential equations, one can obtain that . Accordingly, .

As for any , , it can be proven that , , and are linearly independent if . When , is the corresponding eigenvector; according to the format of the solution to the ordinary differential equations, is the solution of (5); that is, . When , is the corresponding eigenvector, according to the format of the solution to the ordinary differential equations, is the solution to (5), that is, . If , from the format of the fundamental-solution matrix for (4), it can be derived that since any and are linearly independent and is the set of the base vectors of . And this contradicts the fact that . Accordingly, proposition is right.

Simultaneously, if , then both and are the solutions for (4). One can define , where ; then one can obtain that . If and , then it can be drawn that , which is not in accordance with . Accordingly, proposition is positive. This completes the proof.

Theorem 2. The systems before and after variable reconstruction satisfy the pole locations law presented in Theorem 1.

Proof. Note that system matrix determines the pole locations; therefore, it should be obtained at first. Since , their system realization is given by

And the state variables of and are , , respectively.

Similarly the system realization of and is given by

And the state variables of and are, respectively, and .

Also, the system realization of and can be described as

with the state variables of and , respectively.

Define ; then . According to the closed-loop structure diagram of the control system, one has

Then the system matrix of the whole system before variable reconstruction is given by

Define the new state variables ; then, according to (12) and Theorem 1, the corresponding system matrix of can be derived as

As for the SISO system with complex coefficients, the system realization of can be derived from the system realizations of and as follows:

Similarly, the complex-coefficient system matrix of the SISO system after variable reconstruction can be derived as

From (13) and (15), it can be drawn that . Accordingly, the original system and the new system are in line with the same pole locations laws with that of Theorem 1. This concludes the proof.

2.3. Controllability and Observability before and after Variable Reconstruction

Lemma 3. If , then , where , , , .

Proof. When , the conclusion holds obviously. Suppose that the conclusion also holds when ; then when , ; that is, and .

Similarly, . This concludes the proof.

Theorem 4. If the SISO system (2) is controllable and observable, that is, and , where , then the MIMO system (1) is controllable and observable too.

Proof. Denote , where . If the MIMO system (1) is not controllable, then

That is, there exists a set of real numbers which are not all zeros so that

According to Lemma 3, it can be obtained that . Suppose ; then, from (17), it can be drawn that , which is not in accordance with the assumption.

Similarly, it can be proven that if the SISO system (2) is observable, then the MIMO system (1) is observable too. Proof of the proposition is completed.

2.4. Complex-Coefficient Frequency Domain Stability Criterion

According to Theorems 1–4, the following corollary can be achieved.

Corollary 5. As for the complex-coefficient system (15), the sufficient and necessary condition for the system stability is that all the characteristic roots lie in the left half plane (LHP). And the sufficient and necessary condition for the system critical stability is that there is at least a pure imaginary root and there is no right-half plane (RHP) pole.

From the analysis above, it can be drawn that the stability of the cross-coupled antisymmetrical system can be demonstrated by that of the corresponding complex-coefficient system. Note that the argument principle applies to the complex variable function which is defined on the complex plane, and the complex-coefficient open-loop transfer function also belongs to this class of complex variable functions. Since the Nyquist stability criterion is based on the argument principle and the system controllability and observability are the precondition of Nyquist stability criterion, the corresponding SISO system can be employed to analyze the stability of the original control system by using complex-coefficient Nyquist stability criterion. Accordingly, one can achieve the following stability analysis method based on complex-coefficient frequency characteristics. Denote as the number of anticlockwise encirclements of the point by the locus as a representative point traces on the path in the clockwise direction, where is the semicircle contour that encloses the entire right-half plane. Then one can obtain the following stability criterion for two-dimensional cross-coupled antisymmetrical systems.

Criterion 1. As for the cross-coupled two-dimensional antisymmetrical system (12), if its corresponding SISO controlled plant (2) is both controllable and observable, and its number of the RHP poles is , then the sufficient and necessary condition for the stability of system (12) is that .

Compared with the traditional Nyquist stability criterion for the multi variable systems, the proposed one greatly simplified the stability analysis since the complex multivariable system has become a single-variable system.

In fact, as for a special controlled plant, this criterion can be further developed according to its practical physical meaning, for example, the MSR system with heavy gyroscopic effects; it can be proven that the closed-loop right-half-poles and left-half-poles of the complex-coefficient SISO system denote the forward and backward whirling modes of the MSR respectively, which is the theoretical base for the whirling modes stability criterion [41].

3. Simulation and Experimental Results

Note that the large magnetically suspended control moment gyroscope (MSCMG) has a rigid MSR with significant cross-coupling effects which is a typical coupled-coupled antisymmetrical system; thus, it is taken as the example. To demonstrate the effectiveness and superiority of the proposed stability criterion, comparative simulation between the proposed method and the traditional one (root locus with respect to rotor speed) and experiments have been performed.

3.1. Experiment Setup

Figure 3 shows the experimental setup, where the MSR has two 2-degree-of-freedom (DOF) radial MBs and two single-DOF axial MBs, which are driven by five H-bridge unipolar switching power amplifiers independently [42], realizing 5-DOF active control.

Figure 3: Photograph of the experimental setup.

The well-known decentralized PID plus filtered cross feedback control approach [34] is adopted. The main parameters of the MSCMG and the controller coefficients used in simulation and experiments have been listed in Tables 1 and 2.

Table 1: System parameters of the MSCMG.

Parameter
Value

56 kg

0.6032 kg·m²

0.7958 kg·m²

1140 N/A

75926

120 m

0.6032 kg·m²

0.177 m

2.83 N/m

0.113 m

Table 2: Parameters of the controller.

Parameter
Value

2.0

0.167

0.0008

0.00015

280 Hz

10.0

1.7

0.0045

1.0

30 Hz

The experiments are performed in a TMS320C32 digital signal processor (DSP) board. The total digital control system delay can be tested, which is about 150 μs. To facilitate drawing the root locus diagrams [43] and compare fairly, the delay unit is simplified as its second-order Taylor's series approximation when drawing the root locus and the complex-coefficient frequency characteristics plots; that is,

where . In fact, this approximate has little influence on the stability analysis for the MSR system since the interested critical stability frequencies of the whirling modes are relatively low. That is, the phase errors at the critical stability frequencies are tolerant.

The closed-loop rotation dynamic model with the antisymmetrical cross-decoupling controller can be given by [41, 44]

where and are the momenta of inertia of the rotor about the radial and axial directions, respectively; is the rotor speed; and are the distances from the geometry center of the MSR to the radial and the radial sensor; and are the rotor angular displacements about the - and -axes; , , , , and denote the basic controller, the antisymmetrical cross-decoupling controller, the switching power amplifier, the antialias filter, and the control system delay unit, respectively; is the force-displacement factor of the radial magnetic bearings, and is the proportional gain of the displacement sensor.

The open-loop transfer function of the SISO system with complex coefficients can be derived as [44]

According to the definition used by logarithmic frequency stability criterion [45], we make the following definitions: , denote the number of the positive and negative phase cross-over points of in the set of , when changes from 0 to ; , denote the number of the positive and negative phase cross-over points of in , when decreases from to ; , are the number of the positive and negative phase cross-over points of in , when increases from 0 to ; and , denote the number of the positive and negative phase cross-over points of in , when changes from to . According to the proposed Criterion 1 and the whirling modes stability criterion for a MSR with bending modes [41], one can further obtain the following whirling modes stability criterion for the rigid MSR without bending modes.

Criterion 2. As for the rigid MSR system (19), the sufficient and necessary conditions for the nutation and precession modes stability are and , respectively.

3.2. Simulation and Experimental Results

At first, under the conditions that the cross proportion is 48 Hz, the Bode diagram is drawn, respectively, at the rotation speed of 140 Hz and 145 Hz, which is shown in Figure 4. Substituting the system parameters and rotor speed into (20), it can be calculated that and . From Figure 4(a), it can be achieved that and . According to the whirling modes stability criterion, , ; that is, both the forward and backward whirling modes are stable. Simultaneously, from Figure 4(b), it can be achieved that and  − . Therefore, the backward whirling mode is stable, while the forward whirling mode is unstable at the speed of 145 Hz. The root locus with respect to the rotor speed is drawn, which is shown in Figure 5. The speed range is from 0 Hz to 150 Hz and it is plotted every 5 Hz. According to Figure 5, it can be derived that the system is stable at the speed of 140 Hz while unstable at the speed of 145 Hz; that is, the margin speed is located between 140 Hz and 145 Hz, which is in accordance with the conclusion drawn from the proposed stability criterion. It is necessary to note that phase margin of the nutation and precession modes can be obtained from the proposed stability criterion, while it cannot be obtained from the root locus.

(a)

(b)

Figure 4: Complex-coefficient frequency characteristics of the SISO system with 48 Hz cross proportion and different rotation speeds. (a) 140 Hz. (b) 145 Hz.

Figure 5: Root locus with respect to rotor speed from 0 Hz to 150 Hz. Note: the arrowhead denotes the increased direction of rotor speed from 0 to 150 Hz.

Then we set the cross proportion to be 150 Hz, the corresponding simulation and experiments have been developed too, and their results are shown in Figures 6 and 7. Also, it can be obtained that and at the speeds of 200 Hz and 205 Hz and and at the speeds of 55 Hz and 50 Hz. From Figure 6(a), one can obtain that , . According to the proposed criterion, it can be drawn that both the forward and backward whirling modes are stable at the speed of 200 Hz. From Figure 6(b), one can obtain that and ; that is, the backward whirling mode remains stable while the forward whirling mode loses its stability at the speed of 205 Hz. Similarly, from Figures 6(c) and 6(d), the conclusion can be achieved that both the forward and backward modes are stable at the speed of 55 Hz while the backward whirling mode is unstable at the speed of 50 Hz. Also, from Figures 6(c) and 6(d), it can be calculated that the phase margins of the forward whirling mode at the speed of 50 Hz and 55 Hz are 56 degree and 52 degrees, respectively. Usually the rise of the rotation speed will result in the decrease of the phase margins of the forward whirling mode, and consequently the coefficient of the controller should be adjusted with the rotor speed in order to guarantee the adequate phase margins of the different whirling modes in the expected speed range.

(a)

(b)

(c)

(d)

Figure 6: Complex-coefficient Bode diagram of the SISO system with 150 Hz cross proportion and different rotation speeds. (a) 200 Hz. (b) 205 Hz. (c) 55 Hz. (d) 50 Hz.

Figure 7: Root locus with respect to rotor speed from 0 to 220 Hz. Note: the arrowhead denotes the increased direction of rotor speed from 0 to 150 Hz.

Simultaneously, from the root locus with respect to rotor speed from 0 to 220 Hz, the MSR system loses its stability when the rotor speed is lower than 55 Hz or higher than 200 Hz, which also agrees well with the results from the proposed method.

To further testify the correctness of the stability criterion, the corresponding experiments have been performed. For the reason that the MSR is unstable at the speed of 0 Hz, the speed increases from 0 Hz to nearly 60 Hz by setting the cross proportion 48 Hz. And then the MSR decreases its speed slowly with 150 Hz cross proportion. It can be found that the MSR loses its stability and collides with the auxiliary bearing at about 52 Hz. Figures 8(a) and 8(b) present, respectively, its radial displacement responses in -direction (including channels and ) and the corresponding frequency spectrum characteristics before and after the whirling mode loses its stability. Besides, the Lissajous figure of terminal of the MSR, which is generated by and , is also tracked as shown in Figure 8(c). From Figure 8(a), it can be seen that the instability is caused by rotational mode instead of translational mode since the motion directions of channels and are opposite, which can also be testified from Figure 8(c). From Figure 8(b), the critical rotor speed is about 52 Hz, and its corresponding whirling frequency is about 70 Hz, which agrees well with the simulation results. Additionally, from Figure 8(c), it can be obtained that the instability mode is backward whirl instead of forward whirl since the whirling direction is opposite to the rotation.

(a)

(b)

(c)

Figure 8: Experimental results. (a) Radial displacement responses of channels and . (b) Frequency spectrum of displacement . (c) Lissajous figure of terminal of the MSR. Note: the arrowhead denotes the direction of the rotation and the circle drawn by dotted line means the auxiliary bearing.

4. Conclusion

To analyze directly and simply the stability of the cross-coupled two-dimensional antisymmetrical system, this paper presents a stability analysis method based on complex-coefficient frequency characteristics. Through variable reconstruction, the MIMO system is converted into a SISO system with complex coefficients. The closed-loop pole locations law and the controllability as well as observability before and after the variable reconstruction are studied and revealed, which demonstrate that the two systems before and after variable reconstruction possess the same absolute and relative stability. Then the Nyquist stability criterion is extended from a real-coefficient system to a complex-coefficient system according to argument principle.

Combined with the rigid MSR system with significant gyroscopic effects, simulation and experimental results demonstrate that the proposed stability analysis method is feasible, correct, and quite simple to implement. Compared with the existing methods, the presented one can greatly simplify the analysis, design, and debugging of the control system.

Appendix

In Table 1, is the radial clearance of the auxiliary bearings, that is, the protective air gaps of the ball bearings.

In Table 2, and are the proportional and feedback coefficients of the current-loop controller for the MSR system, , , are the proportional, integral, and differential coefficients of the basic PID controller for the MSR system; and are the cross coefficients of the high-pass and low-pass filters, respectively; is the total cross coefficient; and are the cut-off frequencies of the two-order high-pass and low-pass filters, respectively.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors are very grateful to Professor Y. Hong and Professor L. Guo from the Chinese Academy of Sciences and Professor S. S. Ge from the National University of Singapore for their valuable suggestions and instructive comments.

References

T. D. Nguyen, K.-J. Tseng, S. Zhang, and H. T. Nguyen, “A novel axial flux permanent-magnet machine for flywheel energy storage system: design and analysis,” IEEE Transactions on Industrial Electronics, vol. 58, no. 9, pp. 3784–3794, 2011.

S. P. Bhat and P. K. Tiwari, “Controllability of spacecraft attitude using control moment gyroscopes,” IEEE Transactions on Automatic Control, vol. 54, no. 3, pp. 585–590, 2009.

J. Fang and Y. Ren, “High-precision control for a single-gimbal magnetically suspended control moment gyro based on inverse system method,” IEEE Transactions on Industrial Electronics, vol. 58, no. 9, pp. 4331–4342, 2011.

J. Fang and Y. Ren, “Decoupling control of a magnetically suspended rotor in a control moment gyro based on inverse system method,” IEEE/ASME Transactions on Mechatronics, vol. 17, pp. 1133–1144, 2012.

H. Bang and H. D. Choi, “Attitude control of a bias momentum satellite using moment of inertia,” IEEE Transactions on Aerospace and Electronic Systems, vol. 38, no. 1, pp. 243–250, 2002.

Y. Li, Z. Q. Zhu, D. Howe, C. M. Bingham, and D. A. Stone, “Improved rotor-position estimation by signal injection in brushless AC motors, accounting for cross-coupling magnetic saturation,” IEEE Transactions on Industry Applications, vol. 45, no. 5, pp. 1843–1850, 2009.

S. S. Broussev and N. T. Tchamov, “Quality factor analysis for cross-coupled LC oscillators using a time-varying root locus,” IEEE Transactions on Circuits and Systems II, vol. 59, no. 1, pp. 25–29, 2012.

F.-J. Lin, P.-H. Chou, C.-S. Chen, and Y.-S. Lin, “DSP-based cross-coupled synchronous control for dual linear motors via intelligent complementary sliding mode control,” IEEE Transactions on Industrial Electronics, vol. 59, no. 2, pp. 1061–1073, 2012.

H.-J. Lee, S. Jung, and S.-K. Sul, “A current controller design for current source inverter-fed AC machine drive system,” IEEE Transactions on Power Electronics, vol. 28, pp. 1366–1381, 2013.

M. M. Hossain, G. Lu, and Y. Yan, “Optical fiber imaging based tomographic reconstruction of burner flames,” IEEE Transactions on Instrumentation and Measurement, vol. 61, no. 5, pp. 1417–1425, 2012.

F. Heismann, “Improved optical compensator for first- and second-order polarization-mode dispersion,” IEEE Photonics Technology Letters, vol. 17, no. 5, pp. 1016–1018, 2005.

H. Mukaidani and H. Xu, “Pareto optimal strategy for stochastic weakly coupled large scale systems with state dependent system noise,” IEEE Transactions on Automatic Control, vol. 54, no. 9, pp. 2244–2250, 2009.

J. Rootenberg and O. J. Oso, “On cross coupling and stability in nonlinear control systems,” IEEE Transactions on Automatic Control, vol. 16, no. 1, pp. 73–75, 1971.

D. Newman, “The analysis of cross-coupling effects on the stability of two-dimensional, orthogonal, feedback control systems,” IEEE Transactions on Automatic Control, vol. 5, pp. 314–320, 1960.

Z. Chen, “PI and sliding mode control of a Cuk converter,” IEEE Transactions on Power Electronics, vol. 27, no. 8, pp. 3695–3703, 2012.

M. Margaliot and G. Langholz, “The Routh-Hurwitz array and realization of characteristic polynomials,” IEEE Transactions on Automatic Control, vol. 45, no. 12, pp. 2424–2427, 2000.

A. A. El-Marhomy and N. E. Abdel-Sattar, “Stability analysis of rotor-bearing systems via Routh-Hurwitz criterion,” Applied Energy, vol. 77, no. 3, pp. 287–308, 2004.

M.-T. Ho, A. Datta, and S. P. Bhattacharyya, “An elementary derivation of the Routh-Hurwitz criterion,” IEEE Transactions on Automatic Control, vol. 43, no. 3, pp. 405–409, 1998.

C. Yang, Q. Zhang, J. Sun, and T. Chai, “Lur'e Lyapunov function and absolute stability criterion for Lur'e singularly perturbed systems,” IEEE Transactions on Automatic Control, vol. 56, no. 11, pp. 2666–2671, 2011.

K. S. Ray and D. D. Majumder, “Application of circle criteria for stability analysis of linear SISO and MIMO systems associated with fuzzy logic controller,” IEEE Transactions on Systems, Man and Cybernetics, vol. 14, no. 2, pp. 345–349, 1984.

M. Fardad and B. Bamieh, “An extension of the argument principle and nyquist criterion to a class of systems with unbounded,” IEEE Transactions on Automatic Control, vol. 53, no. 1, pp. 379–384, 2008.

R. Turner, S. Walton, and R. Duke, “A case study on the application of the Nyquist stability criterion as applied to interconnected loads and sources on grids,” IEEE Transactions on Industrial Electronics, vol. 60, no. 6, pp. 2740–2749, 2013.

A. G. Yepes, F. D. Freijedo, Ó. López, and J. Doval-Gandoy, “Analysis and design of resonant current controllers for voltage-source converters by means of nyquist diagrams and sensitivity function,” IEEE Transactions on Industrial Electronics, vol. 58, no. 11, pp. 5231–5250, 2011.

J. Kim, D. G. Bates, and I. Postlethwaite, “Evaluation of stochastic effects on biomolecular networks using the generalized nyquist stability criterion,” IEEE Transactions on Automatic Control, vol. 53, no. 8, pp. 1937–1941, 2008.

F. Fontana, “Use of the Nyquist stability criterion in the design of interactive audio digital filters,” IEEE Signal Processing Letters, vol. 18, no. 4, pp. 271–274, 2011.

M. Vidyasagar, R. K. Bertschmann, and C. S. Sallaberger, “Some simplifications of the graphical Nyquist criterion,” IEEE Transactions on Automatic Control, vol. 33, no. 3, pp. 301–305, 1988.

J. Fotakis, M. J. Grimble, and B. Kouvaritakis, “A comparison of characteristic locus and optimal designs for dynamic ship positioning systems,” IEEE Transactions on Automatic Control, vol. 27, no. 6, pp. 1143–1157, 1982.

K. S. Ray and B. A. Kouvaritakis, “An extension of the characteristic sequences method to the case of repeated roots,” IEEE Transactions on Automatic Control, vol. 29, no. 10, pp. 933–934, 1984.

M. N. Saadatzi, J. Poshtan, and M. H. Saadatzi, “Optimal multivariable two-degree-of-freedom control of electric wheelchair using non-dominated sorting genetic algorithm-II,” in Proceedings of the 20th IEEE International Conference on Control Applications (CCA '11), pp. 376–381, September 2011.

Yang, Q. Guan, P. X. Wang, and X. Zhang, “Stability analysis of matrix converter with constant power loads and LC input filter,” in Proceedings of the 7th International Power Electronics and Motion Control Conference (IPEMC '12), pp. 900–904, 2012.

R. S. Baheti, “Multivariable frequency domain controller for magnetic suspension and balance systems,” IEEE Transactions on Automatic Control, vol. 29, no. 8, pp. 725–728, 1984.

O. D. I. Nwokah, R. E. Nordgren, and G. S. Grewal, “Inverse Nyquist array: a quantitative theory,” IEE Proceedings: Control Theory and Applications, vol. 142, no. 1, pp. 23–30, 1995.

A. Postlethwaite and G. J. MacFarlane, A Complex Variable Approach to the Analysis of Linear Multivariable Feedback Systems, Springer, New York, NY, USA, 1979.

Y. Okada, B. Nagai, and T. Shimane, “Cross-feedback stabilization of the digitally controlled magnetic bearing,” Journal of Vibration and Acoustics, vol. 114, no. 1, pp. 54–59, 1992.

M. Ahrens, L. Kučera, and R. Larsonneur, “Performance of a magnetically suspended flywheel energy storage device,” IEEE Transactions on Control Systems Technology, vol. 4, no. 5, pp. 494–502, 1996.

G. V. Brown, A. Kascak, R. H. Jansen, T. P. Dever, and K. P. Duffy, “Stabilizing gyroscopic modes in magnetic-bearing-supported flywheels by using cross-axis proportional gains,” in AIAA Guidance, Navigation, and Control Conference, pp. 1132–1143, San Francisco, Calif, USA, August 2005, AIAA 2005-5955.

A. F. Kascak, G. V. Brown, R. H. Jansen, and T. P. Dever, “Stability limits of a PD controller for a flywheel supported on rigid rotor and magnetic bearings,” in AIAA Guidance, Navigation, and Control Conference, pp. 1144–1155, San Francisco, Calif, USA, August 2005, AIAA 2005-5956.

A. A. El-Marhomy, “Effect of various rotor-bearing system parameters on the whirl stability of system rotational modes of motion,” Energy Conversion and Management, vol. 38, no. 1, pp. 73–82, 1997.

F. Sorge and M. Cammalleri, “Control of hysteretic instability in rotating machinery by elastic suspension systems subject to dry and viscous friction,” Journal of Sound and Vibration, vol. 329, no. 10, pp. 1686–1701, 2010.

S. E. Mushi, Z. Lin, and P. E. Allaire, “Design, construction, and modeling of a flexible rotor active magnetic bearing test rig,” IEEE/ASME Transactions on Mechatronics, vol. 17, pp. 1170–1182, 2008.

Y. Ren, D. Su, and J. Fang, “Whirling modes stability criterion for a magnetically suspended flywheel rotor with significant gyroscopic effects and bending modes,” IEEE Transactions on Power Electronics, vol. 28, pp. 5890–5901, 2013.

Y. Ren and J. Fang, “Current-sensing resistor design to include current derivative in PWM H-bridge unipolar switching power amplifiers for magnetic bearings,” IEEE Transactions on Industrial Electronics, vol. 59, pp. 4590–4600, 2012.

K. S. Yeung and W. T. Wong, “Root-locus plot of systems with time delay,” Electronics Letters, vol. 18, no. 11, pp. 480–481, 1982.

J. Fang, Y. Ren, and Y. Fan, “Nutation and precession stability criterion of magnetically suspended rigid rotors with gyroscopic effects based on positive and negative frequency characteristics,” IEEE Transactions on Industrial Electronics, vol. 61, pp. 2003–2014, 2014.

U. T. Jonsson and C.-Y. Kao, “Consensus of heterogeneous LTI agents,” IEEE Transactions on Automatic Control, vol. 57, no. 8, pp. 2133–2138, 2012.