في نافذة الطباعة اختر «حفظ بصيغة PDF» لحفظ نسخة على جهازك.
1القوانين وطرق الحل
الاستقرار: راوث–هرفيتز والمحل الهندسي للجذور
الاستقرار
جميع جذور كثير الحدود المميِّز في النصف الأيسر من المستوى
جدول راوث، عنصر الصف s¹ في كثير حدود تكعيبي
b1 = (a2a1 − a3a0)/a2
جذور النصف الأيمن من المستوى
عدد تغيّرات الإشارة في العمود الأول
كثير الحدود التكعيبي s³ + a₂s² + a₁s + a₀
مستقر حين a > 0 لجميع المعاملات وa₂·a₁ > a₀
K/(s(s + a)(s + b))
Kmax = ab(a + b), ω = √(ab)
المعادلة المميزة
1 + K·G(s) = 0
الخطوط المقاربة
σa = (Σp − Σz)/(n − m), θ = (2q + 1)·180°/(n − m)
نقطتا الانفصال والالتقاء
dK/ds = 0
شرط المقدار
K = 1/|G(s)|
زاوية المغادرة
θd = 180° − Σ∠(from other poles) + Σ∠(from zeros)
الطريقة
اكتب كثير الحدود المميِّز للحلقة المغلقة: مقام 1 + K·G(s) = 0، بعد فكّه.
ضع المعاملات بالتناوب في الصفين الأولين من جدول راوث وابنِ بقية الصفوف.
عُدّ تغيّرات الإشارة في العمود الأول: كل تغيّر قطب في النصف الأيمن من المستوى. ومع وجود K في الجدول اشترط أن يبقى كل عنصر في العمود الأول موجبًا.
للمحل الهندسي للجذور عيّن الأقطاب والأصفار وقطع المحور الحقيقي والخطوط المقاربة من المركز ونقاط الانفصال حيث dK/ds = 0. واستخدم K = 1/|G(s)| لإيجاد الكسب عند نقطة.
خطأ الحالة المستقرة والتحكم PID
نوع النظام
عدد أقطاب G(s) عند s = 0
ثوابت الخطأ
Kp = lim G(s), Kv = lim s·G(s), Ka = lim s²·G(s) (s → 0)
الخطوة A
ess = A/(1 + Kp)
المنحدر A·t
ess = A/Kv
القطع المكافئ A·t²/2
ess = A/Ka
المتحكم PID
C(s) = Kp + Ki/s + Kd·s
الهدف من الرتبة الثانية
s² + 2ζωns + ωn²
التجاوز والتوقيت
%OS = 100·e−ζπ/√(1−ζ²), tp = π/ωd, ts ≈ 4/(ζωn)
الطريقة
عُدّ أقطاب الحلقة المفتوحة C(s)·G(s) عند s = 0: فذلك هو نوع النظام.
يُتحكَّم في نظام G(s) = 35/(s + 4) في حلقة تغذية راجعة سالبة كسب مستشعرها H = 0,2. أوجد الثابت الزمني τ للحلقة المغلقة.
G(s) = 35/(s + 4)H = 0,2
τ = s
4.المخططات الصندوقيةصعبة
متحكّم تناسبي K يقود النظام G(s) = 1/(s + 7) في حلقة تغذية راجعة سالبة أحادية. تدخل خطوة اضطراب مقدارها 5 عند دخل النظام بينما R = 0. أوجد أصغر كسب K يُبقي خرج الحالة المستقرة عند 0,007143 أو دونه.
G(s) = 1/(s + 7)D = 5css ≤ 0,007143
K =
5.الاستجابة الخطويةمتوسطة
لنظام دالة النقل G(s) = 57/(s² + 6s + 19) ويتلقى دخلًا خطويًا وحدويًا. أوجد زمن الاستقرار ts عند 2 % للاستجابة الخطوية باستخدام ts ≈ 4/(ζωₙ).
G(s) = 57/(s² + 6s + 19)
ts = s
6.الاستجابة الخطويةصعبة
لحلقة تغذية عكسية سالبة وحدوية G(s) = K/(s(s + 2)). أوجد الكسب K الذي يعطي استجابة خطوية بتجاوز 25 %.
يغادر المحل الهندسي للجذور لـ L(s) = K/(s(s + 1)(s + 3)) المحور الحقيقي بين s = 0 وs = −1. أوجد نقطة الانفصال sb.
L(s) = K/(s(s + 1)(s + 3))
sb =
9.خطأ الحالة المستقرة والتحكم PIDمتوسطة
لحلقة تغذية عكسية وحدوية G(s) = K/((s + 5)(s + 10)). أوجد الكسب K الذي يعطي خطأ حالة مستقرة قدره 2 % لدخل خطوة.
G(s) = K/((s + 5)(s + 10))ess = 2 %
K =
10.خطأ الحالة المستقرة والتحكم PIDصعبة
يتحكم متحكم PI بالصيغة C(s) = 1 + 2/s في المحطة G(s) = 5/((s + 2)(s + 4)) ضمن حلقة تغذية عكسية وحدوية. أوجد خطأ الحالة المستقرة ess لدخل المنحدر r(t) = 2t.
C(s) = 1 + 2/sG(s) = 5/((s + 2)(s + 4))r(t) = 2t
ess =
11.مخططات بودي والهوامشمتوسطة
أوجد المقدار الدقيق |G(jω)| لـ G(s) = 196/((s + 7)(s + 14)) عند ω = 8 rad/s، بوحدة dB.
G(s) = 196/((s + 7)(s + 14))ω = 8 rad/s
|G(jω)| = dB
12.مخططات بودي والهوامشصعبة
حلقة تغذية عكسية وحدوية دالة نقلها المفتوحة L(s) = K/(s(s + 3)(s + 4)). أوجد الكسب K الذي يعطي هامش كسب قدره 5 dB.
L(s) = K/(s(s + 3)(s + 4))GM = 5 dB
K =
3المحاليل
1.تحويلات لابلاس
الإجابة: A = 5,833
تحويل لابلاس: F(s) = ∫₀∞ f(t)·e−st dt يحوّل المشتقات إلى ضرب في s: ℒ{f′} = sF(s) − f(0)، ℒ{f″} = s²F(s) − s·f(0) − f′(0). الأزواج القياسية: 1 ↔ 1/s، e−at ↔ 1/(s + a)، t·e−at ↔ 1/(s + a)²، sin ωt ↔ ω/(s² + ω²)، cos ωt ↔ s/(s² + ω²)، والإزاحة s → s + σ تضرب دالة الزمن في e−σt. لعكس F(s) كسرية، جزّئها إلى كسور جزئية: باقي القطب البسيط p هو A = [(s − p)·F(s)] عند s = p. نظرية القيمة النهائية: f(∞) = lims→0 s·F(s)، وتصحّ فقط عندما تقع كل أقطاب s·F(s) في نصف المستوى الأيسر. نظرية القيمة الابتدائية: f(0⁺) = lims→∞ s·F(s).
حلّل المقام
s² + 8s + 7 = (s + 1)(s + 7)
F(s) = A/(s + 1) + B/(s + 7)
قاعدة التغطية
اضرب في (s + 1) وضع s = −1: يختفي حدّ B.
A = (3s + 38)/(s + 7) |s = −1 = 35/6 = 5,833
الإجابة
A = 5,833
(B = −2,833، لذا f(t) = 5,833·e−t − 2,833·e−7t)
2.تحويلات لابلاس
الإجابة: y(0.2) = 3,942
تحويل لابلاس: F(s) = ∫₀∞ f(t)·e−st dt يحوّل المشتقات إلى ضرب في s: ℒ{f′} = sF(s) − f(0)، ℒ{f″} = s²F(s) − s·f(0) − f′(0). الأزواج القياسية: 1 ↔ 1/s، e−at ↔ 1/(s + a)، t·e−at ↔ 1/(s + a)²، sin ωt ↔ ω/(s² + ω²)، cos ωt ↔ s/(s² + ω²)، والإزاحة s → s + σ تضرب دالة الزمن في e−σt. لعكس F(s) كسرية، جزّئها إلى كسور جزئية: باقي القطب البسيط p هو A = [(s − p)·F(s)] عند s = p. نظرية القيمة النهائية: f(∞) = lims→0 s·F(s)، وتصحّ فقط عندما تقع كل أقطاب s·F(s) في نصف المستوى الأيسر. نظرية القيمة الابتدائية: f(0⁺) = lims→∞ s·F(s).
حوّل المعادلة
(s² + 5s + 4)·Y(s) = 4s + 20 + 12/s
Y(s) = (4s² + 20s + 12)/(s(s + 1)(s + 4))
الكسور الجزئية
A = 12/(1·4) = 3
B = −4/((−1)·3) = 1,333
C = −4/((−4)·(−3)) = −0,3333
النتيجة
y(t) = 3 + 1,333·e−t − 0,3333·e−4t
y(0,2) = 3,942
الإجابة
y(0,2) = 3,942
تحقّق: y(0) = A + B + C = 4.
3.المخططات الصندوقية
الإجابة: τ = 0,09091 s
جبر المخططات الصندوقية: الصناديق على التوالي تُضرب (G₁·G₂)، والصناديق على التوازي تُجمع (G₁ + G₂)، والحلقة ذات المسار الأمامي G ومسار التغذية الراجعة H تُختزل إلى T = G/(1 + G·H) للتغذية الراجعة السالبة و G/(1 − G·H) للتغذية الراجعة الموجبة. اختزل الحلقات الداخلية أولًا، ثم اتجه إلى الخارج. أقطاب الحلقة المغلقة هي جذور المعادلة المميّزة 1 + G(s)H(s) = 0. لحلقة مستقرة، دخل خطوة مقداره A يعطي خرج الحالة المستقرة A·T(0)؛ ومع التغذية الراجعة الأحادية يكون خطأ الحالة المستقرة لخطوة A/(1 + Kp) مع Kp = G(0).
دالة تحويل الحلقة المغلقة
T(s) = G/(1 + G·H) = 35/(s + 4 + 7) = 35/(s + 11)
النتيجة
τ = 1/11 = 0,09091 s
الإجابة
τ = 0,09091 s
4.المخططات الصندوقية
الإجابة: K = 693
جبر المخططات الصندوقية: الصناديق على التوالي تُضرب (G₁·G₂)، والصناديق على التوازي تُجمع (G₁ + G₂)، والحلقة ذات المسار الأمامي G ومسار التغذية الراجعة H تُختزل إلى T = G/(1 + G·H) للتغذية الراجعة السالبة و G/(1 − G·H) للتغذية الراجعة الموجبة. اختزل الحلقات الداخلية أولًا، ثم اتجه إلى الخارج. أقطاب الحلقة المغلقة هي جذور المعادلة المميّزة 1 + G(s)H(s) = 0. لحلقة مستقرة، دخل خطوة مقداره A يعطي خرج الحالة المستقرة A·T(0)؛ ومع التغذية الراجعة الأحادية يكون خطأ الحالة المستقرة لخطوة A/(1 + Kp) مع Kp = G(0).
استجابة الاضطراب
css = D·G(0)/(1 + K·G(0)) = 5·1/(7 + 1K)
حلّ لإيجاد K
7 + 1K = 5·1/0,007143 = 700
K = (700 − 7)/1 = 693
الإجابة
K = 693
5.الاستجابة الخطوية
الإجابة: ts = 1,333 s
الاستجابة للدخل الخطوي: النظام من الرتبة الأولى G(s) = K/(τs + 1) يستجيب لخطوة مقدارها U بالعلاقة y(t) = K·U·(1 − e−t/τ): يبلغ 63,2 % من القيمة النهائية K·U بعد ثابت زمني واحد τ و98 % بعد 4τ. والنظام من الرتبة الثانية G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) له القطبان s = −ζωₙ ± jωd مع ωd = ωₙ√(1 − ζ²). وعندما 0 < ζ < 1 تتجاوز الاستجابة الخطوية القيمة النهائية: النسبة المئوية للتجاوز %OS = 100·e−ζπ/√(1−ζ²)، وزمن الذروة tp = π/ωd، وزمن الاستقرار عند 2 % هو ts ≈ 4/(ζωₙ). قارن المقام مع s² + 2ζωₙs + ωₙ² لقراءة ωₙ وζ. وفي نظام كتلة–نابض–مخمّد ωₙ = √(k/m) وζ = c/(2√(km))؛ وفي دائرة RLC على التسلسل ωₙ = 1/√(LC) وζ = (R/2)·√(C/L).
قارن مع الصيغة القياسية
يجب أن يطابق المقام s² + 2ζωₙs + ωₙ².
ωₙ² = 19 → ωₙ = 4,359 rad/s
2ζωₙ = 6 → ζ = 6/(2·4,359) = 0,6882
النتيجة
ζωₙ = 3 → ts ≈ 4/3 = 1,333 s
الإجابة
ts = 1,333 s
6.الاستجابة الخطوية
الإجابة: K = 6,136
الاستجابة للدخل الخطوي: النظام من الرتبة الأولى G(s) = K/(τs + 1) يستجيب لخطوة مقدارها U بالعلاقة y(t) = K·U·(1 − e−t/τ): يبلغ 63,2 % من القيمة النهائية K·U بعد ثابت زمني واحد τ و98 % بعد 4τ. والنظام من الرتبة الثانية G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) له القطبان s = −ζωₙ ± jωd مع ωd = ωₙ√(1 − ζ²). وعندما 0 < ζ < 1 تتجاوز الاستجابة الخطوية القيمة النهائية: النسبة المئوية للتجاوز %OS = 100·e−ζπ/√(1−ζ²)، وزمن الذروة tp = π/ωd، وزمن الاستقرار عند 2 % هو ts ≈ 4/(ζωₙ). قارن المقام مع s² + 2ζωₙs + ωₙ² لقراءة ωₙ وζ. وفي نظام كتلة–نابض–مخمّد ωₙ = √(k/m) وζ = c/(2√(km))؛ وفي دائرة RLC على التسلسل ωₙ = 1/√(LC) وζ = (R/2)·√(C/L).
دالة تحويل الحلقة المغلقة
T(s) = K/(s² + 2s + K)
نسبة التخميد من التجاوز
ζ = −ln(0,25)/√(π² + ln²(0,25)) = 0,4037
الكسب
ωₙ = 2/(2ζ) = 2,477 rad/s
K = ωₙ² = 6,136
الإجابة
K = 6,136
7.الاستقرار: راوث–هرفيتز والمحل الهندسي للجذور
الإجابة: ω = 5,657 rad/s
معيار راوث–هرفيتز: يكون النظام الخطي مستقرًا حين يقع كل جذر من جذور كثير حدوده المميِّز في النصف الأيسر من المستوى s. اكتب المعاملات بالتناوب في الصفين الأولين من جدول راوث. يأتي كل عنصر في صف جديد من الصفين اللذين فوقه؛ وبالنسبة إلى كثير حدود تكعيبي a₃s³ + a₂s² + a₁s + a₀ يكون عنصر الصف s¹ هو b₁ = (a₂·a₁ − a₃·a₀)/a₂. عدد تغيّرات الإشارة في العمود الأول يساوي عدد الجذور في النصف الأيمن من المستوى. وحين يدخل الكسب K في المعاملات يجب أن يبقى كل عنصر في العمود الأول موجبًا: الكسب الذي يجعل عنصر الصف s¹ صفرًا هو حدّ الاستقرار، ثم يعطي الصف s² تردد التذبذب المستمر.
جدول روث
s³ + 12s² + 32s + K = 0
s³: 1, 32
s²: 12, K
s¹: (12·32 − K)/12
s⁰: K
الكسب الحرج
K = 12·32 = 384
المعادلة المساعدة
عند الكسب الحرج يكون الصف s¹ صفرًا؛ ويعطي الصف s² الأقطاب الواقعة على المحور التخيلي.
12s² + 384 = 0 → s = ±j√32
ω = 5,657 rad/s
الإجابة
ω = 5,657 rad/s
8.الاستقرار: راوث–هرفيتز والمحل الهندسي للجذور
الإجابة: sb = −0,4514
المحل الهندسي للجذور: حين يتغير الكسب K من 0 إلى ∞ تتحرك أقطاب الحلقة المغلقة، وهي جذور 1 + K·G(s) = 0، من أقطاب الحلقة المفتوحة إلى أصفار الحلقة المفتوحة أو إلى ما لا نهاية. تنتمي نقطة على المحور الحقيقي إلى المحل الهندسي حين يكون عدد الأقطاب والأصفار الحقيقية على يمينها فرديًا. تتبع الفروع n − m المتجهة إلى ما لا نهاية خطوطًا مقاربة بالزوايا (2q + 1)·180°/(n − m)، وتنطلق من المركز σₐ = (Σ الأقطاب − Σ الأصفار)/(n − m). تغادر الفروع المحور الحقيقي أو تنضم إليه حيث dK/ds = 0. يحدد شرط الزاوية ∠G(s) = 180° الشكل، ويعطي شرط المقدار K = 1/|G(s)| الكسب عند أي نقطة من المحل الهندسي.
الكسب على طول المحور الحقيقي
K(s) = −s(s + 1)(s + 3) = −(s³ + 4s² + 3s)
النقاط الحرجة
dK/ds = −(3s² + 8s + 3) = 0
s = (−8 ± √(64 − 36))/6 = −0,4514, −2,215
اختر النقطة الواقعة على المحل الهندسي
يقع المحل الهندسي على المحور الحقيقي بين 0 و−1؛ أما الجذر الآخر فليس على المحل الهندسي.
sb = −0,4514
الإجابة
sb = −0,4514
9.خطأ الحالة المستقرة والتحكم PID
الإجابة: K = 2450
خطأ الحالة المستقرة: في حلقة تغذية عكسية وحدوية بدالة نقل الحلقة المفتوحة G(s) يكون الخطأ E(s) = R(s)/(1 + G(s))، وبنظرية القيمة النهائية ess = lims→0 s·E(s). نوع النظام هو عدد أقطاب G(s) عند s = 0. ثوابت الخطأ هي Kp = lim G(s) وKv = lim s·G(s) وKa = lim s²·G(s)، كلها عندما s → 0. تعطي دالة خطوة بحجم A الخطأ ess = A/(1 + Kp)، وتعطي دالة منحدرة A·t الخطأ A/Kv، وتعطي دالة قطع مكافئ A·t²/2 الخطأ A/Ka. تتبع حلقة من النوع 0 دالة الخطوة بخطأ منتهٍ، أما حلقة من النوع 1 فتتبعها تمامًا وتتبع دالة المنحدر بخطأ منتهٍ.
ثابت الخطأ المطلوب
1/(1 + Kp) = 0,02 → Kp = 49
الكسب
K/50 = 49 → K = 2450
الإجابة
K = 2450
تحقق أيضًا من استقرار الحلقة المغلقة: في هذه الحلقة من الرتبة الثانية يكون أي K > 0 مستقرًا.
10.خطأ الحالة المستقرة والتحكم PID
الإجابة: ess = 1,6
التحكم P وPI وPID: يعمل المتحكم C(s) = Kp + Ki/s + Kd·s على الخطأ. الكسب التناسبي يجعل الحلقة أسرع ويقلل خطأ محطة من النوع 0، لكنه لا يزيله أبدًا. يضيف الحد التكاملي قطبًا عند s = 0: فيرتفع نوع النظام بواحد، وتُتبع الخطوة دون خطأ. ويضيف الحد التفاضلي تخميدًا. لاختيار الكسوب اكتب كثير الحدود المميِّز للحلقة المغلقة وقارنه بـ s² + 2ζωns + ωn². عندئذٍ %OS = 100·e−ζπ/√(1−ζ²)، وtp = π/ωd مع ωd = ωn√(1 − ζ²)، وts ≈ 4/(ζωn).
ثابت السرعة
s·C(s) = 1s + 2 → 2 عندما s → 0.
Kv = 2·5/8 = 1,25
خطأ المنحدر
ess = 2/1,25 = 1,6
الإجابة
ess = 1,6
يزيل المتحكم PI خطأ الخطوة تمامًا؛ ويُتبع المنحدر بتأخر ثابت.
11.مخططات بودي والهوامش
الإجابة: |G(jω)| = 1,165 dB
الاستجابة الترددية: عوّض s = jω في دالة النقل. المقدار بالديسيبل هو 20·log₁₀|H(jω)| والطور هو مجموع زوايا عوامل البسط ناقص زوايا عوامل المقام. مرشح RC أو RL من الرتبة الأولى له تردد الركن (−3 dB) عند ωc = 1/(RC) أو R/L، وfc = ωc/(2π)؛ ومرشح تمرير الترددات المنخفضة يعطي |H| = 1/√(1 + (ω/ωc)²) وφ = −atan(ω/ωc). وفي مخطط بودي التقريبي (الخطوط المستقيمة) ينحني كل قطب بالمقدار إلى الأسفل بـ 20 dB/decade عند تردد ركنه، وينحني كل صفر إلى الأعلى بـ 20 dB/decade، ويعطي المكامل 1/s ميلًا قدره −20 dB/decade عند كل الترددات. تردد تقاطع الكسب ωc: |L(jω)| = 1 (0 dB)؛ وهامش الطور PM = 180° + ∠L(jωc). تردد تقاطع الطور ω180: ∠L(jω) = −180°؛ وهامش الكسب GM = −20·log₁₀|L(jω180)| dB. تكون الحلقة المغلقة مستقرة حين يكون الهامشان موجبين.
الاستجابة الترددية: عوّض s = jω في دالة النقل. المقدار بالديسيبل هو 20·log₁₀|H(jω)| والطور هو مجموع زوايا عوامل البسط ناقص زوايا عوامل المقام. مرشح RC أو RL من الرتبة الأولى له تردد الركن (−3 dB) عند ωc = 1/(RC) أو R/L، وfc = ωc/(2π)؛ ومرشح تمرير الترددات المنخفضة يعطي |H| = 1/√(1 + (ω/ωc)²) وφ = −atan(ω/ωc). وفي مخطط بودي التقريبي (الخطوط المستقيمة) ينحني كل قطب بالمقدار إلى الأسفل بـ 20 dB/decade عند تردد ركنه، وينحني كل صفر إلى الأعلى بـ 20 dB/decade، ويعطي المكامل 1/s ميلًا قدره −20 dB/decade عند كل الترددات. تردد تقاطع الكسب ωc: |L(jω)| = 1 (0 dB)؛ وهامش الطور PM = 180° + ∠L(jωc). تردد تقاطع الطور ω180: ∠L(jω) = −180°؛ وهامش الكسب GM = −20·log₁₀|L(jω180)| dB. تكون الحلقة المغلقة مستقرة حين يكون الهامشان موجبين.