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

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

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

ในหน้าต่างพิมพ์ ให้เลือก “บันทึกเป็น PDF” เพื่อเก็บสำเนาไว้ในอุปกรณ์ของคุณ

1สูตรและวิธีทำ

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

เสถียรภาพ: รากทั้งหมดของพหุนามลักษณะเฉพาะอยู่ในครึ่งระนาบซ้าย
ตารางราอ์ท ค่าในแถว s¹ ของระบบอันดับสาม: b₁ = (a₂a₁ − a₃a₀)/a₂
รากในครึ่งระนาบขวา: จำนวนครั้งที่เครื่องหมายเปลี่ยนในคอลัมน์แรก
ระบบอันดับสาม s³ + a₂s² + a₁s + a₀: เสถียรเมื่อ a ทุกตัว > 0 และ a₂·a₁ > a₀
K/(s(s + a)(s + b)): K_max = 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
ค่าคงที่ความคลาดเคลื่อน: K_p = lim G(s), K_v = lim s·G(s), K_a = lim s²·G(s) (s → 0)
ขั้น A: e_ss = A/(1 + K_p)
แรมป์ A·t: e_ss = A/K_v
พาราโบลา A·t²/2: e_ss = A/K_a
ตัวควบคุม PID: C(s) = K_p + K_i/s + K_d·s
เป้าหมายอันดับสอง: s² + 2ζω_ns + ω_n²
การเกินเลยและเวลา: %OS = 100·e^{−ζπ/√(1−ζ²)}, t_p = π/ω_d, t_s ≈ 4/(ζω_n)

วิธีการ

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

2โจทย์ฝึก

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

จงหา f(t) สำหรับ F(s) = (5s + 35)/(s² + 6s + 9) และหาค่า f(0.1)

F(s) = (5s + 35)/(s² + 6s + 9)t = 0.1 s

f(0.1) =

2.การแปลงลาปลาซยาก

จงแก้ y″ + 11y′ + 30y = 0 โดย y(0) = 3 และ y′(0) = 3 ด้วยการแปลงลาปลาซ และหา y(0.1)

y″ + 11y′ + 30y = 0y(0) = 3y′(0) = 3t = 0.1 s

y(0.1) =

3.แผนภาพบล็อกปานกลาง

วงป้อนกลับลบหนึ่งหน่วยมี G(s) = 11(s + 2)/((s + 6)(s + 8)) ป้อนฟังก์ชันขั้นขนาด 5 ที่สัญญาณอ้างอิง จงหาค่าผิดพลาดสถานะคงตัว e_ss

G(s) = 11(s + 2)/((s + 6)(s + 8))R = 5

e_ss =

4.แผนภาพบล็อกยาก

ตัวควบคุมสัดส่วน K = 27 ขับพืช G(s) = 8/(s + 9) ในวงป้อนกลับลบหนึ่งหน่วย สิ่งรบกวน D เข้าที่อินพุตของพืช (บวกเข้ากับเอาต์พุตของตัวควบคุม) ป้อนฟังก์ชันขั้นอ้างอิงขนาด 10 และฟังก์ชันขั้นสิ่งรบกวนขนาด 1 พร้อมกัน จงหาเอาต์พุตสถานะคงตัว c_ss

K = 27G(s) = 8/(s + 9)D = 1R = 10

c_ss =

5.การตอบสนองต่อสัญญาณขั้นบันไดปานกลาง

ระบบหนึ่งมีฟังก์ชันถ่ายโอน G(s) = 76/(s² + 8s + 38) และรับอินพุตขั้นบันไดหนึ่งหน่วย จงหาเวลาเข้าสู่ค่าคงตัวที่เกณฑ์ 2 % t_s ของการตอบสนองต่อสัญญาณขั้นบันได โดยใช้ t_s ≈ 4/(ζωₙ)

G(s) = 76/(s² + 8s + 38)

t_s = s

6.การตอบสนองต่อสัญญาณขั้นบันไดยาก

