Random Eigenvalue Problems Revisited S Adhikari Department of Aerospace Engineering, University of Bristol, Bristol, U.K. Email: S.Adhikari@bristol.ac.uk URL: http://www.aer.bris.ac.uk/contact/academic/adhikari/home.html NSSD 2005, IISc Bangalore 21 July, 2005 Random eigenvalue problems revisited – p.1/38
Outline of the Presentation Random eigenvalue problem Existing methods Exact methods Perturbation methods Asymptotic analysis of multidimensional integrals Joint moments and pdf of the natural frequencies Numerical examples & results Conclusions NSSD 2005, IISc Bangalore 21 July, 2005 Random eigenvalue problems revisited – p.2/38
Motivations Random eigenvalue problems is of fundamental interest in various branches of science and engineering: vibration of complex engineering structures (high frequency vibration) high energy physics (energy levels of atomic nuclei, quantum chaos) stability and control of structures (structural buckling with random imperfections) number theory (zeros of Reimann-Zeta function) NSSD 2005, IISc Bangalore 21 July, 2005 Random eigenvalue problems revisited – p.3/38
Random Eigenvalue Problem The random eigenvalue problem of undamped or proportionally damped linear structural systems: K ( x ) φ j = ω 2 j M ( x ) φ j (1) ω j natural frequencies; φ j mode shapes; M ( x ) ∈ R N × N mass matrix and K ( x ) ∈ R N × N stiffness matrix. x ∈ R m is random parameter vector with pdf p x ( x ) = e − L ( x ) − L ( x ) is the log-likelihood function. NSSD 2005, IISc Bangalore 21 July, 2005 Random eigenvalue problems revisited – p.4/38
The Main Issues The aim is to obtain the joint probability density function of the natural frequencies and the eigenvectors in this work we look at the joint statistics of the eigenvalues while several papers are available on the distribution of individual eigenvalues, only first-order perturbation results are available for the joint pdf of the eigenvalues NSSD 2005, IISc Bangalore 21 July, 2005 Random eigenvalue problems revisited – p.5/38
Exact Joint pdf Without any loss of generality the original eigenvalue problem can be expressed by H ( x ) ψ j = ω 2 (2) j ψ j where H ( x ) = M − 1 / 2 ( x ) K ( x ) M − 1 / 2 ( x ) ∈ R N × N ψ j = M 1 / 2 φ j and NSSD 2005, IISc Bangalore 21 July, 2005 Random eigenvalue problems revisited – p.6/38
Exact Joint pdf The joint probability (following Muirhead, 1982) density function of the natural frequencies of an N -dimensional linear positive definite dynamic system is given by π N 2 / 2 � � � ω 2 j − ω 2 p Ω ( ω 1 , ω 2 , · · · , ω N ) = i Γ( N/ 2) i<j ≤ N � � ΨΩ 2 Ψ T � p H ( d Ψ ) (3) O ( N ) where H = M − 1 / 2 KM − 1 / 2 & p H ( H ) is the pdf of H . NSSD 2005, IISc Bangalore 21 July, 2005 Random eigenvalue problems revisited – p.7/38
Eigenvalues of GOE Matrices Suppose the system matrix H is from a Gaussian orthogonal ensemble (GOE). The pdf of H : � � H 2 � � p H ( H ) = exp − θ 2 Trace + θ 1 Trace ( H ) + θ 0 The joint pdf of the natural frequencies: p Ω ( ω 1 , ω 2 , · · · , ω N ) = � � N �� � � � � � ω 2 � θ 2 ω 4 j − θ 1 ω 2 j − ω 2 exp − j − θ 0 i j =1 i<j NSSD 2005, IISc Bangalore 21 July, 2005 Random eigenvalue problems revisited – p.8/38
Eigenvalues of Wishart Matrices If H has a Wishart distribution W N ( N, λ I N ) then the joint pdf of the natural frequencies can be expressed as π N 2 / 2 p Ω ( ω 1 , ω 2 , · · · , ω N ) = (2 λ ) N 2 / 2 (Γ( N/ 2)) 2 � � N N N � � � � � − 1 1 ω 2 ω 2 j − ω 2 exp i i 2 λ ω i i =1 i =1 i<j NSSD 2005, IISc Bangalore 21 July, 2005 Random eigenvalue problems revisited – p.9/38
Limitations of the Exact Method the multidimensional integral over the orthogonal group O ( N ) is difficult to carry out in practice and exact closed-form results can be derived only for few special cases the derivation of an expression of the joint pdf of the system matrix p H ( H ) is non-trivial even if the joint pdf of the random system parameters x is known NSSD 2005, IISc Bangalore 21 July, 2005 Random eigenvalue problems revisited – p.10/38
Limitations of the Exact Method even one can overcome the previous two problems, the joint pdf of the natural frequencies given by Eq. (3) is ‘too much information’ to be useful for practical problems because it is not easy to ‘visualize’ the joint pdf in the space of N natural frequencies, and the derivation of the marginal density functions of the natural frequencies from Eq. (3) is not straightforward, especially when N is large. NSSD 2005, IISc Bangalore 21 July, 2005 Random eigenvalue problems revisited – p.11/38
Perturbation Method Taylor series expansion of ω j ( x ) about the mean x = µ ω j ( x ) ≈ ω j ( µ ) + d T ω j ( µ ) ( x − µ ) + 1 2 ( x − µ ) T D ω j ( µ ) ( x − µ ) Here d ω j ( µ ) ∈ R m and D ω j ( µ ) ∈ R m × m are respec- tively the gradient vector and the Hessian matrix of ω j ( x ) evaluated at x = µ . NSSD 2005, IISc Bangalore 21 July, 2005 Random eigenvalue problems revisited – p.12/38
Joint Statistics Joint statistics of the natural frequencies can be obtained provided it is assumed that the x is Gaussian. Assuming x ∼ N ( µ , Σ ) , first few cumulants can be obtained as � � = E [ ω j ] = ω j + 1 κ (1 , 0) 2Trace , D ω j Σ jk = E [ ω k ] = ω k + 1 κ (0 , 1) 2Trace ( D ω k Σ ) , jk �� � � = Cov ( ω j , ω k ) = 1 κ (1 , 1) + d T 2Trace ( D ω k Σ ) D ω j Σ ω j Σd ω k jk NSSD 2005, IISc Bangalore 21 July, 2005 Random eigenvalue problems revisited – p.13/38
Multidimensional Integrals We want to evaluate an m -dimensional integral over the unbounded domain R m : � m e − f ( x ) d x J = R Assume f ( x ) is smooth and at least twice differentiable The maximum contribution to this integral comes from the neighborhood where f ( x ) reaches its global minimum, say θ ∈ R m NSSD 2005, IISc Bangalore 21 July, 2005 Random eigenvalue problems revisited – p.14/38
Multidimensional Integrals Therefore, at x = θ ∂f ( x ) = 0 , ∀ k or d f ( θ ) = 0 ∂x k Expand f ( x ) in a Taylor series about θ : � � � T D f ( θ )( x − θ ) + ε ( x , θ ) f ( θ ) + 1 2 ( x − θ ) − J = m e d x R � T D f ( θ )( x − θ ) − ε ( x , θ ) d x = e − f ( θ ) 2 ( x − θ ) m e − 1 R NSSD 2005, IISc Bangalore 21 July, 2005 Random eigenvalue problems revisited – p.15/38
Multidimensional Integrals The error ε ( x , θ ) depends on higher derivatives of f ( x ) at x = θ . If they are small compared to f ( θ ) their contribution will negligible to the value of the integral. So we assume that f ( θ ) is large so that � � � � 1 � f ( θ ) D ( j ) ( f ( θ )) � � → 0 for j > 2 � where D ( j ) ( f ( θ )) is j th order derivative of f ( x ) eval- uated at x = θ . Under such assumptions ε ( x , θ ) → 0 . NSSD 2005, IISc Bangalore 21 July, 2005 Random eigenvalue problems revisited – p.16/38
Multidimensional Integrals Use the coordinate transformation: ξ = ( x − θ ) D − 1 / 2 ( θ ) f The Jacobian: � J � = � D f ( θ ) � − 1 / 2 The integral becomes: � � ξ � T ξ m � D f ( θ ) � − 1 / 2 e − 1 J ≈ e − f ( θ ) d ξ 2 R J ≈ (2 π ) m/ 2 e − f ( θ ) � D f ( θ ) � − 1 / 2 or NSSD 2005, IISc Bangalore 21 July, 2005 Random eigenvalue problems revisited – p.17/38
Moments of Single Eigenvalues An arbitrary r th order moment of the natural frequencies can be obtained from � � � µ ( r ) ω r m ω r = E j ( x ) = j ( x ) p x ( x ) d x j R � m e − ( L ( x ) − r ln ω j ( x )) d x , = r = 1 , 2 , 3 · · · R Previous result can be used by choosing f ( x ) = L ( x ) − r ln ω j ( x ) NSSD 2005, IISc Bangalore 21 July, 2005 Random eigenvalue problems revisited – p.18/38
Moments of Single Eigenvalues After some simplifications j ( θ ) e − L ( θ ) µ ( r ) ≈ (2 π ) m/ 2 ω r j � � − 1 / 2 � � � D L ( θ ) + 1 r r d L ( θ ) d L ( θ ) T − � � ω j ( θ ) D ω j ( θ ) � r = 1 , 2 , 3 , · · · θ is obtained from: d ω j ( θ ) r = ω j ( θ ) d L ( θ ) NSSD 2005, IISc Bangalore 21 July, 2005 Random eigenvalue problems revisited – p.19/38
Maximum Entropy pdf Constraints for u ∈ [0 , ∞ ] : � ∞ p ω j ( u ) du = 1 0 � ∞ u r p ω j ( u ) du = µ ( r ) j , r = 1 , 2 , 3 , · · · , n 0 Maximizing Shannon’s measure of entropy � ∞ S = − 0 p ω j ( u ) ln p ω j ( u ) du , the pdf of ω j is p ω j ( u ) = e − { ρ 0 + � n i =1 ρ i u i } = e − ρ 0 e − � n i =1 ρ i u i , u ≥ 0 NSSD 2005, IISc Bangalore 21 July, 2005 Random eigenvalue problems revisited – p.20/38
Maximum Entropy pdf Taking first two moments, the resulting pdf is a truncated Gaussian density function � � ω j ) 2 1 − ( u − � √ p ω j ( u ) = ω j /σ j ) exp 2 σ 2 2 πσ j Φ ( � j j = µ (2) where σ 2 ω 2 − � j j Ensures that the probability of any natural frequencies becoming negative is zero NSSD 2005, IISc Bangalore 21 July, 2005 Random eigenvalue problems revisited – p.21/38
Recommend
More recommend
