Lý thuyết môn kỹ thuật điều khiển tự động - Chương số 2

Chia sẻ: Nguyen Nhi | Ngày: | Loại File: PDF | Số trang:37

Thêm vào BST

Báo xấu

81
lượt xem 5
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

HÀM TRUYỀN ĐẠT Trước tiên ôn tập lại kiến thức về số phức và hàm phức. *Biến phức: s = σ + jw σ: Phần thực (Real part) ω: Phần ảo (Imaginary part) Nếu σ, ω là các số thực thì ta gọi là số phức, còn thay đổi s là biến phức. Biểu diễn biến phức s trên đồ thị như sau:

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: Lý thuyết môn kỹ thuật điều khiển tự động - Chương số 2

CHƯƠNG II HÀM TRUYỀN ĐẠT Trước tiên ôn tập lại kiến thức về số phức và hàm phức. *Biến phức: s = σ + jw σ: Phần thực (Real part) ω: Phần ảo (Imaginary part) Nếu σ, ω là các số thực thì ta gọi là số phức, còn thay đổi s là biến phức. Biểu diễn biến phức s trên đồ thị như sau: jω jω S σ σ O Hình 2.1 * Hàm phức: Là hàm của biến phức S G(s) = Gx + j Gy Cũng bao gồm phần thực và phần ảo. Độ lớn của G ( s ) = G x2 + G y 2 Góc q = tan-1(Gx/Gy), Chiều dương theo chiều kim đồng hồ tính từ trục thực jω -Biểu diễn trên đồ thị: G(s) Gy σ θ O Gx Hình 2.2 Hàm liên hợp của hàm G(s) là: G ( s) = Gx - j Gy http://www.ebook.edu.vn
Một hàm phức, có biến là s = σ + jω . Biến phức S phụ thuộc vào 2 đại lượng độc lập: là phần thực và phần ảo của s. Để biểu diễn hàm G(s) cần có 2 đồ thị, mỗi đồ thị có 2 chiều: Đồ thị của jω ứng với s gọi là phẳng S - - Đồ thị của phần ảo G(S) (ImG) ứng với phần thực của G(S) (ReG) gọi là phẳng G(S). Sự tương ứng giữa các điểm trong hai phẳng đó gọi là một ánh xạ hay biến đổi . Các điểm trong phẳng S được ánh xạ vào các điểm trong phẳng G(S) bằng hàm G. Ví dụ: Hàm phức G(S) = S2 + 1. Điểm S0 = 2 +j 4 được ánh xạ vào điểm G(S0) như sau (S0) = G(2 + j 4) = -11 + j 16 ImG jω ¹G ¸nh x G(S ) 0 16 S j4 0 σ ReG O O 2 -11 Ph¼ng G(S) Ph¼ng S Hình 2.3 * Phẳng S (mặt phẳng phức) Nếu G(S) là hàm hữu tỉ như sau: m bm ∑ ( S + z i ) i =1 G(S) = n ∑ (S + p ) i i =1 - Các giá trị của biến phức S = -zi làm cho G(s) = 0 được gọi là các không của G(s) (Zeros) Các giá trị s = - pi làm cho G(s) → ∞ được gọi là các cực của G(s) ( Poles) - Các cực và các không được xác định bởi: một đại diện phần thực và một đại diện phần ảo của số phức. Biểu diễn các điểm đó trên mặt phẳng phức ( phẳng S) gọi là ánh xạ cực – không của G(s) Ví dụ: 2S 2 − 2S − 4 2( S + 1)( S − 2) = G(s) = S + 5S + 8S + 6 ( S + 3)( S + 1 − j )( S + 1 − j ) 3 2 http://www.ebook.edu.vn
G(s) có các không: s = -1 ; s = 2 và các cực: s = -3; s = -1 – j ; s = -1 +j jω j σ -3 -1 2 Pole -j Zero Ph¼ng S Hình 2.4 *Phẳng G(s): Được biểu diễn trong mặt phẳng với 2 thành phần. Một là phần thực của G(s) – ReG, và một là phần ảo của G(s)- ImG. ánh xạ từ các điểm s0 sang phẳng G(s) là các điểm G(s0). jω ImG ¸nh x¹ G S2 S1 S3 G(S1) G(S2) σ ReG G(S4) S4 Ph¼ng S G(S3) Ph¼ng G(S) Hình 2.5 * Nhận xét: Mối quan hệ giữa phẳng S ( ánh xạ cực – không) *Phép biến đổi Laplace Biến đổi Laplace là cơ sở của một phương pháp giải tích để tìm cả đáp ứng ổn định và đáp ứng quá độ mà các phương trình vi phân tuyến tính hệ số không đổi. Nên phép biến đổi Laplace chỉ dùng biến đổi cho phương trình vi phân tuyến tính. Biến đổi Laplace chuyển phương trình vi phân thành các phương trình đại số nên tìm nghiệm của phương trình đại số đơn giản hơn và từ nghiệm của phương trình đại số tìm được nghiệm của phương trình vi phân. Một ưu điểm là phương pháp này có thể xử lý trực tiếp các điều kiện đầu của hệ thống như một phần của đáp ứng. - Bản chất của phép biến đổi Laplace: Là các phép tính đạo hàm và tích phân gốc được chuyển thành các phép toán đại số thông thường đối với các ảnh, miền xác định rộng. - Hàm gốc: http://www.ebook.edu.vn
Gọi hàm f(t) của biến thực t là hàm gốc nếu nó thoả mãn các điều kiện sau: 1. Hàm f(t) liên tục trên từng đoạn thuộc miền xác định mà t ≥ 0. Giải thích: Lấy [a; b] trên t ≥ 0, luôn chi được trong [a; b] một số hữu hạn khoảng nhỏ [e; x] sao cho trong mỗi khoảng đó f(t) liên tục và tại các mút của mỗi khoảng nhỏ thì f(t) có giới hạn một phía: lim f (t ) < ∞ t →ξ 2. Khi t → +∞ hàm f(t) không tăng nhanh hơn một hàm mũ. Tồn tại M > 0; f (t ) ≤ e α .t ; mọi t >0 a >0 sao cho: a gọi là chỉ số tăng của f(t). 3.f(t) = 0 khi t < 0. Điều kiện này được đưa ra vì trong ứng dụng biến số t thường là thời gian, hàm f(t) biểu diễn một quá trình nào đó mà ta chỉ khảo sát lúc t > 0. Một số ví dụ: a) Hàm h(t) = 0 khi t < 0 1 khi t > 0 Là một hàm gốc : η (t ) ≤ 1 thoả mãn điều kiện hàm f(t) không tăng nhanh hơn một hàm mũ. t ≥ 0 ta lấy t thuộc trong [-1; 1] thì lim η (t ) = 1 ( thoả mãn điều kiện 1) t → +1 η(t) h(t) = 0 khi t < 0 (thoả mãn điều kiện 3) O t Hình 2.6 b) Hàm f(t) = h(t). sint = 0 khi t < 0 sint khi t > 0 η(t).sint η (t ). sin t ≤ 1 = M .e αt ( M = 1; a = 0) η (t ). sin t liên tục trên t ≥ 0 O t η (t ). sin t = 0 khi t 0 http://www.ebook.edu.vn
η (t ).t 2 ≤ 2.e t ( M = 2; a = 1) η(t).t2 O t Hình 2.8 - Toán tử Laplace: Nếu f(t) là một hàm gốc có chỉ số tăng là a thì yêu cầu của f(t) để chuyển đổi được là: ∞ − σ .t dt < ∞ ( a < s < ∞ ) tích phân hội tụ tuyệt đối. ∫ f (t ) .e 0 Biến đổi Laplace là kết quả của một thuật toán chuyển đổi với một hàm thời gian f(t) để cho ta hàm G(s) của biến phức s. ∞ T lim ∫ f (t ).e − st dt = ∫ f (t ).e − st dt ( 0 < e < T ) F(s) = L {f(t)} = T → ∞ε 0+ Biến đổi ngược để tìm gốc f(t): 1 σ + j∞ − st .ds f(t) = ∫ F ( s ).e 2πj σ − j∞ Một số hàm biến đổi Laplace sử dụng phổ biến: Important Laplace Transform Pairs f(t) F(s) Hàm bậc thang h(t) 1 S Hàm xung đơn vị d(t) = 0 nếu t < 0 1 1 nếu 0 ≤ t ≤ t1 0 nếu t > t1 t 1 S2 tn n! S n +1 e-at 1 S+a http://www.ebook.edu.vn
t n −1 .e − at 1 (n − 1)! (S + a) n 1 1 .(1 − e − at ) S(S + a) a sinwt ω S +ω 2 2 coswt S S +ω 2 2 e-at.f(t) F(s + a) skF(s) – sk-1f(0-) – sk-2f’(0-) - ... – f(k-1)(0-) d k f(t) f(k)(t) = dt k t 0 F(s) 1 ∫ ∫ + f(t)dt f(t)dt s s −∞ −∞ * Lưu ý: d Biến s được coi như phép vi phân: s ≡ dt t Và trong tích phân: 1≡ ∫ dt s0 * ứng dụng của toán tử Laplace: - Giải các phương trình vi phân tuyến tính hệ số không đổi - Tìm hàm truyền đạt của hệ thống điều khiển tuyến tính. 2.1. Hàm truyền đạt * Định nghĩa: Hàm truyền đạt (The Transfer function) của một hệ thống tuyến tính được định nghĩa là tỷ số giữa biến đổi Laplace của biếu ra ( đại lượng đáp ứng ra của hệ thống) so với biến đổi Laplace của biến vào ( đại lượng tác động vào hệ thống), Với điều kiện đầu đồng nhất bằng không. Hàm truyền đạt của hệ thống ( phần tử) đặc trưng cho mô tả động lực học của hệ thống. - Một hàm truyền đạt chỉ có thể xác định cho hệ thống tuyến tính, hệ thống bền vững ( tham số không đổi). Một hệ thống không bền vững thường gọi là hệ thống biến thời gian thay đổi, có một hay nhiều tham số thay đổi, và phép biến đổi Laplace không được áp dụng đối với hệ thống này. - Hàm truyền đạt thể hiện tác động vào và đáp ứng ra của trạng thái hệ thống. - Tuy nhiên, hàm truyền đạt không diễn tả thông tin về cấu trúc bên trong của hệ thống và trạng thái hoạt động của hệ thống. http://www.ebook.edu.vn
Output Y(s) l {y(t)} G(s) = = = l {u(t)} Input U(s) Để hiểu về cách xây dựng hàm truyền đạt ta có các ví dụ sau: Ví dụ 1: Cho một hệ thống được mô tả bởi phương trình vi phân sau: d2y dy +4 + 3y = 2r(t) 2 dt dt dy (0) = 0 , và r(t) = 1, t ≥ 0 Điều kiện đầu là: y(0) = 1, dt Biến đổi Laplace: [ s2Y(s) – s y(0) ] + 4[ s Y(s) – y(0) ] + 3 Y(s) = 2 R(s) 1 Thay R(s) = và y(0) = 1 ta được: s 2 s2Y(s) – s + 4 s Y(s) – 4 + 3 Y(s) = s s+4 2 +2 Y(s) = s(s + 4s + 3) (s + 4s + 3) 2 Trong đó: q(s) = s2 +4s + 3 = ( s + 1)(s +3) = 0 là phương trình đặc trưng và d(s) = s − 1/2 −1 3/2 1/3 2/3 + ]+[ + ]+ Y(s) = [ = Y1(s) + Y2(s) + Y3(s) (s + 1) (s + 3) (s + 1) (s + 3) s Biến đổi Laplace ngược: 3 − t 1 − 3t 1 2 .e − .e ] + [ − 1.e − t + .e − 3t ] + y(t) = [ 2 2 3 3 2 Trạng thái ổn định là: lim y(t) = t→∞ 3 Ví dụ 2: Hệ thống cơ khí như hình vẽ ( được mô hình hoá) friction f 2 K M2 V2(t) friction f 1 M1 V1(t) Force r(t) Hình 2.9 http://www.ebook.edu.vn
Trong hình vẽ: K: Độ cứng lò xo ( hằng số lò xo) f1, f2: là các hệ số ma sát V1(t), V2(t): Vận tốc di chuyển của các trọng khối M1 và M2. M1sV1(s) + (f1 + f2)V1(s) – f2V2(s) = R(s) V2 (s) M2sV2(s) + f1(V2(s) – V1(s)) + K =0 s Tương đương với: (M1(s) + (f1 + f2)) V1(s) + (- f1)V2(s) = R(s) K (-f1)V1(s) + (M2(s) + f1 + ) V2(s) = 0 s Hoặc dưới dạng ma trận sau: ⎡(M 1 (s) + f1 + f2 ........(− f1 ) ⎤ ⎡V1 (s) ⎤ ⎡R(s)⎤ ⎢ K ⎥.⎢ = ⎢.........(−f1 )......(M 2 (s) + f1 + )⎥ ⎣V2 (s)⎥ ⎢0 ⎥ ⎦⎣ ⎦ ⎣ s⎦ Vận tốc di chuyển của M1 chính là đại lượng ra, việc tìm V1(s) bởi ma trận nghịch đảo hoặc nguyên tắc Cramer là: (M 2 s + f1 + (K/s)).R(s) V1(s) = (M 1 s + f1 + f2 ).(M 2 s + f1 + (K/s)) − f12 Hàm truyền đạt của hệ thống: (M 2 s + f1 + (K/s)) V1 (s) = G(s) = = (M 1 s + f1 + f2 ).(M 2 s + f1 + (K/s)) − f12 R(s) (M 2 s 2 + f1 s + K).R(s) (M 1 s + f1 + f2 ).(M 2 s 2 + f1 s + K) − f12 s Tại một thời điểm nào đó mà xác định x1(t), thì hàm truyền đạt là: X (s) V1 (s) G(s) = = R(s) sR(s) s Ví dụ 3: Hàm truyền đạt của động cơ dc Động cơ dc là thiết bị phát động mà chuyển từ dạng năng lượng điện sang chuyển động quay. Armature Ra La Rf + Lf Inertia = J Friction = f - http://www.ebook.edu.vn if(t) Load Field
Hình 2.10 Ví dụ 4: Cho hệ cơ học gồm một lò xo có hệ số c, một vật với khối lượng m và bộ giảm chấn có hệ số d được nối với nhau như hình vẽ. Xác định hàm truyền đạt cho hệ cơ đó nếu tín hiệu đầu vào u(t) được định nghĩa là lực bên ngoài tác động lên vật và tín hiệu ra y(t) là quãng đường mà vật đi được. Gọi Fc, Fm, Fd là những lực của lò xo, vật và bộ giảm chấn sinh ra khi vật di chuyển nhằm cản sự dịch chuyển đó thì: Fc = c. y(t) d 2 y(t) Fm = m. c dt 2 Fc Fm u(t) dy(t) Fd = d . m dt Theo tiên đề về cân bằng lực ta được: Fd y(t) 2 d y(t) dy(t) d u(t) = Fc + Fm + Fd = c . y(t) + m. +d. 2 dt dt Biến đổi Laplace: U(s) = ( c + ds + ms2 ). Y(s) Hình 2.11 Hàm truyền đạt của hệ thống là: Y(s) 1 G(s) = = ms + ds + c 2 U(s) Gọi g(t) là hàm gốc của hàm truyền đạt G(s), tức là: g(t) = L-1{G(s)} Theo tính chất của toán tử Laplace ta có: Y(s) = G(s). U(s) +∞ +∞ ∫ g(τ ).u(t − τ )dτ ∫ g(t - τ ).u(τ )dτ ⇔ y(t) = g(t). u(t) = = −∞ −∞ Hàm g(t) được gọi là hàm trọng lượng của hệ thống. Với u(t) = δ (t ) Do U(s) = 1 nên ta có y(t) = g(t) * Hàm truyền đạt trong lĩnh vực Laplace Trên đây mới chỉ giới thiệu hàm truyền đạt giới hạn trong quan hệ tỷ lệ vào – ra đơn giản, đó là một hình thức để mô tả đặc trưng của phần tử hoặc hệ thống. Tuy nhiên có nhiều phần tử có đáp ứng thay đổi theo thời gian. Trong lĩnh vực thời gian đặc tính đó được mô tả bằng phương trình vi phân, phương trình này không trực tiếp dùng làm hàm truyền đạt được. Nếu dùng một hàm truyền đạt với biến số Laplace S, diễn tả được đặc tính động lực của phần tử hoặc hệ thống và phương pháp phân tích trong lĩnh vực thời gian ( tức http://www.ebook.edu.vn
là quá trình quá độ)sẽ tương đối đơn giản giúp ta xác định đáp ứng của phần tử hoặc hệ thống đối với một tín hiệu vào xác định. Đặc trưng của một hệ thống điều khiển, ta có phương trình vi phân tổng quát sau đây: (pn + bn-1pn-1 + ... + b1p + b0). y(t) = ( ampm + am-1pm-1 + ... + a1p + a0 ). x(t) (2.11) a m p m + a m −1 p m −1 + ... + a 1 p + a 0 L (p) .x(t) = m .x(t) y(t) = p + b n −1 p + ... + b 1 p + b 0 n −1 n L n (p) Trong đó: a0, ..., am và b0, ..., bn là những hằng số x(t) hàm kích thích, nó là tín hiệu tác động vào làm kích thích hệ thống y(t) hàm phản ứng. Nó là hàm chuyển tiếp (tín hiệu ra) dưới tác động của tín hiệu vào x(t). Ln(p) = pn + bn-1pn-1 + ... + b1p + b0 Lm(p) = ampm + am-1pm-1 + ... + a1p + a0 Biến đổi Laplace từng số hạng của phương trình (2.11) ta có L[pn y(t)] = sn Y(s) – I(s)n bn-1 L[pn-1 y(t)] =bn-1. sn-1 Y(s) – I(s)n-1 ... am L[pm x(t)] = am sm X(s) – I(s)m am-1 L[pm-1 x(t)] =am-1 sm-1 Y(s) – I(s)m-1 Với I(s)n,...là những điều kiện ban đầu tương ứng với các biến đổi. Thay vào phương trình: (a m s m + a m −1 s m −1 + ... + a 1 s + a 0 ).X(s) + I(s) L m (s).X(s) + I(s) = Y(s) = s n + b n −1 s n −1 + ... + b 1 s + b 0 L n (s) I(s) = I(s)n + I(s)n-1 + ... - I(s)m - I(s)m-1 là tổng những điều kiện đầu. Từ phương trình trên thấy rằng: - các đa thức Lm (s); Ln(s) ở trong miền biến đổi s vẫn giữ nguyên như trong miền toán tử p. - Tử số của chúng cũng có dạng giống nhau, chỉ khác là ở miền s có các điều kiện đầu I(s). - Nếu các điều kiện đầu bằng 0 thì ta có thể biến đổi Laplace của phương trình vi phân bằng cách thay s vào vị trí p, thay Y(s) vào vị trí y(t) và X(s) vào vị trí x(t). Tức là: L m (s) Y(s) = . X ( s) L n (s) L m (s) Và hàm truyền đạt là G(s) = L n (s) http://www.ebook.edu.vn
Vậy có mối quan hệ trong hệ thống điều khiển: “Hàm phản ứng = Hàm truyền đạt x Hàm kích thích” Nếu cho mẫu số của hàm truyền đạt bằng 0 ta sẽ có phương trình đặc trưng: sn + bn-1sn-1 + ... + b1s + b0 = 0 trên cơ sở phương trình đặc trưng ta suy ra các đặc tính chuyển tiếp của hệ thống. - Hàm phản ứng (hàm chuyển tiếp) y(t) có thể xác định với việc biến đổi ngược hàm Y(s) L m (s).X(s) + I(s) y(t) = L-1[ Y(s)] = L-1 [ ] L n (s) Tìm y(t) theo 2 cách: 1) Dùng bảng để xác định các hàm thời gian tương ứng 2) Phân tích hàm đã biến đổi thành tổng những hàm đơn giản hơn và sau đó dùng bảng để biến đổi ngược từng số hạng. Thường dùng phương pháp 2 vì ít khi gặp các hàm đơn giản. Vậy ta tìm hiểu phương pháp 2 như sau: L m (s).X(s) + I(s) Y(s) = L n (s) Hàm kích thích X(s) hay tín hiệu vào có thể viết dưới dạng sau đây: Nx X(s) = Dx L m (s).N x + I(s).D x A( s ) = Y(s) = B( s) L n (s).D x ở đây A(s) và B(s) là những đa thức của s. Để có thể chia Y(s) thành các phân thức, ta phân tích mẫu số B(s). Giả sử các nghiệm của B(s) là r1, r2, ..., rn. Các nghiệm này có thể là nghiệm đơn, nghiệm bội hay là số phức. - Nghiệm đơn: C C C A(s) = 1 + 2 + ... + n Y(s) = B(s) s − r1 s − r2 s − rn Xác định C1, C2, ..., Cn ta dùng phương pháp sau C1 = lim[(s − r1 ).Y(s)] s → r1 C2 = lim[(s − r2 ).Y(s)] s → r2 ... Cn = lim[(s − rn ).Y(s)] s → rn Biết được C1, C2, ..., Cn tìm được biến đổi Laplace ngược trong bảng: http://www.ebook.edu.vn
Cn L-1 [ ] = Cn. e rn t ( t ≥ 0 ) s - rn Vậy, y(t) = L-1[ Y(s)] = C1. e r t + C2. e r t +... + Cn. e r t n 1 2 Hàm chuyển tiếp y(t) mong muốn là hàm tắt dần nên từng phần Cn. e r t là hàm tắt n dần, tức là tất cả các nghiệm r1, r2, ..., rn cần phải là số âm. - Nghiệm bội: B(s) = (s-r)q.(s-r1).(s-r2)...(s-rn) Kq K q −1 K C C A(s) = + + ... + 1 + 1 + ... + n Y(s) = q −1 s − r s − r1 s − rn B(s) (s − r) (s − r) q Xác định Kq: Kq = lim[(s − r) q .Y(s)] s→ r Còn các hệ số còn lại xác định xác định bằng cách: (s − r) q −1 (s − r) q d [(s − r) q Y(s)] = K q −1 + 2K q −2 (s − r) + ... + qC 1 + ... + qC n (s − ri ) (s − rn ) ds d Kq-1 = lim{ [(s − r) q Y(s)]} s→ r ds 1 d2 [(s − r) q Y(s)]} Kq-2 = lim{ 2 s→ r 2 ds .... 1 d (k) [(s − r) q Y(s)]} Kq-k = lim{ (k) s→ r k! ds Vậy: K q .t q −1 .e rt K q −1 .t q −2 .e rt K 2 .t.e rt + + ... + + K 1 .e rt + C1. e r1t + C2. e r2t +... + Cn. y(t) = (q − 1)! (q − 2)! 1! ert n -Nghiệm phức liên hợp: C0 C C A(s) C = + + 1 + ... + n Y(s) = B(s) s − a - jb s − a + jb s − r1 s − rn Xác định các hằng số C, C0 tương ứng với các nghiệm phức liên hợp: C = A(s) A(s) lim [(s − a − jb). ] = lim [ ] (s − a − jb).(s − a + jb).(s − r1 )...(s − rn ) s→a + jb (2jb).(s − r1 )...(s − rn ) s →a + jb 1 .K(a + jb) = 2jb http://www.ebook.edu.vn
4( s ) A/ s = [(s 2 -2as+a 2 +b 2 . K(a+jb) = lim )] S = a + jb ( s − r1)...( S − rn ) B/s s → a + jb A( s ) Co = lim [(s-a+jb) ] ( s − a − jb).( s − a + jb).( s − r1 )....( s − rn ) s → a − jb A( s ) 1 = lim [ ] = - .K.(a-jb) − 2a.( s − r1 ).( s − rn ) 2a s → a − jb A( s ) 1 K(a+jb) = lim [ ] = - .K.(a+jb) ( s − r1 ).( s − rn ) 2a s → a − jb A( s ) = [( 2 - 2as + a 2 +b 2 ) ] s = a − jb B( s) Các trị số k(a+jb) và k(a-jb) là các số phức liên hợp Ta cần thể hiện các số này trên hình vẽ: K(a+jb) = [k(a+jb)]e jα K(a+jb) = [k(a+jb)]e jα [k(a+jb)] = [k(a-jb)] ( Độ dài của véc tơ ) ⇒ C và Co cũng là các số phức liên hợp . 1 .[k(a+jb)].e − jα C= 2 jb Co = -[k(a+jb)]e − jα Từ bảng laplace ta xác định hàm chuyển tiếp y (t ) = c.e ( a + jb ).t +Co.e ( a − jb ).t +C1.e r1t +.....+Cn.e rnt 1 1 [k(a+jb)].e ( a + jb ).t .e − jα + - [k(a-jb).e ( a − jb ).t .e − jα +... ⇒ y (t ) = 2 jb 2 jb 1 [k(a+jb)].e a.t .e j (bt +α ) - e − j (bt +α ) = 2 jb e j (bt +α ) − e − j ( bt +α ) 1 [k(a+jb)].e at . = 2j b 1 [k(a+jb)].e at .sin( (α + bt ) +C1.e r1t +...+Cn.e rnt = b Phương trình trên thể hiện hàm điều hoà sin tắt dần theo hàm mũ, xuất phát từnghiệm phức liên hợp Phần ảo b là tàn số dao động tắt dần . Thời gian của mỗi 2π 1 . Đường bao hình sin là [k(a+jb)].e at . Để hàm mũ giảm dần thì a dao động là b b http://www.ebook.edu.vn
1 [k(a+jb)].e at phải là trị số âm . Trường hợp a = 0, ta sẽ có hàm sin có biên độ b không đổi. j j at (1/b).[K(a+ jb)].e (1/b).[K(a+ jb)] t t O O a0 a0 Hình 2.12 - Nếu các nghiệm nằm ở bên trái trục ảo ( a0) thì dao động sẽ tăng dần. Cách khác xác định đáp ứng thời gian: Đáp ứng thời gian có thể xác định bằng cách tìm các cực của G(s). X(s) vì Y(s) = G(s). X(s) và ước lượng tìm các hệ số của các phân thức của biểu thức Y(s) tại các cực đó. Các hệ số có thể xác định bằng đồ thị nhờ một ánh xạ cực – không của Y(s). ánh xạ này được dựng từ ánh xạ cực – không của G(s) và cộng thêm các cực- không của X(s). Các bước : m b m . ∏(s + z i ) i =1 G(s) = n ∏(s + p i ) i =1 Vì G(s) là một hàm phức nên có thể viết dưới dạng cực như sau: G(s) = P(s) .e jφ = P(s) ∠φ ⎡ Im G ( s ) ⎤ φ = tan −1 ⎢ ⎥ ⎣ Re G ( s ) ⎦ http://www.ebook.edu.vn
Mỗi số phức s, zi, pi, ( s + zi) và ( s + pi) có thể diễn tả bằng một vectơ trong mặt phẳng S. Biểu diễn trên đồ thị: jω jω (s) -zi -p1 -p2 σ σ s+zi -z1 -z2 -p3 s -pi s+pi jω a) b) (s) -p1 -p2 σ -z1 -z2 -p3 c) Hình 2.13 Trong hình a) có một cực –pi và một không – zi và một biến phức S. Vectơ tổng s + zi là vectơ bắt đầu từ không – zi và kết thúc tại s, vectơ s + pi bắt đầu từ cực – pi và kết thúc tại s. m b m . ∏ (s + z i ) § é lín cña vecto (s + zi) i =1 Độ lớn của C = bm. = § é lín cña vecto ( s + pi) n ∏ (s + p i ) i =1 b m . (s + z 1 ).(s + z 2 ) Trường hợp b): C1 = (s + p 1 ).(s + p 2 ) b m . (s + z 1 ).(s + z 2 ) Trường hợp c): C 2 = (s + p 1 ).(s + p 2 ) Diễn tả theo dạng cực thì: Ci = C i .e jφ = C i ∠φi i Hoặc theo toạ độ vuông góc: Ci = Ci . cos φi + j. Ci . sin φi φi = Tổng các góc của các vectơ từ các không đến – pi trừ đi tổng các góc của các vectơ từ các cực tới –pi ( nếu bm > 0) φi = Tổng các góc của các vectơ từ các không đến – pi trừ đi tổng các góc của các vectơ từ các cực tới –pi + 1800 ( nếu bm < 0) -Phương pháp đồ thị này không áp dụng cho trường hợp có các cực trùng nhau ( nghiệm lặp). http://www.ebook.edu.vn
* Hàm truyền đạt trong lĩnh vực tần số Việc phân tích hệ thống nằm trong hai lĩnh vực: Lĩnh vực thời gian và lĩnh vực tần số. -Trong lĩnh vực thời gian: nội dung chủ yếu là các đặc tính động lực của hệ thể hiện trạng thái quá độ (đáp ứng quá độ). Ta đã dùng phương trình vi phân và biến đổi Laplace để nghiên cứu các nghiệm của phương trình ( tức là các đáp ứng của hệ). Áp dụng biến đổi Laplace để giải các phương trình vi phân tuyến tính là phần quan trọng nhất trong nghiên cứu trạng thái quá độ của các hệ tuyến tính thuộc lĩnh vực thời gian. Tuy vậy giải phương trình vi phân để phân tích trạng thái động lực của hệ thống (tức là trong lĩnh vực thời gian) khá phức tạp đối với các hệ không đơn giản. Nhưng phương pháp phân tích đáp ứng tần số ( thuộc lĩnh vực tần số) có thể đánh giá được tính năng của hệ mà không cần giải phương trình vi phân. Phương pháp đáp ứng tần số phân tích các tính năng của hệ xem như một hàm của tần số của tín hiệu vào dạng sin mà không phải là khảo sát đáp ứng thời gian thực tế. Cũng có thể nói phương pháp đáp ứng tần số phân tích đáp ứng dạng sin ổn định của hàm truyền của hệ. Phương pháp này có nhiều ưu điểm: - Cho phép ta ước lượng được dãy tần số ảnh hưởng đến tính năng của hệ - Dễ chỉ cho ta biện pháp thay đổi hệ để đạt các tính năng yêu cầu trong việc thiết kế các hệ thống điều khiển. Bằng đồ thị có thể chỉ cho ta biện pháp phán đoán vấn đề bằng các phương trình vi phân. Nếu các phương trình đã được giải nhưng đáp ứng không đạt yêu cầu thì không dễ quyết định được biện pháp thay đổi hệ thống để đạt chất lượng mong muốn. Phương pháp tần số đã vượt qua được hạn chế đó. - Đáp ứng có thể xác định bằng thực nghiệm cũng tốt không thua kém tính toán giải tích. Ưu điểm này rất quan trọng khi mô tả các phần tử của hệ bằng các phương trình vi phân. 2.2. Đại số sơ đồ khối Sơ đồ khối là một trong các dạng mô hình toán của hệ thống điều khiển, trên sơ đồ thể hiện đại lượng vào – ra của hệ thống và các tính chất của hệ thống. Một số chuyển đổi cơ bản để rút gọn các sơ đồ khối phức tạp. 1. Tổ hợp các khối nối tiếp R C R C G1 G2 G1xG2 Hình 2.14 http://www.ebook.edu.vn
Chứng minh : C= R.G1.G2 = G1.G2.R 2. Tổ hợp các khối song song R C R C + G1 G1+G2 + R G2 Hình 2.15 Tại điểm tụ C = R.G2 + R.G1 = ( G1+G2).R 3. Di chuyển điểm tụ về bên phải một khối : C A + C + A G G + + G B B Hình 2.16 Tại điểm tụ R = A + B Nên C = G. ( A + B) Sơ đồ tương đương là: C = A. G + B. G = G. ( A + B) 4. Di chuyển điểm tụ về bên trái một khối A C C + A G G + + 1/G B B Hình 2.17 5. Di chuyển điểm tán về bên phải một khối C B C B G G C 1/G B Hình 2.18 6. Di chuyển điểm tán về bên trái một khối C B B C G G C C G Hình 2.19 7. Rút gọn hệ thống http://www.ebook.edu.vn
E R+ C R C G G 1+GH B- H Hình 2.20 Chứng minh: C ⇔ C = E. G ; E = R - B ; B = C. H ; Sơ đồ ban đầu: G = E E = R – C. H = R – E. GH E. ( 1 + GH ) = R R E= 1 + G.H C R G G Hàm truyền của hệ thống là: = .= 1 + G.H R 1 + GH R Từ biểu thức ta thấy: Nếu gia lượng tuyến thuận G lớn thì tích GH ≥ 1, lúc này gia C 1 lượng mạch kín còn là = R H - Kết luận: Trạng thái của mạch kín phụ thuộc tính chất tuyến tính của phản hồi H và độc lập với tuyến thuận ( về tính chất ). Nếu tuyến thuận có một vài thay đổi do một vài lí do nào đó thì tuyến phản hồi sẽ trừ khử hiệu quả sự thay đổi của đầu ra. Vì thế không cần điều chỉnh hệ thống, nhưng phải điều chỉnh phần tử phản hồi H. - Nếu mạch kín bị cắt đứt như hình vẽ: E R+ C G G 1 2 + B- H Hình 2.21 Hàm truyền toàn mạch còn G1. G2. H được xem như hàm truyền của mạch hở. B C.H E.G 1 .G 2 .H = = G1. G2. H = E E E * Sơ đồ khối dạng chính tắc: R + E C G + B- H ư Hình 2.22 http://www.ebook.edu.vn
Các đại lượng sau cần xác định rõ: G: Hàm truyền tuyến thuận H: Hàm truyền tuyến phản hồi GH: Hàm truyền mạch hở. E C B G H Hình 2.23 C : Hàm truyền mạch kín (tỷ số điều khiển) R E : Tỷ số tín hiệu tác động ( tỷ số sai lệch ) R B : Tỷ số phản hồi cơ bản R C G Ta có liên hệ sau: = 1 ± GH R * Hệ phản hồi đơn vị: Một hệ phản hồi đơn vị là một hệ trong đó tín hiệu phản hồi cơ bản B bằng đầu ra C. Đây là một trường hợp đặc biệt hay gặp trong thực tế và là sự so sánh trực tiếp giữa đầu ra và đầu vào chuẩn. Vì lúc này khối phản hồi có giá trị đơn vị là 1 nên hàm truyền mạch kín là: C G = 1+ G R Trường hợp này xảy ra khi đầu ra mô phỏng lại đầu vào chuẩn. Bất kỳ hệ phản hồi nào nếu chỉ có các phần tử tuyến tính trong tuyến phản hồi đều có thể đặt dưới dạng một hệ phản hồi đơn vị bằng cách dùng chuyển đổi 4, ta được sơ đồ khối sau: E E R+ C R C + 1 G G.H H + + B- B- H R E=R ±B ±B E= H B = C.h B=C Hình 2.24 C G ⇒ = 1 ± GH R http://www.ebook.edu.vn
C C C R ±C ⇒ E= E= = G G.H G.H H R 1 ± 1) = C( G.H H * Hệ có nhiều tín hiệu vào ra : Nhiều hệ có nhiễu U, hoặc có nhiều tín hiệu vào ( nhiều kích thích ) đồng thời với tín hiệu vào chuẩn R, chúng áp lên hệ tại các điểm khác nhau và mang lại cho hệ những tính năng khác nhau. Khi trong một hệ tuyến tính có mặt nhiều tín hiệu vào ta phải xử lí từng tín hiệu độc lập với nhau, sau đó dựa trên nguyên lí chồng chất cộng đại số các đáp ứng cá biệt của từng tín hiệu với nhau ta sẽ được tín hiệu ra tổng cộng của hệ khi mọi tín hiệu đồng thời tác động lên hệ. Có nghĩa là ta giả thiết từng tín hiệu vào tác dộng riêng biệt đến hệ ( các tín hiệu vào còn lại giả thiết bằng không ) lần lượt làm như vậy với từng tín hiệu vào, sau đó thực hiện một phép cộng đại số các đáp ứng nói trên, để tìm đáp ứng riêng của từng tín hiệu vào, đôi khi cần đến thủ thuật rút gọn sơ đồ khối về dạng chính tắc bằng cách dùng một trong bảy chuyển đổi trên. * MỘT SỐ VÍ DỤ Ví dụ 1 : Xác định dầu ra C của hệ thống: U R+ C ++ G1 G2 - B H Hình 2.25 Cho U = 0 hệ thống đơn giản hoà thành : R+ C(R) G1 G2 - H Hình 2.26 Xác định đầu ra : G1.G2 C (R) = .R 1 + G1.G 2 + Cho R = 0 , chỉ có đầu vào U ta có sơ đồ sau http://www.ebook.edu.vn