Journal of Applied Mathematics

Volume 2014 (2014), Article ID 703178, 9 pages

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

Research Article

Solutions of a Quadratic Inverse Eigenvalue Problem for Damped Gyroscopic Second-Order Systems

Hong-Xiu Zhong, Guo-Liang Chen, and Xiang-Yun Zhang

Department of Mathematics, East China Normal University, Shanghai 200241, China

Received 2 September 2013; Accepted 15 December 2013; Published 21 January 2014

Academic Editor: Zhongxiao Jia

Copyright © 2014 Hong-Xiu Zhong et al. 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.

Given pairs of complex numbers and vectors (closed under conjugation), we consider the inverse quadratic eigenvalue problem of constructing real matrices , , , and , where , and are symmetric, and is skew-symmetric, so that the quadratic pencil has the given pairs as eigenpairs. First, we construct a general solution to this problem with . Then, with the special properties and , we construct a particular solution. Numerical results illustrate these solutions.

1. Introduction

Vibrating structures such as buildings, bridges, highways, and airplanes are distributed parameter systems [1]. Very often a distributed parameter system is first discretized to a matrix second-order using techniques of finite element or finite difference, and then an approximate solution is obtained for the discretized model. Associated with the matrix second-order model is the eigenvalue problem of the quadratic pencil,

where , , , and are, respectively, mass, damping, gyroscopic and stiffness matrices.

The system represented by (1) is called damped gyroscopic system. In general, the gyroscopic matrix is always skew-symmetric, the damping matrix and the stiffness matrix are symmetric, the mass matrix is symmetric positive definite, and they are all real matrices. If , the system is called damped nongyroscopic system, and if , the system is called undamped gyroscopic system.

The damped gyroscopic system has been widely studied in two aspects: the quadratic eigenvalue problem (QEP) and the quadratic inverse eigenvalue problem (QIEP). The QEP involves finding scalars and nonzero vectors , called the eigenvalues and eigenvectors of the system, to satisfy the algebraic equation , when the coefficient matrices are given. Many authors have been devoted to this kind of problems and a series of good results have been made (see, e.g., [2–8]). The QIEP determines or estimates the parameters of the system from observed or expected eigeninformation of . Our main interest in this paper is the corresponding inverse problem: given partially measured information about eigenvalues and eigenvectors, we reconstruct matrices , , , and , satisfied with several conditions, so that has the given eigenpairs. The problem we considered is stated as follows.

Problem 1. Given an eigeninformation pair , where

with

find real matrices , , , and , with being symmetric definite, and being symmetric, and being skew-symmetric, so that

In [9], Gohberg et al. developed a powerful GLR theory to solve the QIEP of the undamped gyroscopic system. In [10], Chu and Xu developed an elegant procedure to obtain a real-valued spectral decomposition of the damped nongyroscopic system. And then, Jia and Wei [11] derived a real-valued spectral decomposition of the undamped gyroscopic system. However, [10, 11] both need all the eigen-information of to obtain the parameters, and it is often impractical or impossible to obtain complete spectral information. Thus, it becomes very interesting to consider a QIEP with only a subset of eigenpairs known.

In [12], Kuo et al. constructed the solutions of the QIEP of the damped nongyroscopic system with given eigenpairs. And for the same system, Cai et al. [13] solved the QIEP with given eigenpairs. Meanwhile, for the damped gyroscopic system, Yuan [14] solved the QIEP with given eigenpairs. In [14], Yuan constructed symmetric positive semidefinite matrix and skew-symmetric matrix for , with and as given matrices. So, it becomes challenging to construct for damped gyroscopic system (1) with given eigenpairs, and this is the goal of this paper.

This paper is organized as follows. In Section 2, we establish the solubility theory of the Problem 1. In Section 3, we develop a simple method to compute a particular solution to the Problem 1 with . Moreover, for and , a simple algorithm is developed to compute a solution in Section 4. Some numerical results are presented in Section 5 to illustrate our main results. In the last section, some conclusions and acknowledgements are given.

Throughout this paper, we use capital letters to denote matrices, and lowercase (bold) letters to denote scalars (vectors). denotes the transpose of the matrix , denotes the identity matrix, and denotes the Moore-Penrose generalized inverse of . We write if is real symmetric positive (semidefinite). The spectrum of is denoted by .

For simplicity, we make the following assumptions.

