رياضيات نظرية التحكم: لابلاس والتغذية الراجعة

ما هي نظرية التحكم؟

تخيّل أنك تشغّل فرن صناعي في معمل زجاج بدمشق — تحتاج لتثبيت الحرارة عند 1200°C بدقة. الفرن يفقد حرارة من الجدران، والمادة الخام الباردة تدخل باستمرار، والوقود يتغير ضغطه. بدون نظام تحكم، الحرارة ستتذبذب بعنف. مع نظام PID محسّن رياضياً، تبقى الحرارة ضمن ±2°C من القيمة المطلوبة.

نظرية التحكم (Control Theory) هي الإطار الرياضي لتصميم أنظمة تحافظ تلقائياً على قيم مطلوبة رغم الاضطرابات. كل نظام تحكم في المصنع — منظم حرارة، متحكم سرعة محرك، ضابط ضغط — مبني على هذه الرياضيات.

تحويل لابلاس: من الزمن إلى التردد

تحويل لابلاس (Laplace Transform) هو الأداة الرياضية الأساسية في نظرية التحكم. يحوّل المعادلات التفاضلية (صعبة الحل) إلى معادلات جبرية (سهلة الحل).

التعريف:

F(s) = L{f(t)} = ∫₀^∞ f(t) × e^(-st) dt

حيث s = σ + jω متغير عقدي (complex variable).

التحويلات الأساسية:

القوة الحقيقية: تحويل المشتقات:

L{df/dt} = sF(s) - f(0)
L{d²f/dt²} = s²F(s) - sf(0) - f'(0)

المشتقة تصبح ضرباً بـ s — المعادلة التفاضلية تتحول لمعادلة جبرية عادية!

مثال: معادلة نابض مع تخميد:

m × d²x/dt² + c × dx/dt + k × x = F(t)

بعد تحويل لابلاس (مع شروط ابتدائية صفرية):

(ms² + cs + k) × X(s) = F(s)

دالة التحويل: بصمة النظام

دالة التحويل (Transfer Function) هي نسبة خرج النظام إلى دخله في مجال لابلاس:

G(s) = Y(s) / U(s)

حيث Y(s) هو الخرج و U(s) هو الدخل. دالة التحويل تصف سلوك النظام بالكامل — هي "بصمته" الرياضية.

أنظمة شائعة:

نظام من الدرجة الأولى (خزان حراري، دارة RC):

G(s) = K / (τs + 1)

حيث K = الربح الساكن (الاستجابة النهائية) و τ = ثابت الزمن (سرعة الاستجابة). بعد 5τ تصل الاستجابة لـ 99% من قيمتها النهائية.

نظام من الدرجة الثانية (نابض-كتلة-مخمّد، محرك مع حمل):

G(s) = ωₙ² / (s² + 2ζωₙs + ωₙ²)

حيث:

• ωₙ = التردد الطبيعي (rad/s) — سرعة الاستجابة
• ζ = نسبة التخميد (بلا وحدة) — مفتاح السلوك:

القيمة المثلى عملياً: ζ ≈ 0.7 — تعطي توازناً بين السرعة وقلة التذبذب (تجاوز أقل من 5%).

المخططات الصندوقية: لغة مهندس التحكم

المخططات الصندوقية (Block Diagrams) تمثل الأنظمة بيانياً. كل صندوق يمثل دالة تحويل، والأسهم تمثل الإشارات.

حلقة التحكم المغلقة الأساسية:

  R(s) --→ [+] --→ [G_c(s)] --→ [G_p(s)] --→ Y(s)
            ↑ -                                |
            |                                  |
            +-------- [H(s)] ←-----------------+

حيث:

• R(s) = الإشارة المرجعية (القيمة المطلوبة)
• G_c(s) = دالة تحويل المتحكم (PID مثلاً)
• G_p(s) = دالة تحويل العملية (المصنع/المحرك/الفرن)
• H(s) = دالة تحويل المستشعر (حساس الحرارة/السرعة)
• Y(s) = الخرج الفعلي

دالة التحويل للحلقة المغلقة:

T(s) = G_c(s) × G_p(s) / (1 + G_c(s) × G_p(s) × H(s))

قواعد تبسيط المخططات:

• تسلسل: G₁(s) × G₂(s) — اضرب دوال التحويل
• توازي: G₁(s) + G₂(s) — اجمع دوال التحويل
• تغذية راجعة: استخدم الصيغة أعلاه

مخطط بود: رؤية التردد

مخطط بود (Bode Plot) يرسم استجابة النظام عند ترددات مختلفة. يتكون من رسمين:

1. مخطط السعة (Magnitude): 20 log₁₀|G(jω)| بوحدة dB مقابل التردد
2. مخطط الطور (Phase): ∠G(jω) بالدرجات مقابل التردد

لنظام من الدرجة الأولى G(s) = K/(τs+1):

• عند ترددات منخفضة (ω << 1/τ): السعة = 20 log K dB، الطور = 0°
• عند تردد القطع (ω = 1/τ): السعة تنخفض 3dB، الطور = -45°
• عند ترددات عالية (ω >> 1/τ): السعة تنخفض بمعدل -20 dB/decade، الطور → -90°

