Stochastic Averaging of Oscillators Excited by Colored Gaussian. Processes. 3 R. Valery Roy
|
|
- Lydia Patrick
- 5 years ago
- Views:
Transcription
1 Stochastic Averaging of Oscillators Excite by Colore Gaussian Processes 3 R. Valery Roy Department of Mechanical Engineering, University of Delaware, Newark, Delaware 976. Abstract: The metho of stochastic averaging has been evelope an applie in the past mainly base on Stratonovich-Khasminskii theorem. We examine in this paper the application of this metho in the case of arbitrary colore Gaussian excitations which can be consiere as output of multiimensional linear lters to white Gaussian noise. The metho use is base on a perturbation theoretic approach of the Fokker-Planck-Kolmogorov equation which governs the response probability ensity function. First, for oscillators with linear elastic forces an non-parametric excitation, it is shown that, to leaing orer of perturbation, the results obtaine matches those erive by application of Stratonovich-Khasminskii theorem in the case of broa-ban excitation. Then, more general results are erive for nearly amiltonian systems perturbe by parametric excitations of uncorrelate colore noises. t is shown that the state probability ensity function is governe by a reuce equation in the \slow" amiltonian variable only, which epens on a number of parameters characterizing the colore noise excitations. Several examples are given for illustration. As a preliminary to these theoretical evelopments, the problem of etermining the eigenfunctions an eigenvalues of the Fokker-Planck operator is aresse for a general class of linear multiimensional systems.. ntrouction n the quest for approximate solution techniques for ranom vibration problems, the metho of stochastic averaging has proven to be a powerful analytic tool. First propose by Stratonovich [] in 963, various mathematical justications, erivations an extensions have appeare in the literature. See for example Khasminskii [2] (966), Roberts [3] (978), hu [4] (983), Lin [5] (986), hu an Lin [6] (99), Re-orse an Spanos [7] (992). Comprehensive reviews attesting to the success of stochastic averaging have been written by Roberts an Spanos [8] (986), an by hu [9] (988). The metho is applicable to systems with light amping forces an excite by broa-ban, quasi-white ranom processes. The stanar stochastic 3 nternational Journal of Non-Linear Mechanics, 29, pp , 994.
2 averaging metho consists of rst transforming the original equations of motion into equations for slowly varying processes an rapily varying processes. Equations for the \slow" variables only can then be obtaine by subsequent averaging, using the Stratonovich-Khasminskii limit theorem; a corresponing Fokker-Planck- Kolmogorov (FPK) equation can be written for the joint probability ensity function (p..f.) of the \slow" variables. ence, the \fast variables" are eliminate. For systems with a single egree of freeom, the approach generally leas to a more tractable one-imensional problem in the energy variable. Another version of stochastic averaging is base on a perturbation theoretic approach to the Fokker-Planck- Kolmogorov equation. A perturbation solution of the bounary value problem governing the state p..f. is sought in the asymptotic limit of small friction forces an small noise excitations. For vanishing value of the perturbation parameter, the system is conservative. Such a perturbation scheme has been use by Blankenship an Papanicolaou [](978), orsthemke an Lefever [] (98), Klosek-Dygas, Matkowsky, an Schuss [2, 3] (988). The purpose of this paper is two-fol: rst to introuce this technique to the engineering mechanics community as an alternative to Stratonovich-Khasminskii theorem, an secon to generalize the results of Klosek-Dygas, Matkowsky, an Schuss [2, 3] to the case of colore Gaussian noise excitations which can be treate as outputs of multiimensional linear lters to white noise. Since the temporal evolution of the system is non-markovian in the case of non-white noise excitation, the response of the combine systemplus-lter must be consiere. Central to the erivations of this paper is the expansion of each perturbation coecient along the eigenfunctions of the Fokker-Planck operator associate with the linear lter. ence in x2, we aress the problem of etermining these eigenfunctions an the corresponing eigenvalues. t is shown that for general multiimensional linear systems satisfying the conition of etaile balance, analytical expressions can be foun in terms of multivariate ermite polynomials. Note that the results of x2 are useful in themselves for purposes other than those of x3 an x4, e. g., for approximate solutions of the FPK equation base on a Bubnov-Galerkin scheme. n x3, the metho is rst applie to the class of oscillators with linear restoring force an non-parametric excitation. The unperturbe system correspons to the classical unampe, unforce linear oscillator. The problem is solve by rst transforming the joint process (x; _x) into amplitue an phase processes (a; ). t is shown that, to leaing orer of perturbation, the results obtaine are ientical to the stanar averaging of the amplitue process by use of Stratonovich-Khasminskii theorem: only the value of the excitation power spectrum at the oscillator's natural frequency is taken into account. n x4, we treat the case of nearly amiltonian systems perturbe by weak friction forces an by weak colore parametric noise excitations. We follow here the methoology of [2] by using a irect perturbation of the FPK equation in the (x; y; u) variables without the use of a transformation to slowly varying processes. The state probability ensity function is shown to be governe, to leaing orer of perturbation, by an FPK equation in the amiltonian variable only. This equation epens on parameters characteristic of the colore noise excitations. Several examples are given. n particular, the results of references [2] an [4, 5] are recovere. 2. Preliminaries: The Eigenvalues an Eigenfunctions of the Fokker-Planck Operator for Linear Systems Consier the general stochastic linear system governe by the equation _u(t) = Au(t) + W(t); () where u = (u ; : : : ; u N ) T. The excitation W(t) is an N-imensional, zero-mean, white Gaussian noise such that hw(t + )W T (t)i = B(). B is assume to be an N 2 N symmetric, non-negative matrix of constant coecients b ij. The N 2 N matrix A of constant, real or complex elements a ij is assume to have a complete set of eigenvalues ; : : : ; N with corresponing eigenvectors e ; : : : ; e N. The eigenvalues ( i ) are assume to all have a negative real part: <( i ) < ; i = ; : : : ; N: (2) 2
3 With these assumptions, Eq.() has a stationary, zero-mean, Gaussian, Markov solution process whose corresponing probability ensity function is given by p s (u) = N exp( 2 ut K u u) (3) where the normalization constant is N = (2) N=2 (et K u ) =2. The matrix K u stans for the stationary response covariance matrix huu T i an is known to be solution of the matrix equation A K u + K u A T = B: (4) We also enote by U the N 2 N matrix whose columns are e ; : : : ; e N an such that U A U = D = iag( ; : : : ; N ): (5) To Eq.() correspons two linear ierential operators, the (forwar) Fokker-Planck operator Lu an its formal ajoint (backwar Kolmogorov operator) L 3 u whose expressions are given by: Lu ui a ij u j + 2 b ij@ 2 u iu j ; (6) L 3 u = a ij u ui + 2 b ij@ 2 u iu j ; (7) where the usual summation convention is use for repeate inices. We seek in this section analytical expressions for the eigenfunctions 9 (u) an 9 3 (u) an the corresponing eigenvalues an of, respectively, Lu an L 3 u ; that is, we seek the solutions of the equations Lu9 = 9 ; L 3 u9 3 = 9 3 : (8) n general, the two sets (9 ) an (9 3 ) are unrelate, except for the conition of biorthogonality u9 (u)9 3 R n (u) = ; if 6= ; (9) in the case of a iscrete spectrum of eigenvalues. There exists however a conition for which the two operators Lu an L 3 u are relate to one another: the conition of etaile balance. A stationary system is sai to satisfy the principle of etaile balance, if the joint probability ensity of the transition from (u; t) to (u ; t + ) is equal to that of the time reverse transition from (u ; t) to (u; t + ). The reverse transition is such that the variable u = (u ; : : : ; u N ) T is transforme into the variable u = ( u ; : : : ; N u N ) T, where i = for an even variable u i (e.g., a isplacement variable), an i = for an o variable u i (e.g., a velocity variable) uner time reversal. enotes the N 2 N iagonal matrix iag( ; : : : ; N ). See Gariner [6] (983) for more etail. Thus, etaile balance is equivalent to, for Markov systems, p(u ; ju; )p s (u) = p(u; ju ; )p s (u ) () an p s (u) = p s (u): () t can be shown [6] that, for the general linear system given by Eq.(), these conitions are equivalent to the Onsager relations given by (AK u ) = (AK u ) T ; (2) B = B: (3) Uner these conitions, if 9 3 is an eigenfunction of L 3 u corresponing to the eigenvalue, then 9(u) = p s (u) 9 3 (u) is an eigenfunction of Lu corresponing to the same eigenvalue. See [6]. Thus, the problem is reuce to that of solving the equation L 3 u9 3 =
4 The solution of the ajoint eigenvalue problem is simplie by using the change of variable u = Uv which iagonalizes the matrix A. The equation L 3 u9 3 (u) = 9 3 (u) becomes U li a ij U jk v vl U li b ij U v kv l 9 3 = 9 3 (4) that is, using U A U = D, k v vk ~ b 2 v kv l 9 3 = 9 3 (5) where the ~ b ij 's are the elements of the symmetric matrix e B = U BU T. The solutions of Eq.(5) are foun by examining the set of multivariate ermite polynomials of the variable v as ene by Appell an Kampe e Feriet (see Erelyi [7], 953): m;:::;mn (v) = () m exp( 2 vt Cv)@ m v m ;:::;v m N exp( N 2 vt Cv) (6) G m;:::;mn (v) = () m exp( 2 vt Cv)@ m w m ;:::;w m N exp( N 2 vt Cv) (7) where w = Cv an C is a positive enite, symmetric N 2 N matrix. Each polynomial G m;:::;mn (v) of orer m = m + + m N is known to satisfy the following system of partial ierential equations (with no summation on the inex k): v vk G l= 2 v kv l G = m k G; k = ; : : : ; N (8) where the ij are the elements of C. ence, G m;:::;m N (v) also satises the following linear combination of these equations: k v vk G k 2 v kv l G = m k k G: (9) Thus, one can reaily see that G m;:::;m N (v) will satisfy equation (5) if k kl + l lk = ~ b kl, that is, if the matrix C is chosen to satisfy D C + C D = e B (2) whose unique solution is C = U T K u U: (2) Thus, C is the covariance matrix K v of the process v = U u whose normalize p..f. is given by 9 (u) = N exp( 2 ut K u u) = N exp( 2 vt K v v) = N exp( 2 wt K v w): (22) with w = K v v = K v U u. Therefore, the eigenfunctions 9(u) are given by with corresponing eigenvalues given by 9 m;:::;m N (u) = 9 (u) 9 3 m ;:::;m N (u) = () m w m ;:::;w m N N 9 (u) (23) m;:::;m N = m + + m N N : (24) For m =, 9 (u) is the stationary p..f. p s (u) corresponing to the eigenvalue =. For m =, the eigenfunctions are reaily foun to be (u) 9 2 (u). 9 N (u) = r w exp( 2 wt K v w) = K v w9 (u) = U u9 (u): (25) 4
5 Since the set of ermite polynomials forms a complete set, the transition probability ensity of the process u(t) can be expane in the form p(u; ju ) = X m =;:::;mn = c m;:::;m N e m;:::;mn 9 m;:::;m N (u)9 3 m ;:::;m N (u ); ( ) (26) The coecients c m;:::;mn are etermine by taking into account the initial conition lim! p(u; ju ) = (u u ) an the biorthogonality conition (9), leaing to: c m;:::;m N = u 9 m;:::;m N (u)9 3 m ;:::;m N (u) : (27) Then, various statistical response functions can be etermine, such as the response stationary correlation matrix R u () = hu(t)u T (t + )i (for ): R u () = where the vector r m;:::;m N = u X u u u T p(u; ju ) p s (u ) m =;:::;mn = c m;:::;m N r m;:::;m N r T m ;:::;m N exp( m;:::;m N ) (28) is ene as r m;:::;m N = u u 9 3 m ;:::;m N (u)9 (u) = u u 9 m;:::;m N (u): (29) The class of linear systems satisfying the conition of etaile balance can be shown to contain the following linear lters of arbitrary orer N: u (N) + a N u (N) + + a 3 u (2) + a 2 u () + a u = W (t) (3) where u (i) = i u=t i, an the constant coecients (a i ) are such that Eq.(3) amits a stationary solution. This equation can be set in the state-space form () by ening u = (u; _u; : : : ; u (N) ) T, with A an B such that, with hw (t)w (t + )i = (), a ij = 8 < : ; for j = i + a j ; for i = N ; otherwise b ij = ; for i = j = N ; otherwise Note that u, u, u (4), : : : are even variables, whereas _u, u (3), : : : are o variables: = iag(; ; : : : ; () N+ ). t can then be verie that the conitions of etaile balance (2) an (3) are satise, thus generalizing the results of Atkinson [8] (973) for the eigenfunctions corresponing to a lter of orer Stochastic Averaging of Oscillators with Non-Parametric Colore Excitation We seek here an asymptotic solution of the probability ensity function of the nonlinear oscillator (3) x(t) + f(x(t); _x(t)) +! 2 o x(t) = p u(t) (32) excite by the colore ranom process u(t), in the limit! of weak noise an weak nonlinearity. The ranom process u(t) is assume to be a stationary, physically realizable Gaussian process, with zero-mean. More specically, it is assume that a Markov realization of this process can be constructe, that is, u(t) is a component of an N-imensional stationary iusion Markov process u(t) = (u = u; u 2 ; : : : ; u N ) T governe by Eq.(). The conition of stability (2), an the conitions of etaile balance (2) an (3) are assume satise, so that the complete sets of eigenvalues an eigenfunctions of the corresponing operator Lu are 5
6 given by Eq.(24) an Eq.(23) respectively. As in x2, 9 (u) enotes the stationary p..f. of u, an correspons to the eigenvalue =. We also enote by R u the correlation function of u(t), that is, the (; ) element of the correlation matrix Ru of u(t). As a rst step, the joint response process (x; _x) is transforme into (a; ) accoring to the relationships x(t) = a(t) cos (t); _x(t) =! o a(t) sin (t): (33) The equations of motion for the amplitue a(t) an phase (t) processes are then given by _a(t) = p f(a cos ;! o a sin ) sin u(t) sin ; (34)! o! o _(t) =! o +! o a f(a cos ;! oa sin )cos p u(t) cos : (35)! o a n the limit!, the process a(t) is clearly a slow variable, whereas the process (t) is a fast variable. The amplitue of the response will not change signicantly over a time interval of the orer O(). To investigate the response of the system over a time interval of the orer O( ), the time variable is scale as t = =, where is a slow time scale. The extene process (a; ; u) is a Markov process, whose probability ensity p (a; ; u; ), assume to exist, is solution of the following FPK equation where the operators L, L an L 2 are given by an Lp = L + p L + L 2 p = (36) L =! + Lu; (37) L = u! o a + u! o cos (38) L 2 =! a (f sin )! o (f cos : (39) A solution of this equation is sought in terms of a perturbation expansion of the parameter in the following form p (a; ; u) = p (a; ; u) + p p (a; ; u) + p 2 (a; ; u) + (4) where, for simplicity, the epenence upon the variable is hereon omitte. nserting (4) into (36) an ientifying terms of equal power of =2 leas to the sequence of equations : L p = ; (4) p : L p = L p ; (42) i=2 : L p i = L p i L 2 p i2 ; i = 2; 3; : : : (43) Each term of expansion (4) is require to be 2-perioic in the variable, integrable over a 2 [; [, an assume to satisfy natural bounary conitions for the variable u. The normalization conition of p is equivalent to 2 a u p (a; ; u) = ; (44) an to, for i, 2 a u p i (a; ; u) = : (45) 6
7 n general, a linear equation of the type L p = q, such as Eq.(42) or Eq.(43), will amit a solution p if its inhomogeneous part q is orthogonal to all elements of the null space ker L 3 of the ajoint of L. nee, for all q 3 2 ker L 3, that is, such that L3 q3 =, the following solvability conition is obtaine: a u q 3 q = a u q 3 L p = a u pl 3 q3 = (46) Thus, to each equation (42) an (43) correspons a solvability conition as given by Eq.(46). Since the operator L 3 =! o@ + L 3 u is the sum of an operator acting on the variable only, an of the operator L 3 u acting on the variable u only, the elements q 3 of ker L 3 are sought in a series expansion along the eigenfunctions 9 3 m ;:::;m N (u) of L 3 u in the following way q 3 (a; ; u) = X m =;:::;mn = q 3 m ;:::;m N (a; )9 3 m ;:::;m N (u): (47) Then, each coecient q 3 m ;:::;m N is solution of Since each eigenvalue satises < m;:::;mn non-zero coecient in expansion (47) is (! + m;:::;m N )q 3 m ;:::;m N = (48) <, except for = corresponing to 9 3 (u) =, the only perioic q 3 (a; ) = q 3 (a): (49) Therefore, the elements of ker L 3 are arbitrary functions of the amplitue a only. This implies that the solvability conition (46) of equation L p = q takes the expression 2 2 u q(a; ; u) = : (5) which imposes that the average of q over the variables, an u is zero. We then procee to n the leaing term p of perturbation expansion (4). p (a; ; u) is solution of Eq.(4) an can be sought in a series expansion of the eigenfunctions of the operator Lu, that is, Then, each term p () m ;:::;m N p (a; ; u) = is solution of X m =;:::;mn = p () = p () 9 + k= m ;:::;m N (a; )9 m;:::;m N (u) p () k 9 k + k;l= p () kl 9 kl + (5) (! + m;:::;m N )p () m ;:::;m N (a; ) = ; (52) an, upon integration over the variable, is a non-zero, perioic function of if an only if m = corresponing to =. t follows that the solution of Eq.(4) is inepenent of the phase variable an takes the form p (a; ; u) = p (a)9 (u); (53) where p (a) remains to be foun from subsequent solvability conitions. To orer O( ), the amplitue process a(t) is uncouple from the phase process (t). At orer p, the term p is solution of the equation (! + Lu)p = L p = u! o 9 (u)(p a p ) sin ; (54) 7
8 whose solvability conition (46) is clearly always satise. As was one for p, p is etermine as the eigenfunction expansion p (a; ; u) = X m =;:::;mn = p () = p () 9 + k= m ;:::;m N (a; )9 m;:::;m N (u) p () k 9 k + k;l= p () kl 9 kl + (55) We rst note that the function u 9 (u) can be expresse as a linear combination of eigenfunctions (9 k ) k=;:::;n of orer m = : u 9 (u) = k= k 9 k (u): (56) where the scalar coecient k is in fact the rst element of the vector r k =c k ene in x2 for the orer m =. Then, the coecients p () m ;:::;m N are solution of (! + m;:::;m N )p () m ;:::;m N (a; ) = k! o (p a p ) sin ; for m = ; ; otherwise (57) The only non-zero perioic solutions are corresponing to m = an =, an p () (a; ) = p() (a) (58) p () k (a; ) = k (p! o a p )(cos + k sin ) (59) k! o corresponing to m = an k = k, for k = ; : : : ; N. The coecient p () (a) correspons to the homogeneous solution of Eq.(54), an is so far unetermine. Finally, p takes the expression p (a; ; u) = p () (a)9 (u)(p a p ) k= k! 2 o + 2 k (cos + k! o sin )9 k (u) (6) R which satises the normality conition (45) if p () (a)a =. The analytical expression for p (a; ) is etermine by application of the solvability conition of equation (43) at the orer, that is, We rst n where F (a) = 2 u (L 2 p + L p ) (a; ; u) = : (6) u L 2 p (a; ; u) p! a (F (a) p ); (62) 2 f(a cos ;! o a sin ) sin : (63) Then, using Eq.(6), an noting that the term p () (a)9 (u) oes not contribute, where 2 2 u L p (a; ; u) = D a ( 4! o 2 a p a p ); (64) D u = k= 2 k k k : (65)! o k 8
9 The scalars k enote the following averages of the process u: k = u u 9 k (u); k = ; : : : N: (66) The scalar D u can in fact be expresse in terms of the excitation power spectral ensity by using the expression of the correlation function (28) an the biorthogonality property of eigenfunctions: S u (!) = R u ()e i! = k= 2 k k k! k ence, D u = S u (! o ). Finally, it is seen that p (a; ) satises the following FPK-type p a F (a) + S u(! o )! o 4! o 2 a (67) p + S u(! o 2 4! o 2 aa p : (68) Thus, it is seen that, to orer O( ), the response probability ensity function is governe by an FPK equation ientical to that erive from Stratonovich-Khasminskii theorem for a broa-ban excitation u(t). To this orer, only the value of the excitation spectral ensity at the oscillator's natural frequency is taken into account, as woul be the case in a white noise approximation of u(t). The joint probability ensity of (x; _x) is given by, to orer O( ), p (x; _x; ) = 2! o a p (a; ; ) = 2! o a p (a; ) (69) an, is therefore inepenent of the phase angle. At the next orer of perturbation the term p () (a) woul be etermine from the solvability conition of the equation L p 3 = L p 2 L 2 p after rst solving for p 2 (a; ; u) up to an unknown term p (2) (a). These higher orer perturbation terms woul take further account of the excitation spectrum an yiel phase-epenent contributions. 4. Stochastic Averaging of Nearly-amiltonian Systems Perturbe by Parametric Colore Excitation We consier in this section systems escribe by a amiltonian function (x; y) an perturbe by weak friction forces an weak parametric colore noise, an whose equations of motion are consiere in the form x + f(x; y) + p (x; y) u(t) y + g(x; y) + p 2 (x; y) u(t): (7) We seek an asymptotic expression for the joint probability ensity of (x(t); y(t)) in the limit!. t is assume that the unperturbe trajectories of the system, that is, for =, are perioic functions of time an are ene in the x; y plane by the close contours of constant values of the amiltonian (x; y). The excitation u(t) is chosen as in x3. For small, the probabilistic behavior of the response is not expecte to change signicantly over time intervals of the orer O(). After scaling the time variable as t = =, the probability ensity of the extene process (x; y; u) is governe by Eq.(36) where the operators L, L, an L 2 are now given L x y + @ y + Lu (7) L = x (x; y) y 2 (x; y) (72) L 2 x f(x; y g(x; : (73) 9
10 Again, we seek an expression for the leaing term p (x; y; u) (omitting the -epenence) of the perturbation expansion (4), solution of the equation L p =. After expaning p on the complete set of eigenfunctions 9 m;:::;m N (u) of Lu, we n that each coecient is solution of the linear partial @ @ y + m;:::;m N p () m ;:::;m N = ; (74) which, accoring to the metho of characteristics, is equivalent to the system of equations x y @x m;:::;m N )p () m ;:::;m N = : (75) The rst two equations ene the contours of constant -values, an their solutions correspon to the unperturbe trajectories, an hence are perioic functions x(; ), y(; ) of the variable. Clearly, the only perioic function p () m ;:::;m N solution of Eq.(75) correspons to the eigenvalue =, an hence, is inepenent of the variable. Thus, the solution of L p = must be of the form p (x; y; u) = p ((x; y)) 9 (u); (76) where p () remains to be etermine from subsequent solvability conitions. ere, plays the role of a slow variable, an that of the fast variable. At the next orer, the solvability conition of equation L p = L p can be state as x y uq 3 L p = uq 3 L p = (77) for all q 3 such that L3 q3 =. We have use the fact that the Jacobian of the transformation from the (x; y) variables to the (; ) variables is equal to. Furthermore, since the elements q 3 are arbitrary functions of the variable only, this conition implies that ul p = @ @ y + m;:::;m N p () m ;:::;m N = where R the contour integral is evaluate on the curves of constant. This conition is always satise since uu 9 (u) =. An expression for p is then sought as an expansion on the set of eigenfunctions (9 m;:::;mn (u)), leaing to the following equation for the expansion coecients: k (@ x ( p ) y ( 2 p )) ; for m = ; (79) ; otherwise which is equivalent to, on each characteristic curve (x(; ), y(; )) (the unperturbe trajectories of the system): m ;:::;m N p () m ;:::;m N = k (@ x y 2 )p k x + y )@ p ; for m = ; ; otherwise ere, the constants k, (k = ; 2; : : : ; N), are ene as in x3. The only non-zero, perioic solutions correspons to m = m + + m N = an. For m =, the solution is function of only, p () (), an remains to be etermine. For m =, Eq.(8) can be viewe as that of a linear system with a perioic input, an whose response is perioic only at steay state. Thus, the perioic solution of Eq.(8) is given by the following (8)
11 convolution integrals (k = ; 2; : : : ; N): where p () k (; ) = p () k (; p J k (; ); (8) k (; ) = k e k (@ x e y e 2 ); (82) J k (; ) = k e k (e + e 2 ); (83) e = (x(; ); y(; )) ; (84) e 2 = 2 (x(; ); y(; )) ; (85) g@ x x (x(; ); y(; )) ; (86) g@ y y (x(; ); y(; )) : (87) Finally, the solvability conition of the next orer equation L p 2 = L 2 p L p is state as u (L 2 p + L p ) = ; (88) which leas successively to, noting that the unetermine term p ( )() oes not contribute, with an nally to with u L 2 p p () F k () = G k () = x + y p ; (89) ul fp () k (; )9 k (u)g = k fp () F k ()g (9) x + y ) e k (@ x e y e 2 ) (9) ul f(@ p )J k (; )9 k (u)g = k fg k ()@ p g ; (92) x + y ) e k (e + e 2 ): (93) ence, this solvability conition leas to the following FPK equation for the joint probability ensity p (x; y) = p ((x; y)) T p fp (f 3 () + F ())g fg()@ p g (94) where T () (the perio of the unperturbe trajectories), f 3 (), F (), an G() are given by f 3 () = F () = G() = T () = (95) x + y ) (96) k= k= k k F k () (97) k k G k (): (98)
12 All these expressions are compute from contour integrals over the noise-free an friction-free trajectory escribe by the functions (x(; ), y(; )), at xe value of the variable. Equation (94) can be easily generalize in the case of multiple, uncorrelate, parametric colore noise excitations: each noise excitation will contribute inepenently to the (noise-inuce) rift F () an iusion G() functions. We illustrate the calculation of Eq.(94) with the example of a general, linear parametric oscillator governe by the equation x + 2! o ( + p 3 (t)) _x +! 2 o( + p 2 (t))x = p (t); (99) where, 2, an 3 are taken to be mutually uncorrelate, colore, Gaussian noise excitations. This oscillator has been investigate by Ariaratnam an Tam [9] (979). The unperturbe trajectories are given by the functions x(; ) = =! o (2) =2 cos! o ; y(; ) = (2) =2 sin! o () of perio T () = T = 2=! o, an where = y 2 =2 +! 2 o x2 =2. Then, using Eq.(67) for the excitation power spectrum, we n f 3 () = 2! o T F () = ; G () = T 2 N X k= F 2 () = ; G 2 () =!2 4 T 2 F 3 () = 2 2! 2 o T N X k= G 3 () = 2! 2 o T 2 N X k= 2 3;k 3;k 3;k 2 ;k ;k ;k = T! o ;k 2 S (! o ) k= 2 2;k 2;k 2;k =!2 4! o ;k 4 T S 2 (2! o ) 2 () = 2 2! 2 o T S 3 () 3;k 3;k 4 3;k + 2 3;k 4! 2 o + 2 3;k! = 2! 2 o T (2S 3 () + S 3 (2! o )) 2 where F i an G i correspon to noise i, i =, 2, an 3 an have etermine accoring to Eqs.(97) an (98). Using these expressions in Eq.(94), the stationary ensity is reaily etermine: p (x; y) = where N is a normalization constant, an where N (c + c 2 (y 2 +! 2 o x2 )) c3=c2 (2) c = S (! o ); c 2 =!2 4 S 2 (2! o ) + 2! 2 o(2s 3 () + S 3 (2! o )); c 3 = 2! o ( +! o S 3 ()): (3) This solution will exist only if c 3 =c 2 >, that is, if 2! o >!2 4 S 2 (2! o ) + 2! 2 o S 3 (2! o ): (4) The results of references [2, 3] are recovere in the case of nonlinear oscillators governe by the equation x + f(x; _x) _x + U (x) = p (x; _x) u(t); (5) in which case, becomes the energy variable E = y 2 =2 + U(x) with y = _x. Note that in references [2, 3] the mean rst-passage time from energy level E to E c, characteristic of the system reliability, is shown to be governe by an \average" equation by aopting a similar perturbation scheme. n the case of non-parametric excitation, that is, where (x; _x) =, the functions F (E) an G(E) of Eq.(36) can be simplie further, after 2
13 expaning the perioic velocity response y(e; ) of the unperturbe system, that is, of 2 x= 2 + U (x) = in Fourier series of perio T (E) = 2=!(E): y(e; ) = _x(e; ) = X y n (E) exp(i n!(e)): (6) n= Then, by substituting this series expansion in Eq.(9) an Eq.(93) we obtain, with =, 2 =, E y = y, an for k = ; 2; : : :; N: F k (E) = ; (7) an ence, F (E) =, an G k () = E = T (E) X y(e; ) e k y(e; ) X n= k k G(E) = T (E) y n (E)y n (E) n= k + in!(e) k= = T (E) 2 jy (E)j 2 S u () + y n (E)y n (E) k + in!(e) : (8) X n= jy n (E)j 2 S u (n!(e))! ; (9) where S u (!) is the excitation power spectrum, an jy n j 2 = y n y n. ence, the results of Roberts [4, 5] are recovere. Similar expressions can be obtaine in the case of parametric excitations. For example, for the oscillator x + f(x; _x) _x + U (x) = p ( (t) + x 2 (t) + _x 3 (t)) ; () where, 2, an 3 are three inepenent Gaussian colore noise excitations, the corresponing functions F i an G i, i = ; 2; 3, are given by F (E) = ; () G (E) = T (E) 2 jy (E)j 2 S () + X n= jy n (E)j 2 S (n!(e))! ; (2) an F 2 (E) = ; (3) G 2 (E) = T (E) 2 jz (E)j 2 S 2 () + X n= jz n (E)j 2 S 2 (n!(e))! ; (4) F 3 (E) = 2 T (E) w (E)S 3 (); (5) G 3 (E) = T (E) 2 jw (E)j 2 S 3 () + X n= jw n (E)j 2 S 3 (n!(e))! : (6) n these expressions, y n, z n, an w n are the Fourier coecients of y(e; ), x(e; )y(e; ), an y 2 (E; ) respectively, corresponing to the unperturbe system. 3
14 Conclusions Moeling excitations by colore noises is justie by the fact that real life uctuations of the environment have a non-vanishing correlation time an a non-white power spectrum. ence, the results of x3 an x4 are useful to stuy the eect of nite correlation time of the excitation on the system's behavior. We have foun that, to leaing orer of perturbation, the state probability ensity function is governe by a reuce equation in the \slow" variable only, that is, the amplitue variable for the harmonic oscillator, or the amiltonian variable for general amiltonian systems. n some cases, we have foun that the results obtaine are ientical to those erive in the past by application of Stratonovich-Khasminskii theorem. n such cases, the assumption of etaile balance for the excitation process u was only neee to the extent that the existence of the eigenfunctions an eigenvalues is then guarantee, but the nal results such as given by Eq. (68) o not explicitly epen on the etaile knowlege of the eigenfunctions an eigenvalues of the noise process. owever, for the more general case of amiltonian systems perturbe by parametric noise, the result of Eq.(94) oes require the assumption of etaile balance in orer to etermine in this equation certain parameters characterizing the ranom excitation. But we surmise that even in this case the results woul be ientical to those obtaine by using Stratonovich-Khasminskii theorem, an shoul only be vali for broa-ban excitations, that is, when the ratio of the relaxation time of the system to the correlation time of the ranom excitation is much larger than one. A rawback of the perturbation scheme is that it oes not impose any restrictions on the shape of the excitation power spectrum. We observe though that narrow-ban excitations corresponing to one or more eigenvalues m;:::;mn close to the natural frequency of the system (near resonance) woul generate nearly singular terms in the present erivations. Thus other asymptotic methos shoul be evelope in the case of narrow-ban ranom excitations. Acknowlegements Many thanks to one of the reviewers for his careful review. References [] Stratonovich, R. L., Topics in the Theory of Ranom Noise, Goron an Breach, New York, (963). [2] Khasminskii, R.., \A Limit Theorem for the Solutions of Dierential Equations with Ranom Right- an Sies", Theory of Probability an Applications, Vol., No. 3, pp , (966). [3] Roberts, J. B., \The Energy Envelope of a Ranomly Excite Non-Linear Oscillator", Journal of Soun an Vibration, Vol. 6, pp , (978). [4] hu, W. Q., \Stochastic Averaging of the Energy Envelope of Nearly Lyapunov Systems", Proceeings of the UTAM Symposium on Ranom Vibrations an Reliability, (Frankfurt/Oer, G. D. R., 982), (K. ennig, E.), pp , Akaemie, Berlin, (983). [5] Lin, Y. K., \Some observations on the Stochastic Averaging Metho", Probabilistic Engineering Mechanics, Vol., No., pp , (986). [6] hu, W. Q., an Lin, Y. K., \On Stochastic Averaging of Energy Envelope", Journal of Engineering Mechanics, ASCE, Vol. 7, No. 8, pp , (99). [7] Re-orse, J. R., an Spanos, P. D., \Generalization to Stochastic Averaging in Ranom Vibration", nternational Journal of Non-Linear Mechanics, Vol. 27, No., pp. 85-, (992). [8] Roberts, J. B., an Spanos, P. D., \Stochastic Averaging: An Approximate Metho for Solving Ranom Vibration Problems", nternational Journal of Non-Linear Mechanics, Vol. 2, No. 2, pp. -34, (986). 4
15 [9] hu, W. Q., \Stochastic Averaging in Ranom Vibration", Applie Mechanics Reviews, Vol. 4, No. 5, 89-99, (988). [] Blankenship, G., an Papanicolaou, G. C., \Stability an Control of Stochastic Systems with Wie-Ban Noise Disturbances", Siam J. Appl. Math., Vol. 34, No. 3, pp , (978). [] orsthemke, W., an Lefever R., \A Perturbation Expansion for External Wie Ban Markovian Noise: Application to Transitions nuce by Ornstein-Uhlenbeck Noise",. Physik B, Vol. 4, pp , (98). [2] Klosek-Dygas, M. M., Matkowsky, B. J., an Schuss,., \Stochastic Stability of Nonlinear Oscillators", Siam J. Appl. Math., Vol. 48, No. 5, pp. 5-27, (988). [3] Klosek-Dygas, M. M., Matkowsky, B. J., an Schuss,., \Colore Noise in Dynamical Systems", Siam J. Appl. Math., 48, No. 2, pp , (988). [4] Roberts, J. B., \A Stochastic Theory for Nonlinear Ship Rolling in rregular Seas", J. Ship Research, Vol. 26, No. 4, pp , (982). [5] Roberts, J. B., \Energy Methos for Non-Linear Systems with Non-white Excitation", Proceeings of the UTAM Symposium on Ranom Vibrations an Reliability, (Frankfurt/Oer, G. D. R., 982), (K. ennig, E.), pp , Akaemie, Berlin, (983). [6] Gariner, C. W., anbook of Stochastic Methos, Springer-Verlag, Berlin, (983). [7] Erelyi, A., igher Transcenental Functions, Vol., McGraw-ill, New York, (953). [8] Atkinson, J. D., \Eigenfunction Expansions for Ranomly Excite Non-Linear Systems", Journal of Soun an Vibration, Vol. 3, No. 2, pp , (973). [9] Ariaratnam, S. T., an Tam, D. S. F., \Ranom Vibration an Stability of a Linear Parametrically Excite Oscillator", AMM, Vol. 59, pp , (979). 5
THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE
Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek
More informationTMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments
Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary
More information19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control
19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior
More informationAdvanced Partial Differential Equations with Applications
MIT OpenCourseWare http://ocw.mit.eu 18.306 Avance Partial Differential Equations with Applications Fall 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.eu/terms.
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More information05 The Continuum Limit and the Wave Equation
Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,
More informationLectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs
Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent
More informationθ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
More informationIntroduction to Markov Processes
Introuction to Markov Processes Connexions moule m44014 Zzis law Gustav) Meglicki, Jr Office of the VP for Information Technology Iniana University RCS: Section-2.tex,v 1.24 2012/12/21 18:03:08 gustav
More information1 Introuction In the past few years there has been renewe interest in the nerson impurity moel. This moel was originally propose by nerson [2], for a
Theory of the nerson impurity moel: The Schrieer{Wol transformation re{examine Stefan K. Kehrein 1 an nreas Mielke 2 Institut fur Theoretische Physik, uprecht{karls{universitat, D{69120 Heielberg, F..
More informationProblems Governed by PDE. Shlomo Ta'asan. Carnegie Mellon University. and. Abstract
Pseuo-Time Methos for Constraine Optimization Problems Governe by PDE Shlomo Ta'asan Carnegie Mellon University an Institute for Computer Applications in Science an Engineering Abstract In this paper we
More informationGLOBAL SOLUTIONS FOR 2D COUPLED BURGERS-COMPLEX-GINZBURG-LANDAU EQUATIONS
Electronic Journal of Differential Equations, Vol. 015 015), No. 99, pp. 1 14. ISSN: 107-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu ftp eje.math.txstate.eu GLOBAL SOLUTIONS FOR D COUPLED
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationRole of parameters in the stochastic dynamics of a stick-slip oscillator
Proceeing Series of the Brazilian Society of Applie an Computational Mathematics, v. 6, n. 1, 218. Trabalho apresentao no XXXVII CNMAC, S.J. os Campos - SP, 217. Proceeing Series of the Brazilian Society
More informationSINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES
Communications on Stochastic Analysis Vol. 2, No. 2 (28) 289-36 Serials Publications www.serialspublications.com SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES
More informationensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y
Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay
More informationand from it produce the action integral whose variation we set to zero:
Lagrange Multipliers Monay, 6 September 01 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine
More informationAccelerate Implementation of Forwaring Control Laws using Composition Methos Yves Moreau an Roolphe Sepulchre June 1997 Abstract We use a metho of int
Katholieke Universiteit Leuven Departement Elektrotechniek ESAT-SISTA/TR 1997-11 Accelerate Implementation of Forwaring Control Laws using Composition Methos 1 Yves Moreau, Roolphe Sepulchre, Joos Vanewalle
More informationG j dq i + G j. q i. = a jt. and
Lagrange Multipliers Wenesay, 8 September 011 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine
More informationCalculus of Variations
16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t
More informationConservation Laws. Chapter Conservation of Energy
20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action
More informationA LIMIT THEOREM FOR RANDOM FIELDS WITH A SINGULARITY IN THE SPECTRUM
Teor Imov r. ta Matem. Statist. Theor. Probability an Math. Statist. Vip. 81, 1 No. 81, 1, Pages 147 158 S 94-911)816- Article electronically publishe on January, 11 UDC 519.1 A LIMIT THEOREM FOR RANDOM
More informationAssignment 1. g i (x 1,..., x n ) dx i = 0. i=1
Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition
More informationSlide10 Haykin Chapter 14: Neurodynamics (3rd Ed. Chapter 13)
Slie10 Haykin Chapter 14: Neuroynamics (3r E. Chapter 13) CPSC 636-600 Instructor: Yoonsuck Choe Spring 2012 Neural Networks with Temporal Behavior Inclusion of feeback gives temporal characteristics to
More informationH. MOUNIER, J. RUDOLPH, M. FLIESS, AND P. ROUCHON. interior mass is considered. It is viewed as a linear delay system. A
ESAIM: Control, Optimisation an Calculus of Variations URL: http://www.emath.fr/cocv/ September 998, Vol. 3, 35{3 TRACKING CONTROL OF A VIBRATING STRING WITH AN INTERIOR MASS VIEWED AS DELAY SYSTEM H.
More informationOptimal Measurement and Control in Quantum Dynamical Systems.
Optimal Measurement an Control in Quantum Dynamical Systems. V P Belavin Institute of Physics, Copernicus University, Polan. (On leave of absence from MIEM, Moscow, USSR) Preprint No 411, Torun, February
More informationThermal conductivity of graded composites: Numerical simulations and an effective medium approximation
JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University
More informationVariational principle for limit cycles of the Rayleigh van der Pol equation
PHYICAL REVIEW E VOLUME 59, NUMBER 5 MAY 999 Variational principle for limit cycles of the Rayleigh van er Pol equation R. D. Benguria an M. C. Depassier Faculta e Física, Pontificia Universia Católica
More informationThe effect of dissipation on solutions of the complex KdV equation
Mathematics an Computers in Simulation 69 (25) 589 599 The effect of issipation on solutions of the complex KV equation Jiahong Wu a,, Juan-Ming Yuan a,b a Department of Mathematics, Oklahoma State University,
More informationSturm-Liouville Theory
LECTURE 5 Sturm-Liouville Theory In the three preceing lectures I emonstrate the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series are just the tip of the iceberg of the theory
More informationConvergence of Random Walks
Chapter 16 Convergence of Ranom Walks This lecture examines the convergence of ranom walks to the Wiener process. This is very important both physically an statistically, an illustrates the utility of
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
More information'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21
Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting
More informationIERCU. Institute of Economic Research, Chuo University 50th Anniversary Special Issues. Discussion Paper No.210
IERCU Institute of Economic Research, Chuo University 50th Anniversary Special Issues Discussion Paper No.210 Discrete an Continuous Dynamics in Nonlinear Monopolies Akio Matsumoto Chuo University Ferenc
More informationLie symmetry and Mei conservation law of continuum system
Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive
More informationConservation laws a simple application to the telegraph equation
J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness
More informationGeneralization of the persistent random walk to dimensions greater than 1
PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,
More informationEntanglement is not very useful for estimating multiple phases
PHYSICAL REVIEW A 70, 032310 (2004) Entanglement is not very useful for estimating multiple phases Manuel A. Ballester* Department of Mathematics, University of Utrecht, Box 80010, 3508 TA Utrecht, The
More informationThe Three-dimensional Schödinger Equation
The Three-imensional Schöinger Equation R. L. Herman November 7, 016 Schröinger Equation in Spherical Coorinates We seek to solve the Schröinger equation with spherical symmetry using the metho of separation
More informationHyperbolic Systems of Equations Posed on Erroneous Curved Domains
Hyperbolic Systems of Equations Pose on Erroneous Curve Domains Jan Norström a, Samira Nikkar b a Department of Mathematics, Computational Mathematics, Linköping University, SE-58 83 Linköping, Sween (
More informationSYMMETRIC KRONECKER PRODUCTS AND SEMICLASSICAL WAVE PACKETS
SYMMETRIC KRONECKER PRODUCTS AND SEMICLASSICAL WAVE PACKETS GEORGE A HAGEDORN AND CAROLINE LASSER Abstract We investigate the iterate Kronecker prouct of a square matrix with itself an prove an invariance
More informationPhysics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1
Physics 5153 Classical Mechanics The Virial Theorem an The Poisson Bracket 1 Introuction In this lecture we will consier two applications of the Hamiltonian. The first, the Virial Theorem, applies to systems
More informationPlacement and tuning of resonance dampers on footbridges
Downloae from orbit.tu.k on: Jan 17, 19 Placement an tuning of resonance ampers on footbriges Krenk, Steen; Brønen, Aners; Kristensen, Aners Publishe in: Footbrige 5 Publication ate: 5 Document Version
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationarxiv: v1 [math-ph] 5 May 2014
DIFFERENTIAL-ALGEBRAIC SOLUTIONS OF THE HEAT EQUATION VICTOR M. BUCHSTABER, ELENA YU. NETAY arxiv:1405.0926v1 [math-ph] 5 May 2014 Abstract. In this work we introuce the notion of ifferential-algebraic
More informationSparse Reconstruction of Systems of Ordinary Differential Equations
Sparse Reconstruction of Systems of Orinary Differential Equations Manuel Mai a, Mark D. Shattuck b,c, Corey S. O Hern c,a,,e, a Department of Physics, Yale University, New Haven, Connecticut 06520, USA
More informationRegularized extremal bounds analysis (REBA): An approach to quantifying uncertainty in nonlinear geophysical inverse problems
GEOPHYSICAL RESEARCH LETTERS, VOL. 36, L03304, oi:10.1029/2008gl036407, 2009 Regularize extremal bouns analysis (REBA): An approach to quantifying uncertainty in nonlinear geophysical inverse problems
More informationEuler equations for multiple integrals
Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................
More informationPartial Differential Equations
Chapter Partial Differential Equations. Introuction Have solve orinary ifferential equations, i.e. ones where there is one inepenent an one epenent variable. Only orinary ifferentiation is therefore involve.
More informationSwitching Time Optimization in Discretized Hybrid Dynamical Systems
Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set
More informationDiagonalization of Matrices Dr. E. Jacobs
Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is
More informationON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM
ON THE OPTIMALITY SYSTEM FOR A D EULER FLOW PROBLEM Eugene M. Cliff Matthias Heinkenschloss y Ajit R. Shenoy z Interisciplinary Center for Applie Mathematics Virginia Tech Blacksburg, Virginia 46 Abstract
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationStable and compact finite difference schemes
Center for Turbulence Research Annual Research Briefs 2006 2 Stable an compact finite ifference schemes By K. Mattsson, M. Svär AND M. Shoeybi. Motivation an objectives Compact secon erivatives have long
More informationarxiv: v2 [cond-mat.stat-mech] 11 Nov 2016
Noname manuscript No. (will be inserte by the eitor) Scaling properties of the number of ranom sequential asorption iterations neee to generate saturate ranom packing arxiv:607.06668v2 [con-mat.stat-mech]
More informationState observers and recursive filters in classical feedback control theory
State observers an recursive filters in classical feeback control theory State-feeback control example: secon-orer system Consier the riven secon-orer system q q q u x q x q x x x x Here u coul represent
More informationLecture 2 - First order linear PDEs and PDEs from physics
18.15 - Introuction to PEs, Fall 004 Prof. Gigliola Staffilani Lecture - First orer linear PEs an PEs from physics I mentione in the first class some basic PEs of first an secon orer. Toay we illustrate
More informationLeast-Squares Regression on Sparse Spaces
Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction
More informationLATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION
The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische
More informationI. INTRODUCTION In this work we analize the slow roll approximation in cosmological moels of nonminimally couple theories of gravity (NMC) starting fr
Slow Roll Ination in Non-Minimally Couple Theories: Hyperextene Gravity Approach Diego F. Torres Departamento e Fsica, Universia Nacional e La Plata C.C. 67, 900, La Plata, Buenos Aires, Argentina Abstract
More informationTopic 7: Convergence of Random Variables
Topic 7: Convergence of Ranom Variables Course 003, 2016 Page 0 The Inference Problem So far, our starting point has been a given probability space (S, F, P). We now look at how to generate information
More informationThe Principle of Least Action
Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of
More informationNOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,
NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which
More informationSystems & Control Letters
Systems & ontrol Letters ( ) ontents lists available at ScienceDirect Systems & ontrol Letters journal homepage: www.elsevier.com/locate/sysconle A converse to the eterministic separation principle Jochen
More informationRank, Trace, Determinant, Transpose an Inverse of a Matrix Let A be an n n square matrix: A = a11 a1 a1n a1 a an a n1 a n a nn nn where is the jth col
Review of Linear Algebra { E18 Hanout Vectors an Their Inner Proucts Let X an Y be two vectors: an Their inner prouct is ene as X =[x1; ;x n ] T Y =[y1; ;y n ] T (X; Y ) = X T Y = x k y k k=1 where T an
More informationChapter 6: Energy-Momentum Tensors
49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.
More informationHarmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method
1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse
More informationarxiv: v1 [hep-lat] 19 Nov 2013
HU-EP-13/69 SFB/CPP-13-98 DESY 13-225 Applicability of Quasi-Monte Carlo for lattice systems arxiv:1311.4726v1 [hep-lat] 19 ov 2013, a,b Tobias Hartung, c Karl Jansen, b Hernan Leovey, Anreas Griewank
More informationMany problems in physics, engineering, and chemistry fall in a general class of equations of the form. d dx. d dx
Math 53 Notes on turm-liouville equations Many problems in physics, engineering, an chemistry fall in a general class of equations of the form w(x)p(x) u ] + (q(x) λ) u = w(x) on an interval a, b], plus
More informationHyperbolic Moment Equations Using Quadrature-Based Projection Methods
Hyperbolic Moment Equations Using Quarature-Base Projection Methos J. Koellermeier an M. Torrilhon Department of Mathematics, RWTH Aachen University, Aachen, Germany Abstract. Kinetic equations like the
More informationA nonlinear inverse problem of the Korteweg-de Vries equation
Bull. Math. Sci. https://oi.org/0.007/s3373-08-025- A nonlinear inverse problem of the Korteweg-e Vries equation Shengqi Lu Miaochao Chen 2 Qilin Liu 3 Receive: 0 March 207 / Revise: 30 April 208 / Accepte:
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationBEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi
BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS Mauro Boccaoro Magnus Egerstet Paolo Valigi Yorai Wari {boccaoro,valigi}@iei.unipg.it Dipartimento i Ingegneria Elettronica
More informationLaplace s Equation in Cylindrical Coordinates and Bessel s Equation (II)
Laplace s Equation in Cylinrical Coorinates an Bessel s Equation (II Qualitative properties of Bessel functions of first an secon kin In the last lecture we foun the expression for the general solution
More informationarxiv: v1 [physics.flu-dyn] 8 May 2014
Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa
More informationANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS
ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS MICHAEL HOLST, EVELYN LUNASIN, AND GANTUMUR TSOGTGEREL ABSTRACT. We consier a general family of regularize Navier-Stokes an Magnetohyroynamics
More informationCOUPLING REQUIREMENTS FOR WELL POSED AND STABLE MULTI-PHYSICS PROBLEMS
VI International Conference on Computational Methos for Couple Problems in Science an Engineering COUPLED PROBLEMS 15 B. Schrefler, E. Oñate an M. Paparakakis(Es) COUPLING REQUIREMENTS FOR WELL POSED AND
More informationCalculus and optimization
Calculus an optimization These notes essentially correspon to mathematical appenix 2 in the text. 1 Functions of a single variable Now that we have e ne functions we turn our attention to calculus. A function
More informationthe solution of ()-(), an ecient numerical treatment requires variable steps. An alternative approach is to apply a time transformation of the form t
Asymptotic Error Analysis of the Aaptive Verlet Metho Stephane Cirilli, Ernst Hairer Beneict Leimkuhler y May 3, 999 Abstract The Aaptive Verlet metho [7] an variants [6] are time-reversible schemes for
More informationApplication of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate
Freun Publishing House Lt., International Journal of Nonlinear Sciences & Numerical Simulation, (9), -, 9 Application of the homotopy perturbation metho to a magneto-elastico-viscous flui along a semi-infinite
More informationarxiv:nlin/ v1 [nlin.cd] 21 Mar 2002
Entropy prouction of iffusion in spatially perioic eterministic systems arxiv:nlin/0203046v [nlin.cd] 2 Mar 2002 J. R. Dorfman, P. Gaspar, 2 an T. Gilbert 3 Department of Physics an Institute for Physical
More informationSemiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom
PHYSICAL REVIEW A 69, 063409 (2004) Semiclassical analysis of long-wavelength multiphoton processes: The Ryberg atom Luz V. Vela-Arevalo* an Ronal F. Fox Center for Nonlinear Sciences an School of Physics,
More informationMath 300 Winter 2011 Advanced Boundary Value Problems I. Bessel s Equation and Bessel Functions
Math 3 Winter 2 Avance Bounary Value Problems I Bessel s Equation an Bessel Functions Department of Mathematical an Statistical Sciences University of Alberta Bessel s Equation an Bessel Functions We use
More informationRobustness and Perturbations of Minimal Bases
Robustness an Perturbations of Minimal Bases Paul Van Dooren an Froilán M Dopico December 9, 2016 Abstract Polynomial minimal bases of rational vector subspaces are a classical concept that plays an important
More informationAbstract A nonlinear partial differential equation of the following form is considered:
M P E J Mathematical Physics Electronic Journal ISSN 86-6655 Volume 2, 26 Paper 5 Receive: May 3, 25, Revise: Sep, 26, Accepte: Oct 6, 26 Eitor: C.E. Wayne A Nonlinear Heat Equation with Temperature-Depenent
More informationParameter estimation: A new approach to weighting a priori information
Parameter estimation: A new approach to weighting a priori information J.L. Mea Department of Mathematics, Boise State University, Boise, ID 83725-555 E-mail: jmea@boisestate.eu Abstract. We propose a
More informationa) Identify the kinematical constraint relating motions Y and X. The cable does NOT slip on the pulley. For items (c) & (e-f-g) use
EAMPLE PROBLEM for MEEN 363 SPRING 6 Objectives: a) To erive EOMS of a DOF system b) To unerstan concept of static equilibrium c) To learn the correct usage of physical units (US system) ) To calculate
More informationF 1 x 1,,x n. F n x 1,,x n
Chapter Four Dynamical Systems 4. Introuction The previous chapter has been evote to the analysis of systems of first orer linear ODE s. In this chapter we will consier systems of first orer ODE s that
More informationSurvey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013
Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing
More informationSTATISTICAL LIKELIHOOD REPRESENTATIONS OF PRIOR KNOWLEDGE IN MACHINE LEARNING
STATISTICAL LIKELIHOOD REPRESENTATIONS OF PRIOR KNOWLEDGE IN MACHINE LEARNING Mark A. Kon Department of Mathematics an Statistics Boston University Boston, MA 02215 email: mkon@bu.eu Anrzej Przybyszewski
More informationJointly continuous distributions and the multivariate Normal
Jointly continuous istributions an the multivariate Normal Márton alázs an álint Tóth October 3, 04 This little write-up is part of important founations of probability that were left out of the unit Probability
More informationFLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS. 1. Introduction
FLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS ALINA BUCUR, CHANTAL DAVID, BROOKE FEIGON, MATILDE LALÍN 1 Introuction In this note, we stuy the fluctuations in the number
More informationMinimum-time constrained velocity planning
7th Meiterranean Conference on Control & Automation Makeonia Palace, Thessaloniki, Greece June 4-6, 9 Minimum-time constraine velocity planning Gabriele Lini, Luca Consolini, Aurelio Piazzi Università
More informationLeft-invariant extended Kalman filter and attitude estimation
Left-invariant extene Kalman filter an attitue estimation Silvere Bonnabel Abstract We consier a left-invariant ynamics on a Lie group. One way to efine riving an observation noises is to make them preserve
More informationWe G Model Reduction Approaches for Solution of Wave Equations for Multiple Frequencies
We G15 5 Moel Reuction Approaches for Solution of Wave Equations for Multiple Frequencies M.Y. Zaslavsky (Schlumberger-Doll Research Center), R.F. Remis* (Delft University) & V.L. Druskin (Schlumberger-Doll
More informationA Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential
Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix
More informationImplicit Differentiation
Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,
More informationINTRODUCTION TO SYMPLECTIC MECHANICS: LECTURE IV. Maurice de Gosson
INTRODUCTION TO SYMPLECTIC MECHANICS: LECTURE IV Maurice e Gosson ii 4 Hamiltonian Mechanics Physically speaking, Hamiltonian mechanics is a paraphrase (an generalization!) of Newton s secon law, popularly
More informationmodel considered before, but the prey obey logistic growth in the absence of predators. In
5.2. First Orer Systems of Differential Equations. Phase Portraits an Linearity. Section Objective(s): Moifie Preator-Prey Moel. Graphical Representations of Solutions. Phase Portraits. Vector Fiels an
More information