การตอบสนองต่อสัญญาณขั้นบันไดของระบบอันดับสองที่ไม่มีศูนย์พุ่งเกิน 8 % และยอดแรกเกิดที่ t = 0.5 s จงหาความถี่ธรรมชาติ ω_n

%OS = 8 %t_p = 0.5 s

ω_n = rad/s

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

วงในแผนภาพทำงานที่ K = 65 และ L(s) = K/(s(s + 5)(s + 8)) เพิ่มอัตราขยายได้กี่เท่า m ก่อนที่วงปิดจะไม่เสถียร

L(s) = K/(s(s + 5)(s + 8))K = 65

m =

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

root locus ของ L(s) = K/((s + 1)(s + 5)(s + 6)) ถูกร่างไว้ จงหาความถี่ ω ที่โลคัสตัดแกนจินตภาพ

L(s) = K/((s + 1)(s + 5)(s + 6))

ω = rad/s

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

วงฟีดแบ็กหนึ่งหน่วยที่มี G(s) = K/(s(s + 2)) ต้องตามแรมป์ r(t) = 10t ด้วยค่าความคลาดเคลื่อนสถานะคงตัว 0.2 จงหาอัตราขยาย K

G(s) = K/(s(s + 2))r(t) = 10te_ss = 0.2

K =

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

ตัวควบคุม PI C(s) = K_p + K_i/s ขับโรงงาน G(s) = 4/(s + 4) ในวงฟีดแบ็กหนึ่งหน่วย วงปิดต้องมี ωₙ = 4 rad/s และ ζ = 0.6 จงหาอัตราขยายสัดส่วน K_p

G(s) = 4/(s + 4)ωₙ = 4 rad/sζ = 0.6

K_p =

11.แผนภาพโบดและขอบเสถียรภาพปานกลาง

ในตัวกรองแถบผ่าน RLC อนุกรมที่แสดง เอาต์พุตอยู่คร่อม R = 10 Ω โดย L = 5 mH และ C = 0.1 µF จงหาแบนด์วิดท์ −3 dB BW ใน Hz

R = 10 ΩL = 5 mHC = 0.1 µF

BW = Hz

12.แผนภาพโบดและขอบเสถียรภาพยาก

วงรอบป้อนกลับหนึ่งหน่วยมีฟังก์ชันถ่ายโอนวงรอบเปิด L(s) = 132/(s(s + 5)(s + 7)) จงหาความถี่จุดตัดเฟส ω₁₈₀ ที่ ∠L(jω) = −180°

L(s) = 132/(s(s + 5)(s + 7))

ω₁₈₀ = rad/s

3สารละลาย

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

คำตอบ: f(0.1) = 5.186

การแปลงลาปลาซ: 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(∞) = lim_s→0 s·F(s) ใช้ได้เฉพาะเมื่อโพลทั้งหมดของ s·F(s) อยู่ในระนาบครึ่งซ้าย ทฤษฎีบทค่าเริ่มต้น: f(0⁺) = lim_s→∞ s·F(s)

แยกตัวประกอบส่วน

s² + 6s + 9 = (s + 3)²

เศษส่วนย่อย

5s + 35 = 5(s + 3) + 20

F(s) = 5/(s + 3) + 20/(s + 3)²

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

f(t) = 5·e^−3t + 20·t·e^−3t

f(0.1) = (5 + 20·0.1)·e^−3·0.1 = 5.186

คำตอบ

f(0.1) = 5.186

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

คำตอบ: y(0.1) = 2.859

แปลงสมการ

s²Y − 3s − 3 + 11(sY − 3) + 30Y = 0

Y(s) = (3s + 36)/((s + 5)(s + 6))

เศษส่วนย่อย

A = 21/1 = 21

B = 18/(−1) = −18

ผลลัพธ์

y(t) = 21·e^−5t − 18·e^−6t

y(0.1) = 2.859

คำตอบ

y(0.1) = 2.859

ตรวจสอบ: y(0) = A + B = 3

3.แผนภาพบล็อก

คำตอบ: e_ss = 3.429

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

ค่าคงที่ความผิดพลาดตำแหน่ง

