المعادلات التفاضلية · مخطط الطور

مخطط الطور

يمثّل مخطط الطور كل حل لنظام من المعادلات التفاضلية بمنحنى في مستوى متغيراته. في النظام الخطي ⁦x′ = Ax⁩، يحدّد أثر ⁦A⁩ ومحدّدها القيمتين الذاتيتين. وما لم تتطابقا، فإنهما تحدّدان نوع الرسم.

مخططات الطور في هذا التصوّر

يمثّل مخطط الطور حلول نظام من المعادلات التفاضلية بمنحنيات في مستوى متغيراته، يمرّ أحدها بكل نقطة بداية. تعرض هذه الصفحة النظام الخطي x′=Axx' = Ax لمصفوفة AA من الرتبة 2×22 \times 2. بتحريك النقطة (τ,Δ)(\tau, \Delta) على خريطة الأثر والمحدّد يتنقل الرسم بين السروج والعقد والبؤر والمراكز، ويُظهر مستوى القيم الذاتية المجاور القيمتين اللتين تحدّدان كل رسم.

النظام وأثره ومحدّده
(xy)′=(abcd)(xy),τ=a+d,Δ=ad−bc\begin{pmatrix} x \\ y \end{pmatrix}' = \begin{pmatrix} a & b \\ c & d \end{pmatrix} \begin{pmatrix} x \\ y \end{pmatrix}, \qquad \tau = a + d, \quad \Delta = ad - bc

يتبع كل منحنى حقل الاتجاهات. أسهمه متساوية الطول، لذا فهي تبيّن الاتجاه فقط. أما السرعة فتظهر في العلامات الموضوعة بفواصل زمنية متساوية، إذ تتقارب حيث تكون الحركة بطيئة. بعض الحلول تسير على مستقيم يمرّ بالأصل. إذا كان Av=λvAv = \lambda v، فإن x(t)=eλtvx(t) = e^{\lambda t} v حل، لأن مشتقته λeλtv\lambda e^{\lambda t} v تساوي تمامًا نتيجة تطبيق AA عليه. إذن هذه المستقيمات هي الاتجاهات الذاتية للمصفوفة AA، ويمثّل λ\lambda معدل التغير على كل منها: تكون الحركة نحو الخارج عندما λ>0\lambda > 0، ونحو الداخل عندما λ<0\lambda < 0، وتنعدم عندما λ=0\lambda = 0.

تحقق القيم الذاتية المعادلة λ2−τλ+Δ=0\lambda^2 - \tau\lambda + \Delta = 0، ولذلك فإن τ\tau هو مجموعهما، وΔ\Delta حاصل ضربهما، وτ2−4Δ\tau^2 - 4\Delta هو المميّز. كل حدّ على الخريطة يقابل حدثًا في مستوى القيم الذاتية. أسفل Δ=0\Delta = 0 تكون القيمتان مختلفتَي الإشارة ويظهر سرج. وعلى هذا الخط تكون إحداهما 0، فيبقى مستقيم كامل من النقاط ساكنًا. على القطع المكافئ τ2=4Δ\tau^2 = 4\Delta تلتقي القيمتان الذاتيتان، وأعلاه تصبحان زوجًا مركّبًا. عندما Δ>0\Delta > 0، تفصل إشارة τ\tau بين الاضمحلال والنمو، وعلى τ=0\tau = 0 يقع الزوج على المحور التخيلي. وهكذا يحدّد τ\tau وΔ\Delta القيمتين الذاتيتين ونوع الرسم، باستثناء حالة القيمة الذاتية المكررة. عندئذ يعتمد النوع أيضًا على ما إذا كانت AA مضاعفًا عدديًا لمصفوفة الوحدة.

عندما تكون القيمتان الذاتيتان المركّبتان α±βi\alpha \pm \beta i مع β>0\beta > 0، يدور كل حل غير صفري حول الأصل مرة واحدة في الزمن 2π/β2\pi/\beta، وخلاله يتغير حجمه بمعامل e2πα/βe^{2\pi\alpha/\beta}. لذلك تقع العلامات التي تفصل بينها دورة واحدة على شعاع واحد من الأصل، ويربط كل علامتين متتاليتين معامل التحجيم نفسه. عندما α=0\alpha = 0، يكون المعامل 1 وينغلق كل مدار غير صفري ليشكّل قطعًا ناقصًا: إنه مركز. مع Δ=1\Delta = 1 وA=(01−1τ)A = \begin{pmatrix} 0 & 1 \\ -1 & \tau \end{pmatrix}، يكون النظام هو المهتز المخمّد x′′+2ζx′+x=0x'' + 2\zeta x' + x = 0 حيث ζ=−τ/2\zeta = -\tau/2.

عند عبور القطع المكافئ، ينطبق المستقيمان أحدهما على الآخر كمقصّ يُغلق، ثم يختفيان حين يصبح الرسم لولبيًا. بعيدًا عن الرأس، لا يكون للمصفوفة النموذجية على القطع المكافئ سوى هذا المستقيم الواحد. وكل حل آخر يدخل الأصل أو يخرج منه مماسًا له. عند (τ,Δ)=(−2,1)(\tau, \Delta) = (-2, 1)، تكون للمصفوفة −I-I القيم الذاتية نفسها، لكن يوجد فيها حل مستقيم في كل اتجاه. تغيير الأساس، A↦PAP−1A \mapsto PAP^{-1}، يغيّر شكل الرسم ويحافظ على τ\tau وΔ\Delta والقيم الذاتية والنوع. أما المضاعف العددي لمصفوفة الوحدة فيبقى كما هو.

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

  • مستقيما العقدة ليسا متعامدين عادةً. الحل الذي لا يقع على أي منهما يدخل الأصل أو يغادره بمحاذاة الأبطأ، ثم يميل نحو الأسرع بعيدًا عن الأصل.
  • المنحنى الأزرق هو الحل الكامل المارّ بنقطة بداية. والمسار الأصفر هو الجزء الموافق للزمن المنقضي منذ وضعها، لذا بعد تغيير AA يعرض الحل الجديد، وليس سجلًا للحل القديم.

التحكم

  • اضغط على موضع فارغ في مستوى الطور لإضافة نقطة بداية. إذا ضغطت قرب مستقيم، توضع النقطة عليه.

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

للمزيد من القراءة: ويكيبيديا: Phase portrait (الإنجليزية).