Solving a Chemical Reaction System by a PDE
|
|
- Esther Sharp
- 5 years ago
- Views:
Transcription
1 Solving a Chemical Reaction System by a PDE François Lemaire, Michel Petitot, Adrien Poteaux, Rosanne Ushi Robira, Samuel Vidal,... University Lille 1 (Symbolic Computation Team) FELIM 2012 Functional Equations in LIMoges 2012 Faculty of Sciences and Techniques of the University of Limoges 6 Mars 2012 Work supported by the ANR LEDA F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
2 Plan 1 Background/Motivation 2 Complete treatment of the degradation 3 General framework 4 Stationary distribution for different examples 5 Conclusion F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
3 Plan 1 Background/Motivation 2 Complete treatment of the degradation 3 General framework 4 Stationary distribution for different examples 5 Conclusion F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
4 Prospective work Initiated by Michel Petitot around Keywords modelization based on chemical reactions systems stochastic simulations statistics, probabilities Application modelization in Biology exact and approximate computation of probabilities and statistics F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
5 Presentation Modelization in Biology one possible approach : chemical reactions system Reaction : list of products list of reactants Example: G + P α β H (1) G k G + P (2) P d 0 (3) Species are G: free gene, H: gene with protein bound, P: protein Difficulties in biology, there are not many molecules reactions do not arise deterministically reactions can arise in (almost) any order randomness is important for understanding cells in biology F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
6 A simple example A k 0 similar to a nuclear desintegration A A Figure: Two different simulations with A(0) = 10 and k=1 Note: the simulations are obtained following the Gillespie algorithm, which repeats the following steps: compute when the next reaction will occur compute which reaction will occur update the numbers of molecules according the reaction occurring F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
7 A simple example A k 0 similar to a nuclear desintegration A A A A A A A Figure: Ten superposed simulations with A(0) = 10 et k=1 F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
8 A simple example A k 0 similar to a nuclear desintegration A Figure: Experimental mean value of 1000 simulations F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
9 Probabilities, statistics Notation/properties one denotes X(t) the (integer) number of molecules of X at time t X(t) is a random variable (i.e. at each time t, X(t) is a random variable) Denote π i (t) = Prob(X(t) = i), i.e. probability that there are i molecules of X at time t. Fundamental questions compute π i (t) compute the mean value, variance, moments... of X(t). avoid Gillespie algorithm (time consuming, and missing rare events) try to get exact or approximate formulaes, or numerical solutions if not possible Models of Stochastic Gene Expression and Weyl Algebra, Vidal et al., ANB2010. F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
10 Plan 1 Background/Motivation 2 Complete treatment of the degradation 3 General framework 4 Stationary distribution for different examples 5 Conclusion F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
11 Complete solution of the degradation A k 0 Initial condition: A(0) = 10 similar to a nuclear desintegration Formal series and the PDE Let us introduce Φ(t, z) = π i (t)z i the chemical master equation (demonstrations skipped) yields : Φ t i=0 = k(1 z) Φ z Maple pdsolve (!) : Φ(t, z) = F ((z 1)e kt ) where F is arbitrary if A(0) = 10, then π 10 (0) = 1, so Φ(0, z) = z 10 easy computation: F(z) = (1 + (z 1)e kt ) 10. One can compute π i (t), and also mean value, variance, moments,... F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
12 Complete solution of the degradation (2) A k 0 similar to a nuclear desintegration Φ(t, z) = π i (t)z i i=0 Mean value (i.e. first moment) by definition: E(X(t)) = i=0 iπ i(t) Φ = z i=1 iπ i(t)z i 1 one obtains the mean value by setting z = 1 E(X(t)) = Φ (t, 1) =... = 10e kt z Second moment similar computations yields (using a second derivative) the second moment where E 2 (X(t)) = i=0 i 2 π i (t) E 2 (X(t)) = 10e t (1 e t ) F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
13 Complete solution of the degradation (2) A k 0 similar to a nuclear desintegration Φ(t, z) = π i (t)z i i= t mean second moment Figure: Mean and second moment, for k = 1 F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
14 Complete solution of the degradation (3) Transformation into system of ODEs: E(X(t)) = ke(x(t)) E 2 (X(t)) = 2kE 2 (X(t)) + ke(t) Infinite If a reaction involves two or more reactants (always the case except in toy examples), an infinite cascad occurs: E(X(t)) = depends on E(X(t)), E 2 (X(t)) E 2 (X(t)) = depends on E(X(t)), E 2 (X(t)), E 3 (X(t))... One cannot truncate the system easily F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
15 Plan 1 Background/Motivation 2 Complete treatment of the degradation 3 General framework 4 Stationary distribution for different examples 5 Conclusion F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
16 General framework Ideal case introduce an indeterminate z i for each species denote Φ(t, z) = π α(t)z α where z = (z 1,..., z n) and α = (α 1,..., α n) transform the system of reactions into Φ t = HΦ where H is a differential operator in z with polynomial coefficients in z (similarities with an Hamiltonian) solve the PDE! How to solve the PDE? symbolically (as in the degradation example) numerically convert the PDE into a ODE system, where unknown functions are mean values, variances,... F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
17 Difficulties Symbolically for exact solving, special functions arise even for easy examples for approximate solving, truncations are necessary where to truncate? Numerically except in toy systems, the operator H always involves second order derivatives the initial condition is given at t = 0 in the form Φ(0, z) = z α 0 (α 0 is the vector of number of species at t = 0) method of characteristics and usual PDE tricks fail Conversion to an ODE system except in toy systems, the ODE system is infinite, and cannot be truncated exactly (infinite cascad) F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
18 A subproblem Stationary distribution when t, the distribution might stabilize, if so, this distribution is time independant, and satisfies HΦ(z) = 0 (no more t) it gives information on the behaviour at t = Degradation example Φ t one looks fo k(1 z) Φ(z) z = 0 yields, Φ(z) = c, where c is a constant = HΦ where H = k(1 z) z since Φ(z) is a distribution, Φ(1) = 1, so Φ(z) = 1 this means that the system stabilizes with π 0 = 1 (no molecules) F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
19 Plan 1 Background/Motivation 2 Complete treatment of the degradation 3 General framework 4 Stationary distribution for different examples 5 Conclusion F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
20 Double degradation with creation 2A µ 0 0 λ A Φ t = HΦ where H = 1 2 µ(1 z2 ) 2 + λ(z 1) z2 For the moment, we have no solution for Φ(t, z), but we could solve the stationary distribution Stationary distribution equations HΦ(z) = 0 where Φ(z) = π i z i, π i 0, Φ(1) = i=0 π i = 1 F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
21 A dimension problem? After simplication (1 + z)φ (z) = 2a Φ(z) where a = λ/µ Φ has nonnegative coefficients, whose sum is 1 Uniqueness of the expected solution the physics tells us the expected solution is unique the ODE yield two independent solutions as formal power series around the origin Looking for the right solution the expected f.p.series solution converges in 1, so in -1 also the expected solution goes "through" the singularity in -1 F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
22 Two solvings Hypergeometric (Petitot) ( ) Φ(z) = K 1 (1 + z) 0 F 1 2 2a(1 + z) where ( ) K = 2 0 F 1 2 4a Expansion at z=-1 By a change of variable y = 1 + z, and φ(y) = b i y i, and using gfun, b 0 = 0 (4) b 1 arbitrary (5) b i+1 = 2a b i for i 1 (6) i(i + 1) Remark: Φ has an infinite radius of convergence, and is extremely convergent (I believe b i = o((i!) 2+ɛ ). F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
23 A generic example p times {}}{ A + A + + A µ 0 0 λ A Φ t Stationary distribution equations = HΦ where H = 1 p! µ(1 zp ) p + λ(z 1) zp (1 + z + + z p 1 )Φ (p) (z) = p!a Φ(z) where a = λ/µ Φ has nonnegative coefficients, whose sum is 1 finding the right solution p-dimensional space of solutions for the PDE the wanted solution goes through the p 1 singularities (p 1 roots of the unity) F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
24 Thoughts / conjectures Stationary distribution equations (1 + z + + z p 1 )Φ (p) (z) = p!a Φ(z) where a = λ/µ Φ has nonnegative coefficients, whose sum is 1 Difficulties "easy" to solve at each singularity, hard to deal with several singularities at the same time a Picard approach does not always converge (by solving the BVP satisfying Φ(w) = 0 for each singularity w, and Φ(1) = 1) Conjectures there is only one solution with an infinite radius of convergence (and the general term is o((i!) p+ɛ ) the expected solution is that one! F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
25 Plan 1 Background/Motivation 2 Complete treatment of the degradation 3 General framework 4 Stationary distribution for different examples 5 Conclusion F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
26 Towards the real examples real examples are in several variables one needs to treat the complete system (not just the stationary distribution) G + P α β H (7) G k G + P (8) P d 0 (9) Φ t = α(h GP) 2 Φ Φ Φ Φ + β(gp H) + k(gp G) + d(1 P) G H H G P where Φ = Φ(t, G, H, P) Remark: the term H GP is the "singularity" for the stationary distribution. F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
27 Software See for a prototype in Maple New version in preparation (Michel Petitot, Kyrill Batmanov, Christian Versari,... ) enough work for a PhD student (proposal made this year in the Computer Algebra team) F. Lemaire (Lille 1) Chemical Reactions and PDE FELIM / 24
Lognormal Moment Closures for Biochemical Reactions
Lognormal Moment Closures for Biochemical Reactions Abhyudai Singh and João Pedro Hespanha Abstract In the stochastic formulation of chemical reactions, the dynamics of the the first M -order moments of
More informationGillespie s Algorithm and its Approximations. Des Higham Department of Mathematics and Statistics University of Strathclyde
Gillespie s Algorithm and its Approximations Des Higham Department of Mathematics and Statistics University of Strathclyde djh@maths.strath.ac.uk The Three Lectures 1 Gillespie s algorithm and its relation
More informationFast Probability Generating Function Method for Stochastic Chemical Reaction Networks
MATCH Communications in Mathematical and in Computer Chemistry MATCH Commun. Math. Comput. Chem. 71 (2014) 57-69 ISSN 0340-6253 Fast Probability Generating Function Method for Stochastic Chemical Reaction
More informationLECTURE 10: REVIEW OF POWER SERIES. 1. Motivation
LECTURE 10: REVIEW OF POWER SERIES By definition, a power series centered at x 0 is a series of the form where a 0, a 1,... and x 0 are constants. For convenience, we shall mostly be concerned with the
More informationComputing the Monodromy Group of a Plane Algebraic Curve Using a New Numerical-modular Newton-Puiseux Algorithm
Computing the Monodromy Group of a Plane Algebraic Curve Using a New Numerical-modular Newton-Puiseux Algorithm Poteaux Adrien XLIM-DMI UMR CNRS 6172 Université de Limoges, France SNC'07 University of
More informationComplex Analysis Slide 9: Power Series
Complex Analysis Slide 9: Power Series MA201 Mathematics III Department of Mathematics IIT Guwahati August 2015 Complex Analysis Slide 9: Power Series 1 / 37 Learning Outcome of this Lecture We learn Sequence
More informationMA22S3 Summary Sheet: Ordinary Differential Equations
MA22S3 Summary Sheet: Ordinary Differential Equations December 14, 2017 Kreyszig s textbook is a suitable guide for this part of the module. Contents 1 Terminology 1 2 First order separable 2 2.1 Separable
More informationter. on Can we get a still better result? Yes, by making the rectangles still smaller. As we make the rectangles smaller and smaller, the
Area and Tangent Problem Calculus is motivated by two main problems. The first is the area problem. It is a well known result that the area of a rectangle with length l and width w is given by A = wl.
More informationMoment Closure Techniques for Stochastic Models in Population Biology
Moment Closure Techniques for Stochastic Models in Population Biology Abhyudai Singh and João Pedro Hespanha Abstract Continuous-time birth-death Markov processes serve as useful models in population biology.
More informationProblem Set 5. 1 Waiting times for chemical reactions (8 points)
Problem Set 5 1 Waiting times for chemical reactions (8 points) In the previous assignment, we saw that for a chemical reaction occurring at rate r, the distribution of waiting times τ between reaction
More informationPolynomial Solutions of the Laguerre Equation and Other Differential Equations Near a Singular
Polynomial Solutions of the Laguerre Equation and Other Differential Equations Near a Singular Point Abstract Lawrence E. Levine Ray Maleh Department of Mathematical Sciences Stevens Institute of Technology
More informationQ-Learning and Stochastic Approximation
MS&E338 Reinforcement Learning Lecture 4-04.11.018 Q-Learning and Stochastic Approximation Lecturer: Ben Van Roy Scribe: Christopher Lazarus Javier Sagastuy In this lecture we study the convergence of
More informationLesson 11: The Special Role of Zero in Factoring
Lesson 11: The Special Role of Zero in Factoring Student Outcomes Students find solutions to polynomial equations where the polynomial expression is not factored into linear factors. Students construct
More informationOne-Parameter Processes, Usually Functions of Time
Chapter 4 One-Parameter Processes, Usually Functions of Time Section 4.1 defines one-parameter processes, and their variations (discrete or continuous parameter, one- or two- sided parameter), including
More informationECON/FIN 250: Forecasting in Finance and Economics: Section 7: Unit Roots & Dickey-Fuller Tests
ECON/FIN 250: Forecasting in Finance and Economics: Section 7: Unit Roots & Dickey-Fuller Tests Patrick Herb Brandeis University Spring 2016 Patrick Herb (Brandeis University) Unit Root Tests ECON/FIN
More information6 Continuous-Time Birth and Death Chains
6 Continuous-Time Birth and Death Chains Angela Peace Biomathematics II MATH 5355 Spring 2017 Lecture notes follow: Allen, Linda JS. An introduction to stochastic processes with applications to biology.
More information2. Mathematical descriptions. (i) the master equation (ii) Langevin theory. 3. Single cell measurements
1. Why stochastic?. Mathematical descriptions (i) the master equation (ii) Langevin theory 3. Single cell measurements 4. Consequences Any chemical reaction is stochastic. k P d φ dp dt = k d P deterministic
More informationStochastic Hybrid Systems: Modeling, analysis, and applications to networks and biology
research supported by NSF Stochastic Hybrid Systems: Modeling, analysis, and applications to networks and biology João P. Hespanha Center for Control Engineering and Computation University of California
More informationTime Series 3. Robert Almgren. Sept. 28, 2009
Time Series 3 Robert Almgren Sept. 28, 2009 Last time we discussed two main categories of linear models, and their combination. Here w t denotes a white noise: a stationary process with E w t ) = 0, E
More informationEfficient Leaping Methods for Stochastic Chemical Systems
Efficient Leaping Methods for Stochastic Chemical Systems Ioana Cipcigan Muruhan Rathinam November 18, 28 Abstract. Well stirred chemical reaction systems which involve small numbers of molecules for some
More informationDiscrete Variables and Gradient Estimators
iscrete Variables and Gradient Estimators This assignment is designed to get you comfortable deriving gradient estimators, and optimizing distributions over discrete random variables. For most questions,
More informationModelling in Systems Biology
Modelling in Systems Biology Maria Grazia Vigliotti thanks to my students Anton Stefanek, Ahmed Guecioueur Imperial College Formal representation of chemical reactions precise qualitative and quantitative
More informationMath 1314 Lesson 4 Limits
Math 1314 Lesson 4 Limits Finding a it amounts to answering the following question: What is happening to the y-value of a function as the x-value approaches a specific target number? If the y-value is
More informationTHE IMAGINARY CUBIC PERTURBATION: A NUMERICAL AND ANALYTIC STUDY JEAN ZINN-JUSTIN
THE IMAGINARY CUBIC PERTURBATION: A NUMERICAL AND ANALYTIC STUDY JEAN ZINN-JUSTIN CEA, IRFU (irfu.cea.fr), IPhT Centre de Saclay 91191 Gif-sur-Yvette Cedex, France and Shanghai University E-mail: jean.zinn-justin@cea.fr
More informationNotes on Computer Systems for Solving Symbolic Equations
Notes on Computer Systems for Solving Symbolic Equations Richard J. Fateman March, 1991, revisited 2005 Abstract Math students learn to solve single equations symbolically beginning in elementary or junior
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationLecture 1: Brief Review on Stochastic Processes
Lecture 1: Brief Review on Stochastic Processes A stochastic process is a collection of random variables {X t (s) : t T, s S}, where T is some index set and S is the common sample space of the random variables.
More informationOn the Bilateral Laplace Transform of the positive even functions and proof of the Riemann Hypothesis. Seong Won Cha Ph.D.
On the Bilateral Laplace Transform of the positive even functions and proof of the Riemann Hypothesis Seong Won Cha Ph.D. Seongwon.cha@gmail.com Remark This is not an official paper, rather a brief report.
More informationPower series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0
Lecture 22 Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) Recall a few facts about power series: a n z n This series in z is centered at z 0. Here z can
More informationIII. Consequences of Cauchy s Theorem
MTH6 Complex Analysis 2009-0 Lecture Notes c Shaun Bullett 2009 III. Consequences of Cauchy s Theorem. Cauchy s formulae. Cauchy s Integral Formula Let f be holomorphic on and everywhere inside a simple
More informationFebruary 26, 2017 COMPLETENESS AND THE LEHMANN-SCHEFFE THEOREM
February 26, 2017 COMPLETENESS AND THE LEHMANN-SCHEFFE THEOREM Abstract. The Rao-Blacwell theorem told us how to improve an estimator. We will discuss conditions on when the Rao-Blacwellization of an estimator
More information7. Kinetics controlled by fluctuations: Kramers theory of activated processes
7. Kinetics controlled by fluctuations: Kramers theory of activated processes Macroscopic kinetic processes (time dependent concentrations) Elementary kinetic process Reaction mechanism Unimolecular processes
More informationIntroduction to Mathematical Modeling
Introduction to Mathematical Modeling - Systems Theory in the Toolbox for Systems Biology The 5 th International Course in Yeast Systems Biology 2011 June 6, 2011, PhD, Assoc Prof Head of Department Systems
More informationThe EG-family of algorithms and procedures for solving linear differential and difference higher-order systems
The EG-family of algorithms and procedures for solving linear differential and difference higher-order systems S. Abramov, D. Khmelnov, A. Ryabenko Computing Centre of the Russian Academy of Sciences Moscow
More informationDamped harmonic motion
Damped harmonic motion March 3, 016 Harmonic motion is studied in the presence of a damping force proportional to the velocity. The complex method is introduced, and the different cases of under-damping,
More informationBasic modeling approaches for biological systems. Mahesh Bule
Basic modeling approaches for biological systems Mahesh Bule The hierarchy of life from atoms to living organisms Modeling biological processes often requires accounting for action and feedback involving
More informationCheng Soon Ong & Christian Walder. Canberra February June 2018
Cheng Soon Ong & Christian Walder Research Group and College of Engineering and Computer Science Canberra February June 2018 Outlines Overview Introduction Linear Algebra Probability Linear Regression
More informationExact Arithmetic on a Computer
Exact Arithmetic on a Computer Symbolic Computation and Computer Algebra William J. Turner Department of Mathematics & Computer Science Wabash College Crawfordsville, IN 47933 Tuesday 21 September 2010
More informationQuantum Mechanics - I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras. Lecture - 9 Introducing Quantum Optics
Quantum Mechanics - I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras Lecture - 9 Introducing Quantum Optics (Refer Slide Time: 00:07) In the last lecture I gave
More informationf(s)ds, i.e. to write down the general solution in
1. First order ODEs. Simplest examples. Existence and uniqueness theorem Example 1.1. Consider the equation (1.1) x (t) = 1 This is an equation because there is an unknown: the function x(t). This is a
More informationPrediction of double gene knockout measurements
Prediction of double gene knockout measurements Sofia Kyriazopoulou-Panagiotopoulou sofiakp@stanford.edu December 12, 2008 Abstract One way to get an insight into the potential interaction between a pair
More informationMACSYMA PROGRAM FOR THE PAINLEVÉ TEST FOR NONLINEAR ORDINARY AND PARTIAL DIFFERENTIAL EQUATIONS
. MACSYMA PROGRAM FOR THE PAINLEVÉ TEST FOR NONLINEAR ORDINARY AND PARTIAL DIFFERENTIAL EQUATIONS Willy Hereman Mathematics Department and Center for the Mathematical Sciences University of Wisconsin at
More informationOn the singularities of non-linear ODEs
On the singularities of non-linear ODEs Galina Filipuk Institute of Mathematics University of Warsaw G.Filipuk@mimuw.edu.pl Collaborators: R. Halburd (London), R. Vidunas (Tokyo), R. Kycia (Kraków) 1 Plan
More informationPartial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation:
Chapter 7 Heat Equation Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation: u t = ku x x, x, t > (7.1) Here k is a constant
More informationStochastic Simulation of Biochemical Reactions
1 / 75 Stochastic Simulation of Biochemical Reactions Jorge Júlvez University of Zaragoza 2 / 75 Outline 1 Biochemical Kinetics 2 Reaction Rate Equation 3 Chemical Master Equation 4 Stochastic Simulation
More informationCh. 14 Stationary ARMA Process
Ch. 14 Stationary ARMA Process A general linear stochastic model is described that suppose a time series to be generated by a linear aggregation of random shock. For practical representation it is desirable
More informationBasic Quantum Mechanics Prof. Ajoy Ghatak Department of Physics Indian Institute of Technology, Delhi
Basic Quantum Mechanics Prof. Ajoy Ghatak Department of Physics Indian Institute of Technology, Delhi Module No. # 05 The Angular Momentum I Lecture No. # 02 The Angular Momentum Problem (Contd.) In the
More informationMATH 1231 MATHEMATICS 1B CALCULUS. Section 5: - Power Series and Taylor Series.
MATH 1231 MATHEMATICS 1B CALCULUS. Section 5: - Power Series and Taylor Series. The objective of this section is to become familiar with the theory and application of power series and Taylor series. By
More information7. MULTIVARATE STATIONARY PROCESSES
7. MULTIVARATE STATIONARY PROCESSES 1 1 Some Preliminary Definitions and Concepts Random Vector: A vector X = (X 1,..., X n ) whose components are scalar-valued random variables on the same probability
More informationExample 4.1 Let X be a random variable and f(t) a given function of time. Then. Y (t) = f(t)x. Y (t) = X sin(ωt + δ)
Chapter 4 Stochastic Processes 4. Definition In the previous chapter we studied random variables as functions on a sample space X(ω), ω Ω, without regard to how these might depend on parameters. We now
More information1 Linear Difference Equations
ARMA Handout Jialin Yu 1 Linear Difference Equations First order systems Let {ε t } t=1 denote an input sequence and {y t} t=1 sequence generated by denote an output y t = φy t 1 + ε t t = 1, 2,... with
More informationGeometric Series and the Ratio and Root Test
Geometric Series and the Ratio and Root Test James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University September 5, 2017 Outline Geometric Series The
More informationPresentation for the Visiting Committee & Graduate Seminar SYMBOLIC COMPUTING IN RESEARCH. Willy Hereman
Presentation for the Visiting Committee & Graduate Seminar SYMBOLIC COMPUTING IN RESEARCH Willy Hereman Dept. of Mathematical and Computer Sciences Colorado School of Mines Golden, CO 80401-1887 Monday,
More information1. (3) Write Gauss Law in differential form. Explain the physical meaning.
Electrodynamics I Midterm Exam - Part A - Closed Book KSU 204/0/23 Name Instructions: Use SI units. Where appropriate, define all variables or symbols you use, in words. Try to tell about the physics involved,
More informationCovariance Stationary Time Series. Example: Independent White Noise (IWN(0,σ 2 )) Y t = ε t, ε t iid N(0,σ 2 )
Covariance Stationary Time Series Stochastic Process: sequence of rv s ordered by time {Y t } {...,Y 1,Y 0,Y 1,...} Defn: {Y t } is covariance stationary if E[Y t ]μ for all t cov(y t,y t j )E[(Y t μ)(y
More informationA Sumudu based algorithm for solving differential equations
Computer Science Journal of Moldova, vol.15, no.3(45), 2007 A Sumudu based algorithm for solving differential equations Jun Zhang Abstract An algorithm based on Sumudu transform is developed. The algorithm
More informationAutomorphisms of varieties and critical density. Jason Bell (U. Michigan)
Automorphisms of varieties and critical density Jason Bell (U. Michigan) belljp@umich.edu 1 Let X be a topological space and let S be an infinite subset of X. We say that S is critically dense if every
More informationIndecomposability parameters in LCFT
Indecomposability parameters in LCFT Romain Vasseur Joint work with J.L. Jacobsen and H. Saleur at IPhT CEA Saclay and LPTENS (Nucl. Phys. B 851, 314-345 (2011), arxiv :1103.3134) ACFTA (Institut Henri
More informationNumerical Methods. King Saud University
Numerical Methods King Saud University Aims In this lecture, we will... Introduce the topic of numerical methods Consider the Error analysis and sources of errors Introduction A numerical method which
More informationMittag-Leffler and Principle of the Argument
Mittag-Leffler and Principle of the Argument Thursday, November 21, 2013 1:54 PM Homework 3 due Friday, November 22 at 5 PM. Homework 4 will be posted tonight, due Wednesday, December 11 at 5 PM. We'll
More informationComplex Variables. Instructions Solve any eight of the following ten problems. Explain your reasoning in complete sentences to maximize credit.
Instructions Solve any eight of the following ten problems. Explain your reasoning in complete sentences to maximize credit. 1. The TI-89 calculator says, reasonably enough, that x 1) 1/3 1 ) 3 = 8. lim
More informationSYMBOLIC SOFTWARE FOR SOLITON THEORY: INTEGRABILITY, SYMMETRIES CONSERVATION LAWS AND EXACT SOLUTIONS. Willy Hereman
. SYMBOLIC SOFTWARE FOR SOLITON THEORY: INTEGRABILITY, SYMMETRIES CONSERVATION LAWS AND EXACT SOLUTIONS Willy Hereman Dept. of Mathematical and Computer Sciences Colorado School of Mines Golden, Colorado
More information2.2 Separable Equations
82 2.2 Separable Equations An equation y = f(x, y) is called separable provided algebraic operations, usually multiplication, division and factorization, allow it to be written in a separable form y =
More informationFinite-Horizon Statistics for Markov chains
Analyzing FSDT Markov chains Friday, September 30, 2011 2:03 PM Simulating FSDT Markov chains, as we have said is very straightforward, either by using probability transition matrix or stochastic update
More informationAC&ST AUTOMATIC CONTROL AND SYSTEM THEORY SYSTEMS AND MODELS. Claudio Melchiorri
C. Melchiorri (DEI) Automatic Control & System Theory 1 AUTOMATIC CONTROL AND SYSTEM THEORY SYSTEMS AND MODELS Claudio Melchiorri Dipartimento di Ingegneria dell Energia Elettrica e dell Informazione (DEI)
More informationPertubation10 Page 1
Pertubation10 Page 1 Singular Perturbation Theory for Algebraic Equations Monday, February 22, 2010 1:58 PM Readings: Holmes, Sec. 1.5 Hinch, Secs. 1.2 & 1.3 Homework 1 due Thursday, February 25. Again
More informationSolving Boundary Value Problems with Homotopy Continuation
Solving Boundary Value Problems with Homotopy Continuation Gene Allgower Colorado State University Andrew Sommese University of Notre Dame Dan Bates University of Notre Dame Charles Wampler GM Research
More informationDynamic Risk Measures and Nonlinear Expectations with Markov Chain noise
Dynamic Risk Measures and Nonlinear Expectations with Markov Chain noise Robert J. Elliott 1 Samuel N. Cohen 2 1 Department of Commerce, University of South Australia 2 Mathematical Insitute, University
More informationPage 404. Lecture 22: Simple Harmonic Oscillator: Energy Basis Date Given: 2008/11/19 Date Revised: 2008/11/19
Page 404 Lecture : Simple Harmonic Oscillator: Energy Basis Date Given: 008/11/19 Date Revised: 008/11/19 Coordinate Basis Section 6. The One-Dimensional Simple Harmonic Oscillator: Coordinate Basis Page
More informationA Structured Approach Part IV. Model Checking in Systems and Synthetic Biology
A Structured Approach Part IV Model Checking in Systems and Synthetic Biology Robin Donaldson Bioinformatics Research Centre University of Glasgow, Glasgow, UK Model Checking - ISMB08 1 Outline Introduction
More informationComplex Analysis Homework 9: Solutions
Complex Analysis Fall 2007 Homework 9: Solutions 3..4 (a) Let z C \ {ni : n Z}. Then /(n 2 + z 2 ) n /n 2 n 2 n n 2 + z 2. According to the it comparison test from calculus, the series n 2 + z 2 converges
More information1 Adeles over Q. 1.1 Absolute values
1 Adeles over Q 1.1 Absolute values Definition 1.1.1 (Absolute value) An absolute value on a field F is a nonnegative real valued function on F which satisfies the conditions: (i) x = 0 if and only if
More informationUNIT-II Z-TRANSFORM. This expression is also called a one sided z-transform. This non causal sequence produces positive powers of z in X (z).
Page no: 1 UNIT-II Z-TRANSFORM The Z-Transform The direct -transform, properties of the -transform, rational -transforms, inversion of the transform, analysis of linear time-invariant systems in the -
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS Basic Concepts Paul Dawkins Table of Contents Preface... Basic Concepts... 1 Introduction... 1 Definitions... Direction Fields... 8 Final Thoughts...19 007 Paul Dawkins i http://tutorial.math.lamar.edu/terms.aspx
More informationComprehensive Examination Quantitative Methods Spring, 2018
Comprehensive Examination Quantitative Methods Spring, 2018 Instruction: This exam consists of three parts. You are required to answer all the questions in all the parts. 1 Grading policy: 1. Each part
More informationNotation Nodes are data points at which functional values are available or at which you wish to compute functional values At the nodes fx i
LECTURE 6 NUMERICAL DIFFERENTIATION To find discrete approximations to differentiation (since computers can only deal with functional values at discrete points) Uses of numerical differentiation To represent
More informationIEOR 4701: Stochastic Models in Financial Engineering. Summer 2007, Professor Whitt. SOLUTIONS to Homework Assignment 9: Brownian motion
IEOR 471: Stochastic Models in Financial Engineering Summer 27, Professor Whitt SOLUTIONS to Homework Assignment 9: Brownian motion In Ross, read Sections 1.1-1.3 and 1.6. (The total required reading there
More informationMATH 566 LECTURE NOTES 6: NORMAL FAMILIES AND THE THEOREMS OF PICARD
MATH 566 LECTURE NOTES 6: NORMAL FAMILIES AND THE THEOREMS OF PICARD TSOGTGEREL GANTUMUR 1. Introduction Suppose that we want to solve the equation f(z) = β where f is a nonconstant entire function and
More informationLecture 7: More Arithmetic and Fun With Primes
IAS/PCMI Summer Session 2000 Clay Mathematics Undergraduate Program Advanced Course on Computational Complexity Lecture 7: More Arithmetic and Fun With Primes David Mix Barrington and Alexis Maciel July
More informationLaplace Transforms Chapter 3
Laplace Transforms Important analytical method for solving linear ordinary differential equations. - Application to nonlinear ODEs? Must linearize first. Laplace transforms play a key role in important
More informationVII.5. The Weierstrass Factorization Theorem
VII.5. The Weierstrass Factorization Theorem 1 VII.5. The Weierstrass Factorization Theorem Note. Conway motivates this section with the following question: Given a sequence {a k } in G which has no limit
More informationThe simplex algorithm
The simplex algorithm The simplex algorithm is the classical method for solving linear programs. Its running time is not polynomial in the worst case. It does yield insight into linear programs, however,
More informationMORE CONSEQUENCES OF CAUCHY S THEOREM
MOE CONSEQUENCES OF CAUCHY S THEOEM Contents. The Mean Value Property and the Maximum-Modulus Principle 2. Morera s Theorem and some applications 3 3. The Schwarz eflection Principle 6 We have stated Cauchy
More informationIntroduction to Computational Stochastic Differential Equations
Introduction to Computational Stochastic Differential Equations Gabriel J. Lord Catherine E. Powell Tony Shardlow Preface Techniques for solving many of the differential equations traditionally used by
More informationWhen do diffusion-limited trajectories become memoryless?
When do diffusion-limited trajectories become memoryless? Maciej Dobrzyński CWI (Center for Mathematics and Computer Science) Kruislaan 413, 1098 SJ Amsterdam, The Netherlands Abstract Stochastic description
More informationSolving the Heat Equation (Sect. 10.5).
Solving the Heat Equation Sect. 1.5. Review: The Stationary Heat Equation. The Heat Equation. The Initial-Boundary Value Problem. The separation of variables method. An example of separation of variables.
More informationTHE LOGIC OF QUANTIFIED STATEMENTS. Predicates and Quantified Statements I. Predicates and Quantified Statements I CHAPTER 3 SECTION 3.
CHAPTER 3 THE LOGIC OF QUANTIFIED STATEMENTS SECTION 3.1 Predicates and Quantified Statements I Copyright Cengage Learning. All rights reserved. Copyright Cengage Learning. All rights reserved. Predicates
More informationExponential Moving Average Based Multiagent Reinforcement Learning Algorithms
Exponential Moving Average Based Multiagent Reinforcement Learning Algorithms Mostafa D. Awheda Department of Systems and Computer Engineering Carleton University Ottawa, Canada KS 5B6 Email: mawheda@sce.carleton.ca
More informationLet's transfer our results for conditional probability for events into conditional probabilities for random variables.
Kolmogorov/Smoluchowski equation approach to Brownian motion Tuesday, February 12, 2013 1:53 PM Readings: Gardiner, Secs. 1.2, 3.8.1, 3.8.2 Einstein Homework 1 due February 22. Conditional probability
More informationQuantum Mechanics- I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras
Quantum Mechanics- I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras Lecture - 4 Postulates of Quantum Mechanics I In today s lecture I will essentially be talking
More information3.3 It All Adds Up. A Develop Understanding Task
3.3 It All Adds Up A Develop Understanding Task Whenever we re thinking about algebra and working with variables, it is useful to consider how it relates to the number system and operations on numbers.
More informationA SUMMARY OF RECURSION SOLVING TECHNIQUES
A SUMMARY OF RECURSION SOLVING TECHNIQUES KIMMO ERIKSSON, KTH These notes are meant to be a complement to the material on recursion solving techniques in the textbook Discrete Mathematics by Biggs. In
More informationStochastic2010 Page 1
Stochastic2010 Page 1 Extinction Probability for Branching Processes Friday, April 02, 2010 2:03 PM Long-time properties for branching processes Clearly state 0 is an absorbing state, forming its own recurrent
More informationLast Update: March 1 2, 201 0
M ath 2 0 1 E S 1 W inter 2 0 1 0 Last Update: March 1 2, 201 0 S eries S olutions of Differential Equations Disclaimer: This lecture note tries to provide an alternative approach to the material in Sections
More informationComputing Monodromy Groups Defined by Plane Algebraic Curves. Adrien Poteaux
Algorithms Seminar 2006 2007, F. Chyzak (ed.), INRIA Available online at the URL http://algo.inria.fr/seminars/. Computing Monodromy Groups Defined by Plane Algebraic Curves Adrien Poteaux Xlim-Dmi, umr
More informationModelling Biochemical Pathways with Stochastic Process Algebra
Modelling Biochemical Pathways with Stochastic Process Algebra Jane Hillston. LFCS, University of Edinburgh 13th April 2007 The PEPA project The PEPA project started in Edinburgh in 1991. The PEPA project
More informationResearch Article On Maslanka s Representation for the Riemann Zeta Function
International Mathematics and Mathematical Sciences Volume 200, Article ID 7447, 9 pages doi:0.55/200/7447 Research Article On Maslana s Representation for the Riemann Zeta Function Luis Báez-Duarte Departamento
More informationON THE APPROXIMATION OF STOCHASTIC CONCURRENT CONSTRAINT PROGRAMMING BY MASTER EQUATION
ON THE APPROXIMATION OF STOCHASTIC CONCURRENT CONSTRAINT PROGRAMMING BY MASTER EQUATION Luca Bortolussi 1,2 1 Department of Mathematics and Informatics University of Trieste luca@dmi.units.it 2 Center
More informationSTOCHASTIC MODELING OF BIOCHEMICAL REACTIONS
STOCHASTIC MODELING OF BIOCHEMICAL REACTIONS Abhyudai Singh and João Pedro Hespanha* Department of Electrical and Computer Engineering University of California, Santa Barbara, CA 93101. Abstract The most
More informationGeometric Series and the Ratio and Root Test
Geometric Series and the Ratio and Root Test James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University September 5, 2018 Outline 1 Geometric Series
More information