Dynamic System In Biology

Size: px
Start display at page:

Download "Dynamic System In Biology"

Transcription

1 Compuaional Siene and Engineering Dnami Ssem In Biolog Yang Cao Deparmen of Compuer Siene hp://ourses.s.v.edu/~s644

2 Ouline Compuaional Siene and Engineering Single Speies opulaion Model Malhus Model Logisi Model Two speies Model Compeiion Model redaor and re Model hase lo and Dnami Ssem Sohasi Model and Simulaion Loka Model Brusselaor Model

3 Malhus Model Compuaional Siene and Engineering I SAID ha populaion, when unheked, inreased in a geomerial raio, and subsisene for man in an arihmeial raio Thomas Malhus

4 Malhus Model Compuaional Siene and Engineering The reproduion rae is proporional o he populaion k Solve i we have e k The populaion in he Unied Saes in ear 6 79 is.9. The orresponding populaion in ear 8 is Wih a daa fiing, we obain:.9 6 e.7 79

5 Compuaional Siene and Engineering Logisi opulaion Model Developed b Belgian mahemaiian ierre Verhuls 88 in 88 The rae of populaion inrease ma be limied, i.e., i ma depend on populaion densi k where k k m k m m m k m k e e e The soluion is

6 Compuaional Siene and Engineering Logisi opulaion Model k m m m k m k e e e The soluion of he Logisi model Wih a daa fiing.4, 97 6 k m

7 Model of wo speies Compeiion Compuaional Siene and Engineering Le he populaion of wo speies be and, and he ompee in he same environmen. If here is no ompeiion, he populaion of X will saisf & r N Wih he ompeiion, & For anoher speies, here is a similar equaion r α N The phsial meaning of and an be undersood as: α he he & r α resoure resoure β eah eah N β X speies Y speies onsume onsume. Thus we have αβ

8 Sae Dnamis lo vs hase lo Compuaional Siene and Engineering Sae Dnamis lo: sae vs ime, hase lo: he sae spae, use arrow o represen he angen veor The phase plo reveals he geomeri proper of a dnami ssem represened b a pair of ODEs. &. &.

9 Sae Dnamis lo vs hase lo Compuaional Siene and Engineering Eample: from differen iniial value, he rajeor follow he direion of he arrows and reahes o is equilibrium sae

10 Sae Dnamis lo vs hase lo Compuaional Siene and Engineering However, a sligh hange of parameers make a big differene in phase plo and lead o a differen onlusion &. &. 9

11 Compuaional Siene and Engineering Sae Dnamis lo vs hase lo 9.. & &

12 Sae Dnamis lo vs hase lo Compuaional Siene and Engineering & & r α N β N The sign of he derivaives are deided b wo values N A dire analsis hrough he phase plo α and α N α r If N > α, X speies will win. N If N < αn, Y speies will win.

13 Model of wo speies redaor and re Compuaional Siene and Engineering Loka-Volerra Model The simples model of predaor-pre ineraions developed independenl b Loka 95 and Volerra 96 Anona s observaion on Shark s populaion during world war I.

14 Model of wo speies redaor and re Compuaional Siene and Engineering Assumpion: The predaor speies is oall dependen on a single pre speies as is onl food suppl, The pre speies has an unlimied food suppl, and here is no hrea o he pre oher han he speifi predaor. Le X represen he pre and Y represen he predaor, wihou he predaor, he Malhus model an be applied However, beause of he predaor, r has o be modified For he predaor, he siuaion is jus he opposie. Thus we ge he ODEs for his model & a & a b & d & a b & d

15 hase lo Analsis & & a b d Compuaional Siene and Engineering There are wo orresponding equilibrium poins:,, or d a b,, a b,,, d

16 Malab Simulaion Resul Compuaional Siene and Engineering Based on eample: &. &.4

17 Effe of arameers Compuaional Siene and Engineering The soluion of he LV predaor-pre model is, d a b where a : b : : d : he naural reproduion rae for he pre he killing rae beause of he predaor he naural deah rae for he predaor he reproduion rae beause of he pre Quesion: Wh he shark raio inreases during world war I?

18 arameer Analsis * *, d a b Compuaional Siene and Engineering When fishing is inrodued in he model, heir effe will be inrease he deah rae of he predaor and derease he reproduion rae for he pre. Thus e, a a e,

19 Compuaional Siene and Engineering Sohasi Modeling Loka reaions: Z Y Y Y X X X A Lead o ODEs A & & The sohasi simulaion generaes ineresing rajeories..,, A

20 Differen Dnami Behavior Compuaional Siene and Engineering