Kₚ = G(0) = 11·2/(6·8) = 0.4583

ค่าผิดพลาดสถานะคงตัว

e_ss = A/(1 + Kₚ) = 5/(1 + 0.4583) = 3.429

คำตอบ

e_ss = 3.429

4.แผนภาพบล็อก

คำตอบ: c_ss = 9.636

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

C/D = G/(1 + K·G) = 8/(s + 9 + 216)

C/R = K·G/(1 + K·G) = 216/(s + 9 + 216)

สภาวะคงตัว

c_ss = 10·216/225 + 1·8/225 = 9.636

คำตอบ

c_ss = 9.636

K ที่มากขึ้นลดผลของสิ่งรบกวน แต่ไม่ถึงศูนย์ด้วยตัวควบคุมสัดส่วน

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

คำตอบ: t_s = 1 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−ζ²)} เวลาถึงจุดสูงสุด t_p = π/ω_d และเวลาเข้าสู่ค่าคงตัวที่เกณฑ์ 2 % t_s ≈ 4/(ζωₙ) ให้เทียบตัวส่วนกับ s² + 2ζωₙs + ωₙ² เพื่ออ่านค่า ωₙ และ ζ สำหรับระบบมวล–สปริง–ตัวหน่วง ωₙ = √(k/m) และ ζ = c/(2√(km)) สำหรับวงจร RLC อนุกรม ωₙ = 1/√(LC) และ ζ = (R/2)·√(C/L)

เทียบกับรูปมาตรฐาน

ตัวส่วนต้องตรงกับ s² + 2ζωₙs + ωₙ²

ωₙ² = 38 → ωₙ = 6.164 rad/s

2ζωₙ = 8 → ζ = 8/(2·6.164) = 0.6489

ผลลัพธ์

ζωₙ = 4 → t_s ≈ 4/4 = 1 s

คำตอบ

t_s = 1 s

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

คำตอบ: ω_n = 8.062 rad/s

อัตราส่วนการหน่วง

ζ = −ln(0.08)/√(π² + ln²(0.08)) = 2.526/√(9.87 + 6.379) = 0.6266

ความถี่ธรรมชาติ

ω_d = π/t_p = π/0.5 = 6.283 rad/s

ωₙ = ω_d/√(1 − ζ²) = 6.283/0.7794 = 8.062 rad/s

คำตอบ

ω_n = 8.062 rad/s

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

คำตอบ: m = 8

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

ตารางเราท์

s³ + 13s² + 40s + K = 0

s³: 1, 40

s²: 13, K

s¹: (13·40 − K)/13

s⁰: K

อัตราขยายสูงสุดที่เสถียร

K_max = 13·40 = 520

ตัวประกอบ

m = 520/65 = 8

คำตอบ

m = 8

ตัวประกอบนี้คือขอบอัตราขยายเป็นอัตราส่วน: 18.06 dB

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

คำตอบ: ω = 6.403 rad/s

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

อัตราขยายวิกฤต

s³ + 12s² + 41s + 30 + K = 0

s¹: (12·41 − 30 − K)/12 = 0 → K = 462

ความถี่ที่ตัด

12s² + 30 + 462 = 12s² + 492 = 0 → s² = −41

ω = √41 = 6.403 rad/s

คำตอบ

ω = 6.403 rad/s

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

คำตอบ: K = 100

ค่าความคลาดเคลื่อนสถานะคงตัว: ในวงฟีดแบ็กหนึ่งหน่วยที่มีฟังก์ชันถ่ายโอนวงเปิด G(s) ค่าความคลาดเคลื่อนคือ E(s) = R(s)/(1 + G(s)) และโดยทฤษฎีบทค่าสุดท้าย e_ss = lim_s→0 s·E(s) ชนิดของระบบคือจำนวนโพลของ G(s) ที่ s = 0 ค่าคงที่ความคลาดเคลื่อนคือ K_p = lim G(s), K_v = lim s·G(s) และ K_a = lim s²·G(s) ทั้งหมดสำหรับ s → 0 ขั้นขนาด A ให้ e_ss = A/(1 + K_p) แรมป์ A·t ให้ A/K_v และพาราโบลา A·t²/2 ให้ A/K_a วงชนิด 0 ตามขั้นด้วยค่าความคลาดเคลื่อนจำกัด วงชนิด 1 ตามขั้นได้อย่างแม่นยำ และตามแรมป์ด้วยค่าความคลาดเคลื่อนจำกัด

