Kỹ thuật điều khiển 1: gói ôn nước rút

Bảng công thức, các bước phương pháp chính cho từng chủ đề và 12 bài tập, kèm lời giải đầy đủ ở cuối.

Bài mớiQuay lại môn học

Trong cửa sổ in, hãy chọn “Lưu dưới dạng PDF” để giữ một bản trên thiết bị của bạn.

1Công thức và phương pháp

Ổn định: Routh–Hurwitz và quỹ đạo nghiệm

Ổn định
mọi nghiệm của đa thức đặc trưng ở nửa mặt phẳng trái
bảng Routh, phần tử s¹ của đa thức bậc ba
b1 = (a2a1 − a3a0)/a2
Nghiệm ở nửa mặt phẳng phải
số lần đổi dấu trong cột đầu tiên
Đa thức bậc ba s³ + a₂s² + a₁s + a₀
ổn định khi mọi a > 0 và a₂·a₁ > a₀
K/(s(s + a)(s + b))
Kmax = ab(a + b), ω = √(ab)
Phương trình đặc trưng
1 + K·G(s) = 0
Đường tiệm cận
σa = (Σp − Σz)/(n − m), θ = (2q + 1)·180°/(n − m)
Điểm tách và điểm nhập
dK/ds = 0
Điều kiện biên độ
K = 1/|G(s)|
Góc xuất phát
θd = 180° − Σ∠(from other poles) + Σ∠(from zeros)

Phương pháp

  1. Viết đa thức đặc trưng vòng kín: mẫu số của 1 + K·G(s) = 0, sau khi nhân ra.
  2. Đặt các hệ số xen kẽ vào hai hàng đầu tiên của bảng Routh và lập các hàng còn lại.
  3. Đếm số lần đổi dấu trong cột đầu tiên: mỗi lần là một cực ở nửa mặt phẳng phải. Khi K có trong bảng, hãy yêu cầu mọi phần tử của cột đầu tiên giữ dương.
  4. Với quỹ đạo nghiệm, đánh dấu các cực và điểm không, các đoạn trên trục thực, các đường tiệm cận từ trọng tâm và các điểm tách/nhập nơi dK/ds = 0. Dùng K = 1/|G(s)| cho độ lợi tại một điểm.

Sai số xác lập và điều khiển PID

Loại của hệ
số cực của G(s) tại s = 0
Hằng số sai số
Kp = lim G(s), Kv = lim s·G(s), Ka = lim s²·G(s) (s → 0)
Bước nhảy A
ess = A/(1 + Kp)
Tín hiệu dốc A·t
ess = A/Kv
Parabol A·t²/2
ess = A/Ka
Bộ điều khiển PID
C(s) = Kp + Ki/s + Kd·s
Mục tiêu bậc hai
s² + 2ζωns + ωn²
Độ vượt quá và thời gian
%OS = 100·e−ζπ/√(1−ζ²), tp = π/ωd, ts ≈ 4/(ζωn)

Phương pháp

  1. Đếm các cực của vòng hở C(s)·G(s) tại s = 0: đó là loại của hệ.
  2. Tính hằng số sai số khớp với đầu vào: Kp cho bước nhảy, Kv cho tín hiệu dốc, Ka cho parabol.
  3. Áp dụng ess = A/(1 + Kp), A/Kv hay A/Ka, hoặc giải nó theo độ lợi bạn cần.
  4. Với yêu cầu về độ cản hoặc tốc độ, viết đa thức đặc trưng vòng kín và so sánh nó với s² + 2ζω_n·s + ω_n².

2Bài tập

1.Biến đổi LaplaceTrung bình
ReIm−1−2−3−4−5−6−7−8−91

F(s) = (3s + 25)/(s² + 6s + 5). Biến đổi ngược của nó có dạng f(t) = A·e−t + B·e−5t. Tìm hệ số A bằng phân thức riêng.

F(s) = (3s + 25)/(s² + 6s + 5)
A =
2.Biến đổi LaplaceKhó
t (s)y(t)t = 0,1y(0,1) = ?

Một hệ thống tuân theo y″ + 9y′ + 20y = 40·u(t). Nó bắt đầu từ trạng thái nghỉ (y(0) = 0, y′(0) = 0) và đầu vào là một bước u(t) = 2 với t ≥ 0. Dùng biến đổi Laplace để tìm đầu ra y(0,1).

