วิศวกรรมระบบควบคุม 1: ชุดสรุปเร่งรัด

สรุปสูตร ขั้นตอนวิธีทำที่สำคัญของทุกหัวข้อ และโจทย์ฝึก 12 ข้อ พร้อมเฉลยละเอียดไว้ท้ายเล่ม

โจทย์ใหม่กลับไปที่รายวิชา

ในหน้าต่างพิมพ์ ให้เลือก “บันทึกเป็น 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. เขียนพหุนามลักษณะเฉพาะของวงปิด: ตัวส่วนของ 1 + K·G(s) = 0 ที่คูณออกแล้ว
  2. ใส่สัมประสิทธิ์สลับกันลงในสองแถวแรกของตารางราอ์ท แล้วสร้างแถวที่เหลือ
  3. นับจำนวนครั้งที่เครื่องหมายเปลี่ยนในคอลัมน์แรก: แต่ละครั้งคือโพลหนึ่งตัวในครึ่งระนาบขวา เมื่อมี K อยู่ในตาราง ให้กำหนดว่าทุกค่าในคอลัมน์แรกต้องเป็นบวก
  4. สำหรับรูตโลกัส ให้ทำเครื่องหมายโพลและศูนย์ ช่วงบนแกนจริง เส้นกำกับจากเซนทรอยด์ และจุดแยกที่ 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)

วิธีการ

  1. นับโพลของวงเปิด C(s)·G(s) ที่ s = 0: นั่นคือชนิดของระบบ
  2. คำนวณค่าคงที่ความคลาดเคลื่อนที่ตรงกับอินพุต: Kp สำหรับขั้น Kv สำหรับแรมป์ Ka สำหรับพาราโบลา
  3. ใช้ ess = A/(1 + Kp), A/Kv หรือ A/Ka หรือแก้สมการหาอัตราขยายที่ต้องการ
  4. สำหรับข้อกำหนดเรื่องการหน่วงหรือความเร็ว ให้เขียนพหุนามลักษณะเฉพาะของวงปิดแล้วเทียบกับ s² + 2ζω_n·s + ω_n²

2โจทย์ฝึก

1.การแปลงลาปลาซปานกลาง
ReIm−2−4−6−8−10

จงหา f(t) สำหรับ F(s) = (3s + 27)/(s² + 2s + 2) ด้วยการทำเป็นกำลังสองสมบูรณ์ และหาค่า f(1)

F(s) = (3s + 27)/(s² + 2s + 2)t = 1 s
f(1) =
2.การแปลงลาปลาซยาก
ReIm−5−10−15−20

ระบบอธิบายด้วยสมการเชิงอนุพันธ์ y″ + 10y′ + 24y = u′ + 18u โดยมีอินพุต u(t) และเอาต์พุต y(t) จงหาฟังก์ชันถ่ายโอน G(s) = Y(s)/U(s) (เงื่อนไขเริ่มต้นเป็นศูนย์) และหาผลตอบสนองอิมพัลส์ g(t) ที่ t = 0.1 s จงหา g(0.1)

y″ + 10y′ + 24y = u′ + 18ut = 0.1 s
g(0.1) =
3.แผนภาพบล็อกปานกลาง
R(s)+3K+8s + 3G(s)C(s)4H₁−−

พืช G(s) = 8/(s + 3) มีวงป้อนกลับลบภายใน (ย่อย) ที่มีอัตราขยาย H₁ = 4 อัตราขยายของตัวควบคุม K = 3 และวงป้อนกลับลบหนึ่งหน่วยภายนอกปิดระบบ จงหาอัตราขยาย DC วงปิด T(0)

G(s) = 8/(s + 3)H₁ = 4K = 3
T(0) =
4.แผนภาพบล็อกยาก
R(s)+Ks(s + 8)(s + 9)G(s)C(s)−

พืช G(s) = K/(s(s + 8)(s + 9)) อยู่ในวงป้อนกลับลบหนึ่งหน่วย ใช้เกณฑ์เราท์–เฮอร์วิตซ์หาอัตราขยายที่มากที่สุด Kmax ที่วงปิดยังเสถียร

G(s) = K/(s(s + 8)(s + 9))ป้อนกลับลบหนึ่งหน่วย
Kmax =
5.การตอบสนองต่อสัญญาณขั้นบันไดปานกลาง
t (s)y(t)2tp​ = ?

ระบบหนึ่งมีฟังก์ชันถ่ายโอน G(s) = 190/(s² + 10s + 95) และรับอินพุตขั้นบันไดหนึ่งหน่วย จงหาเวลาถึงจุดสูงสุด tp ของการตอบสนองต่อสัญญาณขั้นบันได

