วิศวกรรมระบบควบคุม 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.การแปลงลาปลาซปานกลาง
t (s)f(t)t = 0.5f(0.5) = ?

จงหา f(t) สำหรับ F(s) = (4s + 14)/(s² + 2s + 1) และหาค่า f(0.5)

F(s) = (4s + 14)/(s² + 2s + 1)t = 0.5 s
f(0.5) =
2.การแปลงลาปลาซยาก
t (s)y(t)t = 0.25y(0.25) = ?

ระบบเป็นไปตาม y″ + 9y′ + 18y = 72·u(t) โดย y(0) = 3 และ y′(0) = 0 ที่ t = 0 มีอินพุตฟังก์ชันขั้น u(t) = 4 ป้อนเข้า จงใช้การแปลงลาปลาซหา y(0.25)

y″ + 9y′ + 18y = 72·u(t)y(0) = 3y′(0) = 0u(t) = 4t = 0.25 s
y(0.25) =
3.แผนภาพบล็อกปานกลาง
R(s)+17K+10s + 4G(s)C(s)3H₁−−

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

G(s) = 10/(s + 4)H₁ = 3K = 17
T(0) =
4.แผนภาพบล็อกยาก
R(s)+Kตัวควบคุม+D(s)+10s + 6พืชC(s)−

ตัวควบคุมสัดส่วน K ขับพืช G(s) = 10/(s + 6) ในวงป้อนกลับลบหนึ่งหน่วย ฟังก์ชันขั้นสิ่งรบกวนขนาด 10 เข้าที่อินพุตของพืชขณะที่ R = 0 จงหาอัตราขยายที่น้อยที่สุด K ที่ทำให้เอาต์พุตสถานะคงตัวไม่เกิน 0.1667

G(s) = 10/(s + 6)D = 10css ≤ 0.1667
K =
5.การตอบสนองต่อสัญญาณขั้นบันไดปานกลาง
t (s)y(t)1

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

G(s) = 63/(s² + 10s + 63)
ζ =
6.การตอบสนองต่อสัญญาณขั้นบันไดยาก
R(s)+Ks(s + 2)G(s)C(s)−

ลูปป้อนกลับลบแบบหนึ่งหน่วยมี G(s) = K/(s(s + 2)) จงหาอัตราขยาย K ที่ทำให้การตอบสนองต่อสัญญาณขั้นบันไดพุ่งเกิน 30 %

G(s) = K/(s(s + 2))%OS = 30 %
K =
7.เสถียรภาพ: เกณฑ์ราอ์ท–เฮอร์วิตซ์และรูตโลกัสปานกลาง
R(s)+K = 1601s(s + 6)(s + 10)C(s)−

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

L(s) = K/(s(s + 6)(s + 10))K = 160
m =
8.เสถียรภาพ: เกณฑ์ราอ์ท–เฮอร์วิตซ์และรูตโลกัสยาก
ReImζ = 0.8−1−2−3−4−5−6j1j2

root locus ของ L(s) = K/(s(s + 5)) แสดงไว้ จงหาอัตราขยาย K ที่ทำให้โพลวงปิดมีอัตราส่วนการหน่วง ζ = 0.8

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

วงฟีดแบ็กหนึ่งหน่วยมี G(s) = K/((s + 1)(s + 2)) จงหาอัตราขยาย K ที่ให้ค่าความคลาดเคลื่อนสถานะคงตัว 2 % สำหรับอินพุตขั้น

G(s) = K/((s + 1)(s + 2))ess = 2 %
K =
10.ค่าความคลาดเคลื่อนสถานะคงตัวและการควบคุม PIDยาก
R(s)+Kp​ + Ki​/s10s + 5C(s)−

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

G(s) = 10/(s + 5)ωₙ = 8 rad/sζ = 1
Kp =
11.แผนภาพโบดและขอบเสถียรภาพปานกลาง
|H| (dB)0.10.1111010100100−200204060ω (rad/s)φ (°)ω (rad/s)−90−450ω = 5

จงหาเฟส ∠G(jω) ของ G(s) = 20(s + 7)/(s(s + 10)) ที่ ω = 5 rad/s เป็นองศา

G(s) = 20(s + 7)/(s(s + 10))ω = 5 rad/s
∠G(jω) = °
12.แผนภาพโบดและขอบเสถียรภาพยาก
|H| (dB)0.10.1111010100100−80−60−40−2002040ω (rad/s)φ (°)ω (rad/s)−270−180−900ω180​ = ?

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

L(s) = 98/(s(s + 6)(s + 7))
ω180 = rad/s

3สารละลาย

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

คำตอบ: f(0.5) = 5.459

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

