- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Journal of Applied Mathematics
Volume 2014 (2014), Article ID 703178, 9 pages
Solutions of a Quadratic Inverse Eigenvalue Problem for Damped Gyroscopic Second-Order Systems
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.
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.
Vibrating structures such as buildings, bridges, highways, and airplanes are distributed parameter systems . 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 , Gohberg et al. developed a powerful GLR theory to solve the QIEP of the undamped gyroscopic system. In , Chu and Xu developed an elegant procedure to obtain a real-valued spectral decomposition of the damped nongyroscopic system. And then, Jia and Wei  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 , 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.  solved the QIEP with given eigenpairs. Meanwhile, for the damped gyroscopic system, Yuan  solved the QIEP with given eigenpairs. In , 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
Lemma 3 (see ). 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 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.
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 . 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 , 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.
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
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.
4. Particular Solution with and
In practice, the matrix in the Problem 1 with is sometimes required to be symmetric negative definite . In this section, we will apply Theorem 10 to construct such a solution. We first prove the following lemma.
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,
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 :
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.
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.
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.  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.
Guo-Liang Chen is supported by the National Natural Science Foundation of China (no. 11071079).
- 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.