المعادلات التفاضلية · معادلات لوتكا–فولتيرا

معادلات لوتكا–فولتيرا

تمثّل معادلات لوتكا–فولتيرا جماعة من الفرائس وأخرى من المفترسات: تتكاثر الفرائس من تلقاء نفسها وتُؤكل، وتنمو أعداد المفترسات بأكلها وتموت من دونها. ترتفع أعداد الجماعتين وتنخفض في دورات حول نقطة توازن.

معادلات لوتكا–فولتيرا في هذا التصوّر

هنا xx عدد الفرائس وyy عدد المفترسات، وكلاهما بالآلاف. يتغير كل عدد بمعدل لكل فرد: تنمو الفرائس بمعدل α\alpha وتُؤكل بمعدل يتناسب مع عدد المفترسات، وتموت المفترسات بمعدل γ\gamma وتنمو بمعدل يتناسب مع عدد الفرائس. يتتبّع مستوى الطور النقطة (x,y)(x, y) من البداية التي تختارها، ويعرض الرسم الزمني أسفله العددين سنة بعد سنة.

النموذج دون صيد أو حد أعلى
x′=x(α−βy),y′=y(δx−γ)x' = x(\alpha - \beta y), \qquad y' = y(\delta x - \gamma)

معدل الفرائس α−βy\alpha - \beta y يساوي صفرًا على المستقيم y=α/βy = \alpha/\beta، ومعدل المفترسات δx−γ\delta x - \gamma يساوي صفرًا على المستقيم x=γ/δx = \gamma/\delta. يتقاطع خطّا النمو الصفري المتقطعان عند نقطة التوازن (x∗,y∗)=(γ/δ,α/β)(x^*, y^*) = (\gamma/\delta, \alpha/\beta) ويقسمان المستوى إلى أربع مناطق. أسفل y∗y^* تزداد الفرائس، وأعلاه تتناقص. إلى يمين x∗x^* تزداد المفترسات، وإلى يساره تتناقص. لذا يدور المسار عكس عقارب الساعة. تبلغ المفترسات قمتها عندما تتوقف عن الازدياد، أي عندما تعود الفرائس إلى x∗x^* بعد هبوطها من قمتها. هذه المدة هي التأخر المعروض في القراءة.

الدوران وحده لا يعيد المسار إلى البداية؛ فاللولب يدور أيضًا. ما يجعل المسار مغلقًا هو كمية تبقى ثابتة على طول كل حل: V=δx−γln⁡x+βy−αln⁡y,V′=(δx−γ)(α−βy)+(βy−α)(δx−γ)=0.\begin{aligned} V &= \delta x - \gamma \ln x + \beta y - \alpha \ln y, \\ V' &= (\delta x - \gamma)(\alpha - \beta y) + (\beta y - \alpha)(\delta x - \gamma) = 0. \end{aligned} كل مستوى من VV أعلى من قيمتها الصغرى يشكّل منحنى مغلقًا حول نقطة التوازن، فيتبعه المسار إلى بدايته ويكرر الدورة. مع الحد الأعلى KK، يبطئ الازدحام نمو الفرائس ولا تبقى VV ثابتة. في المثال ذي الحد الأعلى، يلتف المسار نحو نقطة التوازن بدلًا من أن ينغلق.

مهما كبرت التقلبات، يساوي متوسط كل دورة كاملة نقطة التوازن تمامًا. بقسمة معادلة المفترسات على yy نحصل على (ln⁡y)′=δx−γ(\ln y)' = \delta x - \gamma. خلال دورة واحدة تعود yy إلى قيمتها الابتدائية، فيكون التكامل 0، ويكون متوسط xx هو γ/δ=x∗\gamma/\delta = x^*. وبالطريقة نفسها تعطي معادلة الفرائس متوسط yy، وهو α/β=y∗\alpha/\beta = y^*. التقلبات الأكبر تغيّر القيم القصوى والدور، لا المتوسطات.

صيد الجماعتين باستمرار بالمعدل نفسه لكل فرد، ε\varepsilon، يحوّل α\alpha إلى α−ε\alpha - \varepsilon وγ\gamma إلى γ+ε\gamma + \varepsilon. دون الحد الأعلى تنتقل نقطة التوازن إلى x∗=(γ+ε)/δx^* = (\gamma + \varepsilon)/\delta، y∗=(α−ε)/βy^* = (\alpha - \varepsilon)/\beta: فرائس أكثر في المتوسط ومفترسات أقل. متوسط الفرائس هو العدد الذي تتعادل عنده زيادة المفترسات وموتها تمامًا، لذا يرتفع عندما يجعل الصيد المفترسات تموت بمعدل أسرع.

أمور تستحق الملاحظة

  • في المعادلات يمكن لعدد أفراد جماعة أن يقترب كثيرًا من 0، لكنه لا يصل إليه أبدًا. أما جماعة حقيقية بهذا العدد الضئيل فقد تنقرض.
  • المعاملات والوحدات للتوضيح فقط، ولم تُضبط لتلائم أي نوع حقيقي.

التحكم

  • لا يمكن أن تقع البداية على محور: الجماعة التي تبدأ بعدد 0 تبقى عند 0.

مواضيع ذات صلة