(A1)The eigenvector matrix in Problem 1 has full column rank, that is, rank .(A2)The eigenvalue matrix in Problem 1 has only simple eigenvalues, that is, .

Remark 2. For the case that , using the assumption of simple eigenvalues, we can partition , where has no zero eigenvalue, and then do discussion with . So, in this paper, we only consider the case that has no zero eigenvalue.

2. General Solution of the Problem

In this section, we will give a general solution to the Problem 1 for given matrix pair as in (2) and (3). At the beginning, we will introduce some lemmas.

Lemma 3 (see [15]). Let and ; then

has a solution if and only if

where is the Moore-Penrose generalized inverse of .

When condition (6) is satisfied, the general solution of (5) is

where is arbitrary and is constrained only by the symmetry requirement that

Lemma 3 directly results in the following lemma.

Lemma 4. Let be a nonsingular matrix and ; then

has a solution if and only if

in which case the general solution is , where is an arbitrary matrix.

Given matrix pair as in Problem 1: let

be the -factorization of , where is orthogonal and is upper triangular. We may require that has positive diagonal entries, since is of full column rank.

Let , so that finding , , , and which satisfy (4) is equivalent to finding , , and which satisfy

and the relations of , , and are

Let ; we can see exists by using . Denoting

where , , , , , , , , , , we will obtain the following main theorem.

Theorem 5. Let , , and be defined as in (14)–(16); then there are real matrices , and satisfy (4) if and only if

(i) is arbitrarily symmetric positive definite,(ii), , are arbitrary,(iii) is arbitrary symmetric,(iv), where is arbitrary symmetric,(v),(vi).

Furthermore, and can be expressed as in (13).

Proof. Necessity. From (14)–(16), we know , , and ; substituting them and (11) into (12), we have

Thus, finding , , and which satisfy (12) is equivalent to finding the submatrices , , , , , and which satisfy (17) and (18). Clearly, it follows from (18) that is determined by

where and are arbitrary.

As and are required to be symmetric positive definite and symmetric, respectively, so are and in (14) and (16). From (17) it follows that

Let be an arbitrary symmetric positive definite matrix. We need to find such that is symmetric; that is, it satisfies

After rearrangement, (21) becomes

Because and is nonsingular, we can get from Lemma 4 that

where is arbitrary symmetric. Substituting (23) into (20) yields (v). Furthermore, and can be expressed as in (13).

Sufficiency. From the description of (i)–(vi), we can obtain that (12) holds; thus (4) holds with

Remark 6. The general solution to the Problem 1 with and prescribed eigenpairs has been given in [12], here, we generalize its solution to the case of .

Remark 7. It is complicated for the more general case , and we will discuss it in our next work. However, here we provide a simple solution. We can select the linear independent columns and the relevant eigenvalues to construct a new and , then do discussion with them.

Remark 8. When , by Theorem 5, the general solution of the Problem 1 is given by

where

with and which can be arbitrarily chosen and is arbitrarily symmetric.

Using Theorem 5, we can construct a solution to the Problem 1 as follows.

Algorithm 9. An algorithm for solving Problem 1 is proposed as follows.

(1)Input and , compute the decomposition of according to (11), and compute .(2) Choose a symmetric positive definite matrix and a symmetric matrix , arbitrarily. Compute and by (iv) and (v) in Theorem 5, respectively.(3) Choose arbitrary and , and compute by (vi) in Theorem 5, and .(4)Choose a symmetric positive definite matrix ; compute .(5)Choose arbitrary matrices and and a symmetric matrix , and form

where is given by (11). Compute and by (13).

3. Particular Solutions with

As we all know, the applications of the undamped gyroscopic system (i.e., ) exist in many fields; for details, see [5]. In this section, we will discuss the particular solutions of Problem 1 with and prescribed eigenpairs. And in this case, in (1) becomes

It is well known that the eigenvalues of (28) have a Hamiltonian structure; that is, they occur in quadruples , possibly collapsing to real or imaginary pairs or single-zero eigenvalues. In [11], Jia and Wei discussed the eigenvalues of in (28) and separated them into four categories. From the assumption (A2), we can know is even. Here we rewrite the given eigeninformation pair of Problem 1 as

with

where

are eigenpairs, and .

In this section, the Problem 1 becomes the following problem.

Given an eigeninformation pair with (29)–(31), find real matrices , , and , with being symmetric definite, being symmetric, and being skew-symmetric, so that,