G(s) = 190/(s² + 10s + 95)
tp = s
6.การตอบสนองต่อสัญญาณขั้นบันไดยาก
k = 500 N/mc = 30 N·s/mm = 5 kgF = 5 Nx(t)

มวล m = 5 kg แขวนอยู่กับสปริง k = 500 N/m และตัวหน่วงแบบหนืด c = 30 N·s/m (m·x″ + c·x′ + k·x = F) เริ่มจากหยุดนิ่งในสมดุล แล้วออกแรงคงที่ F = 5 N กะทันหัน จงหาเวลา tp ที่การกระจัดมีค่ามากที่สุด

m = 5 kgk = 500 N/mc = 30 N·s/mF = 5 N
tp = s
7.เสถียรภาพ: เกณฑ์ราอ์ท–เฮอร์วิตซ์และรูตโลกัสปานกลาง
R(s)+Ks + 7s(s + 1)(s + 4)C(s)−

วงฟีดแบ็กหนึ่งหน่วยมี L(s) = K(s + 7)/(s(s + 1)(s + 4)) โดย K > 0 จงหาอัตราขยายสูงสุด Kmax ที่วงปิดเสถียร

L(s) = K(s + 7)/(s(s + 1)(s + 4))
Kmax =
8.เสถียรภาพ: เกณฑ์ราอ์ท–เฮอร์วิตซ์และรูตโลกัสยาก
ReIm−2−4−6−8−10−12j2sb​

root locus ของ L(s) = K/((s + 5)(s + 11)) แสดงไว้ จงหาอัตราขยาย K ที่จุดแยกออก sb ซึ่งโพลวงปิดสองตัวพบกันบนแกนจริง

L(s) = K/((s + 5)(s + 11))
K =
9.ค่าความคลาดเคลื่อนสถานะคงตัวและการควบคุม PIDปานกลาง
r(t) = 2·t²/2+20(s + 1)s²(s + 1)C(s)−

อินพุตพาราโบลา r(t) = 2·t²/2 ถูกป้อนเข้าวงฟีดแบ็กหนึ่งหน่วยที่มี G(s) = 20(s + 1)/(s²(s + 1)) จงหาค่าความคลาดเคลื่อนสถานะคงตัว ess

G(s) = 20(s + 1)/(s²(s + 1))r(t) = 2·t²/2
ess =
10.ค่าความคลาดเคลื่อนสถานะคงตัวและการควบคุม PIDยาก
R(s)+K = 401s(s + 8)C(s)−

ตัวควบคุมสัดส่วน K = 40 ขับโรงงาน G(s) = 1/(s(s + 8)) ในวงฟีดแบ็กหนึ่งหน่วย จงหาค่าพุ่งเกินร้อยละ %OS ของการตอบสนองต่อสัญญาณขั้นบันได

G(s) = 1/(s(s + 8))K = 40
%OS = %
11.แผนภาพโบดและขอบเสถียรภาพปานกลาง
vin​R = 3.3 kΩC = 100 nFvout​

ตัวกรองความถี่ต่ำผ่าน RCที่แสดงมีR = 3.3 kΩ, C = 100 nF ป้อนอินพุตรูปซายน์ของ 1.4 kHz จงหาอัตราขยาย |H| = |vout/vin| ใน dB

R = 3.3 kΩC = 100 nFf = 1.4 kHz
|H| = dB
12.แผนภาพโบดและขอบเสถียรภาพยาก
|H| (dB)0.10.1111010100100−80−60−40−2002040ω (rad/s)φ (°)ω (rad/s)−180−135−90−450ωc​

วงรอบป้อนกลับหนึ่งหน่วยมีฟังก์ชันถ่ายโอนวงรอบเปิด L(s) = 7/(s(s + 1)) จงหาขอบเฟส PM เป็นองศา

L(s) = 7/(s(s + 1))
PM = °

3สารละลาย

1.การแปลงลาปลาซ

คำตอบ: f(1) = 8.026

การแปลงลาปลาซ: 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² + 2s + 2 = (s + 1)² + 1 = (s + 1)² + 1²

จับคู่เศษ

F(s) = 3·(s + 1)/((s + 1)² + 1²) + (24/1)·1/((s + 1)² + 1²)

การแปลงผกผัน

f(t) = e−t·(3·cos 1t + 24·sin 1t)
f(1) = e−1·(3·cos(1·1) + 24·sin(1·1)) = 8.026

คำตอบ

f(1) = 8.026
2.การแปลงลาปลาซ

คำตอบ: g(0.1) = 1.399

การแปลงลาปลาซ: 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)

ฟังก์ชันถ่ายโอน

