استقرار ليابونوف

يمكن مناقشة أنواع مختلفة من الاستقرار لحلول المعادلات التفاضلية أو معادلات الفرق التي تصف الأنظمة الديناميكية . أهم هذه الأنواع هو ذلك المتعلق باستقرار الحلول بالقرب من نقطة التوازن. ويمكن مناقشة هذا النوع من خلال نظرية ألكسندر ليابونوف . ببساطة، إذا كانت الحلول التي تبدأ بالقرب من نقطة التوازن مستقرة، فإن الحلول تكون مستقرة.xهـ{\displaystyle x_{e}}ابقَ قريباًxهـ{\displaystyle x_{e}}إلى الأبد، إذنxهـ{\displaystyle x_{e}}هل هو مستقر وفقًا لمعيار ليابونوف ؟ وبشكل أقوى، إذاxهـ{\displaystyle x_{e}}هل دالة ليابونوف مستقرة؟ وجميع الحلول التي تبدأ بالقرب منxهـ{\displaystyle x_{e}}تتقارب إلىxهـ{\displaystyle x_{e}}، ثمxهـ{\displaystyle x_{e}}يُقال إنها مستقرة تقاربياً (انظر التحليل التقاربي ). يضمن مفهوم الاستقرار الأسي معدل اضمحلال أدنى، أي تقديرًا لسرعة تقارب الحلول. يمكن توسيع فكرة استقرار ليابونوف لتشمل متعددات الأبعاد اللانهائية، حيث تُعرف بالاستقرار البنيوي ، والذي يهتم بسلوك الحلول المختلفة ولكن "المتقاربة" للمعادلات التفاضلية. يطبق استقرار المدخلات إلى الحالة (ISS) مفاهيم ليابونوف على الأنظمة ذات المدخلات.

تاريخ

سُميت نظرية استقرار ليابونوف نسبةً إلى ألكسندر ميخائيلوفيتش ليابونوف ، عالم الرياضيات الروسي الذي دافع عن أطروحته بعنوان "المسألة العامة لاستقرار الحركة" في جامعة خاركوف ( جامعة خاركوف الوطنية الحالية، في. إن. كارازين ) عام 1892. [ 1 ] كان ليابونوف رائدًا في المساعي الناجحة لتطوير منهج شامل لتحليل استقرار الأنظمة الديناميكية غير الخطية، وذلك بالمقارنة مع الطريقة المحلية الشائعة لتخطيطها حول نقاط التوازن. لم يحظَ عمله، الذي نُشر في البداية باللغة الروسية ثم تُرجم إلى الفرنسية، باهتمام يُذكر لسنوات عديدة. وقد سبقت النظرية الرياضية لاستقرار الحركة، التي وضعها ليابونوف، وقت تطبيقها في العلوم والتكنولوجيا. علاوة على ذلك، لم يُطبّق ليابونوف نفسه نظريته في هذا المجال، إذ انصبّ اهتمامه على استقرار كتل الموائع الدوارة وتطبيقاتها الفلكية. لم يكن لديه طلاب دكتوراه يتابعون أبحاثه في مجال الاستقرار، وكان مصيره مأساويًا للغاية بسبب انتحاره عام ١٩١٨. [ ٢ ] لعدة عقود، غابت نظرية الاستقرار عن الأنظار تمامًا. كان عالم الرياضيات والميكانيكا الروسي السوفيتي نيكولاي غوريفيتش تشيتايف، الذي كان يعمل في معهد قازان للطيران في ثلاثينيات القرن العشرين، أول من أدرك عظمة اكتشاف أ.م. ليابونوف. كانت مساهمة تشيتايف في النظرية [ ٣ ] بالغة الأهمية لدرجة أن العديد من علماء الرياضيات والفيزياء والمهندسين يعتبرونه الوريث المباشر لليابونوف والخليفة العلمي التالي له في ابتكار وتطوير النظرية الرياضية للاستقرار.

ازداد الاهتمام بهذا الموضوع بشكلٍ كبير خلال فترة الحرب الباردة ، عندما تبيّن أن ما يُسمى "طريقة ليابونوف الثانية" (انظر أدناه) قابلة للتطبيق على استقرار أنظمة التوجيه الجوي ، والتي عادةً ما تحتوي على خصائص غير خطية قوية لا يمكن معالجتها بطرق أخرى. وقد نُشر عدد كبير من الأبحاث في ذلك الوقت، وما زال يُنشر لاحقًا، في أدبيات التحكم والأنظمة. [ 4 ] [ 5 ] [ 6 ] [ 7 ] [ 8 ] ومؤخرًا، حظي مفهوم مُعامل ليابونوف (المرتبط بطريقة ليابونوف الأولى في مناقشة الاستقرار) باهتمام واسع في سياق نظرية الفوضى . كما طُبقت طرق استقرار ليابونوف لإيجاد حلول التوازن في مسائل تخصيص حركة المرور. [ 9 ]

تعريف الأنظمة الزمنية المستمرة

لنفترض نظامًا ديناميكيًا غير خطي مستقلًا

x˙=و(x(ت))،x(0)=x0{\displaystyle {\dot {x}}=f(x(t)),\;\;\;\;x(0)=x_{0}}،

