Mathematical Problems in Engineering
Volume 2018 (2018), Article ID 8050132, 8 pages
https://doi.org/10.1155/2018/8050132

Research Article

An Improved Method for Computing Eigenpair Derivatives of Damped System

Pingxin Wang,1,2 Jieer Wu,¹ and Xibei Yang³

¹School of Science, Jiangsu University of Science and Technology, Zhenjiang, Jiangsu 212003, China
²College of Mathematics and Information Science, Hebei Normal University, Shijiazhuang 050024, China
³School of Computer Science, Jiangsu University of Science and Technology, Zhenjiang 212003, China

Correspondence should be addressed to Pingxin Wang; pingxin_wang@hotmail.com

Received 6 November 2017; Revised 31 December 2017; Accepted 23 January 2018; Published 21 February 2018

Academic Editor: Roman Lewandowski

Copyright © 2018 Pingxin Wang 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.

The calculation of eigenpair derivatives plays an important role in vibroengineering. This paper presents an improved algorithm for the eigenvector derivative of the damped systems by dividing it into a particular solution and general solution of the corresponding homogeneous equation. Compared with the existing methods, the proposed algorithm can significantly reduce the condition number of the equation for particular solution. Therefore, the relative errors of the calculated solutions are notably cut down. The results on two numerical examples show that such strategy is effective in reducing the condition numbers for both distinct and repeated eigenvalues.

1. Introduction

It is well-known that eigenvalues and eigenvectors of a structure represent the dynamic characteristics of the structure. The eigenvalues are the natural frequencies of vibration, and the eigenvectors are the shapes of these vibrational modes. According to different vibrational models, eigenvalues problems can be divided into four categories.

Definition 1 (standard eigenvalue problem). Let ; the standard eigenvalue problem is to find scalars and nonzero vectors satisfying

is called the eigenvalue of and is the eigenvectors corresponding to the eigenvalue of , the set of the eigenvalues of .

Definition 2 (generalized eigenvalue problem). Let , and the generalized eigenvalue problem is to find scalars and nonzero vectors satisfying

is called the eigenvalue of matrix pencil and is the eigenvectors corresponding to the eigenvalue of .

Definition 3 (quadratic eigenvalue problem). Let , and the quadratic eigenvalue problem is to find scalars and nonzero vectors satisfying

is called the eigenvalue of and is the eigenvectors corresponding to the eigenvalue of .

Definition 4 (nonlinear eigenvalue problem). Let is a matrix function with parameter ; the nonlinear eigenvalue problem is to find scalars and nonzero vectors satisfying

is called the eigenvalue of a nd is the eigenvectors corresponding to the eigenvalue of .

In particular, when and , nonlinear eigenvalue problem (4) degenerates into standard problem (1), generalized eigenvalue problem (2), and quadratic eigenvalue (3) respectively.

Eigensensitivity analysis of eigenvalue and eigenvector plays an important role in a variety of problems, such as structural optimal design [1, 2], finite model updating [3], structural damage detection [4], and system identification [5]. One of the main tasks of eigensensitivity analysis is to compute the eigenpair derivatives with respect to structural design parameters. Methods for computing the derivatives of eigenvalues and eigenvectors have been studied by many researchers over the past four decades. Although calculating eigenvalue derivative is straightforward, determining eigenvector sensitivity raises several challenges, due in part to the singularity of the coefficient matrix of the linear equation of the eigenvector derivatives. Different numerical methods have been developed for computing derivatives of eigenvectors. Most of the methods for computing eigenpair derivatives of standard eigenvalues, generalized eigenvalues, and quadratic eigenvalue can be divided into five categories: modal method [6–12], algebraic method [6, 13–24], finite difference method [25], iterative method [26–28], and Nelson’s method [6, 29–39]. As for nonlinear eigenvalue problem, Mehrmann and Voss [40, 41] summarize the numerical solution of the general nonlinear eigenvalue problem, which can be divided into Newton method based on matrix decomposition [42–44], Newton method based on nonlinear equations [45], successive approximation method [46], subspace projection method [47–49], perimeter integration method [50–52], and so on.

