Mathematics · الكتاب 5 · Bachelor Year 3

الرياضيات الجامعية — السنة 3

الرياضيات الجامعية — السنة 3 · Bachelor Year 3

19المعادلات التفاضلية العادية

حلّت السنة الجامعية 2 المعادلات التفاضلية الخطية وصاغت مبرهنة كوشي–ليبشيتز؛ ويبرهن هذا الفصل عليها — مرتين: الوجود والوحدانية بنقطة باناخ الصامدة، والبنية الشاملة بنظرية الحلول الأعظمية ومبرهنة الهروب من المتراصات. ثم يُعاد بناء النظرية الخطية على أسس أمينة (المُحلِّل، والفرونسكياني، والأُسّي المصفوفي، ودوهاميل)، ويفتح النصف الثاني من الفصل النظرية النوعية — الانسيابات، ومواضع التوازن، ودوال ليابونوف، والاستقرار بالخطّية: أي كيف نفهم حلولًا لن نحسبها أبدًا. والنواس، في مسألة نهاية الأسبوع، هو دراسة الحالة الأبدية. وفي كل ما يلي، تكون UR×RdU \subseteq \R\times\R^d مفتوحة وتكون f ⁣:URdf \colon U \to \R^d متصلة؛ والحل للمعادلة x=f(t,x)x' = f(t, x) هو تطبيق x ⁣:IRdx \colon I \to \R^d من الصنف C1\mathcal C^1 (حيث II فترة) منحنيه البياني في UU ويحقق المعادلة.

19.1 كوشي–ليبشيتز

تعريف 19.1

تكون ff ليبشيتزية محليًّا في xx إذا كان لكل نقطة من UU جوارٌ VV وثابتٌ LL يحققان f(t,x1)f(t,x2)Lx1x2\norm{f(t, x_1) - f(t, x_2)} \leq L\norm{x_1 - x_2} من أجل (t,x1),(t,x2)V(t, x_1), (t, x_2) \in V. وإذا كانت ff من الصنف C1\mathcal C^1 (أو مجرّد وجود xf\partial_xf واتصالها)، كانت ليبشيتزية محليًّا في xx: أي على جوار محدّب متراص، بمتراجحة القيم المتوسطة مع L=supxfL = \sup\vertiii{\partial_xf}.

مبرهنة 19.2 (كوشي–ليبشيتز، الصيغة المحلية)

لتكن ff متصلة وليبشيتزية محليًّا في xx، ولتكن (t0,x0)U(t_0, x_0) \in U. يوجد T>0T > 0 بحيث يكون لمسألة كوشي

x=f(t,x),x(t0)=x0x' = f(t, x), \qquad x(t_0) = x_0

حلٌّ واحد بالضبط على [t0T,t0+T][t_0 - T, t_0 + T].

برهان. نختار a,b>0a, b > 0 بحيث Q=[t0a,t0+a]×Bˉ(x0,b)UQ = [t_0 - a, t_0 + a]\times\bar B(x_0, b) \subseteq U، ويكون عليه fM\norm f \leq M وتكون ff ليبشيتزية بالثابت LL في xx. والدالة من الصنف C1\mathcal C^1 حلٌّ إذا وفقط إذا حققت المعادلة التكاملية

x(t)=x0+t0tf(s,x(s)) ⁣dsx(t) = x_0 + \int_{t_0}^{t}f\bigl(s, x(s)\bigr)\,\dd s

(بالمبرهنة الأساسية في التفاضل والتكامل، بالاتجاهين). ولتكن T=min(a,bM,12L)T = \min\bigl(a, \frac bM, \frac1{2L}\bigr) و I=[t0T,t0+T]I = [t_0 - T, t_0 + T] و

E={xC(I,Rd):x(t)x0b علىI},\mathcal E = \{x \in \mathcal C(I, \R^d) : \norm{x(t) - x_0} \leq b\ \text{على} I\},

وهي جزء مغلق من فضاء باناخ (C(I,Rd),)(\mathcal C(I, \R^d), \norm\cdot_\infty): أي تام (التعريف 7.1). ونعرّف Φ(x)(t)=x0+t0tf(s,x(s)) ⁣ds\Phi(x)(t) = x_0 + \int_{t_0}^tf(s, x(s))\dd s: فمن أجل xEx \in \mathcal E، Φ(x)(t)x0Mtt0MTb\norm{\Phi(x)(t) - x_0} \leq M\abs{t - t_0} \leq MT \leq b — ومنه يرسل Φ\Phi المجموعةَ E\mathcal E في نفسها — ومن أجل x,yEx, y \in \mathcal E:

Φ(x)(t)Φ(y)(t)t0tLx(s)y(s) ⁣dsLTxy12xy:\norm{\Phi(x)(t) - \Phi(y)(t)} \leq \Bigl|\int_{t_0}^t L\,\norm{x(s) - y(s)}\,\dd s\Bigr| \leq LT\,\norm{x - y}_\infty \leq \tfrac12\norm{x - y}_\infty :

أي تقلّص. وتعطي نقطة باناخ الصامدة (المبرهنة 7.4) نقطةً صامدة وحيدة في E\mathcal E: أي الوجود، والوحدانية بين الحلول التي تبقى في Bˉ(x0,b)\bar B(x_0, b) — لكن كل حل على II يبقى هناك (x(t)x0Mtt0b\norm{x(t) - x_0} \leq M\abs{t - t_0} \leq b ما دام المنحني البياني في QQ، بحجة اتصال): أي الوحدانية على II.

مبرهنة مساعدة 19.3 (غرونوال)

لتكن u ⁣:I[0,)u \colon I \to \intco0\infty متصلة، ولتكن t0It_0 \in I، ولنفترض أن

u(t)α+βt0tu(s) ⁣ds(tI)u(t) \leq \alpha + \beta\,\Bigl|\int_{t_0}^{t}u(s)\,\dd s\Bigr| \qquad (t \in I)

مع α0\alpha \geq 0 و β>0\beta > 0. عندئذٍ u(t)αeβtt0u(t) \leq \alpha\,\eu^{\beta\abs{t - t_0}} على II.

برهان. من أجل tt0t \geq t_0: لتكن v(t)=α+βt0tu(s) ⁣dsv(t) = \alpha + \beta\int_{t_0}^tu(s)\dd s، ومنه uvu \leq v و v=βuβvv' = \beta u \leq \beta v و (veβ(tt0))0(v\eu^{-\beta(t - t_0)})' \leq 0: أي v(t)v(t0)eβ(tt0)=αeβ(tt0)v(t) \leq v(t_0)\eu^{\beta(t-t_0)} = \alpha\eu^{\beta(t - t_0)}. ومن أجل tt0t \leq t_0 نطبّق الشيء نفسه على u~(t)=u(2t0t)\tilde u(t) = u(2t_0 - t).

نتيجة 19.4 (الوحدانية والتعلّق الاتصالي)

تحت فرضيات المبرهنة 19.2، يتوافق حلّان للمعادلة x=f(t,x)x' = f(t,x) يتوافقان عند نقطة واحدة على فترة تعريفهما المشتركة. وكمّيًّا، إذا كان x,yx, y حلّين منحنياهما البيانيان في منطقة تكون فيها ff ليبشيتزية بالثابت LL في xx، فإن

x(t)y(t)x(t0)y(t0)eLtt0.\norm{x(t) - y(t)} \leq \norm{x(t_0) - y(t_0)}\,\eu^{L\abs{t - t_0}} .

برهان. التقدير: يحقق u=xyu = \norm{x - y} أن u(t)u(t0)+Lt0tuu(t) \leq u(t_0) + L\abs{\int_{t_0}^tu} (بطرح المعادلتين التكامليتين)؛ ثم غرونوال. وأما الوحدانية الشاملة: فمجموعة التوافق {t:x(t)=y(t)}\{t : x(t) = y(t)\} مغلقة في الفترة المشتركة، وغير خالية، ومفتوحة — إذ حول أي نقطة توافق، نغطّي قطعةً متراصة من المنحني البياني المشترك بعدد منتهٍ من صناديق ليبشيتز ونطبّق التقدير مع u(t1)=0u(t_1) = 0 على كلٍّ منها: فيكون xyx \equiv y محليًّا. وجزءٌ غير خالٍ مفتوح ومغلق من فترة هو الفترة كلها.

19.2 الحلول الأعظمية

مبرهنة 19.5 (الحلول الأعظمية؛ الهروب من المتراصات)

لنفترض أن ff متصلة وليبشيتزية محليًّا في xx.

  1. لكل مسألة كوشي حلٌّ أعظمي وحيد x ⁣:(T,T+)Rdx \colon \intoo{T_-}{T_+} \to \R^d: وكل حل آخر يمرّ بالنقطة (t0,x0)(t_0, x_0) هو تضييقٌ له. والفترة مفتوحة.
  2. (الهروب) من أجل كل متراصة KUK \subseteq U يوجد ε>0\varepsilon > 0 بحيث (t,x(t))K(t, x(t)) \notin K من أجل كل t(T+ε,T+)t \in \intoo{T_+ - \varepsilon}{T_+} (وبالتناظر عند TT_-): أي إن المنحني البياني لحل أعظمي يغادر في النهاية كل جزء متراص من UU. وعلى وجه الخصوص، من أجل U=R×RdU = \R\times\R^d و T+<+T_+ < +\infty: x(t)+\norm{x(t)} \to +\infty حين tT+t \to T_+^- (أي الانفجار).

برهان. (1) لتكن S\mathcal S مجموعةَ جميع الحلول المارّة بالنقطة (t0,x0)(t_0, x_0)؛ فحسب النتيجة 19.4 يتوافق أي اثنين منها على تقاطع فترتيهما، ومنه تلتحم: فعلى J=ySIyJ = \bigcup_{y \in \mathcal S}I_y، نعرّف x(t)=y(t)x(t) = y(t) من أجل أي yy معرَّف عند tt: وهو حلٌّ معرَّف جيدًا، وأعظمي ووحيد بداهةً. و JJ مفتوح: إذ حلٌّ معرَّف عند طرف كان سيُمدَّد بواسطة المبرهنة 19.2 عند ذلك الطرف.

(2) لنفترض أن الادعاء يخفق عند T+T_+: فتوجد tnT+t_n \to T_+ تحقق (tn,x(tn))K(t_n, x(t_n)) \in K؛ ولاحظ أن هذا يفرض T+<T_+ < \infty، أو، إذا كان T+=T_+ = \infty، فلا شيء نبرهن عليه (KK محدودة زمنيًّا). فليكن T+<T_+ < \infty. وأما التراص: فالثوابت المنتظمة M,L,a,bM, L, a, b تصلح من أجل جميع معطيات كوشي في جوار للمتراصة KK — وعمليًّا، نغطي KK بعدد منتهٍ من الصناديق QiQ_i كما في برهان المبرهنة المحلية وليكن T>0T^* > 0 أصغرَ أزمنة الوجود الموافقة: فكل معطى كوشي في KK يطلق حلًّا يعيش TT^* على الأقل بعد زمنه الابتدائي. وبتطبيق ذلك عند (tn,x(tn))(t_n, x(t_n)) مع tn>T+T/2t_n > T_+ - T^*/2 يمتدّ xx وراء T+T_+ (ويتوافق التمديد مع xx بالوحدانية، ثم يمدّده): وهو تناقض مع الأعظمية. ومنه يغادر المنحني البياني KK نهائيًّا قبل T+T_+. ومن أجل U=R×RdU = \R\times\R^d: إذا كان x(t)↛\norm{x(t)}\not\to\infty، أبقت متتالية tnT+t_n \to T_+ النقاطَ (tn,x(tn))(t_n, x(t_n)) في المتراصة [t0,T+]×Bˉ(0,R)[t_0, T_+]\times\bar B(0, R): وهو مستبعد.

نتيجة 19.6 (الوجود الشامل تحت النمو الخطي)

إذا كان U=I×RdU = I\times\R^d (II فترة مفتوحة) و f(t,x)α(t)x+β(t)\norm{f(t, x)} \leq \alpha(t)\norm x + \beta(t) مع α,β\alpha, \beta متصلتين، فإن كل حل أعظمي معرَّف على II كله.

برهان. على متراصة [t0,T]I[t_0, T] \subseteq I: x(t)x0+t0t(αx+β)\norm{x(t)} \leq \norm{x_0} + \int_{t_0}^t(\alpha\norm x + \beta)، ومنه بغرونوال (مع α,β\alpha, \beta محدودتين بالعددين A,BA, B هناك) x(t)(x0+B(Tt0))eA(Tt0)\norm{x(t)} \leq (\norm{x_0} + B(T - t_0))\eu^{A(T - t_0)}: أي محدود. وإذا كان T+<supIT_+ < \sup I، بقي المنحني البياني في متراصة من I×RdI\times\R^d بجوار T+T_+: وهو يناقض الهروب (المبرهنة 19.5).

19.3 الجمل الخطية

في هذا المقطع تكون A ⁣:IMd(R)A \colon I \to M_d(\R) و b ⁣:IRdb \colon I \to \R^d متصلتين؛ والجملة هي x=A(t)x+b(t)x' = A(t)x + b(t) — أي نموّ خطي: فجميع الحلول الأعظمية تعيش على II كله (النتيجة 19.6).

مبرهنة 19.7 (البنية)

تشكّل حلول الجملة المتجانسة x=A(t)xx' = A(t)x فضاءً متجهيًّا SHS_H بُعده dd؛ ومن أجل كل t0t_0، يكون التقييم xx(t0)x \mapsto x(t_0) تماثلًا SHRdS_H \to \R^d. والمُحلِّل R(t,s)GLd(R)R(t, s) \in GL_d(\R)، المعرَّف بالشرط: أن tR(t,s)vt \mapsto R(t, s)v هو الحل الذي قيمته vv عند ss، يحقق

R(s,s)=I,R(t,u)R(u,s)=R(t,s),tR(t,s)=A(t)R(t,s),R(s,s) = I,\quad R(t, u)R(u, s) = R(t, s),\quad \partial_tR(t,s) = A(t)R(t,s),

وتُحلّ المسألة غير المتجانسة بواسطة صيغة دوهاميل:

x(t)=R(t,t0)x0+t0tR(t,s)b(s) ⁣ds.x(t) = R(t, t_0)\,x_0 + \int_{t_0}^{t}R(t, s)\,b(s)\,\dd s .

وأخيرًا يخضع الفرونسكياني w(t)=detR(t,s)w(t) = \det R(t, s) لصيغة ليوفيل w=trA(t)ww' = \operatorname{tr}A(t)\,w، ومنه w(t)=exp(sttrA)>0w(t) = \exp\bigl(\int_s^t\operatorname{tr}A\bigr) > 0.

برهان. تجعل خطية المعادلة الحلولَ فضاءً متجهيًّا؛ والتقييم خطي، ومتباين (بالوحدانية: إذ حلٌّ ينعدم عند t0t_0 هو 0\equiv 0)، وغامر (بالوجود): فالبُعد dd. وخصائص المُحلِّل تعيد صياغة الوحدانية (إذ يحل طرفا كل متطابقة مسألةَ كوشي نفسها)؛ والقابلية للقلب من R(s,t)R(t,s)=IR(s,t)R(t,s) = I. وأما دوهاميل: فنشتقّ الصيغة — x(t)=A(t)R(t,t0)x0+R(t,t)b(t)+t0tA(t)R(t,s)b(s) ⁣ds=A(t)x(t)+b(t)x'(t) = A(t)R(t,t_0)x_0 + R(t,t)b(t) + \int_{t_0}^tA(t)R(t,s)b(s)\dd s = A(t)x(t) + b(t) (والاشتقاق تحت علامة التكامل مشروع: إذ المقدار المكامَل من الصنف C1\mathcal C^1 في tt بمشتقة متصلة في (t,s)(t,s)؛ أو نتحقق بالمعادلة التكاملية). وأما ليوفيل: w(t+h)=det(R(t+h,t))w(t)w(t + h) = \det\bigl(R(t+h, t)\bigr)w(t) و R(t+h,t)=I+hA(t)+o(h)R(t + h, t) = I + hA(t) + o(h) (من المعادلة التكاملية)، ومنه det=1+htrA(t)+o(h)\det = 1 + h\operatorname{tr}A(t) + o(h) (بنشر det\det عند II): أي w(t)=trA(t)w(t)w'(t) = \operatorname{tr}A(t)\,w(t)؛ ثم نكامل المعادلة التفاضلية الخطية السلّمية.

مبرهنة 19.8 (الأُسّي المصفوفي)

من أجل AMd(C)A \in M_d(\C)، تتقارب المتسلسلة eA=n0Ann!\eu^{A} = \sum_{n\geq0}\frac{A^n}{n!} (تقاربًا مطلقًا، بأي معيار شبه ضربي)، ويكون eA+B=eAeB\eu^{A+B} = \eu^A\eu^B كلما كان AB=BAAB = BA، ويكون tetAt \mapsto \eu^{tA} مُحلِّل الجملة الثابتة: R(t,s)=e(ts)AR(t,s) = \eu^{(t-s)A}؛ وهو من الصنف C\mathcal C^\infty مع  ⁣d ⁣dtetA=AetA\frac{\dd}{\dd t}\eu^{tA} = A\eu^{tA}. وفوق ذلك: إذا كان Reλ<α<0\operatorname{Re}\lambda < -\alpha < 0 من أجل كل قيمة ذاتية λ\lambda للمصفوفة AA، فإن etACeαt\vertiii{\eu^{tA}} \leq C\eu^{-\alpha t} من أجل t0t \geq 0.