y″ + 9y′ + 20y = 40·u(t)u(t) = 2y(0) = 0, y′(0) = 0t = 0,1 s
y(0.1) =
3.Sơ đồ khốiTrung bình
R(s)+9(s + 5)(s + 7)(s + 8)G(s)C(s)−

Một vòng hồi tiếp âm đơn vị có G(s) = 9(s + 5)/((s + 7)(s + 8)). Một bậc thang có độ lớn 1 được đặt vào tín hiệu chuẩn. Tìm sai số xác lập ess.

G(s) = 9(s + 5)/((s + 7)(s + 8))R = 1
ess =
4.Sơ đồ khốiKhó
R(s)+52s(s + 6)G(s)C(s)2H−

Đối tượng G(s) = 52/(s(s + 6)) nằm trong vòng hồi tiếp âm với độ lợi cảm biến H = 2. Các cực vòng kín là số phức. Tìm phần thực σ của chúng.

G(s) = 52/(s(s + 6))H = 2
σ =
5.Đáp ứng bướcTrung bình
t (s)y(t)1

Một hệ có hàm truyền G(s) = 44/(s² + 10s + 44) và nhận tín hiệu vào là bước nhảy đơn vị. Các cực là s = σ ± jωd. Tìm phần thực σ.

G(s) = 44/(s² + 10s + 44)
σ =
6.Đáp ứng bướcKhó
+−5 VR = 7 ΩL = 50 mHC = 47 µFvC​

Một mạch RLC nối tiếp với R = 7 Ω, L = 50 mH và C = 47 µF được nối với nguồn bước nhảy 5 V tại t = 0. Tụ điện chưa tích điện và không có dòng điện chạy trước bước nhảy. Tìm độ vọt lố phần trăm %OS của điện áp tụ điện.

R = 7 ΩL = 50 mHC = 47 µFV = 5 V
%OS = %
7.Ổn định: Routh–Hurwitz và quỹ đạo nghiệmTrung bình
R(s)+13441s(s + 6)(s + 8)C(s)−

Vòng trong sơ đồ có phản hồi âm và L(s) = 1344/(s(s + 6)(s + 8)). Có bao nhiêu cực vòng kín N nằm ở nửa mặt phẳng phải?

L(s) = 1344/(s(s + 6)(s + 8))
N =
8.Ổn định: Routh–Hurwitz và quỹ đạo nghiệmKhó
ReIm−2−4−6−8−10j2sb​

Quỹ đạo nghiệm của L(s) = K/((s + 1)(s + 9)) được cho. Tìm độ lợi K tại điểm tách sb, nơi hai cực vòng kín gặp nhau trên trục thực.

L(s) = K/((s + 1)(s + 9))
K =
9.Sai số xác lập và điều khiển PIDTrung bình
r(t) = 4·t²/2+5(s + 4)s²(s + 1)C(s)−

Đầu vào parabol r(t) = 4·t²/2 được đưa vào một vòng phản hồi đơn vị với G(s) = 5(s + 4)/(s²(s + 1)). Tìm sai số xác lập ess.

G(s) = 5(s + 4)/(s²(s + 1))r(t) = 4·t²/2
ess =
10.Sai số xác lập và điều khiển PIDKhó
R(s)+K = 502s(s + 8)C(s)−

Một bộ điều khiển tỉ lệ K = 50 điều khiển đối tượng G(s) = 2/(s(s + 8)) trong một vòng phản hồi đơn vị. Tìm thời gian đỉnh tp của đáp ứng bước.

G(s) = 2/(s(s + 8))K = 50
tp = s
11.Biểu đồ Bode và độ dự trữTrung bình
vin​C = 22 nFR = 1 kΩvout​

bộ lọc thông cao RC trong hình có R = 1 kΩ, C = 22 nF. Một tín hiệu vào hình sin 3,6 kHz được đặt vào. Tìm độ khuếch đại |H| = |vout/vin| theo dB.

R = 1 kΩC = 22 nFf = 3,6 kHz
|H| = dB
12.Biểu đồ Bode và độ dự trữKhó
|H| (dB)0,10,1111010100100−80−60−40−2002040ω (rad/s)φ (°)ω (rad/s)−270−180−900ω180​ = ?