21 Compuaional Siene and Engineering Brusselaor D Y Y Y X C Y X B X A 4 B B A 4 & & 5..5, 5, 5, 4 B A Lead o ODEs 5.., 5, 5, 4 B A Bifuraion happens around he ondiion: 4 4 A B J. Tson s 97, 974 paper

22 Compuaional Siene and Engineering Oregonaor X Z E D Y Z Y Y C B Y X Y X A 5 4 Ez C z C A Ez A & & & 6.6, 4,.,, 5 4 E C A

23 Oregonaor Compuaional Siene and Engineering A X Y X Y B C Y Y Y 4 D E Z 5 X Z A,., C 4, 4.6, 5E 6

24 Compuaional Siene and Engineering Thanks! Quesions? lao is m friend, Arisole is m friend, bu m bes friend is ruh --- Newon

Outline of Topics. Analysis of ODE models with MATLAB. What will we learn from this lecture. Aim of analysis: Why such analysis matters?

Outline of Topics. Analysis of ODE models with MATLAB. What will we learn from this lecture. Aim of analysis: Why such analysis matters? of Topics wih MATLAB Shan He School for Compuaional Science Universi of Birmingham Module 6-3836: Compuaional Modelling wih MATLAB Wha will we learn from his lecure Aim of analsis: Aim of analsis. Some

More information

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

Simulation-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 information

1 birth rate γ (number of births per time interval) 2 death rate δ proportional to size of population

1 birth rate γ (number of births per time interval) 2 death rate δ proportional to size of population Scienific Comuing I Module : Poulaion Modelling Coninuous Models Michael Bader Par I ODE Models Lehrsuhl Informaik V Winer 7/ Discree vs. Coniuous Models d d = F,,...) ) =? discree model: coninuous model:

More information

Elementary Differential Equations and Boundary Value Problems

Elementary 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 information

Predator - Prey Model Trajectories and the nonlinear conservation law

Predator - 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 information

Comparing Theoretical and Practical Solution of the First Order First Degree Ordinary Differential Equation of Population Model

Comparing Theoretical and Practical Solution of the First Order First Degree Ordinary Differential Equation of Population Model Open Access Journal of Mahemaical and Theoreical Physics Comparing Theoreical and Pracical Soluion of he Firs Order Firs Degree Ordinary Differenial Equaion of Populaion Model Absrac Populaion dynamics

More information

Class Notes 1: Introduction. MAE 82 Engineering Mathematics

Class Notes 1: Introduction. MAE 82 Engineering Mathematics Class Noes 1: Inroducion MAE 82 Engineering Mahemaics CHANGE Rae of Change Basic Mahemaical Models Man of he principles or laws underling he behavior of he naural world are saemens or relaions involving

More information

Math 2214 Solution Test 1B Fall 2017

Math 2214 Solution Test 1B Fall 2017 Mah 14 Soluion Tes 1B Fall 017 Problem 1: A ank has a capaci for 500 gallons and conains 0 gallons of waer wih lbs of sal iniiall. A soluion conaining of 8 lbsgal of sal is pumped ino he ank a 10 galsmin.

More information

Math 2214 Solution Test 1A Spring 2016

Math 2214 Solution Test 1A Spring 2016 Mah 14 Soluion Tes 1A Spring 016 sec Problem 1: Wha is he larges -inerval for which ( 4) = has a guaraneed + unique soluion for iniial value (-1) = 3 according o he Exisence Uniqueness Theorem? Soluion

More information

Stochastic Model for Cancer Cell Growth through Single Forward Mutation

Stochastic 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 information

An Inventory Model for Weibull Time-Dependence. Demand Rate with Completely Backlogged. Shortages

An Inventory Model for Weibull Time-Dependence. Demand Rate with Completely Backlogged. Shortages Inernaional Mahemaial Forum, 5, 00, no. 5, 675-687 An Invenory Model for Weibull Time-Dependene Demand Rae wih Compleely Baklogged Shorages C. K. Tripahy and U. Mishra Deparmen of Saisis, Sambalpur Universiy

More information

Cosumnes River College Principles of Macroeconomics Problem Set 1 Due January 30, 2017

Cosumnes River College Principles of Macroeconomics Problem Set 1 Due January 30, 2017 Spring 0 Cosumnes River College Principles of Macroeconomics Problem Se Due Januar 0, 0 Name: Soluions Prof. Dowell Insrucions: Wrie he answers clearl and concisel on hese shees in he spaces provided.

More information

Mass Transfer Coefficients (MTC) and Correlations I

Mass Transfer Coefficients (MTC) and Correlations I Mass Transfer Mass Transfer Coeffiiens (MTC) and Correlaions I 7- Mass Transfer Coeffiiens and Correlaions I Diffusion an be desribed in wo ways:. Deailed physial desripion based on Fik s laws and he diffusion

More information

ASTR415: Problem Set #5

