गणित · किताब 5 · स्नातक वर्ष 3

विश्वविद्यालय गणित — स्नातक वर्ष 3

विश्वविद्यालय गणित — स्नातक वर्ष 3 · स्नातक वर्ष 3

19साधारण अवकल समीकरण

दूसरे वर्ष ने रैखिक अवकल समीकरण हल किए थे और कोशी–लिप्शिट्ज़ प्रमेय का कथन दिया था; यह अध्याय उसे सिद्ध करता है — और वह भी दो बार: बानाख स्थिर बिंदु से अस्तित्व और अद्वितीयता, तथा महत्तम हलों के सिद्धांत और संहतों-से-पलायन प्रमेय से वैश्विक संरचना। इसके बाद रैखिक सिद्धांत ईमानदार नींव पर फिर से खड़ा किया जाता है (विलायक, व्रोन्स्कियन, आव्यूह चरघातांकी, दुआमेल), और अध्याय का उत्तरार्ध गुणात्मक सिद्धांत खोलता है — प्रवाह, साम्यावस्थाएँ, ल्यापुनोव फलन, और रैखिकीकरण से स्थायित्व: अर्थात् उन हलों को समझना जिन्हें कोई कभी परिकलित नहीं करेगा। सप्ताहांत समस्या में लोलक चिरंतन अध्ययन-विषय है। सर्वत्र UR×RdU \subseteq \R\times\R^d विवृत है और f ⁣:URdf \colon U \to \R^d संतत; x=f(t,x)x' = f(t, x) का हल ऐसा C1\mathcal C^1 प्रतिचित्रण x ⁣:IRdx \colon I \to \R^d (II कोई अंतराल) है जिसका आलेख UU में हो और जो समीकरण संतुष्ट करे।

19.1 कोशी–लिप्शिट्ज़

परिभाषा 19.1

ff xx में स्थानीय लिप्शिट्ज़ कहलाता है यदि UU के प्रत्येक बिंदु का कोई प्रतिवेश VV और कोई नियतांक LL ऐसा हो कि (t,x1),(t,x2)V(t, x_1), (t, x_2) \in V के लिए f(t,x1)f(t,x2)Lx1x2\norm{f(t, x_1) - f(t, x_2)} \leq L\norm{x_1 - x_2}। यदि 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] पर ठीक एक हल है।

उपपत्ति. Q=[t0a,t0+a]×Bˉ(x0,b)UQ = [t_0 - a, t_0 + a]\times\bar B(x_0, b) \subseteq U वाला a,b>0a, b > 0 चुनिए, जिस पर fM\norm f \leq M और ff xx में LL-लिप्शिट्ज़ हो। कोई 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\ I \text{ पर}\},

जो बानाख समष्टि (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 पर कोई भी हल वहीं रहता है (जब तक आलेख QQ में रहता है तब तक x(t)x0Mtt0b\norm{x(t) - x_0} \leq M\abs{t - t_0} \leq b, जो एक सांतत्य तर्क है): अतः 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। तब II पर u(t)αeβtt0u(t) \leq \alpha\,\eu^{\beta\abs{t - t_0}}

उपपत्ति. 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 xx में LL-लिप्शिट्ज़ है, तो

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(T+ε,T+)t \in \intoo{T_+ - \varepsilon}{T_+} के लिए (t,x(t))K(t, x(t)) \notin K (और TT_- पर सममित रूप से): अर्थात् किसी महत्तम हल का आलेख अंततः UU के प्रत्येक संहत उपसमुच्चय को छोड़ देता है। विशेष रूप से U=R×RdU = \R\times\R^d और T+<+T_+ < +\infty के लिए: tT+t \to T_+^- पर x(t)+\norm{x(t)} \to +\infty (विस्फोट)।

उपपत्ति. (1) मान लीजिए S\mathcal S (t0,x0)(t_0, x_0) से होकर जाने वाले समस्त हलों का समुच्चय है; उपप्रमेय 19.4 से कोई भी दो अपने अंतरालों के प्रतिच्छेदन पर सहमत होते हैं, अतः वे जुड़ जाते हैं: J=ySIyJ = \bigcup_{y \in \mathcal S}I_y पर tt पर परिभाषित किसी भी yy के लिए x(t)=y(t)x(t) = y(t) रखिए: यह एक सुपरिभाषित हल है, जो स्पष्टतः महत्तम और अद्वितीय है। JJ विवृत है: किसी सिरे पर परिभाषित हल को उसी सिरे पर प्रमेय 19.2 से आगे बढ़ाया जा सकता।

(2) मान लीजिए दावा T+T_+ पर विफल होता है: तब (tn,x(tn))K(t_n, x(t_n)) \in K वाले tnT+t_n \to T_+ हैं; ध्यान दीजिए कि इससे T+<T_+ < \infty अनिवार्य हो जाता है, अथवा यदि T+=T_+ = \infty हो तो सिद्ध करने को कुछ है ही नहीं (KK समय में परिबद्ध है)। अतः मान लीजिए T+<T_+ < \infty। संहतता: KK के किसी प्रतिवेश के समस्त कोशी आँकड़ों के लिए एकसमान नियतांक M,L,a,bM, L, a, b काम करते हैं — मूर्त रूप में, KK को स्थानीय प्रमेय की उपपत्ति जैसे परिमित संख्या के बक्सों QiQ_i से आच्छादित कीजिए और मान लीजिए T>0T^* > 0 संगत अस्तित्व समयों का न्यूनतम है: तब KK का कोई भी कोशी आँकड़ा ऐसा हल जन्माता है जो अपने आरंभिक समय के आगे कम से कम TT^* तक जीता है। इसे tn>T+T/2t_n > T_+ - T^*/2 के साथ (tn,x(tn))(t_n, x(t_n)) पर लगाने से xx T+T_+ के पार बढ़ जाता है (विस्तार अद्वितीयता से xx के साथ सहमत होता है, फिर उसे आगे बढ़ा देता है): जो महत्तमता का विरोधाभास है। अतः आलेख T+T_+ से पहले KK को निश्चित रूप से छोड़ देता है। 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 विवृत अंतराल) हो और संतत α,β\alpha, \beta के साथ f(t,x)α(t)x+β(t)\norm{f(t, x)} \leq \alpha(t)\norm x + \beta(t), तो प्रत्येक महत्तम हल समूचे 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 हो, तो आलेख T+T_+ के निकट I×RdI\times\R^d के किसी संहत में रहता: जो पलायन (प्रमेय 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 के हल एक dd-विमीय सदिश समष्टि SHS_H बनाते हैं; और प्रत्येक 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 ss पर मान vv वाला हल है,

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) (समाकल के भीतर अवकलन वैध है: समाकल्य tt में C1\mathcal C^1 है और (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) (II पर det\det का प्रसार): इसलिए 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!} अभिसरित होती है (किसी भी उपगुणात्मक मानक में, निरपेक्षतः), और AB=BAAB = BA होने पर eA+B=eAeB\eu^{A+B} = \eu^A\eu^B, तथा 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}। इसके अतिरिक्त: यदि AA के प्रत्येक अभिलक्षणिक मान λ\lambda के लिए Reλ<α<0\operatorname{Re}\lambda < -\alpha < 0 हो, तो t0t \geq 0 पर etACeαt\vertiii{\eu^{tA}} \leq C\eu^{-\alpha t}

