نظرية الاستقرار

مخطط استقرار يصنف خرائط بوانكاريه للنظام الخطي المستقلx=أx،{\displaystyle x'=Ax,}تُصنف هذه العناصر إلى مستقرة أو غير مستقرة وفقًا لخصائصها. يزداد الاستقرار عمومًا باتجاه يسار الرسم البياني. [ 1 ] بعض المصادر أو نقاط التفريغ أو العقد هي نقاط توازن .

في الرياضيات ، تتناول نظرية الاستقرار استقرار حلول المعادلات التفاضلية ومسارات الأنظمة الديناميكية في ظل اضطرابات طفيفة في الشروط الابتدائية. معادلة الحرارة ، على سبيل المثال، هي معادلة تفاضلية جزئية مستقرة لأن الاضطرابات الطفيفة في البيانات الأولية تؤدي إلى تغيرات طفيفة في درجة الحرارة لاحقًا نتيجة لمبدأ القيمة القصوى . في المعادلات التفاضلية الجزئية، يمكن قياس المسافات بين الدوال باستخدام معايير Lp أو معيار القيمة العليا، بينما في الهندسة التفاضلية ، يمكن قياس المسافة بين الفضاءات باستخدام مسافة غروموف-هاوسدورف .

في الأنظمة الديناميكية، يُطلق على المدار اسم "مستقر ليابونوف" إذا كان مساره الأمامي يقع ضمن نطاق صغير كافٍ، أو إذا بقي ضمن نطاق صغير (وربما أكبر). وقد طُوّرت معايير مختلفة لإثبات استقرار المدار أو عدم استقراره. في ظل ظروف مواتية، يمكن اختزال المسألة إلى مشكلة مدروسة جيدًا تتعلق بالقيم الذاتية للمصفوفات . وهناك طريقة أكثر عمومية تعتمد على دوال ليابونوف . عمليًا، يُطبّق أيٌّ من معايير الاستقرار المختلفة.

نظرة عامة في الأنظمة الديناميكية

تتناول أجزاء كثيرة من النظرية النوعية للمعادلات التفاضلية والأنظمة الديناميكية الخصائص التقاربية للحلول والمسارات - أي ما يحدث للنظام بعد فترة زمنية طويلة. أبسط أنواع السلوك هو ما تُظهره نقاط التوازن ، أو النقاط الثابتة، والمدارات الدورية . إذا كان مدار معين مفهومًا جيدًا، فمن الطبيعي أن نتساءل عما إذا كان تغيير طفيف في الشرط الابتدائي سيؤدي إلى سلوك مماثل. تجيب نظرية الاستقرار على الأسئلة التالية: هل سيبقى مدار قريب من مدار معين إلى أجل غير مسمى؟ هل سيتقارب مع المدار المعطى؟ في الحالة الأولى، يُسمى المدار مستقرًا ؛ وفي الحالة الثانية، يُسمى مستقرًا تقاربيًا ، ويُقال إن المدار المعطى جاذب .

محلول متوازنوهـ{\displaystyle f_{e}}يُطلق على النظام المستقل من المعادلات التفاضلية العادية من الدرجة الأولى اسم:

  • مستقر إذا كان لكل (صغير)ϵ>0{\displaystyle \epsilon >0}يوجددلتا>0{\displaystyle \delta >0}بحيث يكون كل حلو(ت){\displaystyle f(t)}وجود شروط أولية ضمن مسافةدلتا{\displaystyle \delta }أيو(ت0)-وهـ<دلتا{\displaystyle \|f(t_{0})-f_{e}\|<\delta }يبقى جزء من التوازن ضمن مسافةϵ{\displaystyle \epsilon }أيو(ت)-وهـ<ϵ{\displaystyle \|f(t)-f_{e}\|<\epsilon }للجميعتت0{\displaystyle t\geq t_{0}}.
  • مستقر تقاربياً إذا كان مستقراً، بالإضافة إلى وجوددلتا0>0{\displaystyle \delta _{0}>0}بحيث كلماو(ت0)-وهـ<دلتا0\displaystyle \|f(t_{0})-f_{e}\|<\delta _{0}}ثمو(ت)وهـ{\displaystyle f(t)\rightarrow f_{e}}مثلت{\displaystyle t\rightarrow \infty }.