Một vòng hồi tiếp đơn vị có hàm truyền vòng hở L(s) = 66/(s(s + 3)(s + 7)). Tìm tần số cắt pha ω180, tại đó ∠L(jω) = −180°.

L(s) = 66/(s(s + 3)(s + 7))
ω180 = rad/s

3Dung dịch

1.Biến đổi Laplace

Đáp số: A = 5,5

Biến đổi Laplace: F(s) = ∫₀∞ f(t)·e−st dt biến đạo hàm thành phép nhân với s: ℒ{f′} = sF(s) − f(0), ℒ{f″} = s²F(s) − s·f(0) − f′(0). Các cặp chuẩn: 1 ↔ 1/s, e−at ↔ 1/(s + a), t·e−at ↔ 1/(s + a)², sin ωt ↔ ω/(s² + ω²), cos ωt ↔ s/(s² + ω²), và phép dịch s → s + σ nhân hàm theo thời gian với e−σt. Để biến đổi ngược một F(s) hữu tỉ, hãy tách thành các phân thức riêng: thặng dư của cực đơn p là A = [(s − p)·F(s)] tại s = p. Định lý giá trị cuối: f(∞) = lims→0 s·F(s), chỉ đúng khi mọi cực của s·F(s) nằm trong nửa mặt phẳng trái. Định lý giá trị đầu: f(0⁺) = lims→∞ s·F(s).

Phân tích mẫu số thành nhân tử

s² + 6s + 5 = (s + 1)(s + 5)
F(s) = A/(s + 1) + B/(s + 5)

Quy tắc che

Nhân với (s + 1) và đặt s = −1: số hạng B biến mất.

A = (3s + 25)/(s + 5) |s = −1 = 22/4 = 5,5

Đáp số

A = 5,5
(B = −2,5, nên f(t) = 5,5·e−t − 2,5·e−5t)
2.Biến đổi Laplace

Đáp số: y(0.1) = 0,2981

Biến đổi Laplace: F(s) = ∫₀∞ f(t)·e−st dt biến đạo hàm thành phép nhân với s: ℒ{f′} = sF(s) − f(0), ℒ{f″} = s²F(s) − s·f(0) − f′(0). Các cặp chuẩn: 1 ↔ 1/s, e−at ↔ 1/(s + a), t·e−at ↔ 1/(s + a)², sin ωt ↔ ω/(s² + ω²), cos ωt ↔ s/(s² + ω²), và phép dịch s → s + σ nhân hàm theo thời gian với e−σt. Để biến đổi ngược một F(s) hữu tỉ, hãy tách thành các phân thức riêng: thặng dư của cực đơn p là A = [(s − p)·F(s)] tại s = p. Định lý giá trị cuối: f(∞) = lims→0 s·F(s), chỉ đúng khi mọi cực của s·F(s) nằm trong nửa mặt phẳng trái. Định lý giá trị đầu: f(0⁺) = lims→∞ s·F(s).

Biến đổi phương trình

Y(s) = 80/(s(s + 4)(s + 5))

Phân thức riêng

A = 80/(4·5) = 4
B = 80/((−4)·1) = −20
C = 80/((−5)·(−1)) = 16

Lượng ra

y(t) = 4 − 20·e−4t + 16·e−5t
y(0,1) = 0,2981

Đáp số

y(0,1) = 0,2981
3.Sơ đồ khối

Đáp số: ess = 0,5545

Đại số sơ đồ khối: các khối nối tiếp thì nhân (G₁·G₂), các khối song song thì cộng (G₁ + G₂), và một vòng có đường thuận G và đường hồi tiếp H rút gọn thành T = G/(1 + G·H) với hồi tiếp âm và G/(1 − G·H) với hồi tiếp dương. Rút gọn các vòng bên trong trước, rồi làm dần ra ngoài. Các cực vòng kín là nghiệm của phương trình đặc trưng 1 + G(s)H(s) = 0. Với vòng ổn định, đầu vào bậc thang có độ lớn A cho đầu ra xác lập A·T(0); với hồi tiếp đơn vị, sai số xác lập đối với bậc thang là A/(1 + Kp) với Kp = G(0).

