الفصل الثاني عشر

العودة إلى الدش

لماذا تذبذبنا، وما مدى قربنا من الحافة؟

الدش بعيون جديدة

جون: ركّبتُ PID الطائرة المسيّرة على الدش (الإحساس السريع نفسه، بعد تحويله إلى وحدات المقبض). فتأرجح: ساخن، بارد، ساخن… تمامًا مثل ميكا في الفصل 0. الطائرة المسيّرة لم تفعل ذلك قط!

ثيو: الطائرة المسيّرة ليس فيها تأخير كبير. أما في الدش، فيضيف الأنبوب تأخيرًا مقداره 2.5 ثانية. دعونا نفهم أخيرًا ما الذي يفعله التأخير بالضبط.

في الدش جزءان بطيئان: الأنبوب، وهو تأخير خالص مقداره 2.5 ثانية، والمُنعِّم، أي رأس الدش وجدران الأنبوب وهي تسخن بالثابت الزمني من الفصل 3، $\tau = 1$ ث (بعد ثانية واحدة يكون قد أغلق 63% من الفجوة نحو درجة حرارة جديدة).

أولاً، اعتراف بخصوص «اليد العادية» في الفصل 0. لقد كانت تدير المقبض بسرعة تتناسب مع مقدار خطأ درجة حرارة الماء: سرعة المقبض $= k\cdot\err{e}$ ، مع $k = 0.8\,\%$ من المقبض لكل ثانية لكل درجة مئوية. السرعة المتناسبة مع الخطأ تعني أن موضع المقبض هو تراكم أخطاء الماضي. هذا متحكم تكاملي، مع $K_i = k$ .

ميكا: لحظة، إذن في الفصل 0 كنتُ متحكم PID؟ حسنًا، متحكم تكاملي. متحكم مفرط الحماس.

التأخير يزحزح الموجة

ثيو: ما زلت أعتقد أن التأخير يجعل كل شيء أبطأ فقط. لا ينبغي أن يكون قادرًا على جعل أي شيء غير مستقر.

لنبدأ بعزل التأخير الزمني الخالص في الأنبوب. تخيّل موجة جيبية تدخل الأنبوب وتصل إلى رأس الدش من دون أن يتغيّر شكلها: موجة الخرج هي نفسها، لكنها تصل بعد 2.5 ثانية. ما مقدار هذا التأخر قياسًا إلى دورة موجية واحدة؟

الفصل 12: العودة إلى الدش ↗

ثيو: لحظة. مع تذبذبة مدتها 5 ثوانٍ يصل الماء الساخن بالضبط حين يطلب المقبض البارد. هذا ليس أبطأ. هذا معكوس.

يبقى تأخير الأنبوب 2.5 ثانية، لكن تأخر الطور الذي يسببه يعتمد على سرعة تذبذب الموجة. سمِّ التأخر $\phi$ والتأخير $L$ . إذا استغرقت التذبذبة الواحدة $T$ ثانية (هنا $T$ زمن، لا درجة حرارة)، فسرعتها الزاوية $\omega=2\pi/T$ راديان في الثانية. التأخر الموجب يعني أن الخرج يلحق بالدخل متأخرًا. يكتب المهندسون التأخر زاويةً سالبة، ولهذا تنحدر المخططات أدناه إلى الأسفل.

التذبذبة الكاملة هي دورة كاملة لدوّار الفصل 5: 360°، أو 2π راديان. إذن المهم هو أي حصة من التذبذبة يغطيها التأخير: $L/T$ .

\phi = \underbrace{\frac{L}{T}}_{\text{حصة من التذبذبة}}\cdot 360^\circ = \omega L \text{ راديان},\qquad \omega = \frac{2\pi}{T}

الفصل 12: العودة إلى الدش ↗

مع تذبذبة مدتها 5 ثوانٍ، تكون 2.5 ثانية نصف تذبذبة: 180°. في نموذج التأخير الخالص هذا، تقع كل قمة في الخرج مقابل قاع في الدخل. ومُنعِّم الدش يصغّر التذبذبات السريعة أيضًا؛ سنضيفه بعد قليل.

ميكا: إذن الـ 2.5 ثانية نفسها دفعة خفيفة لتذبذبة كسولة، وقلب كامل لتذبذبة سريعة. الأنبوب لا يتغير أبدًا. التذبذبة هي التي تتغير!

المسبار يرى التأخير