أينx(ت)دRن{\displaystyle x(t)\in {\mathcal {D}}\subseteq \mathbb {R} ^{n}}يرمز إلى متجه حالة النظام ،د{\displaystyle {\mathcal {D}}}مجموعة مفتوحة تحتوي على نقطة الأصل، وو:دRن{\displaystyle f:{\mathcal {D}}\rightarrow \mathbb {R} ^{n}}هو حقل متجهي متصل علىد{\displaystyle {\mathcal {D}}}. يفترضو{\displaystyle f}يوجد توازن عندxهـ{\displaystyle x_{e}}، لهذا السببو(xهـ)=0{\displaystyle f(x_{e})=0}. ثم:

  1. يُقال إن هذا التوازن مستقر وفقًا لمعيار ليابونوف إذا كان لكلϵ>0{\displaystyle \epsilon >0}يوجددلتا>0{\displaystyle \delta >0}بحيث إذاx(0)-xهـ<دلتا{\displaystyle \|x(0)-x_{e}\|<\delta }ثم لكلت0{\displaystyle t\geq 0}لديناx(ت)-xهـ<ϵ{\displaystyle \|x(t)-x_{e}\|<\epsilon }.
  2. يُقال إن توازن النظام المذكور أعلاه مستقر تقاربياً إذا كان مستقراً وفقاً لمعيار ليابونوف، وكان هناكدلتا>0{\displaystyle \delta >0}بحيث إذاx(0)-xهـ<دلتا{\displaystyle \|x(0)-x_{e}\|<\delta }ثمليمتx(ت)-xهـ=0{\displaystyle \lim _{t\rightarrow \infty }\|x(t)-x_{e}\|=0}.
  3. يُقال إن توازن النظام المذكور أعلاه مستقر أُسّيًا إذا كان مستقرًا تقاربيًا، وُجدتα>0، β>0، دلتا>0{\displaystyle \alpha >0,~\beta >0,~\delta >0}بحيث إذاx(0)-xهـ<دلتا{\displaystyle \|x(0)-x_{e}\|<\delta }ثمx(ت)-xهـαx(0)-xهـهـ-βت{\displaystyle \|x(t)-x_{e}\|\leq \alpha \|x(0)-x_{e}\|e^{-\beta t}}للجميعت0{\displaystyle t\geq 0}.

من الناحية المفاهيمية، فإن معاني المصطلحات المذكورة أعلاه هي كالتالي:

  1. يعني استقرار ليابونوف للتوازن أن الحلول التي تبدأ "قريبة بما فيه الكفاية" من التوازن (ضمن مسافة معينة)دلتا{\displaystyle \delta }(في نطاق مسافة معينة) البقاء "قريبًا بما فيه الكفاية" إلى الأبد (ضمن مسافة معينة)ϵ{\displaystyle \epsilon }منه). لاحظ أن هذا يجب أن يكون صحيحًا لأيϵ{\displaystyle \epsilon }قد يرغب المرء في الاختيار.
  2. يعني الاستقرار التقاربي أن الحلول التي تبدأ قريبة بما فيه الكفاية لا تبقى قريبة بما فيه الكفاية فحسب، بل تتقارب في النهاية إلى حالة التوازن.
  3. يعني الاستقرار الأسي أن الحلول لا تتقارب فحسب، بل تتقارب في الواقع أسرع من أو على الأقل بنفس سرعة معدل معين معروف.αx(0)-xهـهـ-βت{\displaystyle \alpha \|x(0)-x_{e}\|e^{-\beta t}}.

المسارϕ(ت){\displaystyle \phi (t)}جذاب (محلياً) إذا

x(ت)-ϕ(ت)0{\displaystyle \|x(t)-\phi (t)\|\rightarrow 0} مثل ت{\displaystyle t\rightarrow \infty }

لجميع المساراتx(ت){\displaystyle x(t)}التي تبدأ قريبة بما فيه الكفاية منϕ(ت){\displaystyle \phi (t)}وستكون جذابة عالميًا إذا تحققت هذه الخاصية لجميع المسارات.

أي أنه إذا كان x ينتمي إلى داخل مشعبه المستقر ، فإنه يكون مستقرًا تقاربيًا إذا كان جاذبًا ومستقرًا في آنٍ واحد. (توجد أمثلة توضح أن الجاذبية لا تستلزم الاستقرار التقاربي. [ 10 ] [ 11 ] [ 12 ] ومن السهل إنشاء مثل هذه الأمثلة باستخدام الاتصالات المتجانسة .)

إذا كانت مصفوفة جاكوبي للنظام الديناميكي عند نقطة التوازن عبارة عن مصفوفة استقرار (أي إذا كان الجزء الحقيقي من كل قيمة ذاتية سالبًا تمامًا)، فإن نقطة التوازن تكون مستقرة تقاربياً.

نظام الانحرافات

بدلاً من النظر في الاستقرار فقط بالقرب من نقطة التوازن (حل ثابت)x(ت)=xهـ{\displaystyle x(t)=x_{e}}يمكن صياغة تعريفات مماثلة للاستقرار بالقرب من حل عشوائي.x(ت)=ϕ(ت){\displaystyle x(t)=\phi (t)}ومع ذلك، يمكن اختزال الحالة الأكثر عمومية إلى حالة توازن عن طريق تغيير المتغيرات التي تُسمى "نظام الانحرافات". عرّفy=x-ϕ(ت){\displaystyle y=x-\phi (t)}، مع مراعاة المعادلة التفاضلية:

y˙=و(ت،y+ϕ(ت))-ϕ˙(ت)=ز(ت،y){\displaystyle {\dot {y}}=f(t,y+\phi (t))-{\dot {\phi }}(t)=g(t,y)}.

لم يعد هذا نظامًا مستقلاً، ولكنه يمتلك نقطة توازن مضمونة عندy=0{\displaystyle y=0}والتي يكون استقرارها مكافئًا لاستقرار الحل الأصليx(ت)=ϕ(ت){\displaystyle x(t)=\phi (t)}.

طريقة ليابونوف الثانية لتحقيق الاستقرار

اقترح ليابونوف، في عمله الأصلي عام 1892، طريقتين لإثبات الاستقرار . [ 1 ] تعتمد الطريقة الأولى على تطوير الحل في متسلسلة، ثم إثبات تقاربها ضمن حدود معينة. أما الطريقة الثانية، والتي تُعرف الآن بمعيار استقرار ليابونوف أو الطريقة المباشرة، فتستخدم دالة ليابونوف V(x) التي تُشابه دالة الجهد في الديناميكا الكلاسيكية. ويتم تقديمها على النحو التالي لنظامx˙=و(x){\displaystyle {\dot {x}}=f(x)}وجود نقطة توازن عندx=0{\displaystyle x=0}لنفترض دالةV:RنR{\displaystyle V:\mathbb {R} ^{n}\rightarrow \mathbb {R} }بحيث

  • V(x)=0{\displaystyle V(x)=0} إذا وفقط إذاx=0{\displaystyle x=0}
  • V(x)>0{\displaystyle V(x)>0} إذا وفقط إذاx0{\displaystyle x\neq 0}
  • V˙(x)=ددتV(x)=أنا=1نVxأناوأنا(x)=Vو(x)0{\displaystyle {\dot {V}}(x)={\frac {d}{dt}}V(x)=\sum _{i=1}^{n}{\frac {\partial V}{\partial x_{i}}}f_{i}(x)=\nabla V\cdot f(x)\leq 0}لجميع قيمx0{\displaystyle x\neq 0}ملاحظة: من أجل الاستقرار التقاربي،V˙(x)<0{\displaystyle {\dot {V}}(x)<0}لx0{\displaystyle x\neq 0}مطلوب.