Hằng số sai số vị trí

Kₚ = G(0) = 9·5/(7·8) = 0,8036

Sai số xác lập

ess = A/(1 + Kₚ) = 1/(1 + 0,8036) = 0,5545

Đáp số

ess = 0,5545
4.Sơ đồ khối

Đáp số: σ = −3

Đại số sơ đồ khối: các khối nối tiếp thì nhân (G₁·G₂), các khối song song thì cộng (G₁ + G₂), và một vòng có đường thuận G và đường hồi tiếp H rút gọn thành T = G/(1 + G·H) với hồi tiếp âm và G/(1 − G·H) với hồi tiếp dương. Rút gọn các vòng bên trong trước, rồi làm dần ra ngoài. Các cực vòng kín là nghiệm của phương trình đặc trưng 1 + G(s)H(s) = 0. Với vòng ổn định, đầu vào bậc thang có độ lớn A cho đầu ra xác lập A·T(0); với hồi tiếp đơn vị, sai số xác lập đối với bậc thang là A/(1 + Kp) với Kp = G(0).

Phương trình đặc trưng

s(s + 6) + 52·2 = 0 → s² + 6s + 104 = 0

Các nghiệm

s = −3 ± j·√(104 − 9) = −3 ± j9,747
σ = −3

Đáp số

σ = −3
5.Đáp ứng bước

Đáp số: σ = −5

Đáp ứng bước: một hệ bậc nhất G(s) = K/(τs + 1) đáp ứng với một bước nhảy có biên độ U là y(t) = K·U·(1 − e−t/τ): nó đạt 63,2 % giá trị cuối K·U sau một hằng số thời gian τ và 98 % sau 4τ. Một hệ bậc hai G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) có các cực s = −ζωₙ ± jωd với ωd = ωₙ√(1 − ζ²). Với 0 < ζ < 1, đáp ứng bước có độ vọt lố: độ vọt lố phần trăm %OS = 100·e−ζπ/√(1−ζ²), thời gian đỉnh tp = π/ωd, và thời gian xác lập 2 % ts ≈ 4/(ζωₙ). So sánh mẫu số với s² + 2ζωₙs + ωₙ² để suy ra ωₙ và ζ. Với hệ khối lượng–lò xo–giảm chấn ωₙ = √(k/m) và ζ = c/(2√(km)); với mạch RLC nối tiếp ωₙ = 1/√(LC) và ζ = (R/2)·√(C/L).

So sánh với dạng chuẩn

Mẫu số phải khớp với s² + 2ζωₙs + ωₙ².

ωₙ² = 44 → ωₙ = 6,633 rad/s
2ζωₙ = 10 → ζ = 10/(2·6,633) = 0,7538

Kết quả

σ = −ζωₙ = −10/2 = −5

Đáp số

σ = −5
6.Đáp ứng bước

Đáp số: %OS = 71,24 %

Đáp ứng bước: một hệ bậc nhất G(s) = K/(τs + 1) đáp ứng với một bước nhảy có biên độ U là y(t) = K·U·(1 − e−t/τ): nó đạt 63,2 % giá trị cuối K·U sau một hằng số thời gian τ và 98 % sau 4τ. Một hệ bậc hai G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) có các cực s = −ζωₙ ± jωd với ωd = ωₙ√(1 − ζ²). Với 0 < ζ < 1, đáp ứng bước có độ vọt lố: độ vọt lố phần trăm %OS = 100·e−ζπ/√(1−ζ²), thời gian đỉnh tp = π/ωd, và thời gian xác lập 2 % ts ≈ 4/(ζωₙ). So sánh mẫu số với s² + 2ζωₙs + ωₙ² để suy ra ωₙ và ζ. Với hệ khối lượng–lò xo–giảm chấn ωₙ = √(k/m) và ζ = c/(2√(km)); với mạch RLC nối tiếp ωₙ = 1/√(LC) và ζ = (R/2)·√(C/L).

Tần số riêng và hệ số tắt dần

ωₙ = 1/√(LC) = 1/√(0,05·4,7 × 10⁻⁵) = 652,3 rad/s
ζ = (R/2)·√(C/L) = (7/2)·√(4,7 × 10⁻⁵/0,05) = 0,1073