Damped system arises frequently in many areas such as applied mechanics, electrical oscillation, vibroacoustics, fluid dynamics, and signal processing. The equation of motion for the free vibration of a linear damped discrete system with degrees of freedom can be expressed by

where is a design parameter and matrices , and are symmetric, respectively, mass, damping, and stiffness matrices which are analytically dependent on parameter . Suppose that ( does not depend on time ) is a solution of (5). Substituting it into (5) will lead to the following quadratic eigenvalue problem:

If and nonzero vectors satisfy (6), is called the eigenvalue and is the eigenvectors corresponding to the eigenvalue .

Some numerical methods for computing eigenpair derivatives of quadratic eigenvalue problem are presented in [17–24]. Recently, Wang and Dai [39] extended Nelson’s method to compute the eigenvector derivatives of damped system by expressing the eigenvector derivatives as addition of particular solutions and homogeneous solutions. However, when we solve a set of linear equations, we hope the corresponding coefficient matrix has small condition number [34, 53]. The equations of finding particular solutions in [39] are in different levels in their values and it thus has large condition number. In order to overcome this problem, we present an improved algorithm for computing the derivatives of eigenvectors.

The remainder of this paper is arranged as follows. In Section 2, we derive the derivatives of eigenvalues and review some related methods for computing the derivatives eigenvectors. In Section 3, we focus on dealing with the derivatives of eigenvectors and proposed an algorithm which can significantly reduce the condition number of the equation for particular solution. In Section 4, two numerical examples are performed by our proposed method. Finally, we make some concluding remarks in Section 5.

Throughout this paper, we use the following notation. denotes the set of complex matrices, , . is the identity matrix of order . stands for the diagonal matrix with diagonal elements . denotes the transpose of a matrix .

2. Methods for Sensitivity Analysis of Damped System

We assume that is the th eigenvalue of (6) and is the eigenvector corresponding to ; that is,

In this paper, we consider eigenpair derivatives of damped system with distinct or repeated eigenvalues. Without loss of generality, we suppose that are a semisimple eigenvalue with multiplicity of (6). Note that the case of corresponds to a distinct eigenvalue and the following results are also applicable.

Sun [54] has pointed out that the derivatives of repeated eigenvalues are only directionally differentiable and the eigenvectors associated with multiple eigenvalues are not necessarily continuous functions of the design parameters and are generally not derivable. In order to ensure the existence of derivatives, we assume that and are sufficiently differentiable at and eigenvalues are distinct in the neighborhood of . Under these assumptions, the eigenvectors of the quadratic eigenvalue problem (6) are uniquely determined, to within a normalizing condition, for all in the neighborhood of . Often the following normalization is adopted to ensure uniqueness of eigenvectors for simple eigenvalues:

As for repeated eigenvalues , we impose another normalization:

where .

When an eigenvalue is repeated, the associated eigenvectors are only defined up to a subspace with dimension equal to the geometric multiplicity of the eigenvalue. The computed eigenvalues of (6) at are , while the corresponding eigenvectors may not be described as above. To simplify notation, here and henceforth, “” is omitted for variables evaluated at . Let the computed linearly independent eigenvectors corresponding to be the columns of the matrix . It is natural that the eigenvector matrix satisfies (7) and the normalizing condition (9) at too. Then the adjacent eigenvectors can be expressed in terms of as follows:

where is a transformation matrix to be determined and should satisfy .

For convenience, the following notation is adopted in this study. . Differentiating (7) with respect to design parameter and letting yield

Set in (11) and assemble them into a linear system of algebraic equations:

where . Premultiplying each side of (12) by , the eigenvalue derivatives can be obtained by

Substituting (10) into (13) will derive a following eigenvalue problem:

where Equation (14) shows that diagonal elements of matrix are eigenvalues of and matrix are corresponding to right eigenvectors matrix. Hence, solving eigenproblem , yields eigenvalues derivative and matrix .