يعني الاستقرار أن المسارات لا تتغير كثيرًا عند حدوث اضطرابات طفيفة. والحالة المعاكسة، حيث يتنافر مدار قريب مع المدار المحدد، تُعدّ أيضًا ذات أهمية. عمومًا، يؤدي تغيير الحالة الابتدائية في بعض الاتجاهات إلى اقتراب المسار تدريجيًا من المسار المحدد، وفي اتجاهات أخرى إلى ابتعاده عنه. وقد توجد أيضًا اتجاهات يكون فيها سلوك المدار المضطرب أكثر تعقيدًا (لا يتقارب ولا يبتعد تمامًا)، وعندها لا توفر نظرية الاستقرار معلومات كافية حول الديناميكيات.

من الأفكار الأساسية في نظرية الاستقرار إمكانية تحليل السلوك النوعي للمدار تحت تأثير الاضطرابات باستخدام خطية النظام بالقرب من المدار. على وجه الخصوص، عند كل نقطة توازن لنظام ديناميكي سلس ذي فضاء طوري ذي n بُعد ، توجد مصفوفة A من الرتبة n × n ، تُحدد قيمها الذاتية سلوك النقاط المجاورة ( نظرية هارتمان-غروبمان ). بتعبير أدق، إذا كانت جميع القيم الذاتية أعدادًا حقيقية سالبة أو أعدادًا مركبة ذات أجزاء حقيقية سالبة، فإن النقطة تُعد نقطة ثابتة جاذبة مستقرة، وتتقارب النقاط المجاورة إليها بمعدل أُسّي ، انظر استقرار ليابونوف والاستقرار الأُسّي . أما إذا لم تكن أي من القيم الذاتية تخيلية بحتة (أو صفرًا)، فإن اتجاهات الجذب والتنافر ترتبط بالفضاءات الذاتية للمصفوفة A ذات القيم الذاتية التي يكون جزءها الحقيقي سالبًا وموجبًا على التوالي. وتُعرف عبارات مماثلة لاضطرابات المدارات الأكثر تعقيدًا.

استقرار النقاط الثابتة في بعدين

رسم تخطيطي لأربعة من أكثر أنواع النقاط الثابتة شيوعًا

الحالة النموذجية هي استقرار نقطة الأصل في ظل المعادلة التفاضلية الخطية المستقلةX˙=أX{\displaystyle {\dot {X}}=AX}أينX=[xy]{\displaystyle X={\begin{bmatrix}x\\y\end{bmatrix}}}وأ{\displaystyle A}هي مصفوفة 2×2.