ASTR415: 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 information

Age (x) nx lx. Age (x) nx lx dx qx

Age (x) nx lx. Age (x) nx lx dx qx Life Tables Dynamic (horizonal) cohor= cohor followed hrough ime unil all members have died Saic (verical or curren) = one census period (day, season, ec.); only equivalen o dynamic if populaion does no

More information

SIMULATION STUDY OF STOCHASTIC CHANNEL REDISTRIBUTION

SIMULATION STUDY OF STOCHASTIC CHANNEL REDISTRIBUTION Developmens in Business Simulaion and Experienial Learning, Volume 3, 3 SIMULATIO STUDY OF STOCHASTIC CHAEL REDISTRIBUTIO Yao Dong-Qing Towson Universiy dyao@owson.edu ABSTRACT In his paper, we invesigae

More information

Amit Mehra. Indian School of Business, Hyderabad, INDIA Vijay Mookerjee

Amit Mehra. Indian School of Business, Hyderabad, INDIA Vijay Mookerjee RESEARCH ARTICLE HUMAN CAPITAL DEVELOPMENT FOR PROGRAMMERS USING OPEN SOURCE SOFTWARE Ami Mehra Indian Shool of Business, Hyderabad, INDIA {Ami_Mehra@isb.edu} Vijay Mookerjee Shool of Managemen, Uniersiy

More information

USP. Surplus-Production Models

USP. Surplus-Production Models USP Surplus-Producion Models 2 Overview Purpose of slides: Inroducion o he producion model Overview of differen mehods of fiing Go over some criique of he mehod Source: Haddon 2001, Chaper 10 Hilborn and

More information

Chapter 10 INDUCTANCE Recommended Problems:

Chapter 10 INDUCTANCE Recommended Problems: Chaper 0 NDUCTANCE Recommended Problems: 3,5,7,9,5,6,7,8,9,,,3,6,7,9,3,35,47,48,5,5,69, 7,7. Self nducance Consider he circui shown in he Figure. When he swich is closed, he curren, and so he magneic field,

More information

Bifurcation Analysis of a Stage-Structured Prey-Predator System with Discrete and Continuous Delays

Bifurcation Analysis of a Stage-Structured Prey-Predator System with Discrete and Continuous Delays Applied Mahemaics 4 59-64 hp://dx.doi.org/.46/am..4744 Published Online July (hp://www.scirp.org/ournal/am) Bifurcaion Analysis of a Sage-Srucured Prey-Predaor Sysem wih Discree and Coninuous Delays Shunyi

More information

Some Basic Information about M-S-D Systems

Some 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 information

1 Differential Equation Investigations using Customizable

1 Differential Equation Investigations using Customizable Differenial Equaion Invesigaions using Cusomizable Mahles Rober Decker The Universiy of Harford Absrac. The auhor has developed some plaform independen, freely available, ineracive programs (mahles) for

More information

Chapter 4. Truncation Errors

Chapter 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 information

, where P is the number of bears at time t in years. dt (a) If 0 100, lim Pt. Is the solution curve increasing or decreasing?

, where P is the number of bears at time t in years. dt (a) If 0 100, lim Pt. Is the solution curve increasing or decreasing? CALCULUS BC WORKSHEET 1 ON LOGISTIC GROWTH Work he following on noebook paper. Use your calculaor on 4(b) and 4(c) only. 1. Suppose he populaion of bears in a naional park grows according o he logisic

More information

LabQuest 24. Capacitors

LabQuest 24. Capacitors Capaciors LabQues 24 The charge q on a capacior s plae is proporional o he poenial difference V across he capacior. We express his wih q V = C where C is a proporionaliy consan known as he capaciance.

More information

MA 366 Review - Test # 1

MA 366 Review - Test # 1 MA 366 Review - Tes # 1 Fall 5 () Resuls from Calculus: differeniaion formulas, implici differeniaion, Chain Rule; inegraion formulas, inegraion b pars, parial fracions, oher inegraion echniques. (1) Order

More information

5.1 - Logarithms and Their Properties

5.1 - Logarithms and Their Properties Chaper 5 Logarihmic Funcions 5.1 - Logarihms and Their Properies Suppose ha a populaion grows according o he formula P 10, where P is he colony size a ime, in hours. When will he populaion be 2500? We

More information

Section 2.6 Derivatives of products and quotients

Section 2.6 Derivatives of products and quotients Secion 2.6 Derivaives of proucs an quoiens (3/19/08) Overview: In his secion, we erive formulas for erivaives of funcions ha are consruce by aking proucs an quoiens of oher funcions, an we use hese formulas

More information

mywbut.com Lesson 11 Study of DC transients in R-L-C Circuits

