μμ€ν μ€μ predator-prey equations λΌλκ² μμ.
μ΄κΈ°μ λλμ μ«μκ° μ μλ ν λΌλ€μ΄ μ λμ μ«μκ° λ§ λμ΄λλλ°, ν λΌκ° λμ΄λλ©΄ κ·Έλ§νΌ λλλ λμ΄λκΈ° λλ¬Έμ λ€μ ν λΌμ μ«μκ° μ€μ΄λ€κ² λκ³
ν λΌμ μ«μκ° μ€μ΄λ€μμΌλ λ¨Ήμ΄κ° λΆμ‘±ν λλκ° λ€μ μ€μ΄λλ μμ
rμ΄ ν λΌμ μΈκ΅¬μ΄κ³ , wκ° λλμ μΈκ΅¬λΌλ©΄
dr/dt = a*r - b*r*w
dw/dt = c*r*w - d*w
μ΄λ° κ΄κ³λΌκ³ ν¨. μ μμ dr/dt = 0, dw/dt = 0 μ΄ λλ rκ³Ό wλ₯Ό μ°ΎμΌλ©΄ κ· νμ μ΄λ£° μ μλ μ«μκ° λμ΄.
μ λ°©μ μμ Lotka-Volterra equation μ΄λΌκ³ μ΄λ―Έ μ΄λ¦μ΄ μλ λ―ΈλΆλ°©μ μμ΄μ§λ§... μ λ°©μ μμ λ³΄κ³ ν λΌμ λλ μ¬μ΄μ κ΄κ³λ₯Ό μΈκ³Όκ΄κ³λ‘ λνλ΄λ©΄
μ΄λ΄κ±° κ°μ.
μ¬κΈ°μ κ°μ μ λͺ¨λ ν λΌλ μμ°μ¬ νμ§ μκ³ λλνν μλ λΉνλ μ΄λͺ μ κ°μ§κ³ μκ³ Β μμ°μ¬λ λλλ§ ν μ μμ.
κ·Έλμ 보면 μμ b*r*wλ‘ λμ΄ μλ κ²μ. λλκ° μμΌλ©΄ (w=0) ν λΌλ λ€μ§μ§ μμ.
μ λκΈ λκΈν λ€μ΄μ΄κ·Έλ¨μ Causal diagram μ΄λΌκ³ νλ κ²μΈλ°, μ¬μ€ μμλ Causal diagramμ λ¨Όμ κ·Έλ¦¬κ³ , κ±°κΈ°μ λ―ΈλΆλ°©μ μμ λμ΄λ΄μΌν¨.
μμ λ―ΈλΆλ°©μ μμ λμ΄λ΄λ©΄ κ·Έ λ€μμ μ΄μ°νν΄μ νμ΄μ£Όλ©΄ λ¨.
μλ₯Ό λ€λ©΄ dr/dt = a*r - b*rwλ μλ μ½λμ²λΌ νλ¦°λ€.
new_rabbit = rabbit[i] + dt * (a*rabbit[i] - b*rabbit[i]*wolf[i])
a, b, dtμ ν¬κΈ°μ λ°λΌ λ°μ°νκΈ°λ ν¨. μ μ ν μμ μ°Ύμλ λκ³ , μ λ κ² λ§κ³ RK4 κ°μκ±Έ μ¨λ λκ³ , μ μ΄μ μ΄μ°νλ₯Ό implicitνκ² νλ©΄ νμ μλ ΄νλ€.
μμΈν건 μμΉν΄μ ν μ€νΈλΆμ γ±γ±
μ Causal diagramμ λμμ λ§κ² λ°κΏμ€λ€λ©΄ μ¬μν°μ λͺ¨λΈμ ꡬν ν μλ μκ³ ,
κ·Έλ₯ μ λͺ¨λΈ κ·Έλλ‘ μ°λ©΄ λμ μ£Όλ³μ μ΄μ λλ¬Όμ μμ λ°λΌΒ λμ μΈκ·Ό λͺ¬μ€ν°μ μ«μκ° μ€λ₯΄λ½ λ΄λ¦¬λ½νλκ±Έ λΆμν μλ μμλ―.
μ΄λ°κ±Έ feedback loop λΌκ³ ν¨.
μ°λ° 1000λ§λ¦¬ λΉ λ΄λΉμΌμ° 40λ§λ¦¬λ©΄ κ· νμ μ΄λ£¬ λνλ―Όκ΅μ κ³ λ±νκ΅λ₯Ό ꡬνν μ μλ€.
κ²°λ‘ λ΄λΌ μΉ γ γ - dc App
μ°λ°μ λ΄λΉμΌμ§μ μ¬μ΄λ κ·Έκ²λ³΄λ€λ μ€λ¬ν κ΄κ³λΌ μ κ±°λ‘ μΌλ°ννκΈ°λ λΆμ‘±ν¨
μ°λ°μ μκ° μ€λ©΄ λλμ νμκ°μ²΄κ° μ°λ°λ‘ νΈμ λμ΄μΌν¨. λ΄λΉμΌμ§μ΄λΌλ κ°λ μμ²΄κ° λ‘컬μ μΈ λΆλΆμ μν₯μ λΌμΉλ©° λ΄λΆμ κ΄κ³μμ μλ©μ μμκ΄κ³μ λ°λ₯Έ μ¬λ¦¬νμ λ°μ μ νμνλ 컨μ μμ μ λν΄μ μ‘°μ§λλ λ¨Ήμ΄μ¬μ¬λ±μ μκ³ λ¦¬μ¦μΌλ‘ μΌλ°ννκΈ° μ’ μ΄λ ΅λ€κ³ μκ°ν¨
κΈ°μμ체λ ν₯λ―Έλ‘μ°λ μ΄λ‘ μ΄ λ°μ νλκ²μ λ³΄κ³ μΆμ
γ γ γ κ΄μ¬κ°μ Έμ£Όλ μ λ λ μ°Ύμλ³΄κ² μ - dc App
λκΈλ€λ΄λΌ μ°λ°ν λ Όλ¬Έ μ°κ³ μμλ€ γ γ γ γ