उपपत्ति. अभिसरण: An/n!An/n!\vertiii{A^n/n!} \leq \vertiii A^n/n!, जो योग्य है (बानाख बीजगणित MdM_d में अभ्यास 7.1(ख))। क्रमविनिमेय 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.18) से A=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}, जिसमें किसी δ>0\delta > 0 (t0t \geq 0) के लिए etDe(α+δ)t\vertiii{\eu^{tD}} \leq \eu^{-(\alpha + \delta)t} और tt में etN\eu^{tN} बहुपद (शून्यभाविता): अतः गुणनफल Ceαt\leq C\eu^{-\alpha t} है (eδt\eu^{-\delta t} बहुपद को हरा देता है)।

उदाहरण 19.9 (समतल, वर्गीकृत)

व्युत्क्रमणीय AM2(R)A \in M_2(\R) के साथ x=Axx' = Ax के लिए 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; और हल आवर्तकाल 2πω\frac{2\pi}\omega वाले eτt/2×\eu^{\tau t/2}\times घूर्णन हैं।
  • τ=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 पर स्थानीय लिप्शिट्ज़ है। x(0)=x0x(0) = x_0 वाले महत्तम हल के लिए φt(x0)=x(t)\varphi_t(x_0) = x(t) लिखिए (सदिश क्षेत्र का प्रवाह); स्वायत्तता जहाँ परिभाषित हो वहाँ समूह गुण φ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; और स्पर्शोन्मुख रूप से स्थायी यदि इसके अतिरिक्त xˉ\bar x के निकट सभी x0x_0 के लिए φt(x0)xˉ\varphi_t(x_0) \to \bar x

प्रमेय 19.11 (ल्यापुनोव फलन)