mywbut.com Lesson 11 Study of DC transients in R-L-C Circuits mywbu.om esson Sudy of DC ransiens in R--C Ciruis mywbu.om Objeives Be able o wrie differenial equaion for a d iruis onaining wo sorage elemens in presene of a resisane. To develop a horough undersanding

More information

Vehicle Arrival Models : Headway

Vehicle 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 information

Problem Set 7-7. dv V ln V = kt + C. 20. Assume that df/dt still equals = F RF. df dr = =

Problem Set 7-7. dv V ln V = kt + C. 20. Assume that df/dt still equals = F RF. df dr = = 20. Assume ha df/d sill equals = F + 0.02RF. df dr df/ d F+ 0. 02RF = = 2 dr/ d R 0. 04RF 0. 01R 10 df 11. 2 R= 70 and F = 1 = = 0. 362K dr 31 21. 0 F (70, 30) (70, 1) R 100 Noe ha he slope a (70, 1) is

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 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 information

STATE-SPACE MODELLING. A mass balance across the tank gives:

STATE-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 information

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check

More information

) were both constant and we brought them from under the integral.

) 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 information

Numerical Dispersion

Numerical Dispersion eview of Linear Numerical Sabiliy Numerical Dispersion n he previous lecure, we considered he linear numerical sabiliy of boh advecion and diffusion erms when approimaed wih several spaial and emporal

More information

Sterilization D Values

Sterilization D Values Seriliaion D Values Seriliaion by seam consis of he simple observaion ha baceria die over ime during exposure o hea. They do no all live for a finie period of hea exposure and hen suddenly die a once,

More information

An EM algorithm for maximum likelihood estimation given corrupted observations. E. E. Holmes, National Marine Fisheries Service

An EM algorithm for maximum likelihood estimation given corrupted observations. E. E. Holmes, National Marine Fisheries Service An M algorihm maimum likelihood esimaion given corruped observaions... Holmes Naional Marine Fisheries Service Inroducion M algorihms e likelihood esimaion o cases wih hidden saes such as when observaions

More information

Q.1 Define work and its unit?

Q.1 Define work and its unit? CHP # 6 ORK AND ENERGY Q.1 Define work and is uni? A. ORK I can be define as when we applied a force on a body and he body covers a disance in he direcion of force, hen we say ha work is done. I is a scalar

More information

Biol. 356 Lab 8. Mortality, Recruitment, and Migration Rates

Biol. 356 Lab 8. Mortality, Recruitment, and Migration Rates Biol. 356 Lab 8. Moraliy, Recruimen, and Migraion Raes (modified from Cox, 00, General Ecology Lab Manual, McGraw Hill) Las week we esimaed populaion size hrough several mehods. One assumpion of all hese

More information

PHYS-3301 Lecture 5. Chapter 2. Announcement. Sep. 12, Special Relativity. What about y and z coordinates? (x - direction of motion)

PHYS-3301 Lecture 5. Chapter 2. Announcement. Sep. 12, Special Relativity. What about y and z coordinates? (x - direction of motion) Announemen Course webpage hp://www.phys.u.edu/~slee/33/ Tebook PHYS-33 Leure 5 HW (due 9/4) Chaper, 6, 36, 4, 45, 5, 5, 55, 58 Sep., 7 Chaper Speial Relaiiy. Basi Ideas. Consequenes of Einsein s Posulaes

More information

First Order RC and RL Transient Circuits

First Order RC and RL Transient Circuits Firs Order R and RL Transien ircuis Objecives To inroduce he ransiens phenomena. To analyze sep and naural responses of firs order R circuis. To analyze sep and naural responses of firs order RL circuis.

More information

PET467E-Analysis of Well Pressure Tests/2008 Spring Semester/İTÜ Midterm Examination (Duration 3:00 hours) Solutions

PET467E-Analysis of Well Pressure Tests/2008 Spring Semester/İTÜ Midterm Examination (Duration 3:00 hours) Solutions M. Onur 03.04.008 PET467E-Analysis of Well Pressure Tess/008 Spring Semeser/İTÜ Miderm Examinaion (Duraion 3:00 hours) Soluions Name of he Suden: Insrucions: Before saring he exam, wrie your name clearly

More information

Fishing limits and the Logistic Equation. 1

