الفصل التاسع

من التفاضل والتكامل إلى الجبر

هل يمكننا تحويل التفاضل والتكامل إلى جبر؟

الفصل 9: من التفاضل والتكامل إلى الجبر ↗

المشتقات تتحول إلى ضرب

ميكا: أتذكر الورقة اللاصقة من الفصل 4؟ مع eate^{at} كان أخذ الميل مجرد × a\times\, a. راقب الآن ماذا يفعل المسبار بالميل.

كلمات جديدة

f′f' و f′′f''
اختصار لـ dfdt\frac{df}{dt}، ميل ff، ثم ميل ذلك الميل. في الفصل 7 كتبناهما بالنقاط، x˙\dot x و x¨\ddot x.
L{f}\mathcal{L}\{f\}
«تحويل ff»: الشيء نفسه مثل F(s)F(s)، مفيد حين تكون ff تعبيرًا كاملًا.
الحروف الكبيرة
الحرف الكبير هو تحويل حرفه الصغير: FF لـ ff، وHH للارتفاع hh، وXX لـ xx.

هذه هي الخاصية التي تجعل الأمر كله يستحق العناء. ماذا يحدث للتحويل عندما نأخذ مشتقة الإشارة؟ جرّبها مع f=eatf = e^{at}. مشتقتها هي aeata e^{at} (الفصل الرابع)، لذا