मान लीजिए xˉ\bar x कोई साम्यावस्था है और V ⁣:VRV \colon \mathcal V \to \R xˉ\bar x के किसी प्रतिवेश पर ऐसा C1\mathcal C^1 फलन है कि:

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 स्थायी है। और यदि इसके अतिरिक्त xˉ\bar x के बाहर V˙<0\dot V < 0 हो, तो 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 लीजिए (संहतता, धनात्मकता) और B(xˉ,δ)B(\bar x, \delta) पर V<mV < m वाला δ<ε\delta < \varepsilon चुनिए (सांतत्य)। 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 (गोले के भीतर xˉ\bar x से दूरी ρ\geq \rho वाले बिंदुओं के लिए 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 पर स्थानांतरित कीजिए और g(x)=o(x)g(x) = o(\norm x) (C1\mathcal C^1 अवकलनीयता) के साथ F(x)=Ax+g(x)F(x) = Ax + g(x) लिखिए। etACeαt\vertiii{\eu^{tA}} \leq C\eu^{-\alpha t} वाला α>0\alpha > 0 चुनिए (t0t \geq 0; प्रमेय 19.8) और xr\norm x \leq r के लिए g(x)α2Cx\norm{g(x)} \leq \frac{\alpha}{2C}\norm x वाला r>0r > 0b(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 हो, तो यह पूर्व परिबंध सभी tt के लिए x(t)<r\norm{x(t)} < r बनाए रखता है (एक सांतत्य/उन्नयन तर्क: कड़े आकलन को देखते हुए जिन समयों पर 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,x)(x, x')-समतल में लोलक x=sinxx'' = -\sin 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

स्पष्ट रूप से हल कीजिए और महत्तम अंतराल निर्धारित कीजिए: (क) x=x2x' = x^2, x(0)=1x(0) = 1; (ख) x=1+x2x' = 1 + x^2, x(0)=0x(0) = 0; (ग) x=x(1x)x' = x(1-x), x(0)=12x(0) = \frac12। प्रत्येक उत्तर का प्रमेय 19.5(2) और उपप्रमेय 19.6 से मेल कराइए।

हल

हल — अभ्यास 19.1.

(क) चर पृथक करने पर: (,1)\intoo{-\infty}1 पर x(t)=11tx(t) = \frac1{1 - t}: x(t)+x(t) \to +\infty के साथ T+=1T_+ = 1 पर विस्फोट — अर्थात् ठीक प्रमेय 19.5(2)। (ख) (π/2,π/2)\intoo{-\pi/2}{\pi/2} पर x(t)=tantx(t) = \tan t: दोनों सिरों पर विस्फोट। (ग) x(t)=11+etx(t) = \frac{1}{1 + \eu^{-t}}, वैश्विक: हल (0,1)\intoo01 में ही रहता है, जो कोई परिबद्ध समुच्चय है, अतः आलेख परिमित समय में R×R\R\times\R के प्रत्येक संहत से बच नहीं सकता — T±=±T_\pm = \pm\infty। ध्यान दीजिए कि (क), (ख) उपप्रमेय 19.6 का खंडन नहीं करते: x2x^2 और 1+x21 + x^2 की वृद्धि अधिरैखिक है।

अभ्यास 19.2

मान लीजिए x,yx, y x=f(t,x)x' = f(t,x) को हल करते हैं, जहाँ ff R×Rd\R\times\R^d पर xx में वैश्विक रूप से LL-लिप्शिट्ज़ है। (क) 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} प्राप्त हो जाता है। (ख) निष्कर्ष निकालिए कि प्रवाह प्रतिचित्रण x0x(t;x0)x_0 \mapsto x(t; x_0) संतत है, और परिबद्ध समुच्चयों पर वस्तुतः लिप्शिट्ज़।

हल

हल — अभ्यास 19.2.

(क) यह t0=0t_0 = 0 के साथ उपप्रमेय 19.4 का आकलन है। तीक्ष्णता: x=Lxx' = Lx (वैश्विक रूप से LL-लिपशित्ज़) के लिए दो हल ठीक (x0y0)eLt(x_0 - y_0)\eu^{Lt} भर भिन्न होते हैं। (ख) आकलन यह कहता है: समय-tt प्रवाह प्रतिचित्रण प्रारंभिक शर्त में eLt\eu^{L\abs t}-लिपशित्ज़ है — अर्थात् सांतत्य, संहतों में tt के लिए एकसमान रूप से; और अवैश्विक-लिपशित्ज़ ff के परिबद्ध समुच्चयों पर वही तर्क प्रक्षेप-पथों के चारों ओर किसी संहत नली पर स्थानीय अचर के साथ चलाइए।

अभ्यास 19.3 ★★

दर्शाइए कि निम्नलिखित में से प्रत्येक के सभी महत्तम हल R\R पर वैश्विक हैं, और उचित प्रमेय उद्धृत कीजिए: (क) x=sin(tx)x' = \sin(tx); (ख) x=tx1+x2x' = \frac{t\,x}{1 + x^2}; (ग) संतत qq के साथ x+q(t)x=0x'' + q(t)x = 0 (प्रथम-कोटि निकाय में बदलिए); (घ) संतत और परिबद्ध AA के साथ x=A(t)xx' = A(t)x — और x(t)\norm{x(t)} पर ग्रोनवाल परिबंध दीजिए।

हल

हल — अभ्यास 19.3.

(क) sin(tx)1\abs{\sin(tx)} \leq 1: परिबद्ध, अर्थात् α=0\alpha = 0, β=1\beta = 1 के साथ रैखिक वृद्धि: U=R×RU = \R\times\R पर उपप्रमेय 19.6। (ख) tx/(1+x2)t12\abs{tx/(1 + x^2)} \leq \abs t\cdot\frac12: फिर से अवरैखिक (वस्तुतः समय-संहतों पर परिबद्ध): अतः वैश्विक। (ग) 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 की व्यवस्था)। (घ) वैश्विक; और उपप्रमेय 19.6 की भाँति ग्रोनवाल: M=supAM = \sup\vertiii{A} के साथ x(t)x(t0)eMtt0\norm{x(t)} \leq \norm{x(t_0)}\,\eu^{M\abs{t - t_0}}

अभ्यास 19.4 ★★

(क) 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) के लिए etA\eu^{tA} परिकलित कीजिए। (ख) प्रणोदित दोलक x+x=cos(ωt)x'' + x = \cos(\omega t) को (निकाय रूप में) दुआमेल से हल कीजिए, ω1\omega \neq 1 और ω=1\omega = 1 के लिए: अनुनाद धर्मनिरपेक्ष पद tsintt\sin t के रूप में प्रकट होता है।

हल

हल — अभ्यास 19.4.

(क) पहले के लिए 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

(ख) निकाय 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 के लिए: (क) दर्शाइए कि दो हलों का व्रोन्स्कियन w=x1x2x1x2w = x_1x_2' - x_1'x_2 w=pww' = -p\,w (आबेल) को संतुष्ट करता है, और निष्कर्ष निकालिए कि जिन दो हलों के लिए कहीं w0w \neq 0 हो वे एक आधार बनाते हैं। (ख) कोई एक अलुप्त हल x1x_1 दिया हो, तो कोटि अपचयन से सामान्य हल ढूँढ़िए: x2=x1epx12x_2 = x_1\int \frac{\eu^{-\int p}}{x_1^2} रखिए और सत्यापित कीजिए। इसे x1(t)=t2x_1(t) = t^2 के साथ (0,+)\intoo0{+\infty} पर t2x2x=0t^2x'' - 2x = 0 पर लगाइए।

हल

हल — अभ्यास 19.5.

(क) 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) कोई आधार है।

