Generalizing Negative Imaginary Systems Theory to Include Free Body Dynamics:
Control of Highly Resonant Structures with Free Body Motion
Abstract
Negative imaginary (NI) systems play an important role in the robust control of highly resonant flexible structures. In this paper, a generalized NI system framework is presented. A new NI system definition is given, which allows for flexible structure systems with colocated force actuators and position sensors, and with free body motion. This definition extends the existing definitions of NI systems. Also, necessary and sufficient conditions are provided for the stability of positive feedback control systems where the plant is NI according to the new definition and the controller is strictly negative imaginary. The stability conditions in this paper are given purely in terms of properties of the plant and controller transfer function matrices, although the proofs rely on state space techniques. Furthermore, the stability conditions given are independent of the plant and controller system order. As an application of these results, a case study involving the control of a flexible robotic arm with a piezoelectric actuator and sensor is presented.
I Introduction
Flexible structure dynamics arise in many areas such as flexible robot manipulators [1], ground and aerospace vehicles [2], atomic force microscopes (AFMs) [3, 4] and other nanopositioning systems [5, 6, 7, 8]. Flexible structures can be modeled as infinite dimensional distributed parameter systems [9]. However, finite dimensional models are often used for the purpose of designing controllers [10, 11, 9, 12]. In designing controllers for these flexible systems, it is important to consider the effect of highly resonant modes. Such resonant modes are known to adversely affect the stability and performance of flexible structure feedback control systems [12, 13, 14], and are often very sensitive to changes in environmental variables. For instance, a small change in the environment of the system such as changing temperature, can lead to significant changes in the resonant frequencies of such systems. These changes in resonant frequencies can lead to large changes in the gain and phase of the system frequency response at a given frequency, which may lead to instability or poor performance in the corresponding feedback system. In addition, highly resonant modes lead to vibrational effects which limit the ability of control systems to achieve desired levels of performance in many applications such as precision instrumentation, optical systems, precision machine tools, wafer steppers, telescopes, and atomic force microscopes [12]. These issues arising from the presence of highly resonant modes in flexible structures motivate the need for tools to guarantee robust stability and performance in flexible structure control systems.
One common solution to issues of robustness, stability, and performance in the control of highly resonant flexible structures is to use force actuators combined with colocated measurements of velocity, position, or acceleration [12, 13, 14]. Colocated control with velocity measurements, known as negativevelocity feedback, can be used to directly increase the effective damping in the system, thereby facilitating the design of controllers that can guarantee closedloop stability in the presence of parameter variations and unmodeled plant dynamics [12]. Similarly, a class of colocated controllers with position measurements, known as positiveposition feedback controllers, where velocity sensors are replaced with position sensors, can also be used to increase damping in flexible systems as discussed in [13, 15]. Also, positiveposition feedback controllers are robust against uncertainties in resonant frequencies as well as unmodeled plant dynamics, in a similar way to negativevelocity feedback controllers [13, 16, 14].
The properties of negativevelocity feedback has been studied using passivity theory and the theory of positive real (PR) linear time invariant (LTI) systems; e.g., see [17, 18]. However, PR theory cannot be used directly when using position or acceleration measurements [14]. This drawback is important in applications to the field of nanotechnology, especially for nanopositioning systems, where position measurements are widely used; see e.g., [3, 19, 20, 5, 6, 21, 7, 22, 23, 8]. Similar issues also arise in application to the area of robotics where position measurements are also widely used.
Lanzon and Petersen introduced a notion of negative imaginary (NI) systems in [16, 14] for the robust control of flexible structures with force actuators combined with position or acceleration sensors. (SISO) case, NI systems are defined by considering the properties of the imaginary part of the system frequency response and requiring the condition for all . The NI property arises in many practical systems. For example, such systems arise when considering the transfer function from a force actuator to a corresponding colocated position sensor (for instance, a piezoelectric sensor) in a lightly damped structure [13, 14, 3, 24, 25]. Another area where the underlying system dynamics are NI, in the area of nanopositioning systems; see e.g., [3, 19, 20, 5, 6, 21, 7, 22, 23, 8]. Also, the positiveposition feedback control scheme in [13, 26], can be considered using the NI framework. Furthermore, other control methodologies in the literature such as integral resonant control (IRC) [27] and resonant feedback control [28, 29], fit into the NI framework and their stability robustness properties can be explained by NI systems theory.
The stability robustness of interconnected NI systems has been studied in [16, 14]. In these papers, it is shown that a necessary and sufficient condition for the internal stability of a positivefeedback control system (see Fig. 1) consisting of an NI plant with transfer function matrix and a strictly negative imaginary (SNI) controller with transfer function matrix is given by the DC gain condition
(1) 
where the notation denotes the maximum eigenvalue of a matrix with only real eigenvalues. This stability result has been used in a number of practical applications [25, 3, 4, 30, 31, 8]. For example in [25], this stability result is applied to the problem of decentralized control of large vehicle platoons. In [3, 4], the NI stability result is applied to nanopositioning in an atomic force microscope. A positive position feedback control scheme based on the NI stability result provided in [16, 14] is used to design a novel compensation method for a coupled fuselagerotor mode of a rotary wing unmanned aerial vehicle in [30]. In [8], an IRC scheme based on the stability results provided in [16, 14] is used to design an active vibration control system for the mitigation of human induced vibrations in lightweight civil engineering structures, such as floors and footbridges via proofmass actuators. An identification algorithm which enforces the NI constraint is proposed in [31] for estimating model parameters, following which an Integral resonant controller is designed for damping vibrations in flexible structures. In addition, it is shown in [32] that the class of linear systems having NI transfer function matrices is closely related to the class of linear Hamiltonian inputoutput systems. Also, an extension of the NI systems theory to infinitedimensional systems is presented in [33].
The NI framework presented in [16, 14] considers systems with poles in the open left half of the complex plane. This theory has been extended in [34] to include NI systems with poles in the closed left half of the complex plane, except at the origin. Also, further extensions to NI systems theory include the study of NI controller synthesis [35, 36], connections between NI systems analysis and analysis [37], and conditions for robust stability analysis of mixed NI and boundedreal classes of uncertainties [38]. Furthermore, the concept of lossless NI transfer functions is introduced in [39], an algebraic approach to the realization of a lossless NI behavior is presented in [40], and a spectral characterization of NI descriptor systems is discussed in [41]. The NI systems theory can be extended to nonlinear systems using the concept of counterclockwise inputoutput dynamics as presented in [42, 43, 44]. In [44], a sufficient conditions under which a semilinear Duhem model is counterclockwise is given, where the counterclockwise inputoutput system is restricted to periodic input signals. positive feedback interconnection for SISO linear case is provided in [44].
Despite generalizations of the NI systems framework presented in [34], an important class of systems, that cannot be captured by the existing NI systems framework, corresponds to flexible systems with free body motion. These systems arise in areas such as rotating flexible spacecraft [45], rotary cranes [46], robotics and flexible link manipulators [27, 29, 47], and dualstage hard disk drives [48, 49, 50, 51]. Flexible structures with free body motion lead to dynamical models including poles at the origin, which is not covered in earlier work on NI systems theory. In particular, the stability condition (1) is not well defined in the case of flexible structures with free body motion which results in poles at the origin, since in this case, the plant DC gain will be infinite. However, control systems involving flexible structures with free body motion arising in these important application areas still suffer from the stability and performance issues mentioned above. Thus we are motivated to extend the NI robust stability theory developed in [16, 14, 34] so that it can be applied to control systems involving highly resonant flexible structures with free body motion.
Fig. 2 shows a block diagram of a system which includes a flexible structure with free body motion that arises in a problem of diskdrive control; see [52]. Here, a voice coil motor (VCM) is used to actuate the arm of the reader head. The free body motion of the reader head leads to a transfer function from the input to the VCM output which has poles at the origin. Furthermore, the overall system satisfies the NI frequency response property and includes poles at the origin. However, the NI stability results presented in [16, 14, 34] do not allow for poles at the origin and cannot be applied to control systems such as this diskdrive control system.
In this paper, we present a new generalized definition of NI systems which allows for flexible structures with colocated force actuators and position sensors and with free body motion. This definition extends the previous definitions of NI systems presented in [16, 14, 34] to allow for up to two poles at the origin. We also derive new generalized stability conditions for positivefeedback control systems involving an NI plant and an SNI controller.
As in [16, 14, 34], the stability conditions presented in this paper, are given purely in terms of properties of the plant and controller transfer function matrices, although the proofs rely on state space techniques. Furthermore, the stability conditions given are independent of the plant and controller system order and can be stated without using the fact that the plant and the controller transfer function matrices are rational. However, the proofs given in this paper only apply to the rational case.
Preliminary conference versions of the stability results presented in this paper were presented in [53, 54]. However, in this paper, much more general versions of these stability results are presented in Theorems 1  4 and Corollaries 1, 2, which allow for the existence of free body motion in some but not all inputoutput channels. This is important since multivariable control systems involving flexible structures with free body motion usually include free body motion in some but not all inputoutput channels. Also, this paper includes a case study involving the control of a flexible robotic arm, which has not been considered in the previous conference versions of the paper.
This paper is further organized as follows: Section II recalls the existing definition for NI systems and outlines the notation that will be used in the rest of the paper. Section III introduces the new generalized definition for NI systems, which allows for systems with free body dynamics. Also in this section, we present the main stability results in Theorems 1  4 and Corollaries 12. Section IV presents a case study, which involves a flexible robotic arm, as an application of the NI theory presented in this paper. The paper is concluded with a summary and remarks on future work in Section V. All proofs of the presented theorems, lemmas and corollaries are given in the Appendix.
Ii Preliminaries and Notation
In this section, we recall the existing definitions of NI and SNI systems as given in [34] for systems with poles in the closed left half of the complex plane, except at the origin. We will use this existing definition of SNI (first introduced in [16]) systems but in the next section we will present our new definition of generalized NI systems. We also define notation used to describe positive feedback interconnections and internal stability, which will be used to present the main results in this paper.
Consider the following LTI system,
(2)  
(3) 
where and with the square transfer function matrix . The transfer function matrix is said to be strictly proper if . We will use the notation to denote the state space realization (2), (3).
The existing definition of NI systems states that a square transfer function matrix is NI if the following conditions are satisfied [34]:

has no pole at the origin and in .

The corresponding frequency response is such that
for all where is not a pole of .

If with is a pole of , it is at most a simple pole and the residue matrix is positive semidefinite Hermitian.
Definition 1
[34] A square transfer function matrix is SNI if the following conditions are satisfied:

has no pole in .

For all , .
Now, consider a positive feedback interconnection between an NI system with transfer function matrix and an SNI system with transfer function matrix as shown in Fig. 1. Also, suppose that the transfer function matrix has a minimal state space realization and has a minimal state space realization Furthermore, it is assumed that the matrix is nonsingular. Then the closed system has a system matrix given by
Iii Main results
The main contribution of this paper is a generalization of the framework for NI systems presented in [34]. We introduce a new definition of NI systems that will allow for systems with free body dynamics. This generalized definition will be used in a new set of stability conditions that will allow for NI systems with free body motion to be included into the framework of NI systems theory. Henceforth, when a system is said to be NI, we will mean NI as defined below, not NI as defined in earlier papers.
Definition 2
A square transfer function matrix is NI if the following conditions are satisfied:

has no pole in .

For all such that is not a pole of ,

If with is a pole of , then it is a simple pole and the residue matrix is Hermitian and positive semidefinite.

If is a pole of , then for all and is Hermitian and positive semidefinite.
Here, is the frequency response corresponding to the transfer function . Unlike the NI definition presented in [34], Definition 2 allows for poles at the origin. In this case, we cannot use the existing stability results presented in [14, 34, 35], because the stability condition in (1) is not defined. The inclusion of poles at the origin extends the NI systems theory to include flexible systems with free body dynamics. In order to derive a set of stability conditions that allow for NI systems with free body motion, we define the following constant matrices for a given NI transfer function matrix
(5) 
These matrices are the first three coefficients in the Laurent series expansion of the transfer function . These matrices carry information about properties of the free body motion of the system under consideration and will be used in stability conditions for the positive feedback interconnection of NI and SNI systems. Note that the DC gain condition (1) cannot be defined for an NI system with transfer function matrix unless , which reduces to the case where the dynamical system has no free body motion. From Condition 4) in Definition 2, the matrix is required to be Hermitian and positive semidefinite. Hence, it follows (e.g., see [56]) that if , it can be decomposed in the form
(6) 
where is a full column rank matrix.
We now present conditions for the stability of a positive feedback control system involving an NI plant with free body motion. These conditions are stated using the quantities defined in (5). First, we define the Hankel matrix as
(7) 
Suppose that . Using the singular value decomposition (SVD), we can decompose the Hankel matrix as
(8) 
where are unitary matrices, , , and the matrices and each have orthogonal columns. Furthermore, we can decompose the matrix using the SVD as
(9) 
where and are orthogonal matrices, and .
We now introduce some notation which will be used throughout the paper. Given matrices and such that , then the matrix valued function is defined by
(10) 
Using this notation, we define the matrix
(11) 
where the matrix is given by
(12) 
and we will assume that .
We will use the following condition in the theorem which follows:
(13) 
Also, for the case in which is positive semidefinite, we will use the condition
(14) 
Moreover, for the case in which is negative semidefinite, we will use the condition
(15) 
Here, and matrices and are defined in (5), (6), (11), and (12) respectively. Also, denotes the square root of a positive semidefinite matrix.
The following theorem is our first main stability result for the case in which . That is, the system has double poles at the origin.
Theorem 1
Suppose that the square transfer function matrix is strictly proper and NI with , and the transfer function matrix is SNI. Also, suppose that the matrix is nonsingular. If is positive semidefinite, then the closedloop positivefeedback interconnection between and as shown in Fig. 1 is internally stable if and only if conditions (13) and (14) are satisfied. Furthermore, if is negative semidefinite, then the closedloop positivefeedback interconnection between and is internally stable if and only if conditions (13) and (15) are satisfied.
The proof of this and subsequent theorems and corollaries are presented in Appendix B.
We now present a corollary to this theorem which considers the special case in which none of the free body modes of the plant have frictional force present; i.e., . In order to present this corollary, we define the matrix as follows:
(16) 
where we assume that the matrix is nonsingular.
We will use the following condition in the next corollary, which corresponds to condition (13) in Theorem 1:
(17) 
Also, for the case in which is positive semidefinite, we will use the following condition which corresponds to condition (14) in Theorem 1:
(18) 
Moreover, for the case in which is negative semidefinite, we will use the following condition which corresponds to condition (15) in Theorem 1:
(19) 
where .
Corollary 1
Suppose that the transfer function matrix is SNI and the strictly proper transfer function matrix is NI with and . Also, suppose that the matrix is nonsingular. If is positive semidefinite, then the closedloop positivefeedback interconnection between and is internally stable if and only if conditions (17) and (18) are satisfied. Furthermore, if is negative semidefinite, then the closedloop positivefeedback interconnection between and is internally stable if and only if conditions (17) and (19) are satisfied.
The following theorem imposes some extra conditions on the matrix which enables us to relax the sign definiteness condition on the matrix . This then leads to a simplified stability condition.
Theorem 2
Suppose that the transfer function matrix is SNI and the strictly proper transfer function matrix is NI with and . Also, suppose that , where denotes the null space of a matrix. Then the closedloop positivefeedback interconnection between and is internally stable if and only if condition (17) is satisfied.
In Theorem 3, Theorem 4 and Corollary 2, we consider cases which correspond to free body motion with frictional force present. As in Theorem 1, these cases allow for fact that the free body motion may not be present in all inputoutput channels.
In order to present Theorem 3 and Theorem 4, suppose that and . This corresponds to the case when the system has a single pole at the origin. Then we consider the following SVD decomposition of the matrix defined in (5):
(20) 
where , and the matrices and each have orthogonal columns. Also, we define the matrix as follows:
(21) 
where the matrix is assumed to be nonsingular.
We will use the following condition in Theorem 3 and Corollary 2 which corresponds to condition (13) in Theorem 1:
(22) 
For the case in which is positive semidefinite, we also will use the following condition which corresponds to condition (14) in Theorem 1:
(23) 
Moreover, for the case in which is negative semidefinite, we will use the following condition which corresponds to condition (15) in Theorem 1:
(24) 
where .
Theorem 3
Suppose that the transfer function matrix is SNI and the strictly proper transfer function matrix is NI with and . Also, suppose that the matrix nonsingular. If is positive semidefinite, then the closedloop positivefeedback interconnection between and is internally stable if and only if conditions (22) and (23) are satisfied. Furthermore, if is negative semidefinite, then the closedloop positivefeedback interconnection between and is internally stable if and only if conditions (22) and (24) are satisfied.
The following theorem imposes some extra conditions on the matrix which enables us to relax the sign definiteness condition on the matrix . This then leads to a simplified stability condition.
Theorem 4
Suppose that the transfer function matrix is SNI and the strictly proper transfer function matrix is NI with and . Also, suppose that . Then the closedloop positivefeedback interconnection between and is internally stable if and only if condition (22) is satisfied.
Corollary 2
Suppose that the transfer function matrix is SNI and the strictly proper transfer function matrix is NI with either and invertible or and . Then, the closedloop positivefeedback interconnection between and is internally stable if and only if
Remark 1
The case where and corresponds to the existing stability results presented in [14, 34, 35]. In this case, the stability condition reduces to This condition can be obtained from (23) using the fact in this case. Also, we require the assumption . Hence,
Note that using a similar argument to the proof of Theorem 3, we can obtain a similar result under the assumption that
Iv Case Study: Control of Flexible robotic arm
In this section, we present an application of the stability results presented in this paper to the control of a flexible robotic arm system. The robotic arm is pinned to a motor at one end. For the purposes of modeling the flexible robotic arm, we use an equivalent slewing beam model as depicted in Fig. 3; see [57].
The motor allows the robotic arm to traverse in the vertical plane. Two piezoelectric patches are attached to the arm on either side. Here, one piezoelectric patch acts as an actuator while the other is a sensor. The robotic arm system has two inputs and two outputs: the inputs are the voltage applied to the piezoelectric actuator and the torque applied by the motor, whereas the outputs are the voltage produced by the piezoelectric sensor and the motor hub angle . The fact that this system involves colocated “force” actuators and “position” sensors indicates that the system will be NI; e.g., see [14].
Iva Mathematical model for the robotic arm
The beam in Fig. 3 is modeled using the BernoulliEuler equations of motion for a beam with actuating and sensing piezoelectric elements as in [57]:
(25) 
Here, is Young’s modulus and is the second moment of inertia of the beam, is the density of the beam, and is the area of the composite beam. If the thickness of the piezoelectric films are comparable to the thickness of the beam, then the products and would be different in the laminated and nonlaminated areas of the beam. However, since piezoelectric films used in practical applications are often thin compared to the thickness of the beam, these differences will be neglected. Assuming that the products and are uniform over the length of the beam simplifies the modeling procedure.
Now we consider various boundary conditions in modeling the beam. These are given as
(26)  
(27)  
(28)  
(29) 
Here, and are the mass and inertia of the tip, which will be neglected in this paper.Also, (26) represents the inability of the motor joint to undergo transverse motion. As in [57], the time domain beam equation (25) with boundary conditions (26)(29) can be transformed into an equivalent Laplace domain representation as
(30) 
with boundary conditions
(31)  
(32)  
(33)  
(34) 
where the primes indicate spatial derivatives and
(35) 
Note that (30) is the Laplace domain equivalent of the BernoulliEuler beam equation with as a forcing input. Together, (30)(34) represent a set of linear ordinary differential equations with mixed boundary conditions: two at and two at . A state space representation for the system can be formed from equations (30)(34) as in [57]:
(36) 
where represents the Dirac delta function. The equation (IVA) can be written in the general form
(37) 
the solution to which is given by
(38) 
Once the boundary conditions and are known, (IVA) will depend upon three conditions for , namely, . For further details see [57].
IvB Infinite Dimensional Transfer function Model
Here, we present the inputoutput relationship between the two inputs and , and the corresponding collocated outputs and in the form of the transfer function matrix,
(39) 
where and each of the elements of this transfer function matrix is an infinite dimensional transfer function defined in terms of transcendental functions of . Indeed, each of the four transfer functions in (39) can be written as a ratio of numerator and denominator functions computed as
(40)  
(41)  
(42)  
(43) 
Here,
(44) 
where is the hub inertia. Also, the functions , are given by very complicated expressions which can be found in equations (26)(28) in [57].
We now compute the transfer functions in (40)(43) for the case where the piezoelectric actuators and sensors span the entire length of the beam. This corresponds to the substitutions: and . The resulting transfer functions have been verified in [58] for an experimented robotic arm system.
Despite the fact that we have not defined the NI property for infinite dimensional transfer functions, we will provide some calculations which indicate that the infinite dimensional transfer function matrix defined in (39)(43) satisfies the NI conditions given in Definition 2. Since the infinite dimensional transfer function matrix is actually a transcendental function of , , then Condition 2) in Definition 2 is equivalent to the condition
(45) 
for all where is given by (35). Indeed, it is straightforward to verify from the formulas for the transfer function matrix (39)(43) that for all . Also, the function given in (IVB) has an infinite numbers of roots. However, we can check Condition 3) in Definition 2 for a finite number of these roots on the imaginary axis. To do so, we have calculated the first eleven axis roots of numerically. At each of these roots , the corresponding residue matrix is calculated using L’Hopital’s rule as follows:
(46) 
where denotes the first derivative of with respect to .
In this case study, the parameter values for the robotic arm are taken from [58]. These parameter values are shown in the Table I.
Parameter  Value  Unit 

