arxiv:quant-ph/ v1 5 Jul 2004
|
|
- Oliver Barker
- 6 years ago
- Views:
Transcription
1 Numerical Mehods for Sochasic Differenial Equaions Joshua Wilkie Deparmen of Chemisry, Simon Fraser Universiy, Burnaby, Briish Columbia V5A 1S6, Canada Sochasic differenial equaions (sdes) play an imporan role in physics bu exising numerical mehods for solving such equaions are of low accuracy and poor sabiliy. A general sraegy for developing accurae and efficien schemes for solving sochasic equaions in oulined here. High order numerical mehods are developed for inegraion of sochasic differenial equaions wih srong soluions. We demonsrae he accuracy of he resuling inegraion schemes by compuing he errors in approximae soluions for sdes which have known exac soluions. PACS numbers: w, 2.5.-r, 2.7.-c arxiv:quan-ph/4739v1 5 Jul 24 Sochasic differenial equaions (sdes) have a long hisory in physics[1] and play an imporan role in many oher areas of science, engineering and finance[1, 2, 3]. Recenly a number of compuaional echniques have been developed in which high dimensional deerminisic equaions are decomposed ino lower dimensional sochasic equaions. Gisin and Percival[4], for example, reduced a deerminisic maser equaion for he densiy marix ino sochasic equaions for a wavefuncion. Similar approaches are being used o solve he quanum many-body problem for bosons[5], fermions[6] and vibraions[7]. These laer mehods give rise o large ses of coupled sdes which require fas and efficien numerical inegraion schemes. Unforunaely, and in spie of heir widespread use, he available numerical echniques[3] for solving such equaions are far less accurae han comparable mehods for soluion of ordinary differenial equaions (odes). In his manuscrip we show how classical mehods for solving odes, such as Runge-Kua, can be adaped for he soluion of a class of sdes which should include many of he equaions which arise in physical problems. Consider a finie se of sdes, dx j = aj (X, ) d + b j k (X, ) dw k, (1) represened in Iô[1, 2, 3] form, where j = 1,...,n. Here X = (X 1,..., Xn ) and he dw k are independen and normally disribued sochasic differenials wih zero mean and variance d (i.e. sampled N(, d)). The sochasic variables W k are Wiener processes. Now assume ha he coefficiens a j and b j k have regulariy properies which guaranee srong soluions, i.e. ha X j are some fixed funcions of he Wiener processes, and ha hey are differeniable o high order. [Sufficien condiions for srong soluions are discussed in Ref. [3].] We may hen view he soluions of (1) as funcions X j = X j (, W 1,...,W m ) of ime and he Wiener processes. The soluions can herefore be expanded in Taylor series. Keeping erms of order d or less hen gives X j +d = X j + Xj d + X j W k dw k k,l=1 2 X j W k W l dw k dw l. (2) In a mean square sense he produc of differenials dw kdw l is equivalen o δ k,ld in he Iô[1, 2, 3] formulaion of sochasic calculus. Making his replacemen hen yields dx j +d = X j +d Xj = [ Xj X j + W k 2 X j ] d W k2 dw k (3) which when compared o (1) allows us o idenify he firs derivaives X j W k X j = b j k (X, ) (4) = a j (X, ) 1 2 = a j (X, ) X j W k2 n b i k(x, ) bj k (X, ) X i.(5) i=1 Now ha hese firs order derivaives are expressed in erms of a j and b j k, higher order derivaives can be compued. Thus a Taylor expansion of he soluions X j + = X j + Xj + m k,l=1 X j W k W k 2 X j W k W l W k W l +... (6) can be obained for finie displacemens and W k. This Taylor expansion can hen be employed o develop Runge-Kua algorihms and oher inegraion schemes. We illusrae he use of his approach by developing a Runge-Kua mehod for sdes which is closely relaed o he classical Runge-Kua scheme for odes. For given displacemens and W k define f j (X, ) = Xj + m X j W k W k
2 2 = [a j (X, ) n i=1 b i k (X, ) bj k (X, ) X i ] b j k (X, ) W k (7) and consider he following four sage approximaion K 1 j = f j (X i, i ) Figure 1: log 1 X X approximae vs ime for Eq. (9) K 2 j K 3 j = f j (X i K1, i ) = f j (X i K2, i ) K 4 j = f j (X i + K 3, i+1 ) X i+1 = X i (K1 + 2K 2 + 2K 3 + K 4 ) (8) where i is he iniial ime and i+1 = i +. Taylor expansion of his scheme shows ha X i+1 differs from he exac soluion by erms of order higher han 2 (i.e. erms of higher order han 2, ( W k)2, ( W k)4, ( W k ) 2 ( W) l 2, and ( W k ) 2 W W l i ). Thus, his sochasic Runge-Kua algorihm plays a role very similar o is classical counerpar excep ha is order is reduced from four o wo. Generalizaions o higher order Runge-Kua schemes are sraighforward, and we will employ one such scheme in example calculaions, bu deails will no be presened here. While his approach is no compleely general, since i will fail for sdes wih weak soluions or non-differeniable a j and b j k, i should be applicable o a wide range of problems. I can for example be used o solve every one of he equaions wih known soluions abulaed in secion 4.4 of Ref. [3]. To illusrae he accuracy of he mehod and is improvemen over oher known echniques for solving sdes we now consider a number of hese examples. We compare known exac soluions wih numerical soluions obained using he Euler-Maruyama scheme[8], a derivaive free version of he Milsein scheme due o Kloeden and Plaen[9], he classical Runge-Kua scheme (8), and anoher Runge-Kua scheme obained in he manner oulined above from an eighh order welve sep mehod for odes due o Hairer and Wanner[1] (his reproduces he sochasic Taylor expansion up o and including erms of order 4 ). Sochasic differenials were sampled using he rouines gasdev and ran2[11]. As a firs es of hese mehods consider an auonomous nonlinear scalar equaion dx = (1 + X )(1 + X 2 )d + (1 + X2 )dw (9) wih jus one Wiener process. In his example and in all subsequen examples we assume all Wiener processes are iniially zero. The exac soluion o his equaion is[3] X = an( + W + arcan(x )) (1) as can be readily verified using Iô[1, 2, 3] calculus. In Fig. 1 we plo he error log 1 X X approximae vs Figure 2: log 1 X X approximae vs ime for Eq. (11) ime compued wih a ime sep of for a single sochasic rajecory wih iniial condiion X = 1 for he four differen approximaion schemes. The Milsein scheme (long-dashed curve) shows some improvemen over he primiive Euler-Maruyama mehod (solid curve) bu he order wo Runge-Kua scheme (shordashed curve) and order four Runge-Kua scheme (doed curve) perform very much beer. The second example equaion, also from Ref. [3], is an auonomous linear scalar equaion in wo Wiener processes dx = a X d + b 1 X dw 1 + b 2 X dw 2 (11) which has an exac soluion X = X exp{[a 1 2 (b2 1 + b2 2 )] + b 1W 1 + b 2W 2 }. (12) The logarihm base en of he error for he differen schemes, calculaed for iniial condiion X = 1 and ime
3 3 Figure 3: log 1 X 1 X 1 approximae vs ime for Eq. (13) Figure 4: log 1 X X approximae vs ime for Eq. (15) sep.1, is ploed in Fig. 2. Here he Milsein scheme (long-dashed curve) performs no beer han he Euler- Maruyama mehod (solid curve) bu again he order wo Runge-Kua scheme (shor-dashed curve) and order four Runge-Kua scheme (doed curve) show grealy improved accuracy. [Noe ha he apparen improvemen in performance of all schemes a long ime is a resul of he fac ha he soluion decays o zero.] Example 3 is a se of wo coupled linear auonomous sdes Figure 5: log 1 X X approximae vs ime for Eq. (17) dx 1 = 3 2 X1 d + X1 dw 1 X1 dw 2 X2 dw 3 dx 2 = 3 2 X2 d + X2 dw 1 X2 dw 2 + X1 dw 3 (13) wih hree Wiener processes. Here he soluions are X 1 = exp{ 2 + W 1 W 2 } cosw 3 X 2 = exp{ 2 + W 1 W 2 } sinw3. (14) Numerical soluions were calculaed wih a ime sep of.1 and errors in X 1 are represened in Fig. 3. The order wo Runge-Kua scheme (long-dashed curve) and order four Runge-Kua scheme (shor-dashed curve) show improvemen over he Milsein scheme (solid curve). Similar resuls were obained for X 2. The examples we have considered so far have no had explicily ime dependen a j and b j k. Example 4 is a scalar non-auonomous sde 2 dx = [ 1 + X (1 + )2 ]d (1 + )2 dw (15) wih known soluion[3] ( ) X = X (1 + )2 (W + ). (16) Numerical soluions were calculaed using he order wo Runge-Kua scheme and a ime sep of.1, = and X = 1. The error is represened in Fig. 4. As in previous examples a high accuracy is achieved in spie of he rapid growh of he soluion. The comparaive smoohness of he error curve reflecs he fac he he deerminisic par of he soluion dominaes. We now consider an example for which an exac soluion is known bu which is expressed in erms a sochasic inegral. Consider he sochasic Ginzburg-Landau equaion dx = [ X 3 + (α σ2 )X ]d + σx dw (17) wih soluion[3] exp{α + σw } X = X. (18) 1 + 2X 2 exp{2αs + 2σW s}ds We chose α =.1, σ = 4, X = 1 and d = The sochasic inegral was compued using a Riemann sum
4 4 Figure 6: log 1 X X approximae vs ime for Eq. (19) wih exac soluion[3] ) X = arcsinh (e a sinhx + e a e as dw s. (2) Figure 7: log 1 n n approximae vs ime for Eq. (21) We se a =.2, b = 1, X = 1 and d = The sochasic inegral in he exac soluion was calculaed using he Iô[1, 2, 3] inegral formula wih he same ime sep. The error in he soluion calculaed wih he order wo Runge-Kua scheme is ploed in Fig. 6. As in all previous cases considered he accuracy is very good. Thus, he approach o solving sdes advocaed here works very well for he wide range of examples we have considered. The order 4 Runge-Kua mehod is clearly much more accurae han he order 2 Runge- Kua scheme. I also has an embedded lower order Runge-Kua scheme which can be employed o obain an error esimae suiable for sepsize conrol[1]. Hence is should be possible o use variable sepsizes o ensure he accuracy of he soluion. This sor of implemenaion is essenial for solving equaions which do no have known exac soluions. The only subley in developing such a mehod is ensuring ha he correc Wiener pah is mainained even when a sep mus be rejeced. This is achieved[12] by dividing he rejeced differenials d and dw k in wo segmens; d/2 and dw k /2 y followed by d/2 and dw k /2 + y where y is sampled N(, d/2). To illusrae he accuracy of he resuling variable sepsize algorihm we solve he Gisin-Percival[4] sochasic wave equaion for he nonlinear absorber (Eq. 4.2 of Ref. [4]) d ψ =.1(a a) ψ d + (2a 2 a 2 a 2 a 2 a 2 a 2 ) ψ d + 2(a 2 a 2 ) ψ dw (21) wih he same ime sep. Error in he soluion calculaed wih he order wo Runge-Kua scheme is ploed in Fig. 5. Good accuracy is again obained. Finally, we consider an example in which he exac soluion is expressed in erms of a Iô[1, 2, 3] sochasic inegral. Consider he sde dx = anhx (a b2 sech 2 X )d+bsechx dw (19) wih iniial sae ψ() >=. In Fig. 7 we plo he error in mean occupaion number n = M[ ψ a a ψ ] vs ime (Fig. 5 of Ref. [4]) where M[ ] denoes an average over sochasic realisaions. 1, 1, and 2 rajecories were used o calculae he solid curve, dashed curve and doed curve, respecively. Convergence o he exac resul is good. The auhor acknowledges he suppor of he Naural Sciences and Engineering Research Council of Canada. [1] C.W. Gardiner, Handbook of sochasic mehods, (Springer, Berlin, 1983). [2] H. Hasegawa, H. Ezawa, Suppl. Prog. Theor. Phys. 69, 11 (198). [3] P.E. Kloeden and E. Plaen, Numerical soluion of sochasic differenial equaions, (Springer, Berlin, 1995). [4] N. Gisin and I.C. Percival, J. Phys. A 25, 5677 (1992). [5] I. Carusoo and Y. Casin, Laser Physics 13, 59 (23); I. Carusoo, Y. Casin, and J. Dalibard, Phys. Rev. A 63, 2366 (21). [6] O. Juille, Ph. Chomaz, Phys. Rev. Le. 88, (22). [7] J. Wilkie, Phys. Rev. E 67, 1712 (23). [8] See for example secion 1 of chaper 9 in Ref. [3]. [9] See Eq. (1.7) of chaper 11 in Ref. [3]. [1] E. Hairer and G. Wanner, DOP853.f available a
5 5 hp://elib.zib.de/pub/elib/hairer-wanner/nonsiff/. [11] W.H. Press, S.A. Teukolsky, W.T. Veerling and B.P. Flannery, Numerical recipes, (Cambridge Universiy Press, Cambridge, 1992). [12] J.G. Gaines and T.J. Lyons, SIAM J. Appl. Mah. 57, 1455 (1997).
Chapter 2. First Order Scalar Equations
Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.
More informationSingle and Double Pendulum Models
Single and Double Pendulum Models Mah 596 Projec Summary Spring 2016 Jarod Har 1 Overview Differen ypes of pendulums are used o model many phenomena in various disciplines. In paricular, single and double
More informationFractional Method of Characteristics for Fractional Partial Differential Equations
Fracional Mehod of Characerisics for Fracional Parial Differenial Equaions Guo-cheng Wu* Modern Teile Insiue, Donghua Universiy, 188 Yan-an ilu Road, Shanghai 51, PR China Absrac The mehod of characerisics
More informationSimulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010
Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid
More informationPhysics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle
Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,
More informationOn a Fractional Stochastic Landau-Ginzburg Equation
Applied Mahemaical Sciences, Vol. 4, 1, no. 7, 317-35 On a Fracional Sochasic Landau-Ginzburg Equaion Nguyen Tien Dung Deparmen of Mahemaics, FPT Universiy 15B Pham Hung Sree, Hanoi, Vienam dungn@fp.edu.vn
More informationOrdinary differential equations. Phys 750 Lecture 7
Ordinary differenial equaions Phys 750 Lecure 7 Ordinary Differenial Equaions Mos physical laws are expressed as differenial equaions These come in hree flavours: iniial-value problems boundary-value problems
More informationAn Introduction to Malliavin calculus and its applications
An Inroducion o Malliavin calculus and is applicaions Lecure 5: Smoohness of he densiy and Hörmander s heorem David Nualar Deparmen of Mahemaics Kansas Universiy Universiy of Wyoming Summer School 214
More informationt is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...
Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger
More informationA variational radial basis function approximation for diffusion processes.
A variaional radial basis funcion approximaion for diffusion processes. Michail D. Vreas, Dan Cornford and Yuan Shen {vreasm, d.cornford, y.shen}@ason.ac.uk Ason Universiy, Birmingham, UK hp://www.ncrg.ason.ac.uk
More informationChapter 3 Boundary Value Problem
Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le
More informationEXERCISES FOR SECTION 1.5
1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler
More informationRecursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems
8 Froniers in Signal Processing, Vol. 1, No. 1, July 217 hps://dx.doi.org/1.2266/fsp.217.112 Recursive Leas-Squares Fixed-Inerval Smooher Using Covariance Informaion based on Innovaion Approach in Linear
More informationarxiv: v1 [math.ca] 15 Nov 2016
arxiv:6.599v [mah.ca] 5 Nov 26 Counerexamples on Jumarie s hree basic fracional calculus formulae for non-differeniable coninuous funcions Cheng-shi Liu Deparmen of Mahemaics Norheas Peroleum Universiy
More informationMATH 128A, SUMMER 2009, FINAL EXAM SOLUTION
MATH 28A, SUMME 2009, FINAL EXAM SOLUTION BENJAMIN JOHNSON () (8 poins) [Lagrange Inerpolaion] (a) (4 poins) Le f be a funcion defined a some real numbers x 0,..., x n. Give a defining equaion for he Lagrange
More informationVehicle Arrival Models : Headway
Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where
More informationODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004
ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform
More informationSTATE-SPACE MODELLING. A mass balance across the tank gives:
B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing
More informationThe expectation value of the field operator.
The expecaion value of he field operaor. Dan Solomon Universiy of Illinois Chicago, IL dsolom@uic.edu June, 04 Absrac. Much of he mahemaical developmen of quanum field heory has been in suppor of deermining
More informationAnnouncements: Warm-up Exercise:
Fri Apr 13 7.1 Sysems of differenial equaions - o model muli-componen sysems via comparmenal analysis hp//en.wikipedia.org/wiki/muli-comparmen_model Announcemens Warm-up Exercise Here's a relaively simple
More information6.2 Transforms of Derivatives and Integrals.
SEC. 6.2 Transforms of Derivaives and Inegrals. ODEs 2 3 33 39 23. Change of scale. If l( f ()) F(s) and c is any 33 45 APPLICATION OF s-shifting posiive consan, show ha l( f (c)) F(s>c)>c (Hin: In Probs.
More informationOrdinary Differential Equations
Ordinary Differenial Equaions 5. Examples of linear differenial equaions and heir applicaions We consider some examples of sysems of linear differenial equaions wih consan coefficiens y = a y +... + a
More informationSecond quantization and gauge invariance.
1 Second quanizaion and gauge invariance. Dan Solomon Rauland-Borg Corporaion Moun Prospec, IL Email: dsolom@uic.edu June, 1. Absrac. I is well known ha he single paricle Dirac equaion is gauge invarian.
More information(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)
Phase Plane Analysis of Linear Sysems Adaped from Applied Nonlinear Conrol by Sloine and Li The general form of a linear second-order sysem is a c b d From and b bc d a Differeniaion of and hen subsiuion
More information6. Stochastic calculus with jump processes
A) Trading sraegies (1/3) Marke wih d asses S = (S 1,, S d ) A rading sraegy can be modelled wih a vecor φ describing he quaniies invesed in each asse a each insan : φ = (φ 1,, φ d ) The value a of a porfolio
More informationMost Probable Phase Portraits of Stochastic Differential Equations and Its Numerical Simulation
Mos Probable Phase Porrais of Sochasic Differenial Equaions and Is Numerical Simulaion Bing Yang, Zhu Zeng and Ling Wang 3 School of Mahemaics and Saisics, Huazhong Universiy of Science and Technology,
More informationChapter 6. Systems of First Order Linear Differential Equations
Chaper 6 Sysems of Firs Order Linear Differenial Equaions We will only discuss firs order sysems However higher order sysems may be made ino firs order sysems by a rick shown below We will have a sligh
More informationVariational Iteration Method for Solving System of Fractional Order Ordinary Differential Equations
IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 1, Issue 6 Ver. II (Nov - Dec. 214), PP 48-54 Variaional Ieraion Mehod for Solving Sysem of Fracional Order Ordinary Differenial
More informationSection 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients
Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous
More informationA proof of Ito's formula using a di Title formula. Author(s) Fujita, Takahiko; Kawanishi, Yasuhi. Studia scientiarum mathematicarum H Citation
A proof of Io's formula using a di Tile formula Auhor(s) Fujia, Takahiko; Kawanishi, Yasuhi Sudia scieniarum mahemaicarum H Ciaion 15-134 Issue 8-3 Dae Type Journal Aricle Tex Version auhor URL hp://hdl.handle.ne/186/15878
More informationd 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3
and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or
More informationASTR415: Problem Set #5
ASTR45: Problem Se #5 Curran D. Muhlberger Universi of Marland (Daed: April 25, 27) Three ssems of coupled differenial equaions were sudied using inegraors based on Euler s mehod, a fourh-order Runge-Kua
More information10. State Space Methods
. Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he
More informationSome Basic Information about M-S-D Systems
Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,
More informationOrdinary Differential Equations
Lecure 22 Ordinary Differenial Equaions Course Coordinaor: Dr. Suresh A. Karha, Associae Professor, Deparmen of Civil Engineering, IIT Guwahai. In naure, mos of he phenomena ha can be mahemaically described
More informationSolutions to Assignment 1
MA 2326 Differenial Equaions Insrucor: Peronela Radu Friday, February 8, 203 Soluions o Assignmen. Find he general soluions of he following ODEs: (a) 2 x = an x Soluion: I is a separable equaion as we
More informationEfficient Solution of Fractional Initial Value Problems Using Expanding Perturbation Approach
Journal of mahemaics and compuer Science 8 (214) 359-366 Efficien Soluion of Fracional Iniial Value Problems Using Expanding Perurbaion Approach Khosro Sayevand Deparmen of Mahemaics, Faculy of Science,
More information23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes
Half-Range Series 2.5 Inroducion In his Secion we address he following problem: Can we find a Fourier series expansion of a funcion defined over a finie inerval? Of course we recognise ha such a funcion
More informationEXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE
Version April 30, 2004.Submied o CTU Repors. EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Per Krysl Universiy of California, San Diego La Jolla, California 92093-0085,
More informationUnsteady Flow Problems
School of Mechanical Aerospace and Civil Engineering Unseady Flow Problems T. J. Craf George Begg Building, C41 TPFE MSc CFD-1 Reading: J. Ferziger, M. Peric, Compuaional Mehods for Fluid Dynamics H.K.
More informationMath Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.
Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need
More informationImproved Approximate Solutions for Nonlinear Evolutions Equations in Mathematical Physics Using the Reduced Differential Transform Method
Journal of Applied Mahemaics & Bioinformaics, vol., no., 01, 1-14 ISSN: 179-660 (prin), 179-699 (online) Scienpress Ld, 01 Improved Approimae Soluions for Nonlinear Evoluions Equaions in Mahemaical Physics
More information2.160 System Identification, Estimation, and Learning. Lecture Notes No. 8. March 6, 2006
2.160 Sysem Idenificaion, Esimaion, and Learning Lecure Noes No. 8 March 6, 2006 4.9 Eended Kalman Filer In many pracical problems, he process dynamics are nonlinear. w Process Dynamics v y u Model (Linearized)
More informationTitle: On mean-square-stability of two-step Maruyama methods for stochastic delay differential equations
Elsevier Ediorial Sysem(m) for Journal of Compuaional and Applied Mahemaics Manuscrip Draf Manuscrip Number: Tile: On mean-square-sabiliy of wo-sep Maruyama mehods for sochasic delay differenial equaions
More informationNavneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi
Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec
More informationOptimal Path Planning for Flexible Redundant Robot Manipulators
25 WSEAS In. Conf. on DYNAMICAL SYSEMS and CONROL, Venice, Ialy, November 2-4, 25 (pp363-368) Opimal Pah Planning for Flexible Redundan Robo Manipulaors H. HOMAEI, M. KESHMIRI Deparmen of Mechanical Engineering
More informationAn Iterative Method for Solving Two Special Cases of Nonlinear PDEs
Conemporary Engineering Sciences, Vol. 10, 2017, no. 11, 55-553 HIKARI Ld, www.m-hikari.com hps://doi.org/10.12988/ces.2017.7651 An Ieraive Mehod for Solving Two Special Cases of Nonlinear PDEs Carlos
More informationChapter 4. Truncation Errors
Chaper 4. Truncaion Errors and he Taylor Series Truncaion Errors and he Taylor Series Non-elemenary funcions such as rigonomeric, eponenial, and ohers are epressed in an approimae fashion using Taylor
More informationWhat Ties Return Volatilities to Price Valuations and Fundamentals? On-Line Appendix
Wha Ties Reurn Volailiies o Price Valuaions and Fundamenals? On-Line Appendix Alexander David Haskayne School of Business, Universiy of Calgary Piero Veronesi Universiy of Chicago Booh School of Business,
More information23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes
Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals
More informationAn Introduction to Stochastic Programming: The Recourse Problem
An Inroducion o Sochasic Programming: he Recourse Problem George Danzig and Phil Wolfe Ellis Johnson, Roger Wes, Dick Cole, and Me John Birge Where o look in he ex pp. 6-7, Secion.2.: Inroducion o sochasic
More informationInternational Journal of Scientific & Engineering Research, Volume 4, Issue 10, October ISSN
Inernaional Journal of Scienific & Engineering Research, Volume 4, Issue 10, Ocober-2013 900 FUZZY MEAN RESIDUAL LIFE ORDERING OF FUZZY RANDOM VARIABLES J. EARNEST LAZARUS PIRIYAKUMAR 1, A. YAMUNA 2 1.
More informationLinear Response Theory: The connection between QFT and experiments
Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and
More informationFinal Spring 2007
.615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o
More informationConcourse Math Spring 2012 Worked Examples: Matrix Methods for Solving Systems of 1st Order Linear Differential Equations
Concourse Mah 80 Spring 0 Worked Examples: Marix Mehods for Solving Sysems of s Order Linear Differenial Equaions The Main Idea: Given a sysem of s order linear differenial equaions d x d Ax wih iniial
More informationMatrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality
Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]
More informationSolutions from Chapter 9.1 and 9.2
Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is
More informationSimulation of BSDEs and. Wiener Chaos Expansions
Simulaion of BSDEs and Wiener Chaos Expansions Philippe Briand Céline Labar LAMA UMR 5127, Universié de Savoie, France hp://www.lama.univ-savoie.fr/ Workshop on BSDEs Rennes, May 22-24, 213 Inroducion
More informationHomework 10 (Stats 620, Winter 2017) Due Tuesday April 18, in class Questions are derived from problems in Stochastic Processes by S. Ross.
Homework (Sas 6, Winer 7 Due Tuesday April 8, in class Quesions are derived from problems in Sochasic Processes by S. Ross.. A sochasic process {X(, } is said o be saionary if X(,..., X( n has he same
More informationLECTURE 1: GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS
LECTURE : GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS We will work wih a coninuous ime reversible Markov chain X on a finie conneced sae space, wih generaor Lf(x = y q x,yf(y. (Recall ha q
More informationA Note on the Equivalence of Fractional Relaxation Equations to Differential Equations with Varying Coefficients
mahemaics Aricle A Noe on he Equivalence of Fracional Relaxaion Equaions o Differenial Equaions wih Varying Coefficiens Francesco Mainardi Deparmen of Physics and Asronomy, Universiy of Bologna, and he
More informationGeorey E. Hinton. University oftoronto. Technical Report CRG-TR February 22, Abstract
Parameer Esimaion for Linear Dynamical Sysems Zoubin Ghahramani Georey E. Hinon Deparmen of Compuer Science Universiy oftorono 6 King's College Road Torono, Canada M5S A4 Email: zoubin@cs.orono.edu Technical
More informationPredator - Prey Model Trajectories and the nonlinear conservation law
Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories
More informationCERTAIN CLASSES OF SOLUTIONS OF LAGERSTROM EQUATIONS
SARAJEVO JOURNAL OF MATHEMATICS Vol.10 (22 (2014, 67 76 DOI: 10.5644/SJM.10.1.09 CERTAIN CLASSES OF SOLUTIONS OF LAGERSTROM EQUATIONS ALMA OMERSPAHIĆ AND VAHIDIN HADŽIABDIĆ Absrac. This paper presens sufficien
More informationApplication of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing
Applicaion of a Sochasic-Fuzzy Approach o Modeling Opimal Discree Time Dynamical Sysems by Using Large Scale Daa Processing AA WALASZE-BABISZEWSA Deparmen of Compuer Engineering Opole Universiy of Technology
More informationLecture 20: Riccati Equations and Least Squares Feedback Control
34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he
More information14 Autoregressive Moving Average Models
14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class
More informationPENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD
PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD HAN XIAO 1. Penalized Leas Squares Lasso solves he following opimizaion problem, ˆβ lasso = arg max β R p+1 1 N y i β 0 N x ij β j β j (1.1) for some 0.
More informationChapter 8 The Complete Response of RL and RC Circuits
Chaper 8 The Complee Response of RL and RC Circuis Seoul Naional Universiy Deparmen of Elecrical and Compuer Engineering Wha is Firs Order Circuis? Circuis ha conain only one inducor or only one capacior
More informationBook Corrections for Optimal Estimation of Dynamic Systems, 2 nd Edition
Boo Correcions for Opimal Esimaion of Dynamic Sysems, nd Ediion John L. Crassidis and John L. Junins November 17, 017 Chaper 1 This documen provides correcions for he boo: Crassidis, J.L., and Junins,
More informationThe following report makes use of the process from Chapter 2 in Dr. Cumming s thesis.
Zaleski 1 Joseph Zaleski Mah 451H Final Repor Conformal Mapping Mehods and ZST Hele Shaw Flow Inroducion The Hele Shaw problem has been sudied using linear sabiliy analysis and numerical mehods, bu a novel
More informationMorning Time: 1 hour 30 minutes Additional materials (enclosed):
ADVANCED GCE 78/0 MATHEMATICS (MEI) Differenial Equaions THURSDAY JANUARY 008 Morning Time: hour 30 minues Addiional maerials (enclosed): None Addiional maerials (required): Answer Bookle (8 pages) Graph
More informationApplication of He s Variational Iteration Method for Solving Seventh Order Sawada-Kotera Equations
Applied Mahemaical Sciences, Vol. 2, 28, no. 1, 471-477 Applicaion of He s Variaional Ieraion Mehod for Solving Sevenh Order Sawada-Koera Equaions Hossein Jafari a,1, Allahbakhsh Yazdani a, Javad Vahidi
More informationADDITIONAL PROBLEMS (a) Find the Fourier transform of the half-cosine pulse shown in Fig. 2.40(a). Additional Problems 91
ddiional Problems 9 n inverse relaionship exiss beween he ime-domain and freuency-domain descripions of a signal. Whenever an operaion is performed on he waveform of a signal in he ime domain, a corresponding
More informationMaple Tools for Differential Equations A. J. Meir
Maple Tools for Differenial Equaions A. J. Meir Copyrigh (C) A. J. Meir. All righs reserved. This workshee is for educaional use only. No par of his publicaion may be reproduced or ransmied for profi in
More information) were both constant and we brought them from under the integral.
YIELD-PER-RECRUIT (coninued The yield-per-recrui model applies o a cohor, bu we saw in he Age Disribuions lecure ha he properies of a cohor do no apply in general o a collecion of cohors, which is wha
More informationMath 334 Fall 2011 Homework 11 Solutions
Dec. 2, 2 Mah 334 Fall 2 Homework Soluions Basic Problem. Transform he following iniial value problem ino an iniial value problem for a sysem: u + p()u + q() u g(), u() u, u () v. () Soluion. Le v u. Then
More informationR.#W.#Erickson# Department#of#Electrical,#Computer,#and#Energy#Engineering# University#of#Colorado,#Boulder#
.#W.#Erickson# Deparmen#of#Elecrical,#Compuer,#and#Energy#Engineering# Universiy#of#Colorado,#Boulder# Chaper 2 Principles of Seady-Sae Converer Analysis 2.1. Inroducion 2.2. Inducor vol-second balance,
More informationCHAPTER 6: FIRST-ORDER CIRCUITS
EEE5: CI CUI T THEOY CHAPTE 6: FIST-ODE CICUITS 6. Inroducion This chaper considers L and C circuis. Applying he Kirshoff s law o C and L circuis produces differenial equaions. The differenial equaions
More informationSome Ramsey results for the n-cube
Some Ramsey resuls for he n-cube Ron Graham Universiy of California, San Diego Jozsef Solymosi Universiy of Briish Columbia, Vancouver, Canada Absrac In his noe we esablish a Ramsey-ype resul for cerain
More informationSolving a System of Nonlinear Functional Equations Using Revised New Iterative Method
Solving a Sysem of Nonlinear Funcional Equaions Using Revised New Ieraive Mehod Sachin Bhalekar and Varsha Dafardar-Gejji Absrac In he presen paper, we presen a modificaion of he New Ieraive Mehod (NIM
More informationTesting the Random Walk Model. i.i.d. ( ) r
he random walk heory saes: esing he Random Walk Model µ ε () np = + np + Momen Condiions where where ε ~ i.i.d he idea here is o es direcly he resricions imposed by momen condiions. lnp lnp µ ( lnp lnp
More informationStochastic Model for Cancer Cell Growth through Single Forward Mutation
Journal of Modern Applied Saisical Mehods Volume 16 Issue 1 Aricle 31 5-1-2017 Sochasic Model for Cancer Cell Growh hrough Single Forward Muaion Jayabharahiraj Jayabalan Pondicherry Universiy, jayabharahi8@gmail.com
More informationCSE 3802 / ECE Numerical Methods in Scientific Computation. Jinbo Bi. Department of Computer Science & Engineering
CSE 3802 / ECE 3431 Numerical Mehods in Scienific Compuaion Jinbo Bi Deparmen of Compuer Science & Engineering hp://www.engr.uconn.edu/~jinbo 1 Ph.D in Mahemaics The Insrucor Previous professional experience:
More informationLecture 1 Overview. course mechanics. outline & topics. what is a linear dynamical system? why study linear systems? some examples
EE263 Auumn 27-8 Sephen Boyd Lecure 1 Overview course mechanics ouline & opics wha is a linear dynamical sysem? why sudy linear sysems? some examples 1 1 Course mechanics all class info, lecures, homeworks,
More informationPhysics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution
Physics 7b: Saisical Mechanics Fokker-Planck Equaion The Langevin equaion approach o he evoluion of he velociy disribuion for he Brownian paricle migh leave you uncomforable. A more formal reamen of his
More informationCHAPTER 2 Signals And Spectra
CHAPER Signals And Specra Properies of Signals and Noise In communicaion sysems he received waveform is usually caegorized ino he desired par conaining he informaion, and he undesired par. he desired par
More informationModal identification of structures from roving input data by means of maximum likelihood estimation of the state space model
Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix
More informationDifferential Equations
Mah 21 (Fall 29) Differenial Equaions Soluion #3 1. Find he paricular soluion of he following differenial equaion by variaion of parameer (a) y + y = csc (b) 2 y + y y = ln, > Soluion: (a) The corresponding
More informationarxiv: v1 [math.na] 23 Feb 2016
EPJ Web of Conferences will be se by he publisher DOI: will be se by he publisher c Owned by he auhors, published by EDP Sciences, 16 arxiv:163.67v1 [mah.na] 3 Feb 16 Numerical Soluion of a Nonlinear Inegro-Differenial
More informationT L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB
Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal
More informationHaar Wavelet Operational Matrix Method for Solving Fractional Partial Differential Equations
Copyrigh 22 Tech Science Press CMES, vol.88, no.3, pp.229-243, 22 Haar Wavele Operaional Mari Mehod for Solving Fracional Parial Differenial Equaions Mingu Yi and Yiming Chen Absrac: In his paper, Haar
More informationElementary Differential Equations and Boundary Value Problems
Elemenar Differenial Equaions and Boundar Value Problems Boce. & DiPrima 9 h Ediion Chaper 1: Inroducion 1006003 คณ ตศาสตร ว ศวกรรม 3 สาขาว ชาว ศวกรรมคอมพ วเตอร ป การศ กษา 1/2555 ผศ.ดร.อร ญญา ผศ.ดร.สมศ
More informationHamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t
M ah 5 2 7 Fall 2 0 0 9 L ecure 1 0 O c. 7, 2 0 0 9 Hamilon- J acobi Equaion: Explici Formulas In his lecure we ry o apply he mehod of characerisics o he Hamilon-Jacobi equaion: u + H D u, x = 0 in R n
More informationGMM - Generalized Method of Moments
GMM - Generalized Mehod of Momens Conens GMM esimaion, shor inroducion 2 GMM inuiion: Maching momens 2 3 General overview of GMM esimaion. 3 3. Weighing marix...........................................
More informationMA 214 Calculus IV (Spring 2016) Section 2. Homework Assignment 1 Solutions
MA 14 Calculus IV (Spring 016) Secion Homework Assignmen 1 Soluions 1 Boyce and DiPrima, p 40, Problem 10 (c) Soluion: In sandard form he given firs-order linear ODE is: An inegraing facor is given by
More informationNumerical Solution of Fractional Variational Problems Using Direct Haar Wavelet Method
ISSN: 39-8753 Engineering and echnology (An ISO 397: 7 Cerified Organizaion) Vol. 3, Issue 5, May 4 Numerical Soluion of Fracional Variaional Problems Using Direc Haar Wavele Mehod Osama H. M., Fadhel
More informationDEPARTMENT OF ECONOMICS AND FINANCE COLLEGE OF BUSINESS AND ECONOMICS UNIVERSITY OF CANTERBURY CHRISTCHURCH, NEW ZEALAND
DEPARTMENT OF ECONOMICS AND FINANCE COLLEGE OF BUSINESS AND ECONOMICS UNIVERSITY OF CANTERBURY CHRISTCHURCH, NEW ZEALAND Asymmery and Leverage in Condiional Volailiy Models Michael McAleer WORKING PAPER
More informationEcon107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)
I. Definiions and Problems A. Perfec Mulicollineariy Econ7 Applied Economerics Topic 7: Mulicollineariy (Sudenmund, Chaper 8) Definiion: Perfec mulicollineariy exiss in a following K-variable regression
More informationGeneralized Chebyshev polynomials
Generalized Chebyshev polynomials Clemene Cesarano Faculy of Engineering, Inernaional Telemaic Universiy UNINETTUNO Corso Viorio Emanuele II, 39 86 Roma, Ialy email: c.cesarano@unineunouniversiy.ne ABSTRACT
More information