(ख) 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 के लिए, अर्थात् x1=t2x_1 = t^2 के साथ x2t2x=0x'' - \frac2{t^2}x = 0 (p=0p = 0): 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 वाला प्रत्येक हल बढ़ता हुआ और वैश्विक है, जिसकी \mp\infty पर सीमाएँ 00 और 11 हैं; और स्पष्ट रूप से हल करके पुष्टि कीजिए। अधिक व्यापक रूप से दर्शाइए कि अदिश स्वायत्त हल एकदिष्ट होते हैं, और निष्कर्ष निकालिए: विमा 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: अनंतस्पर्शी रूप से स्थायी (विमा 11 में प्रमेय 19.12)। 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; और सममित रूप से -\infty पर x0x \to 0। स्पष्ट रूप से 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 पर विचार कीजिए। (क) दर्शाइए कि HH हलों के अनुदिश अचर है। (ख) 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.

(क)  ⁣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। (ख) H=y22+x44H = \frac{y^2}2 + \frac{x^4}4 के स्तर समुच्चय संहत हैं (HH प्रबलकारी), अतः हल संहतों में फँस जाते हैं: वैश्विक और परिबद्ध (प्रमेय 19.5)। (0,0)(0,0) की स्थायित्व: V=HV = H V˙=0\dot V = 0 के साथ धनात्मक निश्चित है (H=0H = 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। (क) दर्शाइए कि V(x,y)=y22+1cosxV(x, y) = \frac{y^2}2 + 1 - \cos x V˙=cy20\dot V = -cy^2 \leq 0 को संतुष्ट करता है: अतः मूल बिंदु स्थायी है। (ख) V˙\dot V पूरे अक्ष y=0y = 0 पर लुप्त हो जाता है: अतः ल्यापुनोव की कड़ी कसौटी विफल हो जाती है। फिर भी रैखिकीकरण (प्रमेय 19.12) से स्पर्शोन्मुख स्थायित्व सिद्ध कीजिए: (0,0)(0,0) पर रैखिकीकृत आव्यूह के अभिलक्षणिक मान परिकलित कीजिए और प्रत्येक c>0c > 0 के लिए Re<0\operatorname{Re} < 0 जाँचिए। (ग) साम्यावस्था (π,0)(\pi, 0) पर क्या होता है? रैखिकीकरण परिकलित कीजिए और निष्कर्ष निकालिए (एक अभिलक्षणिक मान धनात्मक: अस्थायित्व — आप अस्थायित्व के कथन का अनौपचारिक उपयोग कर सकते हैं अथवा रैखिक निकाय का कोई स्पष्ट भागता हुआ हल प्रस्तुत कर सकते हैं)।

हल

हल — अभ्यास 19.8.

(क) 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)। (ख) (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}}2c2c \geq 2 होने पर दोनों वास्तविक ऋणात्मक, और 0<c<20 < c < 2 होने पर वास्तविक भाग c2<0-\frac c2 < 0 वाले सम्मिश्र। सभी दशाओं में Reλ<0\operatorname{Re}\lambda < 0: अतः प्रमेय 19.12 अनंतस्पर्शी स्थायित्व दे देती है (अपह्रासी V˙\dot V के बावजूद)। (ग) (π,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)। अस्थिर अभिलक्षणिक सदिश के अनुदिश रैखिक निकाय का λ+>0\lambda_+ > 0 वाला स्पष्ट रूप से भागता हुआ हल eλ+tv+\eu^{\lambda_+t}v_+ है: अर्थात् उलटा लोलक प्रत्येक अवमंदन के लिए अस्थिर है।

अभ्यास 19.9 ★★

(अद्वितीयता की सीमा) α(0,1)\alpha \in \intoo01 के लिए दर्शाइए कि समस्या x=xαx' = \abs x^{\alpha}, x(0)=0x(0) = 0 के अनंत हल हैं (समस्या 7.1, भाग III को अनुकूलित कीजिए)। इसके विपरीत दर्शाइए कि α=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 पर स्थानीय लिपशित्ज़ शर्त है: α<1\alpha < 1 के लिए xα\abs x^\alpha के वहाँ अपरिबद्ध अंतर भागफल हैं।

अभ्यास 19.10 ★★★