ค่าคงที่ความเร็วที่ต้องการ

K_v = A/e_ss = 10/0.2 = 50

อัตราขยาย

K = 2·K_v = 100

คำตอบ

K = 100

อัตราขยายที่มากทำให้วงนี้หน่วงน้อยมาก: ζ = 0.1 ความขัดแย้งนั้นคือเหตุผลที่ใช้ตัวชดเชย

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

คำตอบ: K_p = 0.2

การควบคุม P, PI และ PID: ตัวควบคุม C(s) = K_p + K_i/s + K_d·s กระทำต่อค่าความคลาดเคลื่อน อัตราขยายสัดส่วนทำให้วงเร็วขึ้นและลดค่าความคลาดเคลื่อนของโรงงานชนิด 0 แต่ไม่เคยกำจัดมัน พจน์อินทิกรัลเพิ่มโพลที่ s = 0: ชนิดของระบบเพิ่มขึ้นหนึ่ง และตามขั้นได้โดยไม่มีค่าความคลาดเคลื่อน พจน์อนุพันธ์เพิ่มการหน่วง ในการเลือกอัตราขยาย เขียนพหุนามลักษณะเฉพาะวงปิดและเปรียบเทียบกับ s² + 2ζω_ns + ω_n² จากนั้น %OS = 100·e^{−ζπ/√(1−ζ²)}, t_p = π/ω_d โดย ω_d = ω_n√(1 − ζ²) และ t_s ≈ 4/(ζω_n)

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

1 + C(s)G(s) = 0 คูณด้วย s(s + 4)

s² + (4 + 4K_p)s + 4K_i

พหุนามเป้าหมาย

s² + 2·0.6·4s + 4² = s² + 4.8s + 16

จับคู่สัมประสิทธิ์

4 + 4K_p = 4.8 → K_p = 0.2

คำตอบ

K_p = 0.2

พจน์อินทิกรัลทำให้วงเป็นชนิด 1 ดังนั้นตามขั้นได้โดยไม่มีค่าความคลาดเคลื่อนสถานะคงตัว

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

คำตอบ: BW = 318.3 Hz

การตอบสนองเชิงความถี่ (frequency response): แทน s = jω ในฟังก์ชันถ่ายโอน ขนาดเป็นเดซิเบลคือ 20·log₁₀|H(jω)| และเฟสคือผลรวมของมุมของตัวประกอบเศษลบด้วยมุมของตัวประกอบส่วน ตัวกรอง RC หรือ RL อันดับหนึ่งมีความถี่มุม (−3 dB) ที่ ω_c = 1/(RC) หรือ R/L และ f_c = ω_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) จุดตัดเฟส ω₁₈₀: ∠L(jω) = −180° ขอบอัตราขยาย GM = −20·log₁₀|L(jω₁₈₀)| dB วงรอบปิดเสถียรเมื่อขอบทั้งสองเป็นบวก

ความถี่เรโซแนนซ์

ω₀ = 1/√(LC) = 1/√(0.005·1 × 10⁻⁷) = 44721 rad/s

f₀ = ω₀/(2π) = 7118 Hz

แบนด์วิดท์

BW = R/L = 10/0.005 = 2000 rad/s = 318.3 Hz

คำตอบ

BW = 318.3 Hz

12.แผนภาพโบดและขอบเสถียรภาพ

คำตอบ: ω₁₈₀ = 5.916 rad/s

จุดตัดเฟส

(ω/5)·(ω/7) = 1 → ω₁₈₀ = √(5·7) = 5.916 rad/s

คำตอบ

ω₁₈₀ = 5.916 rad/s

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