مسبار الفصل 8 يعطي العامل نفسه لأي إشارة، لا للدوّار وحده. الإشارة التي تصل متأخرة $L$ ثانية تبدأ بتكديس المساحة متأخرة $L$ ثانية، وحينها يكون المسبار قد تقلّص بالفعل بمقدار $e^{-sL}$ . وقبل وصولها تكون الإشارة المتأخرة صفرًا، لذا ضع $u = t - L$ :

\int_0^\infty f(t-L)\,e^{-st}\,dt = \int_0^\infty f(u)\,e^{-s(u+L)}\,du = e^{-sL}\,F(s)

إذن التأخير الخالص يضرب كل تحويل في $e^{-sL}$ . هذه هي دالة النقل الخاصة به، وعلى الخط $s = i\omega$ هي بالضبط الإدارة إلى الخلف بمقدار $\omega L$ التي رأيناها أعلاه.

ثيو: لكن $e^{-sL}$ ليست كسرًا من كثيرات حدود. أين أقطابها؟

جون: لا يمكنك حصرها في قائمة. ضعها في حلقة فيصبح لـ $1 + C\,P = 0$ عدد لا نهائي من الحلول. لهذا يتوقف هذا الفصل عن رسم خرائط الأقطاب ويقيس التذبذبات بدلًا منها.

قياس الدش تذبذبًا بعد آخر

الآن أضف المُنعِّم إلى تأخير الأنبوب. رأس الدش وجدران الأنبوب تحتاج وقتًا لتسخن، فيؤخّر الدش الكامل تغيّرات المقبض وينعّمها معًا. الخطة: حرّك المقبض بتذبذبة بسرعة واحدة، وانتظر حتى يثبت الإيقاع، ثم قِس رقمين:

الكسب (الكلمة نفسها التي نستخدمها لكسب المتحكم: الحجم الخارج ÷ الحجم الداخل): كم تبلغ تذبذبة درجة الحرارة، مقارنةً بتذبذبة خليط الماء الذي طلبناه.
الطور: كم تتأخر، على شكل زاوية: حصتها من تذبذبة كاملة، مضروبة في 360°.

كل قياس يضع نقطة على مخططين. قس خمس سرعات على الأقل، بطيئة وسريعة.

الفصل 12: العودة إلى الدش ↗

المسطرة المضغوطة

في هذه المخططات، كل ×10 خطوة متساوية على المحور: من 0.1 إلى 1 يأخذ المساحة نفسها التي يأخذها من 1 إلى 10. منتصف المسافة بين 0.1 و1 هو نحو 0.3، لا 0.55. المسطرة مضغوطة لكي تتسع تذبذبة مدتها 100 ثانية وتذبذبة مدتها 3 ثوانٍ في المخطط نفسه. انظر إلى كسب المُنعِّم على هذه المسطرة: مستوٍ للتذبذبات البطيئة، ثم ينزلق على خط شبه مستقيم للسريعة، مع انحناءة قرب $\omega = 1/\tau$ .

جون: الخطوات المتساوية أضعاف متساوية، لا إضافات متساوية. مثل الأوكتافات على البيانو!

للمسطرة المضغوطة هدية ثانية: عليها يصبح الضرب جمعًا. فالضرب ×10 خطوة واحدة للأعلى دائمًا، أيًّا كانت نقطة البداية، و×0.1 خطوة واحدة للأسفل. لهذا ضُغط محور الكسب أيضًا. حول الحلقة تُضرب المكاسب بعضها في بعض (اليد × المُنعِّم × الأنبوب)، لذا تُجمع ارتفاعاتها على هذه المسطرة ببساطة، تمامًا مثل تأخرات الطور. إنه لوغاريتم الفصل 4 وهو يعمل: $\ln(a\,b) = \ln a + \ln b$ .

المُنعِّم يصغّر التذبذبات السريعة، بينما يواصل الأنبوب إضافة تأخر الطور. منحنيا الكسب والطور هذان مقابل سرعة التذبذب يشكّلان مخطط بود. ومعًا هما الاستجابة الترددية للدش: كيف يعامل كل سرعة تذبذب. الدخل هنا هو درجة حرارة الماء المخلوط، لذا يكون كسب التذبذبات البطيئة 1. لو استخدمنا نسبة المقبض دخلًا، لكان الكسب أكبر 45 مرة. للحصول على الاستجابة الترددية من الصيغة، ضع $s = i\omega$ : الحيلة نفسها التي وجدت حافة الفصل 11. (وهذا ينجح لأن الدش وحده مستقر، فتستقر تذبذباته على إيقاع ثابت يمكننا قياسه.)