Hub inertia,  0.0348  
Beam length,  2  
Volumetric mass density,  2712.6  
Cross sectional area,  483.87  
Young’s Modulus,  69.0  
Area moment of inertia,  1.63  
Coupling Coefficient  0.340   
Capacitance,  68.35  
Thickness  3.05 
.
Table II shows the calculated roots of and the minimum eigenvalue of the corresponding residue matrix given in (IVB). Also, the matrix is found to be which is positive semidefinite.
n  Root (rad/s)  Minimum eigenvalue 

0  0  0 
1  3.395326441  0.1434 
2  9.501801884  0.2553 
3  17.08210071  0.1320 
4  29.32863976  0.0361 
5  47.01240951  0.0142 
6  96.84550724  0.0049 
7  128.7332003  0.0034 
8  165.2195349  0.0025 
9  206.2898971  0.0019 
10  251.9420283  0.0015 
.
These results show that the infinite dimensional transfer function matrix satisfies the conditions of Definition 2, at least for the first ten resonant modes.
IvC Approximate Finitedimensional Transfer Function Matrix
The transfer functions in (40)(43) are irrational functions of . We now approximate these transfer functions by rational functions in in order to design a suitable controller for the robotic arm system and to simulate its performance. Various methods such as the Maclaurin series expansion presented in [59], the RayleighRitz method [60], and the assumed modes method [60] are available in literature for the finite dimensional approximation of such an infinite dimensional model. Here, we adopt a partial fraction approach to obtain a finite dimensional approximation of . This method is similar to the assumed modes technique described in [60]. The finite dimensional model can be written as
(47) 
Here in the infinite dimensional model is approximated by
(48) 
where, are the first n axis roots of . Also, the coefficient matrices are computed using a partial fraction expansion method. That is,
(49) 
The constant is chosen so that
(50) 
where is such that is not a root of . We consider the first resonant mode; i.e., for the controller design. The corresponding coefficient matrices were computed and were found to be and . Also, the poles were computed to be .
The finite dimensional model in (IVC) is NI, since for all , where is not a pole for . This follows because in this example, is real and symmetric for all such that is not a pole of . Also, the coefficient matrices are positive semidefinite which implies that Condition 3) in Definition 2 is satisfied. Moreover, which implies that Condition 4) in Definition 2 is satisfied.
IvD Controller design
According to Theorem 1 if a plant is NI, any SNI controller which satisfies the conditions of Theorem 1 will stabilize the system. The fact that the robotic arm plant involves colocated “force” actuators and “position” sensors indicates that this plant should be NI. In particular, the finite dimensional approximation to the robotic arm model derived in Subsection IVC was shown to be NI. We will now use a finite dimensional model of the form (IVC) to design a controller for the system. First, we compute the matrices , and in (5), for the finite dimension approximate system where n=1 in (IVC) to obtain
(51) 
This implies that we can use Corollary 1 to guarantee the stability of the positive feedback interconnection between the plant and an SNI controller.
In this case study, an integral resonant controller (IRC) is chosen to stabilize the system; e.g., see [14]. An IRC is a first order controller which takes the form
(52) 
This controller is SNI if and is a symmetric matrix [14]. Now, we chose the controller matrices and such that the conditions of Corollary 1 are satisfied. We choose the controller matrices as follows:
(53) 
This leads to a controller DC gain matrix of . To check the stability conditions in Corollary 1, we first compute the matrix in (6) using in (IVD). This yields Also, the matrix in (16) is calculated as , which is negative semidefinite. Then we conclude where Also, Thus, the conditions of Corollary 1 are satisfied.
To verify the performance of the closed loop system, we simulate the response of this system corresponding to a step change in the reference position of the robotic arm; see Fig. 5. This step response is shown in Fig. 6. Also, the corresponding response of the piezo sensor output is shown in Fig. 7. Here, the step responses were calculated using finite dimensional plant models defined in (IVC) for different numbers of modes, n=2,3…7.
To this end, we have used the proposed controller which is designed for the finite dimensional model with n=1 when applied to the plant with finite dimensional model where n=2,3…7 in order to check the performance and robustness of the proposed controller. In fact, the performance of the closed loop system is found to improve by increasing the number of modes; see Fig 6 and Fig 7.
Note that the controller parameters in (53) were chosen by process of trial and error to obtain good closed loop performance the case of the nominal plant model, . An alternative approach, which would be useful in the case of a more complicated SNI controller structure, would be to use an optimization procedure to obtain the controller parameters; e.g., see [61].
V Conclusion
In this paper, new stability results for the positivefeedback interconnection of negative imaginary systems have been derived. A new NI definition is presented, which allows for systems having free body dynamics to be considered as NI systems. This work can be used in controller design to allow for a broader class of NI systems than considered in previous work. The application of the main results in this paper has been illustrated via a case study involving the control of a flexible robotic arm.
References
 [1] D. G. Wilson, R. D. Robinett, G. G. Parker, and G. P. Starr, “Augmented sliding mode control for flexible link manipulators,” Journal of Intelligent and Robotic Systems, vol. 34, no. 4, pp. 415–430, 2002.
 [2] M. Harigae, I. Yamaguchi, T. Kasai, H. Igawa, and T. Suzuki, “Control of large space structures using GPS modal parameter identification and attitude and deformation estimation,” Electronics and Communications in Japan, vol. 86, no. 4, pp. 63–71, 2003.
 [3] B. Bhikkaji and S. Moheimani, “Fast scanning using piezoelectric tube nanopositioners: A negative imaginary approach,” in Proc. IEEE/ASME Int. Conf. Advanced Intelligent Mechatronics AIM, Singapore, July 2009, pp. 274–279.
 [4] I. A. Mahmood, S. O. R. Moheimani, and B. Bhikkaji, “A new scanning method for fast atomic force microscopy,” IEEE Transactions on Nanotechnology, vol. 10, no. 2, pp. 203–216, 2011.
 [5] S. Salapaka, A. Sebastian, J. P. Cleveland, and M. V. Salapaka, “High bandwidth nanopositioner: A robust control approach,” Review of Scientific Instruments, vol. 73, no. 9, pp. 3232–3241, 2002.
 [6] J. R. van Hulzen, G. Schitter, P. M. J. Van den Hof, and J. van Eijk, “Modal actuation for high bandwidth nanopositioning,” in Proc. American Control Conference, Baltimore, Maryland, USA, July 2010, pp. 6525–6530.
 [7] S. Devasia, E. Eleftheriou, and S. O. R. Moheimani, “A survey of control issues in nanopositioning,” IEEE Transactions on Control Systems Technology, vol. 15, no. 5, pp. 802–823, 2007.
 [8] I. M. Diaz, E. Pereira, and P. Reynolds, “Integral resonant control scheme for cancelling humaninduced vibrations in lightweight pedestrian structures,” Structural Control and Health Monitoring, vol. 19, no. 1, pp. 55–69, 2012.
 [9] W. Ray, “Some recent applications of distributed parameter systems theory âa survey,” Automatica, vol. 14, no. 3, pp. 281 – 287, 1978.
 [10] R. Curtain and K. Morris, “Transfer functions of distributed parameter systems: A tutorial,” Automatica, vol. 45, no. 5, pp. 1101 – 1116, 2009.
 [11] M. Demetriou and R. Smith, Research Directions in Distributed Parameter Systems, ser. Frontiers in applied mathematics. Society for Industrial and Applied Mathematics (SIAM, 3600 Market Street, Floor 6, Philadelphia, PA 19104), 2003.
 [12] A. Preumont, Vibration Control of Active Structures: An Introduction, N. Jacobs, Ed. Springer, 2011.
 [13] J. L. Fanson and T. K. Caughley, “Positive position feedback control for large space structures,” AIAA Journal, vol. 28, no. 4, pp. 717–724, Apr. 1990.
 [14] I. R. Petersen and A. Lanzon, “Feedback control of negative imaginary systems,” IEEE Control System Magazine, vol. 30, no. 5, pp. 54–72, 2010.
 [15] M. J. Balas, “Direct velocity feedback control of large space structures1,” AIAA Journal of Guidance and Control, vol. 2, no. 3, pp. 252–253, MayJune 1979.
 [16] A. Lanzon and I. R. Petersen, “Stability robustness of a feedback interconnection of systems with negative imaginary frequency response,” IEEE Transactions on Automatic Control, vol. 53, no. 4, pp. 1042–1046, 2008.
 [17] B. D. O. Anderson and S. Vongpanitlerd, Network Analysis and Synthesis: A Modern Systems Approach. Englewood Cliffs, N.J., USA: PrenticeHall, 1973.
 [18] B. Brogliato, R. Lozano, B. Maschke, and O. Egeland, Dissipative Systems Analysis and Control, 2nd ed., ser. Communications and Control Engineering. London, UK: Springer, 2007.
 [19] J. Dong, S. M. Salapaka, and P. M. Ferreira, “Robust MIMO control of a parallel kinematics nanopositioner for high resolution high bandwidth tracking and repetitive tasks,” in Proc. 46th IEEE Conf. Decision and Control, New Orleans, LA, Dec 2007, pp. 4495–4500.
 [20] A. Sebastian and S. M. Salapaka, “Design methodologies for robust nanopositioning,” IEEE Transactions on Control Systems Technology, vol. 13, no. 6, pp. 868–876, 2005.
 [21] Y. Michellod, P. Mullhaupt, and D. Gillet, “Strategy for the control of a dualstage nanopositioning system with a single metrology,” in Proc. IEEE Conf. Robotics, Automation and Mechatronics, Bangkok, June 2006, pp. 1–8.
 [22] R. K. Messenger, Q. T. Aten, T. W. McLain, and L. L. Howell, “Piezoresistive feedback control of a MEMS thermal actuator,” Journal of Microelectromechanical Systems, vol. 18, no. 6, pp. 1267–1278, 2009.
 [23] K. El Rifai, O. El Rifai, and K. YoucefToumi, “On dual actuation in atomic force microscopes,” in Proc. American Control Conference, vol. 4, 2004, pp. 3128–3133.
 [24] Y. K. Yong, B. Ahmed, and S. O. R. Moheimani, “Atomic force microscopy with a 12electrode piezoelectric tube scanner,” Review of Scientific Instruments, vol. 81, no. 3, p. 033701, March 2010.
 [25] C. Cai and G. Hagen, “Stability analysis for a string of coupled stable subsystems with negative imaginary frequency response,” IEEE Transactions on Automatic Control, vol. 55, no. 8, pp. 1958–1963, Aug. 2010.
 [26] C. J. Goh and T. K. Caughley, “On the stability problem caused by finite actuator dynamics in the collocated control of large space structures,” International Journal of Control, vol. 41, no. 3, pp. 787–802, Mar. 1995.
 [27] E. Pereira, S. S. Aphale, V. Feliu, and S. O. R. Moheimani, “Integral resonant control for vibration damping and precise tippositioning of a singlelink flexible manipulator,” IEEE/ASME Transactions on Mechatronics, vol. 16, no. 2, pp. 232–240, 2011.
 [28] D. Halim and S. O. R. Moheimani, “Spatial resonant control of flexible structuresapplication to a piezoelectric laminate beam,” IEEE Transactions on Control Systems Technology, vol. 9, no. 1, pp. 37–53, 2001.
 [29] I. A. Mahmood, S. O. R. Moheimani, and B. Bhikkaji, “Precise tip positioning of a flexible manipulator using resonant control,” IEEE /ASME Transactions on Mechatronics, vol. 13, no. 2, pp. 180–186, 2008.
 [30] B. Ahmed and H. Pota, “Dynamic compensation for control of a rotary wing UAV using positive position feedback,” Journal of Intelligent and Robotic Systems, vol. 61, no. 14, pp. 43–56, 2011.
 [31] B. Bhikkaji, S. O. R. Moheimani, and I. R. Petersen, “A negative imaginary approach to modeling and control of a collocated structure,” IEEE/ASME Transactions on Mechatronics, vol. 17, no. 4, pp. 717–727, 2012.
 [32] A. van der Schaft, “Positive feedback interconnection of Hamiltonian systems,” in Proceedings of the 50th IEEE Conference on Decision and Control and European Control Conference (CDCECC), Orlando, FL, USA, Dec 2011.
 [33] M. R. Opmeer, “Infinitedimensional negative imaginary systems,” IEEE Transactions on Automatic Control, vol. 56, no. 12, pp. 2973–2976, 2011.
 [34] J. Xiong, I. R. Petersen, and A. Lanzon, “A negative imaginary lemma and the stability of interconnections of linear negative imaginary systems,” IEEE Transactions on Automatic Control, vol. 55, no. 10, pp. 2342–2347, 2010.
 [35] Z. Song, A. Lanzon, S. Patra, and I. Petersen, “Towards controller synthesis for systems with negative imaginary frequency response,” IEEE Transactions on Automatic Control, vol. 55, no. 6, pp. 1506–1511, June 2010.
 [36] Z. Song, A. Lanzon, S. Patra, and I. R. Petersen, “A negativeimaginary lemma without minimality assumptions and robust statefeedback synthesis for uncertain negativeimaginary systems,” Systems & Control Letters, vol. 61, no. 12, pp. 1269 – 1276, 2012. [Online]. Available: http://www.sciencedirect.com/science/article/pii/S0167691112001557
 [37] S. Engelken, S. Patra, A. Lanzon, and I. R. Petersen, “Stability analysis of negative imaginary systems with real parametric uncertainty  the SISO case,” IET Control Theory and Applications, vol. 4, no. 11, pp. 2631 – 2638, 2010.
 [38] S. Patra and A. Lanzon, “Stability analysis of interconnected systems with mixed negativeimaginary and smallgain properties,” IEEE Transactions on Automatic Control, vol. 56, no. 6, pp. 1395 – 1400, 2011.
 [39] J. Xiong, I. R. Petersen, and A. Lanzon, “Finite frequency negative imaginary systems,,” IEEE Transactions on Automatic Control, vol. 57, pp. 2917–2922, 2012.
 [40] S. Rao, “An algebraic approach to the realization of lossless negative imaginary behaviors,” SIAM Journal on Control and Optimization, vol. 50, no. 3, pp. 1700–1720, 2012.
 [41] P. Benner and M. Voigt, “Spectral characterization and enforcement of negative imaginariness for descriptor systems,” April 2012, not. [Online]. Available: {}{}}{http://www.mpimagdeburg.mpg.de/preprints/}{cmtt}
 [42] D.~Angeli, ``Systems with counterclockwise inputoutput dynamics,'' IEEE Transactions on Automatic Control, vol.~51, no.~7, pp. 11301143, 2006.
 [43] , ``Multistability in systems with counterclockwise input output dynamics,'' IEEE Transactions on Automatic Control, vol.~52, no.~4, pp. 596 609, 2007.
 [44] A.~K. Padthe, J.~H. Oh, and D.~S. Bernstein, ``Counterclockwise dynamics of a rateindependent semilinear Duhem model,'' in Proceedings of the 44th IEEE Conference on Decision and Control, Seville, Spain, December 2005, pp. 80008005.
 [45] P.~Hughes, Spacecraft Attitude Dynamics. New York: Wiley, 1986.
 [46] T.~Gustafsson, ``On the design and implementation of a rotary crane controller,'' European Journal of Control, vol.~2, no.~2, pp. 166175, 1996.
 [47] S.~B. Choi, S.~S. Cho, H.~C. Shin, and H.~K. Kim, ``Quantitative feedback theory control of a singlelink flexible manipulator featuring piezoelectric actuator and sensor,'' Smart Materials and Structures, vol.~8, no.~3, pp. 338349, 1999.
 [48] Y.~Li, R.~Horowitz, and R.~Evans, ``Vibration control of a PZT actuated suspension dualstage servo system using a PZT sensor,'' IEEE Transactions on Magnetics, vol.~39, no.~2, pp. 932  937, 2003.
 [49] S.M. Suh, C.~C. Chung, and S.H. Lee, ``Discretetime LQG/LTR dualstage controller design in magnetic disk drives,'' IEEE Transactions on Magnetics, vol.~37, no.~4, pp. 1891  1895, 2001.
 [50] T.~B. Goh, Z.~Li, B.~M. Chen, T.~H. Lee, and T.~Huang, ``Design and implementation of a hard disk drive servo system using robust and perfect tracking approach,'' IEEE Transactions on Control Systems Technology, vol.~9, no.~2, pp. 221233, 2001.
 [51] S.~Devasia, ``Timeoptimal control with pre/post actuation for dualstage systems,'' IEEE Transactions on Control Systems Technology, vol.~20, no.~2, pp. 323334, 2012.
 [52] C.~Laorpacharapan and L.~Y. Pao, ``Shaped timeoptimal feedback control for diskdrive systems with backelectromotive force,'' IEEE Transactions on Magnetics, vol.~40, no.~1, pp. 8596, 2004.
 [53] M.~A. Mabrok, A.~G. Kallapur, I.~R. Petersen, and A.~Lanzon, ``A stability result on the feedback interconnection of negative imaginary systems with poles at the origins,'' in Proceedings of the 2012 Australian Control Conference, Sydney, Australia, 2012.
 [54] , ``A new stability result for the feedback interconnection of negative imaginary systems with a