The derivatives of the eigenvectors cannot be solved directly by (11) since the coefficient matrix is singular. There are many different methods to address the singularity problem such as modal method, algebraic method, Nelson’s method, iterative method, and finite difference method. Here we briefly introduce modal method, algebraic method, and Nelson’s method.

2.1. Modal Method

The basic idea of modal method is to express each eigenvector derivative as a linear combination of all the eigenvectors:

where are the complete eigenvalues space and are sets of complex constants which are given by

However, it is difficult to obtain all eigenvectors, in practice, especially for those large and complex engineering structures. So we can only use a part of the low-order modal as the basis vector and ignore the contribution of those unknown high-order modes, which is the so-called truncated modal method. The corrections to the problem of modal truncation error for damped system have been studied by several authors [14, 15]. However, most of the existing methods obtain the derivatives of complex mode shapes of viscously damped systems by considering the state-space forms which linearize the quadratic eigenvalue problem to standard or generalized eigenvalue problem. These procedures will complicate the derivations and the practical computations since the size of the problem is doubled.

2.2. Algebraic Method

In algebraic method, the eigenvector derivatives are calculated by assembling the derivatives of eigenproblems and the additional constraints obtained from the derivative of normalization condition into a linear system of algebraic equations. In this subsection, we briefly introduce the process of algebraic method for distinct eigenvalue [6].

Suppose that is a distinct eigenvalue of (6) and is the eigenvector corresponding to . Differentiating (8) with respect to yields

Assembling (17) and the derivatives of (6) into a linear system of algebraic equation,

where the order of coefficient matrix on the left side of (18) is and the matrix on the right side is . The derivatives can be found by solving (18). The coefficient matrix can be decomposed into upper and lower triangular forms and then a forward and backward substitution scheme may be used to compute . This method has been extended to repeated eigenvalue of symmetric systems with viscous damping by [13]. However, Wu et al. [34] pointed out that there was a mistake in Lee’s extension. Later, Xu et al. [19] extended the algebraic method to compute the eigensolution sensitivities of asymmetric viscously damped systems. Li et al. [20, 21] extended the algebraic method to the more general nonviscous damped systems.

2.3. Nelson’s Method

By expressing the eigenvector derivatives as a particular solution and a homogeneous solution, Nelson [29] proposed an efficient method for computing the eigenpair derivatives for undamped systems, where only eigenpair under consideration is required. Friswell and Adhikari [31] extended Nelson’s method to symmetric and asymmetric systems with viscous damping by expressing the derivative of each eigenvector as a particular solution and a homogeneous solution of the singularity problem. Recently, Wang and Dai [39] extended Nelson’s method to compute the eigenvector derivatives of damped system with repeated eigenvalues. In [39], the solution of (12) can be expressed by

where is the general solution of the corresponding homogeneous equations. is a constant matrix. is a particular solution of (12), that is,

The particular solution may be required to be -orthogonal with respect to , that is,

The particular solution can be obtained by combining (20) and (21):

Using nonsingularity of the coefficient matrix in (22), we may solve (22) to achieve the particular solution to (12). The computation of the constant matrix in (19) employs the second-order differential information of (7). The method is an exact method and only requires the eigenvector of interest. Similar to algebraic method, Nelson’s method also needs matrix decomposition to obtain the particular solution for each eigensolution and calculate the particular solution by forward and backward substitutions.

3. The Proposed Algorithm for Eigenvector Derivatives

The algorithm presented in [39] gives a simple method for computing the eigenvector derivatives of damped system with repeated eigenvalues. As we know, we hope that the calculated solution is a close representation of the true solution when we solve a set of linear equations. This requires that the corresponding coefficient matrix has small condition number [34, 53]. However, elements in the coefficient matrix of (22) are in different levels in their values and it thus has large condition number. In order to overcome this problem, we will propose an improved method to reduce condition number of the coefficient matrix in this section.

Note that the coefficient matrix in (11) is of rank ; it cannot be inverted. We assume that the eigenvector derivative have the following form:

where is the general solution of the corresponding homogeneous equations and is a particular solution of (11); that is, satisfies

Particular solutions may be required to be -orthogonal with respect to ; that is,

Combining (24) and (25), we obtain the following linear system:

The nonsingularity of the coefficient matrix in (26) has been proved in [39]. From above discussions, we know that the part of the solution to (26) is the particular solution to (11). Using nonsingularity of the coefficient matrix in (26), we can obtain the particular solution . When we solve a set of linear equations, we hope that the calculated solution is a close representation of the true solution. This requires that the corresponding coefficient matrix has small condition number [34, 53]. However, elements in the coefficient matrix of (26) are in different levels in their values and it thus has large condition number. In order to overcome this problem, we revise (26) into the following form:

where is the th diagonal element of .

It is easy to show that (26) and (27) have the same solution and given by that solution is the particular solution of (11). However, coefficient matrices in (26) and (27) have different condition numbers. We will show that condition number of (27) is much smaller than that of (26) in numerical examples in Section 4.

Once the particular solution is obtained, coefficient can be determined by the following [39]:

where

4. Numerical Examples

In order to test the effectiveness of the proposed method, two numerical examples are employed in this section. The first one is a three-degree-of-freedom damped system with distinct eigenvalues and the second one is four-degree-of-freedom mass-spring-damper system with repeated eigenvalues. We use maximum row sum matrix norm to compute the condition numbers of (26) and (27). All codes are run in Matlab R2016a with machine precision on a personal computer.

Example 1. Consider a three-degree-of-freedom damped system shown in Figure 1. From Figure 1, we have the mass, stiffness, and damping matrices as

Figure 1: A three-degree-of-freedom damped system.

Set and is chosen as the parameter. The system has six distinct eigenvalues at , which are and and the corresponding normalized eigenvectors are

Substituting the above eigenvectors to (26) and (27), respectively, we compute the condition numbers of (26) and (27) and list the results in Table 1.

Table 1: Condition numbers of coefficient matrices in (26) and (27).

Eigenvalue
Condition number of (26)
Condition number of (27)

59.7389
9.6172

59.7389
9.6172

89.2379
11.2299

89.2379
11.2299

90.8835
2.5304

90.8835
2.5304

From Table 1, we can see that the condition numbers of (27) are considerably reduced. Therefore, the relative errors of the calculated solutions in the proposed method are notably cut down. Example 1 shows that the proposed algorithm is effective when the eigenvalues of damped are distinct. Next, we will give another example with repeated eigenvalues.

Example 2. Figure 2 shows a four-degree-of-freedom mass-spring-damper system, which is an example of [39]. We recalculate the example by using our proposed method and compare condition number of (27) with the result of (26). Let . The damping is chosen as the parameter. It is easy to verify that the system has two 2-repeated eigenvalues: and one 4-repeated eigenvalue: at . We do not consider eigenvalue because it is defective in this problem. For the eigenvalues and , the corresponding normalized eigenvectors are

Substituting the above eigenvectors to (26) and (27), we can get the condition numbers of coefficient matrices as in Table 2.

Table 2: Condition numbers of coefficient matrices in (26) and (27).

Eigenvalue
Condition number of (26)
Condition number of (27)

223.3214
22.3665

223.3214
22.3665

Figure 2: A four-degree-of-freedom mass-spring-damper system.

From Table 2, we can see that the condition numbers of coefficient matrices in (27) decrease significantly. Example 2 shows that our proposed method is effective in reducing the condition numbers for repeated eigenvalues.

5. Concluding Remarks