قراءة مخطط بود — المعلومات الحرجة:

• هامش الربح (Gain Margin): كم يمكن زيادة الربح قبل أن يفقد النظام استقراره — يُقاس عند التردد الذي يكون فيه الطور = -180°
• هامش الطور (Phase Margin): كم يبعد الطور عن -180° عند التردد الذي تكون فيه السعة = 0 dB
• القاعدة: هامش ربح > 6 dB وهامش طور > 30° يضمنان استقراراً مريحاً

مخطط نايكويست وموضع الجذور

مخطط نايكويست (Nyquist Plot) يرسم G(jω) في المستوي العقدي (الجزء الحقيقي مقابل الجزء التخيلي) مع تغيّر التردد من 0 إلى ∞.

معيار نايكويست للاستقرار: النظام المغلق مستقر إذا وفقط إذا كان عدد مرات إحاطة منحنى نايكويست للنقطة (-1, 0) يساوي عدد أقطاب دالة الحلقة المفتوحة الواقعة في نصف المستوي الأيمن.

ببساطة: إذا كانت الحلقة المفتوحة مستقرة (لا أقطاب في نصف المستوي الأيمن)، فالنظام المغلق مستقر إذا لم يُحِط المنحنى بالنقطة (-1, 0).

مخطط موضع الجذور (Root Locus):

يرسم مسار أقطاب النظام المغلق مع تغيّر ربح المتحكم K من 0 إلى ∞.

القواعد الأساسية:

• يبدأ من أقطاب الحلقة المفتوحة (عند K=0)
• ينتهي عند أصفار الحلقة المفتوحة (عند K→∞)
• النظام مستقر طالما بقيت جميع الأقطاب في نصف المستوي الأيسر (Re(s) < 0)
• عندما يعبر قطب المحور التخيلي — يفقد النظام استقراره

معايير الاستقرار

معيار راوث-هورويتز (Routh-Hurwitz)

طريقة جبرية لتحديد الاستقرار بدون حساب الجذور فعلياً. من المعادلة المميزة:

s⁴ + a₃s³ + a₂s² + a₁s + a₀ = 0

نبني جدول راوث — إذا كانت جميع عناصر العمود الأول موجبة، فالنظام مستقر. أي تغيير في الإشارة يعني وجود قطب في نصف المستوي الأيمن (عدم استقرار).

شرط ضروري (لكن غير كافٍ): جميع معاملات المعادلة المميزة يجب أن تكون موجبة وموجودة. إذا كان أي معامل سالباً أو صفراً — النظام غير مستقر حتماً.

رياضيات التحكم PID

متحكم PID هو الأكثر استخداماً في الصناعة (أكثر من 95% من حلقات التحكم). دالة تحويله:

G_c(s) = K_p + K_i/s + K_d × s

أو بالصيغة المعيارية:

G_c(s) = K_p × (1 + 1/(T_i × s) + T_d × s)

حيث:

• K_p = الربح التناسبي — يتناسب مع الخطأ الحالي
• K_i = K_p/T_i = الربح التكاملي — يزيل الخطأ الثابت بالتراكم
• K_d = K_p × T_d = الربح التفاضلي — يتنبأ بالخطأ المستقبلي

تأثير كل حد:

ضبط معاملات PID

طريقة زيغلر-نيكولز (Ziegler-Nichols)

أشهر طريقة تجريبية لضبط PID — لا تحتاج نموذجاً رياضياً، فقط تجربة على النظام الحقيقي:

الطريقة الأولى — استجابة الخطوة:

1. ضع المتحكم في وضع يدوي
2. طبّق تغييراً مفاجئاً (خطوة) في الدخل
3. من منحنى الاستجابة، قِس زمن التأخير L وثابت الزمن T

الطريقة الثانية — التذبذب الحرج:

1. ابدأ بتحكم تناسبي فقط (K_i = K_d = 0)
2. زِد K_p تدريجياً حتى يبدأ النظام بالتذبذب المستمر
3. سجّل الربح الحرج K_u ودور التذبذب T_u

ملاحظة مهمة: زيغلر-نيكولز تعطي نقطة بداية جيدة لكن تحتاج ضبط دقيق يدوي — التجاوز الناتج عادة 25% وهو كثير لبعض التطبيقات. في المصانع الحديثة، تُستخدم طرق أكثر تطوراً مثل التحكم الأمثل (Optimal Control) و التحكم التنبؤي (Model Predictive Control - MPC).

من النظرية إلى المصنع

في الواقع الصناعي السوري، معظم أنظمة التحكم تعمل بمتحكمات PLC مع خوارزميات PID مضمّنة. فهم الرياضيات ليس ترفاً — بل ضرورة لـ:

• تشخيص المشاكل: لماذا يتذبذب الخرج؟ هل المشكلة في الربح أم في التخميد؟
• ضبط الأداء: تقليل زمن الاستقرار وتقليل التجاوز
• اختيار المتحكم المناسب: متى يكفي P؟ متى نحتاج PI؟ ومتى نحتاج PID كاملاً؟
• تصميم أنظمة جديدة: اختيار المستشعرات والمشغلات بناءً على متطلبات الأداء الرياضية