TRAVELLING WAVES. Morteza Fotouhi Sharif Univ. of Technology

Size: px
Start display at page:

Download "TRAVELLING WAVES. Morteza Fotouhi Sharif Univ. of Technology"

Transcription

1 TRAVELLING WAVES Morteza Fotohi Sharif Univ. of Technology Mini Math NeroScience Mini Math NeroScience Agst 28

2 REACTION DIFFUSION EQUATIONS U = DU + f ( U ) t xx x t > U n D d 1 = d j > d n 2

3 Travelling wave is a soltion of the form Uxt (,) = Vx ( ct) = x ct cv = DV + f ( V ) V V = W W D 1 = [ cw + f( V)] 3

4 REST STATE W x eqilibrim V 4

5 PULSE WAVE W c x V Homoclinic Orbit 5

6 FRONT WAVE W c x Hetroclinic Orbit V 6

7 PERIODIC WAVE W c x Periodic Orbit V 7

8 EXAMPLE FitzHgh-Nagmo g v = v + f () v t t xx = βv fv () = v ( α v )( v 1) 1), < α < 1, 1 < β V = W W = cw f( V) + U U β = V c 8

9 9

10 STABILITY If U(x) is a stationary soltion of Reaction- Diffsion eqation, U t = DU xx + f ( U ) we call it stable when for every initial vale close to U in some norm, i.e. the soltion satisfies U < ε X (., t) U( + h) X < δ Frthermore, it is called asymptotically stable, if it is stable and tends towards Ux ( + h), for a constant h. (., t) U( + h) X 1

11 STABLITY OF TRAVELLING WAVES Uxt (,) = Vx ( ctt,) V = DV + cv + f ( V ) t Vx ( ct ) Travelling wave is a stationary soltion of this PDE. We mean its (asymptotically) y) stability as a stable soltion of this PDE. 11

12 LINEAR STABILITY Linearize the eqation V = DV + cv + f ( V ) t arond the stationary soltion V() V t = DV + cv + f U ( V) V L = D + c + f U ( V) 12

13 SPECTRUM Resolvent set ρ( L) { λ L λi L : D ( L ) X X = has a bonded inverse } K > : h X! U X,( L λi) U = h U X K h X 13

14 Spec ( L ) = \ ρ ( L ) = Σ Σ pt ess Σ pt : point spectrm or eigenvale defined as the kernel of L λi is nontrivial. Σ = Spec ( L )\ Σ ess pt is essential spectrm Example: L : a a1 a a1 (,, ) (,,, ) λ = Spec( L) is not an eigenvale bt becase L = (1,,, ) doesn t have 14 soltion in

15 SPECTRUM OF LINEAR EQUATION Proposition: If () is a travelling wave soltion and then V V, Spec( L). L = D + c + f U ( V) Proof. Differentiate We fid find DV + cv + f ( V ) = LV = 15

16 NONLINEAR STABILITY If is a travelling wave soltion of U t = DU xx + f ( U ) with, is a simple eigenvale of L and the other V () λ = {Re( λ) α < } spectrm are located in, then V is asymptotically stable. 16

17 CASE1: REST STATE, ν V() = V () = e n ν, \{} L = λ Sbstitte for some in the eqation to find eigenvales. We find 2 = [ + + U ( )] λ ν D cνi f V 2 d(, λν): = det[ νd + ( cν λ) I + f ( V )] U 17

18 CASE1: REST STATE, V() = V Theorem: Spec( L) = { λ d(, ik) = Proof. λ for some k } d(, λ ik), k Assme, we will show that D + c + ( f ( V ) I ) = h( ) h U λ X for every, has soltion and X K h A( λ) 18 = + v v h X

19 A ( λ ) I = 1 1 D [ λi fu ( V)] cd A( λ) 1 det[ A( λ) νi] = d( λ, ν) det D hyperbolic since d(, λ ik), k is hyperbolic E s ( λ ) = E ( λ ) = stable sbspace, nstable sbspace E s n ( λ) E ( λ) = 19

20 STRUCTURE OF SPEC(L) Spec ( L ) = { λ d ( λ, ik ) = for some k } = { λ SpecA( λ) i } f U ( V ) Spec ( L ) d ( λ, ik ) =, d λ( λ, ik ) ( ik) k k All eigenvales of lie in If, then there λ is a crve defined for sch that λ ( ik ) = λ, λ ( ik ) SpecA ( λ ) k As, we have Reλ 2

21 STRUCTURE OF SPEC(L) Spectrm lies in sector { λ < R} {arg λ > π δ} Assme that where i e ϕ λ = SpecL 2 ε < ε 1, ϕ π δ then we show that roots of i ϕ (, ) det[ 2 ( e d λν = νd + cν ) I + f ( )] 2 U V = ε ν are far from i. Let ν = ε 2 i ϕ 2 det[ ν D + ( cνε e ) I + ε f ( V )] = U 21

22 ε = e ν =± iϕ/2 For, we have is far from d j imaginary axis. And for < ε 1, ν = ν + O () ε is far too. 22

23 EXAMPLE = + a t xx v + cv + av = 2 d (, λ ik ) = k + ick + a λ = For negative parameter a< the rest wave (x,t)= is stable. 23