عندئذٍ تُسمى V(x) دالة ليابونوف ، ويكون النظام مستقرًا وفقًا لمعيار ليابونوف. (لاحظ أنV(0)=0{\displaystyle V(0)=0}مطلوب؛ وإلا على سبيل المثالV(x)=1/(1+|x|){\displaystyle V(x)=1/(1+|x|)}سيثبت ذلكx˙(ت)=x{\displaystyle {\dot {x}}(t)=x}(مستقر محليًا). يلزم شرط إضافي يُسمى "الملاءمة" أو "عدم التقييد القطري" لاستنتاج الاستقرار العالمي. وينطبق الأمر نفسه على الاستقرار التقاربي العالمي (GAS).

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

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

تعريف لأنظمة الزمن المنفصل

يكاد تعريف أنظمة الزمن المتقطع أن يكون مطابقاً لتعريف أنظمة الزمن المستمر. ويُقدّم التعريف أدناه ذلك باستخدام لغة بديلة شائعة الاستخدام في النصوص الرياضية.

ليكن ( X , d ) فضاءً متريًا ، و f  : XX دالة متصلة . يُقال إن النقطة x في X مستقرة وفقًا لمعيار ليابونوف ، إذا،

ϵ>0 دلتا>0 yX [د(x،y)<دلتانشمال د(ون(x)،ون(y))<ϵ].{\displaystyle \forall \epsilon >0\ \exists \delta >0\ \forall y\in X\ \left[d(x,y)<\delta \Rightarrow \forall n\in \mathbf {N} \ d\left(f^{n}(x),f^{n}(y)\right)<\epsilon \right].}

نقول إن x مستقر تقاربياً إذا كان ينتمي إلى داخل مجموعته المستقرة ، أي إذا،

دلتا>0[د(x،y)<دلتاليمند(ون(x)،ون(y))=0].{\displaystyle \exists \delta >0\left[d(x,y)<\delta \Rightarrow \lim _{n\to \infty }d\left(f^{n}(x),f^{n}(y)\right)=0\right].}

استقرار نماذج فضاء الحالة الخطية

نموذج فضاء الحالة الخطي

x˙=أx{\displaystyle {\dot {\textbf {x}}}=A{\textbf {x}}}،

أينأ{\displaystyle A}تكون المصفوفة منتهية، وتكون مستقرة تقاربياً (في الواقع، مستقرة أُسّياً ) إذا كانت جميع الأجزاء الحقيقية للقيم الذاتية لـأ{\displaystyle A}سالبة. هذا الشرط يعادل الشرط التالي: [ 13 ]

أتيم+مأ{\displaystyle A^{\textsf {T}}M+MA}

تكون سالبة التحديد لبعض المصفوفات موجبة التحديدم=متي{\displaystyle M=M^{\textsf {T}}}(دالة ليابونوف ذات الصلة هيV(x)=xتيمx{\displaystyle V(x)=x^{\textsf {T}}Mx}.)

وبناءً على ذلك، نموذج فضاء الحالة الخطي المتقطع زمنيًا

xت+1=أxت{\displaystyle {\textbf {x}}_{t+1}=A{\textbf {x}}_{t}}

تكون مستقرة تقاربياً (في الواقع، مستقرة أُسّياً) إذا كانت جميع القيم الذاتية لـأ{\displaystyle A}لها معامل أصغر من واحد.

تم تعميم هذا الشرط الأخير ليشمل الأنظمة المبدلة: نظام زمني منفصل مبدل خطي (يحكمه مجموعة من المصفوفات){أ1،...،أم}{\displaystyle \{A_{1},\dots ,A_{m}\}})

xت+1=أأناتxت،أأنات{أ1،...،أم}{\displaystyle {{\textbf {x}}_{t+1}}=A_{i_{t}}{\textbf {x}}_{t},\quad A_{i_{t}}\in \{A_{1},\dots ,A_{m}\}}

تكون مستقرة تقاربياً (في الواقع، مستقرة أُسّياً) إذا كان نصف القطر الطيفي المشترك للمجموعة{أ1،...،أم}{\displaystyle \{A_{1},\dots ,A_{m}\}}أصغر من واحد.

استقرار الأنظمة ذات المدخلات

النظام الذي يحتوي على مدخلات (أو عناصر تحكم) يكون على الشكل التالي:

x˙=و(x،u){\displaystyle {\dot {\textbf {x}}}={\textbf {f}}({\textbf {x}},{\textbf {u}})}

حيث يمكن اعتبار المدخل u(t) (الذي يعتمد عادةً على الزمن) بمثابة دالة تحكم ، أو مدخل خارجي ، أو محفز ، أو اضطراب ، أو إجبار . وقد ثبت [ 14 ] أنه بالقرب من نقطة التوازن التي تكون مستقرة وفقًا لمعيار ليابونوف، يظل النظام مستقرًا في ظل اضطرابات صغيرة. أما بالنسبة لاضطرابات المدخلات الأكبر، فإن دراسة هذه الأنظمة تخضع لنظرية التحكم وتُطبق في هندسة التحكم . بالنسبة للأنظمة ذات المدخلات، يجب تحديد تأثير هذه المدخلات على استقرار النظام. النهجان الرئيسيان لهذا التحليل هما استقرار BIBO ( للأنظمة الخطية ) واستقرار المدخلات بالنسبة للحالة (ISS) ( للأنظمة غير الخطية ).