As well as in Section 2, let

where is partitioned conforming with that of in (14), and we can easily calculate that also satisfies the action of in Section 2, except that has an additional property, that is, . In the following theorem, we will discuss the solubility of Problem 1 with and .

Theorem 10. Let , , and be defined as in (14), (16) and (34); then there are real matrices , , and which satisfy (33) if and only if

(i) is arbitrarily symmetric positive definite,(ii) and are arbitrary,(iii) is arbitrary symmetric,(iv),(v),(vi),

in which

with and being arbitrary real numbers.

Proof. Necessity. Same as the proof of Theorem 5, we can get

where and are arbitrary. We also have

and satisfies

After rearrangement, (38) becomes

It is easily seen that (39) has a particular solution

Next, we consider the homogeneous equation

Substituting into (41), we get

Write , where is partitioned conforming with that of in (29). Then we observe that

When , (43) can be rewritten as

in which stands for the Kronecker product and vec stands for the column vectorization of a matrix. Because , and assumption (A2), we know that is nonsingular; therefore , so .

Now we discuss the structures of matrices , , which are skew-symmetric. For simplicity, we denote by . Then we need to solve

Since has the form in (30), we can easily compute that the general solution of (45) has the form

Thus, the general solution of the homogeneous equation (41) has the form

with defined in (35). This, together with (40), gives rise to the general solution of (39)

Substituting (48) into (37) yields (v).

Sufficiency. From the description of (i)–(vi), we can obtain that (33) holds.

Remark 11. When , by Theorem 10, the solution of the Problem 1 with is given by

where

with and which can be arbitrarily chosen.

4. Particular Solution with and

In practice, the matrix in the Problem 1 with is sometimes required to be symmetric negative definite [5]. In this section, we will apply Theorem 10 to construct such a solution. We first prove the following lemma.

Lemma 12. For any given matrix defined in (35), we can construct a symmetric positive definite matrix so that defined in Theorem 10 is symmetric negative definite.

Proof. Since , it is easy to see that in Theorem 10 is symmetric negative definite if and only if the matrix

is symmetric negative definite.

By the assumption (A2), we can first construct a symmetric positive definite matrix so that . Then we use to construct the desired .

From (35) and (30), we denote

with

Here and are arbitrary real numbers. Take

with

Using (53), if we choose , , , , , and such that

then and . Obviously, such real numbers , , , , , and can be easily chosen. Once is determined, the required can be chosen by

Furthermore,

Using Lemma 12, we can construct a particular solution to the Problem 1 with , as follows.

Algorithm 13. An algorithm for solving the Problem 1 with , is proposed as follows.

(1)Input and , compute the decomposition of according to (11), and compute .(2)Choose as in (35) arbitrarily and compute by (52) and (53).(3)Construct a symmetric positive definite matrix by (54)–(57), compute by (iv), and compute by (v) in Theorem 10 or by (58).(4)Choose arbitrary and , and compute by (vi) in Theorem 10, and .(5)Choose a symmetric positive definite matrix and a symmetric negative definite matrix ; compute , .(6)Choose an arbitrary skew-symmetric matrix , and form

where is given by (11).

Remark 14. When , we only need to choose by (54)–(56), and compute by (57), that is, is the whole ; then use the same method described in Remark 8, we can obtain the particular solution of the Problem 1 with and .

5. Numerical Examples

In this section, we present two numerical examples to illustrate the solutions constructed in Sections 2 and 4, respectively. For presentation, we report all numbers in 5 significant digits only, though all calculations are carried out in full precision.

Example 1. In this example, we use Algorithm 9 to construct the general solution of the Problem 1. The partially prescribed eigeninformation as in (2)-(3) is given by the following eigenvalues and eigenvectors, which are from [12]:

It is easy to check that the matrix pair satisfy the assumptions (A1) and (A2). According to Algorithm 9, by randomly choosing

we can figure out

It is easy to check that is symmetric positive definite, and are symmetric, and is skew-symmetric. We define the residual as

and the numerical results are shown in Table 1. This shows that Algorithm 9 to construct the general solution of the Problem 1 is effective.

Table 1: Example 1.

Example 2. In this example, we use Algorithm 13 to construct the general solution of the Problem 1 with and . The partially prescribed eigeninformation as in (29)–(31) is given by randomly generated eigenvalues and eigenvectors