Fishing limits and the Logistic Equation. 1 Fishing limis and he Logisic Equaion. 1 1. The Logisic Equaion. The logisic equaion is an equaion governing populaion growh for populaions in an environmen wih a limied amoun of resources (for insance,

More information

Single and Double Pendulum Models

Single 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 information

CH.7. PLANE LINEAR ELASTICITY. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

CH.7. PLANE LINEAR ELASTICITY. Continuum Mechanics Course (MMC) - ETSECCPB - UPC CH.7. PLANE LINEAR ELASTICITY Coninuum Mechanics Course (MMC) - ETSECCPB - UPC Overview Plane Linear Elasici Theor Plane Sress Simplifing Hpohesis Srain Field Consiuive Equaion Displacemen Field The Linear

More information

4.1 - Logarithms and Their Properties

4.1 - Logarithms and Their Properties Chaper 4 Logarihmic Funcions 4.1 - Logarihms and Their Properies Wha is a Logarihm? We define he common logarihm funcion, simply he log funcion, wrien log 10 x log x, as follows: If x is a posiive number,

More information

Math 36. Rumbos Spring Solutions to Assignment #6. 1. Suppose the growth of a population is governed by the differential equation.

Math 36. Rumbos Spring Solutions to Assignment #6. 1. Suppose the growth of a population is governed by the differential equation. Mah 36. Rumbos Spring 1 1 Soluions o Assignmen #6 1. Suppose he growh of a populaion is governed by he differenial equaion where k is a posiive consan. d d = k (a Explain why his model predics ha he populaion

More information

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still. Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in

More information

15. Bicycle Wheel. Graph of height y (cm) above the axle against time t (s) over a 6-second interval. 15 bike wheel

15. Bicycle Wheel. Graph of height y (cm) above the axle against time t (s) over a 6-second interval. 15 bike wheel 15. Biccle Wheel The graph We moun a biccle wheel so ha i is free o roae in a verical plane. In fac, wha works easil is o pu an exension on one of he axles, and ge a suden o sand on one side and hold he

More information

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)

Econ107 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 information

Introduction to Physical Oceanography Homework 5 - Solutions

Introduction to Physical Oceanography Homework 5 - Solutions Laure Zanna //5 Inroducion o Phsical Oceanograph Homework 5 - Soluions. Inerial oscillaions wih boom fricion non-selecive scale: The governing equaions for his problem are This ssem can be wrien as where

More information

3. Differential Equations

3. Differential Equations 3. Differenial Equaions 3.. inear Differenial Equaions of Firs rder A firs order differenial equaion is an equaion of he form d() d ( ) = F ( (),) (3.) As noed above, here will in general be a whole la

More information

Dynamic Analysis of Damped Driven Pendulum using Laplace Transform Method

Dynamic Analysis of Damped Driven Pendulum using Laplace Transform Method , ISSN 0974-570X (Online), ISSN 0974-578 (Prin), Vol. 6; Issue No. 3; Year 05, Copyrigh 05 by CESER PUBLICATIONS Dynamic Analysis of Damped Driven Pendulum using Laplace Transform Mehod M.C. Agarana and

More information

1.6. Slopes of Tangents and Instantaneous Rate of Change

1.6. Slopes of Tangents and Instantaneous Rate of Change 1.6 Slopes of Tangens and Insananeous Rae of Change When you hi or kick a ball, he heigh, h, in meres, of he ball can be modelled by he equaion h() 4.9 2 v c. In his equaion, is he ime, in seconds; c represens

More information