Kết quả

%OS = 100·e−ζπ/√(1−ζ²) = 71,24 %

Đáp số

%OS = 71,24 %
7.Ổn định: Routh–Hurwitz và quỹ đạo nghiệm

Đáp số: N = 2

Tiêu chuẩn Routh–Hurwitz: một hệ tuyến tính ổn định khi mọi nghiệm của đa thức đặc trưng của nó nằm ở nửa trái của mặt phẳng s. Viết các hệ số xen kẽ vào hai hàng đầu tiên của bảng Routh. Mỗi phần tử của một hàng mới đến từ hai hàng phía trên nó; với đa thức bậc ba a₃s³ + a₂s² + a₁s + a₀ phần tử s¹ là b₁ = (a₂·a₁ − a₃·a₀)/a₂. Số lần đổi dấu trong cột đầu tiên bằng số nghiệm ở nửa mặt phẳng phải. Với độ lợi K trong các hệ số, mọi phần tử của cột đầu phải dương: độ lợi làm phần tử s¹ bằng không là giới hạn ổn định, và hàng s² khi đó cho tần số của dao động duy trì.

Đa thức đặc trưng

s³ + 14s² + 48s + 1344 = 0

Bảng Routh

s³: 1, 48
s²: 14, 1344
s¹: −48
s⁰: 1344

Số lần đổi dấu trong cột đầu tiên

1 → 14 → −48 → 1344
N = 2

Đáp số

N = 2

Độ lợi vượt quá giới hạn K = 672, nên một cặp phức đã vượt sang nửa mặt phẳng phải.

8.Ổn định: Routh–Hurwitz và quỹ đạo nghiệm

Đáp số: K = 16

Quỹ đạo nghiệm số: khi độ lợi K đi từ 0 đến ∞, các cực vòng kín, tức các nghiệm của 1 + K·G(s) = 0, di chuyển từ các cực vòng hở đến các điểm không vòng hở hoặc đến vô cực. Một điểm trên trục thực thuộc quỹ đạo khi số cực và điểm không thực ở bên phải nó là số lẻ. n − m nhánh đi đến vô cực theo các đường tiệm cận ở các góc (2q + 1)·180°/(n − m), xuất phát từ trọng tâm σₐ = (Σ cực − Σ điểm không)/(n − m). Các nhánh rời hoặc nhập vào trục thực tại nơi dK/ds = 0. Điều kiện góc ∠G(s) = 180° quyết định hình dạng, và điều kiện biên độ K = 1/|G(s)| cho độ lợi tại bất kỳ điểm nào của quỹ đạo.

Điểm tách

dK/ds = 0 cho K = −(s + 1)(s + 9) cho 2s + 10 = 0.

sb = −5

Điều kiện biên độ

K = |−5 + 1|·|−5 + 9| = 4·4 = 16

Đáp số

K = 16

Với độ lợi này, vòng kín bị cản tới hạn (ζ = 1).

9.Sai số xác lập và điều khiển PID

Đáp số: ess = 0,2

Sai số xác lập: trong một vòng phản hồi đơn vị với hàm truyền vòng hở G(s), sai số là E(s) = R(s)/(1 + G(s)), và theo định lý giá trị cuối ess = lims→0 s·E(s). Loại của hệ là số cực của G(s) tại s = 0. Các hằng số sai số là Kp = lim G(s), Kv = lim s·G(s) và Ka = lim s²·G(s), đều với s → 0. Một bước nhảy có độ lớn A cho ess = A/(1 + Kp), một tín hiệu dốc A·t cho A/Kv và một parabol A·t²/2 cho A/Ka. Một vòng loại 0 bám theo bước nhảy với sai số hữu hạn, một vòng loại 1 bám theo nó chính xác và bám theo tín hiệu dốc với sai số hữu hạn.

Hằng số gia tốc

Ka = 5·4/1 = 20

Sai số parabol

ess = 4/20 = 0,2

Đáp số

ess = 0,2
10.Sai số xác lập và điều khiển PID

Đáp số: tp = 0,3428 s