G(s) = Y(s)/U(s) = (s + 18)/(s² + 10s + 24)

เศษส่วนย่อย

s² + 10s + 24 = (s + 4)(s + 6)
A = 7, B = −6

ผลตอบสนองอิมพัลส์

g(t) = 7·e−4t − 6·e−6t
g(0.1) = 1.399

คำตอบ

g(0.1) = 1.399
3.แผนภาพบล็อก

คำตอบ: T(0) = 0.4068

พีชคณิตแผนภาพบล็อก: บล็อกอนุกรมคูณกัน (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)

วงใน

Gᵢ = G/(1 + G·H₁) = 8/(s + 3 + 8·4) = 8/(s + 35)

วงนอก

T(s) = K·Gᵢ/(1 + K·Gᵢ) = 24/(s + 35 + 24) = 24/(s + 59)

ผลลัพธ์

T(0) = 24/59 = 0.4068

คำตอบ

T(0) = 0.4068
4.แผนภาพบล็อก

คำตอบ: Kmax = 1224

พีชคณิตแผนภาพบล็อก: บล็อกอนุกรมคูณกัน (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)

สมการลักษณะเฉพาะ

s³ + 17s² + 72s + K = 0

ตารางเราท์

s³ | 1 72
s² | 17 K
s¹ | (17·72 − K)/17
s⁰ | K

ขีดจำกัดเสถียรภาพ

ทุกค่าในคอลัมน์แรกต้องเป็นบวก

Kmax = 17·72 = 1224

คำตอบ

Kmax = 1224
5.การตอบสนองต่อสัญญาณขั้นบันได

คำตอบ: tp = 0.3755 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 + ωₙ²

ωₙ² = 95 → ωₙ = 9.747 rad/s
2ζωₙ = 10 → ζ = 10/(2·9.747) = 0.513

ผลลัพธ์

ωd = 9.747·√(1 − 0.513²) = 8.367 rad/s → tp = π/ωd = 0.3755 s

คำตอบ

tp = 0.3755 s
6.การตอบสนองต่อสัญญาณขั้นบันได

คำตอบ: tp = 0.3293 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)

ความถี่ธรรมชาติและอัตราส่วนความหน่วง

ωₙ = √(k/m) = √(500/5) = 10 rad/s
ζ = c/(2√(km)) = 30/(2·√(500·5)) = 0.3

ผลลัพธ์

ωd = 9.539 rad/s
tp = π/ωd = 0.3293 s

คำตอบ

tp = 0.3293 s
7.เสถียรภาพ: เกณฑ์ราอ์ท–เฮอร์วิตซ์และรูตโลกัส

คำตอบ: Kmax = 10

เกณฑ์ Routh–Hurwitz: ระบบเชิงเส้นเสถียรเมื่อรากทุกตัวของพหุนามลักษณะเฉพาะอยู่ในครึ่งซ้ายของระนาบ s เขียนสัมประสิทธิ์สลับกันลงในสองแถวแรกของตาราง Routh แต่ละช่องของแถวใหม่มาจากสองแถวด้านบน สำหรับกำลังสามเชิงลูกบาศก์ a₃s³ + a₂s² + a₁s + a₀ ช่อง s¹ คือ b₁ = (a₂·a₁ − a₃·a₀)/a₂ จำนวนการเปลี่ยนเครื่องหมายในคอลัมน์แรกเท่ากับจำนวนรากในครึ่งระนาบขวา เมื่อมีอัตราขยาย K ในสัมประสิทธิ์ ทุกช่องในคอลัมน์แรกต้องเป็นบวก: อัตราขยายที่ทำให้ช่อง s¹ เป็นศูนย์คือขีดจำกัดเสถียรภาพ และแถว s² ให้ความถี่ของการสั่นที่คงอยู่

พหุนามลักษณะเฉพาะ

s(s + 1)(s + 4) + K(s + 7) = s³ + 5s² + (4 + K)s + 7K

ตารางเราท์

s³: 1, 4 + K
s²: 5, 7K
s¹: (5·(4 + K) − 7K)/5
s⁰: 7K

เงื่อนไข

20 − 2K > 0
K < 20/2 = 10

คำตอบ

Kmax = 10

ถ้าศูนย์อยู่ทางขวาของ s = −5 วงจะคงเสถียรสำหรับทุก K > 0

8.เสถียรภาพ: เกณฑ์ราอ์ท–เฮอร์วิตซ์และรูตโลกัส

คำตอบ: K = 9

root locus: เมื่ออัตราขยาย 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)| ให้อัตราขยาย ณ จุดใดก็ได้ของโลคัส

จุดแยกออก

