Mathematics · الكتاب 2 · Grades 10–12

رياضيات المرحلة الثانوية

رياضيات المرحلة الثانوية · Grades 10–12

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

تربط المعادلة التفاضلية دالة بمشتقاتها. وتعبّر الفيزياء والكيمياء والبيولوجيا والاقتصاد عن قوانينها بهذه الصورة: فمعدل تبرّد جسم، وتفكك نواة مشعّة، ونمو مجموعة سكانية، كلها عبارات عن yy'. ويحل هذا الفصل حلًا كاملًا المعادلات الخطية من المرتبة الأولى ذات المعاملات الثابتة.

26.1 المعادلة y=ayy' = ay

تعريف 26.1 (المعادلة التفاضلية)

المعادلة التفاضلية هي معادلة مجهولها دالة yy، وتتضمن yy ومشتقاتها. و الحل على فترة II هو دالة قابلة للاشتقاق تحقق المعادلة عند كل نقطة من II.

مبرهنة 26.2 (حلول y=ayy' = ay)

ليكن aRa \in \R. حلول المعادلة y=ayy' = ay على R\R هي بالضبط الدوال

y(x)=Ceax,CR.y(x) = C\,\eu^{ax}, \qquad C \in \R .

ومن أجل كل زوج (x0,y0)(x_0, y_0) يوجد حل وحيد حيث y(x0)=y0y(x_0) = y_0.

برهان. كل دالة y=Ceaxy = C\eu^{ax} تحقق y=Caeax=ayy' = Ca\,\eu^{ax} = ay. وبالعكس، لتكن yy أي حل ولنضع z(x)=y(x)eaxz(x) = y(x)\,\eu^{-ax}. عندئذٍ

z(x)=y(x)eaxay(x)eax=(y(x)ay(x))eax=0,z'(x) = y'(x)\,\eu^{-ax} - a\,y(x)\,\eu^{-ax} = \bigl(y'(x) - a y(x)\bigr)\eu^{-ax} = 0,

إذن zz ثابتة، ولتكن CC، ويكون y(x)=Ceaxy(x) = C\eu^{ax}. والشرط الابتدائي y(x0)=y0y(x_0) = y_0 يفرض C=y0eax0C = y_0 \eu^{-a x_0}، على نحو وحيد.

ملاحظة 26.3

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

تفرض المعادلة y' = y ميلًا عند كل نقطة من المستوي (القطع الرمادية). والحلول C x — وهنا C = 1 (بالأزرق) و C = 0.3 (بالأحمر) و C = -0.5 (بالبرتقالي) — هي بالضبط المنحنيات التي تتبع حقل الميول هذا.
تفرض المعادلة y=yy' = y ميلًا عند كل نقطة من المستوي (القطع الرمادية). والحلول CexC\eu^{x} — وهنا C=1C = 1 (بالأزرق) و C=0.3C = 0.3 (بالأحمر) و C=0.5C = -0.5 (بالبرتقالي) — هي بالضبط المنحنيات التي تتبع حقل الميول هذا.

مثال 26.4 (التفكك الإشعاعي)

تحقق كمية مشعّة العلاقة N=λNN' = -\lambda N حيث λ>0\lambda > 0، إذن N(t)=N0eλtN(t) = N_0\,\eu^{-\lambda t}؛ وانظر المثال 23.9 من أجل نصف العمر.

26.2 المعادلة y=ay+by' = ay + b

مبرهنة 26.5 (حلول y=ay+by' = ay + b)

ليكن a0a \neq 0 و bRb \in \R. حلول المعادلة y=ay+by' = ay + b على R\R هي بالضبط الدوال

y(x)=Ceaxba,CR.y(x) = C\,\eu^{ax} - \frac{b}{a}, \qquad C \in \R .

والدالة الثابتة yp=bay_p = -\frac ba هي الحل التوازني. ومن أجل كل (x0,y0)(x_0, y_0) يوجد حل وحيد حيث y(x0)=y0y(x_0) = y_0؛ وإذا كان a<0a < 0، آل كل حل إلى التوازن ba-\frac ba لما x+x \to +\infty.

برهان. المقدار الثابت yp=bay_p = -\frac{b}{a} يحقق yp=0=ayp+by_p' = 0 = a y_p + b. والآن تكون yy حلًا إذا وفقط إذا كان

(yyp)=y=ay+b=a(yyp)+ayp+b=0=a(yyp),(y - y_p)' = y' = ay + b = a(y - y_p) + \underbrace{a y_p + b}_{=\,0} = a (y - y_p),

أي إذا وفقط إذا كانت z=yypz = y - y_p حلًا للمعادلة z=azz' = az. وحسب المبرهنة 26.2، z=Ceaxz = C\eu^{ax}، ومنه الصيغة والوجود والوحدانية. وإذا كان a<0a < 0، فإن eax0\eu^{ax} \to 0 لما x+x \to +\infty، إذن y(x)bay(x) \to -\frac ba.

طريقة 26.6 (حل y=ay+by' = ay + b مع شرط ابتدائي)

  1. أوجد الحل التوازني yp=bay_p = -\frac{b}{a} (بحل y=0y' = 0).
  2. اكتب الحل العام y=Ceax+ypy = C\eu^{ax} + y_p.
  3. حدّد CC من الشرط الابتدائي.
  4. تحقق من السلوك بعيد المدى في مواجهة الحدس الفيزيائي (فهل يتقارب الحل إلى التوازن؟).

والاستراتيجية نفسها — حل خاص + الحل العام للمعادلة المتجانسة — تمتد إلى y=ay+f(x)y' = ay + f(x): انظر القضية 26.8.

مثال 26.7 (قانون نيوتن في التبرّد)

فنجان قهوة عند 8080\,^\circC يقف في غرفة عند 2020\,^\circC. وينص قانون نيوتن على أن درجة الحرارة TT تحقق T=k(T20)T' = -k(T - 20) من أجل k>0k > 0 ما، أي T=kT+20kT' = -kT + 20k. والتوازن هو 2020، و T(t)=20+60ektT(t) = 20 + 60\,\eu^{-kt}: فتبرد القهوة بسرعة أسية إلى درجة حرارة الغرفة.

مهما كانت درجة حرارتها الابتدائية، تتقارب كل حلول T' = -k(T - 20) تقاربًا أسيًا إلى التوازن T = 20.
مهما كانت درجة حرارتها الابتدائية، تتقارب كل حلول T=k(T20)T' = -k(T - 20) تقاربًا أسيًا إلى التوازن T=20T = 20.

26.3 المعادلة y=ay+f(x)y' = ay + f(x)

قضية 26.8 (بنية مجموعة الحلول)

لتكن ff متصلة على فترة II وليكن ypy_p حلًا خاصًا واحدًا للمعادلة

y=ay+f(x)y' = ay + f(x)

على II. عندئذٍ تكون الحلول على II هي بالضبط الدوال y=Ceax+ypy = C\eu^{ax} + y_p، حيث CRC \in \R.

برهان. كما في المبرهنة 26.5: تكون yy حلًا إذا وفقط إذا كانت z=yypz = y - y_p تحقق z=(ay+f)(ayp+f)=azz' = (ay + f) - (ay_p + f) = az، أي إذا وفقط إذا كانت z=Ceaxz = C\eu^{ax}.

طريقة 26.9 (تخمين حل خاص)

ابحث عن حل خاص من شكل ff نفسه:

  • إذا كانت ff كثير حدود درجته nn: فجرّب كثير حدود درجته nn؛
  • f(x)=αekxf(x) = \alpha\,\eu^{kx} حيث kak \neq a: فجرّب yp=βekxy_p = \beta\,\eu^{kx}؛
  • f(x)=αeaxf(x) = \alpha\,\eu^{ax} (حالة الرنين): فجرّب yp=βxeaxy_p = \beta x\,\eu^{ax}.

ثم عوّض في المعادلة وطابق المعاملات.

مثال 26.10

حل y=2y+4xy' = 2y + 4x. جرّب yp=αx+βy_p = \alpha x + \beta: فإن α=2(αx+β)+4x\alpha = 2(\alpha x + \beta) + 4x من أجل كل xx يفرض 2α+4=02\alpha + 4 = 0 و α=2β\alpha = 2\beta، إذن α=2\alpha = -2 و β=1\beta = -1. والحل العام: y=Ce2x2x1y = C\eu^{2x} - 2x - 1.

26.4 تمارين

تمرين 26.1

حل على R\R: (أ) y=3yy' = 3y حيث y(0)=2y(0) = 2؛ (ب) 2y+y=02y' + y = 0 حيث y(0)=1y(0) = -1؛ (ج) y=y+5y' = -y + 5 حيث y(0)=0y(0) = 0.

حل

حل التمرين 26.1.

(أ) y=Ce3xy = C\eu^{3x}؛ و y(0)=2y(0) = 2 يعطي y=2e3xy = 2\eu^{3x}.

(ب) y=12yy' = -\frac12 y، إذن y=Cex/2y = C\eu^{-x/2}؛ و y(0)=1y(0) = -1 يعطي y=ex/2y = -\eu^{-x/2}.

(ج) التوازن yp=5y_p = 5؛ و y=Cex+5y = C\eu^{-x} + 5؛ و y(0)=0y(0) = 0 يعطي C=5C = -5: y=5(1ex)y = 5\left(1 - \eu^{-x}\right).

تمرين 26.2

تنمو مجموعة بكتيرية بمعدل يتناسب مع حجمها، فتتضاعف كل 33 ساعات. اكتب المعادلة التفاضلية التي يحققها التعداد N(t)N(t) وحدّد ثابت التناسب.

حل

حل التمرين 26.2.

N=kNN' = kN، إذن N(t)=N0ektN(t) = N_0 \eu^{kt}. والتضاعف في 33 ساعات يعني e3k=2\eu^{3k} = 2، أي

k=ln230.231 h1.k = \frac{\ln 2}{3} \approx 0.231\ \text{h}^{-1}.

تمرين 26.3

تحقق من أن yp(x)=xexy_p(x) = x\,\eu^{x} حل للمعادلة y=y+exy' = y + \eu^x، و أعطِ كل الحلول على R\R.

حل

حل التمرين 26.3.

yp(x)=ex+xex=yp(x)+exy_p'(x) = \eu^x + x\eu^x = y_p(x) + \eu^x: إذن ypy_p حل خاص (وهذه هي حالة الرنين في الطريقة 26.9). وحسب القضية 26.8، تكون الحلول y=Cex+xex=(C+x)exy = C\eu^{x} + x\eu^{x} = (C + x)\,\eu^x، حيث CRC \in \R.

تمرين 26.4 ★★

حل y=2y+exy' = -2y + \eu^{x} حيث y(0)=1y(0) = 1. (إرشاد: ابحث عن حل خاص من الشكل βex\beta\,\eu^{x}.)

حل

حل التمرين 26.4.

جرّب yp=βexy_p = \beta\eu^x: فإن βex=2βex+ex\beta\eu^x = -2\beta\eu^x + \eu^x يعطي 3β=13\beta = 1، إذن yp=13exy_p = \frac13\eu^x. والحل العام y=Ce2x+13exy = C\eu^{-2x} + \frac13\eu^{x}؛ والشرط y(0)=1y(0) = 1 يعطي C=23C = \frac23:

y(x)=23e2x+13ex.y(x) = \frac{2}{3}\,\eu^{-2x} + \frac{1}{3}\,\eu^{x}.

تمرين 26.5 ★★

يتفكك الكربون 14 بنصف عمر قدره 57305730 سنة. وتحوي عينة أثرية 60%60\% من كربون 14 الذي في كائن حي. قدّر عمرها.

حل

حل التمرين 26.5.

N(t)=N0eλtN(t) = N_0\,\eu^{-\lambda t} مع λ=ln25730\lambda = \frac{\ln 2}{5730} (المثال 23.9). ونحل eλt=0.6\eu^{-\lambda t} = 0.6:

t=ln(1/0.6)λ=5730ln(5/3)ln25730×0.51080.69314220 سنة.t = \frac{\ln(1/0.6)}{\lambda} = 5730\,\frac{\ln(5/3)}{\ln 2} \approx 5730 \times \frac{0.5108}{0.6931} \approx 4220 \text{ سنة}.

تمرين 26.6 ★★

خزان يحوي 100100 L من الماء النقي. ويدخله محلول ملحي فيه 0.20.2 kg من الملح في كل لتر بمعدل 55 L/min، ويخرج المزيج (المحفوظ متجانسًا) بالمعدل نفسه. وليكن m(t)m(t) كتلة الملح في الخزان عند اللحظة tt (بالدقائق).

  1. برّر أن m=1m20m' = 1 - \dfrac{m}{20}.
  2. حل، وحدّد نهاية m(t)m(t) لما t+t \to +\infty. وفسّر.
حل

حل التمرين 26.6.

1. يدخل الملح بمعدل 0.2×5=10.2 \times 5 = 1 kg/min. ويحمل التدفق الخارج تركيزًا m100\frac{m}{100} kg/L بمعدل 55 L/min، أي m20\frac{m}{20} kg/min. ومنه m=1m20m' = 1 - \frac{m}{20}.

2. التوازن mp=20m_p = 20؛ و m(t)=20+Cet/20m(t) = 20 + C\eu^{-t/20}، و m(0)=0m(0) = 0 يعطي C=20C = -20:

m(t)=20(1et/20)t+20 kg.m(t) = 20\left(1 - \eu^{-t/20}\right) \xrightarrow[t\to+\infty]{} 20 \text{ kg}.

وعلى المدى الطويل يساوي تركيز الخزان تركيز المحلول الداخل: 0.20.2 kg/L ×\times 100100 L =20= 20 kg.

تمرين 26.7 ★★

قافزة بالمظلة كتلتها 8080 kg تسقط تحت تأثير الثقالة (g=9.8 m/s2g = 9.8\ \text{m/s}^2) ومقاومة الهواء المتناسبة مع السرعة، بحيث تحقق سرعتها v=gkmvv' = g - \frac{k}{m}v حيث k=16k = 16 kg/s.

  1. حل المعادلة حيث v(0)=0v(0) = 0.
  2. احسب السرعة الحدية limt+v(t)\lim_{t\to+\infty} v(t)، ثم الزمن اللازم لبلوغ 95%95\% منها.
حل

حل التمرين 26.7.

1. v=gkmvv' = g - \frac{k}{m} v مع km=1680=0.2\frac km = \frac{16}{80} = 0.2. والتوازن v=mgk=9.80.2=49v_\infty = \frac{mg}{k} = \frac{9.8}{0.2} = 49 m/s؛ و v(t)=49+Ce0.2tv(t) = 49 + C\eu^{-0.2t}، و v(0)=0v(0) = 0 يعطي

v(t)=49(1e0.2t).v(t) = 49\left(1 - \eu^{-0.2 t}\right).

2. السرعة الحدية 4949 m/s (176\approx 176 km/h). ونريد 1e0.2t=0.951 - \eu^{-0.2t} = 0.95، أي e0.2t=0.05\eu^{-0.2t} = 0.05:

t=ln200.23.000.215 s.t = \frac{\ln 20}{0.2} \approx \frac{3.00}{0.2} \approx 15 \text{ s}.

تمرين 26.8 ★★★

(المعادلة اللوجستية.) تحقق مجموعة سكانية y(t)(0,1)y(t) \in \intoo{0}{1} (منسوبةً إلى التعداد الأقصى) العلاقة

y=y(1y).y' = y(1 - y).
  1. ليكن z=1yz = \dfrac{1}{y}. بيّن أن zz تحقق المعادلة الخطية z=z+1z' = -z + 1.
  2. حل بالنسبة إلى zz، ثم بالنسبة إلى yy، حيث y(0)=110y(0) = \frac{1}{10}.
  3. بيّن أن y(t)1y(t) \to 1 لما t+t \to +\infty وارسم تقريبيًا شكل منحنى الحل.
حل

حل التمرين 26.8.

1. z=1yz = \frac1y يعطي z=yy2=y(1y)y2=1yy=1y+1=z+1z' = -\frac{y'}{y^2} = -\frac{y(1-y)}{y^2} = -\frac{1-y}{y} = -\frac1y + 1 = -z + 1.

2. التوازن zp=1z_p = 1، إذن z(t)=1+Cetz(t) = 1 + C\eu^{-t}. ومن y(0)=110y(0) = \frac{1}{10}، يكون z(0)=10z(0) = 10، إذن C=9C = 9 و

y(t)=11+9et.y(t) = \frac{1}{1 + 9\,\eu^{-t}} .

3. لما t+t \to +\infty، يكون 9et09\eu^{-t} \to 0 و y(t)1y(t) \to 1. والمنحنى هو المنحنى اللوجستي الكلاسيكي على شكل S (السيني): نمو بطيء في البداية (yy صغير و yyy' \approx y)، وأسرع نمو عند y=12y = \frac12 (حيث يكون y=y(1y)y' = y(1-y) أعظميًا)، ثم إشباع نحو سعة الاحتمال 11.

تمرين 26.9 ★★★

لتكن yy حلًا للمعادلة y=ay+by' = ay + b وليكن un=y(n)u_n = y(n) من أجل nNn \in \N. بيّن أن (un)(u_n) تحقق علاقة تراجعية حسابية هندسية un+1=qun+ru_{n+1} = q u_n + r، وعبّر عن qq و rr بدلالة aa و bb. وما الذي يوافقه الشرط q<1\abs{q} < 1 بالنسبة إلى المعادلة التفاضلية؟

حل

حل التمرين 26.9.

حسب المبرهنة 26.5، y(x)=Ceaxbay(x) = C\eu^{ax} - \frac ba. ومنه

un+1=Cea(n+1)ba=ea(Ceanba)+ba(ea1)=qun+ru_{n+1} = C\eu^{a(n+1)} - \frac ba = \eu^{a}\left(C\eu^{an} - \frac ba\right) + \frac ba\left(\eu^a - 1\right) = q\,u_n + r

مع q=eaq = \eu^{a} و r=ba(ea1)r = \frac{b}{a}\left(\eu^{a} - 1\right). والشرط q<1\abs q < 1 يعني ea<1\eu^a < 1، أي a<0a < 0: وهو بالضبط الشرط الذي تتقارب تحته حلول المعادلة التفاضلية إلى التوازن ba-\frac ba — وفعلًا فإن النقطة الثابتة للعلاقة التراجعية هي r1q=ba\frac{r}{1 - q} = -\frac{b}{a}.

26.5 مسألة: الساعة التي داخل الأشياء

مسألة 26.1

مسألة نهاية الأسبوع — الكربون 14 يؤرّخ الكهوف، وقانون نيوتن في التبرّد يوقّت جريمة، ومعادلة صغيرة واحدة ترتدي أربعة أثواب

المعادلة التفاضلية قانون تغيّر؛ وحلّها يحوّل القانون إلى ساعة. والمعادلة الصغيرة نفسها y=ay+by' = ay + b (المبرهنة 26.5) تدقّ داخل الفحم ما قبل التاريخي، والقهوة المتبرّدة، وقطرات المطر الساقطة، و المحاليل الوريدية في المستشفى — وقراءة تلك الساعات هي شغل هذه المسألة: فهي تؤرّخ رسوم لاسكو، وتحدّد ساعة وفاة، وتدقّق فرضياتها كما يفعل عالِم جيد.

الجزء الأول — التمكّن.

  1. حل y=3yy' = 3y حيث y(0)=2y(0) = 2؛ ثم y=2y+6y' = -2y + 6 حيث y(0)=0y(0) = 0 (الطريقة 26.6).
  2. أعد استخراج حيلة الوحدانية وراء المبرهنة 26.2: إذا كان y=ayy' = ay، فاحسب مشتقة y(t)eaty(t)\,\eu^{-at} واستنتج أن كل حل هو CeatC\eu^{at}.
  3. من أجل y=2y+6y' = -2y + 6: أوجد الحل التوازني، ثم صِف مصير كل حل آخر لما t+t \to +\infty. (وهي النقط الثابتة لعلاقات المسألة 13.1 التراجعية، وقد صارت متصلة.)
  4. من أجل قانون التفكك y=yτy' = -\frac{y}{\tau}: بيّن أن نصف العمر هو t1/2=τln2t_{1/2} = \tau \ln 2.
  5. ارسم تقريبيًا (أو صِف) عائلة حلول y=2y+6y' = -2y + 6: فماذا تفعل الحلول التي تبدأ فوق 33؟ وتحت 33؟ وعند 33؟

الجزء الثاني — الكربون 14. تحافظ الأنسجة الحية على نسبة ثابتة من الكربون 14 المشعّ؛ وعند الموت يتوقف الاستهلاك ويتفكك المخزون: N=λNN' = -\lambda N، ونصف العمر 57305\,730 سنة.

  1. احسب λ\lambda (في السنة).
  2. عظمة تحتفظ بنسبة 20%20\,\% من نسبة الحي: فكم عمرها؟
  3. يحتفظ فحم من كهوف لاسكو المرسومة بنحو 15%15\,\%: فأرّخ الرسوم.
  4. رجل الجليد أوتزي، الذي وُجد في نهر جليدي في الألب، بلغ قياسه نحو 53%53\,\%: فأرّخه (ويقول علماء الآثار نحو 53005\,300 سنة — فكيف كان أداؤك؟).
  5. لماذا لا يستطيع الكربون 14 تأريخ الديناصورات؟ احسب الجزء الباقي بعد عشرة أنصاف أعمار، وصُغ الجزء بعد 6565 مليون سنة على صورة قوة للعدد 22، ثم اختم.
  6. حدود الخطأ: إذا لم تكن نسبة أوتزي 53%53\,\% معلومة إلا بدقة ±1%\pm 1\,\%، فاحسب مجال العمر. فأي دقة في المخبر تشتري أي دقة في التأريخ؟

الجزء الثالث — تبرّد نيوتن، وجريمة. جسم درجة حرارته TT في محيط درجة حرارته ثابتة TaT_a يبرد وفق T=k(TTa)T' = -k\,(T - T_a).

  1. حل المعادلة (بالتعويض z=TTaz = T - T_a):

    T(t)=Ta+(T0Ta)ekt.T(t) = T_a + (T_0 - T_a)\,\eu^{-kt} .
  2. قهوة تُصبّ عند 9090\,^\circC في غرفة عند 2020\,^\circC تقيس 7070\,^\circC بعد 55 دقائق. حدّد kk، ثم زمن الانتظار حتى تصير قابلة للشرب عند 5555\,^\circC.
  3. مسألة الحليب: لشرب أدفأ قهوة ممكنة خلال عشر دقائق، أينبغي وضع الحليب البارد الآن أم في اللحظة الأخيرة؟ أجب بالمعادلة (فماذا تفعل إضافة الحليب بالفجوة TTaT - T_a، وكيف تقود الفجوة الفقد؟).
  4. الطب الشرعي: يُعثر على جثة عند منتصف الليل عند 3030\,^\circC في غرفة عند 2020\,^\circC؛ وبعد ساعة تقيس 2828\,^\circC. وبافتراض 3737\,^\circC عند لحظة الوفاة، استعد kk من القياسين، ثم احسب ساعة الوفاة.
  5. دقّق ساعة الطبيب الشرعي: سمِّ ثلاثة انتهاكات واقعية لفرضيات النموذج، وسمِّ الجهة التي يحرف نحوها كل منها ساعة الوفاة المقدَّرة.

الجزء الرابع — أربعة أثواب ومنحنى S.

  1. المعادلة اللوجستية في التمرين 26.8، y=y(1y)y' = y(1 - y) حيث y(0)=0.1y(0) = 0.1: باستعمال تعويض التمرين، استخرج y(t)=11+9ety(t) = \dfrac{1}{1 + 9\eu^{-t}}، ثم أوجد لحظة نقطة الانعطاف y=12y = \frac12 — وهو تاريخ عالِم الأوبئة في المسألة 22.1، وقد صار قابلًا للحساب.
  2. قطرة مطر تطيع v=10v2v' = 10 - \frac v2 (الثقالة ناقص قوة السحب)، حيث v(0)=0v(0) = 0. أوجد السرعة الحدية، ثم حل لإيجاد v(t)v(t)، ثم احسب متى تبلغ القطرة 95%95\,\% من السرعة الحدية.
  3. محلول وريدي في مستشفى يسلّم دواءً بمعدل ثابت بينما يطرحه الجسم تناسبيًا: c=40.5cc' = 4 - 0.5c حيث c(0)=0c(0) = 0. أوجد التركيز المستقر والزمن اللازم لبلوغ نصفه. وصُغ في جملة واحدة التوازي مع علاقات القرض التراجعية في التمرين 26.9.
  4. الخاتمة — معادلة واحدة، وأربعة أثواب: التفكك، والتبرّد، والسقوط، والجرعات، وكلها y=ay+by' = ay + b بإشارات وأسماء مختلفة. اسرد حلقة النمذجة التي شغّلتها هذه المسألة أربع مرات (قانون \to معادلة \to حل \to معايرة \to تنبؤ \to تدقيق)، وسمِّ الظاهرة التي تحتاج إلى yy'' — وأي مسألة لقيتها من قبل.
حل

حل المسألة 26.1.

1. y=2e3ty = 2\eu^{3t}. ومن أجل الثانية: التوازن 33، إذن y=3+Ce2ty = 3 + C\eu^{-2t} مع y(0)=0y(0) = 0: أي C=3C = -3: y=3(1e2t)y = 3\left(1 - \eu^{-2t}\right).

2. (yeat)=yeatayeat=(yay)eat=0\left(y\eu^{-at}\right)' = y'\eu^{-at} - ay\eu^{-at} = (y' - ay)\eu^{-at} = 0: فالجداء مقدار ثابت CC، إذن y=Ceaty = C\eu^{at} — فلا حل يفلت.

3. التوازن: y3y \equiv 3. وكل حل آخر هو 3+Ce2t3 + C\eu^{-2t} حيث C0C \neq 0: فتموت الدالة الأسية وينزلق الحل إلى 33 — وهي نقطة ثابتة مستقرة، التوأم المتصل للنقط الثابتة التراجعية في المسألة 13.1.

4. y=y0et/τy = y_0\eu^{-t/\tau} ينتصف عندما يكون et/τ=12\eu^{-t/\tau} = \frac12: أي t=τln2t = \tau\ln 2.

5. كل الحلول انسحابات شاقولية للتفكك نحو 33: فالتي تبدأ فوقه تهبط إلى 33، والتي تبدأ تحته ترتفع إلى 33، والحل الثابت 33 يقف ساكنًا — إنه قِمع حول التوازن.

6. λ=ln257301.21×104\lambda = \frac{\ln 2}{5730} \approx 1.21 \times 10^{-4} في السنة.

7. eλt=0.2\eu^{-\lambda t} = 0.2: t=ln5λ13300t = \frac{\ln 5}{\lambda} \approx 13\,300 سنة.

8. t=ln(1/0.15)λ15700t = \frac{\ln(1/0.15)}{\lambda} \approx 15\,700 سنة: فثيران لاسكو من أواخر العصر الجليدي.

9. t=ln(1/0.53)λ5250t = \frac{\ln(1/0.53)}{\lambda} \approx 5\,250 سنة: أي ضمن عمر إنسان من قيمة علماء الآثار — فالساعة تعمل.

10. عشرة أنصاف أعمار تترك 2100.1%2^{-10} \approx 0.1\,\%: أي على حافة القياس. وبعد 6565 مليون سنة يكون الجزء 265000000/57302113442^{-65\,000\,000/5\,730} \approx 2^{-11\,344}: فلا تبقى ذرة من المخزون الأصلي في أي أحفورة — وتُؤرَّخ الديناصورات بساعات أبطأ (البوتاسيوم–الأرغون وأمثاله).

11. 52%52\,\% تعطي 5405\approx 5\,405 سنة، و 54%54\,\% تعطي 5095\approx 5\,095: أي إن ±1%\pm 1\,\% من الكيمياء تصير نحو ±155\pm 155 سنة من التاريخ — فالدقة في المخبر دقة في لوحة المتحف.

12. z=TTaz = T - T_a تحقق z=kzz' = -kz: z=(T0Ta)ektz = (T_0 - T_a)\eu^{-kt}، و T=Ta+zT = T_a + z.

13. الفجوة الابتدائية T0Ta=70T_0 - T_a = 70؛ وبعد خمس دقائق تكون الفجوة 7020=5070 - 20 = 50، إذن 50=70e5k50 = 70\eu^{-5k}: أي k=ln(70/50)50.067k = \frac{\ln(70/50)}{5} \approx 0.067 في الدقيقة. وتصير قابلة للشرب عند: 5520=35=70ekt55 - 20 = 35 = 70\eu^{-kt}: t=ln2k10.3t = \frac{\ln 2}{k} \approx 10.3 دقيقة.

14. معدل الفقد يتناسب مع الفجوة TTaT - T_a: فالقهوة الحارقة تنزف الحرارة أسرع. وإضافة الحليب الآن تقلّص الفجوة فورًا، فيُفقد حرّ أقل خلال العشر دقائق؛ أما إضافته في النهاية فتترك القهوة تبرد بأقصى سرعة أولًا. فمن أجل أدفأ فنجان: الحليب أولًا. (والمضيفون العجولون يقلبون الترموديناميك.)

15. الفجوات عن المحيط: عند منتصف الليل 1010، وبعد ساعة 88: إذن ek=0.8\eu^{-k} = 0.8: k=ln1080.223k = \ln\frac{10}{8} \approx 0.223 في الساعة. وعند الوفاة كانت الفجوة 1717؛ وقد تفككت إلى 1010 عند منتصف الليل: s=ln(17/10)k2.4s = \frac{\ln(17/10)}{k} \approx 2.4 ساعة. فساعة الوفاة: نحو 21:4021{:}40.

16. غرفة مدفَّأة أو ذات تيار هواء (TaT_a غير ثابتة) تحرف الساعة في أي من الجهتين؛ والملابس أو كتلة الجسم تغيّر kk (وجداول الطبيب الشرعي تصحّح لذلك) — و kk الخاطئ يضرب الفترة كلها بمعامل؛ وجثة نُقلت من مكان آخر تعيد ضبط TaT_a في منتصف التفكك، فتزيّف وفاة أبكر أو أمتأخر. فالمعادلة أمينة؛ والفرضيات هي التي تحمل الخطر.

17. z=1yz = \frac1y تطيع z=1zz' = 1 - z: z=1+9etz = 1 + 9\eu^{-t} (من z(0)=10z(0) = 10)، إذن y=11+9ety = \frac{1}{1 + 9\eu^{-t}}. والانعطاف عند y=12y = \frac12: 9et=19\eu^{-t} = 1: t=ln92.2t = \ln 9 \approx 2.2 — وهو تاريخ انقلاب الوباء، مباشرة من الصيغة.

18. السرعة الحدية: v=0v' = 0 عند v=20v = 20 m/s. والحل: v(t)=20(1et/2)v(t) = 20\left(1 - \eu^{-t/2}\right). ثم 0.950.95: et/2=0.05\eu^{-t/2} = 0.05: t=2ln206t = 2\ln 20 \approx 6 s — فقطرات المطر تبلغ سرعتها المستقرة في ثوانٍ، ولهذا لا يقتل المطر.

19. الحالة المستقرة: c=8c = 8. ونصفها: 4=8(1et/2)4 = 8\left(1 - \eu^{-t/2}\right): et/2=12\eu^{-t/2} = \frac12: t=2ln21.4t = 2\ln 2 \approx 1.4 ساعة. فالمحلول الوريدي هو التوأم المتصل للقرض: دخل ثابت، وخروج تناسبي — و التمرين 26.9 يجعل القاموس مضبوطًا.

20. التفكك (a<0a < 0 و b=0b = 0)، والتبرّد (y=TTay = T - T_a)، والسقوط (a<0a < 0 و b>0b > 0: اقتراب من التوازن من تحت)، والجرعات (المثل نفسه، في وريد): معادلة واحدة، وأربعة عوالم. والحلقة: صُغ قانون التغيّر؛ واكتب المعادلة؛ وحلّ (الطريقة 26.6)؛ وعاير الثوابت على القياسات؛ وتنبّأ؛ ثم دقّق الفرضيات (السؤال 16). وأما التذبذب فيحتاج إلى y=ω2yy'' = -\omega^2 y — وقد لُقي وحُلّ و أُطلق متأرجحًا في المسألة 24.1.

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

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