ماذا يفعل المُنعِّم بالدوّار

ثيو: حيلة الفصل 4 مرة أخرى: خمّن أن الجواب هو الدوّار نفسه، لكن بحجم وزاوية مختلفين.

\tau\,\dot y + y = u,\quad u = e^{i\omega t},\ y = Y\,e^{i\omega t} \;\Rightarrow\; (1 + i\omega\tau)\,Y = 1 \;\Rightarrow\; Y = \frac{1}{1+i\omega\tau}

ضع الدوّار $e^{i\omega t}$ في قاعدة المُنعِّم، فيجب أن يحقق السهم $Y$ أن $(1 + i\omega\tau)\,Y = 1$ . أسهم الفصل 7: السهم $1 + i\omega\tau$ طوله $\sqrt{1+\omega^2\tau^2}$ ويشير للأعلى بالزاوية $\arctan(\omega\tau)$ . القسمة عليه تصغّر بذلك الطول وتدير إلى الخلف بتلك الزاوية. التذبذبات البطيئة بالكاد تلاحظ. أما السريعة فتنضغط ويقترب تأخرها من ربع دورة، لكنه لا يتجاوزه أبدًا.

الفصل 12: العودة إلى الدش ↗

G(i\omega) = \frac{e^{-i\omega L}}{1 + i\omega\tau}, \qquad |G| = \frac{1}{\sqrt{1+\omega^2\tau^2}}, \qquad \angle G = -\arctan(\omega\tau) - \omega L

الفصل 12: العودة إلى الدش ↗

جون: وجدتُها! قرب 0.95 راديان/ثانية يقلب الدش وحده التذبذبة رأسًا على عقب. تذكّروا هذه السرعة.

ما مدى قربنا من الهاوية؟

الآن أغلق الحلقة: اليد ← المقبض ← الأنبوب ← الجلد ← اليد. تتبّع تذبذبة طوال الدورة. بالنسبة ليد الفصل 0، رحلة الدوران حول الحلقة هي $G_{\circ}$ (دائرة صغيرة ترمز إلى الحلقة، لأن $L$ هو التأخير بالفعل):

G_{\circ}(s) = \underbrace{\frac{k}{s}}_{\text{اليد}}\cdot\underbrace{\frac{45\,e^{-2.5s}}{s+1}}_{\text{الدش}}

كومة التذبذبة

ميكا: يدي تكوّم الخطأ. ماذا تفعل الكومة بالتذبذبة؟

\int \cos(\omega t)\,dt = \frac{\sin(\omega t)}{\omega} \qquad\Longleftrightarrow\qquad \frac{1}{i\omega}:\ \text{ربع دورة إلى الخلف، والحجم } \tfrac{1}{\omega}\text{ من الأصل}

المساحة تحت جيب التمام هي جيب، متأخر ربع تذبذبة: تأخر 90°، عند كل سرعة. والكومة أكبر للتذبذبات البطيئة، لأن كل حدبة تدوم أطول قبل أن تُلغى: يتناسب الحجم مع $1/\omega$ . قال الفصل 5 إن السرعة المتجهة هي $i\omega$ مضروبًا في الموضع. التكويم هو العكس: القسمة على $i\omega$ . هذا هو $k/s$ في الحلقة.

التغذية الراجعة تطرح (الخطأ = الهدف − درجة الحرارة). هذا ضرب في (−1)، نصف دورة الفصل 5، فتنقلب الموجة مرة. إذا كانت رحلة الدوران حول الحلقة تؤخرها أيضًا بمقدار 180° تمامًا، يلغي الانقلابان أحدهما الآخر، فيصل كل تصحيح متزامنًا تمامًا مع التأرجح، كمن يدفع أرجوحة الملعب في اللحظة المناسبة بالضبط. إذا أعادت الحلقة الموجة بحجمها الكامل (كسب 1)، يستمر التأرجح من تلقاء نفسه؛ وإذا كان أكبر من 1، يكبر. هذه هي $1 + C\,P = 0$ من الفصل 10 مجددًا، على الخط $s = i\omega$ : اللفّة التي تعود مقلوبة وبحجمها الكامل هي $G_{\circ}(i\omega) = -1$ .

جون: بالنسبة لحلقة الدش البسيطة هذه، ينخفض الطور عبر −180° عند سرعة تذبذب واحدة ذات صلة. ابحث عن هذه السرعة على مخطط الطور، ثم اقرأ كسب الحلقة عند السرعة نفسها على مخطط الكسب. هل هو أقل من 1 أم أعلى منه؟