L{f′}=as−a=ss−a−1=s F(s)−f(0)\begin{aligned} \mathcal{L}\{f'\} &= \frac{a}{s-a} = \frac{s}{s-a} - 1 \\ &= s\,F(s) - f(0) \end{aligned}

تحولت المشتقة إلى ضرب في ss، ناقص القيمة الابتدائية f(0)f(0). التفاضل والتكامل أصبح جبرًا!

الفصل 9: من التفاضل والتكامل إلى الجبر ↗

ثيو: مهلاً. هذه حيلة تعمل مع eate^{at} لأن eate^{at} حالة خاصة. بعد قليل ستقول لي إنها قاعدة جاهزة علينا حفظها. لا أصدق أنها تعمل مع أي منحنى.

جون: إذن فلنجربها على أبشع منحنى يمكنك رسمه.

ارسم أي إشارة تريدها في المخطط الأول (أو اختر إعدادًا مسبقًا). نحسب ميلها، ثم نحسب طرفي القاعدة كليهما بآلة المساحة: تحويل الميل، و s F(s)−f(0)s\,F(s) - f(0). يُظهر المخطط الثاني هذين الطرفين لكل ss.

الفصل 9: من التفاضل والتكامل إلى الجبر ↗

ثيو: إنهما متطابقان، حتى مع خربشتي المتعرجة. إذن ليست قاعدة للحفظ في النهاية: إنها ما تفعله آلة المساحة بأي ميل.

إليك الجبر الناقص، على شكل صورة. يتغيّر حاصل الضرب لسببين. تتغيّر ff بينما يبقى المسبار ثابتًا، ويتغيّر المسبار بينما تبقى ff ثابتة. تخيّل مستطيلًا ضلعاه ff و e−ste^{-st}: يكبر على طول ضلع، ثم على طول الآخر. لهذه الخطوة اسم: التكامل بالتجزئة.

(f e−st)′=f′ e−st−s f e−st\big(f\,e^{-st}\big)' = f'\,e^{-st} - s\,f\,e^{-st}

اجمع تلك التغيّرات من 0 إلى ∞ (الفصل 3: المساحة تُلغي الميل) فتحصل على القيمة عند النهاية ناقص القيمة عند البداية. هذا ما يعنيه القوسان [ ⋅ ]0∞\big[\,\cdot\,\big]_0^\infty:

[f e−st]0∞=∫0∞f′ e−st dt−s F(s)\big[f\,e^{-st}\big]_0^\infty = \int_0^\infty f'\,e^{-st}\,dt - s\,F(s)

عند الطرف البعيد يكون المسبار قد سحق كل شيء إلى الصفر (ما دام ss يتغلب على نمو الإشارة)، فالطرف الأيسر هو ببساطة 0−f(0)0 - f(0). انقل s F(s)s\,F(s) إلى الطرف الآخر وها هي القاعدة.

ثيو: إذن f(0)f(0) ليست ترقيعًا. إنها المستطيل الذي بدأنا به.

طبّق القاعدة نفسها مرة أخرى، هذه المرة على f′f'، التي قيمتها الابتدائية f′(0)f'(0)، أي السرعة الابتدائية:

L{f′′}=s L{f′}−f′(0)=s(s F(s)−f(0))−f′(0)=s2F(s)−s f(0)−f′(0)\begin{aligned} \mathcal{L}\{f''\} &= s\,\mathcal{L}\{f'\} - f'(0) \\ &= s\big(s\,F(s) - f(0)\big) - f'(0) \\ &= s^2 F(s) - s\,f(0) - f'(0) \end{aligned}

جدولنا الخاص

تطبع الكتب المدرسية «جدول تحويلات لابلاس» وكأنه نزل من السماء. سنكسب جدولنا نحن، صفًّا تلو الآخر. كل صف هو مجرد آلة الفحص والمساحة مطبَّقة على إشارة واحدة، مع التحقق من أن الصيغة تتفق مع المساحة الفعلية.

الفصل 9: من التفاضل والتكامل إلى الجبر ↗

ميكا: صف الجيب هو المفضّل عندي. فالجيب ليس سوى دوّارَين توأمَين، وتحويل كل دوّار هو الصف الأسّي مع a=±iωa = \pm i\omega. في الفصل 6 جمعنا التوأمين لنحصل على جيب التمام. اطرحهما بدلًا من ذلك فتُلغى الأجزاء الأفقية، ويبقى 2isin⁡ωt2i\sin\omega t. من هنا تأتي الـ 2i2i! وضرب مقامَي التوأمين هو حاصل الضرب المرآتي من الفصل 7، (s−iω)(s+iω)=s2+ω2(s - i\omega)(s + i\omega) = s^2 + \omega^2.

(مع aa مركّب يظل الصف الأسّي صالحًا ما دام الجزء الأفقي من ss يتغلب على الجزء الأفقي من aa. بالنسبة للتوأمين هذا يعني ببساطة σ>0\sigma > 0.)

والصف الأخير لم يحتج إلى أي عمل جديد إطلاقًا: فالضرب في e−σte^{-\sigma t} ينقل ببساطة ss إلى s+σs+\sigma. أربعة صفوف، كلها اشتُقّت وكلها خضعت للتحقق.

قاعدة أخرى كنا نستخدمها بهدوء. آلة المساحة عادلة: اجعل الإشارة أطول 3 مرات فتصبح كل شريحة من مساحتها أطول 3 مرات، وكذلك المساحة كلها. اجمع إشارتين فتُجمع مساحتاهما. يسمّي الرياضيون هذا الخطية.

الفصل 9: من التفاضل والتكامل إلى الجبر ↗

جون: هذا ما يتيح لنا تحويل معادلة حدًّا بعد حد، مع الأعداد التي أمامها وكل شيء. ومعادلة الطائرة المسيّرة هي بالضبط مجموع من هذا النوع. هيا نستخدم الجدول!

عكس توحيد المقامات

ثيو: كل إجابة سنحصل عليها تخرج كسرًا واحدًا كبيرًا. لكن جدولنا لا يعرف إلا الكسور الصغيرة: 1s\frac1s، و1s−a\frac{1}{s-a}، والجيوب. كيف نكسر كسرًا كبيرًا إلى أجزاء من الجدول؟

في المدرسة كنت تجمع الكسور بإيجاد مقام مشترك. هنا نسير بالعكس: نبدأ من كسر واحد ونبحث عن الأجزاء الصغيرة التي مجموعها يساويه. هذا يُسمّى التفكيك إلى كسور جزئية:

1s+1s+2=2s+2s (s+2)إذن2s+2s (s+2)=As+Bs+2\frac{1}{s} + \frac{1}{s+2} = \frac{2s + 2}{s\,(s+2)} \qquad\text{إذن}\qquad \frac{2s+2}{s\,(s+2)} = \frac{A}{s} + \frac{B}{s+2}

ماذا نضع في بسط كل جزء؟ **مجهول واحد لكل قوة من ss في مقامه.** المقام مثل ss أو s+2s + 2 يأخذ عددًا عاديًا. أما المقام التربيعي مثل s2+2s+5s^2 + 2s + 5 فيأخذ B s+CB\,s + C: مجهولين، لأن في التربيعي متسعًا لاثنين.

الفصل 9: من التفاضل والتكامل إلى الجبر ↗

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

لإيجاد الأعداد طريق مختصر لأجزاء القوة الأولى: حيلة التغطية. اضرب في مقام الجزء، ثم اجعل ss القيمة التي تجعله صفرًا: تختفي كل الأجزاء الأخرى ولا يبقى إلا بسطه.

A=[s⋅2s+2s (s+2)]s=0=22=1,B=[(s+2)⋅2s+2s (s+2)]s=−2=−2−2=1A = \Big[s \cdot \frac{2s+2}{s\,(s+2)}\Big]_{s=0} = \frac{2}{2} = 1, \qquad B = \Big[(s+2) \cdot \frac{2s+2}{s\,(s+2)}\Big]_{s=-2} = \frac{-2}{-2} = 1

المقام التربيعي يحتاج إلى خطوة أخرى قبل أن يستطيع الجدول قراءته: ادفعه إلى الشكل (s+σ)2+ω2(s + \sigma)^2 + \omega^2 عن طريق إكمال المربع، بتنصيف العدد الأوسط وتربيعه. إذن s2+2s+5=(s+1)2+4s^2 + 2s + 5 = (s + 1)^2 + 4: جيب تمام وجيب متضائلان مع σ=1\sigma = 1 و ω=2\omega = 2.

حلّ الطائرة المسيّرة بالجبر

هذه هي طائرة الفصل 2 المسيّرة تحت تحكم P، هذه المرة بلا دفع تحويم مجاني، تبدأ من حافة ارتفاعها متر واحد. عند t=0t=0 نطلقها ونطلب منها الوصول إلى مترين. في الحساب من تلك اللحظة فصاعدًا، تتحول القيمة المرجعية الجديدة وقوة الجاذبية الثابتة كلتاهما مثل الخطوة:

m h′′+c h′=Kp (2−h)−mgm\,\out{h}'' + c\,\out{h}' = \eff{K_p}\,(\sp{2} - \out{h}) - \dis{m g}

الخطوة 1: حوّل كل حدّ. لأن المساحات تتضاعف وتُجمع (قسم الجدول)، نستطيع تحويل كل حدّ وحده. الخطوة بمقدار rr تصبح r/sr/s، والميول تصبح ضربًا في ss. إليك ما كتبناه:

m s2H+c s H=Kp(2s−H)−mgsm\,s^2\out{H} + c\,s\,\out{H} = \eff{K_p}\Big(\frac{\sp{2}}{s} - \out{H}\Big) - \frac{\dis{m g}}{s}

الخطوة 2: صار الأمر جبرًا الآن. اجمع H\out{H} في طرف واحد، ثم اقسم:

(ms2+cs+Kp) H=Kp 2−mgsH(s)=N(s)s (ms2+cs+Kp)N(s)=Kp 2−mg\begin{gathered} \begin{aligned} &(m s^2 + c s + \eff{K_p})\,\out{H} \\ &\quad = \frac{\eff{K_p}\,\sp{2} - \dis{m g}}{s} \end{aligned} \\[8pt] \out{H}(s) = \frac{N(s)}{s\,\big(m s^2 + c s + \eff{K_p}\big)} \\[8pt] N(s) = \eff{K_p}\,\sp{2} - \dis{m g} \end{gathered}

ميكا: انظروا إلى هذا. لم تبقَ أي ميول، مجرد حروف يُعاد ترتيبها. ثيو، دفترك توقف عن الدخان.

الخطوة 3: قسّمها إلى أجزاء من الجدول، كما في القسم السابق: عدد واحد في بسط الجزء ss، و Bs+CB s + C في بسط الجزء التربيعي:

H(s)=As+Bs+Cms2+cs+Kp\out{H}(s) = \frac{A}{s} + \frac{B s + C}{m s^2 + c s + \eff{K_p}}

كيف نجد AA و BB و CC؟ اضرب الطرفين في المقام المشترك s (ms2+cs+Kp)s\,(m s^2 + c s + \eff{K_p}):

N(s)=A (ms2+cs+Kp)+s (Bs+C)\begin{aligned} N(s) &= A\,(ms^2+cs+\eff{K_p}) \\ &\quad + s\,(Bs+C) \end{aligned}

ثم قارن الطرفين، قوةً واحدة من ss في كل مرة. الأعداد الثابتة تعطي AA، وحدود s2s^2 تعطي BB، وحدود ss تعطي CC:

A=2−mgKpB=−m A,C=−c A\begin{gathered} A = \sp{2} - \frac{\dis{mg}}{\eff{K_p}} \\ B = -m\,A, \qquad C = -c\,A \end{gathered}

لا تظهر فيزياء جديدة هنا: لقد فككنا ببساطة مقامًا مشتركًا لنصل إلى أجزاء يستطيع جدولنا إعادتها إلى الزمن. لاحظ أن AA لا يهتم بالمكان الذي بدأت منه الطائرة المسيّرة.

ثيو: ولم نضطر إلى تخمين este^{st} ولو مرة واحدة. الجدول خمّن عنا.

الخطوة 4: ابحث عن كل جزء في الجدول الذي اشتققناه للتو: A/sA/s خطوة، والجزء الثاني جيب تمام متضائل زائد جيب متضائل (الاستطراد أدناه يُظهر كيف يطابق هذين الصفين). اقرأ h(t)h(t).

جون: لنتحقق منها بالمحاكاة!

توقّع أولًا! نبدأ الطائرة المسيّرة على حافة ارتفاعها متر واحد ونشغّل صيغتنا بجانب المحاكاة. ماذا سنرى؟

الفصل 9: من التفاضل والتكامل إلى الجبر ↗

ميكا: أوه. الصيغة تبدأ من صفر، لكن الطائرة المسيّرة تبدأ على الحافة عند 1 متر. هل الجبر يكذب؟

ثيو: الجبر سليم. نحن كذبنا. انظروا إلى الخطوة 1: كتبنا L{h′′}=s2H\mathcal{L}\{h''\} = s^2 H و L{h′}=s H\mathcal{L}\{h'\} = s\,H، وأسقطنا كل h(0)h(0). هذا هو بالضبط الجزء الذي يتذكر الحافة.

اقلب المفتاح لإعادة حدود الارتفاع الابتدائي، وراقب الخطوتين 1 و2 والمعاملات أعلاه تتغيّر معه. الصيغة المتقطعة والمحاكاة الزرقاء تنطبقان الآن إحداهما على الأخرى تمامًا: أكبر فرق أقل من جزء من مليون من المتر، مجرد تقريب. غيّر Kp\eff{K_p} أو ارتفاع الحافة وستظلان متطابقتين.

الفصل 9: من التفاضل والتكامل إلى الجبر ↗

ثيو: في الفصل 2 اضطررنا إلى تطيير الطائرة المسيّرة لنعرف انخفاضها المستمر. هنا سقط ببساطة من كسر.

ما نعرفه الآن

ما توصلنا إليه

  • تصبح الميول ضربًا: L{f′}=sF(s)−f(0)\mathcal{L}\{f'\} = sF(s) - f(0). الحد f(0)f(0) يتذكر أين بدأت الإشارة. انسَه فتبدأ إجابتك من المكان الخطأ.
  • اشتققنا جدولنا الخاص: 1/s1/s، و1/(s−a)1/(s-a)، وω/(s2+ω2)\omega/(s^2+\omega^2)، والنسخة المُزاحة للجيب المتضائل. المساحات تتضاعف وتُجمع، فنستطيع تحويل معادلة حدًّا بعد حد.
  • الكسور الكبيرة تتفكك إلى كسور جزئية: مجهول واحد لكل قوة من ss في كل مقام، نجده بحيلة التغطية. والمقام التربيعي يصبح جيب تمام وجيبًا متقلصين.
  • حوّل، وقم بالجبر، وقسّم إلى أجزاء من الجدول، واقرأ الإجابة. الجزء A/sA/s هو المكان الذي تستقر فيه الطائرة المسيّرة، مع انخفاضها المستمر. إنها تطابق المحاكاة تمامًا، ما دمت تحتفظ بالقيم الابتدائية.

جون: "كل مرة حللنا فيها شيئاً، كانت الإجابة (شيء) ÷ (كثيرة حدود في s)، وكان سلوك الطائرة المسيّرة يأتي من حيث يكون المقام صفراً. ماذا لو كان ذلك الكسر هو الطائرة المسيّرة؟"

الفصل 9: من التفاضل والتكامل إلى الجبر ↗