برهان. التقارب: An/n!An/n!\vertiii{A^n/n!} \leq \vertiii A^n/n!، وهي قابلة للجمع (التمرين 7.1(b) في جبر باناخ MdM_d). ومن أجل A,BA, B متبادلتين: يُعاد ترتيب جداء كوشي للمتسلسلتين المتقاربتين تقاربًا مطلقًا، عبر مبرهنة ذات الحدين (الصحيحة حين AB=BAAB = BA)، إلى n(A+B)nn!\sum_n\frac{(A+B)^n}{n!}. وأما القابلية للاشتقاق فمباشرةً: e(t+h)AetAh=etAehAIhetAA\frac{\eu^{(t+h)A} - \eu^{tA}}h = \eu^{tA}\frac{\eu^{hA} - I}{h} \to \eu^{tA}A لأن ehAIhAn2hn1Ann!=O(h)\norm{\frac{\eu^{hA} - I}h - A} \leq \sum_{n\geq2}\frac{\abs h^{n-1}\vertiii A^n}{n!} = O(h). ومنه يحل te(ts)Avt \mapsto \eu^{(t - s)A}v مسألةَ كوشي المعرِّفة للمقدار R(t,s)vR(t,s)v. وأما الحدّ الطيفي: فبصيغة جوردان (المبرهنة 3.18A=P(D+N)P1A = P(D + N)P^{-1} حيث DD قطرية تحمل القيم الذاتية، و NN معدومة القوة، و DN=NDDN = ND. عندئذٍ etA=PetDetNP1\eu^{tA} = P\,\eu^{tD}\eu^{tN}P^{-1} مع etDe(α+δ)t\vertiii{\eu^{tD}} \leq \eu^{-(\alpha + \delta)t} من أجل δ>0\delta > 0 ما (حيث t0t \geq 0)، ويكون etN\eu^{tN} كثير حدود في tt (بانعدام القوة): فيكون الجداء Ceαt\leq C\eu^{-\alpha t} (إذ يغلب eδt\eu^{-\delta t} كثيرَ الحدود).

مثال 19.9 (المستوي، مصنَّفًا)

من أجل x=Axx' = Ax مع AM2(R)A \in M_2(\R) قابلة للقلب، تتحدّد صورة الطور بجوار 00 بالعددين τ=trA\tau = \operatorname{tr}A و δ=detA\delta = \det A، عبر القيمتين الذاتيتين λ±=τ±τ24δ2\lambda_\pm = \frac{\tau \pm \sqrt{\tau^2 - 4\delta}}2:

  • δ<0\delta < 0: قيمتان ذاتيتان حقيقيتان مختلفتا الإشارة — أي سرج؛ فيدخل مساران ويخرج مساران، وتمرّ البقية مرورًا. وهو غير مستقر دائمًا.
  • δ>0\delta > 0 و τ24δ\tau^2 \geq 4\delta: قيمتان ذاتيتان حقيقيتان بالإشارة نفسها (=signτ= \operatorname{sign}\tau) — أي عقدة، مستقرة إذا وفقط إذا كان τ<0\tau < 0؛ والمسارات مماسّة للاتجاه الذاتي البطيء.
  • δ>0\delta > 0 و τ2<4δ\tau^2 < 4\delta و τ0\tau \neq 0: قيمتان ذاتيتان عقديتان مرافقتان τ2±iω\frac\tau2 \pm \iu\omega — أي حلزون (بؤرة)، مستقر إذا وفقط إذا كان τ<0\tau < 0؛ والحلول دورانات eτt/2×\eu^{\tau t/2}\times بالدور 2πω\frac{2\pi}\omega.
  • τ=0\tau = 0 و δ>0\delta > 0: قيمتان ذاتيتان تخيّليتان محضتان — أي مركز: بمدارات مغلقة (قطوع ناقصة)، واستقرارٌ دون استقرار مقارب، وهو بالضبط الحدّ الفاصل الذي لا تستطيع المبرهنة 19.12 البتّ فيه من أجل الجمل غير الخطية (وهنا يجلس موضع توازن النواس السفلي، المسألة 19.1).

ويحمل القطع المكافئ الحدّي τ2=4δ\tau^2 = 4\delta العقدَ المنحلّة (أي كتل جوردان: بمسارات ذات اتجاه مماس واحد). وكل شيء يُقرأ من عددين — ولهذا يكون أول ردّ فعل قبل أي صورة طور مستوية هو حساب tr\operatorname{tr} و det\det؛ فمثلًا A=(011c)A = \bigl(\begin{smallmatrix}0 & 1\\ -1 & -c\end{smallmatrix}\bigr) (المذبذب المخمَّد): δ=1>0\delta = 1 > 0 و τ=c\tau = -c: أي حلزون مستقر من أجل 0<c<20 < c < 2، وعقدة مستقرة من أجل c2c \geq 2 — أي التخميد الناقص مقابل التخميد الزائد، في لمحة واحدة.

19.4 الانسيابات ومواضع التوازن والاستقرار

ننظر الآن في المعادلة الذاتية x=F(x)x' = F(x)، حيث F ⁣:ΩRdF \colon \Omega \to \R^d ليبشيتزية محليًّا على المفتوحة ΩRd\Omega \subseteq \R^d. ونكتب φt(x0)=x(t)\varphi_t(x_0) = x(t) للحل الأعظمي الذي يحقق x(0)=x0x(0) = x_0 (أي انسياب الحقل المتجهي)؛ وتعطي الذاتية خاصية الزمرة φt+s=φtφs\varphi_{t+s} = \varphi_t\circ\varphi_s حيثما عُرِّفت (إذ يحل الطرفان المسألة نفسها عند الزمن ss).

تعريف 19.10

موضع التوازن نقطةٌ xˉ\bar x تحقق F(xˉ)=0F(\bar x) = 0 (ومنه φt(xˉ)=xˉ\varphi_t(\bar x) = \bar x). ويكون مستقرًّا إذا وُجد من أجل كل ε>0\varepsilon > 0 عددٌ δ>0\delta > 0 بحيث يستلزم x0xˉ<δ\norm{x_0 - \bar x} < \delta وجود الحل من أجل كل t0t \geq 0 مع φt(x0)xˉ<ε\norm{\varphi_t(x_0) - \bar x} < \varepsilon؛ ويكون مستقرًّا مقاربًا إذا كان فوق ذلك φt(x0)xˉ\varphi_t(x_0) \to \bar x من أجل كل x0x_0 بجوار xˉ\bar x.

مبرهنة 19.11 (دوال ليابونوف)

ليكن xˉ\bar x موضع توازن ولتكن V ⁣:VRV \colon \mathcal V \to \R دالةً من الصنف C1\mathcal C^1 على جوار للنقطة xˉ\bar x تحقق:

V(xˉ)=0,V(x)>0 من أجلxxˉ,V˙(x)=V(x)F(x)0.V(\bar x) = 0,\qquad V(x) > 0 \text{ من أجل} x \neq \bar x, \qquad \dot V(x) = \nabla V(x)\cdot F(x) \leq 0 .

عندئذٍ يكون xˉ\bar x مستقرًّا. وإذا كان فوق ذلك V˙<0\dot V < 0 خارج xˉ\bar x، كان xˉ\bar x مستقرًّا مقاربًا.

برهان. على امتداد حل،  ⁣d ⁣dtV(x(t))=V˙(x(t))0\frac{\dd}{\dd t}V(x(t)) = \dot V(x(t)) \leq 0: أي إن VV متناقصة. وبإعطاء ε\varepsilon (صغير بما يكفي بحيث Bˉ(xˉ,ε)V\bar B(\bar x, \varepsilon) \subseteq \mathcal V)، لتكن m=min{V(x):xxˉ=ε}>0m = \min\{V(x) : \norm{x - \bar x} = \varepsilon\} > 0 (بالتراص والإيجابية) ونختار δ<ε\delta < \varepsilon بحيث V<mV < m على B(xˉ,δ)B(\bar x, \delta) (بالاتصال). فحلٌّ ينطلق من B(xˉ,δ)B(\bar x, \delta) يحقق V(x(t))<mV(x(t)) < m من أجل كل الأزمنة اللاحقة، ومنه لا يمكنه أن يبلغ الكرة xxˉ=ε\norm{x - \bar x} = \varepsilon (حيث VmV \geq m): فيبقى في الكرة — وعندئذٍ يوجد من أجل كل t0t \geq 0: إذ يبقى الحل في المتراصة Bˉ\bar B، ومنه تفرض المبرهنة 19.5(2) (الهروب من المتراصات) أن T+=+T_+ = +\infty. وهذا هو الاستقرار.

وأما الحالة المقاربة: فليكن x(t)x(t) منطلقًا من B(xˉ,δ)B(\bar x, \delta)؛ تتناقص V(x(t))V(x(t)) إلى c0c \geq 0 ما. فإذا كان c>0c > 0: بقي المسار في K={xBˉ(xˉ,ε):V(x)c}K = \{x \in \bar B(\bar x, \varepsilon): V(x) \geq c\}، وهي مجموعة متراصة تستبعد جوارًا للنقطة xˉ\bar x (إذ VV متصلة مع V(xˉ)=0<cV(\bar x) = 0 < c). وعلى KK تكون الدالة V˙\dot V متصلة وسالبة تمامًا، ومنه μ=maxKV˙<0\mu = \max_K\dot V < 0 (بالتراص)؛ عندئذٍ V(x(t))V(x(0))+μtV(x(t)) \leq V(x(0)) + \mu t \to -\infty: وهو محال لأن V0V \geq 0. ومنه c=0c = 0، و x(t)xˉx(t) \to \bar x (إذ للنقاط التي تبعد ρ\geq \rho عن xˉ\bar x داخل الكرة Vmρ>0V \geq m_\rho > 0).

مبرهنة 19.12 (الاستقرار بالخطّية)

لتكن FF من الصنف C1\mathcal C^1 مع F(xˉ)=0F(\bar x) = 0 و A=DF(xˉ)A = DF(\bar x). إذا كان لكل قيمة ذاتية للمصفوفة AA الشرطُ Reλ<0\operatorname{Re}\lambda < 0، فإن xˉ\bar x مستقر مقاربًا.

برهان. ننقل xˉ\bar x إلى 00 ونكتب F(x)=Ax+g(x)F(x) = Ax + g(x) مع g(x)=o(x)g(x) = o(\norm x) (بالقابلية للاشتقاق من الصنف C1\mathcal C^1). ونختار α>0\alpha > 0 يحقق etACeαt\vertiii{\eu^{tA}} \leq C\eu^{-\alpha t} (حيث t0t \geq 0؛ المبرهنة 19.8) و r>0r > 0 يحقق g(x)α2Cx\norm{g(x)} \leq \frac{\alpha}{2C}\norm x من أجل xr\norm x \leq r. وبدوهاميل مع b(s)=g(x(s))b(s) = g(x(s)):

x(t)=etAx0+0te(ts)Ag(x(s)) ⁣ds,x(t) = \eu^{tA}x_0 + \int_0^t\eu^{(t-s)A}g(x(s))\,\dd s,

وهي صحيحة ما دام x(s)r\norm{x(s)} \leq r. عندئذٍ يحقق u(t)=eαtx(t)u(t) = \eu^{\alpha t}\norm{x(t)}

u(t)Cx0+0tCα2Cu(s) ⁣ds,u(t) \leq C\norm{x_0} + \int_0^{t}C\,\frac{\alpha}{2C}\,u(s)\,\dd s ,

ومنه تعطي غرونوال أن u(t)Cx0eαt/2u(t) \leq C\norm{x_0}\eu^{\alpha t/2}، أي x(t)Cx0eαt/2\norm{x(t)} \leq C\norm{x_0}\eu^{-\alpha t/2}. وإذا كان x0<r/C\norm{x_0} < r/C، أبقى الحدّ القبلي x(t)<r\norm{x(t)} < r من أجل كل tt (بحجة اتصال أو ترقية: إذ مجموعة الأزمنة التي يكون فيها xr\norm x \leq r مفتوحة ومغلقة في [0,T+)\intco0{T_+} بحكم التقدير الأكيد)، ويكون الحل شاملًا، ويتقارب إلى 00 أُسّيًّا: أي استقرار مقارب.

طريقة 19.13

حين تواجه معادلة تفاضلية: (1) الوجود والوحدانية — تحقق من الليبشيتزية المحلية (وهي عادةً C1\mathcal C^1)؛ (2) الشمول — بالنمو الخطي، أو الحدّية، أو متراصة صامدة عبر دالة ليابونوف أو تكامل أول؛ وإن لم يكن، فاشتبه بالانفجار واختبر على الكاريكاتير السلّمي x=x2x' = x^2؛ (3) الجمل الخطية — المُحلِّل، ودوهاميل، وبنيةُ القيم الذاتية للمصفوفة AA من أجل المعاملات الثابتة؛ (4) الأسئلة النوعية — مواضع التوازن، والخطّية، وصيدُ دالة ليابونوف (أي الطاقة حين تكون الجملة ميكانيكية) أو تكاملٍ أول تحبس مجموعاتُ مستوياته المسارات. وتمشي مسألة نهاية الأسبوع بالطريقة كلها عبر النواس.

صورة الطور للنواس x'' = - x في المستوي (x, x'): أي المنحنيات المستوية للطاقة E = x'22 - x. والمنحنيات المغلقة (بالأزرق): تذبذبات، E < 1؛ والمنحنيات الجارية (بالبرتقالي): دورات كاملة، E > 1؛ وبينهما الفاصلة (بالأحمر)، E = 1، الواصلة بين موضعَي التوازن غير المستقرين (±π, 0). وتبرهن مسألة نهاية الأسبوع على كل ما توحي به هذه الصورة.
صورة الطور للنواس x=sinxx'' = -\sin x في المستوي (x,x)(x, x'): أي المنحنيات المستوية للطاقة E=x22cosxE = \frac{x'^2}2 - \cos x. والمنحنيات المغلقة (بالأزرق): تذبذبات، E<1E < 1؛ والمنحنيات الجارية (بالبرتقالي): دورات كاملة، E>1E > 1؛ وبينهما الفاصلة (بالأحمر)، E=1E = 1، الواصلة بين موضعَي التوازن غير المستقرين (±π,0)(\pm\pi, 0). وتبرهن مسألة نهاية الأسبوع على كل ما توحي به هذه الصورة.

19.5 تمارين

تمرين 19.1

حُلّ صراحةً وعيّن الفترة الأعظمية: (a) x=x2x' = x^2 مع x(0)=1x(0) = 1؛ (b) x=1+x2x' = 1 + x^2 مع x(0)=0x(0) = 0؛ (c) x=x(1x)x' = x(1-x) مع x(0)=12x(0) = \frac12. ووفّق بين كل جواب والمبرهنة 19.5(2) والنتيجة 19.6.

حل

حل التمرين 19.1.

(a) بفصل المتغيّرات: x(t)=11tx(t) = \frac1{1 - t} على (,1)\intoo{-\infty}1: أي انفجار عند T+=1T_+ = 1 مع x(t)+x(t) \to +\infty — وهو بالضبط المبرهنة 19.5(2). (b) x(t)=tantx(t) = \tan t على (π/2,π/2)\intoo{-\pi/2}{\pi/2}: أي انفجار عند الطرفين. (c) x(t)=11+etx(t) = \frac{1}{1 + \eu^{-t}}، وهو شامل: إذ يبقى الحل في (0,1)\intoo01، وهي مجموعة محدودة، ومنه لا يمكن للمنحني البياني أن يهرب من كل متراصة من R×R\R\times\R في زمن منتهٍ — T±=±T_\pm = \pm\infty. ولاحظ أن (a) و(b) لا تناقضان النتيجة 19.6: إذ للمقدارين x2x^2 و 1+x21 + x^2 نموٌّ فوق خطي.

تمرين 19.2

ليكن x,yx, y حلّين للمعادلة x=f(t,x)x' = f(t,x) حيث ff ليبشيتزية شاملًا بالثابت LL في xx على R×Rd\R\times\R^d. (a) برهن على x(t)y(t)x(0)y(0)eLt\norm{x(t) - y(t)} \leq \norm{x(0) - y(0)}\eu^{L\abs t}، وبيّن بمثال (خطي!) أن العامل eLt\eu^{L\abs t} مبلوغ. (b) استنتج أن تطبيق الانسياب x0x(t;x0)x_0 \mapsto x(t; x_0) متصل، بل ليبشيتزي على المجموعات المحدودة.

حل

حل التمرين 19.2.

(a) هذا هو تقدير النتيجة 19.4 مع t0=0t_0 = 0. وأما الضبط: فمن أجل x=Lxx' = Lx (وهي ليبشيتزية شاملًا بالثابت LL)، يختلف حلّان بالمقدار (x0y0)eLt(x_0 - y_0)\eu^{Lt} بالضبط. (b) ويُقرأ التقدير: أن تطبيق الانسياب عند الزمن tt ليبشيتزي بالثابت eLt\eu^{L\abs t} في الشرط الابتدائي — أي الاتصال، بانتظام من أجل tt في المتراصات؛ وعلى المجموعات المحدودة من أجل ff غير الليبشيتزية شاملًا، نجري الشيء نفسه على أنبوب متراص حول المسارات بالثابت المحلي.

تمرين 19.3 ★★

برهن على أن لكلٍّ مما يلي حلولًا أعظمية شاملة على R\R، مستشهدًا بالمبرهنة الصحيحة: (a) x=sin(tx)x' = \sin(tx)؛ (b) x=tx1+x2x' = \frac{t\,x}{1 + x^2}؛ (c) x+q(t)x=0x'' + q(t)x = 0 مع qq متصلة (بالتحويل إلى جملة من الرتبة الأولى)؛ (d) x=A(t)xx' = A(t)x مع AA متصلة ومحدودة — وأعطِ حدّ غرونوال للمقدار x(t)\norm{x(t)}.

حل

حل التمرين 19.3.

(a) sin(tx)1\abs{\sin(tx)} \leq 1: أي محدود، وهو نمو خطي مع α=0\alpha = 0 و β=1\beta = 1: النتيجة 19.6 على U=R×RU = \R\times\R. (b) tx/(1+x2)t12\abs{tx/(1 + x^2)} \leq \abs t\cdot\frac12: وهو تحت خطي مرة أخرى (بل محدود على المتراصات الزمنية): فشامل. (c) X=(x,x)X = (x, x'): X=(01q(t)0)XX' = \bigl(\begin{smallmatrix}0 & 1\\ -q(t) & 0\end{smallmatrix}\bigr)X: أي خطية بمعاملات متصلة: فشاملة (في إطار المبرهنة 19.7). (d) شامل؛ وغرونوال كما في النتيجة 19.6: x(t)x(t0)eMtt0\norm{x(t)} \leq \norm{x(t_0)}\,\eu^{M\abs{t - t_0}} مع M=supAM = \sup\vertiii{A}.

تمرين 19.4 ★★

(a) احسب etA\eu^{tA} من أجل A=(0110)A = \bigl(\begin{smallmatrix}0 & -1\\ 1 & 0\end{smallmatrix}\bigr)، و (λ10λ)\bigl(\begin{smallmatrix}\lambda & 1\\ 0 & \lambda\end{smallmatrix}\bigr)، و (0110)\bigl(\begin{smallmatrix}0 & 1\\ 1 & 0\end{smallmatrix}\bigr). (b) حُلّ المذبذب المُقحَم x+x=cos(ωt)x'' + x = \cos(\omega t) بدوهاميل (بصيغة الجملة)، من أجل ω1\omega \neq 1 ومن أجل ω=1\omega = 1: فيظهر الرنين بوصفه الحدّ الدنيوي tsintt\sin t.

حل

حل التمرين 19.4.

(a) A2=IA^2 = -I من أجل الأولى: etA=costI+sintA=(costsintsintcost)\eu^{tA} = \cos t\,I + \sin t\,A = \bigl(\begin{smallmatrix}\cos t & -\sin t\\ \sin t & \cos t\end{smallmatrix}\bigr). وكتلة جوردان: إذ تتبادل λI\lambda I و NN: etA=eλt(1t01)\eu^{tA} = \eu^{\lambda t}\bigl(\begin{smallmatrix}1 & t\\ 0 & 1\end{smallmatrix}\bigr). والثالثة: A2=IA^2 = I: etA=coshtI+sinhtA\eu^{tA} = \cosh t\,I + \sinh t\,A.

(b) الجملة X=(0110)X+(0cosωt)X' = \bigl(\begin{smallmatrix}0&1\\-1&0 \end{smallmatrix}\bigr)X + \bigl(\begin{smallmatrix}0\\ \cos\omega t\end{smallmatrix}\bigr)؛ ويعطي دوهاميل مع مُحلِّل الدوران الحلولَ الخاصة: فمن أجل ω1\omega \neq 1، xp(t)=cos(ωt)1ω2x_p(t) = \frac{\cos(\omega t)}{1 - \omega^2} (تحقق مباشرةً)؛ ومن أجل ω=1\omega = 1 ينتج التكامل 0tsin(ts)coss ⁣ds=t2sint\int_0^t\sin(t - s)\cos s\,\dd s = \frac t2\sin t النموَّ الدنيوي xp=t2sintx_p = \frac t2\sin t: أي الرنين — إذ يضخّ الإقحام الطاقة عند التردد الطبيعي فتنمو السعة خطيًّا.

تمرين 19.5 ★★

من أجل المعادلة السلّمية x+p(t)x+q(t)x=0x'' + p(t)x' + q(t)x = 0: (a) برهن على أن الفرونسكياني w=x1x2x1x2w = x_1x_2' - x_1'x_2 لحلّين يحقق w=pww' = -p\,w (أبيل)، واستنتج أن حلّين يكون w0w \neq 0 من أجلهما في موضع ما يشكّلان أساسًا. (b) وبإعطاء حل واحد غير منعدم x1x_1، جد الحل العام بتخفيض الرتبة: ضع x2=x1epx12x_2 = x_1\int \frac{\eu^{-\int p}}{x_1^2} وتحقق. وطبّق ذلك على t2x2x=0t^2x'' - 2x = 0 على (0,+)\intoo0{+\infty} مع x1(t)=t2x_1(t) = t^2.

حل

حل التمرين 19.5.

(a) w=x1x2x1x2=x1(px2qx2)(px1qx1)x2=pww' = x_1x_2'' - x_1''x_2 = x_1(-px_2' - qx_2) - (-px_1' - qx_1)x_2 = -p\,w: w(t)=w(t0)et0tpw(t) = w(t_0)\eu^{-\int_{t_0}^tp}، وهو لا ينعدم أبدًا ولا ينعدم انعدامًا تامًّا. وإذا كان w0w \neq 0، كانت المتجهتان (xi,xi)(t0)(x_i, x_i')(t_0) مستقلتين في R2\R^2، وبما أن فضاء الحلول ذو البُعد 22 (المبرهنة 19.7 من أجل الجملة)، تكون (x1,x2)(x_1, x_2) أساسًا.

(b) مع u=epx12u = \int\frac{\eu^{-\int p}}{x_1^2}: x2=x1ux_2 = x_1u، و x2=x1u+epx1x_2' = x_1'u + \frac{\eu^{-\int p}}{x_1}، و

x2+px2+qx2=u(x1+px1+qx1)+(pepx1+pepx1)=0x_2'' + px_2' + qx_2 = u\,(x_1'' + px_1' + qx_1) + \Bigl(-\,p\frac{\eu^{-\int p}}{x_1} + p\frac{\eu^{-\int p}}{x_1}\Bigr) = 0

(إذ تتلاشى الحدود المتقاطعة بالضبط؛ فانشر بعناية). ومن أجل t2x2x=0t^2x'' - 2x = 0، أي x2t2x=0x'' - \frac2{t^2}x = 0 (حيث p=0p = 0) مع x1=t2x_1 = t^2: u=t4=13t3u = \int t^{-4} = -\frac1{3t^3}، ومنه x2=13tx_2 = -\frac1{3t}: فالحل العام at2+btat^2 + \frac bt.

تمرين 19.6 ★★

(اللوجستية) من أجل x=x(1x)x' = x(1 - x): عيّن جميع مواضع التوازن واستقرارها (بواسطة المبرهنة 19.12 ومباشرةً)؛ وبرهن على أن كل حل يحقق x(0)(0,1)x(0) \in \intoo01 متزايدٌ وشامل، بنهايتين 00 و 11 عند \mp\infty؛ وحُلّ صراحةً للتأكيد. وبرهن عمومًا على أن الحلول الذاتية السلّمية رتيبة، واخلص: أنه لا حلول دورية غير ثابتة في البُعد 11.

حل

حل التمرين 19.6.

مواضع التوازن 0,10, 1؛ و f(x)=12xf'(x) = 1 - 2x: فإن f(0)=1>0f'(0) = 1 > 0 (غير مستقر — إذ تبتعد الحلول المجاورة، كما تبيّن الصيغة الصريحة)، و f(1)=1<0f'(1) = -1 < 0: أي مستقر مقاربًا (المبرهنة 19.12 في البُعد 11). ومن أجل x(0)(0,1)x(0) \in \intoo01: يكون f>0f > 0 هناك، ومنه يتزايد الحل ما دام في (0,1)\intoo01؛ ولا يمكنه أبدًا أن يبلغ 00 أو 11 (بالوحدانية: إذ هما مساران)، ومنه يبقى، ويكون محدودًا — ومن ثَمّ شاملًا — ويتزايد إلى نهاية L(x(0),1]L \in \intoc{x(0)}1. فإذا كان f(L)0f(L) \neq 0، كان xc>0x' \geq c > 0 بجوار النهاية، فيدفع xx وراء LL: ومنه f(L)=0f(L) = 0، أي L=1L = 1؛ وبالتناظر x0x \to 0 عند -\infty. ويؤكد الحل الصريح x(t)=11+Cetx(t) = \frac1{1 + C\eu^{-t}} كلَّ ذلك. وعمومًا: لو كان لحلٍّ ذاتي سلّمي x(t0)=0x'(t_0) = 0، لكان x(t0)x(t_0) موضعَ توازن ولجعلت الوحدانية xx ثابتة؛ وإلا احتفظ f(x(t))f(x(t)) بإشارة ثابتة (إذ لا ينعدم أبدًا، و tf(x(t))t \mapsto f(x(t)) متصلة): فتكون xx رتيبة تمامًا — ومنه يستحيل وجود حل دوري غير ثابت.

تمرين 19.7 ★★

(التكاملات الأولى) لتكن H ⁣:ΩRH \colon \Omega \to \R من الصنف C1\mathcal C^1 ولننظر في الجملة الهاملتونية المستوية x=yHx' = \partial_yH، y=xHy' = -\partial_xH. (a) برهن على أن HH ثابتة على امتداد الحلول. (b) ومن أجل H=y22+x44H = \frac{y^2}2 + \frac{x^4}4: برهن على أن جميع الحلول شاملة ومحدودة، وعلى أن المبدأ مستقر (بليابونوف: HH) رغم أن الخطّية ((0100)\bigl(\begin{smallmatrix}0&1\\0&0\end{smallmatrix}\bigr)) ليست مستقرة مقاربًا: فالخطّية قد تكون غير حاسمة.

حل

حل التمرين 19.7.

(a)  ⁣d ⁣dtH(x,y)=Hxx+Hyy=HxHy+Hy(Hx)=0\frac{\dd}{\dd t}H(x,y) = H_xx' + H_yy' = H_xH_y + H_y(-H_x) = 0. (b) مجموعات المستويات للمقدار H=y22+x44H = \frac{y^2}2 + \frac{x^4}4 متراصة (إذ HH قسرية)، ومنه تُحبس الحلول في متراصات: فهي شاملة ومحدودة (المبرهنة 19.5). وأما استقرار (0,0)(0,0): فإن V=HV = H معيَّنة موجبةً (إذ H=0H = 0 عند المبدأ وحده) مع V˙=0\dot V = 0: المبرهنة 19.11. والخطّية x=yx' = y و y=0y' = 0 لها المصفوفة المعدومة القوة غير القابلة للقُطرنة بقيمة ذاتية 00: فتصمت المبرهنة 19.12 (إذ تخفق فرضيتها Reλ<0\operatorname{Re}\lambda < 0)، وبالفعل تكون الجملة المخطَّطة غير مستقرة (إذ ينجرف y00y_0 \ne 0) بينما تكون غير الخطية مستقرة: فالخطّية عند موضع توازن غير زائدي لا تبرهن على شيء.

تمرين 19.8 ★★★

(النواس المخمَّد) x+cx+sinx=0x'' + cx' + \sin x = 0 مع c>0c > 0؛ وبصيغة الجملة: x=yx' = y و y=sinxcyy' = -\sin x - cy. (a) برهن على أن V(x,y)=y22+1cosxV(x, y) = \frac{y^2}2 + 1 - \cos x يحقق V˙=cy20\dot V = -cy^2 \leq 0: أي إن المبدأ مستقر. (b) وينعدم V˙\dot V على المحور y=0y = 0 كله: فيخفق محك ليابونوف الأكيد. برهن على الاستقرار المقارب رغم ذلك، بالخطّية (المبرهنة 19.12): احسب القيم الذاتية للمصفوفة المخطَّطة عند (0,0)(0,0) وتحقق من Re<0\operatorname{Re} < 0 من أجل كل c>0c > 0. (c) وماذا يحدث عند موضع التوازن (π,0)(\pi, 0)؟ احسب الخطّية واخلص (إذ قيمة ذاتية واحدة موجبة: أي عدم استقرار — ولك أن تستعمل عبارة عدم الاستقرار بصورة غير رسمية أو تنتج حلًّا هاربًا صريحًا للجملة الخطية).

حل

حل التمرين 19.8.

(a) V˙=yy+sinxx=y(sinxcy)+ysinx=cy20\dot V = y\,y' + \sin x\cdot x' = y(-\sin x - cy) + y\sin x = -cy^2 \leq 0، و V=y22+(1cosx)V = \frac{y^2}2 + (1 - \cos x) معيَّنة موجبةً على {x<2π}\{\abs x < 2\pi\} حول المبدأ: أي مستقر (المبرهنة 19.11). (b) المصفوفة المخطَّطة عند (0,0)(0,0) هي (011c)\bigl(\begin{smallmatrix}0 & 1\\ -1 & -c\end{smallmatrix}\bigr)، وكثير حدودها المميِّز λ2+cλ+1\lambda^2 + c\lambda + 1: بجذرين c±c242\frac{-c \pm \sqrt{c^2 - 4}}2 — وكلاهما حقيقي سالب إذا كان c2c \geq 2، وعقديان بجزء حقيقي c2<0-\frac c2 < 0 إذا كان 0<c<20 < c < 2. وفي جميع الحالات Reλ<0\operatorname{Re}\lambda < 0: فتعطي المبرهنة 19.12 الاستقرار المقارب (رغم انحلال V˙\dot V). (c) وعند (π,0)(\pi, 0): sin(π+u)=sinu\sin(\pi + u) = -\sin u، والخطّية (011c)\bigl(\begin{smallmatrix}0&1\\1&-c\end{smallmatrix}\bigr)، والمميِّز λ2+cλ1\lambda^2 + c\lambda - 1: بجذرين مختلفَي الإشارة (λ+λ=1\lambda_+\lambda_- = -1). وعلى امتداد المتجهة الذاتية غير المستقرة، للجملة الخطية الحلُّ الهارب صراحةً eλ+tv+\eu^{\lambda_+t}v_+ مع λ+>0\lambda_+ > 0: فالنواس المقلوب غير مستقر من أجل كل تخميد.

تمرين 19.9 ★★

(حدود الوحدانية) من أجل α(0,1)\alpha \in \intoo01، برهن على أن للمسألة x=xαx' = \abs x^{\alpha} مع x(0)=0x(0) = 0 عددًا لانهائيًّا من الحلول (بتكييف المسألة 7.1، الجزء الثالث). وبرهن بالعكس على أنه من أجل α=1\alpha = 1 (أي x=xx' = \abs x) يكون الحل المارّ بالنقطة 00 وحيدًا، وعيّن بالضبط أي فرضية من المبرهنة 19.2 تميّز بين الحالتين.