مثال

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

y¨+y-ε(y˙33-y˙)=0.{\displaystyle {\ddot {y}}+y-\varepsilon \left({\frac {{\dot {y}}^{3}}{3}}-{\dot {y}}\right)=0.}

يترك

x1=y،x2=y˙{\displaystyle x_{1}=y,x_{2}={\dot {y}}}

بحيث يكون النظام المقابل هو

x˙1=x2،x˙2=-x1+ε(x233-x2).{\displaystyle {\begin{aligned}&{\dot {x}}_{1}=x_{2},\\&{\dot {x}}_{2}=-x_{1}+\varepsilon \left({\frac {x_{2}^{3}}{3}}-{x_{2}}\right).\end{aligned}}}

الأصلx1=0، x2=0{\displaystyle x_{1}=0,\ x_{2}=0}هي نقطة التوازن الوحيدة. لنختر كدالة ليابونوف

V=12(x12+x22){\displaystyle V={\frac {1}{2}}\left(x_{1}^{2}+x_{2}^{2}\right)}

وهي موجبة تمامًا بشكل واضح . ومشتقتها هي

V˙=x1x˙1+x2x˙2=x1x2-x1x2+εx243-εx22=εx243-εx22.{\displaystyle {\dot {V}}=x_{1}{\dot {x}}_{1}+x_{2}{\dot {x}}_{2}=x_{1}x_{2}-x_{1}x_{2}+\varepsilon {\frac {x_{2}^{4}}{3}}-\varepsilon {x_{2}^{2}}=\varepsilon {\frac {x_{2}^{4}}{3}}-\varepsilon {x_{2}^{2}}.}

يبدو أنه إذا كانت المعلمةε{\displaystyle \varepsilon }إذا كانت القيمة موجبة، فإن الاستقرار يكون تقاربياً لـx22<3.{\displaystyle x_{2}^{2}<3.}لكن هذا خطأ، لأنV˙{\displaystyle {\dot {V}}}لا يعتمد علىx1{\displaystyle x_{1}}وستكون قيمتها صفرًا في كل مكان علىx1{\displaystyle x_{1}}المحور. التوازن مستقر وفقًا لمعيار ليابونوف ولكنه ليس مستقرًا تقاربيًا.

معضلة باربالات واستقرار الأنظمة المتغيرة مع الزمن

قد يكون من الصعب إيجاد دالة ليابونوف ذات مشتقة سالبة محددة كما هو مطلوب بموجب معيار استقرار ليابونوف، ومع ذلك، فإن الدالةV{\displaystyle V}معV˙{\displaystyle {\dot {V}}}قد يكون ذلك متاحًا فقط في حالة شبه التحديد السالب. في الأنظمة المستقلة، يمكن تطبيق نظرية المجموعة الثابتة لإثبات الاستقرار التقاربي، لكن هذه النظرية لا تنطبق عندما تكون الديناميكيات دالة للزمن. [ 15 ]

