This preview has intentionally blurred sections. Sign up to view the full version.
View Full Document
Unformatted text preview: Lab 5: Nonlinear Systems Goals In this lab you will use the pplane8 program to study two nonlinear systems by direct numerical simulation. The first model, from population biology, displays interesting nonlinear oscillations (so called limit cycles ). They are isolated closed orbits like the red circular one in Figure 1. Figure 1: Limit cycle If one starts on the red circle, one keeps going around on it counterclockwise. The limit cycle in the picture is an attractor , that is, if the initial point is sufficiently close to the circle then its orbit will approach the circle. The second model is a system whose solutions depend on a parameter. Neither of these systems is described by exactly solvable systems of differential equations. Although much may be learned from strictly theoretical analyses, we must ultimately rely on computational methods to extract quantitative predictions from these systems. 1 Predatorprey interactions In class we considered a model of predatorprey species interactions known as the Lotka Volterra model (referred to in Section 6.3 as the predatorprey system ). If x describes the size of a population of rabbits and y describes a population of foxes (which feed on the population of rabbits) then the LotkaVolterra model of their interactions says that there are positive constants a, b, c, d so that dx dt = x ( a by ) , (1) dy dt = y ( c + dx ) . (2) Unfortunately, this model predicts some unlikely behavior. In the absence of foxes ( y = 0), equation (1) becomes dx/dt = ax . In other words, without any foxes the rabbits will always grow exponentially without bound. And even if the predator population is small, they will always eat the prey at a rate proportional to their product. In other words, 10 foxes surrounded by 100,000 rabbits would each have to eat ten times more than 10 foxes surrounded by 10,000 rabbits. If the rabbit population could be held at a fixed level x > c/d , equation (2) becomes dy/dt = Cy where C = c + dx > 0. In other words, if the rabbit population is maintained at a given level, above some threshold, the fox population will always grow exponentially without bound. None of these predictions are ecologically reasonable. The following model addresses these problems. For positive values of r , a more reasonable model of the two populations is the system dx dt = x (1 x ) 5 xy 5 x + 1 , (3) dy dt = ry 1 y x . (4) In the absence of predators, the prey satisfies the logistic equation with equilibrium population x = 1. In the presence of predators, prey is consumed at a rate 5 xy/ (5 x + 1)....
View
Full
Document
This note was uploaded on 02/08/2012 for the course MATH 216 taught by Professor Gavinlarose during the Winter '02 term at University of MichiganDearborn.
 Winter '02
 gavinlarose
 Differential Equations, Equations, Linear Systems

Click to edit the document details