Điều khiển P, PI và PID: bộ điều khiển C(s) = Kp + Ki/s + Kd·s tác động lên sai số. Một độ lợi tỉ lệ làm vòng nhanh hơn và giảm sai số của đối tượng loại 0, nhưng không bao giờ loại bỏ nó. Thành phần tích phân thêm một cực tại s = 0: loại của hệ tăng thêm một, và một bước nhảy được bám theo không có sai số. Thành phần vi phân thêm sự cản. Để chọn các độ lợi, viết đa thức đặc trưng vòng kín và so sánh nó với s² + 2ζωns + ωn². Khi đó %OS = 100·e−ζπ/√(1−ζ²), tp = π/ωd với ωd = ωn√(1 − ζ²), và ts ≈ 4/(ζωn).

Đa thức đặc trưng

s(s + 8) + 50·2 = s² + 8s + 100

ωₙ và ζ

ωₙ = √100 = 10 rad/s
ζ = 8/(2·10) = 0,4

Thời gian đỉnh

ωd = 10·√(1 − 0,4²) = 9,165 rad/s
tp = π/9,165 = 0,3428 s

Đáp số

tp = 0,3428 s
11.Biểu đồ Bode và độ dự trữ

Đáp số: |H| = −7,023 dB

Đáp ứng tần số: thay s = jω vào hàm truyền. Độ lớn theo decibel là 20·log₁₀|H(jω)| và pha là tổng các góc của các nhân tử ở tử trừ đi các góc của các nhân tử ở mẫu. Một bộ lọc RC hoặc RL bậc nhất có tần số góc (−3 dB) tại ωc = 1/(RC) hoặc R/L, và fc = ωc/(2π); bộ lọc thông thấp cho |H| = 1/√(1 + (ω/ωc)²) và φ = −atan(ω/ωc). Trong biểu đồ Bode tiệm cận (đường thẳng), mỗi cực bẻ độ lớn xuống 20 dB/decade tại tần số góc của nó, mỗi điểm không bẻ nó lên 20 dB/decade, và một bộ tích phân 1/s cho −20 dB/decade ở mọi tần số. Tần số cắt biên ωc: |L(jω)| = 1 (0 dB); độ dự trữ pha PM = 180° + ∠L(jωc). Tần số cắt pha ω180: ∠L(jω) = −180°; độ dự trữ biên GM = −20·log₁₀|L(jω180)| dB. Vòng kín ổn định khi cả hai độ dự trữ đều dương.

Tần số góc

ωc = 1/(RC) = 1/(1000·2,2 × 10⁻⁸) = 45455 rad/s
fc = 7234 Hz → f/fc = 3600/7234 = 0,4976

Độ khuếch đại

|H| = 0,4976/√(1 + 0,4976²) = 0,4455
20·log₁₀(0,4455) = −7,023 dB

Đáp số

|H| = −7,023 dB
12.Biểu đồ Bode và độ dự trữ

Đáp số: ω180 = 4,583 rad/s

Đáp ứng tần số: thay s = jω vào hàm truyền. Độ lớn theo decibel là 20·log₁₀|H(jω)| và pha là tổng các góc của các nhân tử ở tử trừ đi các góc của các nhân tử ở mẫu. Một bộ lọc RC hoặc RL bậc nhất có tần số góc (−3 dB) tại ωc = 1/(RC) hoặc R/L, và fc = ωc/(2π); bộ lọc thông thấp cho |H| = 1/√(1 + (ω/ωc)²) và φ = −atan(ω/ωc). Trong biểu đồ Bode tiệm cận (đường thẳng), mỗi cực bẻ độ lớn xuống 20 dB/decade tại tần số góc của nó, mỗi điểm không bẻ nó lên 20 dB/decade, và một bộ tích phân 1/s cho −20 dB/decade ở mọi tần số. Tần số cắt biên ωc: |L(jω)| = 1 (0 dB); độ dự trữ pha PM = 180° + ∠L(jωc). Tần số cắt pha ω180: ∠L(jω) = −180°; độ dự trữ biên GM = −20·log₁₀|L(jω180)| dB. Vòng kín ổn định khi cả hai độ dự trữ đều dương.

Tần số cắt pha

(ω/3)·(ω/7) = 1 → ω180 = √(3·7) = 4,583 rad/s

Đáp số

ω180 = 4,583 rad/s

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