วิศวกรรมระบบควบคุม 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.การแปลงลาปลาซปานกลาง
ReIm5

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

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

จงแก้ y″ + 7y′ + 12y = 0 โดย y(0) = 1 และ y′(0) = −12 ด้วยการแปลงลาปลาซ และหา y(0.3)

y″ + 7y′ + 12y = 0y(0) = 1y′(0) = −12t = 0.3 s
y(0.3) =
3.แผนภาพบล็อกปานกลาง
R(s)+10K+5s + 3G(s)C(s)5H₁−−

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

G(s) = 5/(s + 3)H₁ = 5K = 10
T(0) =
4.แผนภาพบล็อกยาก
R(s)+74s(s + 6)G(s)C(s)0.5H−

พืช G(s) = 74/(s(s + 6)) อยู่ในวงป้อนกลับลบที่ตัวตรวจจับมีอัตราขยาย H = 0.5 โพลวงปิดเป็นจำนวนเชิงซ้อน จงหาส่วนจริง σ ของพวกมัน

G(s) = 74/(s(s + 6))H = 0.5
σ =
5.การตอบสนองต่อสัญญาณขั้นบันไดปานกลาง
t (s)y(t)1?

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

G(s) = 26/(s² + 7s + 26)
%OS = %
6.การตอบสนองต่อสัญญาณขั้นบันไดยาก
t (s)y(t)1tp​ = 0.5 s25 %

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

%OS = 25 %tp = 0.5 s
ωn = rad/s
7.เสถียรภาพ: เกณฑ์ราอ์ท–เฮอร์วิตซ์และรูตโลกัสปานกลาง
R(s)+K1(s + 1)(s + 2)(s + 3)C(s)−

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

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

root locus ของ L(s) = K/(s(s + 5)(s + 10)) ออกจากแกนจริงระหว่าง s = 0 และ s = −5 จงหาจุดแยกออก sb

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

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

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

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

G(s) = 1/(s + 4)ωₙ = 8 rad/sζ = 1
Kp =
11.แผนภาพโบดและขอบเสถียรภาพปานกลาง
vin​R = 3.3 kΩC = 1 nFvout​

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

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

วงรอบป้อนกลับหนึ่งหน่วยมีฟังก์ชันถ่ายโอนวงรอบเปิด L(s) = K/(s(s + 1)(s + 7)) จงหาอัตราขยาย K ที่ให้ขอบอัตราขยาย 8 dB

L(s) = K/(s(s + 1)(s + 7))GM = 8 dB
K =

3สารละลาย

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

คำตอบ: f(0.75) = −0.4063

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

จับคู่เศษ

F(s) = 4·(s + 2)/((s + 2)² + 6²) + (6/6)·6/((s + 2)² + 6²)

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

f(t) = e−2t·(4·cos 6t + sin 6t)
f(0.75) = e−2·0.75·(4·cos(6·0.75) + sin(6·0.75)) = −0.4063

คำตอบ

f(0.75) = −0.4063
2.การแปลงลาปลาซ

คำตอบ: y(0.3) = −0.5418

การแปลงลาปลาซ: 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²Y − 1s − (−12) + 7(sY − 1) + 12Y = 0
Y(s) = (s − 5)/((s + 3)(s + 4))

เศษส่วนย่อย

A = −8/1 = −8
B = −9/(−1) = 9

ผลลัพธ์

y(t) = −8·e−3t + 9·e−4t
y(0.3) = −0.5418

คำตอบ

y(0.3) = −0.5418

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

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

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

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

วงนอก

T(s) = K·Gᵢ/(1 + K·Gᵢ) = 50/(s + 28 + 50) = 50/(s + 78)

ผลลัพธ์

T(0) = 50/78 = 0.641

คำตอบ

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

คำตอบ: σ = −3

พีชคณิตแผนภาพบล็อก: บล็อกอนุกรมคูณกัน (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(s + 6) + 74·0.5 = 0 → s² + 6s + 37 = 0

ราก

s = −3 ± j·√(37 − 9) = −3 ± j5.292
σ = −3

คำตอบ

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

คำตอบ: %OS = 5.155 %

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

ωₙ² = 26 → ωₙ = 5.099 rad/s
2ζωₙ = 7 → ζ = 7/(2·5.099) = 0.6864

ผลลัพธ์

%OS = 100·e−ζπ/√(1−ζ²) = 100·e−2.965 = 5.155 %

คำตอบ

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

คำตอบ: ωn = 6.868 rad/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)

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

ζ = −ln(0.25)/√(π² + ln²(0.25)) = 1.386/√(9.87 + 1.922) = 0.4037

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

ωd = π/tp = π/0.5 = 6.283 rad/s
ωₙ = ωd/√(1 − ζ²) = 6.283/0.9149 = 6.868 rad/s

คำตอบ

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

คำตอบ: Kmax = 60

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

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

s³ + 6s² + 11s + 6 + K = 0

ตารางเราท์

s³: 1, 11
s²: 6, 6 + K
s¹: (6·11 − 6 − K)/6
s⁰: 6 + K

เงื่อนไข

ช่อง s¹ ต้องคงเป็นบวก

K < 6·11 − 6 = 60

คำตอบ

Kmax = 60

ที่ขีดจำกัด วงสั่นที่ ω = √11 = 3.317 rad/s

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

คำตอบ: sb = −2.113

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

อัตราขยายตามแกนจริง

K(s) = −s(s + 5)(s + 10) = −(s³ + 15s² + 50s)

จุดนิ่ง

dK/ds = −(3s² + 30s + 50) = 0
s = (−30 ± √(900 − 600))/6 = −2.113, −7.887

เลือกจุดบนโลคัส

โลคัสอยู่บนแกนจริงระหว่าง 0 และ −5 รากอีกตัวไม่อยู่บนโลคัส

sb = −2.113

คำตอบ

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

คำตอบ: K = 180

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

อัตราขยาย

K/20 = 9 → K = 180

คำตอบ

K = 180

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

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

คำตอบ: Kp = 12

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

s² + (4 + 1Kp)s + 1Ki

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

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

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

4 + 1Kp = 16 → Kp = 12

คำตอบ

Kp = 12

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

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⁻⁹) = 303030 rad/s
fc = 48229 Hz → f/fc = 140000/48229 = 2.903

อัตราขยาย

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

คำตอบ

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

คำตอบ: K = 22.29

การตอบสนองเชิงความถี่ (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)·(ω/7) = 1 → ω180 = √(1·7) = 2.646 rad/s

ขนาดที่จุดตัดเฟส

|L(jω180)| = K/(ω·√(ω² + 1)·√(ω² + 49)) = K/(1·7·(1 + 7)) = K/56

อัตราขยายที่ต้องการ

56/K = 108/20 = 2.512
K = 56/2.512 = 22.29

คำตอบ

K = 22.29

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