It is easy to check that the matrix pair satisfy the assumptions (A1) and (A2). According to Algorithm 13, by randomly choosing

and choosing

we can figure out

It is easy to check that is symmetric positive definite, is symmetric negative definite, and is skew-symmetric. We define the residual as

and the numerical results are shown in Table 2. This shows that Algorithm 13 to construct the particular solution of the Problem 1 with and is effective.

Table 2: Example 2.

6. Conclusions

In this paper, we first use techniques involving matrix decompositions to derive an expression of the general solution to the question, for a set of given pairs of complex numbers and vectors (closed under conjugation), under assumptions (A1) and (A2). Then, with the special properties and , we construct a particular solution. Numerical results illustrate these solutions.

For another case of , it is rather complex, and the proof method in Theorem 5 seems not to be used directly to find a solution of Problem 1. Fortunately, for the damped nongyroscopic system, Cai et al. [13] solved the QIEP with given eigenpairs. However, case has never been discussed in the literature for damped gyroscopic system. It might be an interesting research and needs further investigation.

Conflict of Interests

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

Acknowledgment

Guo-Liang Chen is supported by the National Natural Science Foundation of China (no. 11071079).

References

M. J. Balas, “Trends in large space structure control theory: fondest hopes, wildest dreams,” IEEE Transactions on Automatic Control, vol. 27, no. 3, pp. 522–535, 1982.

Z. Bai and Y. Su, “Soar: a second-order arnoldi method for the solution of the quadratic eigenvalue problem,” SIAM Journal on Matrix Analysis and Applications, vol. 26, no. 3, pp. 640–659, 2005.

C. Guo, “Numerical solution of a quadratic eigenvalue problem,” Linear Algebra and Its Applications, vol. 385, no. 1–3, pp. 391–406, 2004.

Z. Jia and Y. Sun, “A refined Jacobi-Davidson method for the quadratic eigenvalue problem,” in Proceedings of the 10th WSEAS International Confenrence on APPLIED MATHEMATICS, pp. 1150–3155, Dallas, Tex, USA, November 2006.

K. Meerbergen, “The quadratic arnoldi method for the solution of the quadratic eigenvalue problem,” SIAM Journal on Matrix Analysis and Applications, vol. 30, no. 4, pp. 1463–1482, 2008.

J. Qian and W. Lin, “A numerical method for quadratic eigenvalue problems of gyroscopic systems,” Journal of Sound and Vibration, vol. 306, no. 1-2, pp. 284–296, 2007.

F. Tisseur and K. Meerbergen, “The quadratic eigenvalue problem,” SIAM Review, vol. 43, no. 2, pp. 235–286, 2001.

L. Zhou, L. Bao, Y. Lin, Y. Wei, and Q. Wu, “Restarted generalized second-order krylov subspace methods for solving quadratic eigenvalue problems,” World Academy of Science, Engineering and Technology, vol. 67, pp. 429–436, 2010.

I. Gohberg, P. Lancaster, and L. Rodman, “On selfadjoint matrix polynomials,” Integral Equations and Operator Theory, vol. 2, no. 3, pp. 434–439, 1979.

M. Chu and S. Xu, “Spectral decomposition of real symmetric quadratic λ-matrices and its applications,” Mathematics of Computation, vol. 78, no. 265, pp. 293–313, 2009.

Z. Jia and M. Wei, “A real-valued spectral decomposition of the undamped gyroscopic system with applications,” SIAM Journal on Matrix Analysis and Applications, vol. 32, no. 2, pp. 584–604, 2011.

Y. Kuo, W. Lin, and S. Xu, “Solutions of the partially described inverse quadratic eigenvalue problem,” SIAM Journal on Matrix Analysis and Applications, vol. 29, no. 1, pp. 33–53, 2006.

Y. Cai, Y. Kuo, W. Lin, and S. Xu, “Solutions to a quadratic inverse eigenvalue problem,” Linear Algebra and Its Applications, vol. 430, no. 5-6, pp. 1590–1606, 2009.

Y. Yuan, “An inverse eigenvalue problem for damped gyroscopic second-order systems,” Mathematical Problems in Engineering, vol. 2009, Article ID 725616, 10 pages, 2009.

H. W. Braden, “The equations ,” Journal on Matrix Analysis and Applications, vol. 20, no. 2, pp. 295–302, 1999.