Matemática · किताब 3 · Graduação — Ano 1

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

विश्वविद्यालय गणित — स्नातक वर्ष 1 · Graduação — Ano 1

5रैखिक अवकल समीकरण

अवकल समीकरण पहले उच्चतर माध्यमिक खंड में मिले थे; यहाँ उनका उपचार पूरी उपपत्तियों के साथ और अधिक व्यापकता में है: चर गुणांकों वाले प्रथम कोटि के रैखिक समीकरण (जिन्हें अचरों के विचरण की विधि पूरी तरह हल कर देती है) और अचर गुणांकों वाले द्वितीय कोटि के रैखिक समीकरण, जो दोलनों के आदर्श हैं। दोनों स्थितियाँ एक ही संरचना दिखाती हैं: व्यापक हल == एक विशिष्ट हल ++ समांगी समीकरण का व्यापक हल

5.1 प्रथम कोटि के रैखिक समीकरण

परिभाषा 5.1

मान लीजिए II एक अंतराल है और a,b ⁣:IRa, b \colon I \to \R (या C\C) संतत है। समीकरण

(E) ⁣:y+a(x)y=b(x),(E)\colon\quad y' + a(x)\,y = b(x),

अज्ञात अवकलनीय फलन y ⁣:IRy \colon I \to \R (या C\C) में प्रथम कोटि का रैखिक अवकल समीकरण कहलाता है। समीकरण (H) ⁣:y+a(x)y=0(H)\colon y' + a(x) y = 0 उसका समांगी समीकरण है।

प्रमेय 5.2 (समांगी समीकरण हल करना)

मान लीजिए AA II पर aa का प्रतिअवकलज है (वह विद्यमान है: अध्याय 15)। II पर (H)(H) के हल ठीक ये फलन हैं

y(x)=λeA(x),λR (या C).y(x) = \lambda\, \eu^{-A(x)}, \qquad \lambda \in \R \text{ (या } \C).

उपपत्ति. ये फलन हल हैं: y=λAeA=ayy' = -\lambda A' \eu^{-A} = -a y। विलोमतः, मान लीजिए yy (H)(H) को हल करता है और z(x)=y(x)eA(x)z(x) = y(x)\, \eu^{A(x)} रखिए। तब

z=yeA+yaeA=(y+ay)eA=0,z' = y' \eu^{A} + y\, a\, \eu^{A} = (y' + ay)\, \eu^{A} = 0 ,

अतः zz अंतराल II पर अचर है, मान लीजिए z=λz = \lambda: y=λeAy = \lambda \eu^{-A}। (तर्क पर ध्यान दीजिए: कोई हल नहीं छूटता, क्योंकि प्रत्येक हल घोषित रूप में लिखा जा चुका है।)

उदाहरण 5.3 (चर गुणांक वाला समांगी समीकरण)

R\R पर y+(cosx)y=0y' + (\cos x)\,y = 0 हल कीजिए। a(x)=cosxa(x) = \cos x का एक प्रतिअवकलज A(x)=sinxA(x) = \sin x है, अतः हल हैं

y(x)=λesinx,λR.y(x) = \lambda\,\eu^{-\sin x}, \qquad \lambda \in \R .

दो पाठ। प्रत्येक हल आवर्ती है (आवर्तकाल 2π2\pi) और λ=0\lambda = 0 के अतिरिक्त कभी लुप्त नहीं होता — λ\lambda का चिह्न सदा के लिए yy का चिह्न रहता है, क्योंकि चरघातांकी शून्य को पार नहीं कर सकता। और y(0)=y0y(0) = y_0 से होकर जाने वाला हल y0esinxy_0\eu^{-\sin x} है: प्रत्येक आरंभिक बिंदु से होकर इस कुल का ठीक एक वक्र, जो प्रमेय 5.4 (2) का एकविमीय चित्र है।

प्रमेय 5.4 (अचरों का विचरण; कोशी समस्या)

ऊपर के संकेतन के साथ:

  1. II पर (E)(E) के हल ठीक ये हैं

    y(x)=(λ+x0xb(t)eA(t) ⁣dt)eA(x),λR,y(x) = \Bigl(\lambda + \int_{x_0}^{x} b(t)\, \eu^{A(t)} \dd t \Bigr)\, \eu^{-A(x)}, \qquad \lambda \in \R,

    जहाँ x0Ix_0 \in I नियत है। समतुल्य रूप से: (H)(H) का व्यापक हल और (E)(E) का एक विशिष्ट हल।

  2. प्रत्येक x0Ix_0 \in I और y0y_0 के लिए कोशी समस्या(E)(E) और y(x0)=y0y(x_0) = y_0” का II पर ठीक एक हल है।

उपपत्ति. (1) अचरों का विचरण नामक विधि के अनुसार y=μ(x)eA(x)y = \mu(x)\, \eu^{-A(x)} रूप के हल खोजिए, जहाँ μ\mu अवकलनीय है — इससे व्यापकता नहीं खोती, क्योंकि II पर प्रत्येक फलन इस रूप में लिखा जा सकता है (μ=yeA\mu = y\,\eu^{A})। प्रतिस्थापित करने पर,

y+ay=μeAμaeA+aμeA=μeA,y' + ay = \mu' \eu^{-A} - \mu a \eu^{-A} + a\mu\eu^{-A} = \mu'\, \eu^{-A},

अतः yy (E)(E) को हल करता है यदि और केवल यदि μ(x)=b(x)eA(x)\mu'(x) = b(x)\,\eu^{A(x)}, अर्थात् यदि और केवल यदि किसी अचर λ\lambda के लिए μ(x)=λ+x0xb(t)eA(t) ⁣dt\mu(x) = \lambda + \int_{x_0}^x b(t)\eu^{A(t)}\dd t (अंतराल पर एक ही संतत फलन के दो प्रतिअवकलज अचर से भिन्न होते हैं)।

(2) सूत्र में y(x0)=λeA(x0)y(x_0) = \lambda\,\eu^{-A(x_0)}: शर्त y(x0)=y0y(x_0) = y_0 λ=y0eA(x0)\lambda = y_0 \eu^{A(x_0)} को अद्वितीय रूप से निर्धारित कर देती है।

उदाहरण 5.5

I=(0,+)I = \intoo{0}{+\infty} पर y+yx=x2y' + \dfrac{y}{x} = x^2 हल कीजिए। यहाँ a(x)=1xa(x) = \frac 1x, A(x)=lnxA(x) = \ln x, eA(x)=1x\eu^{-A(x)} = \frac 1x। समांगी हल: λx\frac{\lambda}{x}अचरों का विचरण: μ(x)=x2x=x3\mu'(x) = x^2 \cdot x = x^3, अतः μ=x44+λ\mu = \frac{x^4}{4} + \lambda, और

y(x)=x34+λx,λR.y(x) = \frac{x^3}{4} + \frac{\lambda}{x}, \qquad \lambda \in \R .

आरंभिक शर्त y(1)=0y(1) = 0 के साथ: λ=14\lambda = -\frac14जाँच: y+yx=3x24λx2+x24+λx2=x2y' + \frac yx = \frac{3x^2}{4} - \frac{\lambda}{x^2} + \frac{x^2}{4} + \frac{\lambda}{x^2} = x^2

उदाहरण 5.6 (अनुमान लगाना समाकलन से अच्छा है)

R\R पर y+2xy=xy' + 2x\,y = x हल कीजिए। अचरों का विचरण काम करता है (A=x2A = x^2, μ=xex2\mu' = x\,\eu^{x^2}, μ=12ex2+λ\mu = \frac12\eu^{x^2} + \lambda), पर यह देख लेना कि अचर yp=12y_p = \frac12 समीकरण को हल करता है (0+2x12=x0 + 2x\cdot\frac12 = x), अधिक तेज़ है। समांगी हलों λex2\lambda\,\eu^{-x^2} के साथ:

y(x)=12+λex2,λR.y(x) = \frac12 + \lambda\,\eu^{-x^2}, \qquad \lambda \in \R .

जब x±x \to \pm\infty, तब प्रत्येक हल अत्यंत तेज़ी से 12\frac12 की ओर अभिसरित होता है: अचर विशिष्ट हल एक साम्यावस्था है जिससे सभी हल आ मिलते हैं। सार: सामान्य विधि चलाने से पहले दस सेकंड कोई स्पष्ट विशिष्ट हल (अचर, एकपदी, दाएँ पक्ष का गुणज) ढूँढ़ने में लगाइए; संरचना प्रमेय फिर काम पूरा कर देती है।

टिप्पणी 5.7 (अंतराल महत्त्व रखते हैं)

यह प्रमेय ऐसे अंतराल पर रहती है जहाँ aa और bb संतत हों। R\R^* पर y+yx=0y' + \frac yx = 0 के लिए हल (0,+)\intoo{0}{+\infty} पर λx\frac{\lambda}{x} और (,0)\intoo{-\infty}{0} पर μx\frac{\mu}{x} हैं, जिनके अचर स्वतंत्र हैं: 00 की विचित्रता को पार करके एक ही सूत्र के जुड़ने का कोई कारण नहीं है।

उदाहरण 5.8 (सम्मिश्र दायाँ पक्ष, दो वास्तविक उत्तर)

yy=cosxy' - y = \cos x और yy=sinxy' - y = \sin x को एक ही चाल में हल कीजिए। दाएँ पक्ष eix\eu^{\iu x} के साथ C\C में काम कीजिए: yp=ceixy_p = c\,\eu^{\iu x} आज़माने पर c(i1)eix=eixc(\iu - 1)\eu^{\iu x} = \eu^{\iu x} मिलता है, अतः

c=1i1=1i2,yp=(1+i)(cosx+isinx)2=sinxcosx2+isinxcosx2.c = \frac1{\iu - 1} = \frac{-1 - \iu}2, \qquad y_p = -\frac{(1 + \iu)(\cos x + \iu\sin x)}2 = \frac{\sin x - \cos x}2 + \iu\,\frac{-\sin x - \cos x}2 .

चूँकि समीकरण के गुणांक वास्तविक हैं, वास्तविक और काल्पनिक भाग अलग हो जाते हैं: sinxcosx2\frac{\sin x - \cos x}2 yy=cosxy' - y = \cos x को हल करता है, और sinx+cosx2-\frac{\sin x + \cos x}2 yy=sinxy' - y = \sin x को (पहले की जाँच: अवकलज cosx+sinx2\frac{\cos x + \sin x}2, में से फलन घटाने पर cosx\cos x मिलता है)। एक सम्मिश्र पंक्ति ने अचरों के विचरण की दो दौड़ों की जगह ले ली — वही मितव्ययिता जिसे विधि 5.13 द्वितीय कोटि के लिए व्यवस्थित करता है, और अध्याय 3 का बार-बार मिलने वाला लाभांश।

5.2 अचर गुणांकों वाले द्वितीय कोटि के रैखिक समीकरण

परिभाषा 5.9

मान लीजिए a,bRa, b \in \R और f ⁣:IRf \colon I \to \R संतत है। समीकरण

(E) ⁣:y+ay+by=f(x)(E)\colon\quad y'' + a\,y' + b\,y = f(x)

अचर गुणांकों वाला द्वितीय कोटि का रैखिक समीकरण कहलाता है; (H) ⁣:y+ay+by=0(H)\colon y'' + ay' + by = 0 उसका समांगी समीकरण है, और χ(r)=r2+ar+b\chi(r) = r^2 + ar + b उसका अभिलक्षणिक बहुपद

प्रमेय 5.10 (समांगी हल)

मान लीजिए Δ=a24b\Delta = a^2 - 4b χ\chi का विविक्तकर है। R\R पर (H)(H) के वास्तविक हल ये हैं:

  1. यदि Δ>0\Delta > 0, तो दो वास्तविक मूलों r1r2r_1 \neq r_2 के साथ:   y=λer1x+μer2x\;y = \lambda\, \eu^{r_1 x} + \mu\, \eu^{r_2 x};
  2. यदि Δ=0\Delta = 0, तो द्विक मूल r0r_0 के साथ:   y=(λ+μx)er0x\;y = (\lambda + \mu x)\, \eu^{r_0 x};
  3. यदि Δ<0\Delta < 0, तो मूलों α±iω\alpha \pm \iu\omega (ω>0\omega > 0) के साथ:   y=eαx(λcosωx+μsinωx)\;y = \eu^{\alpha x} (\lambda \cos\omega x + \mu \sin\omega x);

और हर स्थिति में (λ,μ)(\lambda, \mu) R2\R^2 पर चलते हैं।