الفصل 12: العودة إلى الدش ↗

حرّك المنزلق حتى يتجاوز المجموع 180°، بين 0.45 و0.46 راديان/ثانية. عند العبور بالضبط، نحو 0.457 راديان/ثانية، يعود التأرجح بحجم 0.72 من حجمه، فكان يمكن لليد أن تكون أقوى بنحو 1 ÷ 0.72 ≈ 1.4 مرة قبل أن تبلغ الحافة.

كلمات جديدة

هامش الكسب: كم مرة يمكنك مضاعفة كسب المتحكم قبل أن يصل كسب الحلقة إلى 1 عند سرعة −180°. أقل من 1 = فوق الهاوية.
هامش الطور: عند السرعة التي يكون فيها كسب الحلقة 1 تمامًا، كم من التأخر الإضافي يلزم للوصول إلى −180°. الأكثر أكثر أمانًا؛ و45°–60° هدف شائع (الاستطراد أدناه يشرح السبب).

هامش الطور تخميدٌ متنكر

جون: لماذا من 45° إلى 60°؟ من اختار ذلك؟

في الحلقات البسيطة من الفصل 7، يتحرك هامش الطور ونسبة التخميد معًا. نحو 43° يرافق ζ = 0.4 (تجاوز 25%)، و52° يرافق ζ = 0.5 (16%)، و59° يرافق ζ = 0.6 (أقل من 10%). هذه الأزواج دقيقة لأبسط حلقة في الفصل 7، $\omega_n^2/\big(s(s+2\zeta\omega_n)\big)$ . إذن هامش 45°–60° يعني ζ ≈ 0.4–0.6: أكثر ارتدادًا من 0.7–1 في الفصل 7. ومع ذلك فهو هدف شائع، لأن الهامش أيضًا مسافة أمان من الجرف (مقدار ما قد يخطئ النموذج)، لا مجرد وعد بشأن التجاوز. ومع وجود أنبوب في الحلقة يبقى دليلًا جيدًا، لا قانونًا دقيقًا.

في دشّنا، يبلغ طور $G_{\circ}$ القيمة −180° عند $\omega \approx 0.457$ راديان/ثانية، أي فترة تذبذب نحو 13.8 ثانية. اليد العادية ( $k = 0.008$ ) هامش كسبها نحو 1.4 فقط، وهامش طورها نحو 22° (حيث كسب حلقتها 1، عند 0.34 راديان/ثانية، يكون الطور −158°). قريبة من الحافة! اليد الأقوى تضاعف الكسب: 1.4 ÷ 2 = 0.7، فوق الجرف. لهذا بالضبط لم تستقر أبدًا.

توقّع أولًا! افترض أن الأنبوب أطول بمرتين (تأخير 5 ثوانٍ). ماذا يحدث لأكبر كسب لليد يظل مستقرًا؟

الفصل 12: العودة إلى الدش ↗

ميكا: ضعف الأنبوب، ونحو نصف اليد. والتأرجحات عند الحافة تصبح أبطأ أيضًا: قرابة 24 ثانية!

ثيو: حسنًا، أتراجع عن كلامي. التأخير لا يجعل الأمور أبطأ فحسب، بل يلتهم هامش الطور، وهامش الطور هو الشيء الوحيد بينك وبين التذبذب.

شرح دشّ الفصل 0 الخاص بك

إليك تجربتك الخاصة من الفصل 0. إذا كانت تأرجحاتك قريبة من حد الاستقرار، فينبغي أن تكون فترتها قريبةً مما تتنبأ به الحلقة. لكن هناك فرقًا مهمًا: لا يتصرف جميع الناس بالطريقة نفسها عند التحكم بالمقبض.

الفصل 12: العودة إلى الدش ↗

يد السرعة (تدير المقبض ما دام الإحساس خاطئًا) متحكم I: حافتها عند تأرجحات مدتها 13.8 ثانية. يد الموضع («ساخن جدًا، ضع المقبض هنا») تتصرف مثل متحكم P: $G_{\circ}(s) = 45K\,e^{-2.5s}/(s+1)$ . يد الموضع ليس لها كومة، فلا 90° خاصة بها: على الأنبوب والمُنعِّم وحدهما أن يبلغا 180°، وهذا يحدث عند تذبذبة أسرع، $\omega \approx 0.95$ راديان/ثانية (سرعة القلب التي وجدتها جون). هذا يعني تأرجحات مدتها 6.6 ثانية، لأي $K$ أكبر من نحو 0.031 مقبض لكل درجة مئوية. معظم الناس مزيج من الاثنين. (اليد الآلية العادية تتأرجح أبطأ قليلًا من 13.8 ثانية، نحو 15.5 ثانية، لأنها تقع داخل الحافة مباشرة؛ أما تأرجحات اليد الأقوى، نحو 11.6 ثانية، فتنضغط لأن المقبض يصطدم بطرفيه.)