(पूर्व परिबंध हलों को फँसा लेते हैं) मान लीजिए F ⁣:RdRdF \colon \R^d \to \R^d स्थानीय लिप्शिट्ज़ है और xR\norm x \geq R होने पर F(x),x0\langle F(x), x\rangle \leq 0। (क) दर्शाइए कि संवृत गोला Bˉ(0,R)\bar B(0, R) धनात्मक रूप से अपरिवर्ती है: भीतर से आरंभ होने वाले हल t0t \geq 0 के लिए भीतर ही रहते हैं। (यदि x(t2)>R\norm{x(t_2)} > R हो, तो x(t1)=R\norm{x(t_1)} = R वाला अंतिम समय t1<t2t_1 < t_2 लीजिए और [t1,t2]\intcc{t_1}{t_2} पर  ⁣d ⁣dtx(t)2\frac{\dd}{\dd t}\norm{x(t)}^2 का अध्ययन कीजिए।) (ख) गोले के आँकड़ों के लिए वैश्विक अग्रगामी अस्तित्व निष्कर्ष के रूप में निकालिए। फिर प्रवणक निकाय x=G(x)x' = -\nabla G(x), GC2G \in \mathcal C^2 लीजिए, जहाँ x\norm x \to \infty पर G(x)+G(x) \to +\infty: दर्शाइए कि हलों के अनुदिश GG घटता है, कि प्रत्येक हल (परिबद्ध) अधोस्तर समुच्चय {GG(x0)}\{G \leq G(x_0)\} में रहता है, और वैश्विक अग्रगामी अस्तित्व निष्कर्ष निकालिए।

हल

हल — अभ्यास 19.10.

(क) मान लीजिए x(0)R\norm{x(0)} \leq R वाले किसी t2>0t_2 > 0 के लिए x(t2)>R\norm{x(t_2)} > 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 और (t1,t2]\intoc{t_1}{t_2} पर x(t)>R\norm{x(t)} > R। उस अंतराल पर 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: विरोधाभास। गोला धनात्मक रूप से निश्चर है। (ख) संहत गोले में फँसा कोई हल 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 पर विचार कीजिए। (क) दर्शाइए कि महत्तम हल किसी [0,T+)\intco0{T_+} पर विद्यमान है जिसमें T+<T_+ < \infty: y=y2y' = y^2, y(0)=1y(0) = 1 से तुलना कीजिए (अभीष्ट तुलना प्रमेयिका सिद्ध कीजिए: यदि xF(x)x' \geq F(x) और x(0)y(0)x(0) \geq y(0) के साथ y=F(y)y' = F(y) हो, तो जहाँ दोनों जीते हैं वहाँ xyx \geq y), और T+1T_+ \leq 1 निष्कर्ष निकालिए। (ख) T+T_+ को नीचे से परिबद्ध कीजिए: [0,1]\intcc01 पर xx2+1x' \leq x^2 + 1; z(0)=1z(0) = 1 वाले अधिहल z=z2+1z' = z^2 + 1 से तुलना कीजिए, जो z(t)=tan(t+π4)z(t) = \tan\bigl(t + \frac\pi4\bigr) से हल हो जाता है, और T+π4T_+ \geq \frac\pi4 निष्कर्ष निकालिए। (ग) π4T+1\frac\pi4 \leq T_+ \leq 1 संकलित कीजिए और सार बाँधिए: दायाँ पक्ष अधिरैखिक रूप से बढ़ता हो तो वैश्विक अस्तित्व मर जाता है (अभ्यास 19.3 इसका प्रतिबिंदु है), और सीमा  ⁣dsF(s)\int^{\infty}\frac{\dd s}{F(s)} के अभिसरण पर पड़ती है।

हल

हल — अभ्यास 19.11.

(क) तुलना प्रमेयिका: मान लीजिए उभयनिष्ठ अंतराल पर 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 0F(x)=x2F(x) = x^2 के साथ: y(t)=11ty(t) = \frac1{1 - t} 11 पर विस्फोटित होता है, और जब तक xx जीवित है तब तक xyx \geq y; यदि T+>1T_+ > 1 हो, तो xx t=1t = 1 पर परिमित होता जबकि yy \to \infty पर हावी रहता: असंगत। अतः T+1T_+ \leq 1

