로트카-볼테라 모델이란?
로트카-볼테라 방정식은 포식자와 피식자라는 두 종이 상호작용하는 생물학적 시스템의 동역학을 기술하는 두 개의 1차 비선형 미분방정식입니다. 1920년대에 알프레드 로트카(Alfred Lotka)와 비토 볼테라(Vito Volterra)가 각각 독립적으로 도출한 이 모델은 개체수가 주기적으로 늘었다 줄어드는 현상을 잘 보여줍니다. 피식자가 늘어나면 포식자가 먹이를 충분히 얻어 번성하고, 그 결과 피식자가 급감하면 포식자도 굶어 줄어들며, 이 순환이 끊임없이 반복됩니다.
계산기 사용법
초기 피식자 개체수(\(x_0\))와 포식자 개체수(\(y_0\)), 그리고 네 가지 비율 파라미터(\(\alpha\) 피식자 증식률, \(\beta\) 포식률, \(\delta\) 먹이 섭취에 따른 포식자 증가율, \(\gamma\) 포식자 사망률)를 입력하고, 시간 간격(\(dt\))과 적분할 단계 수를 지정하세요. 계산기는 고전적인 4차 룽게-쿠타(RK4) 기법으로 시스템을 단계별로 전진시키며, 최종 개체수와 각 종이 도달한 최대 개체수, 그리고 해석적으로 구한 평형점을 함께 보여줍니다.
공식 풀이
피식자 방정식 \(\frac{dx}{dt} = \alpha x - \beta x y\)는 피식자가 지수적으로 증가하지만 포식자와 마주칠 때마다 줄어든다는 의미입니다. 포식자 방정식 \(\frac{dy}{dt} = \delta x y - \gamma y\)는 사냥에 성공할수록 포식자가 늘고 자연적으로는 사망한다는 뜻입니다.
$$\frac{dx}{dt} = \alpha\,x - \beta\,x\,y \qquad \frac{dy}{dt} = \delta\,x\,y - \gamma\,y$$0이 아닌 평형점은 \(x^{*} = \frac{\gamma}{\delta}\), \(y^{*} = \frac{\alpha}{\beta}\)에 위치하며, 이 지점에서는 두 도함수가 모두 0이 되어 개체수가 일정하게 유지됩니다.
$$x^{*} = \frac{\gamma}{\delta} \qquad y^{*} = \frac{\alpha}{\beta}$$
계산 예시
\(x_0 = 40\), \(y_0 = 9\), \(\alpha = 0.1\), \(\beta = 0.02\), \(\delta = 0.01\), \(\gamma = 0.1\), \(dt = 0.1\)로 1000단계(\(t\)가 100까지)를 진행하면 개체수가 진동합니다. 평형점은 피식자 $$x^{*} = \frac{0.1}{0.01} = 10,$$ 포식자 $$y^{*} = \frac{0.1}{0.02} = 5$$이며, 이 값이 바로 순환이 그 주위를 도는 중심이 됩니다.
자주 묻는 질문
어떤 적분법을 사용하나요? 4차 룽게-쿠타(RK4) 기법을 사용합니다. 같은 시간 간격에서도 단순한 오일러 방식보다 훨씬 정확합니다.
개체수가 폭발하거나 0이 되는 이유는 무엇인가요? 시간 간격이 지나치게 크거나 비율 값이 극단적이면 수치해가 불안정해질 수 있습니다. \(dt\)를 줄이면 더 매끄럽고 정확한 순환을 얻을 수 있습니다.
평형점에서 시작하면 어떻게 되나요? \(x_0 = \frac{\gamma}{\delta}\)이고 \(y_0 = \frac{\alpha}{\beta}\)이면 두 도함수가 모두 0이므로 개체수가 일정하게 유지됩니다. 이것이 바로 시스템의 고정점입니다.