เศษส่วนย่อย

4s + 14 = 4(s + 1) + 10
F(s) = 4/(s + 1) + 10/(s + 1)²

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

f(t) = 4·e−t + 10·t·e−t
f(0.5) = (4 + 10·0.5)·e−0.5 = 5.459

คำตอบ

f(0.5) = 5.459
2.การแปลงลาปลาซ

คำตอบ: y(0.25) = 6.619

การแปลงลาปลาซ: 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² + 9s + 18)·Y(s) = 3s + 27 + 288/s
Y(s) = (3s² + 27s + 288)/(s(s + 3)(s + 6))

เศษส่วนย่อย

A = 288/(3·6) = 16
B = 234/((−3)·3) = −26
C = 234/((−6)·(−3)) = 13

ผลลัพธ์

y(t) = 16 − 26·e−3t + 13·e−6t
y(0.25) = 6.619

คำตอบ

y(0.25) = 6.619

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

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

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

พีชคณิตแผนภาพบล็อก: บล็อกอนุกรมคูณกัน (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₁) = 10/(s + 4 + 10·3) = 10/(s + 34)

วงนอก

T(s) = K·Gᵢ/(1 + K·Gᵢ) = 170/(s + 34 + 170) = 170/(s + 204)

ผลลัพธ์

T(0) = 170/204 = 0.8333

คำตอบ

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

คำตอบ: K = 59.4

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

แก้หา K

6 + 10K = 10·10/0.1667 = 600
K = (600 − 6)/10 = 59.4

คำตอบ

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

คำตอบ: ζ = 0.6299

การตอบสนองต่อสัญญาณขั้นบันได: ระบบอันดับหนึ่ง 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 + ωₙ²

ωₙ² = 63 → ωₙ = 7.937 rad/s
2ζωₙ = 10 → ζ = 10/(2·7.937) = 0.6299

ผลลัพธ์

ζ = 0.6299

คำตอบ

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

คำตอบ: K = 7.809

การตอบสนองต่อสัญญาณขั้นบันได: ระบบอันดับหนึ่ง 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.3)/√(π² + ln²(0.3)) = 0.3579

อัตราขยาย

ωₙ = 2/(2ζ) = 2.794 rad/s
K = ωₙ² = 7.809

คำตอบ

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

คำตอบ: m = 6

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

ตารางเราท์

s³ + 16s² + 60s + K = 0
s³: 1, 60
s²: 16, K
s¹: (16·60 − K)/16
s⁰: K

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

Kmax = 16·60 = 960

ตัวประกอบ

m = 960/160 = 6

คำตอบ

m = 6

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

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

คำตอบ: K = 9.766

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² + 5s + K = s² + 2ζωₙs + ωₙ²
2·0.8·ωₙ = 5 → ωₙ = 3.125 rad/s

อัตราขยาย

K = ωₙ² = 3.125² = 9.766

คำตอบ

K = 9.766

โพลวงปิดอยู่ที่เส้น ζ ตัดกิ่งแนวตั้งที่ s = −2.5

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

คำตอบ: K = 98

ค่าความคลาดเคลื่อนสถานะคงตัว: ในวงฟีดแบ็กหนึ่งหน่วยที่มีฟังก์ชันถ่ายโอนวงเปิด 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/2 = 49 → K = 98

คำตอบ

K = 98

ตรวจสอบเสถียรภาพวงปิดด้วย: สำหรับวงอันดับสองนี้ K > 0 ใด ๆ เสถียร

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

คำตอบ: Kp = 1.1

การควบคุม 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)

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

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

s² + (5 + 10Kp)s + 10Ki

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

s² + 2·1·8s + 8² = s² + 16s + 64

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

5 + 10Kp = 16 → Kp = 1.1

คำตอบ

Kp = 1.1

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

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

คำตอบ: ∠G(jω) = −81.03 °

การตอบสนองเชิงความถี่ (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 วงรอบปิดเสถียรเมื่อขอบทั้งสองเป็นบวก

บวกมุม

∠G = atan(5/7) − 90° − atan(5/10) = 35.54° − 90° − 26.57°

ผลลัพธ์

∠G(j5) = −81.03°

คำตอบ

∠G(jω) = −81.03°
12.แผนภาพโบดและขอบเสถียรภาพ

คำตอบ: ω180 = 6.481 rad/s

การตอบสนองเชิงความถี่ (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 วงรอบปิดเสถียรเมื่อขอบทั้งสองเป็นบวก

จุดตัดเฟส

(ω/6)·(ω/7) = 1 → ω180 = √(6·7) = 6.481 rad/s

คำตอบ

ω180 = 6.481 rad/s

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