This paper has outlined a numerical method for computing the derivatives of eigenvalues and corresponding eigenvectors of symmetric damped system by dividing eigenvector derivatives into a particular solution and general solution of the corresponding homogeneous equation. The proposed algorithm requires only the information of those eigenvectors corresponding to the repeated eigenvalues to extend the governing equations of particular solutions. Compared with the existing methods, the proposed algorithm can significantly reduce the condition number of the governing equations of particular solutions. Therefore, the relative errors of the calculated solutions are notably cut down. Numerical examples have demonstrated the validity of the proposed method for distinct eigenvalues and repeated eigenvalues. The present study is the first step for the research of eigensensitivity. The high-order derivatives of eigenvalues and eigenvectors are specially interesting topic that remains to be further investigated.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (nos. 61503160, 11571171, 61572242, and 61773012), the Natural Science Foundation of the Jiangsu Higher Education Institutions of China (no. 15KJB110004), and the Fundamental Research Funds of Jiangsu University of Science and Technology.

References

E. J. Haug, K. K. Choi, and V. Komkov, Design Sensitivity Analysis of Structural Systems, Academic Press, Orlando, Fla, USA, 1986.

S. Adhikari, Structural Dynamics with Generalized Damping Models: Analysis, Wiley-ISTE, 2013.

Y. Y. Yuan and H. Liu, “An iterative updating method for undamped structural systems,” Meccanica, vol. 47, no. 3, pp. 699–706, 2012.

A. Messina, E. J. Williams, and T. Contursi, “Structural damage detection by a sensitivity and statistical-based method,” Journal of Sound and Vibration, vol. 216, no. 5, pp. 791–808, 1998.

S. Adhikari, Structural Dynamics with Generalized Damping Models: Identification, Wiley-ISTE, 2013.

R. L. Fox and M. P. Kapoor, “Rates of change of eigenvalues and eigenvectors,” AIAA Journal, vol. 6, no. 12, pp. 2426–2429, 1968.

S. Adhikari, “Rates of change of eigenvalues and eigenvectors in damped dynamic system,” AIAA Journal, vol. 37, no. 11, pp. 1452–1458, 1999.

S. Adhikari, “Calculation of derivative of complex modes using classical normal modes,” Computers & Structures, vol. 77, no. 6, pp. 625–633, 2000.

S. Adhikari and M. I. Friswell, “Eigenderivative analysis of asymmetric non-conservative systems,” International Journal for Numerical Methods in Engineering, vol. 51, no. 6, pp. 709–733, 2001.

S. Adhikari, “Derivative of eigensolutions of nonviscously damped linear systems,” AIAA Journal, vol. 40, no. 10, pp. 2061–2069, 2002.

I.-W. Lee and G.-H. Jung, “An efficient algebraic method for the computation of natural frequency and mode shape sensitivities. Part I. Distinct natural frequencies,” Computers & Structures, vol. 62, no. 3, pp. 429–435, 1997.

I.-W. Lee, D.-O. Kim, and G.-H. Jung, “Natural frequency and mode shape sensitivities of damped systems: Part I. Distinct natural frequencies,” Journal of Sound and Vibration, vol. 223, no. 3, pp. 399–412, 1999.

I.-W. Lee, D.-O. Kim, and G.-H. Jung, “Natural frequency and mode shape sensitivities of damped systems: part II, multiple natural frequencies,” Journal of Sound and Vibration, vol. 223, no. 3, pp. 413–424, 1999.

Y. Q. Zhao, Z. S. Liu, S. H. Chen, and G. Y. Zhang, “An accurate modal truncation method for eigenvector derivatives,” Computers & Structures, vol. 73, pp. 609–614, 1999.

Y.-J. Moon, B.-W. Kim, M.-G. Ko, and I.-W. Lee, “Modified modal methods for calculating eigenpair sensitivity of asymmetric damped system,” International Journal for Numerical Methods in Engineering, vol. 60, no. 11, pp. 1847–1860, 2004.

K.-M. Choi, H.-K. Jo, W.-H. Kim, and I.-W. Lee, “Sensitivity analysis of non-conservative eigensystems,” Journal of Sound and Vibration, vol. 274, no. 3-5, pp. 997–1011, 2004.

N. Guedria, H. Smaoui, and M. Chouchane, “A direct algebraic method for eigensolution sensitivity computation of damped asymmetric systems,” International Journal for Numerical Methods in Engineering, vol. 68, no. 6, pp. 674–689, 2006.