بدلاً من ذلك، تسمح مبرهنة باربالات بإجراء تحليل شبيه بتحليل ليابونوف لهذه الأنظمة غير المستقلة. وتستند هذه المبرهنة إلى الملاحظات التالية: بافتراض أن f دالة للزمن فقط:

  • امتلاكو˙(ت)0{\displaystyle {\dot {f}}(t)\to 0}لا يعني ذلك أنو(ت){\displaystyle f(t)}له حد أقصى عندت{\displaystyle t\to \infty }. على سبيل المثال،و(ت)=الخطيئة(ln(ت))،ت>0{\displaystyle f(t)=\sin(\ln(t)),\;t>0}.
  • امتلاكو(ت){\displaystyle f(t)}يقترب من الحد الأقصى معت{\displaystyle t\to \infty }لا يعني ذلك أنو˙(ت)0{\displaystyle {\dot {f}}(t)\to 0}. على سبيل المثال،و(ت)=الخطيئة(ت2)/ت،ت>0{\displaystyle f(t)=\sin \left(t^{2}\right)/t,\;t>0}.
  • امتلاكو(ت){\displaystyle f(t)}ذات حد أدنى ومتناقصة (و˙0{\displaystyle {\dot {f}}\leq 0}يشير ذلك إلى أنها تتقارب إلى حد معين. لكنه لا يحدد ما إذا كانت كذلك أم لا.و˙0{\displaystyle {\dot {f}}\to 0}مثلت{\displaystyle t\to \infty }.

تقول ليمّة باربالات : إذا و(ت){\displaystyle f(t)}له حد نهائي لأنت{\displaystyle t\to \infty }وإذاو˙{\displaystyle {\dot {f}}}تكون متصلة بانتظام (الشرط الكافي للاستمرارية المنتظمة هو أنو¨{\displaystyle {\ddot {f}}}(محدود)، ثمو˙(ت)0{\displaystyle {\dot {f}}(t)\to 0}مثلت{\displaystyle t\to \infty }[ 16 ]

وهناك صيغة بديلة كما يلي: ليكنص[1،){\displaystyle p\in [1,\infty )}وq(1،]{\displaystyle q\in (1,\infty ]}. لوولص(0،){\displaystyle f\in L^{p}(0,\infty )}وو˙لq(0،){\displaystyle {\dot {f}}\in L^{q}(0,\infty )}، ثمو(ت)0{\displaystyle f(t)\to 0}مثلت.{\displaystyle t\to \infty .}[ 17 ]

تكون اللمة صحيحة أيضاً في الحالة المتجهة، وذلك بالصيغة التالية: ليكنو(ت){\displaystyle f(t)}لتكن دالة متصلة بانتظام ذات قيم في فضاء باناخهـ{\displaystyle E}وافترض أن0تو(τ)دτ{\displaystyle \textstyle \int _{0}^{t}f(\tau )\mathrm {d} \tau }له حد نهائي لأنت{\displaystyle t\to \infty }. ثمو(ت)0{\displaystyle f(t)\to 0}مثلت{\displaystyle t\to \infty }[ 18 ]

المثال التالي مأخوذ من الصفحة 125 من كتاب سلوتين ولي بعنوان "التحكم غير الخطي التطبيقي" . [ 15 ]

لنفترض نظامًا غير مستقل

هـ˙=-هـ+زw(ت){\displaystyle {\dot {e}}=-e+g\cdot w(t)}
ز˙=-هـw(ت).{\displaystyle {\dot {g}}=-e\cdot w(t).}

هذا غير مستقل لأن المدخلw{\displaystyle w}هي دالة للزمن. افترض أن المدخلw(ت){\displaystyle w(t)}محدود.

أخذV=هـ2+ز2{\displaystyle V=e^{2}+g^{2}}أعطِV˙=-2هـ20.{\displaystyle {\dot {V}}=-2e^{2}\leq 0.}

هذا يعني أنV(ت)V(0){\displaystyle V(t)\leq V(0)}بحسب الشرطين الأولين، وبالتاليهـ{\displaystyle e}وز{\displaystyle g}محدودة. لكنها لا تقول شيئًا عن تقاربهـ{\displaystyle e}إلى الصفر، كماV˙{\displaystyle {\dot {V}}}هي فقط شبه محددة سلبية (ملاحظة)ز{\displaystyle g}يمكن أن تكون قيمتها غير صفرية عندماV˙{\displaystyle {\dot {V}}}=0) والديناميكيات غير مستقلة.

باستخدام ليمّة باربالات:

V¨=-4هـ(-هـ+زw){\displaystyle {\ddot {V}}=-4e(-e+g\cdot w)}.

هذا محدود لأنهـ{\displaystyle e}،ز{\displaystyle g}وw{\displaystyle w}محدودة. وهذا يعنيV˙0{\displaystyle {\dot {V}}\to 0}مثلت{\displaystyle t\to \infty }وبالتاليهـ0{\displaystyle e\to 0}وهذا يثبت أن الخطأ يتقارب.

استقرار الأنظمة المتغيرة مع الزمن ذات الاضطرابات المتلاشية والمحدودة

لنفترض المعادلة التفاضلية المساعدةv˙(ت)=-q(ت)β(v(ت))+هـ(ت)،{\displaystyle {\dot {v}}(t)=-q(t)\beta (v(t))+e(t),}للجميعتت0{\displaystyle t\geq t_{0}}، مع الدولةvR{\displaystyle v\in \mathbb {R} }والحالة الابتدائيةv(ت0)0{\displaystyle v(t_{0})\geq 0}الوظيفةβج0(R،R){\displaystyle \beta \in C^{0}(\mathbb {R} ,\mathbb {R} )}يزداد بشكل مطرد ويرضيβ(0)=0{\displaystyle \beta (0)=0}الوظائفهـ{\displaystyle e}وq{\displaystyle q}ينتمي إلىج0(R،R+){\displaystyle C^{0}(\mathbb {R} ,\mathbb {R} _{+})}تكمن أهمية المعادلة التفاضلية السابقة في أنه يمكن ربط فئة واسعة من متباينات ليابونوف بها من خلال مبدأ المقارنة.

افترض أن ذلك لجميعتت0{\displaystyle t\geq t_{0}}،q(ت)>0{\displaystyle q(t)>0}ومعت0q(ت)دت={\displaystyle \int _{t_{0}}^{\infty }q(t)\,dt=\infty }وليمتهـ(ت)q(ت)=لR+{}.{\displaystyle \lim _{t\to \infty }{\frac {e(t)}{q(t)}}=L\in \mathbb {R} _{+}\cup \{\infty \}.}تشير الخاصية السابقة إلى اضطراب محدود عندمال>0{\displaystyle L>0}تكون محدودة واضطرابًا متلاشيًا عندمال=0{\displaystyle L=0}.

لكل شرط ابتدائيv(ت0)0{\displaystyle v(t_{0})\geq 0}وكل حلv(ت){\displaystyle v(t)}مع أقصى فترة وجود[ت0،ω){\displaystyle [t_{0},\omega )}، أينت0<ω{\displaystyle t_{0}<\omega \leq \infty }، تنطبق الخصائص التالية: [ 19 ]

  1. v(ت)0{\displaystyle v(t)\geq 0}للجميعت[ت0،ω){\displaystyle t\in [t_{0},\omega )}.
  2. لول[0،){\displaystyle L\in [0,\infty )}ولRأنزهـ{β}{\displaystyle L\in \mathrm {Range} \{\beta \}}، ثم ω={\displaystyle \omega =\infty }،v<{\displaystyle \|v\|_{\infty }<\infty }وليمتv(ت)=β-1(ل){\displaystyle \lim _{t\to \infty }v(t)=\beta ^{-1}(L)}.
  3. لول={\displaystyle L=\infty }،v{\displaystyle v}لا يساوي الصفر بشكل منتظم، وليمsβ(s)={\displaystyle \lim _{s\to \infty }\beta (s)=\infty }، ثم ω={\displaystyle \omega =\infty }وليمتv(ت)={\displaystyle \lim _{t\to \infty }v(t)=\infty }.

وقد تم استخلاص النتائج السابقة أيضًا في الأدبيات ضمن سياقات مختلفة؛ انظر على سبيل المثال [ 20 ] [ 21 ] [ 22 ] [ 23 ] [ 24 ] .

استقرار ليابونوف للأنظمة المتغيرة مع الزمن ذات الاضطرابات غير المحدودة

ضع في اعتبارك النظام التالي  :ζ˙(ت)=ز(ت،ζ(ت))؛تت0،{\displaystyle {\dot {\zeta }}(t)=g{\big (}t,\zeta (t){\big )};t\geq t_{0},}ζ(ت0)=ζ0،{\displaystyle \zeta (t_{0})=\zeta _{0},}أين(ت0،ζ0)R×Rم{\displaystyle (t_{0},\zeta _{0})\in \mathbb {R} \times \mathbb {R} ^{m}}، حلζ(ت){\displaystyle \zeta (t)}فيRم{\displaystyle \mathbb {R} ^{m}}(م{\displaystyle m}(عدد صحيح موجب تمامًا)، ودالة محددة جيدًاز:[ت0،)×RمRم{\displaystyle g:[t_{0},\infty )\times \mathbb {R} ^{m}\rightarrow \mathbb {R} ^{m}}معز(ت،0)=0{\displaystyle g(t,0)=0}،تت0{\displaystyle \forall t\geq t_{0}}.

افترض أن النظام يحقق شروط كاراثيودوري؛ أي التعيينتز(ت،ζ){\displaystyle t\mapsto g(t,\zeta )}محصورة محليًا بشكل أساسي على[ت0،)×Rم{\displaystyle [t_{0},\infty )\times \mathbb {R} ^{m}}، قابلة للقياس لكلζRم{\displaystyle \zeta \in \mathbb {R} ^{m}}وهي مستمرة تقريبًا لكلتت0{\displaystyle t\geq t_{0}}يقبل النظام حلاً محلياً متصلاً تماماً وفقاً لمعادلة كاراثيودوري، وهو حل معرف على فترة قصوى.[ت0،ω){\displaystyle [t_{0},\omega )}.

افترض وجود ثوابتα>0{\displaystyle \alpha >0}،β>0{\displaystyle \beta >0}معα<β{\displaystyle \alpha <\beta }الدوال المتصلة بشكل مطلق محليًا ر1ج0(R،R){\displaystyle r_{1}\in C^{0}(\mathbb {R} ,\mathbb {R} )}،ر2ج0(R،R){\displaystyle r_{2}\in C^{0}(\mathbb {R} ,\mathbb {R} )}ودالة قابلة للقياس وفقًا لمقياس لوبيغح:RR{\displaystyle h:\mathbb {R} \rightarrow \mathbb {R} }بما يفي بما يلي:

(أنا) (-1)α=-1{\displaystyle (-1)^{\alpha }=-1}و(-1)β{\displaystyle (-1)^{\beta }}محدد جيدًا (لسنا بحاجة إلى(-1)β{\displaystyle (-1)^{\beta }}(أن تكون محددة جيدًا إذا تم ضمان إيجابية الحلول).

(ii)ر1(ت)>0{\displaystyle r_{1}(t)>0}،ر2(ت)>0{\displaystyle r_{2}(t)>0}،ت>ت0{\displaystyle \forall t>t_{0}}وح(ت)>0{\displaystyle h(t)>0}بالنسبة للجميع تقريباًت>ت0{\displaystyle t>t_{0}}.

(ثالثاً)ليمتر2(ت)ر1(ت)={\displaystyle \lim _{t\rightarrow \infty }{\frac {r_{2}(t)}{r_{1}(t)}}=\infty }(ومن ثم يترتب على ذلك أن الاضطراب غير محدود، كما يتضح من متباينة ليابونوف الواردة لاحقًا).

(رابعاً) ليمتΛ(ت)=0{\displaystyle \lim _{t\rightarrow \infty }\Lambda (t)=0}، أين Λ(ت):=ر1(ت)ر˙2(ت)-ر˙1(ت)ر2(ت)ح(ت)(ر1(ت))2β-α-1β-α(ر2(ت))β-2α+1β-α، بالنسبة للجميع تقريباً ت>ت0.{\displaystyle \Lambda (t):={\frac {r_{1}(t){\dot {r}}_{2}(t)-{\dot {r}}_{1}(t)r_{2}(t)}{h(t)(r_{1}(t))^{\frac {2\beta -\alpha -1}{\beta -\alpha }}(r_{2}(t))^{\frac {\beta -2\alpha +1}{\beta -\alpha }}}},{\mbox{ for almost all }}t>t_{0}.}

(خامساً) لكل حلζ(ت){\displaystyle \zeta (t)}النظام ذو أقصى فترة وجود[ت0،ω){\displaystyle [t_{0},\omega )}توجد ثوابت موجبةدلتا{\displaystyle \delta }،σ{\displaystyle \sigma }،ج1{\displaystyle c_{1}}،ج2{\displaystyle c_{2}}ودالة ليابونوف Vج1(R×Rم،R+){\displaystyle V\in C^{1}(\mathbb {R} \times \mathbb {R} ^{m},\mathbb {R} _{+})}مُرضٍ (لكل شخص تقريبًا)ت(ت0،ω){\displaystyle t\in (t_{0},\omega )}ذلك يرضيV(ت،ζ(ت))<دلتا{\displaystyle V(t,\zeta (t))<\delta }):

ج1|κ|σV(ت،κ)ج2|κ|σ،تR،κRم،{\displaystyle c_{1}|\kappa |^{\sigma }\leq V(t,\kappa )\leq c_{2}|\kappa |^{\sigma },\forall t\in \mathbb {R} ,\forall \kappa \in \mathbb {R} ^{m},}

V(ت،ζ(ت))ز(ت،ζ(ت))(-ر1(ت)Vα(ت،ζ(ت))+ر2(ت)Vβ(ت،ζ(ت)))ح(ت).{\displaystyle \nabla V(t,\zeta (t))\cdot g(t,\zeta (t))\leq {\big (}-r_{1}(t)V^{\alpha }(t,\zeta (t))+r_{2}(t)V^{\beta }(t,\zeta (t)){\big )}h(t).}

ثم، يوجدج3>0{\displaystyle c_{3}>0}بحيث يكون لكل|ζ0|<ج3{\displaystyle |\zeta _{0}|<c_{3}}، واحد يحصلω={\displaystyle \omega =\infty }و|ζ(ت)|ج2ج1σ|ζ0|،تت0،{\displaystyle |\zeta (t)|\leq {\sqrt[{\sigma }]{\frac {c_{2}}{c_{1}}}}|\zeta _{0}|,\forall t\geq t_{0},}لهذا السبب.ζ=0{\displaystyle \zeta =0}مستقر بشكل منتظم. علاوة على ذلك، فإن نقطة الأصل مستقرة تقاربياً. [ 25 ]

مثال

يمكن تمثيل نموذج نمو السكان المتغير مع الزمن مع تأثير ألي بواسطة [ 25 ]

شمال˙(ت)=R(ت)شمال(ت)(شمال(ت)أ(ت)-1)(1-شمال(ت)ك(ت))،{\displaystyle {\dot {N}}(t)=R(t)N(t)\left({\frac {N(t)}{A(t)}}-1\right)\left(1-{\frac {N(t)}{K(t)}}\right),}

أينتت0{\displaystyle t\geq t_{0}}، ولايةشمال(ت)R{\displaystyle N(t)\in \mathbb {R} }، دالة قابلة للقياس وفقًا لمقياس ليبيغ R:RR{\displaystyle R:\mathbb {R} \rightarrow \mathbb {R} }معR(ت)>0{\displaystyle R(t)>0}بالنسبة للجميع تقريباًت>ت0{\displaystyle t>t_{0}}والدوال المتصلة بشكل مطلق محليًا أج0(R،R){\displaystyle A\in C^{0}(\mathbb {R} ,\mathbb {R} )}و كج0(R،R){\displaystyle K\in C^{0}(\mathbb {R} ,\mathbb {R} )}لهذا السبب. أ(ت)>0{\displaystyle A(t)>0}وك(ت)>0{\displaystyle K(t)>0}للجميعت>ت0{\displaystyle t>t_{0}}.

الجانب الأيمن من المعادلة يكون محليًا ليبشيتز فيشمال{\displaystyle N}وبالتالي يوجد حل فريد بفترة وجود قصوى[ت0،ω){\displaystyle [t_{0},\omega )}الأصلشمال=0{\displaystyle N=0}هي نقطة توازن. نهدف إلى استنباط الشروط التي تجعلشمال=0{\displaystyle N=0}توازن انقراض مستقر بشكل منتظم وتوازن انقراض مستقر تقاربياً. ولتحقيق هذه الغاية، افترض أن

ت0R(ت)دت=،{\displaystyle \int _{t_{0}}^{\infty }R(t)dt=\infty ,}ليمت(1أ(ت)+1ك(ت))=،{\displaystyle \lim _{t\rightarrow \infty }\left({\frac {1}{A(t)}}+{\frac {1}{K(t)}}\right)=\infty ,} و ليمت(أ˙(ت)(ك(ت))2+ك˙(ت)(أ(ت))2R(ت)أ(ت)ك(ت)(أ(ت)+ك(ت)))=0.{\displaystyle \lim _{t\rightarrow \infty }\left({\frac {{\dot {A}}(t)(K(t))^{2}+{\dot {K}}(t)(A(t))^{2}}{R(t)A(t)K(t)(A(t)+K(t))}}\right)=0.}

يتركV=شمال2{\displaystyle V=N^{2}}تتحقق المتباينة معج1=1{\displaystyle c_{1}=1}،ج2=1{\displaystyle c_{2}=1}وσ=2{\displaystyle \sigma =2}. تؤدي الحسابات البسيطة إلى

V˙(ت)2R(ت)(-V(ت)+(1أ(ت)+1ك(ت))V32(ت)) بالنسبة للجميع تقريباً ت>ت0،{\displaystyle {\dot {V}}(t)\leq 2R(t)\left(-V(t)+\left({\frac {1}{A(t)}}+{\frac {1}{K(t)}}\right)V^{\frac {3}{2}}(t)\right){\mbox{ for almost all }}t>t_{0},}

والتي تأخذ شكل المتباينة معα=1{\displaystyle \alpha =1}،β=32{\displaystyle \beta ={\frac {3}{2}}}،ح(ت)=2R(ت){\displaystyle h(t)=2R(t)}،ر1(ت)=1{\displaystyle r_{1}(t)=1}،ر2(ت)=1أ(ت)+1ك(ت){\displaystyle r_{2}(t)={\frac {1}{A(t)}}+{\frac {1}{K(t)}}} ودلتا{\displaystyle \delta }هو أمر اعتباطي في(ت0،){\displaystyle (t_{0},\infty )}علاوة على ذلك، يمكن للمرء بسهولة إظهار الوظيفةΛ(ت){\displaystyle \Lambda (t)}يؤول إلى الصفر عندمات{\displaystyle t}يؤول إلى اللانهاية. وبالتالي، تتحقق جميع الشروط، ومن ثم يوجدج3>0{\displaystyle c_{3}>0}بحيث يكون لكل|شمال0|<ج3{\displaystyle |N_{0}|<c_{3}}، واحد يحصلω={\displaystyle \omega =\infty }و

|شمال(ت)|ج2ج1σ|شمال0|،تت0،{\displaystyle |N(t)|\leq {\sqrt[{\sigma }]{\frac {c_{2}}{c_{1}}}}|N_{0}|,\forall t\geq t_{0},}

وبالتاليشمال=0{\displaystyle N=0}مستقر بشكل موحد. في الواقع،شمال=0{\displaystyle N=0}مستقر بشكل منتظم ومستقر تقاربياً.

انظر أيضاً

مراجع

  1. 1 2 ليابونوف، أ.م. المشكلة العامة لاستقرار الحركة (باللغة الروسية)، أطروحة دكتوراه، جامعة خاركيف 1892. الترجمات الإنجليزية: (1) استقرار الحركة ، دار النشر الأكاديمية، نيويورك ولندن، 1966 (2) المشكلة العامة لاستقرار الحركة ، (ترجمة أ.ت. فولر) تايلور وفرانسيس، لندن 1992. يتضمن الكتاب سيرة ذاتية بقلم سميرنوف وقائمة مراجع شاملة لأعمال ليابونوف.
  2. شيرباكوف 1992 .
  3. ^ Chetaev، NG حول المسارات المستقرة للديناميكيات، Kazan Univ Sci Notes، vol.4 no.1 1936؛ استقرار الحركة، نُشر في الأصل باللغة الروسية عام 1946 بواسطة ОГИЗ. جوس. هناك نظرية تقنية. مضاء، موسكو-لينينغراد. ترجمة مورتون نادلر، أكسفورد، 1961، 200 صفحة.
  4. ^ ليتوف، آم (1955). Устойчивость нелинейный regулируемых систем [ استقرار أنظمة التحكم غير الخطية ] (بالروسية). موسكو: جوستيخيزدات.ترجمة إنجليزية. برينستون 1961
  5. كالمان، ر . إي .؛ بيرترام، ج. ف. (1960). "تحليل وتصميم أنظمة التحكم باستخدام "الطريقة الثانية" لليابونوف: الجزء الأول - أنظمة الزمن المستمر". مجلة الهندسة الأساسية . 82 (2): 371-393 . doi : 10.1115/1.3662604 .
  6. لاسال، جيه بي ؛ ليفشيتز، إس. (1961). الاستقرار باستخدام طريقة ليابونوف الثانية مع التطبيقات . نيويورك: أكاديميك برس.
  7. باركس، بي سي (1962). "طريقة ليابونوف في نظرية التحكم الآلي". التحكم . 1 نوفمبر 1962 2 ديسمبر 1962.
  8. كالمان، ر. إي. (1963). "دوال ليابونوف لمسألة لور في التحكم الآلي" . وقائع الأكاديمية الوطنية للعلوم في الولايات المتحدة الأمريكية . 49 ( 2): 201-205 . Bibcode : 1963PNAS...49..201K . doi : 10.1073/pnas.49.2.201 . PMC 299777. PMID 16591048 .  
  9. سميث، إم جيه؛ ويستن، إم بي (1995). "نموذج مستمر لتخصيص حركة المرور يوميًا ووجود توازن ديناميكي مستمر للمستخدم". حوليات بحوث العمليات . 60 (1): 59-79 . doi : 10.1007/BF02031940 . S2CID 14034490 . 
  10. هان، وولفغانغ (1967). استقرار الحركة . سبرينغر. الصفحات 191-194 ، القسم 40. doi : 10.1007/978-3-642-50085-5 . ISBN  978-3-642-50087-9.
  11. براون، فيليب؛ غرون، لارس؛ كيلت، كريستوفر م. (2021). (عدم) استقرار التضمينات التفاضلية: مفاهيم، ومكافئات، وخصائص شبيهة بليابونوف . سبرينغر. ص 19-20 ، مثال 2.18. doi : 10.1007/978-3-030-76317-6 . ISBN  978-3-030-76316-9. S2CID 237964551 . 
  12. فينوجراد، ر. إي. (1957). "قصور طريقة الأسس المميزة لدراسة المعادلات التفاضلية غير الخطية" . دوكلادي أكاديميي ناوك (باللغة الروسية). 114 (2): 239-240 .
  13. جوه، بي إس (1977). "الاستقرار العالمي في الأنظمة متعددة الأنواع". المجلة الأمريكية لعلم الطبيعة . 111 (977): 135-143 . Bibcode : 1977ANat..111..135G . doi : 10.1086/283144 . S2CID 84826590 . 
  14. مالكين، آي جي. نظرية استقرار الحركة، موسكو 1952 (غوستخيزدات)، الفصل الثاني، الفقرة 4 (بالروسية). الترجمة الإنجليزية، مكتب الخدمات اللغوية، واشنطن AEC -tr-3352؛ الأصل: حول الاستقرار في ظل الاضطرابات المؤثرة باستمرار، الرياضيات التطبيقية 1944، المجلد 8، العدد 3، 241-245 (بالروسية)؛ ​​الجمعية الأمريكية للرياضيات، الترجمة رقم 8
  15. 1 2 سلوتين، جان جاك إي؛ ويبينج لي (1991). تطبيق التحكم غير الخطي . نيوجيرسي: برنتيس هول.
  16. I. Barbălat، أنظمة المعادلات التفاضلية للتذبذبات غير الخطية، القس الرياضيات. تطبيق بيور. 4 (1959) 267–270، ص. 269.
  17. ب. فاركاس وآخرون، اختلافات في ليمّة باربالات، المجلة الأمريكية للرياضيات الشهرية (2016) 128، العدد 8، 825-830، DOI: 10.4169/amer.math.monthly.123.8.825، ص 827.
  18. ب. فاركاس وآخرون، اختلافات في ليمّة باربالات، المجلة الأمريكية للرياضيات الشهرية (2016) 128، العدد 8، 825-830، DOI: 10.4169/amer.math.monthly.123.8.825، ص 826.
  19. ناصر، إم إف إم (2020). تقارب الحالة لفئة من المعادلات التفاضلية المتغيرة مع الزمن . مجلة IMA للتحكم الرياضي والمعلومات، 37 ، 27-38.
  20. ناصر، م.ف.م. وإخوان، ف. (2019). استقرار الأنظمة المتغيرة مع الزمن في غياب دوال ليابونوف الصارمة . مجلة IMA للتحكم الرياضي والمعلومات . DOI: 10.1093/imamci/dnx056.
  21. جيانغ، ز.، لين، ي. ووانغ، ي. (2009). تثبيت الأنظمة غير الخطية المتغيرة مع الزمن: نهج دالة ليابونوف للتحكم . مجلة علوم الأنظمة والتعقيد ، 22 ، 683-696.
  22. ماليسوف، م. ومازينك، ف. (2005). ملاحظات إضافية حول دوال ليابونوف المستقرة تمامًا من المدخلات إلى الحالة للأنظمة المتغيرة مع الزمن . أوتوماتيكا ، 41 ، 1973-1978.
  23. مازينك، ف. (2003). دوال ليابونوف الصارمة للأنظمة المتغيرة مع الزمن . أوتوماتيكا ، 39 ، 349-353.
  24. مو، إكس. وتشنغ، دي. (2005). حول استقرار وتثبيت أنظمة التحكم غير الخطية المتغيرة مع الزمن . المجلة الآسيوية للتحكم ،
  25. 1 2 ناصر، إم إف إم، سلوك حلول الأنظمة غير المستقلة ذات الاضطرابات غير المحدودة بالقرب من اللانهاية الزمنية ، مجلة IMA للتحكم الرياضي والمعلومات، 39، ص 60-79، النظرية 6.2، 2022.

للمزيد من القراءة

تتضمن هذه المقالة مواد من موقع asymptotically stable على PlanetMath ، وهو مرخص بموجب رخصة Creative Commons Attribution/Share-Alike .