الفصل 12: العودة إلى الدش ↗

ميكا: إذن حين كنتُ أنتر المقبض مباشرة إلى موضع ما، كنتُ يد موضع. تأرجحاتي ذات الست ثوانٍ كانت من الأنبوب والمُنعِّم، لا من أصابع خرقاء.

صمّم دشًا روبوتيًا

لنتفوّق الآن على أيدينا. يستخدم الروبوت PI: يتفاعل P مع خطأ درجة الحرارة الحالي، ويتذكر I الخطأ المتراكم. معًا يجعلان المقبض يستجيب بسرعة بينما يجد I الوضع الدائم الصحيح.

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

الفصل 12: العودة إلى الدش ↗

جون: …إنه يتذبذب ذهابًا وإيابًا، ساخن–بارد–ساخن، كل 6 ثوانٍ، مثل يد الموضع في الفصل 0. المكاسب الكبيرة دفعت سرعة الكسب 1 إلى ما بعد سرعة −180°. لا هامش طور على الإطلاق.

ثيو: الطائرة المسيّرة احتملت مكاسب كبيرة لأنها لم يكن فيها تأخر تقريبًا. هنا، المتانة تعني اللطف: اترك هامشًا للتأخير. جرّبوا مكاسبي اللطيفة، نحو $K_p \approx 1$ و $K_i \approx 0.6$ .

P يتفاعل مع الآن، وI مع الكومة. إضافة قليل من الآن تحرّك المقبض أبكر قليلًا، وهذا يعيد بعض الطور: تقدّم الطور، نحو 37° عند سرعة −180° القديمة مع مكاسب ثيو.

متحكم PI بهامش جيد يرتاح في حوالي 9.4 ثانية، دون تأرجحات كبيرة متكررة. هذا يتفوق على اليد العادية (التي أمضت معظم الدقيقة وهي تتأرجح قبل أن تهدأ) وعلى معظم البشر. لكنه لا يستطيع التغلب على الفيزياء: لا شيء يمكنه التفاعل مع ماء لم يصل بعد. لدى المهندسين حيلة للتأخيرات الطويلة، هي متنبئ سميث، الذي يشغّل محاكاة صغيرة للأنبوب داخل المتحكم. هذه قصة لدفتر آخر.

المتانة هي الدرس الحقيقي. نموذجك للدش ليس مثاليًا أبدًا: قد يكون الأنبوب أطول، والسخان أضعف. الهامش يقول كم يمكن أن يخطئ النموذج، برقم تستطيع أن تحسّه. هامش الطور (بالراديان) مقسومًا على سرعة الكسب 1 هو هامش التأخير: كم ثانية من الأنبوب الإضافي تستطيع الحلقة أن تتحمل.

الفصل 12: العودة إلى الدش ↗

ثيو: اليد العادية لا تتحمل إلا نحو 1.1 ثانية من الأنبوب الإضافي. روبوتي اللطيف يتحمل 3.5 ثانية. هذا هو معنى المتانة.

ما فهمناه الآن

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

تأخير $L$ يؤخر التذبذبة بمقدار $\omega L$ : لا شيء للتذبذبات البطيئة، والكثير للتذبذبات السريعة.
الكومة (I) تؤخر كل تذبذبة 90°؛ والمُنعِّم حتى 90°؛ والأنبوب $\omega L$ ، بلا حد.
إدخال موجات جيبية وقياس الكسب والطور يعطي الاستجابة الترددية؛ وعند رسمها، تكون مخطط بود، وهي مجرد $G(i\omega)$ .
بالنسبة لحلقة الدش البسيطة هذه، يبدأ التذبذب الذاتي المستمر عندما تصل الحلقة إلى طور −180° مع كسب 1؛ والكسب الأكبر يجعل التأرجحات تنمو حتى تتدخل الحدود.
يقيس هامش الكسب والطور المسافة إلى تلك الهاوية. التأخير الأطول يعني هامشًا أقل، لذا تحتاج إلى تحكم ألطف.
كانت تأرجحات الفصل 0 عبارة عن متحكم (نحن) بهامش صغير جدًا.

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

الفصل 12: العودة إلى الدش ↗