M. Chouchane, N. Guedria, and H. Smaoui, “Eigensensitivity computation of asymmetric damped systems using an algebraic approach,” Mechanical Systems and Signal Processing, vol. 21, no. 7, pp. 2761–2776, 2007.

Z. H. Xu, H. X. Zhong, X. W. Zhu, and B. S. Wu, “An efficient algebraic method for computing eigensolution sensitivity of asymmetric damped systems,” Journal of Sound and Vibration, vol. 327, no. 3-5, pp. 584–592, 2009.

L. Li, Y. J. Hu, X. L. Wang, and L. Ling, “Computation of eigensolution derivatives for nonviscously damped systems using the algebraic method,” AIAA Journal, vol. 50, no. 10, pp. 2282–2284, 2012.

L. Li, Y. J. Hu, X. L. Wang, and L. Ling, “Eigensensitivity analysis for asymmetric nonviscous systems,” AIAA Journal, vol. 51, no. 3, pp. 738–741, 2013.

L. Li, Y. J. Hu, and X. L. Wang, “A study on design sensitivity analysis for general nonlinear eigenproblems,” in Mechanical Systems & Signal Processing, pp. 88–105, 34, 2013.

L. Li, Y. J. Hu, and X. L. Wang, “Numerical methods for evaluating the sensitivity of element modal strain energy,” Finite Elements in Analysis and Design, vol. 64, pp. 13–23, 2013.

J. Qian, A. L. Andrew, D. L. Chu, and R. C. E. Tan, “Computing derivatives of repeated eigenvalues and corresponding eigenvectors of quadratic eigenvalue problems,” SIAM Journal on Matrix Analysis and Applications, vol. 34, no. 3, pp. 1089–1111, 2013.

T. R. Sutter, C. J. Camarda, J. L. Walsh, and H. M. Adelman, “Comparison of several methods for calculating vibration mode shape derivatives,” AIAA Journal, vol. 26, no. 12, pp. 1506–1511, 1988.

A. L. Andrew and R. C. E. Tan, “Computation of derivatives of repeated eigenvalues and corresponding eigenvectors by simultaneous iteration,” AIAA Journal, vol. 34, no. 10, pp. 2214–2216, 1996.

H. Q. Xie and H. Dai, “Simultaneous iterative method for eigenpair derivatives of damped systems,” AIAA Journal, vol. 43, no. 1, pp. 1403–1405, 2013.

H. Q. Xie, “Simultaneous iterative method for eigenpair derivatives of damped systems,” AIAA Journal, vol. 51, no. 1, pp. 236–239, 2013.

R. B. Nelson, “Simplified calculation of eigenvector derivatives,” American Institute of Aeronautics and Astronautics. Journal, vol. 14, no. 9, pp. 1201–1205, 1976.

A. L. Andrew and R. C. Tan, “Computation of derivatives of repeated eigenvalues and the corresponding eigenvectors of symmetric matrix pencils,” SIAM Journal on Matrix Analysis and Applications, vol. 20, no. 1, pp. 78–100, 1999.

M. I. Friswell and S. Adhikari, “Derivatives of complex eigenvectors using Nelson's method,” AIAA Journal, vol. 38, no. 12, pp. 2355–2357, 2000.

N. Guedria, M. Chouchane, and H. Smaoui, “Second-order eigensensitivity analysis of asymmetric damped systems using Nelson's method,” Journal of Sound and Vibration, vol. 300, no. 5, pp. 974–992, 2007.

S. Adhikari and M. I. Friswell, “Calculation of eigensolution derivatives for nonviscously damped systems using Nelson's method,” AIAA Journal, vol. 44, no. 8, pp. 1799–1806, 2006.

B. S. Wu, Z. H. Xu, and Z. G. Li, “A note on computing eigenvector derivatives with distinct and repeated eigenvalues,” Communications in Numerical Methods in Engineering, vol. 23, no. 3, pp. 241–251, 2007.