उपपत्ति. पहले यह देखिए कि rCr \in \C के लिए xerxx \mapsto \eu^{rx} (H)(H) को हल करता है यदि और केवल यदि χ(r)=0\chi(r) = 0 (प्रतिस्थापित कीजिए: (r2+ar+b)erx=0(r^2 + ar + b)\eu^{rx} = 0)। इसीलिए चरघातांकी स्वाभाविक पहला अनुमान हैं: अवकलन erx\eu^{rx} पर संख्या rr से गुणन की भाँति काम करता है, अतः अवकल समीकरण संख्यात्मक समीकरण χ(r)=0\chi(r) = 0 बन जाता है — पूरी विश्लेषणात्मक समस्या एक ही द्विघात के मूल खोजने में सिमट जाती है।

मुख्य पद अज्ञात का ऐसा परिवर्तन है जो कोटि घटा देता है। मान लीजिए rr χ\chi का (संभवतः सम्मिश्र) मूल है और y=zerxy = z\, \eu^{rx} लिखिए, जिससे व्यापकता नहीं खोती। तब

y+ay+by=(z+(2r+a)z+χ(r)z)erx=(z+(2r+a)z)erx,y'' + ay' + by = \bigl(z'' + (2r + a) z' + \chi(r) z\bigr)\eu^{rx} = \bigl(z'' + (2r + a) z'\bigr)\eu^{rx},

अतः (H)(H) u=zu = z' के लिए प्रथम कोटि का समीकरण u+(2r+a)u=0u' + (2r + a) u = 0 बन जाता है।

स्थिति Δ0\Delta \neq 0: r=r1r = r_1 चुनिए; तब 2r1+a=r1r22r_1 + a = r_1 - r_2 (क्योंकि r1+r2=ar_1 + r_2 = -a)। प्रमेय 5.2 से, किसी अचर cc के लिए z=ce(r2r1)xz' = c\,\eu^{(r_2 - r_1)x}; R\R पर समाकलन करने पर μ=cr2r1\mu = \frac{c}{r_2 - r_1} के साथ z=μe(r2r1)x+λz = \mu\, \eu^{(r_2 - r_1)x} + \lambda, और इसलिए y=zer1x=λer1x+μer2xy = z\, \eu^{r_1 x} = \lambda\,\eu^{r_1x} + \mu\,\eu^{r_2x}। जब Δ<0\Delta < 0, तब मूल α±iω\alpha \pm \iu\omega हैं और सम्मिश्र हल c1,c2Cc_1, c_2 \in \C के साथ y=c1e(α+iω)x+c2e(αiω)xy = c_1\eu^{(\alpha+\iu\omega)x} + c_2\eu^{(\alpha-\iu\omega)x} हैं। इनमें से कौन वास्तविक-मान वाले हैं? चूँकि e(α+iω)x=e(αiω)x\conj{\eu^{(\alpha+\iu\omega)x}} = \eu^{(\alpha-\iu\omega)x}, yy का संयुग्मी c2e(α+iω)x+c1e(αiω)x\conj{c_2}\, \eu^{(\alpha+\iu\omega)x} + \conj{c_1}\,\eu^{(\alpha-\iu\omega)x} है, और सभी xx के लिए y=yy = \conj y c2=c1c_2 = \conj{c_1} पर बाध्य कर देता है (दोनों चरघातांकी रैखिकतः स्वतंत्र हैं: दो बिंदुओं पर मान लीजिए, अथवा eαx\eu^{\alpha x} से भाग देकर x=0x=0 पर तुलना कीजिए)। वास्तविक λ,μ\lambda, \mu के साथ c1=λiμ2c_1 = \frac{\lambda - \iu\mu}2 लिखने पर:

y=2(c1eαx(cosωx+isinωx))=eαx(λcosωx+μsinωx),y = 2\,\Re\Bigl(c_1\,\eu^{\alpha x}(\cos\omega x + \iu\sin\omega x)\Bigr) = \eu^{\alpha x}\bigl(\lambda\cos\omega x + \mu\sin\omega x\bigr),

और विलोमतः ऐसा हर फलन एक हल है (वास्तविक समीकरण के सम्मिश्र हल का वास्तविक भाग): अतः वास्तविक हल-समष्टि वही है जो घोषित की गई थी।

स्थिति Δ=0\Delta = 0: r=r0r = r_0, 2r0+a=02r_0 + a = 0, अतः z=0z'' = 0: z=λ+μxz = \lambda + \mu x और y=(λ+μx)er0xy = (\lambda + \mu x)\eu^{r_0 x}

उदाहरण 5.11 (एक कोशी समस्या, आरंभ से अंत तक)

y(0)=0y(0) = 0, y(0)=1y'(0) = 1 के साथ y3y+2y=0y'' - 3y' + 2y = 0 हल कीजिए। अभिलक्षणिक बहुपद r23r+2=(r1)(r2)r^2 - 3r + 2 = (r - 1)(r - 2) के वास्तविक मूल 11 और 22 हैं: व्यापक हल y=λex+μe2xy = \lambda\eu^{x} + \mu\eu^{2x}। दोनों शर्तें रैखिक निकाय देती हैं

λ+μ=0,λ+2μ=1,\lambda + \mu = 0, \qquad \lambda + 2\mu = 1 ,

अतः μ=1\mu = 1, λ=1\lambda = -1:

y(x)=e2xex.y(x) = \eu^{2x} - \eu^{x} .

जाँच: y(0)=0y(0) = 0; y=2e2xexy' = 2\eu^{2x} - \eu^x के लिए y(0)=1y'(0) = 1; और y3y+2y=(46+2)e2x+(1+32)ex=0y'' - 3y' + 2y = (4 - 6 + 2)\eu^{2x} + (-1 + 3 - 2)\eu^x = 0। उत्तर के आकार पर ध्यान दीजिए: -\infty के पास धीमी विधा ex-\eu^x प्रभावी है; ++\infty के पास तेज़ विधा e2x\eu^{2x}। हलों को भिन्न क्षय या वृद्धि दरों वाली विधाओं के अध्यारोपण के रूप में पढ़ना लाभदायक आदत है — सप्ताहांत समस्या का क्षणिक/स्थायी विभाजन इसी प्रकार सँभाला गया है।

क्षयशील हलों वाले y'' + ay' + by = 0 की तीन अवस्थाएँ: अवमंदित दोलन (सम्मिश्र मूल), क्रांतिक वापसी (द्विक मूल), अति-अवमंदित क्षय (दो वास्तविक मूल)। कौन-सी अवस्था होगी, यह केवल = a2 - 4b के चिह्न से पढ़ ली जाती है — कुछ भी हल करने से पहले।
क्षयशील हलों वाले y+ay+by=0y'' + ay' + by = 0 की तीन अवस्थाएँ: अवमंदित दोलन (सम्मिश्र मूल), क्रांतिक वापसी (द्विक मूल), अति-अवमंदित क्षय (दो वास्तविक मूल)। कौन-सी अवस्था होगी, यह केवल Δ=a24b\Delta = a^2 - 4b के चिह्न से पढ़ ली जाती है — कुछ भी हल करने से पहले।

प्रमेय 5.12 (संरचना और कोशी समस्या)

  1. यदि ypy_p (E)(E) का एक विशिष्ट हल है, तो (E)(E) के हल ठीक yp+yhy_p + y_h हैं, जहाँ yhy_h (H)(H) के हलों पर चलता है।
  2. (अध्यारोपण) यदि y1y_1 y+ay+by=f1y'' + ay' + by = f_1 को हल करता है और y2y_2 y+ay+by=f2y'' + ay' + by = f_2 को, तो y1+y2y_1 + y_2 दाएँ पक्ष f1+f2f_1 + f_2 वाले समीकरण को हल करता है।
  3. सभी x0Ix_0 \in I और (y0,y0)(y_0, y_0') के लिए कोशी समस्या(E)(E), y(x0)=y0y(x_0) = y_0, y(x0)=y0y'(x_0) = y_0'” का II पर ठीक एक हल है। (अस्तित्व, एक विशिष्ट हल मान लेने पर; अद्वितीयता पूरी तरह।)

उपपत्ति. (1) yy+ay+byy \mapsto y'' + ay' + by की रैखिकता से yy (E)(E) को हल करता है यदि और केवल यदि yypy - y_p (H)(H) को हल करता है। (2) वही रैखिकता है।

(3) (1) से यह सिद्ध करना पर्याप्त है कि अचर (λ,μ)(\lambda, \mu) किसी भी आँकड़े (y0,y0)(y_0, y_0') के अनुसार सदा, और अद्वितीय रूप से, समायोजित किए जा सकते हैं। चर का स्थानांतरण करके x0=0x_0 = 0 मान लीजिए। प्रमेय 5.10 की स्थिति (1) में y=λer1x+μer2xy = \lambda\eu^{r_1x} + \mu\eu^{r_2x} देता है

y(0)=λ+μ,y(0)=r1λ+r2μ:y(0) = \lambda + \mu, \qquad y'(0) = r_1\lambda + r_2\mu :

जो (λ,μ)(\lambda, \mu) में ऐसा रैखिक निकाय है जिसका सारणिक r2r10r_2 - r_1 \neq 0 है; उसे स्पष्ट रूप से हल करने पर μ=y0r1y0r2r1\mu = \frac{y_0' - r_1y_0}{r_2 - r_1} और λ=y0μ\lambda = y_0 - \mu: ठीक एक हल। स्थिति (2) में y(0)=λy(0) = \lambda और y(0)=r0λ+μy'(0) = r_0\lambda + \mu: निकाय त्रिभुजीय है जिसका सारणिक 11 है, और λ=y0\lambda = y_0, μ=y0r0y0\mu = y_0' - r_0y_0 से हल हो जाता है। स्थिति (3) में y(0)=λy(0) = \lambda और y(0)=αλ+ωμy'(0) = \alpha\lambda + \omega\mu: सारणिक ω0\omega \neq 0, और λ=y0\lambda = y_0, μ=y0αy0ω\mu = \frac{y_0' - \alpha y_0}\omega से हल। हर स्थिति में प्रतिचित्रण (λ,μ)(y(x0),y(x0))(\lambda, \mu) \mapsto (y(x_0), y'(x_0)) एक रैखिक एकैकी आच्छादन है — अध्याय 20 की भाषा इस स्थिति-दर-स्थिति जाँच को एक वाक्य में सिकोड़ देगी।

विधि 5.13 (f(x)=P(x)eγxf(x) = P(x)\,\eu^{\gamma x} के लिए विशिष्ट हल)

जब दायाँ पक्ष P(x)eγxP(x)\,\eu^{\gamma x} हो, जहाँ PP बहुपद है और γR\gamma \in \R (यह बहुपदों, चरघातांकियों और सम्मिश्र γ\gamma या अध्यारोपण के द्वारा cos\cos तथा sin\sin को समेट लेता है), तब इस रूप का विशिष्ट हल खोजिए

yp(x)=xmQ(x)eγx,m=बहुपद χ के मूल के रूप में γ की बहुलता (m=0,1 या 2),y_p(x) = x^{m}\, Q(x)\, \eu^{\gamma x}, \qquad m = \text{बहुपद } \chi \text{ के मूल के रूप में } \gamma \text{ की बहुलता} \ (m = 0, 1 \text{ या } 2),

जहाँ QQ PP के समान घात का बहुपद है, जिसके गुणांक प्रतिस्थापन और पहचान से मिलते हैं। f=Kcosωxf = K\cos\omega x (या sin\sin) के लिए दाएँ पक्ष KeiωxK\eu^{\iu\omega x} के साथ हल कीजिए और वास्तविक (क्रमशः काल्पनिक) भाग लीजिए।

उदाहरण 5.14 (अध्यारोपण काम पर)

R\R पर yy=ex+4y'' - y = \eu^{x} + 4 हल कीजिए। समांगी: χ(r)=r21\chi(r) = r^2 - 1, मूल ±1\pm1, अतः yh=λex+μexy_h = \lambda\eu^x + \mu\eu^{-x}। दाएँ पक्ष को तोड़िए और हर टुकड़े का उपचार विधि-बॉक्स से कीजिए। टुकड़ा ex\eu^x: यहाँ γ=1\gamma = 1 χ\chi का सरल मूल है, अतः y1=cxexy_1 = c\,x\,\eu^x आज़माइए: तब y1y1=c(x+2)excxex=2cexy_1'' - y_1 = c(x + 2)\eu^x - cx\eu^x = 2c\,\eu^x, जिससे c=12c = \frac12 मिलता है। टुकड़ा 44: γ=0\gamma = 0 मूल नहीं है; अचर y2=4y_2 = -4 काम करता है। अध्यारोपण (प्रमेय 5.12 (2)) से:

y=xex24+λex+μex,(λ,μ)R2.y = \frac{x\,\eu^x}2 - 4 + \lambda\,\eu^x + \mu\,\eu^{-x}, \qquad (\lambda, \mu) \in \R^2 .

ध्यान दीजिए कि दोनों टुकड़ों ने भिन्न आकार माँगे (m=1m = 1 बनाम m=0m = 0): बहुलता की जाँच प्रत्येक घातांक पर अलग-अलग लगती है, और अनुमान लगाने से पहले दाएँ पक्ष को तोड़ने का यही पूरा प्रयोजन है।

उदाहरण 5.15 (बहुलता का नियम काम पर)

R\R पर y+y=xy'' + y' = x हल कीजिए। दायाँ पक्ष P(x)=xP(x) = x के साथ P(x)e0xP(x)\eu^{0 \cdot x} है, और γ=0\gamma = 0 χ(r)=r2+r=r(r+1)\chi(r) = r^2 + r = r(r + 1) का सरल मूल है: अतः m=1m = 1, और सही अनुमान yp=x(αx+β)=αx2+βxy_p = x\,(\alpha x + \beta) = \alpha x^2 + \beta x है, जो PP से एक घात ऊपर है। प्रतिस्थापित करने पर:

yp+yp=2α+(2αx+β)=2αx+(2α+β),y_p'' + y_p' = 2\alpha + (2\alpha x + \beta) = 2\alpha x + (2\alpha + \beta) ,

और xx से पहचान करने पर α=12\alpha = \frac12, β=1\beta = -1 मिलते हैं: yp=x22xy_p = \frac{x^2}2 - x। व्यापक हल: y=x22x+λ+μexy = \frac{x^2}2 - x + \lambda + \mu\,\eu^{-x}। यदि हमने yp=αx+βy_p = \alpha x + \beta का अनुमान लगाया होता (बहुलता की उपेक्षा करके), तो प्रतिस्थापन yp+yp=αy_p'' + y_p' = \alpha देता, जो अचर है — α,β\alpha, \beta का कोई चयन xx से मेल नहीं खा सकता, और यह विफलता संरचनात्मक है: अचर तो पहले से ही समांगी समीकरण को हल करते हैं, अतः वे बाएँ पक्ष के लिए अदृश्य हैं। गुणनखंड xmx^m ठीक इसीलिए है कि समांगी हल-समष्टि से बाहर निकला जा सके।

टिप्पणी 5.16 (तीस सेकंड की बीमा-पॉलिसी)

इस अध्याय का हर हल किया गया समीकरण प्रतिस्थापन की जाँच पर समाप्त होता है, और यह सजावट नहीं है। अवकल समीकरण की संगणना अनेक छोटे पदों को जोड़ती है (एक प्रतिअवकलज, एक गुणनफल नियम, दो अचर), और एक ही चिह्न की भूल अदृश्य रूप से आगे बढ़ती जाती है; अंतिम सूत्र को समीकरण में वापस रखने पर एक अवकलन के मूल्य पर लगभग सारी भूलें पकड़ में आ जाती हैं। यह प्रतिवर्त तीन परतों में विकसित कीजिए: अकेले विशिष्ट हल की जाँच कीजिए (समांगी भाग तो वैसे भी कट जाता है), पूरे हल पर आरंभिक शर्तों की जाँच कीजिए, और जब कोई प्राचल हो तब किसी अपकर्षी मान की जाँच कीजिए (क्या व्यापक Ω\Omega का सूत्र Ω=0\Omega = 0 पर ज्ञात उत्तर दुबारा दे देता है?)। इस आदत की कीमत आधा मिनट है; और वह “शायद ठीक” को “सत्यापित” में बदल देती है।

टिप्पणी 5.17 (सामान्य भूलें)

  1. पहले मानकीकरण कीजिए। सूत्र मानकर चलते हैं कि समीकरण y+a(x)y=b(x)y' + a(x)y = b(x) है — yy' पर गुणांक 11 के साथ। xy2y=x3xy' - 2y = x^3 के लिए aa और bb की पहचान करने से पहले xx से भाग दीजिए (00 से बचने वाले अंतराल पर), जैसा अभ्यास 5.2 में है।
  2. प्रति विमा एक अचर, और वह अंत में नियत। प्रथम कोटि का व्यापक हल एक अचर रखता है, द्वितीय कोटि का दो; आरंभिक शर्तें पूरे हल yp+yhy_p + y_h पर लगाई जाती हैं, कभी अकेले yhy_h पर नहीं — ypy_p जोड़ने से पहले उन्हें लगा देना सबसे बार-बार होने वाली संरचनात्मक भूल है।
  3. बहुलता का ध्यान रखिए। जो विशिष्ट-हल अनुमान समांगी समीकरण को हल करता हो, वह बाएँ पक्ष के लिए अदृश्य है; विधि 5.13 का गुणनखंड xmx^m वैकल्पिक नहीं है (उदाहरण 5.15)।
  4. अंतराल उत्तर का हिस्सा हैं। हल उन्हीं अंतरालों पर रहते हैं जहाँ गुणांक संतत हों; किसी विचित्रता को पार करके जोड़ने से मिथ्या अचर बन सकते हैं (अभ्यास 5.12) या अद्वितीयता नष्ट हो सकती है। “R\R^* पर हल कीजिए” का अर्थ है दो स्वतंत्र समस्याएँ।

उदाहरण 5.18 (अनुनाद से हटकर प्रणोदन)

R\R पर y+4y=sinxy'' + 4y = \sin x हल कीजिए। स्वाभाविक आवृत्ति 22, प्रणोदन आवृत्ति 11: चूँकि i\iu χ(r)=r2+4\chi(r) = r^2 + 4 का मूल नहीं है, बहुलता m=0m = 0 है और एक सादा ज्यावक्रीय फलन पर्याप्त है। yp=αsinxy_p = \alpha\sin x आज़माइए (कोज्या की आवश्यकता नहीं: समीकरण में yy' वाला पद नहीं है, और sin\sin sin\sin को फिर से उत्पन्न कर देता है):

yp+4yp=αsinx+4αsinx=3αsinx,y_p'' + 4y_p = -\alpha\sin x + 4\alpha\sin x = 3\alpha\sin x ,

अतः α=13\alpha = \frac13 और व्यापक हल है

y=sinx3+λcos2x+μsin2x.y = \frac{\sin x}3 + \lambda\cos 2x + \mu\sin 2x .

प्रत्येक हल परिबद्ध रहता है: आवृत्तियों 11 (प्रणोदित) और 22 (स्वाभाविक) के दो दोलनों का अध्यारोपण। अगले उदाहरण से तुलना कीजिए, जहाँ स्वाभाविक आवृत्ति पर प्रणोदन उत्तर का आकार ही बदल देता है।

उदाहरण 5.19 (एक प्रणोदित दोलन)

y+y=cosxy'' + y = \cos x, y(0)=0y(0) = 0, y(0)=0y'(0) = 0 हल कीजिए।

समांगी: χ(r)=r2+1\chi(r) = r^2 + 1, मूल ±i\pm\iu: yh=λcosx+μsinxy_h = \lambda\cos x + \mu \sin x

विशिष्ट: दायाँ पक्ष (eix)\Re(\eu^{\iu x}) है, जहाँ γ=i\gamma = \iu χ\chi का सरल मूल है: zp=cxeixz_p = c\, x\, \eu^{\iu x} आज़माइए (cCc \in \C)। तब zp+zp=c(2i)eixz_p'' + z_p = c\,(2\iu)\eu^{\iu x}, जो c=12i=i2c = \frac{1}{2\iu} = -\frac\iu2 के लिए eix\eu^{\iu x} के बराबर है। अतः zp=i2x(cosx+isinx)z_p = -\frac{\iu}{2} x (\cos x + \iu \sin x) और yp=(zp)=xsinx2y_p = \Re(z_p) = \frac{x \sin x}{2}

व्यापक हल: y=xsinx2+λcosx+μsinxy = \frac{x\sin x}{2} + \lambda\cos x + \mu\sin x। शर्तें: y(0)=λ=0y(0) = \lambda = 0; y=sinx+xcosx2+μcosxy' = \frac{\sin x + x\cos x}{2} + \mu\cos x, अतः y(0)=μ=0y'(0) = \mu = 0। उत्तर: y=xsinx2y = \frac{x\sin x}{2} — ऐसा दोलन जिसका आयाम रैखिक रूप से बढ़ता है: यही अनुनाद की परिघटना है, जो निकाय को उसकी स्वाभाविक आवृत्ति पर प्रणोदित करने से होती है।

अनुनाद: y'' + y = x का हल y = x x/2 रेखाओं y = ± x2 (बिंदुकित) के बीच, निरंतर बढ़ते आयाम के साथ दोलन करता है।
अनुनाद: y+y=cosxy'' + y = \cos x का हल y=xsinx2y = \frac{x \sin x}{2} रेखाओं y=±x2y = \pm\frac x2 (बिंदुकित) के बीच, निरंतर बढ़ते आयाम के साथ दोलन करता है।

टिप्पणी 5.20 (मध्यांतर: रैखिकता एक ज्यामिति है)

इस अध्याय के हर हल-समुच्चय के आकार पर फिर दृष्टि डालिए: एक विशेष हल, और उसके साथ समांगी हलों की ऐसी समष्टि जिसमें एक मुक्त अचर (प्रथम कोटि) या दो (द्वितीय कोटि) हैं। रैखिक बीजगणित के अध्याय (अध्याय 18, 19 और 20) ठीक-ठीक शब्दावली देंगे: प्रतिचित्रण L(y)=y+ay+byL(y) = y'' + ay' + by रैखिक है, उसके समांगी हल LL की अष्टि बनाते हैं, जो ऐसी सदिश समष्टि है जिसकी विमा समीकरण की कोटि के बराबर है — “प्रति कोटि एक अचर” की ईमानदार सामग्री यही है — और L(y)=fL(y) = f का हल-समुच्चय एक ऐफ़ीन उपसमष्टि है, अर्थात् अष्टि का स्थानांतरण। यहाँ तक कि प्रमेय 5.12 का कोशी प्रतिचित्रण (λ,μ)(y(x0),y(x0))(\lambda, \mu) \mapsto (y(x_0), y'(x_0)) भी दो समतलों के बीच रैखिक एकैकी आच्छादन है, अर्थात् एक व्युत्क्रमणीय 2×22 \times 2 निकाय (अध्याय 21)। इस अध्याय में कुछ भी दोबारा करने की आवश्यकता नहीं होगी — केवल उसका नाम बदलना होगा, और यह नामकरण रैखिक बीजगणित के लिए सर्वोत्तम संभव तैयारी है: वहाँ की हर अमूर्त परिभाषा यहाँ पहले ही अपनी रोटी कमा चुकी है।

टिप्पणी 5.21 (यह अध्याय कहाँ काम आता है)

संरचना प्रमेय — (E)(E) के हल “एक विशिष्ट हल और (H)(H) के हल” बनाते हैं — उस प्रतिरूप की पहली उपस्थिति है जिसे अध्याय 18 और 20 नाम देगा: (H)(H) का हल-समुच्चय रैखिक प्रतिचित्रण yy+ay+byy \mapsto y'' + ay' + by की अष्टि है, और (E)(E) का हल-समुच्चय उसका ऐफ़ीन स्थानांतरण। अभिलक्षणिक बहुपद अध्याय 21 में आव्यूह के अभिलक्षणिक बहुपद के रूप में फिर प्रकट होता है: द्वितीय कोटि का समीकरण वेश बदला हुआ 2×22 \times 2 प्रथम कोटि का निकाय है, और स्नातक वर्ष 2 का खंड इस दृष्टिकोण को व्यवस्थित करता है। अचरों के विचरण की माँगी हुई समाकलों की आपूर्ति अध्याय 15 करता है, और नीचे दी गई सप्ताहांत समस्या — प्रणोदित अवमंदित दोलक — विज्ञानों के हर दोलन-प्रश्न का आदर्श स्थिति है, परिपथों से लेकर झूलते पुलों तक।

5.3 अभ्यास

अभ्यास 5.1

R\R पर हल कीजिए:   y+2y=e3x\;y' + 2y = \eu^{3x}; फिर कोशी समस्या y(0)=1y(0) = 1

हल

हल — अभ्यास 5.1.

समांगी: yh=λe2xy_h = \lambda\,\eu^{-2x}। विशिष्ट: yp=ce3xy_p = c\,\eu^{3x} आज़माइए (33 r+2r + 2 का मूल नहीं है): 3c+2c=13c + 2c = 1, c=15c = \frac15। व्यापक हल: y=e3x5+λe2xy = \frac{\eu^{3x}}{5} + \lambda\,\eu^{-2x}y(0)=1y(0) = 1 के साथ: 15+λ=1\frac15 + \lambda = 1, λ=45\lambda = \frac45, अतः y=e3x+4e2x5y = \frac{\eu^{3x} + 4\,\eu^{-2x}}{5}

अभ्यास 5.2

(0,+)\intoo{0}{+\infty} पर हल कीजिए:   xy2y=x3\;x y' - 2y = x^3 (पहले समीकरण को मानकीकृत रूप में रखिए)

हल

हल — अभ्यास 5.2.

(0,+)\intoo{0}{+\infty} पर xx से भाग दीजिए: y2xy=x2y' - \frac{2}{x}\,y = x^2। यहाँ A(x)=2lnxA(x) = -2\ln x, eA(x)=x2\eu^{-A(x)} = x^2: समांगी हल λx2\lambda x^2अचरों का विचरण: μ(x)=x2x2=1\mu'(x) = x^2 \cdot x^{-2} = 1, अतः μ=x+λ\mu = x + \lambda और

y(x)=x3+λx2,λR.y(x) = x^3 + \lambda x^2, \qquad \lambda \in \R .

जाँच: x(3x2+2λx)2(x3+λx2)=x3x(3x^2 + 2\lambda x) - 2(x^3 + \lambda x^2) = x^3

अभ्यास 5.3

R\R पर वास्तविक व्यापक हल देते हुए हल कीजिए:   y3y+2y=0\;y'' - 3y' + 2y = 0;   y+4y+4y=0\;y'' + 4y' + 4y = 0;   y2y+5y=0\;y'' - 2y' + 5y = 0

हल

हल — अभ्यास 5.3.

y3y+2y=0y'' - 3y' + 2y = 0: मूल 11 और 22;   y=λex+μe2x\;y = \lambda\,\eu^{x} + \mu\,\eu^{2x}

y+4y+4y=0y'' + 4y' + 4y = 0: द्विक मूल 2-2;   y=(λ+μx)e2x\;y = (\lambda + \mu x)\,\eu^{-2x}

y2y+5y=0y'' - 2y' + 5y = 0: मूल 1±2i1 \pm 2\iu;   y=ex(λcos2x+μsin2x)\;y = \eu^{x}(\lambda\cos 2x + \mu\sin 2x)

अभ्यास 5.4

R\R पर yy=x2y'' - y = x^2 हल कीजिए, फिर कोशी समस्या y(0)=0y(0) = 0, y(0)=1y'(0) = 1

हल

हल — अभ्यास 5.4.

समांगी: मूल ±1\pm 1, yh=λex+μexy_h = \lambda\,\eu^x + \mu\,\eu^{-x}। बहुपदीय दाएँ पक्ष के साथ विशिष्ट (γ=0\gamma = 0 मूल नहीं है): yp=ax2+bx+cy_p = ax^2 + bx + c; प्रतिस्थापित करने पर 2a(ax2+bx+c)=x22a - (ax^2 + bx + c) = x^2 से a=1a = -1, b=0b = 0, c=2a=2c = 2a = -2 मिलते हैं: yp=x22y_p = -x^2 - 2। व्यापक हल y=x22+λex+μexy = -x^2 - 2 + \lambda\eu^x + \mu\eu^{-x}

कोशी: y(0)=2+λ+μ=0y(0) = -2 + \lambda + \mu = 0 और y(0)=λμ=1y'(0) = \lambda - \mu = 1: λ=32\lambda = \frac32, μ=12\mu = \frac12। अतः y=x22+3ex+ex2y = -x^2 - 2 + \frac{3\eu^x + \eu^{-x}}{2}

अभ्यास 5.5 ★★

(π2,π2)\intoo{-\frac\pi2}{\frac\pi2} पर हल कीजिए:   y+ytanx=sin2x\;y' + y\tan x = \sin 2x

हल

हल — अभ्यास 5.5.

a(x)=tanxa(x) = \tan x, A(x)=ln(cosx)A(x) = -\ln(\cos x) (वैध: अंतराल पर cos>0\cos > 0), eA=cosx\eu^{-A} = \cos x: समांगी हल λcosx\lambda\cos xअचरों का विचरण: μ(x)=sin2x1cosx=2sinx\mu'(x) = \sin 2x \cdot \frac{1}{\cos x} = 2\sin x, अतः μ=2cosx+λ\mu = -2\cos x + \lambda और

y(x)=2cos2x+λcosx.y(x) = -2\cos^2 x + \lambda \cos x .

जाँच: y=4cosxsinxλsinxy' = 4\cos x \sin x - \lambda\sin x और ytanx=2cosxsinx+λsinxy\tan x = -2\cos x\sin x + \lambda \sin x; उनका योग 2cosxsinx=sin2x2\cos x\sin x = \sin 2x है, जैसा अपेक्षित था।

अभ्यास 5.6 ★★

R\R पर y4y+3y=(2x+1)exy'' - 4y' + 3y = (2x + 1)\,\eu^{x} हल कीजिए। (बहुलता का ध्यान रखिए: क्या 11 अभिलक्षणिक बहुपद का मूल है?)

हल

हल — अभ्यास 5.6.

χ(r)=r24r+3=(r1)(r3)\chi(r) = r^2 - 4r + 3 = (r-1)(r-3): γ=1\gamma = 1 सरल मूल है (m=1m = 1)। yp=x(ax+b)exy_p = x(ax + b)\,\eu^x आज़माइए। u=ax2+bxu = ax^2 + bx के साथ,

yp4yp+3yp=(u+(24)u+χ(1)u)ex=(2a2(2ax+b))ex.y_p'' - 4y_p' + 3y_p = \bigl(u'' + (2 - 4)u' + \chi(1) u\bigr)\eu^x = \bigl(2a - 2(2ax + b)\bigr)\eu^x .

(2x+1)ex(2x + 1)\eu^x से पहचान कीजिए: 4a=2-4a = 2 और 2a2b=12a - 2b = 1, अतः a=12a = -\frac12, b=1b = -1। व्यापक हल:

y=(x22+x)ex+λex+μe3x,(λ,μ)R2.y = -\Bigl(\frac{x^2}{2} + x\Bigr)\eu^{x} + \lambda\,\eu^{x} + \mu\,\eu^{3x}, \qquad (\lambda, \mu) \in \R^2 .

अभ्यास 5.7 ★★

R\R पर y+4y=sin2x+xy'' + 4y = \sin 2x + x हल कीजिए (अध्यारोपण; प्रत्येक दाएँ पक्ष का उपचार अलग-अलग कीजिए)

हल

हल — अभ्यास 5.7.

समांगी: yh=λcos2x+μsin2xy_h = \lambda\cos 2x + \mu\sin 2x

दायाँ पक्ष xx (γ=0\gamma = 0 मूल नहीं है): 4(ax+b)=x4(ax + b) = x के साथ y1=ax+by_1 = ax + b: y1=x4y_1 = \frac x4

दायाँ पक्ष sin2x=(e2ix)\sin 2x = \Im(\eu^{2\iu x}), जहाँ 2i2\iu r2+4r^2 + 4 का सरल मूल है: z=cxe2ixz = c\,x\,\eu^{2\iu x} आज़माइए; तब z+4z=4ice2ixz'' + 4z = 4\iu c\,\eu^{2\iu x}, जो c=14i=i4c = \frac{1}{4\iu} = -\frac{\iu}{4} के लिए e2ix\eu^{2\iu x} के बराबर है। अतः z=ix4(cos2x+isin2x)z = -\frac{\iu x}{4}(\cos 2x + \iu \sin 2x) और y2=(z)=xcos2x4y_2 = \Im(z) = -\frac{x\cos 2x}{4}

अध्यारोपण से:

y=x4xcos2x4+λcos2x+μsin2x.y = \frac{x}{4} - \frac{x\cos 2x}{4} + \lambda\cos 2x + \mu\sin 2x .

अभ्यास 5.8 ★★

T0=80T_0 = 80\,^\circC तापमान वाली कॉफ़ी का प्याला 2020\,^\circC के कमरे में रखा है। न्यूटन का शीतलन नियम k>0k > 0 के साथ T=k(T20)T' = -k\,(T - 20) कहता है। T(t)T(t) के लिए हल कीजिए, और यह देखते हुए कि 1010 मिनट बाद कॉफ़ी 5050\,^\circC पर है, बताइए कि वह 2525\,^\circC पर कब पहुँचेगी।

हल

हल — अभ्यास 5.8.

समीकरण T+kT=20kT' + kT = 20k का अचर विशिष्ट हल 2020 है और समांगी हल λekt\lambda\eu^{-kt} हैं: T(t)=20+λektT(t) = 20 + \lambda\,\eu^{-kt}, और T(0)=80T(0) = 80 से λ=60\lambda = 60 मिलता है:

T(t)=20+60ekt.T(t) = 20 + 60\,\eu^{-kt} .

T(10)=50T(10) = 50: e10k=12\eu^{-10k} = \frac12, अतः k=ln210k = \frac{\ln 2}{10}। फिर T(t)=25T(t) = 25 के लिए ekt=560=112\eu^{-kt} = \frac{5}{60} = \frac{1}{12} आवश्यक है, अर्थात्

t=ln12k=10ln12ln235.8 मिनट।t = \frac{\ln 12}{k} = 10\,\frac{\ln 12}{\ln 2} \approx 35.8 \text{ मिनट।}

अभ्यास 5.9 ★★★

(अवमंदित दोलक) ε0\varepsilon \geq 0 के लिए y+2εy+y=0y'' + 2\varepsilon y' + y = 0 पर विचार कीजिए।

  1. ε[0,1)\varepsilon \in \intco{0}{1}, ε=1\varepsilon = 1 और ε>1\varepsilon > 1 के लिए हल कीजिए।
  2. दिखाइए कि ε>0\varepsilon > 0 के लिए प्रत्येक हल ++\infty पर 00 की ओर जाता है, और ε=0\varepsilon = 0 के लिए अशून्य हल ऐसा नहीं करते।
  3. ε(0,1)\varepsilon \in \intoo{0}{1} के लिए दिखाइए कि किसी अशून्य हल के शून्य नियमित अंतराल पर हैं, जिनका अंतराल π1ε2\frac{\pi}{\sqrt{1 - \varepsilon^2}} है।
हल

हल — अभ्यास 5.9.

  1. χ(r)=r2+2εr+1\chi(r) = r^2 + 2\varepsilon r + 1, Δ=4(ε21)\Delta = 4(\varepsilon^2 - 1)ε[0,1)\varepsilon \in \intco{0}{1} के लिए: मूल ε±i1ε2-\varepsilon \pm \iu\sqrt{1 - \varepsilon^2}, अतः ω=1ε2\omega = \sqrt{1 - \varepsilon^2} के साथ y=eεt(λcosωt+μsinωt)y = \eu^{-\varepsilon t}\bigl(\lambda\cos\omega t + \mu\sin\omega t\bigr)ε=1\varepsilon = 1 के लिए: द्विक मूल 1-1, y=(λ+μt)ety = (\lambda + \mu t)\,\eu^{-t}ε>1\varepsilon > 1 के लिए: वास्तविक मूल r±=ε±ε21r_\pm = -\varepsilon \pm \sqrt{\varepsilon^2 - 1}, दोनों <0< 0, और y=λer+t+μerty = \lambda\eu^{r_+t} + \mu\eu^{r_-t}
  2. ε(0,1)\varepsilon \in \intoo{0}{1} के लिए: yeεt(λ+μ)0\abs y \leq \eu^{-\varepsilon t}(\abs\lambda + \abs\mu) \to 0ε=1\varepsilon = 1 के लिए: (λ+μt)et0(\lambda + \mu t)\eu^{-t} \to 0 (चरघातांकी बहुपद को हरा देता है, प्रतिज्ञप्ति 4.6)। ε>1\varepsilon > 1 के लिए: दोनों चरघातांकी क्षय होते हैं, क्योंकि r±<0r_\pm < 0 (वस्तुतः ε21<ε\sqrt{\varepsilon^2 - 1} < \varepsilon)। ε=0\varepsilon = 0 के लिए: y=λcost+μsinty = \lambda\cos t + \mu\sin t का आयाम अचर λ2+μ20\sqrt{\lambda^2 + \mu^2} \neq 0 रहता है, जब तक y=0y = 0 न हो।
  3. R=λ2+μ2>0R = \sqrt{\lambda^2 + \mu^2} > 0 के साथ λcosωt+μsinωt=Rcos(ωtφ)\lambda\cos\omega t + \mu\sin\omega t = R\cos(\omega t - \varphi) लिखिए। yy के शून्य वही हैं जो cos(ωtφ)\cos(\omega t - \varphi) के हैं (गुणनखंड eεt\eu^{-\varepsilon t} कभी लुप्त नहीं होता): ωtφπ2(modπ)\omega t - \varphi \equiv \frac\pi2 \pmod \pi, जो अंतराल πω=π1ε2\frac{\pi}{\omega} = \frac{\pi}{\sqrt{1 - \varepsilon^2}} वाली समांतर श्रेणी है।

अभ्यास 5.10 ★★★

ऐसे सभी दो बार अवकलनीय फलन f ⁣:RRf \colon \R \to \R खोजिए कि

x,yR,f(x+y)+f(xy)=2f(x)f(y),\forall x, y \in \R, \qquad f(x + y) + f(x - y) = 2 f(x) f(y),

जहाँ f(0)0f(0) \ne 0 और ff अचर नहीं है। संकेत: yy नियत कीजिए, 00 पर xx के सापेक्ष दो बार अवकलन कीजिए; दिखाइए कि किसी अचर cc के लिए f(0)=1f(0) = 1 और f=cff'' = c f; फिर cc के चिह्न के अनुसार हल कीजिए और जाँचिए कि कौन-से हल फलनीय समीकरण को संतुष्ट करते हैं।

हल

हल — अभ्यास 5.10.

x=y=0x = y = 0 रखिए: 2f(0)=2f(0)22f(0) = 2f(0)^2, और f(0)0f(0) \neq 0 से f(0)=1f(0) = 1 मिलता है। xx नियत कीजिए और समीकरण का yy के सापेक्ष दो बार अवकलन कीजिए:

f(x+y)+f(xy)=2f(x)f(y).f''(x+y) + f''(x-y) = 2 f(x) f''(y) .

y=0y = 0 रखने पर:   2f(x)=2f(x)f(0)\;2f''(x) = 2 f(x) f''(0), अर्थात्

f(x)=cf(x),c=f(0).f''(x) = c\,f(x), \qquad c = f''(0).

स्थिति c=ω2>0c = \omega^2 > 0: f(x)=λcoshωx+μsinhωxf(x) = \lambda\cosh\omega x + \mu\sinh\omega x; f(0)=1f(0) = 1 से λ=1\lambda = 1 मिलता है। फलनीय समीकरण में रखने पर और योग-सूत्रों (प्रतिज्ञप्ति 4.18) का प्रयोग करने पर समीकरण μ=0\mu = 0 पर बाध्य कर देता है (sinhωxsinhωy\sinh\omega x \sinh\omega y के गुणांकों की तुलना कीजिए या x=yx = y पर मान लीजिए): f=coshωxf = \cosh\omega x, जो cosh(x+y)+cosh(xy)=2coshxcoshy\cosh(x+y) + \cosh(x-y) = 2\cosh x\cosh y को सचमुच संतुष्ट करता है।

स्थिति c=ω2<0c = -\omega^2 < 0: इसी प्रकार f(x)=cosωxf(x) = \cos\omega x (ω0\omega \neq 0), जो समीकरण को संतुष्ट करता है।

स्थिति c=0c = 0: f(0)=1f(0) = 1 के साथ ff ऐफ़ीन है: f(x)=1+μxf(x) = 1 + \mu x; समीकरण μ=0\mu = 0 पर बाध्य करता है, जो वर्जित है (ff अचर नहीं है)।

निष्कर्ष: हल f(x)=cosωxf(x) = \cos\omega x और f(x)=coshωxf(x) = \cosh\omega x, ω>0\omega > 0 हैं।

अभ्यास 5.11 ★★

(ऑयलर समीकरण) (0,+)\intoo{0}{+\infty} पर x2yxy+y=0x^2 y'' - x y' + y = 0 हल कीजिए। संकेत: z(t)=y(et)z(t) = y(\eu^t) रखिए, अर्थात् x=etx = \eu^t प्रतिस्थापित कीजिए, और दिखाइए कि zz अचर गुणांकों वाले रैखिक समीकरण को संतुष्ट करता है।

हल

हल — अभ्यास 5.11.

z(t)=y(et)z(t) = y(\eu^t) रखिए, जिससे x>0x > 0 के लिए y(x)=z(lnx)y(x) = z(\ln x)। तब

y(x)=z(lnx)x,y(x)=z(lnx)z(lnx)x2,y'(x) = \frac{z'(\ln x)}x, \qquad y''(x) = \frac{z''(\ln x) - z'(\ln x)}{x^2} ,

और समीकरण में प्रतिस्थापित करने पर:

x2yxy+y=(zz)z+z=z2z+z=0.x^2y'' - xy' + y = \bigl(z'' - z'\bigr) - z' + z = z'' - 2z' + z = 0 .

अभिलक्षणिक बहुपद (r1)2(r - 1)^2: द्विक मूल 11, अतः z(t)=(λ+μt)etz(t) = (\lambda + \mu t)\,\eu^t, और चर x=etx = \eu^t में लौटकर:

y(x)=(λ+μlnx)x,λ,μR.y(x) = (\lambda + \mu\ln x)\,x, \qquad \lambda, \mu \in \R .

अभ्यास 5.12 ★★★

अज्ञात अवकलनीय फलन y ⁣:RRy \colon \R \to \R के साथ पूरी वास्तविक रेखा पर समीकरण xy=2yx\,y' = 2y पर विचार कीजिए।

  1. (0,+)\intoo{0}{+\infty} पर और (,0)\intoo{-\infty}{0} पर हल कीजिए।
  2. दिखाइए कि किन्हीं भी अचरों a,bRa, b \in \R के लिए वह फलन जो x0x \geq 0 के लिए ax2ax^2 के बराबर हो और x<0x < 0 के लिए bx2bx^2 के बराबर, R\R पर अवकलनीय है और सर्वत्र समीकरण को हल करता है।
  3. निष्कर्ष निकालिए कि R\R पर हल-समुच्चय दो-प्राचल वाला कुल है, और समझाइए कि यह प्रमेय 5.4 की अद्वितीयता का विरोध क्यों नहीं करता।
हल

हल — अभ्यास 5.12.

  1. प्रत्येक अंतराल पर मानकीकृत रूप y2xy=0y' - \frac2x\,y = 0 में: A(x)=2lnxA(x) = -2\ln\abs x, अतः हल (0,+)\intoo0{+\infty} पर y=ax2y = a x^2 और (,0)\intoo{-\infty}0 पर y=bx2y = b x^2 हैं, जिनके अचर स्वतंत्र हैं (प्रमेय 5.2)।
  2. मान लीजिए x0x \geq 0 के लिए y=ax2y = ax^2 और x<0x < 0 के लिए bx2bx^2। प्रत्येक खुली अर्धरेखा पर yy अवकलनीय है और xy=2yxy' = 2y00 पर: अंतर-भागफल y(h)y(0)h=ah\frac{y(h) - y(0)}h = ah या bhbh 00 की ओर जाते हैं, अतः y(0)=0y'(0) = 0 विद्यमान है, और x=0x = 0 पर समीकरण 0y(0)=2y(0)=00 \cdot y'(0) = 2y(0) = 0 कहता है: जो संतुष्ट है। अतः yy पूरे R\R पर समीकरण को हल करता है।
  3. R\R पर हल ठीक यही जोड़े हुए फलन हैं: प्रथम कोटि के समीकरण के लिए दो-प्राचल वाला कुल। प्रमेय 5.4 से कोई विरोध नहीं है, क्योंकि यहाँ उसकी परिकल्पनाएँ विफल हैं: y+a(x)y=0y' + a(x)y = 0 के रूप में लिखने पर गुणांक a(x)=2xa(x) = -\frac2x 00 पर संतत नहीं है — वस्तुतः परिभाषित ही नहीं है — अतः R\R ऐसा अंतराल नहीं है जिस पर वह प्रमेय लागू हो। 00 की विचित्रता दोनों अर्धरेखाओं को अलग कर देती है, और मान y(0)=0y(0) = 0 अनिवार्य हो जाता है, जो एक ओर से दूसरी ओर कोई सूचना नहीं ले जाता। x00x_0 \neq 0 पर हर कोशी आँकड़ा हल को केवल उसी अर्धरेखा पर निर्धारित करता है जिसमें x0x_0 है।

5.4 समस्या: प्रणोदित अवमंदित दोलक

समस्या 5.1

एक ही समीकरण श्यान माध्यम में स्प्रिंग पर लगे द्रव्यमान, RLC परिपथ के आवेश, और हवा में झूलती इमारत — तीनों को नियंत्रित करता है:

(EΩ) ⁣:x+2λx+ω02x=Acos(Ωt),(E_\Omega)\colon\quad x'' + 2\lambda x' + \omega_0^2\,x = A\cos(\Omega t),

जहाँ λ0\lambda \geq 0 अवमंदन है, ω0>0\omega_0 > 0 स्वाभाविक आवृत्ति है, और A>0A > 0, Ω>0\Omega > 0 प्रणोदन के आयाम तथा आवृत्ति हैं। यह समस्या उसका पूरा व्यवहार निकाल लेती है: क्षणिक भागों का क्षय, अद्वितीय आवर्ती स्थायी अवस्था, अनुनाद वक्र और उसकी तीक्ष्णता (गुणता कारक), अनवमंदित स्थिति के विस्पंद, और वह ऊर्जा-संतुलन जो दोलन को बनाए रखता है। जब तक अन्यथा न कहा जाए, 0<λ<ω00 < \lambda < \omega_0 (अल्प-अवमंदित अवस्था) और हम ωd=ω02λ2\omega_d = \sqrt{\omega_0^2 - \lambda^2} लिखते हैं।

भाग I — मुक्त दोलक। यहाँ A=0A = 0

  1. 0<λ<ω00 < \lambda < \omega_0 के लिए और λ=0\lambda = 0 के लिए समांगी समीकरण (H)(H) हल कीजिए। (अवस्थाएँ λω0\lambda \geq \omega_0 अभ्यास 5.9 में सँभाली गई थीं; उन्हें उद्धृत कीजिए।)
  2. दिखाइए कि प्रत्येक λ>0\lambda > 0 के लिए (H)(H) के सभी हल ++\infty पर 00 की ओर जाते हैं — तीनों अवस्थाओं में।
  3. (H)(H) के किसी हल पर ऊर्जा E(t)=12x(t)2+12ω02x(t)2\mathcal E(t) = \frac12 x'(t)^2 + \frac12\omega_0^2\,x(t)^2 परिभाषित कीजिए। दिखाइए E(t)=2λx(t)20\mathcal E'(t) = -2\lambda\,x'(t)^2 \leq 0, और उससे (बिना कुछ हल किए) निकालिए कि प्रत्येक λ0\lambda \geq 0 के लिए कोशी समस्या(H)(H), x(t0)=x(t0)=0x(t_0) = x'(t_0) = 0” का केवल शून्य हल है।
  4. 0<λ<ω00 < \lambda < \omega_0 के लिए अशून्य हल को x(t)=Reλtcos(ωdtφ)x(t) = R\,\eu^{-\lambda t}\cos(\omega_d t - \varphi) लिखिए और मान लीजिए Td=2πωdT_d = \frac{2\pi}{\omega_d} छद्म-आवर्तकाल है। दिखाइए कि x(t+Td)=eλTdx(t)x(t + T_d) = \eu^{-\lambda T_d}\,x(t): प्रत्येक झूला पिछले झूले का अचर गुणक eδ\eu^{-\delta} से सिकुड़ा हुआ रूप है, δ=2πλωd\delta = \frac{2\pi\lambda}{\omega_d} (लघुगणकीय ह्रासांक)। ω0=1\omega_0 = 1, λ=0.1\lambda = 0.1 के लिए δ\delta संगणित कीजिए।
  5. गुणता कारक Q=ω02λQ = \dfrac{\omega_0}{2\lambda} परिभाषित कीजिए। दिखाइए कि समय 1λ\frac1\lambda (एक आयाम-e\eu-गुणन) के बाद दोलक ωd2πλ\frac{\omega_d}{2\pi\lambda} छद्म-आवर्तकाल पूरे कर चुका होता है, जो दुर्बल अवमंदन (λω0\lambda \ll \omega_0) के लिए लगभग Qπ\frac Q\pi है: गुणता कारक π\pi तक यही गिनता है कि आयाम के e\eu से घटने से पहले कितने दोलन बच रहे।

भाग II — स्थायी अवस्था। अब A>0A > 0 और λ>0\lambda > 0

  1. zCz \in \C (विधि 5.13) के साथ zeiΩtz\,\eu^{\iu\Omega t} के वास्तविक भाग के रूप में विशिष्ट हल खोजिए। दिखाइए कि यह इनके साथ काम करता है

    z=Aω02Ω2+2iλΩ.z = \frac{A}{\omega_0^2 - \Omega^2 + 2\iu\lambda\Omega} .
  2. उससे स्थायी अवस्था को आयाम–कला रूप में निकालिए: xp(t)=R(Ω)cos(Ωtφ(Ω))x_p(t) = R(\Omega)\cos\bigl(\Omega t - \varphi(\Omega)\bigr), जहाँ

    R(Ω)=A(ω02Ω2)2+4λ2Ω2,tanφ=2λΩω02Ω2,φ(0,π).R(\Omega) = \frac{A}{\sqrt{(\omega_0^2 - \Omega^2)^2 + 4\lambda^2\Omega^2}}, \qquad \tan\varphi = \frac{2\lambda\Omega}{\omega_0^2 - \Omega^2}, \quad \varphi \in \intoo0\pi .
  3. दोनों चरम अवस्थाओं की व्याख्या कीजिए: Ω0+\Omega \to 0^+ पर RR और φ\varphi की सीमाएँ संगणित कीजिए (अर्ध-स्थैतिक अनुक्रिया A/ω02A/\omega_0^2, कला 00) और Ω+\Omega \to +\infty पर भी (RA/Ω20R \sim A/\Omega^2 \to 0, कला π\to \pi: द्रव्यमान अत्यधिक तेज़ प्रणोदन के विपरीत चलता है)।
  4. दिखाइए कि (EΩ)(E_\Omega) का प्रत्येक हल xpx_p और (H)(H) के किसी हल का योग है, अतः आरंभिक शर्तें जो भी हों, वह t+t \to +\infty पर स्थायी अवस्था xpx_p की ओर अभिसरित होता है: क्षणिक भाग के मर जाने पर दोलक को यह स्मरण नहीं रहता कि उसका आरंभ कैसे हुआ था।
  5. दिखाइए कि (EΩ)(E_\Omega) का एकमात्र आवर्ती हल xpx_p है।
  6. एक कोशी समस्या अंत तक कीजिए: x(0)=x(0)=0x(0) = x'(0) = 0 के साथ x+2x+2x=costx'' + 2x' + 2x = \cos t के लिए दिखाइए कि हल है

    x(t)=cost+2sint5etcost+3sint5,x(t) = \frac{\cos t + 2\sin t}5 - \eu^{-t}\,\frac{\cos t + 3\sin t}5 ,

    और क्षणिक तथा स्थायी भागों की पहचान कीजिए।

भाग III — अनुनाद वक्र। (0,+)\intoo0{+\infty} पर ΩR(Ω)\Omega \mapsto R(\Omega) का अध्ययन।

  1. u=Ω2u = \Omega^2 और g(u)=(ω02u)2+4λ2ug(u) = (\omega_0^2 - u)^2 + 4\lambda^2 u रखकर दिखाइए: यदि 2λ2<ω022\lambda^2 < \omega_0^2, तो RR अनुनाद आवृत्ति Ωr=ω022λ2\Omega_r = \sqrt{\omega_0^2 - 2\lambda^2} पर कठोर अधिकतम प्राप्त करता है, जहाँ

    Rmax=R(Ωr)=A2λω02λ2.R_{\max} = R(\Omega_r) = \frac{A}{2\lambda\sqrt{\omega_0^2 - \lambda^2}} .
  2. दिखाइए कि RmaxR(0)=Q(1λ2ω02)1/2Q\dfrac{R_{\max}}{R(0)} = Q\,\bigl(1 - \tfrac{\lambda^2}{\omega_0^2}\bigr)^{-1/2} \geq Q: अनुनाद पर प्रणोदन (मूलतः) गुणता कारक से प्रवर्धित होता है।
  3. सिद्ध कीजिए कि वेग का आयाम V(Ω)=ΩR(Ω)V(\Omega) = \Omega\,R(\Omega) ठीक Ω=ω0\Omega = \omega_0 पर अधिकतम है (Ωr\Omega_r पर नहीं), और यह कि वहाँ कला φ(ω0)=π2\varphi(\omega_0) = \frac\pi2 है: Ω=ω0\Omega = \omega_0 पर वेग बल के ठीक समान कला में है।
  4. (पट्ट-चौड़ाई) g(u)=2g(ur)g(u) = 2\,g(u_r) यथार्थ रूप से हल कीजिए, जहाँ ur=ω022λ2u_r = \omega_0^2 - 2\lambda^2, और उससे निकालिए कि R=Rmax/2R = R_{\max}/\sqrt2 वाली दोनों आवृत्तियाँ Ω±\Omega_\pm Ω+2Ω2=4λω02λ2\Omega_+^2 - \Omega_-^2 = 4\lambda\sqrt{\omega_0^2 - \lambda^2} को संतुष्ट करती हैं; निष्कर्ष निकालिए कि दुर्बल अवमंदन के लिए पट्ट-चौड़ाई Ω+Ω2λ\Omega_+ - \Omega_- \approx 2\lambda है, अर्थात् Qω0Ω+ΩQ \approx \frac{\omega_0}{\Omega_+ - \Omega_-}: तीक्ष्ण अनुनाद-शिखर उच्च-QQ निकायों के होते हैं।
  5. ω0=1\omega_0 = 1, λ=0.05\lambda = 0.05 (Q=10Q = 10), A=1A = 1 के लिए संख्यात्मक चित्र: Ωr\Omega_r, RmaxR_{\max}, स्थैतिक अनुक्रिया R(0)R(0) और सन्निकट पट्ट-चौड़ाई संगणित कीजिए।
  6. दिखाइए कि यदि 2λ2ω022\lambda^2 \geq \omega_0^2, तो RR (0,+)\intoo0{+\infty} पर कठोरतः ह्रासमान है: भारी अवमंदन वाले निकायों में अनुनाद-शिखर होता ही नहीं।

भाग IV — अवमंदन के बिना: विस्पंद और अनुनाद यहाँ λ=0\lambda = 0

  1. Ωω0\Omega \neq \omega_0 के लिए x+ω02x=Acos(Ωt)x'' + \omega_0^2 x = A\cos(\Omega t) का व्यापक हल खोजिए।
  2. कोशी समस्या x(0)=x(0)=0x(0) = x'(0) = 0 हल कीजिए और उत्तर को गुणनफल-रूप में बदलिए

    x(t)=2Aω02Ω2sin((ω0Ω)t2)sin((ω0+Ω)t2).x(t) = \frac{2A}{\omega_0^2 - \Omega^2}\, \sin\Bigl(\frac{(\omega_0 - \Omega)t}2\Bigr) \sin\Bigl(\frac{(\omega_0 + \Omega)t}2\Bigr) .
  3. ω0\omega_0 के निकट Ω\Omega के लिए इस गुणनफल को आवृत्ति ω0+Ω2\frac{\omega_0 + \Omega}2 के तीव्र दोलन के रूप में पढ़िए, जो आवृत्ति ω0Ω2\frac{\abs{\omega_0 - \Omega}}2 के मंद अन्वालोप से मॉडुलित है: यही विस्पंद हैं। अन्वालोप का आवर्तकाल और अधिकतम आयाम दीजिए, और ध्यान दीजिए कि Ωω0\Omega \to \omega_0 होने पर दोनों कैसे फट पड़ते हैं।
  4. tt नियत कीजिए और प्रश्न 19 के सूत्र में Ωω0\Omega \to \omega_0 लीजिए: दिखाइए कि सीमा है

    x(t)=Atsin(ω0t)2ω0,x_\infty(t) = \frac{A\,t\,\sin(\omega_0 t)}{2\omega_0} ,

    और सीधे जाँचिए कि xx_\infty x(0)=x(0)=0x(0) = x'(0) = 0 के साथ अनुनादी समीकरण x+ω02x=Acos(ω0t)x'' + \omega_0^2 x = A\cos(\omega_0 t) को हल करता है (उदाहरण 5.19 से तुलना कीजिए): अनुनाद निरंतर धीमे और निरंतर बड़े विस्पंदों की सीमा है।

  5. अनुनाद की दोनों नियतियों की तुलना कीजिए: अवमंदन के बिना रैखिक वृद्धि At2ω0\frac{At}{2\omega_0}, और दुर्बल अवमंदन के साथ RmaxQAω02R_{\max} \approx Q\,\frac{A}{\omega_0^2} पर संतृप्ति। एक वाक्य में: कौन-सी भौतिक क्रियाविधि पहली को दूसरी में बदल देती है?

भाग V — ऊर्जा-संतुलन और संश्लेषण।

  1. भाग II की स्थायी अवस्था में एक आवर्तकाल 2πΩ\frac{2\pi}\Omega पर (क) प्रणोदन द्वारा प्रविष्ट शक्ति Pin(t)=Acos(Ωt)xp(t)P_{\mathrm{in}}(t) = A\cos(\Omega t) \cdot x_p'(t) और (ख) अवमंदन द्वारा क्षयित शक्ति Pdiss(t)=2λxp(t)2P_{\mathrm{diss}}(t) = 2\lambda\,x_p'(t)^2 का औसत संगणित कीजिए। दिखाइए कि दोनों औसत λR2Ω2\lambda\,R^2\Omega^2 के बराबर हैं: प्रणोदन ठीक उतना ही देता है जितना अवमंदन जलाता है — इसीलिए स्थायी अवस्था स्थायी है।
  2. इस समस्या में ठीक कहाँ प्रयोग हुआ: (क) संरचना प्रमेय प्रमेय 5.12; (ख) सम्मिश्र चरघातांकी विधि; (ग) अध्याय 4 की शैली में वास्तविक चर वाले फलन का अध्ययन? प्रत्येक के लिए एक वाक्य।
  3. संश्लेषण: (EΩ)(E_\Omega) का पूरा व्यवहार-मानचित्र वर्णित कीजिए — मुक्त बनाम प्रणोदित, अवमंदित बनाम अनवमंदित, और शिखर की ऊँचाई, पट्ट-चौड़ाई तथा क्षणिक भाग के जीवनकाल को एक साथ नियंत्रित करने वाले एकमात्र विमाहीन घुंडी QQ की भूमिका — और बताइए कि कहानी आगे कहाँ चलती है: प्रथम कोटि के 2×22 \times 2 निकाय (अध्याय 21 और स्नातक वर्ष 2 का खंड) तथा किसी व्यापक आवर्ती प्रणोदन का ज्यावक्रीय फलनों में वियोजन (फूरिये श्रेणी, स्नातक वर्ष 3 के खंड में), जिसके लिए इस समस्या की ज्यावक्रीय स्थिति मूल निर्माण-खंड है।
हल

हल — समस्या 5.1.

1. 0<λ<ω00 < \lambda < \omega_0 के लिए χ(r)=r2+2λr+ω02\chi(r) = r^2 + 2\lambda r + \omega_0^2, Δ=4(λ2ω02)<0\Delta = 4(\lambda^2 - \omega_0^2) < 0: मूल λ±iωd-\lambda \pm \iu\omega_d, अतः प्रमेय 5.10 से

x(t)=eλt(λ1cosωdt+μ1sinωdt),(λ1,μ1)R2.x(t) = \eu^{-\lambda t}\bigl(\lambda_1\cos\omega_d t + \mu_1\sin\omega_d t\bigr), \qquad (\lambda_1, \mu_1) \in \R^2 .

λ=0\lambda = 0 के लिए: x=λ1cosω0t+μ1sinω0tx = \lambda_1\cos\omega_0 t + \mu_1\sin\omega_0 t। क्रांतिक (λ=ω0\lambda = \omega_0) और अति-अवमंदित (λ>ω0\lambda > \omega_0) अवस्थाएँ वही हैं जो अभ्यास 5.9 की हैं (समय के पुनर्मापन के बाद): (λ1+μ1t)eλt(\lambda_1 + \mu_1 t)\eu^{-\lambda t}, अथवा r±=λ±λ2ω02r_\pm = -\lambda \pm \sqrt{\lambda^2 - \omega_0^2} के साथ er±t\eu^{r_\pm t} के संयोजन।

2. अल्प-अवमंदित: xeλt(λ1+μ1)0\abs x \leq \eu^{-\lambda t}(\abs{\lambda_1} + \abs{\mu_1}) \to 0। क्रांतिक: (λ1+μ1t)eλt0(\lambda_1 + \mu_1 t)\eu^{-\lambda t} \to 0, क्योंकि चरघातांकी बहुपदों को हरा देते हैं (प्रतिज्ञप्ति 4.6)। अति-अवमंदित: r<r+=λ+λ2ω02<0r_- < r_+ = -\lambda + \sqrt{\lambda^2 - \omega_0^2} < 0, क्योंकि λ2ω02<λ\sqrt{\lambda^2 - \omega_0^2} < \lambda; दोनों चरघातांकी क्षय होते हैं।

3. x=2λxω02xx'' = -2\lambda x' - \omega_0^2 x का प्रयोग करते हुए (H)(H) के किसी हल पर:

E(t)=xx+ω02xx=x(2λxω02x)+ω02xx=2λx20.\mathcal E'(t) = x'x'' + \omega_0^2 x x' = x'\bigl(-2\lambda x' - \omega_0^2 x\bigr) + \omega_0^2 xx' = -2\lambda\,x'^2 \leq 0 .

यदि x(t0)=x(t0)=0x(t_0) = x'(t_0) = 0, तो E(t0)=0\mathcal E(t_0) = 0; E\mathcal E ऋणेतर और अवर्धमान है, अतः [t0,+)\intco{t_0}{+\infty} पर E0\mathcal E \equiv 0, जिससे वहाँ x0x \equiv 0 अनिवार्य हो जाता है; tt0t \leq t_0 के लिए वही तर्क x~(t)=x(2t0t)\tilde x(t) = x(2t_0 - t) पर चलाइए, जो अवमंदन λ-\lambda वाले समीकरण को हल करता है पर उसमें भी E~(t0)=0\tilde{\mathcal E}(t_0) = 0 और E~0\tilde{\mathcal E} \geq 0 के साथ E~=+2λx~20\tilde{\mathcal E}' = +2\lambda\tilde x'^2 \geq 0 है; ऋणेतर, अह्रासमान और (,t0]\intoc{-\infty}{t_0} के दाएँ सिरे पर शून्य होने का अर्थ है सर्वत्र शून्य। अतः R\R पर x0x \equiv 0 — अद्वितीयता की एक ऊर्जा-उपपत्ति, जो सभी λ0\lambda \geq 0 के लिए वैध है।

4. x(t+Td)=ReλteλTdcos(ωdt+2πφ)=eλTdx(t)x(t + T_d) = R\,\eu^{-\lambda t}\eu^{-\lambda T_d}\cos(\omega_d t + 2\pi - \varphi) = \eu^{-\lambda T_d}x(t)। प्रति छद्म-आवर्तकाल संकुचन-गुणक δ=λTd=2πλωd\delta = \lambda T_d = \frac{2\pi\lambda}{\omega_d} के साथ eδ\eu^{-\delta} है। ω0=1\omega_0 = 1, λ=0.1\lambda = 0.1 के लिए: ωd=0.99=0.99499\omega_d = \sqrt{0.99} = 0.99499, अतः δ=0.628320.99499=0.6315\delta = \frac{0.62832}{0.99499} = 0.6315: हर झूला अपने आयाम का e0.6353%\eu^{-0.63} \approx 53\% बनाए रखता है।

5. आयाम-गुणक eλt\eu^{-\lambda t} है, जो t=1λt = \frac1\lambda में e\eu से घटता है। उस अंतराल में 1/λTd=ωd2πλ\frac{1/\lambda}{T_d} = \frac{\omega_d}{2\pi\lambda} छद्म-आवर्तकाल आते हैं। λω0\lambda \ll \omega_0 के लिए ωdω0\omega_d \approx \omega_0, और यह ω02πλ=Qπ\approx \frac{\omega_0}{2\pi\lambda} = \frac Q\pi है। Q=300Q = 300 वाली गिटार की तार लगभग सौ आवर्तकाल तक गूँजती है; Q=1Q = 1 वाला दरवाज़े का अवमंदक एक भी पूरा नहीं करता।

6. बाएँ पक्ष में zeiΩtz\,\eu^{\iu\Omega t} प्रतिस्थापित करने पर z(Ω2+2iλΩ+ω02)eiΩtz\,(-\Omega^2 + 2\iu\lambda\Omega + \omega_0^2)\,\eu^{\iu\Omega t} मिलता है, जो ठीक z=Aω02Ω2+2iλΩz = \frac{A}{\omega_0^2 - \Omega^2 + 2\iu\lambda\Omega} के लिए AeiΩtA\,\eu^{\iu\Omega t} के बराबर है (हर अशून्य है: उसका काल्पनिक भाग 2λΩ>02\lambda\Omega > 0 है)। चूँकि गुणांक वास्तविक हैं, वास्तविक भाग xp=(zeiΩt)x_p = \Re\bigl(z\eu^{\iu\Omega t}\bigr) दाएँ पक्ष (AeiΩt)=AcosΩt\Re\bigl(A\eu^{\iu\Omega t}\bigr) = A\cos\Omega t वाले समीकरण को हल करता है।

7. D=(ω02Ω2)2+4λ2Ω2D = (\omega_0^2 - \Omega^2)^2 + 4\lambda^2\Omega^2 और φ(0,π)\varphi \in \intoo0\pi के साथ ω02Ω2+2iλΩ=Deiφ\omega_0^2 - \Omega^2 + 2\iu\lambda\Omega = \sqrt{D}\,\eu^{\iu\varphi} लिखिए (काल्पनिक भाग 2λΩ2\lambda\Omega धनात्मक है), अतः tanφ=2λΩω02Ω2\tan\varphi = \frac{2\lambda\Omega}{\omega_0^2 - \Omega^2}। तब z=ADeiφz = \frac{A}{\sqrt D}\eu^{-\iu\varphi} और

xp(t)=(ADei(Ωtφ))=Rcos(Ωtφ),R=AD.x_p(t) = \Re\Bigl(\frac A{\sqrt D}\,\eu^{\iu(\Omega t - \varphi)}\Bigr) = R\cos(\Omega t - \varphi), \qquad R = \frac A{\sqrt D} .

8. जब Ω0+\Omega \to 0^+: Dω04D \to \omega_0^4, अतः RA/ω02R \to A/\omega_0^2 और φ(0,π2)\varphi \in \intoo0{\frac\pi2} के साथ tanφ0+\tan\varphi \to 0^+: φ0\varphi \to 0। द्रव्यमान बल का अनुसरण अर्ध-स्थैतिक रूप से करता है, और उसका विस्थापन बल/दृढ़ता होता है। जब Ω+\Omega \to +\infty: DΩ4D \sim \Omega^4, अतः RA/Ω20R \sim A/\Omega^2 \to 0, और φπ\varphi \to \pi (सम्मिश्र संख्या ω02Ω2+2iλΩ\omega_0^2 - \Omega^2 + 2\iu\lambda\Omega कोणांक π\to \pi के साथ द्वितीय चतुर्थांश में चली जाती है): द्रव्यमान मुश्किल से हिलता है, और कला-विरोध में — जड़त्व प्रभावी है।

9. प्रमेय 5.12 (1) से प्रत्येक हल x=xp+xhx = x_p + x_h है, जहाँ xhx_h (H)(H) को हल करता है; प्रश्न 2 से xh(t)0x_h(t) \to 0, अतः x(t)xp(t)0x(t) - x_p(t) \to 0: सभी हल उसी स्थायी अवस्था की ओर अभिसरित होते हैं। आरंभिक शर्तें केवल क्षणिक भाग को आकार देती हैं।

10. यदि xx आवर्ती हल है, तो xxp=xhx - x_p = x_h (H)(H) का ऐसा आवर्ती हल है जो ++\infty पर 00 की ओर जाता है; सीमा 00 वाला आवर्ती फलन तत्समक रूप से 00 होता है (एक आवर्तकाल के उसके मान सदा दोहराते हैं, अतः हर मान 00 की ओर जाते किसी अनुक्रम की सीमा है)। अतः x=xpx = x_p

11. यहाँ λ=1\lambda = 1, ω02=2\omega_0^2 = 2, Ω=1\Omega = 1, A=1A = 1: z=121+2i=12i5z = \frac1{2 - 1 + 2\iu} = \frac{1 - 2\iu}5, अतः

xp=((12i)(cost+isint)5)=cost+2sint5.x_p = \Re\Bigl(\frac{(1 - 2\iu)(\cos t + \iu\sin t)}5\Bigr) = \frac{\cos t + 2\sin t}5 .

समांगी: r2+2r+2r^2 + 2r + 2 के मूल 1±i-1 \pm \iu हैं: xh=et(Ccost+Ssint)x_h = \eu^{-t}(C\cos t + S\sin t)। शर्तें: x(0)=15+C=0x(0) = \frac15 + C = 0 से C=15C = -\frac15 मिलता है; अवकलन करने पर x(0)=25C+S=0x'(0) = \frac25 - C + S = 0 से S=C25=35S = C - \frac25 = -\frac35 मिलता है। अतः

x(t)=cost+2sint5स्थायीetcost+3sint5क्षणिक,x(t) = \underbrace{\frac{\cos t + 2\sin t}5}_{\text{स्थायी}} - \underbrace{\eu^{-t}\,\frac{\cos t + 3\sin t}5}_{\text{क्षणिक}} ,

जहाँ क्षणिक भाग et\eu^{-t} की तरह मरता है।

12. प्रसार करने पर ur=ω022λ2u_r = \omega_0^2 - 2\lambda^2 के साथ g(u)=u22(ω022λ2)u+ω04=(uur)2+g(ur)g(u) = u^2 - 2(\omega_0^2 - 2\lambda^2)u + \omega_0^4 = (u - u_r)^2 + g(u_r), और

g(ur)=ω04ur2=(ω02ur)(ω02+ur)=2λ2(2ω022λ2)=4λ2(ω02λ2).g(u_r) = \omega_0^4 - u_r^2 = (\omega_0^2 - u_r)(\omega_0^2 + u_r) = 2\lambda^2\bigl(2\omega_0^2 - 2\lambda^2\bigr) = 4\lambda^2(\omega_0^2 - \lambda^2) .

यदि 2λ2<ω022\lambda^2 < \omega_0^2, तो ur>0u_r > 0 स्वीकार्य वर्ग-आवृत्ति है: gg का वहाँ कठोर निम्नतम है, अतः R=A/gR = A/\sqrt g का Ωr=ur=ω022λ2\Omega_r = \sqrt{u_r} = \sqrt{\omega_0^2 - 2\lambda^2} पर कठोर अधिकतम है, जहाँ Rmax=A/g(ur)=A2λω02λ2R_{\max} = A/\sqrt{g(u_r)} = \frac{A}{2\lambda\sqrt{\omega_0^2 - \lambda^2}}

13. R(0)=A/ω02R(0) = A/\omega_0^2, अतः

RmaxR(0)=ω022λω02λ2=ω02λω0ω02λ2=Q(1λ2ω02)1/2Q.\frac{R_{\max}}{R(0)} = \frac{\omega_0^2}{2\lambda\sqrt{\omega_0^2 - \lambda^2}} = \frac{\omega_0}{2\lambda}\cdot \frac{\omega_0}{\sqrt{\omega_0^2 - \lambda^2}} = Q\,\Bigl(1 - \frac{\lambda^2}{\omega_0^2}\Bigr)^{-1/2} \geq Q .

दुर्बल अवमंदन के लिए संशोधन-गुणक 11 के निकट है: अनुनाद स्थैतिक विस्थापन को मूलतः QQ से गुणा कर देता है।

14. u=Ω2u = \Omega^2 के साथ V(Ω)2=A2ug(u)V(\Omega)^2 = \frac{A^2 u}{g(u)}। उसके अवकलज का चिह्न g(u)ug(u)=(ω02u)2+4λ2uu(2(uω02)+4λ2)=(ω02u)2+2u(ω02u)=(ω02u)(ω02+u)g(u) - u\,g'(u) = (\omega_0^2 - u)^2 + 4\lambda^2 u - u\bigl(2(u - \omega_0^2) + 4\lambda^2\bigr) = (\omega_0^2 - u)^2 + 2u(\omega_0^2 - u) = (\omega_0^2 - u)(\omega_0^2 + u) का चिह्न है, जो u<ω02u < \omega_0^2 के लिए धनात्मक और उससे आगे ऋणात्मक है: अतः प्रत्येक अवमंदन के लिए कठोर अधिकतम ठीक Ω=ω0\Omega = \omega_0 पर है। वहाँ φ(0,π)\varphi \in \intoo0\pi के साथ tanφ\tan\varphi फट पड़ता है: φ=π2\varphi = \frac\pi2, और xp(t)=RΩsin(Ωtπ2)=RΩcos(Ωt)x_p'(t) = -R\Omega\sin(\Omega t - \frac\pi2) = R\Omega\cos(\Omega t) बल के ठीक समान कला में है: इष्टतम शक्ति-अंतरण।

15. R=Rmax/2    g(u)=2g(ur)    (uur)2=g(ur)=4λ2(ω02λ2)R = R_{\max}/\sqrt2 \iff g(u) = 2g(u_r) \iff (u - u_r)^2 = g(u_r) = 4\lambda^2(\omega_0^2 - \lambda^2), जिससे

u±=ur±2λω02λ2,Ω+2Ω2=4λω02λ2.u_\pm = u_r \pm 2\lambda\sqrt{\omega_0^2 - \lambda^2}, \qquad \Omega_+^2 - \Omega_-^2 = 4\lambda\sqrt{\omega_0^2 - \lambda^2} .

फिर Ω+Ω=Ω+2Ω2Ω++Ω\Omega_+ - \Omega_- = \frac{\Omega_+^2 - \Omega_-^2}{\Omega_+ + \Omega_-}, और λω0\lambda \ll \omega_0 के लिए दोनों Ω±ω0\Omega_\pm \approx \omega_0: Ω+Ω4λω02ω0=2λ\Omega_+ - \Omega_- \approx \frac{4\lambda\omega_0}{2\omega_0} = 2\lambda, अतः ω0Ω+Ωω02λ=Q\frac{\omega_0}{\Omega_+ - \Omega_-} \approx \frac{\omega_0}{2\lambda} = Q। अनुनाद-शिखर की चौड़ाई नापना उसका गुणता कारक नापना है।

16. Q=10Q = 10; Ωr=12(0.05)2=0.995=0.9975\Omega_r = \sqrt{1 - 2(0.05)^2} = \sqrt{0.995} = 0.9975; Rmax=12×0.0510.0025=10.1×0.99875=10.01R_{\max} = \frac1{2 \times 0.05 \sqrt{1 - 0.0025}} = \frac1{0.1 \times 0.99875} = 10.01; स्थैतिक अनुक्रिया R(0)=1R(0) = 1; पट्ट-चौड़ाई 2λ=0.1\approx 2\lambda = 0.111 ऊँचाई के पठार पर Q\approx Q ऊँचाई की एक ऊँची पतली नोक।

17. यदि 2λ2ω022\lambda^2 \geq \omega_0^2, तो सभी u>0u > 0 के लिए ur0u_r \leq 0 और g(u)=2(uur)>0g'(u) = 2(u - u_r) > 0: (0,+)\intoo0{+\infty} पर gg कठोरतः बढ़ता है, अतः R=A/gR = A/\sqrt g R(0)=A/ω02R(0) = A/\omega_0^2 से कठोरतः घटता है: अनुक्रिया शून्य आवृत्ति पर सबसे बड़ी है और कोई शिखर नहीं है।

18. γ=iΩ\gamma = \iu\Omega r2+ω02r^2 + \omega_0^2 का मूल नहीं है (क्योंकि Ωω0\Omega \neq \omega_0), अतः m=0m = 0 के साथ विधि 5.13 xp=AcosΩtω02Ω2x_p = \frac{A\cos\Omega t}{\omega_0^2 - \Omega^2} देता है (प्रतिस्थापित करके जाँचिए: कोज्या का Ω2+ω02-\Omega^2 + \omega_0^2 गुना)। व्यापक हल:

x(t)=AcosΩtω02Ω2+λ1cosω0t+μ1sinω0t.x(t) = \frac{A\cos\Omega t}{\omega_0^2 - \Omega^2} + \lambda_1\cos\omega_0 t + \mu_1\sin\omega_0 t .

19. x(0)=0x(0) = 0 λ1=Aω02Ω2\lambda_1 = -\frac A{\omega_0^2 - \Omega^2} पर बाध्य करता है, और x(0)=0x'(0) = 0 μ1=0\mu_1 = 0 पर:

x(t)=Aω02Ω2(cosΩtcosω0t)=2Aω02Ω2sin((ω0Ω)t2)sin((ω0+Ω)t2),x(t) = \frac{A}{\omega_0^2 - \Omega^2}\, \bigl(\cos\Omega t - \cos\omega_0 t\bigr) = \frac{2A}{\omega_0^2 - \Omega^2}\, \sin\Bigl(\frac{(\omega_0 - \Omega)t}2\Bigr) \sin\Bigl(\frac{(\omega_0 + \Omega)t}2\Bigr) ,

यह गुणनफल-सूत्र cosacosb=2sinba2sinb+a2\cos a - \cos b = 2\sin\frac{b - a}2\sin\frac{b + a}2 को a=Ωta = \Omega t, b=ω0tb = \omega_0 t के साथ लगाने से मिलता है।

20. ω0\omega_0 के निकट Ω\Omega के लिए दूसरी ज्या तीव्र आवृत्ति ω0+Ω2ω0\frac{\omega_0 + \Omega}2 \approx \omega_0 पर दोलन करती है, जबकि पहली आवृत्ति ω0Ω2\frac{\abs{\omega_0 - \Omega}}2 का मंद अन्वालोप है: तीव्र दोलन का आयाम अन्वालोप-आवर्तकाल 2πω0Ω\frac{2\pi}{\abs{\omega_0 - \Omega}} के साथ बढ़ता-घटता है (प्रति अन्वालोप-आवर्तकाल दो विस्पंद) और अधिकतम 2Aω02Ω2\frac{2A}{\abs{\omega_0^2 - \Omega^2}} तक पहुँचता है। जब Ωω0\Omega \to \omega_0, तब विस्पंद और धीमे (आवर्तकाल \to \infty) तथा और ऊँचे (आयाम \to \infty) दोनों हो जाते हैं।

21. tt नियत कीजिए। जब Ωω0\Omega \to \omega_0:

2Aω02Ω2sin((ω0Ω)t2)=2Aω0+Ωsin((ω0Ω)t2)ω0Ω2A2ω0t2=At2ω0,\frac{2A}{\omega_0^2 - \Omega^2}\sin\Bigl(\frac{(\omega_0 - \Omega)t}2\Bigr) = \frac{2A}{\omega_0 + \Omega}\cdot \frac{\sin\bigl(\frac{(\omega_0 - \Omega)t}2\bigr)} {\omega_0 - \Omega} \longrightarrow \frac{2A}{2\omega_0}\cdot\frac t2 = \frac{At}{2\omega_0} ,

जबकि sin((ω0+Ω)t2)sin(ω0t)\sin\bigl(\frac{(\omega_0 + \Omega)t}2\bigr) \to \sin(\omega_0 t): सीमा x(t)=Atsinω0t2ω0x_\infty(t) = \frac{At\sin\omega_0 t}{2\omega_0} है। सीधी जाँच: C=A2ω0C = \frac A{2\omega_0} के साथ x=Ctsinω0tx_\infty = Ct\sin\omega_0 t में x=2Cω0cosω0tCω02tsinω0tx_\infty'' = 2C\omega_0\cos\omega_0 t - C\omega_0^2 t \sin\omega_0 t है, अतः x+ω02x=2Cω0cosω0t=Acosω0tx_\infty'' + \omega_0^2 x_\infty = 2C\omega_0\cos\omega_0 t = A\cos\omega_0 t, जहाँ x(0)=0x_\infty(0) = 0 और x(0)=Csin0+Cω00cos0=0x_\infty'(0) = C\sin 0 + C\omega_0 \cdot 0 \cdot \cos 0 = 0 — और ω0=A=1\omega_0 = A = 1 के लिए यह ठीक उदाहरण 5.19 है। अनुनाद विस्पंदों का अपकर्ष है: अन्वालोप की पहली उभार, अनंत लंबाई तक खिंची हुई।

22. अवमंदन के बिना अनुनादी आयाम रैखिक रूप से और बिना परिबंध बढ़ता है; अवमंदन λ>0\lambda > 0 के साथ यह वृद्धि RmaxQAω02R_{\max} \approx Q\,\frac A{\omega_0^2} पर संतृप्त हो जाती है। क्रियाविधि: क्षय ऊर्जा को ऐसी दर से हटाता है जो आयाम के साथ बढ़ती है (प्रश्न 23), अतः संचय ठीक तब रुक जाता है जब अवमंदन उतनी ही तेज़ी से ऊर्जा जलाने लगता है जितनी तेज़ी से प्रणोदन उसे देता है।

23. xp=RΩsin(Ωtφ)x_p' = -R\Omega\sin(\Omega t - \varphi) के साथ, एक आवर्तकाल पर औसत cos2=sin2=12\langle\cos^2\rangle = \langle\sin^2\rangle = \frac12 और sincos=0\langle\sin\cos\rangle = 0 देते हैं:

Pdiss=2λR2Ω2sin2(Ωtφ)=λR2Ω2;\langle P_{\mathrm{diss}}\rangle = 2\lambda\,R^2\Omega^2\,\langle\sin^2(\Omega t - \varphi)\rangle = \lambda R^2\Omega^2 ;

और sin(Ωtφ)=sinΩtcosφcosΩtsinφ\sin(\Omega t - \varphi) = \sin\Omega t\cos\varphi - \cos\Omega t\sin\varphi का प्रसार करने पर:

Pin=ARΩcosΩtsin(Ωtφ)=ARΩsinφ2.\langle P_{\mathrm{in}}\rangle = -AR\Omega\,\bigl\langle\cos\Omega t\,\sin(\Omega t - \varphi)\bigr\rangle = AR\Omega\,\frac{\sin\varphi}2 .

चूँकि sinφ=2λΩD=2λΩRA\sin\varphi = \frac{2\lambda\Omega}{\sqrt D} = \frac{2\lambda\Omega R}A, यह ARΩ22λΩRA=λR2Ω2\frac{AR\Omega}2 \cdot \frac{2\lambda\Omega R}A = \lambda R^2\Omega^2 है: प्रविष्ट और क्षयित शक्तियाँ ठीक-ठीक संतुलित हैं — यही स्थायी अवस्था का परिभाषक गुणधर्म है।

24. (क) संरचना प्रमेय ने प्रत्येक हल को स्थायी अवस्था और क्षणिक भाग में बाँट दिया (प्रश्न 9–11) और अद्वितीयता को समांगी समस्या तक सीमित कर दिया। (ख) सम्मिश्र विधि ने विशिष्ट हल की खोज को सम्मिश्र संख्याओं के एक ही भाग में बदल दिया (प्रश्न 6), जहाँ आयाम और कला एक मापांक तथा एक कोणांक से पढ़ लिए जाते हैं। (ग) अनुनाद वक्र शुद्ध फलन-अध्ययन है — u=Ω2u = \Omega^2 में एक द्विघात, उसका निम्नतम, उसके स्तर-समुच्चय — अध्याय 4 की शैली में (प्रश्न 12–17)।

25. मुक्त और अवमंदित: क्षयशील छद्म-दोलन, जीवनकाल 1λ\frac1\lambda, लगभग Qπ\frac Q\pi झूले। प्रणोदित और अवमंदित: क्षणिक भाग मर जाते हैं और प्रणोदन-आवृत्ति पर एक अद्वितीय ज्यावक्रीय स्थायी अवस्था बची रहती है, जिसका आयाम ω0\omega_0 के पास शिखर पर होता है (ऊँचाई Q×\approx Q \times स्थैतिक, चौड़ाई ω0Q\approx \frac{\omega_0}Q) और कला 00 से ω0\omega_0 पर π2\frac\pi2 से होती हुई π\pi तक बहती है। मुक्त और अनवमंदित: शाश्वत दोलन। प्रणोदित और अनवमंदित: विस्पंद, जो यथार्थ अनुरणन पर रैखिक रूप से बढ़ते अनुनाद में अपकर्षित हो जाते हैं। एक ही विमाहीन संख्या Q=ω02λQ = \frac{\omega_0}{2\lambda} सब कुछ साध देती है — शिखर की ऊँचाई, पट्ट-चौड़ाई और क्षणिक जीवनकाल उसी घुंडी के तीन पाठ हैं। आगे की कड़ी: x+2λx+ω02xx'' + 2\lambda x' + \omega_0^2x को प्रथम कोटि के निकाय के रूप में फिर से लिखना अध्याय 21 तथा स्नातक वर्ष 2 के खंड की आव्यूह विधियाँ खोल देता है, और किसी भी आवर्ती प्रणोदन का ज्यावक्रीय फलनों में वियोजन (फूरिये श्रेणी, स्नातक वर्ष 3 का खंड) इस समस्या के एक-आवृत्ति विश्लेषण को सार्वभौमिक निर्माण-खंड बना देता है: प्रत्येक आवृत्ति के लिए हल कीजिए और अध्यारोपित कर दीजिए।