Let us start with a two dimensional case. We consider a vector ( x,

Let us start with a two dimensional case. We consider a vector ( x, Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our

More information

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x WEEK-3 Reciaion PHYS 131 Ch. 3: FOC 1, 3, 4, 6, 14. Problems 9, 37, 41 & 71 and Ch. 4: FOC 1, 3, 5, 8. Problems 3, 5 & 16. Feb 8, 018 Ch. 3: FOC 1, 3, 4, 6, 14. 1. (a) The horizonal componen of he projecile

More information

Kalman filtering for maximum likelihood estimation given corrupted observations.

Kalman filtering for maximum likelihood estimation given corrupted observations. alman filering maimum likelihood esimaion given corruped observaions... Holmes Naional Marine isheries Service Inroducion he alman filer is used o eend likelihood esimaion o cases wih hidden saes such

More information

Lecture 5. Time series: ECM. Bernardina Algieri Department Economics, Statistics and Finance

Lecture 5. Time series: ECM. Bernardina Algieri Department Economics, Statistics and Finance Lecure 5 Time series: ECM Bernardina Algieri Deparmen Economics, Saisics and Finance Conens Time Series Modelling Coinegraion Error Correcion Model Two Seps, Engle-Granger procedure Error Correcion Model

More information

Jang-Ting Guo Lecture 1-1. Introduction and Some Basics. The building blocks of modern macroeconomics are

Jang-Ting Guo Lecture 1-1. Introduction and Some Basics. The building blocks of modern macroeconomics are Jang-Ting Guo Leure - Inroduion and Some Basis The building bloks of modern maroeonomis are () Solow (Neolassial) growh model Opimal (Ramse) growh model Real business le (RBC) model () Overlapping generaions

More information

Chapter 12: Velocity, acceleration, and forces

Chapter 12: Velocity, acceleration, and forces To Feel a Force Chaper Spring, Chaper : A. Saes of moion For moion on or near he surface of he earh, i is naural o measure moion wih respec o objecs fixed o he earh. The 4 hr. roaion of he earh has a measurable

More information

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

More information

( ) ( ) ( ) + W aa. ( t) n t ( ) ( )

( ) ( ) ( ) + W aa. ( t) n t ( ) ( ) Supplemenary maerial o: Princeon Universiy Press, all righs reserved From: Chaper 3: Deriving Classic Models in Ecology and Evoluionary Biology A Biologis s Guide o Mahemaical Modeling in Ecology and Evoluion

More information

Economics 8105 Macroeconomic Theory Recitation 6

Economics 8105 Macroeconomic Theory Recitation 6 Economics 8105 Macroeconomic Theory Reciaion 6 Conor Ryan Ocober 11h, 2016 Ouline: Opimal Taxaion wih Governmen Invesmen 1 Governmen Expendiure in Producion In hese noes we will examine a model in which

More information

Object tracking: Using HMMs to estimate the geographical location of fish

Object tracking: Using HMMs to estimate the geographical location of fish Objec racking: Using HMMs o esimae he geographical locaion of fish 02433 - Hidden Markov Models Marin Wæver Pedersen, Henrik Madsen Course week 13 MWP, compiled June 8, 2011 Objecive: Locae fish from agging

More information

Nonlinear spring-mass system

Nonlinear spring-mass system Basic phenomenology of simple nonlinear vibraion (free and forced) Manoj Srinivasan (6) Nonlinear spring-mass sysem Damper Spring Mass Hardening Sofening O No damping graviy g lengh l A mass m Do.8.6..

More information

CHEMICAL KINETICS: 1. Rate Order Rate law Rate constant Half-life Temperature Dependence

CHEMICAL KINETICS: 1. Rate Order Rate law Rate constant Half-life Temperature Dependence CHEMICL KINETICS: Rae Order Rae law Rae consan Half-life Temperaure Dependence Chemical Reacions Kineics Chemical ineics is he sudy of ime dependence of he change in he concenraion of reacans and producs.

More information

The average rate of change between two points on a function is d t

The average rate of change between two points on a function is d t SM Dae: Secion: Objecive: The average rae of change beween wo poins on a funcion is d. For example, if he funcion ( ) represens he disance in miles ha a car has raveled afer hours, hen finding he slope

More information

Mathematical Foundations -1- Choice over Time. Choice over time. A. The model 2. B. Analysis of period 1 and period 2 3

Mathematical Foundations -1- Choice over Time. Choice over time. A. The model 2. B. Analysis of period 1 and period 2 3 Mahemaial Foundaions -- Choie over Time Choie over ime A. The model B. Analysis of period and period 3 C. Analysis of period and period + 6 D. The wealh equaion 0 E. The soluion for large T 5 F. Fuure

More information

18 Biological models with discrete time

18 Biological models with discrete time 8 Biological models wih discree ime The mos imporan applicaions, however, may be pedagogical. The elegan body of mahemaical heory peraining o linear sysems (Fourier analysis, orhogonal funcions, and so

More information

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes

23.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 information

Number of modes per unit volume of the cavity per unit frequency interval is given by: Mode Density, N

Number of modes per unit volume of the cavity per unit frequency interval is given by: Mode Density, N SMES404 - LASER PHYSCS (LECTURE 5 on /07/07) Number of modes per uni volume of he aviy per uni frequeny inerval is given by: 8 Mode Densiy, N (.) Therefore, energy densiy (per uni freq. inerval); U 8h

More information

Optimal Transform: The Karhunen-Loeve Transform (KLT)

Optimal Transform: The Karhunen-Loeve Transform (KLT) Opimal ransform: he Karhunen-Loeve ransform (KL) Reall: We are ineresed in uniary ransforms beause of heir nie properies: energy onservaion, energy ompaion, deorrelaion oivaion: τ (D ransform; assume separable)

More information

Outline. Math Elementary Differential Equations. Solution of Linear Growth and Decay Models. Example: Linear Decay Model

Outline. Math Elementary Differential Equations. Solution of Linear Growth and Decay Models. Example: Linear Decay Model Mah 337 - Elemenary Differenial Equaions Joseph M. Mahaffy, jmahaffy@sdsu.edu Deparmen of Mahemaics and Saisics Dynamical Sysems Group Compuaional Sciences Research Cener San Diego Sae Universiy San Diego,

More information

Problem Set 9 Due December, 7

Problem Set 9 Due December, 7 EE226: Random Proesses in Sysems Leurer: Jean C. Walrand Problem Se 9 Due Deember, 7 Fall 6 GSI: Assane Gueye his problem se essenially reviews Convergene and Renewal proesses. No all exerises are o be

More information

MOMENTUM CONSERVATION LAW

MOMENTUM CONSERVATION LAW 1 AAST/AEDT AP PHYSICS B: Impulse and Momenum Le us run an experimen: The ball is moving wih a velociy of V o and a force of F is applied on i for he ime inerval of. As he resul he ball s velociy changes

More information

Section 3.8, Mechanical and Electrical Vibrations

Section 3.8, Mechanical and Electrical Vibrations Secion 3.8, Mechanical and Elecrical Vibraions Mechanical Unis in he U.S. Cusomary and Meric Sysems Disance Mass Time Force g (Earh) Uni U.S. Cusomary MKS Sysem CGS Sysem fee f slugs seconds sec pounds

More information

Module 3: The Damped Oscillator-II Lecture 3: The Damped Oscillator-II

Module 3: The Damped Oscillator-II Lecture 3: The Damped Oscillator-II Module 3: The Damped Oscillaor-II Lecure 3: The Damped Oscillaor-II 3. Over-damped Oscillaions. This refers o he siuaion where β > ω (3.) The wo roos are and α = β + α 2 = β β 2 ω 2 = (3.2) β 2 ω 2 = 2

More information

SPH3U: Projectiles. Recorder: Manager: Speaker:

SPH3U: Projectiles. Recorder: Manager: Speaker: SPH3U: Projeciles Now i s ime o use our new skills o analyze he moion of a golf ball ha was ossed hrough he air. Le s find ou wha is special abou he moion of a projecile. Recorder: Manager: Speaker: 0

More information

Note: For all questions, answer (E) NOTA means none of the above answers is correct.

Note: For all questions, answer (E) NOTA means none of the above answers is correct. Thea Logarihms & Eponens 0 ΜΑΘ Naional Convenion Noe: For all quesions, answer means none of he above answers is correc.. The elemen C 4 has a half life of 70 ears. There is grams of C 4 in a paricular

More information

10.6 Parametric Equations

10.6 Parametric Equations 0_006.qd /8/05 9:05 AM Page 77 Secion 0.6 77 Parameric Equaions 0.6 Parameric Equaions Wha ou should learn Evaluae ses of parameric equaions for given values of he parameer. Skech curves ha are represened

More information

( ) is the stretch factor, and x the

( ) is the stretch factor, and x the (Lecures 7-8) Liddle, Chaper 5 Simple cosmological models (i) Hubble s Law revisied Self-similar srech of he universe All universe models have his characerisic v r ; v = Hr since only his conserves homogeneiy

More information

Basic Circuit Elements Professor J R Lucas November 2001

Basic Circuit Elements Professor J R Lucas November 2001 Basic Circui Elemens - J ucas An elecrical circui is an inerconnecion of circui elemens. These circui elemens can be caegorised ino wo ypes, namely acive and passive elemens. Some Definiions/explanaions

More information

Lawrence Livermore National Laboratory. 1 LLNL-PRES-xxxxxx

Lawrence Livermore National Laboratory. 1 LLNL-PRES-xxxxxx Lawrene Livermore Naional Laboraor 1 LLNL-PRES- RTE as a TSO makes deisions using ools based on modeling of he power ssem simulaion and opimizaion. Reliabili rieria N-1 => Wha if analsis Requiremens: full

More information

κt π = (5) T surrface k BASELINE CASE

κt π = (5) T surrface k BASELINE CASE II. BASELINE CASE PRACICAL CONSIDERAIONS FOR HERMAL SRESSES INDUCED BY SURFACE HEAING James P. Blanhard Universi of Wisonsin Madison 15 Engineering Dr. Madison, WI 5376-169 68-63-391 blanhard@engr.is.edu

More information

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow 1 KEY Mah 4 Miderm I Fall 8 secions 1 and Insrucor: Sco Glasgow Please do NOT wrie on his eam. No credi will be given for such work. Raher wrie in a blue book, or on our own paper, preferabl engineering

More information

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan Ground Rules PC11 Fundamenals of Physics I Lecures 3 and 4 Moion in One Dimension A/Prof Tay Seng Chuan 1 Swich off your handphone and pager Swich off your lapop compuer and keep i No alking while lecure

More information

ME425/525: Advanced Topics in Building Science

ME425/525: Advanced Topics in Building Science ME425/525: Advanced Topics in Building Science Indoor environmenal qualiy for susainable buildings: Lecure 6 Dr. Ellio T. Gall, Ph.D. Lecure 6 Today s objecives o Error propagaion Apply o SS soluion (venilaion,

More information

CS376 Computer Vision Lecture 6: Optical Flow

CS376 Computer Vision Lecture 6: Optical Flow CS376 Compuer Vision Lecure 6: Opical Flow Qiing Huang Feb. 11 h 2019 Slides Credi: Krisen Grauman and Sebasian Thrun, Michael Black, Marc Pollefeys Opical Flow mage racking 3D compuaion mage sequence

More information

GINI MEAN DIFFERENCE AND EWMA CHARTS. Muhammad Riaz, Department of Statistics, Quaid-e-Azam University Islamabad,

GINI MEAN DIFFERENCE AND EWMA CHARTS. Muhammad Riaz, Department of Statistics, Quaid-e-Azam University Islamabad, GINI MEAN DIFFEENCE AND EWMA CHATS Muhammad iaz, Deparmen of Saisics, Quaid-e-Azam Universiy Islamabad, Pakisan. E-Mail: riaz76qau@yahoo.com Saddam Akbar Abbasi, Deparmen of Saisics, Quaid-e-Azam Universiy

More information

Linear Quadratic Regulator (LQR) - State Feedback Design

Linear Quadratic Regulator (LQR) - State Feedback Design Linear Quadrai Regulaor (LQR) - Sae Feedbak Design A sysem is expressed in sae variable form as x = Ax + Bu n m wih x( ) R, u( ) R and he iniial ondiion x() = x A he sabilizaion problem using sae variable

More information

Chapter 15 Lasers, Laser Spectroscopy, and Photochemistry

Chapter 15 Lasers, Laser Spectroscopy, and Photochemistry Chaper 15 Lasers, Laser Specroscopy, and Phoochemisry ackground: In his chaper we will alk abou ligh amplificaion by simulaed emission of radiaion (LASER), heir impac on specroscopy and ligh-iniiaed reacions

More information

Vector autoregression VAR. Case 1

Vector autoregression VAR. Case 1 Vecor auoregression VAR So far we have focused mosl on models where deends onl on as. More generall we migh wan o consider oin models ha involve more han one variable. There are wo reasons: Firs, we migh

More information

References are appeared in the last slide. Last update: (1393/08/19)

References are appeared in the last slide. Last update: (1393/08/19) SYSEM IDEIFICAIO Ali Karimpour Associae Professor Ferdowsi Universi of Mashhad References are appeared in he las slide. Las updae: 0..204 393/08/9 Lecure 5 lecure 5 Parameer Esimaion Mehods opics o be

More information

Y 0.4Y 0.45Y Y to a proper ARMA specification.

Y 0.4Y 0.45Y Y to a proper ARMA specification. HG Jan 04 ECON 50 Exercises II - 0 Feb 04 (wih answers Exercise. Read secion 8 in lecure noes 3 (LN3 on he common facor problem in ARMA-processes. Consider he following process Y 0.4Y 0.45Y 0.5 ( where

More information

3 at MAC 1140 TEST 3 NOTES. 5.1 and 5.2. Exponential Functions. Form I: P is the y-intercept. (0, P) When a > 1: a = growth factor = 1 + growth rate

3 at MAC 1140 TEST 3 NOTES. 5.1 and 5.2. Exponential Functions. Form I: P is the y-intercept. (0, P) When a > 1: a = growth factor = 1 + growth rate 1 5.1 and 5. Eponenial Funcions Form I: Y Pa, a 1, a > 0 P is he y-inercep. (0, P) When a > 1: a = growh facor = 1 + growh rae The equaion can be wrien as The larger a is, he seeper he graph is. Y P( 1

More information

INVERSE RESPONSE COMPENSATION BY ESTIMATING PARAMETERS OF A PROCESS COMPRISING OF TWO FIRST ORDER SYSTEMS

INVERSE RESPONSE COMPENSATION BY ESTIMATING PARAMETERS OF A PROCESS COMPRISING OF TWO FIRST ORDER SYSTEMS Inernaional Journal of Informaion Technology and nowledge Managemen July-December 0, Volume 5, No., pp. 433-438 INVERSE RESPONSE COMPENSATION BY ESTIMATING PARAMETERS OF A PROCESS COMPRISING OF TWO FIRST

More information

Second-Order Differential Equations

Second-Order Differential Equations WWW Problems and Soluions 3.1 Chaper 3 Second-Order Differenial Equaions Secion 3.1 Springs: Linear and Nonlinear Models www m Problem 3. (NonlinearSprings). A bod of mass m is aached o a wall b means

More information

Module 4: Time Response of discrete time systems Lecture Note 2

Module 4: Time Response of discrete time systems Lecture Note 2 Module 4: Time Response of discree ime sysems Lecure Noe 2 1 Prooype second order sysem The sudy of a second order sysem is imporan because many higher order sysem can be approimaed by a second order model

More information