P. Omenzetter, “Sensitivity analysis of the eigenvalue problem for general dynamic systems with application to bridge deck flutter,” Journal of Engineering Mechanics, vol. 138, no. 6, pp. 675–682, 2012.

P. X. Wang and H. Dai, “Eigensensitivity of symmetric damped systems with repeated eigenvalues by generalized inverse,” Journal of Engineering Mathematics, vol. 96, pp. 201–210, 2016.

P. X. Wang and H. Dai, “Calculation of eigenpair derivatives for asymmetric damped systems with distinct and repeated eigenvalues,” International Journal for Numerical Methods in Engineering, vol. 103, no. 7, pp. 501–515, 2015.

P. X. Wang, X. Yang, and J. Mi, “Computing eigenpair derivatives of asymmetric damped system by generalized inverse,” Journal of Vibroengineering, vol. 19, no. 7, pp. 5290–5300, 2017.

P. X. Wang and H. Dai, “Calculation of eigenpair derivatives for symmetric quadratic eigenvalue problem with repeated eigenvalues,” Computational & Applied Mathematics, vol. 35, no. 1, pp. 17–28, 2016.

V. Mehrmann and H. Voss, “Nonlinear eigenvalue problems: a challenge for modern eigenvalue methods,” GAMM-Mitteilungen, vol. 27, no. 2, pp. 121–152, 2005.

H. Voss, “Projection methods for nonlinear sparse eigenvalue problems,” In Annals of the European Academy of Sciences, pp. 152–183, 2005.

R. C. Li, “Compute multiply nonlinear eigenvalues,” Journal of Computational Mathematics, vol. 10, no. 1, pp. 1–20, 1992.

H. Dai and Z.-Z. Bai, “On smooth LU decompositions with applications to solutions of nonlinear eigenvalue problems,” Journal of Computational Mathematics, vol. 28, no. 6, pp. 745–766, 2010.

J.-S. Guo, W. W. Lin, and C. S. Wang, “Nonequivalence deflation for the solution of matrix latent value problems,” Linear Algebra and its Applications, vol. 231, pp. 15–45, 1995.

A. Neumaier, “Residual inverse iteration for the nonlinear eigenvalue problem,” SIAM Journal on Numerical Analysis, vol. 22, no. 5, pp. 914–923, 1985.

A. Ruhe, “Algorithms for the nonlinear eigenvalue problem,” SIAM Journal on Numerical Analysis, vol. 10, pp. 674–689, 1973.

H. Voss, “An Arnoldi method for nonlinear eigenvalue problems,” BIT Numerical Mathematics, vol. 44, no. 2, pp. 387–401, 2004.

H. Voss, “A Jacobi-Davidson method for nonlinear and nonsymmetric eigenproblems,” Computers & Structures, vol. 85, no. 17-18, pp. 1284–1292, 2007.

R. Van Beeumen, E. Jarlebring, and W. Michiels, “A rank-exploiting infinite Arnoldi algorithm for nonlinear eigenvalue problems,” Numerical Linear Algebra with Applications, vol. 23, no. 4, pp. 607–628, 2016.

J. Asakura, T. Sakurai, H. Tadano, T. Ikegami, and K. Kimura, “A numerical method for nonlinear eigenvalue problems using contour integrals,” JSIAM Letters, vol. 1, pp. 52–55, 2009.

W.-J. Beyn, “An integral method for solving nonlinear eigenvalue problems,” Linear Algebra and its Applications, vol. 436, no. 10, pp. 3839–3863, 2012.

X.-P. Chen and H. Dai, “Stability analysis of time-delay systems using a contour integral method,” Applied Mathematics and Computation, vol. 273, pp. 390–397, 2016.

G. H. Golub and C. F. van Loan, Matrix Computations, Johns Hopkins University Press, London, UK, 3rd edition, 1996.

J. G. Sun, “Multiple eigenvalue sensitivity analysis,” Linear Algebra and its Applications, vol. 137/138, pp. 183–211, 1990.