24 CASE 2: PERIODIC WAVE, V( + q) = V( ) = ( L λ I ) = D + c + ( f ( V ( )) λ I ) I A(, λ) = = 1 1 v v D [ λi f ( ( ))] U V cd v U A ( + q, λ ) = A (, λ ) Floqet representation B( ( ) () (, ) λ R λ e v = v 24

25 The point spectrm is empty. Spec ( L ) = { λ det( B ( λ ) ik ) = for somek } for somek = { λ Spec( B( λ)) i } Eigenfnctions are of the form () = () e ik where ( + q) = ( ) per per Sectoriality of spectrm is also tre in this case. per 25

26 CASE 3: FRONT, V() V AS ± ± L = D + c + fu ( V ( )) A(, λ) I = 1 1 D [ λi fu ( V( ))] cd E s I lim A ( λ, ) = A ( ) ± λ = 1 1 D [ λi f ( )] U V± cd ( V +, λ ) = : ( ) 2n + v v, v v V E (, λ ) = : () as where + ( ) = 2n as, where v v ( ) = v v

27 λ is in the resolvent set of L if and only if, A± ( λ) are both hyperbolic with the same Morse index, and + dim E ( λ) = dim E ( λ) E (, λ) E s (, λ) n + = λ is in the point spectrm pt, if and only if, A± ( λ) are both hyperbolic with the same Morse index, bt + dim E ( λ) = dim E ( λ) E s + (, λ) E (, λ) 27

28 λ is in the essential spectrm ess, if either at least one of the two asymptotic matrices A ( λ) is not hyperbolic, or else if it does, bt the Morse indices are different. ± 28

29 CASE 4: PULSE, V() V AS ± Special case of front wave with this different that the Morse indices are always the same. lim A( λ, ) = A ( λ ) ± 29

30 EVANS FUNCTION Choose analytic bases s { V j ( λ )} j = 1,, kand { V ( λ )} = 1,, for s E (, λ ) and E (, λ), respectively. s s 1 k 1 n k j j n k E ( λ ) = det[ V ( λ ),, V ( λ ), V ( λ ),, V ( λ )] Reslt: I. ( ) s E λ = E (, λ) E (, λ) λis an eigenvale. II. λ The order of as a zero of the Evans fnction is eqal to the algebraic mltiplicity of as an eigenvale of L. 3 λ

31 EXAMPLE Stationary soltion: = + t xx 3 qx ( ) = 2Sechx Linear eqation: v = v + (3 q ( x ) 1) v t xx 2 L = + (3 q( x) 1) xx Ax (, λ) = A () λ 2 ± = 3() qx 1 λ 1 λ Σ ess = (, 1) 31

32 1+λ λ x λ 2 ( x, λ ) = e [1+ 1+ λ tanh( x) Sech ( x )] λ x λ 2 + ( x, λ) = e [ λ tanh( x ) Sech ( x )] 3 (, λ) (, λ) + 2 E( λ) = det = λ( λ 3) 1+ λ (, λ) + (, λ) 9 Σ pt = { 1,,3} 32

33 NEURAL FIELD: INTEGRO-DIFFERENTIAL EQUATION 1 xt (, ) α t + = xt (,) + γ wyfx ()(( yt,)) dy Rest state: xt (,) = Linear Eqation: + = γ f ( ) w ( y ) dy 1 xt (, ) α t + = xt (,) + β wyx ()( ytdy,) β = γ f ( ) 33

34 L = + w()( y x y) dy β + Eigenfnctions: ikx x ( ) = e λ + 1 = βwk ˆ( ) If we assme that wy () = w( y), then wk ˆ( ) is a real even fnction of k and the stability condition is β wˆ < 1 max 34

35 REFERENCE B. Sandstede, Stability of travelling waves, Handbok of Dynamical Systems II, 22. A. I. Volpert, V. A. Volpert, and V. A. Volpert, Traveling wave soltions of parabolic systems, Translations of Mathematical Monographs 14, American Mathematical Society, L. Zhang S. Coombes C.R. Laing 35

Stability and instability of nonlinear waves:

Stability and instability of nonlinear waves: Stability and instability of nonlinear waves: Introduction 1. Nonlinear Waves 2. Stability problems 3. Stability of pulses and fronts 4. Stability of periodic waves 1 Nonlinear waves particular solutions

More information

Review of Dynamic complexity in predator-prey models framed in difference equations

Review of Dynamic complexity in predator-prey models framed in difference equations Review of Dynamic complexity in predator-prey models framed in difference eqations J. Robert Bchanan November 10, 005 Millersville University of Pennsylvania email: Robert.Bchanan@millersville.ed Review

More information

The Evans function and the stability of travelling waves

The Evans function and the stability of travelling waves The Evans function and the stability of travelling waves Jitse Niesen (University of Leeds) Collaborators: Veerle Ledoux (Ghent) Simon Malham (Heriot Watt) Vera Thümmler (Bielefeld) PANDA meeting, University

More information

The Scalar Conservation Law

The Scalar Conservation Law The Scalar Conservation Law t + f() = 0 = conserved qantity, f() =fl d dt Z b a (t, ) d = Z b a t (t, ) d = Z b a f (t, ) d = f (t, a) f (t, b) = [inflow at a] [otflow at b] f((a)) f((b)) a b Alberto Bressan

More information

Subcritical bifurcation to innitely many rotating waves. Arnd Scheel. Freie Universitat Berlin. Arnimallee Berlin, Germany

Subcritical bifurcation to innitely many rotating waves. Arnd Scheel. Freie Universitat Berlin. Arnimallee Berlin, Germany Sbcritical bifrcation to innitely many rotating waves Arnd Scheel Institt fr Mathematik I Freie Universitat Berlin Arnimallee 2-6 14195 Berlin, Germany 1 Abstract We consider the eqation 00 + 1 r 0 k2

More information

Conditions for Approaching the Origin without Intersecting the x-axis in the Liénard Plane

Conditions for Approaching the Origin without Intersecting the x-axis in the Liénard Plane Filomat 3:2 (27), 376 377 https://doi.org/.2298/fil7276a Pblished by Faclty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Conditions for Approaching

More information

Nonlinear stability of time-periodic viscous shocks. Margaret Beck Brown University

Nonlinear stability of time-periodic viscous shocks. Margaret Beck Brown University Nonlinear stability of time-periodic viscous shocks Margaret Beck Brown University Motivation Time-periodic patterns in reaction-diffusion systems: t x Experiment: chemical reaction chlorite-iodite-malonic-acid

More information

LINEAR COMBINATIONS AND SUBSPACES

LINEAR COMBINATIONS AND SUBSPACES CS131 Part II, Linear Algebra and Matrices CS131 Mathematics for Compter Scientists II Note 5 LINEAR COMBINATIONS AND SUBSPACES Linear combinations. In R 2 the vector (5, 3) can be written in the form

More information

A Single Species in One Spatial Dimension

A Single Species in One Spatial Dimension Lectre 6 A Single Species in One Spatial Dimension Reading: Material similar to that in this section of the corse appears in Sections 1. and 13.5 of James D. Mrray (), Mathematical Biology I: An introction,

More information

INPUT-OUTPUT APPROACH NUMERICAL EXAMPLES

INPUT-OUTPUT APPROACH NUMERICAL EXAMPLES INPUT-OUTPUT APPROACH NUMERICAL EXAMPLES EXERCISE s consider the linear dnamical sstem of order 2 with transfer fnction with Determine the gain 2 (H) of the inpt-otpt operator H associated with this sstem.

More information

MAXIMUM AND ANTI-MAXIMUM PRINCIPLES FOR THE P-LAPLACIAN WITH A NONLINEAR BOUNDARY CONDITION. 1. Introduction. ν = λ u p 2 u.

MAXIMUM AND ANTI-MAXIMUM PRINCIPLES FOR THE P-LAPLACIAN WITH A NONLINEAR BOUNDARY CONDITION. 1. Introduction. ν = λ u p 2 u. 2005-Ojda International Conference on Nonlinear Analysis. Electronic Jornal of Differential Eqations, Conference 14, 2006, pp. 95 107. ISSN: 1072-6691. URL: http://ejde.math.txstate.ed or http://ejde.math.nt.ed

More information

Evans function review

Evans function review Evans function review Part I: History, construction, properties and applications Strathclyde, April 18th 2005 Simon J.A. Malham http://www.ma.hw.ac.uk/ simonm/talks/ Acknowledgments T.J. Bridges, C.K.R.T.

More information

The Linear Quadratic Regulator

The Linear Quadratic Regulator 10 The Linear Qadratic Reglator 10.1 Problem formlation This chapter concerns optimal control of dynamical systems. Most of this development concerns linear models with a particlarly simple notion of optimality.

More information

Second-Order Wave Equation

Second-Order Wave Equation Second-Order Wave Eqation A. Salih Department of Aerospace Engineering Indian Institte of Space Science and Technology, Thirvananthapram 3 December 016 1 Introdction The classical wave eqation is a second-order

More information

Nonlinear stability of semidiscrete shocks for two-sided schemes

Nonlinear stability of semidiscrete shocks for two-sided schemes Nonlinear stability of semidiscrete shocks for two-sided schemes Margaret Beck Boston University Joint work with Hermen Jan Hupkes, Björn Sandstede, and Kevin Zumbrun Setting: semi-discrete conservation

More information

Control Systems

Control Systems 6.5 Control Systems Last Time: Introdction Motivation Corse Overview Project Math. Descriptions of Systems ~ Review Classification of Systems Linear Systems LTI Systems The notion of state and state variables

More information

Absolute and convective instabilities of waves on unbounded and large bounded domains

Absolute and convective instabilities of waves on unbounded and large bounded domains Physica D 145 (2000) 233 277 Absolute and convective instabilities of waves on unbounded and large bounded domains Björn Sandstede a, Arnd Scheel b, a Department of Mathematics, Ohio State University,

More information

4 Exact laminar boundary layer solutions

4 Exact laminar boundary layer solutions 4 Eact laminar bondary layer soltions 4.1 Bondary layer on a flat plate (Blasis 1908 In Sec. 3, we derived the bondary layer eqations for 2D incompressible flow of constant viscosity past a weakly crved

More information

Spectral stability of periodic waves in dispersive models

Spectral stability of periodic waves in dispersive models in dispersive models Collaborators: Th. Gallay, E. Lombardi T. Kapitula, A. Scheel in dispersive models One-dimensional nonlinear waves Standing and travelling waves u(x ct) with c = 0 or c 0 periodic

More information

Math 273b: Calculus of Variations

Math 273b: Calculus of Variations Math 273b: Calcls of Variations Yacob Kreh Homework #3 [1] Consier the 1D length fnctional minimization problem min F 1 1 L, or min 1 + 2, for twice ifferentiable fnctions : [, 1] R with bonary conitions,

More information

Homework 5 Solutions

Homework 5 Solutions Q Homework Soltions We know that the colmn space is the same as span{a & a ( a * } bt we want the basis Ths we need to make a & a ( a * linearly independent So in each of the following problems we row

More information

Applied Mathematics Letters

Applied Mathematics Letters Applied Mathematics Letters 23 (21 49 416 Contents lists available at ScienceDirect Applied Mathematics Letters jornal homepage: www.elsevier.com/locate/aml Exponential trichotomy and homoclinic bifrcation

More information

Elements of Coordinate System Transformations

Elements of Coordinate System Transformations B Elements of Coordinate System Transformations Coordinate system transformation is a powerfl tool for solving many geometrical and kinematic problems that pertain to the design of gear ctting tools and

More information

ECON3120/4120 Mathematics 2, spring 2009

ECON3120/4120 Mathematics 2, spring 2009 University of Oslo Department of Economics Arne Strøm ECON3/4 Mathematics, spring 9 Problem soltions for Seminar 4, 6 Febrary 9 (For practical reasons some of the soltions may inclde problem parts that

More information

Study Guide for Linear Algebra Exam 2

Study Guide for Linear Algebra Exam 2 Study Guide for Linear Algebra Exam 2 Term Vector Space Definition A Vector Space is a nonempty set V of objects, on which are defined two operations, called addition and multiplication by scalars (real

More information

Figure 1 Probability density function of Wedge copula for c = (best fit to Nominal skew of DRAM case study).

Figure 1 Probability density function of Wedge copula for c = (best fit to Nominal skew of DRAM case study). Wedge Copla This docment explains the constrction and properties o a particlar geometrical copla sed to it dependency data rom the edram case stdy done at Portland State University. The probability density

More information

Exponential Stability of the Traveling Fronts for a Pseudo-Para. Pseudo-Parabolic Fisher-KPP Equation

Exponential Stability of the Traveling Fronts for a Pseudo-Para. Pseudo-Parabolic Fisher-KPP Equation Exponential Stability of the Traveling Fronts for a Pseudo-Parabolic Fisher-KPP Equation Based on joint work with Xueli Bai (Center for PDE, East China Normal Univ.) and Yang Cao (Dalian University of

More information

Formal Methods for Deriving Element Equations

Formal Methods for Deriving Element Equations Formal Methods for Deriving Element Eqations And the importance of Shape Fnctions Formal Methods In previos lectres we obtained a bar element s stiffness eqations sing the Direct Method to obtain eact

More information

Math 263 Assignment #3 Solutions. 1. A function z = f(x, y) is called harmonic if it satisfies Laplace s equation:

Math 263 Assignment #3 Solutions. 1. A function z = f(x, y) is called harmonic if it satisfies Laplace s equation: Math 263 Assignment #3 Soltions 1. A fnction z f(x, ) is called harmonic if it satisfies Laplace s eqation: 2 + 2 z 2 0 Determine whether or not the following are harmonic. (a) z x 2 + 2. We se the one-variable

More information

Mean Value Formulae for Laplace and Heat Equation

Mean Value Formulae for Laplace and Heat Equation Mean Vale Formlae for Laplace and Heat Eqation Abhinav Parihar December 7, 03 Abstract Here I discss a method to constrct the mean vale theorem for the heat eqation. To constrct sch a formla ab initio,

More information

TRANSONIC EVAPORATION WAVES IN A SPHERICALLY SYMMETRIC NOZZLE

TRANSONIC EVAPORATION WAVES IN A SPHERICALLY SYMMETRIC NOZZLE TRANSONIC EVAPORATION WAVES IN A SPHERICALLY SYMMETRIC NOZZLE XIAOBIAO LIN AND MARTIN WECHSELBERGER Abstract. This paper stdies the liqid to vapor phase transition in a cone shaped nozzle. Using the geometric

More information

OPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIELD OF A POLYHEDRAL BODY WITH LINEARLY INCREASING DENSITY 1

OPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIELD OF A POLYHEDRAL BODY WITH LINEARLY INCREASING DENSITY 1 OPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIEL OF A POLYHERAL BOY WITH LINEARLY INCREASING ENSITY 1 V. POHÁNKA2 Abstract The formla for the comptation of the gravity field of a polyhedral body

More information

Reduction of the HH-Model to two dimensions (general)

Reduction of the HH-Model to two dimensions (general) Spiking models 1 Neronal codes Spiking models: Hodgkin Hxley Model (brief repetition) Redction of the HH-Model to two dimensions (general) FitzHgh-Nagmo Model Integrate and Fire Model 2 Varying firing

More information

Homogeneous Liner Systems with Constant Coefficients

Homogeneous Liner Systems with Constant Coefficients Homogeneos Liner Systems with Constant Coefficients Jly, 06 The object of stdy in this section is where A is a d d constant matrix whose entries are real nmbers. As before, we will look to the exponential

More information

Partial Differential Equations with Applications

Partial Differential Equations with Applications Universit of Leeds MATH 33 Partial Differential Eqations with Applications Eamples to spplement Chapter on First Order PDEs Eample (Simple linear eqation, k + = 0, (, 0) = ϕ(), k a constant.) The characteristic

More information

MATH 304 Linear Algebra Lecture 23: Diagonalization. Review for Test 2.

MATH 304 Linear Algebra Lecture 23: Diagonalization. Review for Test 2. MATH 304 Linear Algebra Lecture 23: Diagonalization. Review for Test 2. Diagonalization Let L be a linear operator on a finite-dimensional vector space V. Then the following conditions are equivalent:

More information

Stability Analysis of Stationary Solutions for the Cahn Hilliard Equation

Stability Analysis of Stationary Solutions for the Cahn Hilliard Equation Stability Analysis of Stationary Solutions for the Cahn Hilliard Equation Peter Howard, Texas A&M University University of Louisville, Oct. 19, 2007 References d = 1: Commun. Math. Phys. 269 (2007) 765

More information

Travelling waves. Chapter 8. 1 Introduction

Travelling waves. Chapter 8. 1 Introduction Chapter 8 Travelling waves 1 Introduction One of the cornerstones in the study of both linear and nonlinear PDEs is the wave propagation. A wave is a recognizable signal which is transferred from one part

More information

AM 205: lecture 14. Last time: Boundary value problems Today: Numerical solution of PDEs

AM 205: lecture 14. Last time: Boundary value problems Today: Numerical solution of PDEs AM 205: lecture 14 Last time: Boundary value problems Today: Numerical solution of PDEs ODE BVPs A more general approach is to formulate a coupled system of equations for the BVP based on a finite difference

More information

Optimization via the Hamilton-Jacobi-Bellman Method: Theory and Applications

Optimization via the Hamilton-Jacobi-Bellman Method: Theory and Applications Optimization via the Hamilton-Jacobi-Bellman Method: Theory and Applications Navin Khaneja lectre notes taken by Christiane Koch Jne 24, 29 1 Variation yields a classical Hamiltonian system Sppose that

More information

Numerical Evans Function Computation Jeffrey Humpherys Brigham Young University INTRODUCTION

Numerical Evans Function Computation Jeffrey Humpherys Brigham Young University INTRODUCTION Numerical Evans Function Computation Jeffrey Humpherys Brigham Young University INTRODUCTION 1 Traveling Wave Consider the evolution equation u t = F(u, u x,u xx,...). A traveling wave with speed is a

More information

STEP Support Programme. STEP III Hyperbolic Functions: Solutions

STEP Support Programme. STEP III Hyperbolic Functions: Solutions STEP Spport Programme STEP III Hyperbolic Fnctions: Soltions Start by sing the sbstittion t cosh x. This gives: sinh x cosh a cosh x cosh a sinh x t sinh x dt t dt t + ln t ln t + ln cosh a ln ln cosh

More information

HADAMARD-PERRON THEOREM

HADAMARD-PERRON THEOREM HADAMARD-PERRON THEOREM CARLANGELO LIVERANI. Invariant manifold of a fixed point He we will discss the simplest possible case in which the existence of invariant manifolds arises: the Hadamard-Perron theorem.

More information

MATH 425, HOMEWORK 3 SOLUTIONS

MATH 425, HOMEWORK 3 SOLUTIONS MATH 425, HOMEWORK 3 SOLUTIONS Exercise. (The differentiation property of the heat equation In this exercise, we will use the fact that the derivative of a solution to the heat equation again solves the

More information

Linear System Theory (Fall 2011): Homework 1. Solutions

Linear System Theory (Fall 2011): Homework 1. Solutions Linear System Theory (Fall 20): Homework Soltions De Sep. 29, 20 Exercise (C.T. Chen: Ex.3-8). Consider a linear system with inpt and otpt y. Three experiments are performed on this system sing the inpts

More information

Coherent structures near the boundary between excitable and oscillatory media

Coherent structures near the boundary between excitable and oscillatory media Coherent structures near the boundary between excitable and oscillatory media Jeremy Bellay University of Minnesota Department of Computer Science 200 Union St. S.E. Minneapolis, MN 55455, USA Arnd Scheel

More information

KdV PRESERVES WHITE NOISE

KdV PRESERVES WHITE NOISE KdV PESEVES WHIE NOISE JEEMY QUASEL AND BENEDEK VALKÓ Abstract. It is shown that white noise is an invariant measre for the Korteweg-deVries eqation on. his is a conseqence of recent reslts of Kappeler

More information

Linear Algebra II Lecture 8

Linear Algebra II Lecture 8 Linear Algebra II Lecture 8 Xi Chen 1 1 University of Alberta October 10, 2014 Outline 1 2 Definition Let T 1 : V W and T 2 : V W be linear transformations between two vector spaces V and W over R. Then

More information

Solving First Order PDEs

Solving First Order PDEs Solving Ryan C. Trinity University Partial Differential Equations January 21, 2014 Solving the transport equation Goal: Determine every function u(x, t) that solves u t +v u x = 0, where v is a fixed constant.

More information

ON THE EIGENVALUE OF INFINITE MATRICES WITH NONNEGATIVE OFF-DIAGONAL ELEMENTS

ON THE EIGENVALUE OF INFINITE MATRICES WITH NONNEGATIVE OFF-DIAGONAL ELEMENTS ON THE EIGENVALUE OF INFINITE MATRICES WITH NONNEGATIVE OFF-DIAGONAL ELEMENTS N. APREUTESEI AND V. VOLPERT The paper is devoted to infinite-dimensional difference operators. Some spectral properties of

More information

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 4 (Products and quotients) & (Logarithmic differentiation) A.J.Hobson

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 4 (Products and quotients) & (Logarithmic differentiation) A.J.Hobson JUST THE MATHS UNIT NUMBER 104 DIFFERENTIATION 4 (Products and quotients) & (Logarithmic differentiation) by AJHobson 1041 Products 1042 Quotients 1043 Logarithmic differentiation 1044 Exercises 1045 Answers

More information

Solving a System of Equations

Solving a System of Equations Solving a System of Eqations Objectives Understand how to solve a system of eqations with: - Gass Elimination Method - LU Decomposition Method - Gass-Seidel Method - Jacobi Method A system of linear algebraic

More information

Lecture 9: 3.4 The Geometry of Linear Systems

Lecture 9: 3.4 The Geometry of Linear Systems Lectre 9: 3.4 The Geometry of Linear Systems Wei-Ta Ch 200/0/5 Dot Prodct Form of a Linear System Recall that a linear eqation has the form a x +a 2 x 2 + +a n x n = b (a,a 2,, a n not all zero) The corresponding

More information

CRITERIA FOR TOEPLITZ OPERATORS ON THE SPHERE. Jingbo Xia

CRITERIA FOR TOEPLITZ OPERATORS ON THE SPHERE. Jingbo Xia CRITERIA FOR TOEPLITZ OPERATORS ON THE SPHERE Jingbo Xia Abstract. Let H 2 (S) be the Hardy space on the nit sphere S in C n. We show that a set of inner fnctions Λ is sfficient for the prpose of determining

More information

Chapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS

Chapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS Chapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS 3. System Modeling Mathematical Modeling In designing control systems we mst be able to model engineered system dynamics. The model of a dynamic system

More information

Asymptotics of dissipative nonlinear evolution equations with ellipticity: different end states

Asymptotics of dissipative nonlinear evolution equations with ellipticity: different end states J. Math. Anal. Appl. 33 5) 5 35 www.elsevier.com/locate/jmaa Asymptotics of dissipative nonlinear evoltion eqations with ellipticity: different end states enjn Dan, Changjiang Zh Laboratory of Nonlinear

More information

THE HOHENBERG-KOHN THEOREM FOR MARKOV SEMIGROUPS

THE HOHENBERG-KOHN THEOREM FOR MARKOV SEMIGROUPS THE HOHENBERG-KOHN THEOREM FOR MARKOV SEMIGROUPS OMAR HIJAB Abstract. At the basis of mch of comptational chemistry is density fnctional theory, as initiated by the Hohenberg-Kohn theorem. The theorem

More information

Advanced topics in Finite Element Method 3D truss structures. Jerzy Podgórski

Advanced topics in Finite Element Method 3D truss structures. Jerzy Podgórski Advanced topics in Finite Element Method 3D trss strctres Jerzy Podgórski Introdction Althogh 3D trss strctres have been arond for a long time, they have been sed very rarely ntil now. They are difficlt

More information

Integration of Basic Functions. Session 7 : 9/23 1

Integration of Basic Functions. Session 7 : 9/23 1 Integration o Basic Fnctions Session 7 : 9/3 Antiderivation Integration Deinition: Taking the antiderivative, or integral, o some nction F(), reslts in the nction () i ()F() Pt simply: i yo take the integral

More information

1 Undiscounted Problem (Deterministic)

1 Undiscounted Problem (Deterministic) Lectre 9: Linear Qadratic Control Problems 1 Undisconted Problem (Deterministic) Choose ( t ) 0 to Minimize (x trx t + tq t ) t=0 sbject to x t+1 = Ax t + B t, x 0 given. x t is an n-vector state, t a

More information

Solutions for Math 225 Assignment #5 1

Solutions for Math 225 Assignment #5 1 Solutions for Math 225 Assignment #5 1 (1) Find a polynomial f(x) of degree at most 3 satisfying that f(0) = 2, f( 1) = 1, f(1) = 3 and f(3) = 1. Solution. By Lagrange Interpolation, ( ) (x + 1)(x 1)(x

More information

Spatial decay of rotating waves in parabolic systems

Spatial decay of rotating waves in parabolic systems Spatial decay of rotating waves in parabolic systems Nonlinear Waves, CRC 701, Bielefeld, June 19, 2013 Denny Otten Department of Mathematics Bielefeld University Germany June 19, 2013 CRC 701 Denny Otten

More information

1 Existence of Travelling Wave Fronts for a Reaction-Diffusion Equation with Quadratic- Type Kinetics

1 Existence of Travelling Wave Fronts for a Reaction-Diffusion Equation with Quadratic- Type Kinetics 1 Existence of Travelling Wave Fronts for a Reaction-Diffusion Equation with Quadratic- Type Kinetics Theorem. Consider the equation u t = Du xx + f(u) with f(0) = f(1) = 0, f(u) > 0 on 0 < u < 1, f (0)

More information

Dispersion relations, stability and linearization

Dispersion relations, stability and linearization Dispersion relations, stability and linearization 1 Dispersion relations Suppose that u(x, t) is a function with domain { < x 0}, and it satisfies a linear, constant coefficient partial differential

More information

Nonlinear convective stability of travelling fronts near Turing and Hopf instabilities

Nonlinear convective stability of travelling fronts near Turing and Hopf instabilities Nonlinear convective stability of travelling fronts near Turing and Hopf instabilities Margaret Beck Joint work with Anna Ghazaryan, University of Kansas and Björn Sandstede, Brown University September

More information

(i) [7 points] Compute the determinant of the following matrix using cofactor expansion.

(i) [7 points] Compute the determinant of the following matrix using cofactor expansion. Question (i) 7 points] Compute the determinant of the following matrix using cofactor expansion 2 4 2 4 2 Solution: Expand down the second column, since it has the most zeros We get 2 4 determinant = +det

More information

Geometric Image Manipulation. Lecture #4 Wednesday, January 24, 2018

Geometric Image Manipulation. Lecture #4 Wednesday, January 24, 2018 Geometric Image Maniplation Lectre 4 Wednesda, Janar 4, 08 Programming Assignment Image Maniplation: Contet To start with the obvios, an image is a D arra of piels Piel locations represent points on the

More information

Vectors in Rn un. This definition of norm is an extension of the Pythagorean Theorem. Consider the vector u = (5, 8) in R 2

Vectors in Rn un. This definition of norm is an extension of the Pythagorean Theorem. Consider the vector u = (5, 8) in R 2 MATH 307 Vectors in Rn Dr. Neal, WKU Matrices of dimension 1 n can be thoght of as coordinates, or ectors, in n- dimensional space R n. We can perform special calclations on these ectors. In particlar,

More information

Section 7.4: Integration of Rational Functions by Partial Fractions

Section 7.4: Integration of Rational Functions by Partial Fractions Section 7.4: Integration of Rational Fnctions by Partial Fractions This is abot as complicated as it gets. The Method of Partial Fractions Ecept for a few very special cases, crrently we have no way to

More information

Analyticity of semigroups generated by Fleming-Viot type operators

Analyticity of semigroups generated by Fleming-Viot type operators Analyticity of semigroups generated by Fleming-Viot type operators Elisabetta Mangino, in collaboration with A. Albanese Università del Salento, Lecce, Italy s Au(x) = x i (δ ij x j )D ij u + b i (x)d

More information

Diagonalization of the Coupled-Mode System.

Diagonalization of the Coupled-Mode System. Diagonalization of the Coupled-Mode System. Marina Chugunova joint work with Dmitry Pelinovsky Department of Mathematics, McMaster University, Canada Collaborators: Mason A. Porter, California Institute

More information

Solutions to Math 152 Review Problems for Exam 1

Solutions to Math 152 Review Problems for Exam 1 Soltions to Math 5 Review Problems for Eam () If A() is the area of the rectangle formed when the solid is sliced at perpendiclar to the -ais, then A() = ( ), becase the height of the rectangle is and

More information

Applied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.

Applied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R

More information

FEA Solution Procedure

FEA Solution Procedure EA Soltion Procedre (demonstrated with a -D bar element problem) EA Procedre for Static Analysis. Prepare the E model a. discretize (mesh) the strctre b. prescribe loads c. prescribe spports. Perform calclations

More information

Continuous dependence estimates for the ergodic problem with an application to homogenization

Continuous dependence estimates for the ergodic problem with an application to homogenization Continuous dependence estimates for the ergodic problem with an application to homogenization Claudio Marchi Bayreuth, September 12 th, 2013 C. Marchi (Università di Padova) Continuous dependence Bayreuth,

More information

Characterizations of probability distributions via bivariate regression of record values

Characterizations of probability distributions via bivariate regression of record values Metrika (2008) 68:51 64 DOI 10.1007/s00184-007-0142-7 Characterizations of probability distribtions via bivariate regression of record vales George P. Yanev M. Ahsanllah M. I. Beg Received: 4 October 2006

More information

On Multiobjective Duality For Variational Problems

On Multiobjective Duality For Variational Problems The Open Operational Research Jornal, 202, 6, -8 On Mltiobjective Dality For Variational Problems. Hsain *,, Bilal Ahmad 2 and Z. Jabeen 3 Open Access Department of Mathematics, Jaypee University of Engineering

More information

Math 116 First Midterm October 14, 2009

Math 116 First Midterm October 14, 2009 Math 116 First Midterm October 14, 9 Name: EXAM SOLUTIONS Instrctor: Section: 1. Do not open this exam ntil yo are told to do so.. This exam has 1 pages inclding this cover. There are 9 problems. Note

More information

SUBORDINATION RESULTS FOR A CERTAIN SUBCLASS OF NON-BAZILEVIC ANALYTIC FUNCTIONS DEFINED BY LINEAR OPERATOR

SUBORDINATION RESULTS FOR A CERTAIN SUBCLASS OF NON-BAZILEVIC ANALYTIC FUNCTIONS DEFINED BY LINEAR OPERATOR italian jornal of pre and applied mathematics n. 34 215 375 388) 375 SUBORDINATION RESULTS FOR A CERTAIN SUBCLASS OF NON-BAZILEVIC ANALYTIC FUNCTIONS DEFINED BY LINEAR OPERATOR Adnan G. Alamosh Maslina

More information

Math Ordinary Differential Equations

Math Ordinary Differential Equations Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x

More information

On continuous time contract theory

On continuous time contract theory Ecole Polytechnique, France Journée de rentrée du CMAP, 3 octobre, 218 Outline 1 2 Semimartingale measures on the canonical space Random horizon 2nd order backward SDEs (Static) Principal-Agent Problem

More information

1 Introduction Travelling-wave solutions of parabolic equations on the real line arise in a variety of applications. An important issue is their stabi

1 Introduction Travelling-wave solutions of parabolic equations on the real line arise in a variety of applications. An important issue is their stabi Essential instability of pulses, and bifurcations to modulated travelling waves Bjorn Sandstede Department of Mathematics The Ohio State University 231 West 18th Avenue Columbus, OH 4321, USA Arnd Scheel

More information

Ground Rules. PC1221 Fundamentals of Physics I. Position and Displacement. Average Velocity. Lectures 7 and 8 Motion in Two Dimensions

Ground Rules. PC1221 Fundamentals of Physics I. Position and Displacement. Average Velocity. Lectures 7 and 8 Motion in Two Dimensions PC11 Fndamentals of Physics I Lectres 7 and 8 Motion in Two Dimensions A/Prof Tay Sen Chan 1 Grond Rles Switch off yor handphone and paer Switch off yor laptop compter and keep it No talkin while lectre

More information

New approach for tanh and extended-tanh methods with applications on Hirota-Satsuma equations

New approach for tanh and extended-tanh methods with applications on Hirota-Satsuma equations Volume 28, N. 1, pp. 1 14, 2009 Copyright 2009 SBMAC ISSN 0101-8205 www.scielo.br/cam New approach for tanh and extended-tanh methods with applications on Hirota-Satsuma equations HASSAN A. ZEDAN Mathematics

More information

Simplified Identification Scheme for Structures on a Flexible Base

Simplified Identification Scheme for Structures on a Flexible Base Simplified Identification Scheme for Strctres on a Flexible Base L.M. Star California State University, Long Beach G. Mylonais University of Patras, Greece J.P. Stewart University of California, Los Angeles

More information

Solving First Order PDEs

Solving First Order PDEs Solving Ryan C. Trinity University Partial Differential Equations Lecture 2 Solving the transport equation Goal: Determine every function u(x, t) that solves u t +v u x = 0, where v is a fixed constant.

More information

Minimal periods of semilinear evolution equations with Lipschitz nonlinearity

Minimal periods of semilinear evolution equations with Lipschitz nonlinearity Minimal periods of semilinear evolution equations with Lipschitz nonlinearity James C. Robinson a Alejandro Vidal-López b a Mathematics Institute, University of Warwick, Coventry, CV4 7AL, U.K. b Departamento

More information

Solutions to the Calculus and Linear Algebra problems on the Comprehensive Examination of January 28, 2011

Solutions to the Calculus and Linear Algebra problems on the Comprehensive Examination of January 28, 2011 Solutions to the Calculus and Linear Algebra problems on the Comprehensive Examination of January 8, Solutions to Problems 5 are omitted since they involve topics no longer covered on the Comprehensive

More information

Partial Differential Equations, Winter 2015

Partial Differential Equations, Winter 2015 Partial Differential Equations, Winter 215 Homework #2 Due: Thursday, February 12th, 215 1. (Chapter 2.1) Solve u xx + u xt 2u tt =, u(x, ) = φ(x), u t (x, ) = ψ(x). Hint: Factor the operator as we did

More information

Solving the Lienard equation by differential transform method

Solving the Lienard equation by differential transform method ISSN 1 746-7233, England, U World Jornal of Modelling and Simlation Vol. 8 (2012) No. 2, pp. 142-146 Solving the Lienard eqation by differential transform method Mashallah Matinfar, Saber Rakhshan Bahar,

More information

Turning points and traveling waves in FitzHugh-Nagumo type equations

Turning points and traveling waves in FitzHugh-Nagumo type equations Turning points and traveling waves in FitzHugh-Nagumo type equations Weishi Liu and Erik Van Vleck Department of Mathematics University of Kansas, Lawrence, KS 66045 E-mail: wliu@math.ku.edu, evanvleck@math.ku.edu

More information

(1.1) g(x) = f'ek(xy)fy)dy,

(1.1) g(x) = f'ek(xy)fy)dy, THE G-FÜNCTIONS AS UNSYMMETRICAL FOURIER KERNELS. Ill ROOP NARAIN 1. The fnctions k(x) and h(x) are said to form a pair of Forier kernels if the reciprocal eqations (1.1) g(x) = f'ek(xy)fy)dy, J o (1.1')

More information

Convergence to diffusion waves for nonlinear evolution equations with different end states

Convergence to diffusion waves for nonlinear evolution equations with different end states J. Math. Anal. Appl. 338 8) 44 63 www.elsevier.com/locate/jmaa Convergence to diffsion waves for nonlinear evoltion eqations with different end states Walter Allegretto, Yanping Lin, Zhiyong Zhang Department

More information

Theorem 1. ẋ = Ax is globally exponentially stable (GES) iff A is Hurwitz (i.e., max(re(σ(a))) < 0).

Theorem 1. ẋ = Ax is globally exponentially stable (GES) iff A is Hurwitz (i.e., max(re(σ(a))) < 0). Linear Systems Notes Lecture Proposition. A M n (R) is positive definite iff all nested minors are greater than or equal to zero. n Proof. ( ): Positive definite iff λ i >. Let det(a) = λj and H = {x D

More information

Analytical Investigation of Hyperbolic Equations via He s Methods

Analytical Investigation of Hyperbolic Equations via He s Methods American J. of Engineering and Applied Sciences (4): 399-47, 8 ISSN 94-7 8 Science Pblications Analytical Investigation of Hyperbolic Eqations via He s Methods D.D. Ganji, M. Amini and A. Kolahdooz Department

More information

1 The space of linear transformations from R n to R m :

1 The space of linear transformations from R n to R m : Math 540 Spring 20 Notes #4 Higher deriaties, Taylor s theorem The space of linear transformations from R n to R m We hae discssed linear transformations mapping R n to R m We can add sch linear transformations

More information

Lecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.

Lecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University. Lecture 4 Chapter 4: Lyapunov Stability Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 4 p. 1/86 Autonomous Systems Consider the autonomous system ẋ

More information

Group Method. December 16, Oberwolfach workshop Dynamics of Patterns

Group Method. December 16, Oberwolfach workshop Dynamics of Patterns CWI, Amsterdam heijster@cwi.nl December 6, 28 Oberwolfach workshop Dynamics of Patterns Joint work: A. Doelman (CWI/UvA), T.J. Kaper (BU), K. Promislow (MSU) Outline 2 3 4 Interactions of localized structures

More information

Linear Algebra, part 2 Eigenvalues, eigenvectors and least squares solutions

Linear Algebra, part 2 Eigenvalues, eigenvectors and least squares solutions Linear Algebra, part 2 Eigenvalues, eigenvectors and least squares solutions Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Main problem of linear algebra 2: Given

More information

Slow-fast dynamics in Josephson junctions

Slow-fast dynamics in Josephson junctions Er. Phys.. B 3, 93 33 (3) DOI:./epjb/e3-5- THE EUROPEAN PHYSICAL OURNAL B Slow-fast dynamics in osephson jnctions E. Nemann a and A. Pikovsky Department of Physics, University of Potsdam, Am Neen Palais,

More information