15.08.2015 Views

Introduction to the Modeling and Analysis of Complex Systems

introduction-to-the-modeling-and-analysis-of-complex-systems-sayama-pdf

introduction-to-the-modeling-and-analysis-of-complex-systems-sayama-pdf

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

4.5. BUILDING YOUR OWN MODEL EQUATION 53• Exponential growth• Convergence <strong>to</strong> a certain population limitWe should check <strong>the</strong> first one first. The original model already shows exponential growthby itself, so this is already done. So we move on <strong>to</strong> <strong>the</strong> second one. Obviously, <strong>the</strong> originalmodel doesn’t show such convergence, so this is what we will need <strong>to</strong> implement in <strong>the</strong>model.The third tip says you need <strong>to</strong> focus on a specific component <strong>to</strong> be revised. There aremany options here. You could revise a, x t−1 , or you could even add ano<strong>the</strong>r term <strong>to</strong> <strong>the</strong>right h<strong>and</strong> side. But in this particular case, <strong>the</strong> convergence <strong>to</strong> a certain limit means that<strong>the</strong> growth ratio a should go <strong>to</strong> 1 (i.e., no net growth). So, we can focus on <strong>the</strong> a part, <strong>and</strong>replace it by an unknown function <strong>of</strong> <strong>the</strong> population size f(x t−1 ). The model equation nowlooks like this:x t = f(x t−1 )x t−1 (4.20)Now your task just got simpler: just <strong>to</strong> design function f(x). Think about constraints it has<strong>to</strong> satisfy. f(x) should be close <strong>to</strong> <strong>the</strong> original constant a when <strong>the</strong> population is small,i.e., when <strong>the</strong>re are enough environmental resources, <strong>to</strong> show exponential growth. In <strong>the</strong>meantime, f(x) should approach 1 when <strong>the</strong> population approaches a carrying capacity <strong>of</strong><strong>the</strong> environment (let’s call it K for now). Ma<strong>the</strong>matically speaking, <strong>the</strong>se constraints meanthat <strong>the</strong> function f(x) needs <strong>to</strong> go through <strong>the</strong> following two points: (x, f(x)) = (0, a) <strong>and</strong>(K, 1).And this is where <strong>the</strong> fourth tip comes in. If you have no additional information aboutwhat <strong>the</strong> model should look like, you should choose <strong>the</strong> simplest possible form that satisfies<strong>the</strong> requirements. In this particular case, a straight line that connects <strong>the</strong> two pointsabove is <strong>the</strong> simplest one (Fig. 4.5), which is given byf(x) = − a − 1 x + a. (4.21)KYou can plug this form in<strong>to</strong> <strong>the</strong> original equation, <strong>to</strong> complete a new ma<strong>the</strong>matical equation:(x t = − a − 1 )Kx t−1 + a x t−1 (4.22)Now it seems your model building is complete. Following <strong>the</strong> fifth tip, let’s check if<strong>the</strong> new model behaves <strong>the</strong> way you intended. As <strong>the</strong> tip suggests, testing with extremevalues <strong>of</strong>ten helps find out possible issues in <strong>the</strong> model. What happens when x t−1 = 0? In

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!