(ख) उलटी तुलना (वही प्रमेयिका, भूमिकाएँ बदली हुई): [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

(ग) मिलाकर: π4T+1\frac\pi4 \leq T_+ \leq 1 (संख्यात्मक रूप से T+0.96T_+ \approx 0.96)। शिक्षा: xx में अधिरैखिक FF वाले x=F(t,x)x' = F(t, x) के लिए हल परिमित समय में तब-तब फट जाते हैं जब  ⁣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 को हल करते हैं। (क) व्रोन्स्कियन सर्वसमिका स्थापित कीजिए: W=uvuvW = uv' - u'v के साथ W=(q1q2)uvW' = (q_1 - q_2)\,uv। (ख) (श्टुर्म) दर्शाइए कि uu के दो क्रमागत शून्यों a<ba < b के बीच या तो vv (a,b)\intoo ab में कहीं लुप्त होता है, या वहाँ q1=q2q_1 = q_2 और vuv \propto u ((a,b)\intoo ab पर u>0u > 0 और v>0v > 0 भी मान लीजिए; (क) का aa से bb तक समाकलन कीजिए और सीमांत पदों W(a),W(b)W(a), W(b) के चिह्न देखिए)। (ग) निष्कर्ष निकालिए: qm>0q \geq m > 0 वाले u+qu=0u'' + qu = 0 के हल लंबाई π/m\pi/\sqrt m के प्रत्येक अंतराल में कम से कम एक बार लुप्त होते हैं (v+mv=0v'' + mv = 0 से तुलना कीजिए); और q0q \leq 0 वाले हल R\R पर अधिकतम एक बार। दोनों को q±1q \equiv \pm1 पर परखिए।

हल

हल — अभ्यास 19.12.

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

(ख) मान लीजिए a<ba < b uu के क्रमागत शून्य हैं; (a,b)\intoo ab पर u>0u > 0 मानकीकृत कीजिए (अतः u(a)>0u'(a) > 0, u(b)<0u'(b) < 0 — जो अद्वितीयता से शून्येतर है, क्योंकि u(a)=u(a)=0u(a) = u'(a) = 0 u0u \equiv 0 अनिवार्य कर देता)। मान लीजिए (a,b)\intoo ab में vv का कोई शून्य नहीं; वहाँ v>0v > 0 मानकीकृत कीजिए (अतः सांतत्य से v(a),v(b)0v(a), v(b) \geq 0)। (क) का समाकलन कीजिए:

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। पूरी शृंखला में बराबरी: विवृत अंतराल पर uv>0uv > 0 के साथ (q1q2)uv=0\int(q_1 - q_2)uv = 0 वहाँ q1=q2q_1 = q_2 अनिवार्य कर देता है; और W(a)=W(b)=0W(a) = W(b) = 0 v(a)=v(b)=0v(a) = v(b) = 0 अनिवार्य कर देता है; तब [a,b]\intcc ab पर W0W \equiv 0 (उसका अवकलज लुप्त होता है), अर्थात् (a,b)\intoo ab पर (v/u)=W/u2=0(v/u)' = -W/u^2 = 0: vuv \propto u

(ग) 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: तब (ख) से v+qv=0v'' + qv = 0 का प्रत्येक हल vv लंबाई π/m\pi/\sqrt m के प्रत्येक विवृत अंतराल में लुप्त होता है (अपह्रासी विकल्प qmq \equiv m में vuv \propto u भी लुप्त होता है)। और यदि q0q \leq 0 हो: तो q1=qq_1 = q, u=vu = v के साथ (ख) लगाइए, तथा v=0v'' = 0 के शून्य-रहित हल 1\mathbf 1 के साथ q2=0q_2 = 0। यदि vv के दो क्रमागत शून्य होते, तो (ख) या तो उनके बीच 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 — निकाय के रूप में: R2\R^2 पर x=yx' = y, y=sinxy' = -\sin x — गतिकी की ड्रोसोफिला है: लिखने में सरल, प्रारंभिक सूत्रों से हल करने में असंभव, फिर भी गुणात्मक विधि से पूर्णतः बोधगम्य। मान लीजिए E(x,y)=y22cosxE(x, y) = \frac{y^2}2 - \cos x (ऊर्जा)।

भाग I — वैश्विक संरचना।

  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 के साथ ल्यापुनोव; अस्पर्शोन्मुखता: ऊर्जा संरक्षण कक्षाओं को मूल बिंदु से दूर स्तर वक्रों पर फँसा देता है।)

भाग II — दोलन और उनका आवर्तकाल। 1<E0<1-1 < E_0 < 1 स्थिर कीजिए और a(0,π)a \in \intoo0\pi (आयाम) के साथ E0=cosaE_0 = -\cos a लिखिए।

  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

    दर्शाइए (एक पूर्ण दीर्घवृत्तीय समाकल)। प्रभावी अभिसरण से निष्कर्ष निकालिए कि a0+a \to 0^+ पर T(a)2πT(a) \to 2\pi: अर्थात् संनादी सीमा, जो आयाम से स्वतंत्र है — गैलीलियो की सन्निकट समकालिकता, और उसका यथार्थ संशोधन T(a)=2π(1+k24+O(k4))T(a) = 2\pi\bigl(1 + \frac{k^2}4 + O(k^4)\bigr) (समाकल्य का प्रसार कीजिए और प्रसामान्य अभिसरण से उचित ठहराते हुए पदशः समाकलन कीजिए)।

  4. दर्शाइए कि aπa \to \pi^- पर T(a)+T(a) \to +\infty (k1k \to 1 होने पर φ=π2\varphi = \frac\pi2 के निकट समाकल्य को नीचे से परिबद्ध कीजिए, अथवा एकदिष्ट अभिसरण लगाइए): अर्थात् विभाजक के पास पहुँचते हुए लोलक बिना किसी सीमा के धीमा होता जाता है।

भाग III — विभाजक।

  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 के साथ)। सीधे सत्यापित कीजिए कि वह लोलक समीकरण हल करता है, और t±t \to \pm\infty पर उसकी तथा y(t)y(t) की सीमाएँ परिकलित कीजिए।

  2. निष्कर्ष निकालिए: विभाजक कक्षा काठी (π,0)(-\pi, 0) (tt \to -\infty पर) को काठी (π,0)(\pi, 0) (t+t \to +\infty पर) से जोड़ती है, पर किसी तक भी परिमित समय में नहीं पहुँचती — जो अद्वितीयता के अनुरूप है (किसी काठी तक परिमित समय में पहुँचना उपप्रमेय 19.4 का खंडन क्यों करता?)।

भाग IV — घूर्णन, और पूरा चित्र।

  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. पूरा प्रावस्था चित्र (अध्याय का चित्र) संकलित कीजिए और उसकी प्रत्येक विशेषता को पूरी तरह उचित ठहराइए, तथा विधि का दस-पंक्ति सार लिखिए: ऊर्जा, स्तर समुच्चय, संहतता, अद्वितीयता — अध्याय की प्रत्येक प्रमेय कहाँ आई? हमें कहीं भी सामान्य हल के किसी सूत्र की आवश्यकता पड़ी?