dK/ds = 0 สำหรับ K = −(s + 5)(s + 11) ให้ 2s + 16 = 0

sb = −8

เงื่อนไขขนาด

K = |−8 + 5|·|−8 + 11| = 3·3 = 9

คำตอบ

K = 9

สำหรับอัตราขยายนี้วงปิดหน่วงวิกฤต (ζ = 1)

9.ค่าความคลาดเคลื่อนสถานะคงตัวและการควบคุม PID

คำตอบ: ess = 0.1

ค่าความคลาดเคลื่อนสถานะคงตัว: ในวงฟีดแบ็กหนึ่งหน่วยที่มีฟังก์ชันถ่ายโอนวงเปิด 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 ตามขั้นได้อย่างแม่นยำ และตามแรมป์ด้วยค่าความคลาดเคลื่อนจำกัด

ค่าคงที่ความเร่ง

Ka = 20·1/1 = 20

ค่าความคลาดเคลื่อนของพาราโบลา

ess = 2/20 = 0.1

คำตอบ

ess = 0.1
10.ค่าความคลาดเคลื่อนสถานะคงตัวและการควบคุม PID

คำตอบ: %OS = 7.691 %

การควบคุม 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(s + 8) + 40·1 = s² + 8s + 40

ωₙ และ ζ

ωₙ = √40 = 6.325 rad/s
ζ = 8/(2·6.325) = 0.6325

การเกินเลย

%OS = 100·e−0.6325π/√(1 − 0.6325²) = 7.691 %

คำตอบ

%OS = 7.691 %
11.แผนภาพโบดและขอบเสถียรภาพ

คำตอบ: |H| = −9.743 dB

การตอบสนองเชิงความถี่ (frequency response): แทน s = jω ในฟังก์ชันถ่ายโอน ขนาดเป็นเดซิเบลคือ 20·log₁₀|H(jω)| และเฟสคือผลรวมของมุมของตัวประกอบเศษลบด้วยมุมของตัวประกอบส่วน ตัวกรอง RC หรือ RL อันดับหนึ่งมีความถี่มุม (−3 dB) ที่ ωc = 1/(RC) หรือ R/L และ fc = ωc/(2π) ตัวกรองความถี่ต่ำผ่านให้ |H| = 1/√(1 + (ω/ωc)²) และ φ = −atan(ω/ωc) ในแผนภาพโบดแบบเส้นกำกับ (เส้นตรง) โพล (pole) แต่ละตัวหักขนาดลง 20 dB/decade ที่ความถี่มุมของมัน ศูนย์ (zero) แต่ละตัวหักขึ้น 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 วงรอบปิดเสถียรเมื่อขอบทั้งสองเป็นบวก

ความถี่มุม

ωc = 1/(RC) = 1/(3300·1 × 10⁻⁷) = 3030 rad/s
fc = 482.3 Hz → f/fc = 1400/482.3 = 2.903

อัตราขยาย

|H| = 1/√(1 + 2.903²) = 0.3257
20·log₁₀(0.3257) = −9.743 dB

คำตอบ

|H| = −9.743 dB
12.แผนภาพโบดและขอบเสถียรภาพ

คำตอบ: PM = 21.39 °

การตอบสนองเชิงความถี่ (frequency response): แทน s = jω ในฟังก์ชันถ่ายโอน ขนาดเป็นเดซิเบลคือ 20·log₁₀|H(jω)| และเฟสคือผลรวมของมุมของตัวประกอบเศษลบด้วยมุมของตัวประกอบส่วน ตัวกรอง RC หรือ RL อันดับหนึ่งมีความถี่มุม (−3 dB) ที่ ωc = 1/(RC) หรือ R/L และ fc = ωc/(2π) ตัวกรองความถี่ต่ำผ่านให้ |H| = 1/√(1 + (ω/ωc)²) และ φ = −atan(ω/ωc) ในแผนภาพโบดแบบเส้นกำกับ (เส้นตรง) โพล (pole) แต่ละตัวหักขนาดลง 20 dB/decade ที่ความถี่มุมของมัน ศูนย์ (zero) แต่ละตัวหักขึ้น 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 วงรอบปิดเสถียรเมื่อขอบทั้งสองเป็นบวก

จุดตัดอัตราขยาย

ω⁴ + 1ω² − 49 = 0
ω² = (−1 + √(1 + 4·49))/2 = 6.518
ωc = 2.553 rad/s

ขอบเฟส

∠L(jωc) = −90° − atan(2.553/1) = −158.6°
PM = 180° + ∠L = 21.39°

คำตอบ

PM = 21.39°

ExamProblems · examproblems.com/courses/control-engineering-1