User Tools

Site Tools


tutorial_15-implementing_using_one_compartment

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
tutorial_15-implementing_using_one_compartment [2020/02/19 18:03]
ecoloco
tutorial_15-implementing_using_one_compartment [2023/04/03 17:38] (current)
dmytroh
Line 3: Line 3:
 ===== 3/3 ===== ===== 3/3 =====
  
-When implementing the Lotka-Volterra model using one [[Compartment|Compartment]] block (one for both the predator and prey), the compartment needs to be vectorized. I.e. the rabbit and fox species needs to be added as modelled materials. Then either one can use the same approach as in the two Compartment solution, using Source/Sink and Transfer blocks as means to model the dynamic behaviour or just explicitly define the differential equation inside the Compartment.+When implementing the Lotka-Volterra model using one [[Compartment|Compartment]] block (one for both the predator and prey), the compartment needs to be vectorized. I.e. the rabbit and fox species needs to be added as modelled materials or created as indexes from Index-list. Then either one can use the same approach as in the two Compartment solution, using Source/Sink and Transfer blocks as means to model the dynamic behaviour or just explicitly define the differential equation inside the Compartment.
  
-|Ecosystem| | | +====Option 1==== 
-        | | | + 
-        | | |+Create parameters according to table below: 
 + 
 +^Parameter          ^Unit      ^Value ^Description                                                     ^ 
 +|a                  |1/year    |0.04  |the natural growth rate of rabbits in the absence of predation  
 +|c                  |1/year    |0.2   |the natural death rate of foxes in the absence of food (rabbits)
 +|b                  |1/(year F)|5.0E-4|the death rate per encounter of rabbits due to predation        | 
 +|e                  |F/R       |0.1   |the efficiency of turning predated rabbits into foxes           | 
 +|Init_population                  |n         |Initial animal population        | 
 + 
 +Values for parameter Initial population: 
 +  * //Init_population[Rabbits]= 1000// 
 +  * //Init_population[Foxes]= 4// 
 + 
 +Create source, sink, compartments and transfer according to picture below: 
 +{{rabbitfoxmodel3.png?600}} 
 + 
 +Equations for Transfers: 
 +  * //AnimalBirth[Rabbits]=Animals·a// 
 +  * //AnimalBirth[Foxes]=Animals[Rabbits]·Animals[Foxes]·e·b// 
 +  * //AnimalDeath[Rabbits]=b·Animals[Rabbits]·Animals[Foxes]// 
 +  * //AnimalDeath[Foxes]=c·Animals// 
 + 
 +**Initial Conditions** for compartment: 
 +  * For //Animals[Rabbits] = Init_population// 
 +  * For //Animals[Foxes] = Init_population// 
 + 
 + 
 + 
 + 
 +====Option 2==== 
 +Create parameters according to table below: 
 + 
 +^Parameter          ^Unit      ^Value ^Description                                                     ^ 
 +|a                  |1/year    |0.04  |the natural growth rate of rabbits in the absence of predation 
 +|c                  |1/year    |0.2   |the natural death rate of foxes in the absence of food (rabbits)| 
 +|b                  |1/(year F)|5.0E-4|the death rate per encounter of rabbits due to predation        | 
 +|e                  |F/R       |0.1   |the efficiency of turning predated rabbits into foxes           | 
 + 
 +Create only one compartment with name //Animals//.  
 + 
 +{{rabbitfoxmodel4.png}} 
 + 
 +**Initial Conditions** for compartment: 
 +  * For //Animals[Rabbits] = 1000// 
 +  * For //Animals[Foxes] = 4// 
 + 
 +**dy/dt** for compartment: 
 +  * Rabbits: //dAnimals/dt=Animals[Rabbits]·a-b·Animals[Rabbits]·Animals[Foxes]// 
 +  * Foxes: //dAnimals/dt=Animals[Rabbits]·Animals[Foxes]·e·b-c·Animals[Foxes]//
  
-//''%%Figure 1 Ecolego matrix of one Compartment approach%%''//\\ 
  
  
Line 15: Line 62:
  
 **Solution to the excersise** **Solution to the excersise**
- +  * {{lotka-volterra.eco}} 
-  * {{tutorial:Lotka-Volterra.eco?linkonly}}+  * old one {{tutorial:Lotka-Volterra.eco?linkonly}}
  
  
tutorial_15-implementing_using_one_compartment.1582131812.txt.gz · Last modified: 2020/02/19 18:03 by ecoloco