भाग V — सूक्ष्मदर्शी के नीचे आवर्तकाल फलन।

  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. प्रश्न 10 के घूर्णन आवर्तकाल τ\tau पर लौटिए: दर्शाइए कि τ\tau (1,+)\intoo1{+\infty} पर कड़ाई से घटता हुआ है, कि E01+E_0 \to 1^+ पर τ(E0)+\tau(E_0) \to +\infty (एकदिष्ट अभिसरण), और यह कि E0+E_0 \to +\infty पर 2E0τ(E0)2π\sqrt{2E_0}\,\tau(E_0) \to 2\pi (प्रभावी अभिसरण): अर्थात् तेज चक्करदार घूर्णन स्पर्शोन्मुख रूप से कोणीय चाल 2E0\sqrt{2E_0} का मुक्त घूर्णन है।

भाग VI — विधि का निर्यात: लोत्का–वोल्तेरा। लोलक का नुस्खा — प्रथम समाकल, संहत स्तर वक्र, अद्वितीयता — किसी पारितंत्र को भी हल कर देता है। α,β,γ,δ>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 साम्यावस्था के चारों ओर एक संवृत वक्र है: किसी संहत अंतराल [x,x+]x\intcc{x_-}{x_+} \ni x_* पर दो संतत शाखाएँ y±(x)y_\pm(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 और γ\gamma के स्थान पर αε\alpha - \varepsilon और γ+ε\gamma + \varepsilon आ जाते हैं। औसत जनसंख्याओं का क्या होता है? दांकोना के प्रेक्षण (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)} है, और यह कि E0E_0 \to \infty पर मॉडुलन अनुपात E0+1E011\sqrt{\frac{E_0 + 1}{E_0 - 1}} \to 1: अर्थात् तेज घूर्णन स्पर्शोन्मुख रूप से एकसमान होता है।
  9. (घूर्णन आवर्तकाल की एकदिष्टता) दर्शाइए कि τ(E0)\tau(E_0) C1\mathcal C^1 है और (1,+)\intoo1{+\infty} पर कड़ाई से घटता हुआ (समाकल चिह्न के भीतर अवकलन कीजिए, प्रत्येक [1+δ,)\intco{1 + \delta}\infty पर प्रभुत्व के साथ), जिसमें E01+E_0 \to 1^+ पर τ\tau \to \infty और E0E_0 \to \infty पर τ0\tau \to 0। ऊर्जा अक्ष के अनुदिश लोलक का पूरा द्विशाखन चित्र संकलित कीजिए: E0=1E_0 = -1 पर साम्यावस्थाएँ, 1<E0<1-1 < E_0 < 1 पर 2π2\pi से \infty तक बढ़ते आवर्तकाल वाले दोलन, 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); रैखिकीकरण cos(kπ)=1-\cos(k\pi) = \mp1 के साथ (0110)\bigl(\begin{smallmatrix}0&1\\ \mp1&0\end{smallmatrix}\bigr): सम kk के लिए अभिलक्षणिक मान ±i\pm\iu (केंद्र प्रकार, अपने आप में अनिर्णायक) और विषम kk के लिए ±1\pm1 (काठी)।

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 लिखने पर वह समुच्चय x[2kπa,2kπ+a]x \in [2k\pi - a, 2k\pi + a] पर संवृत वक्रों y=±2(cosxcosa)y = \pm\sqrt{2(\cos x - \cos 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 पर बना रहता है, जिसकी मूल बिंदु से दूरी धनात्मक है (वक्र xx-अक्ष से केवल ±a\pm a पर मिलता है): φt(a,0)↛(0,0)\varphi_t(a, 0) \not\to (0,0)

4. स्तर वक्र CaC_a पर: कोई साम्यावस्था नहीं (y=0y = 0 0<a<π0 < a < \pi के लिए sin(±a)0\sin(\pm a) \neq 0 के साथ x=±ax = \pm a अनिवार्य कर देता है), अतः चाल (y,sinx)\norm{(y, -\sin x)} का संहत CaC_a पर कोई धनात्मक निम्निष्ठ mm है। (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 में वापस पार होता है: अतः हल समय TT पर (a,0)(a, 0) पर लौट आता है। तब अद्वितीयता (उपप्रमेय 19.4) प्रसारित हो जाती है: सभी tt के लिए x(t+T)=x(t)x(t + T) = x(t): अर्थात् आवर्ती।

5. उस शाखा पर जहाँ y>0y > 0: y=2(cosxcosa)y = \sqrt{2(\cos x - \cos a)} और  ⁣dt= ⁣dxy\dd t = \frac{\dd x}{y}; a-a से aa तक xx के समाकलन से आधा आवर्त मिलता है, और सममिति 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^- पर अभिसरण: sina>0\sin a > 0 के साथ cosxcosa=sin(a)(ax)+O((ax)2)\cos x - \cos a = \sin(a)(a - x) + O((a-x)^2): अतः समाकल्य (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 दे देती है। u=k2sin2φu = k^2\sin^2\varphi (k<1k < 1 के लिए प्रसामान्य अभिसरण) और 0π/2sin2=π4\int_0^{\pi/2}\sin^2 = \frac\pi4 के साथ (1u)1/2=1+u2+3u28+(1 - u)^{-1/2} = 1 + \frac u2 + \frac{3u^2}8 + \cdots का प्रसार करने पर:

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} तक जाते हैं, जिसका समाकल अपसारित होता है: अतः एकदिष्ट अभिसरण से aπa \to \pi^- पर T(a)+T(a) \to +\infty

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: अर्थात् प्रक्षेप-पथ विभाजक पर है, और जहाँ y>0y > 0 हो वहाँ y2=2(1+cosx)y^2 = 2(1 + \cos x) का अवकलन y=sinxy' = -\sin x पुनरुत्पन्न कर देता है। सीमाएँ: t±t \to \pm\infty पर x±πx \to \pm\pi और y0y \to 0

9. कक्षा अग्र दिशा में काठी (π,0)(\pi, 0) की ओर और पश्च दिशा में (π,0)(-\pi, 0) की ओर जाती है पर कभी पहुँचती नहीं: यदि वह किसी परिमित समय tt^* पर (π,0)(\pi, 0) तक पहुँच जाती, तो दो भिन्न महत्तम हल — विभाजक हल और काठी पर अचर हल — एक ही बिंदु (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-आवर्ती है, अतः जब-जब xx 2π2\pi भर आगे बढ़ता है तब-तब yy अपने मान पर लौट आता है; और इसके लिए आवश्यक समय

τ(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) विभाजित कीजिए) के साथ आगमन दिखाया गया मान दे देता है। द्विपद श्रेणी: cn=(2n)!4n(n!)2(0,1]c_n = \frac{(2n)!}{4^n(n!)^2} \in \intoc01 के साथ (1u)1/2=n0cnun(1 - u)^{-1/2} = \sum_{n\geq0}c_nu^n, त्रिज्या 11kk0<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 के अनुदिश प्रभावी अभिसरण लागू होता है। a0+a \to 0^+ पर T2πT \to 2\pi (प्रश्न 6) और aπa \to \pi^- पर T+T \to +\infty (प्रश्न 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_* पर लुप्त होता है, और 0+0^+ पर तथा ++\infty पर दोनों जगह f+f \to +\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 में संवृत है क्योंकि QQ की परिसीमा पर H+H \to +\infty। कोई महत्तम हल अपने संहत स्तर समुच्चय पर बना रहता है, अतः वह परिमित समय में प्रत्येक संहत को नहीं छोड़ सकता: इसलिए प्रमेय 19.5 उसे वैश्विक बना देती है।

19. h>hh > h_* स्थिर कीजिए और c=hg(y)>f(x)c = h - g(y_*) > f(x_*) रखिए। चूँकि ff (0,x]\intoc0{x_*} पर ++\infty से f(x)f(x_*) तक कड़ाई से घटता है और [x,+)\intco{x_*}{+\infty} पर कड़ाई से बढ़कर ++\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 में संतत हैं (yy_* के दोनों ओर gg के कड़ाई से एकदिष्ट संतत प्रतिबंधों के प्रतिलोम), जिनमें xx±x \to x_\pm पर y±(x)yy_\pm(x) \to y_*; और x=x±x = x_\pm पर अद्वितीय हल y=yy = y_* है। अतः {H=h}Q\{H = h\}\cap Q [x,x+]\intcc{x_-}{x_+} पर y+y_+ और yy_- के आलेखों का संघ है, जो (x±,y)(x_\pm, y_*) पर चिपके हुए हैं: अर्थात् (x,y)(x_*, y_*) के चारों ओर कोई संवृत वक्र — लोलक के अंडाकारों का सादृश्य।

20. h>hh > h_* के साथ Ch={H=h}QC_h = \{H = h\}\cap Q लीजिए: 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_* के आगे वर्धमान और धनात्मक है), जबकि स्तर संबंध और टेलर असमिका ChC_h के संहत yy-परास पर g(y(x))g(y)12(maxg)(yy(x))2g(y_-(x)) - g(y_*) \leq \frac12\,\bigl(\max g''\bigr)\,(y_* - y_-(x))^2 दे देते हैं: अतः c>0c > 0 के साथ yy(x)cx+xy_* - y_-(x) \geq c\,\sqrt{x_+ - x}, और समाकल्य 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_* में प्रवेश करती है, सममित आकलन से (ff और gg की भूमिकाएँ बदलकर) ऊपरी कोने (x,y+)(x_*, y_+) तक चढ़ती है, और इसी प्रकार चारों चापों के चारों ओर: अतः परिमित समय T>0T > 0 के बाद हल अपने आरंभ बिंदु पर लौट आता है। उपप्रमेय 19.4 से वह TT-आवर्ती है — प्रश्न 4 का तर्क, अक्षरशः।

21. QQ में किसी TT-आवर्ती कक्षा पर 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: अतः cosx=±1\cos x = \pm1 पर y2y^2 के चरम, जो कहे गए परिबंध देते हैं और 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}): अतः विभाजक के पास आवर्त फट जाता है, जो दोलन की ओर से भाग II के T(a)T(a) \to \infty से मेल खाता है। ऊर्जा अक्ष यों पढ़ी जाती है: E0=1E_0 = -1 पर विश्राम; (1,1)\intoo{-1}1 पर दोलन, 2π2\pi \nearrow \infty; E0=1E_0 = 1 पर अनंत रूप से धीमा विभाजक; और उसके परे घूर्णन, 0\infty \searrow 0। एक ही समाकल, और लोलक का समूचा जीवन।