حل

حل التمرين 19.9.

من أجل α(0,1)\alpha \in \intoo01: إضافةً إلى x0x \equiv 0، تكون كل

xc(t)={0tc,((1α)(tc))1/(1α)tc,x_c(t) = \begin{cases}0 & t \leq c,\\ \bigl((1-\alpha)(t - c)\bigr)^{1/(1-\alpha)} & t \geq c, \end{cases}

من الصنف C1\mathcal C^1 وتحل المعادلة (إذ يجعل الأُسّ 11α>1\frac1{1-\alpha} > 1 المشتقةَ منعدمة عند cc): أي متصلٌ من الحلول المارّة بالنقطة (0,0)(0,0). ومن أجل α=1\alpha = 1: تكون xxx \mapsto \abs x ليبشيتزية شاملًا بالثابت 11 (abab\abs{\abs a - \abs b} \leq \abs{a - b})، ومنه تنطبق المبرهنة 19.2 ويكون الحل الوحيد المارّ بالنقطة 00 هو x0x \equiv 0. والحدّ الفاصل هو بالضبط شرط الليبشيتزية المحلية عند 00: إذ لنسب فروق xα\abs x^\alpha هناك عدمُ حدّية من أجل α<1\alpha < 1.

تمرين 19.10 ★★★

(الحدود القبلية تحبس الحلول) لتكن F ⁣:RdRdF \colon \R^d \to \R^d ليبشيتزية محليًّا مع F(x),x0\langle F(x), x\rangle \leq 0 كلما كان xR\norm x \geq R. (a) برهن على أن الكرة المغلقة Bˉ(0,R)\bar B(0, R) صامدة إيجابيًّا: أي إن الحلول التي تنطلق من الداخل تبقى فيه من أجل t0t \geq 0. (إذا كان x(t2)>R\norm{x(t_2)} > R، فانظر في آخر زمن t1<t2t_1 < t_2 يحقق x(t1)=R\norm{x(t_1)} = R وادرس  ⁣d ⁣dtx(t)2\frac{\dd}{\dd t}\norm{x(t)}^2 على [t1,t2]\intcc{t_1}{t_2}.) (b) استنتج الوجود الأمامي الشامل من أجل المعطيات في الكرة. ثم عالج الجملة التدرّجية x=G(x)x' = -\nabla G(x)، GC2G \in \mathcal C^2 مع G(x)+G(x) \to +\infty حين x\norm x \to \infty: برهن على أن GG متناقصة على امتداد الحلول، وعلى أن كل حل يبقى في مجموعة المستوى الأدنى (المحدودة) {GG(x0)}\{G \leq G(x_0)\}، واخلص إلى الوجود الأمامي الشامل.

حل

حل التمرين 19.10.

(a) لنفترض أن x(t2)>R\norm{x(t_2)} > R من أجل t2>0t_2 > 0 ما مع x(0)R\norm{x(0)} \leq R، ولتكن t1=sup{tt2:x(t)R}t_1 = \sup\{t \leq t_2 : \norm{x(t)} \leq R\}: عندئذٍ x(t1)=R\norm{x(t_1)} = R و x(t)>R\norm{x(t)} > R على (t1,t2]\intoc{t_1}{t_2}. وعلى تلك الفترة يحقق g(t)=x(t)2g(t) = \norm{x(t)}^2 أن g(t)=2x(t),F(x(t))0g'(t) = 2\langle x(t), F(x(t))\rangle \leq 0 (إذ تنطبق الفرضية: x(t)R\norm{x(t)} \geq R)، ومنه g(t2)g(t1)=R2g(t_2) \leq g(t_1) = R^2: وهو تناقض. فالكرة صامدة إيجابيًّا. (b) ولا يمكن لحلٍّ محبوس في الكرة المتراصة أن يحقق T+<T_+ < \infty (المبرهنة 19.5(2)): أي شامل أماميًّا. وأما الجملة التدرّجية:  ⁣d ⁣dtG(x(t))=G,G=G(x(t))20\frac{\dd}{\dd t}G(x(t)) = \langle\nabla G, -\nabla G\rangle = -\norm{\nabla G(x(t))}^2 \leq 0: فتتناقص GG، ومنه يبقى الحل في {GG(x0)}\{G \leq G(x_0)\}، وهي محدودة (بالقسرية: إذ خارج كرة كبيرة يكون G>G(x0)G > G(x_0)) ومغلقة: أي متراصة. فالهروب مستحيل: أي إن كل حل لجملة تدرّجية قسرية شامل أماميًّا، منزلقًا إلى أسفل إلى الأبد.

تمرين 19.11 ★★

(الانفجار بالمقارنة) لننظر في x=x2+t2x' = x^2 + t^2 مع x(0)=1x(0) = 1. (a) برهن على أن الحل الأعظمي موجود على [0,T+)\intco0{T_+} ما مع T+<T_+ < \infty: بالمقارنة مع y=y2y' = y^2 و y(0)=1y(0) = 1 (وبرهن على المبرهنة المساعدة للمقارنة التي تحتاجها: أنه إذا كان xF(x)x' \geq F(x) وكان y=F(y)y' = F(y) مع x(0)y(0)x(0) \geq y(0)، فإن xyx \geq y حيثما يعيش كلاهما)، واستنتج T+1T_+ \leq 1. (b) وحُدّ T+T_+ من الأسفل: فعلى [0,1]\intcc01، xx2+1x' \leq x^2 + 1؛ وقارن مع الحل الأعلى z=z2+1z' = z^2 + 1 مع z(0)=1z(0) = 1، المحلول بواسطة z(t)=tan(t+π4)z(t) = \tan\bigl(t + \frac\pi4\bigr)، واخلص إلى T+π4T_+ \geq \frac\pi4. (c) اجمع π4T+1\frac\pi4 \leq T_+ \leq 1 وصُغ العبرة: أن النمو الفوق خطي للطرف الأيمن هو ما يقتل الوجود الشامل (التمرين 19.3 بوصفه النقيض)، والحدّ الفاصل هو تقارب  ⁣dsF(s)\int^{\infty}\frac{\dd s}{F(s)}.

حل

حل التمرين 19.11.

(a) المبرهنة المساعدة للمقارنة: لتكن w=xyw = x - y على الفترة المشتركة؛ فإن w(0)0w(0) \geq 0 و w=xyF(x)F(y)=c(t)ww' = x' - y' \geq F(x) - F(y) = c(t)w مع c(t)=F(x)F(y)xyc(t) = \frac{F(x) - F(y)}{x - y} محدودة على الفترات الزمنية المتراصة (لأن FF ليبشيتزية محليًّا)؛ عندئذٍ (wec)0(w\eu^{-\int c})' \geq 0، ومنه w0w \geq 0 في كل مكان. ومع F(x)=x2F(x) = x^2: ينفجر y(t)=11ty(t) = \frac1{1 - t} عند 11، و xyx \geq y ما دام xx حيًّا؛ فإذا كان T+>1T_+ > 1، لكان xx منتهيًا عند t=1t = 1 بينما يهيمن على yy \to \infty: وهو محال. ومنه T+1T_+ \leq 1.

(b) والمقارنة المعكوسة (بالمبرهنة نفسها، بتبادل الدورين): على [0,1][0,T+)\intcc01\cap\intco0{T_+}، يعطي t21t^2 \leq 1 أن xx2+1x' \leq x^2 + 1، بينما يحقق z(t)=tan(t+π4)z(t) = \tan(t + \frac\pi4) أن z=z2+1z' = z^2 + 1 و z(0)=1=x(0)z(0) = 1 = x(0): ومنه xzx \leq z ما دام كلاهما معرَّفًا. وبما أن zz منتهٍ على [0,π4)\intco0{\frac\pi4}، لا يمكن للحل xx أن ينفجر قبل π4\frac\pi4: أي T+π4T_+ \geq \frac\pi4.

(c) ومعًا: π4T+1\frac\pi4 \leq T_+ \leq 1 (وعدديًّا T+0.96T_+ \approx 0.96). والعبرة: أنه من أجل x=F(t,x)x' = F(t, x) مع FF فوق خطية في xx، تنفجر الحلول في زمن منتهٍ كلما كان  ⁣dsF(s)<\int^\infty\frac{\dd s}{F(s)} < \infty (إذ يبلغ حل المقارنة اللانهاية في ذلك الزمن المنتهي)؛ وأما النمو الخطي، حيث يتباعد التكامل، فيفرض الوجودَ الشامل (التمرين 19.3). وهو تكامل أوسغود نفسه كما في المسألة 7.1، لكنه يحكم الآن الهروبَ إلى اللانهاية بدل الهروب من الصفر.

تمرين 19.12 ★★★

(مبرهنة شتورم في المقارنة) ليكن q1q2q_1 \leq q_2 متصلتين على فترة II، ولتكن u0u \neq 0 حلًّا للمعادلة u+q1u=0u'' + q_1u = 0، ولتكن v0v \neq 0 حلًّا للمعادلة v+q2v=0v'' + q_2v = 0. (a) أثبت متطابقة الفرونسكياني: أنه مع W=uvuvW = uv' - u'v، W=(q1q2)uvW' = (q_1 - q_2)\,uv. (b) (شتورم) برهن على أنه بين جذرين متتاليين a<ba < b للدالة uu، إما أن تنعدم vv في موضع ما من (a,b)\intoo ab، وإما أن يكون q1=q2q_1 = q_2 و vuv \propto u هناك (افترض u>0u > 0 على (a,b)\intoo ab و v>0v > 0 أيضًا؛ وكامل (a) من aa إلى bb وافحص إشارتَي الحدّين الحدّيين W(a),W(b)W(a), W(b)). (c) استنتج: أن حلول u+qu=0u'' + qu = 0 مع qm>0q \geq m > 0 تنعدم مرة واحدة على الأقل في كل فترة طولها π/m\pi/\sqrt m (بالمقارنة مع v+mv=0v'' + mv = 0)؛ وأن الحلول مع q0q \leq 0 تنعدم مرة واحدة على الأكثر على R\R. واختبر كليهما على q±1q \equiv \pm1.

حل

حل التمرين 19.12.

(a) W=uvuv=u(q2v)(q1u)v=(q1q2)uvW' = uv'' - u''v = u(-q_2v) - (-q_1u)v = (q_1 - q_2)\,uv.

(b) ليكن a<ba < b جذرين متتاليين للدالة uu؛ ونعيِّر بحيث u>0u > 0 على (a,b)\intoo ab (ومنه u(a)>0u'(a) > 0 و u(b)<0u'(b) < 0 — وهما غير معدومين بالوحدانية، إذ كان u(a)=u(a)=0u(a) = u'(a) = 0 سيفرض u0u \equiv 0). ولنفترض أن vv لا تقبل أي جذر في (a,b)\intoo ab؛ ونعيِّر v>0v > 0 هناك (ومنه v(a),v(b)0v(a), v(b) \geq 0 بالاتصال). ونكامل (a):

W(b)W(a)=ab(q1q2)uv    0.W(b) - W(a) = \int_a^b(q_1 - q_2)\,uv \;\leq\; 0 .

لكن W(a)=u(a)v(a)u(a)v(a)=u(a)v(a)0W(a) = u(a)v'(a) - u'(a)v(a) = -u'(a)v(a) \leq 0 و W(b)=u(b)v(b)0W(b) = -u'(b)v(b) \geq 0: ومنه W(b)W(a)0W(b) - W(a) \geq 0. والتساوي في كل المراحل: فإن (q1q2)uv=0\int(q_1 - q_2)uv = 0 مع uv>0uv > 0 على الفترة المفتوحة يفرض q1=q2q_1 = q_2 هناك؛ ويفرض W(a)=W(b)=0W(a) = W(b) = 0 أن v(a)=v(b)=0v(a) = v(b) = 0؛ عندئذٍ W0W \equiv 0 على [a,b]\intcc ab (لأن مشتقته منعدمة)، أي (v/u)=W/u2=0(v/u)' = -W/u^2 = 0 على (a,b)\intoo ab: أي vuv \propto u.

(c) نأخذ q1=mq_1 = m و u=sin(m(tt0))u = \sin(\sqrt m(t - t_0))، وجذوره المتتالية تتباعد بمقدار π/m\pi/\sqrt m، و q2=qmq_2 = q \geq m: فحسب (b) تنعدم كل حل vv للمعادلة v+qv=0v'' + qv = 0 في كل فترة مفتوحة طولها π/m\pi/\sqrt m (وفي البديل المنحلّ qmq \equiv m، تنعدم vuv \propto u أيضًا). وإذا كان q0q \leq 0 بدلًا من ذلك: طبّقنا (b) مع q1=qq_1 = q و u=vu = v، ومع q2=0q_2 = 0 والحل العديم الجذور 1\mathbf 1 للمعادلة v=0v'' = 0. فلو قبلت vv جذرين متتاليين، لفرضت (b) إما جذرًا للمقدار 1\mathbf 1 بينهما وإما الحالة المنحلّة 1v\mathbf 1 \propto v — وكلاهما محال: ومنه تنعدم vv مرة واحدة على الأكثر. والفحوص: من أجل q=1q = 1، ينعدم sint\sin t كل π=π/1\pi = \pi/\sqrt1، كما تنبّأنا؛ ومن أجل q=1q = -1، ينعدم sinht\sinh t مرة واحدة بالضبط ولا ينعدم et\eu^t أبدًا — أي جذر واحد على الأكثر، كما تنبّأنا.

19.6 مسألة: النواس، محلولًا تمامًا

مسألة 19.1

مسألة نهاية الأسبوع — التذبذبات والدورانات والفاصلة والدور

معادلة النواس x=sinxx'' = -\sin x — وبصيغة الجملة: x=yx' = y، y=sinxy' = -\sin x على R2\R^2 — هي ذبابة الفاكهة في الديناميكا: سهلة الكتابة، مستحيلة الحل بالصيغ الابتدائية، ومع ذلك مفهومة تمامًا بالطريقة النوعية. ولتكن E(x,y)=y22cosxE(x, y) = \frac{y^2}2 - \cos x (أي الطاقة).

الجزء الأول — البنية الشاملة.

  1. برهن على أن جميع الحلول الأعظمية شاملة (معرَّفة على R\R): باستعمال E˙=0\dot E = 0 والمبرهنة 19.5. ومواضع التوازن: (kπ,0)(k\pi, 0)؛ صنّف خطّياتها (من نمط المركز من أجل kk زوجي، ومن نمط السرج من أجل kk فردي).
  2. برهن على أن المسارات محتواة في مجموعات المستويات {E=E0}\{E = E_0\}، وارسم أو صف هذه المسارات حسب قيمة E0[1,+)E_0 \in \intco{-1}{+\infty}: E0=1E_0 = -1 (مواضع التوازن)، و 1<E0<1-1 < E_0 < 1 (منحنيات مغلقة حول (2kπ,0)(2k\pi, 0))، و E0=1E_0 = 1 (الفاصلة المارّة بالنقطتين (±π,0)(\pm\pi, 0))، و E0>1E_0 > 1 (منحنيات بيانية فوق xx: أي دورانات).
  3. برهن على أن موضع التوازن السفلي (0,0)(0,0) مستقر لكنه ليس مستقرًّا مقاربًا. (بليابونوف مع V=E+1V = E + 1؛ وعدم الاستقرار المقارب: إذ يحبس حفظُ الطاقة المداراتِ على منحنيات مستوية بعيدة عن المبدأ.)

الجزء الثاني — التذبذبات ودورها. نثبّت 1<E0<1-1 < E_0 < 1 ونكتب E0=cosaE_0 = -\cos a مع a(0,π)a \in \intoo0\pi (أي السعة).

  1. برهن على أن الحل الذي يحقق x(0)=ax(0) = a و y(0)=0y(0) = 0 يتذبذب: x(t)[a,a]x(t) \in \intcc{-a}a، ويكون المدار المنحني المغلق y2=2(cosxcosa)y^2 = 2(\cos x - \cos a). وبرّر أن الحل دوري: فالمدار منحنٍ متراص بلا مواضع توازن، ويُجرى بسرعة محدودة من الأسفل — اجعل هذا حجةً (إذ يعود الحل إلى نقطته الابتدائية في زمن منتهٍ، ثم تفرض الوحدانية الدورية).
  2. أثبت صيغة الدور

    T(a)=40a ⁣dx2(cosxcosa)T(a) = 4\int_0^{a}\frac{\dd x}{\sqrt{2(\cos x - \cos a)}}

    (على ربع مدار، y= ⁣dx ⁣dt>0y = \frac{\dd x}{\dd t} > 0 ثم افصل المتغيّرات؛ وبرّر التقارب المعتلّ عند x=ax = a).

  3. (التذبذبات الصغيرة) عوّض sinx2=sina2sinφ\sin\frac x2 = \sin\frac a2\,\sin\varphi وبرهن على

    T(a)=40π/2 ⁣dφ1k2sin2φ,k=sina2T(a) = 4\int_0^{\pi/2}\frac{\dd\varphi}{\sqrt{1 - k^2\sin^2\varphi}}, \qquad k = \sin\frac a2

    (وهو تكامل إهليلجي تام). واستنتج بالتقارب المهيمن أن T(a)2πT(a) \to 2\pi حين a0+a \to 0^+: أي النهاية التوافقية، المستقلة عن السعة — أي تساوي الأزمنة التقريبي عند غاليليو، مع تصحيحه المضبوط T(a)=2π(1+k24+O(k4))T(a) = 2\pi\bigl(1 + \frac{k^2}4 + O(k^4)\bigr) (بنشر المقدار المكامَل والمكاملة حدًّا حدًّا، مبرّرًا ذلك بالتقارب الناظمي).

  4. برهن على T(a)+T(a) \to +\infty حين aπa \to \pi^- (بحدّ المقدار المكامَل من الأسفل بجوار φ=π2\varphi = \frac\pi2 حين k1k \to 1، أو بتطبيق التقارب الرتيب): فباقتراب النواس من الفاصلة، يتباطأ دون حدّ.

الجزء الثالث — الفاصلة.

  1. من أجل E0=1E_0 = 1، y=2cosx2y = 2\cos\frac x2 على الفرع العلوي: افصل المتغيّرات وكامل لتجد الحل الصريح

    x(t)=4arctan(et)πx(t) = 4\arctan\bigl(\eu^{t}\bigr) - \pi

    (مع x(0)=0x(0) = 0 و y(0)=2y(0) = 2). وتحقق مباشرةً من أنه يحل معادلة النواس، واحسب نهاياته ونهايات y(t)y(t) حين t±t \to \pm\infty.

  2. اخلص: أن مدار الفاصلة يصل بين السرج (π,0)(-\pi, 0) (حين tt \to -\infty) والسرج (π,0)(\pi, 0) (حين t+t \to +\infty) لكنه لا يبلغ أيًّا منهما في زمن منتهٍ — وهو متسق مع الوحدانية (ولماذا كان بلوغُ سرجٍ في زمن منتهٍ سيناقض النتيجة 19.4؟).

الجزء الرابع — الدورانات، والصورة الكاملة.

  1. من أجل E0>1E_0 > 1: برهن على أن yy لا ينعدم أبدًا، وأن xx رتيبة تمامًا وشاملة مع x(t)±x(t) \to \pm\infty، وأن ty(t)t \mapsto y(t) دورية بالدور

    τ(E0)=ππ ⁣dx2(E0+cosx).\tau(E_0) = \int_{-\pi}^{\pi} \frac{\dd x}{\sqrt{2(E_0 + \cos x)}} .
  2. ركّب صورة الطور الكاملة (أي شكل الفصل) مع تبرير كامل لكل سمة، واكتب ملخصًا من عشرة أسطر للطريقة: الطاقة، ومجموعات المستويات، والتراص، والوحدانية — وكيف دخلت كل مبرهنة من مبرهنات الفصل. وأين احتجنا يومًا إلى صيغة للحل العام؟

الجزء الخامس — دالة الدور تحت المجهر.

  1. برهن على عزوم واليس

    Wn=0π/2sin2nφ ⁣dφ=π2(2n)!4n(n!)2W_n = \int_0^{\pi/2}\sin^{2n}\varphi\,\dd\varphi = \frac\pi2\cdot\frac{(2n)!}{4^n\,(n!)^2}

    بالتراجع (بالمكاملة بالتجزئة)، وانشر المقدار المكامَل في السؤال 6 بمتسلسلة ذات الحدين، وبرّر المكاملة حدًّا حدًّا لتحصل على المتسلسلة الكاملة

    T(a)=2πn0((2n)!4n(n!)2) ⁣2k2n=2π(1+k24+9k464+O(k6)),k=sina2.T(a) = 2\pi\sum_{n\geq0} \Bigl(\frac{(2n)!}{4^n(n!)^2}\Bigr)^{\!2}k^{2n} = 2\pi\Bigl(1 + \frac{k^2}4 + \frac{9k^4}{64} + O(k^6)\Bigr), \qquad k = \sin\frac a2 .
  2. حوّل إلى السعة:

    T(a)=2π(1+a216+11a43072+O(a6))T(a) = 2\pi\Bigl(1 + \frac{a^2}{16} + \frac{11\,a^4}{3072} + O(a^6)\Bigr)

    (بتعويض نشر sina2\sin\frac a2 والتجميع). فيخفق تساوي الأزمنة من الرتبة a2a^2، وقد صار الإخفاق الآن مكمَّمًا حتى الرتبة a4a^4.

  3. برهن على أن aT(a)a \mapsto T(a) متصلة ومتزايدة تمامًا على (0,π)\intoo0\pi، واخلص مع السؤالين 6 و7 إلى أن TT تقابل من (0,π)\intoo0\pi على (2π,+)\intoo{2\pi}{+\infty}: أي إن كل دور فوق حرج تحققه سعةٌ واحدة بالضبط.
  4. (حساب صانع الساعات) نواسٌ مضبوط عند سعة متلاشية يحفظ الزمن المثالي؛ برهن على أنه إذا شُغّل عند السعة aa تأخّر بالكسر a216+O(a4)\frac{a^2}{16} + O(a^4) من الزمن المثالي، واحسب الانحراف من أجل a=0.2a = 0.2 راديان: نحو 216216 ثانية في اليوم. (وخدّا هويغنس الدويريان والسعات الثابتة الصغيرة في الميزانات جوابان على هذا العدد.)
  5. عُد إلى دور الدوران τ\tau في السؤال 10: برهن على أن τ\tau متناقص تمامًا على (1,+)\intoo1{+\infty}، وأن τ(E0)+\tau(E_0) \to +\infty حين E01+E_0 \to 1^+ (بالتقارب الرتيب)، وأن 2E0τ(E0)2π\sqrt{2E_0}\,\tau(E_0) \to 2\pi حين E0+E_0 \to +\infty (بالتقارب المهيمن): أي إن الدوران السريع هو مقاربًا دورانٌ حر بالسرعة الزاوية 2E0\sqrt{2E_0}.

الجزء السادس — الطريقة مصدَّرةً: لوتكا–فولتيرا. وصفة النواس — تكاملٌ أول، ومنحنيات مستوية متراصة، ووحدانية — تحل نظامًا بيئيًّا. نثبّت α,β,γ,δ>0\alpha, \beta, \gamma, \delta > 0 وننظر، على الربع المفتوح Q=(0,+)×(0,+)Q = \intoo0{+\infty}\times\intoo0{+\infty}، في

x=x(αβy),y=y(δxγ)x' = x\,(\alpha - \beta y), \qquad y' = y\,(\delta x - \gamma)

(حيث xx الفريسة و yy المفترسات).

  1. برهن على أن QQ صامدة — إذ المحوران اتحادان لمدارات قابلة للحساب صراحةً، لا يمكن لأي حل أن يعبرها (النتيجة 19.4) — وعلى أن موضع التوازن الوحيد في QQ هو (x,y)=(γδ,αβ)(x_*, y_*) = \bigl(\frac\gamma\delta, \frac\alpha\beta\bigr).
  2. برهن على أن

    H(x,y)=δxγlnx+βyαlnyH(x, y) = \delta x - \gamma\ln x + \beta y - \alpha\ln y

    تكاملٌ أول، وعلى أن H=f(x)+g(y)H = f(x) + g(y) مع f,gf, g محدّبتين تمامًا وسويّتين على (0,+)\intoo0{+\infty} بأصغريتين عند xx_* و yy_*، واستنتج أن جميع الحلول الأعظمية في QQ شاملة.

  3. برهن على أنه من أجل h>h=H(x,y)h > h_* = H(x_*, y_*) تكون مجموعة المستوى {H=h}Q\{H = h\}\cap Q منحنيًا مغلقًا حول موضع التوازن: أي فرعان متصلان y±(x)y_\pm(x) فوق فترة متراصة [x,x+]x\intcc{x_-}{x_+} \ni x_*، ملتحمان عند الطرفين — وهو نظير بيضاويات النواس.
  4. برهن على أن كل مدار غير توازني في QQ دوري: أثبت الدوران عكس عقارب الساعة عبر المناطق الأربع التي يقطعها المستقيمان x=xx = x_* و y=yy = y_*، وحُدّ زمن عبور كل قوس بتكامل ذي شذوذ جذري تربيعي متقارب (كما في السؤال 5)، ثم أغلق بالوحدانية (كما في السؤال 4).
  5. (قانون فولتيرا في المتوسطات) إذا كان TT دورَ مدار كهذا، برهن على أن

    1T0Tx(t) ⁣dt=γδ,1T0Ty(t) ⁣dt=αβ:\frac1T\int_0^Tx(t)\,\dd t = \frac\gamma\delta, \qquad \frac1T\int_0^Ty(t)\,\dd t = \frac\alpha\beta :

    أي إن المتوسطات الزمنية تساوي قيم التوازن، مهما كانت السعة (بمكاملة (lnx)=αβy(\ln x)' = \alpha - \beta y على دور واحد).

  6. (مفارقة الصيد) اصطد النوعين بالمعدل ε(0,α)\varepsilon \in \intoo0\alpha: فتحتفظ الجملة بشكلها مع αε\alpha - \varepsilon و γ+ε\gamma + \varepsilon بدل α\alpha و γ\gamma. فماذا يحدث لمتوسطَي التعداد؟ واشرح ملاحظة دانكونا (1914–1918): أنه لما تناقص الصيد في الأدرياتيك أثناء الحرب، ازدادت نسبة المفترسات (أسماك القرش) في المصيد — ولماذا يفضّل الصيدُ المعتدل الفريسةَ.
  7. اكتب العبرة في عشرة أسطر: أي مبرهنات الفصل تشغّل كل خطوة، وما الذي يحلّ محل طاقة النواس، ولماذا لم تحتج أي من الجملتين — ولا تقبل — حلًّا ابتدائيًّا بصيغة مغلقة.
  8. (تضمين السرعة) في نظام الدوران E0>1E_0 > 1، برهن على أن y=xy = x' يتذبذب بين 2(E01)\sqrt{2(E_0 - 1)} (عند xπmod2πx \equiv \pi \bmod 2\pi) و 2(E0+1)\sqrt{2(E_0 + 1)} (عند x0x \equiv 0)، وعلى أن المتوسط الزمني للمقدار yy على دور واحد يساوي 2πτ(E0)\frac{2\pi}{\tau(E_0)} بالضبط، وعلى أن نسبة التضمين E0+1E011\sqrt{\frac{E_0 + 1}{E_0 - 1}} \to 1 حين E0E_0 \to \infty: أي إن الدوران السريع منتظمٌ مقاربًا.
  9. (رتابة دور الدوران) برهن على أن τ(E0)\tau(E_0) من الصنف C1\mathcal C^1 ومتناقص تمامًا على (1,+)\intoo1{+\infty} (بالاشتقاق تحت علامة التكامل، مع الهيمنة على كل [1+δ,)\intco{1 + \delta}\infty)، مع τ\tau \to \infty حين E01+E_0 \to 1^+ و τ0\tau \to 0 حين E0E_0 \to \infty. وركّب صورة التشعّب الكاملة للنواس على امتداد محور الطاقة: مواضع التوازن عند E0=1E_0 = -1، والتأرجحات بدور يتزايد من 2π2\pi إلى \infty على 1<E0<1-1 < E_0 < 1، والفاصلة عند E0=1E_0 = 1، والدورانات بدور يتناقص من \infty إلى 00 بعدها.
حل

حل المسألة 19.1.

1. E˙=yy+sinxx=ysinx+ysinx=0\dot E = yy' + \sin x\cdot x' = -y\sin x + y\sin x = 0: أي إن الطاقة تكامل أول. وعلى حل أعظمي، y2=2(E0+cosx)2(E0+1)y^2 = 2(E_0 + \cos x) \leq 2(E_0 + 1): فيكون yy محدودًا؛ عندئذٍ ينمو x(t)x(0)+tsupy\abs{x(t)} \leq \abs{x(0)} + t\sup\abs y خطيًّا على الأكثر: فعلى أي فترة زمنية منتهية يبقى المسار في متراصة من R2\R^2، ومنه تفرض المبرهنة 19.5(2) أن T±=±T_\pm = \pm\infty. ومواضع التوازن (kπ,0)(k\pi, 0)؛ والخطّية (0110)\bigl(\begin{smallmatrix}0&1\\ \mp1&0\end{smallmatrix}\bigr) مع cos(kπ)=1-\cos(k\pi) = \mp1: بقيمتين ذاتيتين ±i\pm\iu من أجل kk زوجي (من نمط المركز، وهو غير حاسم بذاته) و ±1\pm1 من أجل kk فردي (سرج).

2. ثبات EE على امتداد الحلول يحصر كل مسار في مجموعة مستوى {y2=2(E0+cosx)}\{y^2 = 2(E_0 + \cos x)\}. فمن أجل E0=1E_0 = -1: النقاط (2kπ,0)(2k\pi, 0) وحدها. ومن أجل 1<E0<1-1 < E_0 < 1: بكتابة E0=cosaE_0 = -\cos a، تكون المجموعة اتحادًا منفصلًا لمنحنيات مغلقة y=±2(cosxcosa)y = \pm\sqrt{2(\cos x - \cos a)} فوق x[2kπa,2kπ+a]x \in [2k\pi - a, 2k\pi + a]، واحدٍ حول كل موضع توازن مستقر. ومن أجل E0=1E_0 = 1: المنحنيات y=±2cosx2y = \pm2\cos\frac x2 الواصلة بين السروج المتتالية — أي الفاصلة — مع السروج نفسها. ومن أجل E0>1E_0 > 1: منحنيان بيانيان y=±2(E0+cosx)y = \pm\sqrt{2(E_0 + \cos x)}، معرَّفان من أجل كل xx، ولا يمسّان y=0y = 0 أبدًا.

3. ينعدم V=E+1=y22+(1cosx)V = E + 1 = \frac{y^2}2 + (1 - \cos x) عند (0,0)(0,0)، وهو موجب على جوار مثقوب (حيث x<2π\abs x < 2\pi)، و V˙=00\dot V = 0 \leq 0: فتعطي المبرهنة 19.11 الاستقرار. وليس مقاربًا: إذ يبقى الحل المارّ بالنقطة (a,0)(a, 0) (حيث 0<a0 < a صغير) على المنحني المستوي E=cosaE = -\cos a، ومسافته إلى المبدأ موجبة (إذ لا يلاقي المنحني محور السينات إلا عند ±a\pm a): أي φt(a,0)↛(0,0)\varphi_t(a, 0) \not\to (0,0).

4. على المنحني المستوي CaC_a: لا مواضع توازن (إذ y=0y = 0 يفرض x=±ax = \pm a مع sin(±a)0\sin(\pm a) \neq 0 من أجل 0<a<π0 < a < \pi)، ومنه للسرعة (y,sinx)\norm{(y, -\sin x)} أصغريةٌ موجبة mm على المتراصة CaC_a. ونتتبّع الحل من (a,0)(a, 0): ففي نصف المستوي السفلي x=y<0x' = y < 0، ومنه يتناقص xx من aa إلى a-a في زمن منتهٍ (إذ تتقارب تكاملات ربع الدور ونصفه: بتحليل السؤال 5)، فيبلغ (a,0)(-a, 0)؛ وبالتماثل (x,y)(x,y)(x, y) \mapsto (x, -y)، ttt \mapsto -t للمعادلة، يُجرى النصف العلوي رجوعًا في الزمن نفسه T2\frac T2: فيرجع الحل إلى (a,0)(a, 0) عند الزمن TT. ثم تنتشر الوحدانية (النتيجة 19.4): أي x(t+T)=x(t)x(t + T) = x(t) من أجل كل tt: فهو دوري.

5. على الفرع الذي يكون فيه y>0y > 0: y=2(cosxcosa)y = \sqrt{2(\cos x - \cos a)} و  ⁣dt= ⁣dxy\dd t = \frac{\dd x}{y}؛ وبمكاملة xx من a-a إلى aa نجد نصف الدور، ويناصف التماثل xxx \mapsto -x التكاملَ مرة أخرى:

T(a)=40a ⁣dx2(cosxcosa).T(a) = 4\int_0^a\frac{\dd x}{\sqrt{2(\cos x - \cos a)}} .

وأما التقارب عند x=ax = a^-: فإن cosxcosa=sin(a)(ax)+O((ax)2)\cos x - \cos a = \sin(a)(a - x) + O((a-x)^2) مع sina>0\sin a > 0: ومنه يسلك المقدار المكامَل مثل (2sina(ax))1/2\bigl(2\sin a\,(a - x)\bigr)^{-1/2}، وهو قابل للمكاملة.

6. مع sinx2=ksinφ\sin\frac x2 = k\sin\varphi و k=sina2k = \sin\frac a2: cosxcosa=2(k2sin2x2)=2k2cos2φ\cos x - \cos a = 2(k^2 - \sin^2\frac x2) = 2k^2\cos^2\varphi و 12cosx2 ⁣dx=kcosφ ⁣dφ\frac12\cos\frac x2\,\dd x = k\cos\varphi\,\dd\varphi، ومنه

T(a)=40π/2 ⁣dφ1k2sin2φ.T(a) = 4\int_0^{\pi/2} \frac{\dd\varphi}{\sqrt{1 - k^2\sin^2\varphi}} .

وحين a0+a \to 0^+ يكون k0k \to 0: فمن أجل kk0<1k \leq k_0 < 1 يُهيمَن على المقدار المكامَل بالمقدار (1k02sin2φ)1/2(1 - k_0^2\sin^2\varphi)^{-1/2}، وهو متصل على [0,π/2]\intcc0{\pi/2}: فيعطي التقارب المهيمن T4π2=2πT \to 4\cdot\frac\pi2 = 2\pi. وبنشر (1u)1/2=1+u2+3u28+(1 - u)^{-1/2} = 1 + \frac u2 + \frac{3u^2}8 + \cdots مع u=k2sin2φu = k^2\sin^2\varphi (بتقارب ناظمي من أجل k<1k < 1) ومع 0π/2sin2=π4\int_0^{\pi/2}\sin^2 = \frac\pi4:

T(a)=2π(1+k24+O(k4)):T(a) = 2\pi\Bigl(1 + \frac{k^2}{4} + O(k^4)\Bigr) :

فلا يصح تساوي الأزمنة إلا من الرتبة الأولى؛ ويتزايد الدور مع السعة.

7. حين k1k \uparrow 1 تتزايد المقادير المكامَلة إلى (1sin2φ)1/2=1cosφ(1 - \sin^2\varphi)^{-1/2} = \frac1{\cos\varphi}، وتكامله متباعد: فبالتقارب الرتيب، T(a)+T(a) \to +\infty حين aπa \to \pi^-.

8. على الفرع y=2cosx2y = 2\cos\frac x2 (x<π\abs x < \pi):  ⁣dx2cos(x/2)= ⁣dt\frac{\dd x}{2\cos(x/2)} = \dd t؛ ومع u=x2u = \frac x2 و  ⁣ducosu=lntan(u2+π4)\int\frac{\dd u}{\cos u} = \ln\tan\bigl(\frac u2 + \frac\pi4\bigr)، يكون t=lntan(x4+π4)t = \ln\tan\bigl(\frac x4 + \frac\pi4\bigr)، أي

x(t)=4arctan(et)π,y(t)=x(t)=4et1+e2t=2cosht.x(t) = 4\arctan(\eu^t) - \pi, \qquad y(t) = x'(t) = \frac{4\eu^t}{1 + \eu^{2t}} = \frac{2}{\cosh t} .

وأما التحقق بالطاقة: فمع θ=arctanet\theta = \arctan\eu^t، sin2θ=1cosht\sin2\theta = \frac1{\cosh t}، ومنه cosx=cos4θ=1+2cosh2t\cos x = -\cos4\theta = -1 + \frac{2}{\cosh^2t} و E=y22cosx=2cosh2t+12cosh2t=1E = \frac{y^2}2 - \cos x = \frac2{\cosh^2t} + 1 - \frac2{\cosh^2t} = 1: فيقع المسار على الفاصلة، ويعيد اشتقاق y2=2(1+cosx)y^2 = 2(1 + \cos x) حيث y>0y > 0 إنتاجَ y=sinxy' = -\sin x. والنهايات: x±πx \to \pm\pi و y0y \to 0 حين t±t \to \pm\infty.

9. يؤول المدار إلى السرج (π,0)(\pi, 0) أمامًا وإلى (π,0)(-\pi, 0) خلفًا لكنه لا يصل أبدًا: فلو بلغ (π,0)(\pi, 0) عند زمن منتهٍ tt^*، لمرّ حلّان أعظميان متمايزان — حل الفاصلة والحل الثابت عند السرج — بالنقطة نفسها (t,(π,0))(t^*, (\pi, 0))، وهو يناقض النتيجة 19.4. فلا تُقارَب السروج إلا مقاربًا.

10. من أجل E0>1E_0 > 1: y2=2(E0+cosx)2(E01)>0y^2 = 2(E_0 + \cos x) \geq 2(E_0 - 1) > 0: فيحتفظ yy بإشارته، و x=y2(E01)\abs{x'} = \abs y \geq \sqrt{2(E_0 - 1)}: فتكون xx رتيبة تمامًا وشاملة (السؤال 1) مع x(t)±x(t) \to \pm\infty. وبما أن y(t)=±2(E0+cosx(t))y(t) = \pm\sqrt{2(E_0 + \cos x(t))} و cos\cos دورية بالدور 2π2\pi، يرجع yy إلى قيمته كلما تقدّم xx بمقدار 2π2\pi؛ والزمن اللازم هو

τ(E0)=x0x0+2π ⁣dx2(E0+cosx)=ππ ⁣dx2(E0+cosx)\tau(E_0) = \int_{x_0}^{x_0 + 2\pi}\frac{\dd x}{\sqrt{2(E_0 + \cos x)}} = \int_{-\pi}^{\pi}\frac{\dd x}{\sqrt{2(E_0 + \cos x)}}

(بالتعويض؛ وبالدورية): فيكون yy دوريًّا بالدور τ\tau — أي إن النواس يدور بمعدل دوران ثابت مقاربًا 2π/τ2E02\pi/\tau \approx \sqrt{2E_0} من أجل الطاقات الكبيرة.

11. الطريقة، بالترتيب: الطاقة (E˙=0\dot E = 0) تختزل الانسياب ذا البُعدين إلى منحنيات مستوية ذات بُعد واحد؛ وحدّية yy على كل مستوى مع الهروب من المتراصات تعطي الوجود الشامل؛ وتراص المستويات المغلقة يعطي حدود السرعة ومن ثَمّ الدورية؛ والوحدانية تحوّل العودة الأولى إلى دورية مضبوطة، وتمنع الوصول إلى السروج في زمن منتهٍ، وتفصل أنماط المدارات؛ والخطّية وليابونوف تصنّفان مواضع التوازن؛ وتكامل الدور يُحلَّل بمبرهنات التقارب في الفصل 10. ولم نمتلك في أي لحظة — ولم نحتج — إلى حل عام بصيغة مغلقة: إذ استخرجت النظرية النوعية كل سمة من سمات الحركة من المعادلة نفسها.

12. بالتجزئة: Wn=0π/2sin2n1φsinφ ⁣dφ=(2n1)0π/2sin2n2φcos2φ ⁣dφ=(2n1)(Wn1Wn)W_n = \int_0^{\pi/2}\sin^{2n-1}\varphi \cdot\sin\varphi\,\dd\varphi = (2n-1)\int_0^{\pi/2} \sin^{2n-2}\varphi\cos^2\varphi\,\dd\varphi = (2n-1)(W_{n-1} - W_n)، ومنه Wn=2n12nWn1W_n = \frac{2n-1}{2n}W_{n-1}؛ ومع W0=π2W_0 = \frac\pi2 و j=1n2j12j=(2n)!4n(n!)2\prod_{j=1}^n\frac{2j-1}{2j} = \frac{(2n)!}{4^n(n!)^2} (بشطر (2n)!=2nn!(2j1)(2n)! = 2^nn!\prod(2j-1))، يعطي التراجع القيمةَ المعروضة. ومتسلسلة ذات الحدين: (1u)1/2=n0cnun(1 - u)^{-1/2} = \sum_{n\geq0}c_nu^n مع cn=(2n)!4n(n!)2(0,1]c_n = \frac{(2n)!}{4^n(n!)^2} \in \intoc01، ونصف القطر 11. ومن أجل kk0<1k \leq k_0 < 1 و u=k2sin2φu = k^2\sin^2\varphi، تتقارب المتسلسلة cnk2nsin2nφ\sum c_nk^{2n}\sin^{2n}\varphi تقاربًا ناظميًّا في φ\varphi (cnk2nk02nc_nk^{2n} \leq k_0^{2n})، ومنه تكون المكاملة حدًّا حدًّا في صيغة السؤال 6 مشروعة:

T(a)=4n0cnk2nWn=4π2n0cn2k2n=2πn0((2n)!4n(n!)2) ⁣2k2n.T(a) = 4\sum_{n\geq0}c_nk^{2n}W_n = 4\cdot\frac\pi2\sum_{n\geq0}c_n^2\,k^{2n} = 2\pi\sum_{n\geq0} \Bigl(\frac{(2n)!}{4^n(n!)^2}\Bigr)^{\!2}k^{2n} .

ومع c0=1c_0 = 1 و c1=12c_1 = \frac12 و c2=38c_2 = \frac38: T(a)=2π(1+k24+9k464+O(k6))T(a) = 2\pi\bigl(1 + \frac{k^2}4 + \frac{9k^4}{64} + O(k^6)\bigr)، والباقي منتظم من أجل kk0k \leq k_0 (إذ يُهيمَن على الذيل بمتسلسلة هندسية).

13. k=sina2=a2a348+O(a5)k = \sin\frac a2 = \frac a2 - \frac{a^3}{48} + O(a^5)، ومنه

k2=a24a448+O(a6),k4=a416+O(a6),k^2 = \frac{a^2}4 - \frac{a^4}{48} + O(a^6), \qquad k^4 = \frac{a^4}{16} + O(a^6) ,

و

T(a)2π=1+14(a24a448)+964a416+O(a6)=1+a216+11a43072+O(a6),\frac{T(a)}{2\pi} = 1 + \frac14\Bigl(\frac{a^2}4 - \frac{a^4}{48}\Bigr) + \frac9{64}\cdot\frac{a^4}{16} + O(a^6) = 1 + \frac{a^2}{16} + \frac{11\,a^4}{3072} + O(a^6) ,

لأن 1192+91024=16+273072=113072-\frac1{192} + \frac9{1024} = \frac{-16 + 27}{3072} = \frac{11}{3072}.

14. في الصيغة الإهليلجية للسؤال 6، يكون ak=sina2a \mapsto k = \sin\frac a2 تقابلًا متصلًا متزايدًا تمامًا من (0,π)\intoo0\pi على (0,1)\intoo01، ومن أجل كل φ(0,π/2]\varphi \in \intoc0{\pi/2} يكون المقدار المكامَل (1k2sin2φ)1/2(1 - k^2\sin^2\varphi)^{-1/2} متزايدًا تمامًا في kk: فيكون TT متزايدًا تمامًا. وأما الاتصال: فعلى kk0<1k \leq k_0 < 1 يُهيمَن على المقدار المكامَل بالدالة المتصلة (1k02sin2φ)1/2(1 - k_0^2\sin^2\varphi)^{-1/2}، ومنه ينطبق التقارب المهيمن على امتداد kkk' \to k. ومع النهايتين T2πT \to 2\pi حين a0+a \to 0^+ (السؤال 6) و T+T \to +\infty حين aπa \to \pi^- (السؤال 7)، تجعل الرتابة الأكيدة ومبرهنة القيم الوسطى TT تقابلًا من (0,π)\intoo0\pi على (2π,+)\intoo{2\pi}{+\infty}.

15. تعدّ الساعة التأرجحات؛ وحين تُضبط عند سعة متلاشية، تقيّد الدور التوافقي 2π2\pi لكل تأرجحة (بوحدة زمن النواس). وحين تُشغَّل عند السعة aa، يكون الدور الحقيقي T(a)=2π(1+a216+O(a4))T(a) = 2\pi\bigl(1 + \frac{a^2}{16} + O(a^4)\bigr): فتقيّد الساعة 2π2\pi بينما ينقضي فعلًا T(a)T(a)، ومنه تتأخّر بالكسر

T(a)2πT(a)=a216+O(a4).\frac{T(a) - 2\pi}{T(a)} = \frac{a^2}{16} + O(a^4) .

ومن أجل a=0.2a = 0.2 راديان (نحو 11.511.5 درجة): a216=0.0416=1400\frac{a^2}{16} = \frac{0.04}{16} = \frac1{400}، ولليوم 8640086400 ثانية: فتفقد الساعة 86400/400=21686400/400 = 216 ثانية — أي نحو ثلاث دقائق ونصف — في اليوم. ومن هنا العلاجان التاريخيان: فرضُ سعة ثابتة ضئيلة (بالميزان)، أو ثنيُ القيد بحيث يكون الدور مستقلًّا عن السعة تمامًا (خدّا هويغنس الدويريان، 1657).

16. في τ(E0)=ππ ⁣dx2(E0+cosx)\tau(E_0) = \int_{-\pi}^{\pi}\frac{\dd x}{\sqrt{2(E_0 + \cos x)}} يكون المقدار المكامَل، من أجل كل xx مثبَّت، متناقصًا تمامًا في E0E_0: فيكون τ\tau متناقصًا تمامًا. وحين E01E_0 \downarrow 1 تتزايد المقادير المكامَلة نقطةً نقطة إلى (2(1+cosx))1/2=12cosx2\bigl(2(1 + \cos x)\bigr)^{-1/2} = \frac1{2\abs{\cos\frac x2}}، وتكامله على (π,π)\intoo{-\pi}\pi متباعد (إذ ينعدم cosx2\cos\frac x2 من الرتبة الأولى عند ±π\pm\pi): فيعطي التقارب الرتيب τ(E0)+\tau(E_0) \to +\infty. وحين E0+E_0 \to +\infty:

2E0τ(E0)=ππ ⁣dx1+cosx/E02π\sqrt{2E_0}\,\tau(E_0) = \int_{-\pi}^{\pi} \frac{\dd x}{\sqrt{1 + \cos x/E_0}} \longrightarrow 2\pi

بالتقارب المهيمن (إذ من أجل E02E_0 \geq 2 يكون المقدار المكامَل على الأكثر 2\sqrt2). ومنه τ2π/2E0\tau \approx 2\pi/\sqrt{2E_0}: أي إن الدورة الواحدة تستغرق زمن الدوران الحر بالسرعة 2E0\sqrt{2E_0}، وقد اختُزل الكمون إلى تموّج — وهو يوافق معدل الدوران في السؤال 10.

17. يحمل المحوران الحلين الصريحين t(x0eαt,0)t \mapsto (x_0\eu^{\alpha t}, 0) و t(0,y0eγt)t \mapsto (0, y_0\eu^{-\gamma t})، مع موضع التوازن (0,0)(0, 0): فهما اتحادان لمدارات. والحقل C1\mathcal C^1، ومنه ليبشيتزي محليًّا؛ وحلٌّ ينطلق من QQ ويلمس محورًا كان سيمرّ بنقطة من أحد تلك المدارات، فيتوافق معه حسب النتيجة 19.4 — وهو مستحيل، إذ يعيش أحدهما على المحور والآخر لا. ومنه فإن QQ صامدة في اتجاهَي الزمن. وأما مواضع التوازن في QQ: فإن x>0x > 0 يفرض αβy=0\alpha - \beta y = 0 و y>0y > 0 يفرض δxγ=0\delta x - \gamma = 0: أي النقطة الوحيدة (x,y)=(γδ,αβ)(x_*, y_*) = \bigl(\frac\gamma\delta, \frac\alpha\beta\bigr).

18. على امتداد حل،

H˙=(δγx)x+(βαy)y=(δxγ)(αβy)+(βyα)(δxγ)=0.\dot H = \Bigl(\delta - \frac\gamma x\Bigr)x' + \Bigl(\beta - \frac\alpha y\Bigr)y' = (\delta x - \gamma)(\alpha - \beta y) + (\beta y - \alpha)(\delta x - \gamma) = 0 .

وللمقدار f(x)=δxγlnxf(x) = \delta x - \gamma\ln x المشتقةُ f=γ/x2>0f'' = \gamma/x^2 > 0، ولا تنعدم ff' إلا عند xx_*، و f+f \to +\infty عند 0+0^+ وعند ++\infty معًا: أي محدّبة تمامًا وسويّة، بأصغرية f(x)f(x_*)؛ وبالمثل g(y)=βyαlnyg(y) = \beta y - \alpha\ln y، بأصغرية g(y)g(y_*). ومنه HhH \geq h_* مع التساوي عند (x,y)(x_*, y_*) وحدها، وتكون كل مجموعة مستوى أدنى {Hh}Q\{H \leq h\}\cap Q متراصة: إذ يحصر f(x)hg(y)f(x) \leq h - g(y_*) المتغيّرَ xx في فترة متراصة من (0,+)\intoo0{+\infty} بالسويّة، وبالمثل yy، والمجموعة مغلقة في R2\R^2 لأن H+H \to +\infty عند حافة QQ. ويبقى الحل الأعظمي على مجموعة مستواه المتراصة، ومنه لا يمكنه أن يغادر كل متراصة في زمن منتهٍ: فتجعله المبرهنة 19.5 شاملًا.

19. نثبّت h>hh > h_* ونضع c=hg(y)>f(x)c = h - g(y_*) > f(x_*). وبما أن ff تتناقص تمامًا من ++\infty إلى f(x)f(x_*) على (0,x]\intoc0{x_*} وتتزايد تمامًا راجعةً إلى ++\infty على [x,+)\intco{x_*}{+\infty}، يكون للمعادلة f(x)=cf(x) = c جذران بالضبط x<x<x+x_- < x_* < x_+، مع {fc}=[x,x+]\{f \leq c\} = \intcc{x_-}{x_+}. ومن أجل x(x,x+)x \in \intoo{x_-}{x_+}: للمعادلة g(y)=hf(x)>g(y)g(y) = h - f(x) > g(y_*) جذران بالضبط y(x)<y<y+(x)y_-(x) < y_* < y_+(x)، متصلان في xx (بوصفهما مقلوبَي التضييقين المتصلين الرتيبين تمامًا للدالة gg على جانبَي yy_*)، مع y±(x)yy_\pm(x) \to y_* حين xx±x \to x_\pm؛ وعند x=x±x = x_\pm يكون الحل الوحيد y=yy = y_*. ومنه فإن {H=h}Q\{H = h\}\cap Q اتحادُ منحنيَي y+y_+ و yy_- فوق [x,x+]\intcc{x_-}{x_+}، ملتحمين عند (x±,y)(x_\pm, y_*): أي منحنٍ مغلق حول (x,y)(x_*, y_*) — وهو نظير بيضاويات النواس.

20. لتكن Ch={H=h}QC_h = \{H = h\}\cap Q مع h>hh > h_*: فموضع التوازن الوحيد في QQ خارج ChC_h، ومنه لا ينعدم الحقل عليه أبدًا. وأما الإشارات: x=βx(yy)x' = \beta x(y_* - y)، y=δy(xx)y' = \delta y(x - x_*): فتسير الحركة يمينًا تحت المستقيم y=yy = y_*، وصعودًا يمين x=xx = x_*، ويسارًا فوقه، ونزولًا على اليسار — أي دوران عكس عقارب الساعة. ونتتبّع الحل من نقطة (x0,y(x0))(x_0, y_-(x_0)) من الفرع السفلي المفتوح، حيث x>0x' > 0: فيكون زمن بلوغ الزاوية اليمنى B=(x+,y)B = (x_+, y_*)

x0x+ ⁣dxβx(yy(x)).\int_{x_0}^{x_+}\frac{\dd x}{\beta x\,\bigl(y_* - y_-(x)\bigr)} .

وبجوار x+x_+ نختار x(x,x+)x' \in \intoo{x_*}{x_+}؛ فمن أجل x[x,x+]x \in \intcc{x'}{x_+}، f(x+)f(x)f(x)(x+x)f(x_+) - f(x) \geq f'(x')\,(x_+ - x) (إذ ff' متزايدة وموجبة بعد xx_*)، بينما تعطي علاقة المستوى ومتراجحة تايلور أن g(y(x))g(y)12(maxg)(yy(x))2g(y_-(x)) - g(y_*) \leq \frac12\,\bigl(\max g''\bigr)\,(y_* - y_-(x))^2 على مجال yy المتراص للمنحني ChC_h: ومنه yy(x)cx+xy_* - y_-(x) \geq c\,\sqrt{x_+ - x} مع c>0c > 0، ويكون المقدار المكامَل O((x+x)1/2)O\bigl((x_+ - x)^{-1/2}\bigr): أي قابلًا للمكاملة — وهو تقارب السؤال 5، منقولًا. وفي غير ذلك من الفرع يكون المقدار المكامَل متصلًا. ومنه تُبلَغ BB في زمن منتهٍ؛ وهناك y=δy(x+x)>0y' = \delta y_*(x_+ - x_*) > 0، فيدخل المدار المنطقة x>xx > x_* و y>yy > y_*، ويصعد إلى الزاوية العليا (x,y+)(x_*, y_+) بالتقدير المتناظر (بتبادل دورَي ff و gg)، وهكذا حول الأقواس الأربعة: فبعد زمن منتهٍ T>0T > 0 يرجع الحل إلى نقطة انطلاقه. وحسب النتيجة 19.4 يكون دوريًّا بالدور TT — أي حجة السؤال 4 حرفيًّا.

21. على مدار دوري بالدور TT في QQ، تكون tlnx(t)t \mapsto \ln x(t) من الصنف C1\mathcal C^1 ودورية بالدور TT، ومنه

0=0T(lnx) ⁣dt=0T(αβy) ⁣dt=αTβ0Ty ⁣dt,0 = \int_0^T(\ln x)'\,\dd t = \int_0^T(\alpha - \beta y)\,\dd t = \alpha T - \beta\int_0^Ty\,\dd t ,

فنجد 1T0Ty=αβ\frac1T\int_0^Ty = \frac\alpha\beta؛ وبالمثل يعطي 0=0T(lny)=δ0Tx ⁣dtγT0 = \int_0^T(\ln y)' = \delta\int_0^Tx\,\dd t - \gamma T أن 1T0Tx=γδ\frac1T\int_0^Tx = \frac\gamma\delta. فالمتوسطات الزمنية هي قيم التوازن، من أجل كل مدار مهما كانت السعة — أي قانون حفظ لم يضعه أحد باليد.

22. مع الصيد تعود الجملة من نمط لوتكا–فولتيرا، بالوسائط αε\alpha - \varepsilon و β\beta و γ+ε\gamma + \varepsilon و δ\delta (ويستمر موضع التوازن الداخلي لأن ε<α\varepsilon < \alpha). وبتطبيق السؤال 21 على الجملة الجديدة:

xˉ=γ+εδ  (متوسط الفريسة يرتفع),yˉ=αεβ  (متوسط المفترس ينخفض):\bar x = \frac{\gamma + \varepsilon}{\delta} \ \ (\text{متوسط الفريسة يرتفع}), \qquad \bar y = \frac{\alpha - \varepsilon}{\beta} \ \ (\text{متوسط المفترس ينخفض}) :

فيزيح الصيدُ العشوائي التوازنَ لصالح الفريسة. وتقرأ معطيات دانكونا ذلك بالعكس: فقد قلّصت الحرب الصيد، فانخفض ε\varepsilon، ومنه ارتفع متوسط المفترس (αε)/β(\alpha - \varepsilon)/\beta وانخفض متوسط الفريسة — أي نسبة أكبر من أسماك القرش في المصيد، وهو بالضبط ما سجّلته أسواق سمك الأدرياتيك. وهذا مبدأ فولتيرا، وهو الآلية نفسها وراء مفارقات المبيدات: فقتلُ المستويين الغذائيين معًا يفيد المستوى المأكول.

23. الوصفة، في الحالتين: (أ) تكاملٌ أول — EE من أجل النواس و HH هنا، يُوجد بفصل  ⁣dy/ ⁣dx\dd y/\dd x — يطوي المستوي إلى منحنيات؛ (ب) والسويّة وتراص مجموعات المستويات يعطيان الوجود الشامل عبر المبرهنة 19.5؛ (ج) وهندسة المستويات — بيضاويات، من شكل cos\cos هناك ومن التحدّب الأكيد للدالتين ff و gg هنا — تُقرأ من التكامل لا من الانسياب؛ (د) وحقلٌ غير منعدم على بيضاوي متراص مع شذوذات زاوية قابلة للمكاملة يفرض زمن عودة منتهيًا؛ (ه) والوحدانية (النتيجة 19.4) تحوّل العودة إلى دورية وتمنع الوصول إلى مواضع التوازن في زمن منتهٍ؛ (و) والعوائد — نشرات الدور، وقوانين المتوسطات — تأتي من مبرهنات التقارب مطبَّقةً على التكاملات الناتجة. ولا x=sinxx'' = -\sin x ولا لوتكا–فولتيرا يقبل حلًّا ابتدائيًّا بصيغة مغلقة (تكاملات إهليلجية في حالة، ومنحنيات مستوية متسامية في الأخرى)، ولم يُحتَج إلى ذلك في أي لحظة: إذ سلّمت المعادلةُ نفسها، مستجوَبةً نوعيًّا، الحركةَ بكاملها.

24. على مستوى طاقة E0>1E_0 > 1: y2=2(E0+cosx)>0y^2 = 2(E_0 + \cos x) > 0: فتقع القيم القصوى للمقدار y2y^2 عند cosx=±1\cos x = \pm1، فتعطي الحدود المذكورة، وهي مبلوغة عند x0,πx \equiv 0, \pi. والمتوسط الزمني على دور واحد τ=τ(E0)\tau = \tau(E_0): إذ يتقدّم xx بمقدار 2π2\pi بالضبط، ومنه

1τ0τy ⁣dt=x(τ)x(0)τ=2πτ.\frac1\tau\int_0^\tau y\,\dd t = \frac{x(\tau) - x(0)}{\tau} = \frac{2\pi}\tau .

ونسبة السرعتين القصويين هي E0+1E01=1+O(E01)1\sqrt{\frac{E_0 + 1}{E_0 - 1}} = 1 + O(E_0^{-1}) \to 1: أي إنه عند الطاقة العالية يكون تموّج الكمون ±1\pm1 مهملًا أمام E0E_0، فيدور النواس دورانًا شبه منتظم — فتستوي لوحة الغسيل.

25. على E01+δE_0 \geq 1 + \delta، يُهيمَن على المقدار المكامَل في τ(E0)=ππ ⁣dx2(E0+cosx)\tau(E_0) = \int_{-\pi}^{\pi}\frac{\dd x}{\sqrt{2(E_0 + \cos x)}} بالمقدار (2δ)1/2(2\delta)^{-1/2}، وعلى مشتقته بالنسبة إلى E0E_0

E012(E0+cosx)=1(2(E0+cosx))3/2\partial_{E_0}\frac1{\sqrt{2(E_0 + \cos x)}} = -\frac{1}{\bigl(2(E_0 + \cos x)\bigr)^{3/2}}

بالمقدار (2δ)3/2(2\delta)^{-3/2}، وكلاهما قابل للمكاملة على المتراصة [π,π]\intcc{-\pi}\pi: فينطبق الاشتقاق تحت علامة التكامل (المبرهنة 10.15) ويعطي τ(E0)<0\tau'(E_0) < 0 (إذ المقدار المكامَل سالب تمامًا): أي متناقص تمامًا ومن الصنف C1\mathcal C^1. وأما النهايات: فحين E0E_0 \to \infty، τ2π2(E01)0\tau \leq \frac{2\pi}{\sqrt{2(E_0 - 1)}} \to 0؛ وحين E01+E_0 \to 1^+: نكتب E0+cosx=(E01)+2cos2x2E_0 + \cos x = (E_0 - 1) + 2\cos^2\frac x2؛ فيتزايد المقدار المكامَل حين يتناقص E0E_0، ومنه بالتقارب الرتيب

τ(E0)    ππ ⁣dx2cosx2=+,\tau(E_0) \;\nearrow\; \int_{-\pi}^{\pi}\frac{\dd x}{2\,\abs{\cos\frac x2}} = +\infty,

والتكامل الحدّي متباعد عند x=±πx = \pm\pi (إذ هناك cosx2xπ2\abs{\cos\frac x2} \sim \frac{\abs{x \mp \pi}}2، وهو 1\frac1{\abs\cdot} غير قابل للمكاملة): فينفجر الدور باقتراب الفاصلة، موافقًا T(a)T(a) \to \infty في الجزء الثاني من جهة التأرجح. ويُقرأ محور الطاقة: سكونٌ عند E0=1E_0 = -1؛ وتأرجحات، 2π2\pi \nearrow \infty، على (1,1)\intoo{-1}1؛ والفاصلة البطيئة لانهائيًّا عند E0=1E_0 = 1؛ ودورانات، 0\infty \searrow 0، بعدها. تكاملٌ واحد، وحياة النواس بأسرها.

المصطلحات المعرَّفة في هذا الفصل

عرض كل المصطلحات (395) في المسرد