كنا نقوم أحيانًا بتغيير الأساس بواسطةX=جX{\displaystyle X'=CX}بالنسبة لبعض المصفوفات القابلة للعكسج{\displaystyle C}، مما يعطيX˙=ج-1أجX{\displaystyle {\dot {X}}'=C^{-1}ACX'}نقولج-1أج{\displaystyle C^{-1}AC}يكون "أ{\displaystyle A}"على الأساس الجديد". منذالمحققأ=المحققج-1أج{\displaystyle \det A=\det C^{-1}AC}وtrأ=trج-1أج{\displaystyle \operatorname {tr} A=\operatorname {tr} C^{-1}AC}يمكننا تصنيف استقرار الأصل باستخدامالمحققأ{\displaystyle \det A}وtrأ{\displaystyle \operatorname {tr} A}، مع استخدام تغيير الأساس بحرية.

تصنيف أنواع الاستقرار

لوالمحققأ=0{\displaystyle \det A=0}ثم رتبةأ{\displaystyle A}هو صفر أو واحد.

  • إذا كانت الرتبة صفرًا، فإنأ=0{\displaystyle A=0}ولا يوجد تدفق.
  • إذا كانت الرتبة واحدًا، فـكيرأ{\displaystyle \ker A}وأناأ{\displaystyle \operatorname {im} A}كلاهما أحادي البعد.
    • لوكيرأ=أناأ{\displaystyle \ker A=\operatorname {im} A}ثم دعv{\displaystyle v}فترةكيرأ{\displaystyle \ker A}ودعw{\displaystyle w}أن تكون صورة معكوسة لـv{\displaystyle v}ثم في{v،w}{\displaystyle \{v,w\}}أساس،أ=[0100]{\displaystyle A={\begin{bmatrix}0&1\\0&0\end{bmatrix}}}وبالتالي فإن التدفق عبارة عن قص على طولv{\displaystyle v}الاتجاه. في هذه الحالة،trأ=0{\displaystyle \operatorname {tr} A=0}.
    • لوكيرأأناأ{\displaystyle \ker A\neq \operatorname {im} A}ثم دعv{\displaystyle v}فترةكيرأ{\displaystyle \ker A}ودعw{\displaystyle w}فترةأناأ{\displaystyle \operatorname {im} A}ثم في{v،w}{\displaystyle \{v,w\}}أساس،أ=[000أ]{\displaystyle A={\begin{bmatrix}0&0\\0&a\end{bmatrix}}}لبعض الأعداد الحقيقية غير الصفريةأ{\displaystyle a}.
      • لوtrأ>0{\displaystyle \operatorname {tr} A>0}ثم يصبح غير مستقر، ويتباعد بمعدلأ{\displaystyle a}منكيرأ{\displaystyle \ker A}على طول الترجمات المتوازية لـأناأ{\displaystyle \operatorname {im} A}.
      • لوtrأ<0{\displaystyle \operatorname {tr} A<0}ثم يصبح مستقراً، ويتقارب بمعدلأ{\displaystyle a}لكيرأ{\displaystyle \ker A}على طول الترجمات المتوازية لـأناأ{\displaystyle \operatorname {im} A}.

لوالمحققأ0{\displaystyle \det A\neq 0}، نجد أولاً الشكل الطبيعي لجوردان للمصفوفة، للحصول على أساس{v،w}{\displaystyle \{v,w\}}في أيأ{\displaystyle A}وهو أحد الأشكال الثلاثة الممكنة:

  • [أ00ب]{\displaystyle {\begin{bmatrix}a&0\\0&b\end{bmatrix}}}أينأ،ب0{\displaystyle a,b\neq 0}.
    • لوأ،ب>0{\displaystyle a,b>0}، ثم{4المحققأ-(trأ)2=-(أ-ب)20المحققأ=أب>0{\displaystyle {\begin{cases}4\det A-(\operatorname {tr} A)^{2}=-(ab)^{2}\leq 0\\\det A=ab>0\end{cases}}}الأصل هو مصدر ، بمنحنيات متكاملة من الشكلy=جxب/أ{\displaystyle y=cx^{b/a}}
    • وبالمثل بالنسبة لـأ،ب<0{\displaystyle a,b<0}. الأصل هو حوض .
    • لوأ>0>ب{\displaystyle a>0>b}أوأ<0<ب{\displaystyle a<0<b}، ثمالمحققأ<0{\displaystyle \det A<0}والأصل هو نقطة سرجية . مع منحنيات تكاملية على شكلy=جx-|ب/أ|{\displaystyle y=cx^{-|b/a|}}.
  • [أ10أ]{\displaystyle {\begin{bmatrix}a&1\\0&a\end{bmatrix}}}أينأ0{\displaystyle a\neq 0}ويمكن تبسيط ذلك أكثر عن طريق تغيير الأساس معج=[1/أ001]{\displaystyle C={\begin{bmatrix}1/a&0\\0&1\end{bmatrix}}}وبعد ذلكأ=أ[1101]{\displaystyle A=a{\begin{bmatrix}1&1\\0&1\end{bmatrix}}}يمكننا حل المعادلة بشكل صريح لـX˙=أX{\displaystyle {\dot {X}}=AX}معأ=أ[1101]{\displaystyle A=a{\begin{bmatrix}1&1\\0&1\end{bmatrix}}}الحل هوX(ت)=هـأتX(0){\displaystyle X(t)=e^{At}X(0)}معهـأت=هـأت[1أت01]{\displaystyle e^{At}=e^{at}{\begin{bmatrix}1&at\\0&1\end{bmatrix}}}تُسمى هذه الحالة " العقدة المنحلة ". المنحنيات التكاملية في هذا الأساس هي تمددات مركزية لـx=ylny{\displaystyle x=y\ln y}بالإضافة إلى المحور السيني.
    • لوtrأ>0{\displaystyle \operatorname {tr} A>0}إذا كان المصدر متدهوراً، فهو مصدر متدهور . وإلا فهو مصب متدهور .
    • في كلتا الحالتين،4المحققأ-(trأ)2=0{\displaystyle 4\det A-(\operatorname {tr} A)^{2}=0}
  • أ[كوسθالخطيئةθ-الخطيئةθكوسθ]{\displaystyle a{\begin{bmatrix}\cos \theta &\sin \theta \\-\sin \theta &\cos \theta \end{bmatrix}}}أينأ>0،θ(-π،π]{\displaystyle a>0,\theta \in (-\pi ,\pi ]}في هذه الحالة،4المحققأ-(trأ)2=(2أالخطيئةθ)20{\displaystyle 4\det A-(\operatorname {tr} A)^{2}=(2a\sin \theta )^{2}\geq 0}.
    • لوθ(-π،-π/2)(π/2،π]{\displaystyle \theta \in (-\pi ,-\pi /2)\cup (\pi /2,\pi ]}إذن، هذا حوض حلزوني . في هذه الحالة، {4المحققأ-(trأ)2>0trأ<0{\displaystyle {\begin{cases}4\det A-(\operatorname {tr} A)^{2}>0\\\operatorname {tr} A<0\end{cases}}}الخطوط التكاملية عبارة عن حلزونات لوغاريتمية .
    • لوθ(-π/2،π/2){\displaystyle \theta \in (-\pi /2,\pi /2)}إذن، هذا مصدر حلزوني . في هذه الحالة، {4المحققأ-(trأ)2>0trأ>0{\displaystyle {\begin{cases}4\det A-(\operatorname {tr} A)^{2}>0\\\operatorname {tr} A>0\end{cases}}}الخطوط التكاملية عبارة عن حلزونات لوغاريتمية .
    • لوθ=-π/2،π/2{\displaystyle \theta =-\pi /2,\pi /2}إذن، هذا دوران (" استقرار محايد ") بمعدلأ{\displaystyle a}لا يتحرك باتجاه الأصل ولا بعيدًا عنه. في هذه الحالة، trأ=0{\displaystyle \operatorname {tr} A=0}الخطوط المتكاملة عبارة عن دوائر.

يُظهر مخطط الاستقرار على اليمين ملخص النتائج. في كل حالة، باستثناء حالة4المحققأ-(trأ)2=0{\displaystyle 4\det A-(\operatorname {tr} A)^{2}=0}القيم(trأ،المحققأ){\displaystyle (\operatorname {tr} A,\det A)}يسمح بتصنيف فريد لنوع التدفق.

في الحالة الخاصة لـ4المحققأ-(trأ)2=0{\displaystyle 4\det A-(\operatorname {tr} A)^{2}=0}هناك حالتان لا يمكن التمييز بينهما بواسطة(trأ،المحققأ){\displaystyle (\operatorname {tr} A,\det A)}في كلتا الحالتين،أ{\displaystyle A}لها قيمة ذاتية واحدة فقط، بتعدد جبري 2.

  • إذا كانت القيمة الذاتية تمتلك فضاءً ذاتيًا ثنائي الأبعاد ( تعدد هندسي 2)، فإن النظام يكون عقدة مركزية (تسمى أحيانًا " نجمة " أو " عقدة حرجة ") والتي تكون إما مصدرًا (عندماtrأ>0{\displaystyle \operatorname {tr} A>0}) أو حوض (عندماtrأ<0{\displaystyle \operatorname {tr} A<0}). [ 2 ]
  • إذا كان للنظام فضاء ذاتي أحادي البعد ( تعدد هندسي 1)، فإن النظام يكون عقدة متدهورة (إذاالمحققأ>0{\displaystyle \det A>0}) أو تدفق القص (إذاالمحققأ=0{\displaystyle \det A=0}).

تدفق يحافظ على المساحة

متىtrأ=0{\displaystyle \operatorname {tr} A=0}لديناالمحققهـأت=هـtr(أ)ت=1{\displaystyle \det e^{At}=e^{\operatorname {tr} (A)t}=1}لذا فإن التدفق يحافظ على المساحة. في هذه الحالة، يتم تصنيف نوع التدفق بواسطةالمحققأ{\displaystyle \det A}.

  • لوالمحققأ>0{\displaystyle \det A>0}إذن، فهو دوران ("استقرار محايد") حول نقطة الأصل.
  • لوالمحققأ=0{\displaystyle \det A=0}إذن فهو تدفق قصي.
  • لوالمحققأ<0{\displaystyle \det A<0}إذن، الأصل هو نقطة سرجية.

استقرار النقاط الثابتة

أبسط أنواع المدارات هو النقطة الثابتة، أو حالة التوازن. إذا كان النظام الميكانيكي في حالة توازن مستقر، فإن دفعة صغيرة ستؤدي إلى حركة موضعية، مثل تذبذبات صغيرة كما في حالة البندول . في نظام ذي تخميد ، تكون حالة التوازن المستقر مستقرة تقاربياً. من ناحية أخرى، بالنسبة لحالة توازن غير مستقرة، مثل كرة مستقرة على قمة تل، فإن بعض الدفعات الصغيرة ستؤدي إلى حركة ذات سعة كبيرة قد تتقارب أو لا تتقارب إلى الحالة الأصلية.

توجد اختبارات مفيدة للاستقرار في حالة النظام الخطي. ويمكن في كثير من الأحيان استنتاج استقرار النظام غير الخطي من استقرار عملية تحويله إلى نظام خطي .

خرائط

لتكن f : RR دالة قابلة للتفاضل باستمرار ولها نقطة ثابتة a ، f ( a ) = a . لنعتبر النظام الديناميكي الناتج عن تكرار الدالة f :

xن+1=و(xن)،ن=0،1،2،....{\displaystyle x_{n+1}=f(x_{n}),\quad n=0,1,2,\ldots .}

تكون النقطة الثابتة a مستقرة إذا كانت القيمة المطلقة لمشتقة f عند a أقل من 1، وغير مستقرة إذا كانت أكبر من 1. وذلك لأن الدالة f بالقرب من النقطة a لها تقريب خطي بميل f' ( a ) :

و(x)و(أ)+و(أ)(x-أ).{\displaystyle f(x)\approx f(a)+f'(a)\left(x-a\right).}

هكذا

xن+1=و(xن)و(أ)+و(أ)(xن-أ)=أ+و(أ)(xن-أ){\displaystyle {\begin{aligned}x_{n+1}=f(x_{n})&\approx f(a)+f'(a)\left(x_{n}-a\right)\\&=a+f'(a)\left(x_{n}-a\right)\end{aligned}}}و(أ)xن+1-أxن-أ{\displaystyle \Rightarrow f'(a)\approx {\frac {x_{n+1}-a}{x_{n}-a}}}

وهذا يعني أن المشتقة تقيس معدل اقتراب التكرارات المتتالية من النقطة الثابتة a أو ابتعادها عنها. إذا كانت المشتقة عند a تساوي 1 أو -1 بالضبط، فسنحتاج إلى مزيد من المعلومات لتحديد الاستقرار.

يوجد معيار مماثل للدالة f : R nR n القابلة للتفاضل باستمرار والتي لها نقطة ثابتة a ، ويُعبَّر عنه بدلالة مصفوفة جاكوبي الخاصة بها عند a ، J a ( f ) . إذا كانت جميع القيم الذاتية لـ J أعدادًا حقيقية أو مركبة بقيم مطلقة أقل من 1، فإن a نقطة ثابتة مستقرة؛ أما إذا كانت قيمة مطلقة واحدة منها على الأقل أكبر من 1، فإن a غير مستقرة. وكما هو الحال بالنسبة لـ n = 1، فإن حالة كون أكبر قيمة مطلقة تساوي 1 تتطلب مزيدًا من البحث  ، إذ أن اختبار مصفوفة جاكوبي غير حاسم. وينطبق المعيار نفسه بشكل أعم على التشاكلات التفاضلية لمتشعب أملس .

الأنظمة الخطية المستقلة

يمكن تحليل استقرار النقاط الثابتة لنظام من المعادلات التفاضلية الخطية ذات المعاملات الثابتة من الدرجة الأولى باستخدام القيم الذاتية للمصفوفة المقابلة.

نظام مستقل

x=أx،{\displaystyle x'=Ax,}

حيث x ( t ) ∈ Rⁿ و A مصفوفة من الرتبة n × n ذات عناصر حقيقية، ولها حل ثابت

x(ت)=0.{\displaystyle x(t)=0.}

(بمعنى آخر، نقطة الأصل 0 ∈ R n هي نقطة توازن للنظام الديناميكي المقابل). يكون هذا الحل مستقرًا تقاربيًا عندما t → ∞ ("في المستقبل") إذا وفقط إذا كان Re ( λ ) < 0 لجميع القيم الذاتية λ للمصفوفة A. وبالمثل، يكون مستقرًا تقاربيًا عندما t → −∞ ("في الماضي") إذا وفقط إذا كان Re( λ ) > 0 لجميع القيم الذاتية λ للمصفوفة A. إذا وُجدت قيمة ذاتية λ للمصفوفة A بحيث Re( λ ) > 0، فإن الحل يكون غير مستقر عندما t → ∞ .

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

الأنظمة المستقلة غير الخطية

يمكن في كثير من الأحيان إثبات الاستقرار التقاربي للنقاط الثابتة لنظام غير خطي باستخدام نظرية هارتمان-جروبمان .

لنفترض أن v حقل متجه من الرتبة C1 في R^ n والذي ينعدم عند النقطة p ، أي v ( p ) = 0. عندئذٍ يكون النظام المستقل المقابل

x=v(x){\displaystyle x'=v(x)}

له حل ثابت

x(ت)=ص.{\displaystyle x(t)=p.}

لتكن J <sub>p</sub> ( v ) مصفوفة جاكوبية من الرتبة n × n لحقل المتجهات v عند النقطة p . إذا كانت جميع القيم الذاتية لـ J<sub>p </sub> لها جزء حقيقي سالب تمامًا، فإن الحل يكون مستقرًا تقاربيًا. يمكن اختبار هذا الشرط باستخدام معيار راوث-هرويتز .

دالة ليابونوف للأنظمة الديناميكية العامة

تتمثل إحدى الطرق العامة لإثبات استقرار ليابونوف أو الاستقرار التقاربي لنظام ديناميكي في استخدام دوال ليابونوف .

انظر أيضاً

مراجع

  1. الرياضيات إيغوالد - الجبر الخطي: أنظمة المعادلات التفاضلية الخطية: تحليل الاستقرار الخطي تم الوصول إليه في 10 أكتوبر 2019.
  2. "موقع نود - موسوعة الرياضيات" . encyclopediaofmath.org . مؤرشف من الأصل بتاريخ 26-10-2020 . تم الاطلاع عليه بتاريخ 30-03-2023 .