Bài giảng môn học Kỹ thuật điều khiển tự động (Chuẩn kiến thức)

các phần tử điều khiển cơ bản 1.4 Phân loại hệ thống điều khiển 1.4.1. Hệ thống điều khiển hở và hệ thống điều khiển kín. 1.4.2. Hệ thống điều khiển liên tục và gián đoạn 1.5 Tuyến tính hóa các hệ thống phi tuyến 1.6 Ứng dụng MatLab Chương 2. HÀM TRUYỀN ĐẠT 2.1 Hàm truyền đạt 2.2 Sơ đồ khối - Đại số sơ đồ khối 2.3 Graph tín hiệu và qui tắc Mason 2.4. Các hệ thống lấy mẫu dữ liệu 2.5 Hàm truyền đạt của hệ thống rời rạc 2.6 Ứng dụng MatLab Chương 3. KHÔNG GIAN TRẠNG THÁI. 3.1 Các mô hình không gian trạng thái. 3.2 Mô hình không gian trạng thái và các phương trình vi phân 3.3 Xác định biến trạng thái từ hàm truyền 3.4 Xác định hàm đáp ứng từ phương trình trạng thái 3.5 Ứng dụng MatLab Chương 4. ỔN ĐỊNH CỦA HỆ THỐNG ĐIỀU KHIỂN TUYẾN TÍNH. 4.1 Khái niệm chung 4.2 Khái niệm ổn định và các định nghĩa chính 4.3 Trị riêng và tính ổn định của hệ thống 4.4 Các tiêu chuẩn ổn định 4.5 Ứng dụng MatLab Chương 5. TÍNH ĐIỀU KHIỂN VÀ QUAN SÁT ĐƯỢC CỦA HỆ THỐNG ĐIỀU KHIỂN. 5.1 Tính điều khiển được của các hệ thống liên tục. 5.2 Tính quan sát được của các hệ thống liên tục. 5.3 Tính điều khiển được của các hệ thống gián đoạn. 5.4 Tính quan sát được của các hệ thống gián đoạn. 5.5 Ứng dụng MATLAB. Chương 6. THIẾT KẾ HỆ THỐNG ĐIỀU KHIỂN. 6.1 Mở đầu. 6.2 Các khâu động học của hệ thống điều khiển. Chương 7. THIẾT KẾ HỆ THỐNG ĐIỀU KHIỂN BẰNG THUỶ LỰC. 7.1. Các phần tử cơ bản 7.1.1. Bơm dầu. 7.1.2. Van tràn, van an toàn. 7.1.3. Van giảm áp 7.1.4. Bộ điều chỉnh và ổn định tốc độ. 7.1.5. Van điều khiển. 7.1.6. Cơ cấu chấp hành. I. Phần 1: Phần lý thuyết I.1. Yêu cầu đối với sinh viên - Mục tiêu: Nội dung cơ bản của hệ thống điều khiển tự động, Phân tích và tổng hợp được một hệ thống điều khiển. - Nhiệm vụ của sinh viên: Dự học lý thuyết: đầy đủ Thảo luận: đầy đủ. - Đánh giá: Chấm điểm Thảo luận : 20% Kiểm tra giữa kỳ: 20% Thi kết thúc học phần : 60% I.2. Các nội dung cụ thể Chương 1 CÁC VẤN ĐỀ CƠ BẢN CỦA HỆ THỐNG ĐIỀU KHIỂN TỰ ĐỘNG 1.1- Các nội dung cơ bản của hệ thống điều khiển. * Điều khiển: Là tác động lên đối tượng để đối tượng làm việc theo một mục đích nào đó. * Hệ thống điều khiển: Là một tập hợp các thành phần vật lý có liên hệ tác động qua lại với nhau để chỉ huy hoặc hiệu chỉnh bản thân đối tượng hay một hệ thống khác. * Xung quanh ta có rất nhiều hệ thống điều khiển nhưng có thể phân chia thành 3 dạng hệ thống điều khiển cơ bản. - Hệ thống điều khiển nhân tạo. - Hệ thống điều khiển tự nhiên (bao gồm điều khiển sinh vật). - Hệ thống điều khiển tự nhiên và nhân tạo. Trong các hệ thống đó đối tượng điều khiển có thể là hệ thống vật lý, thiết bị kỹ thuật, cơ chế sinh vật, hệ thống kinh tế, quá trình v.v... đối tượng nghiên cứu là các thiết bị kỹ thuật gọi là điều khiển học kỹ thuật. Mỗi hệ thống (hoặc phần tử của hệ thống) kỹ thuật, đều chịu tác động của bên ngoài và cho ta các đáp ứng. Gọi tác động vào là đầu vào, tác động ra là đầu ra ( hoặc tín hiệu vào, tín hiệu ra). Hình 1-1 * Nhiệm vụ của lý thuyết điều khiển tự động Lý thuyết điều khiển tự động giải quyết 2 nhiệm vụ chính: - Phân tích hệ thống - Tổng hợp hệ thống Phân tích hệ thống: Nhiệm vụ này nhằm xác định đặc tính đầu ra của hệ sau đó đem so sánh với những chỉ tiêu yêu cầu để đánh giá chất lượng điều khiển của hệ thống đó. Muốn phân tích hệ thống điều khiển tự động người ta dùng phương pháp trực tiếp hoặc gián tiếp để giải quyết 2 vấn đề cơ bản. - Tính ổn định của hệ thống Hệ thống (hoặc phần tử của hệ thống) Các tác động vào Các đáp ứng - Chất lượng của quá trình điều khiển- quá trình xác lập trạng thái tĩnh và trạng thái động (trạng thái quá độ). Để giải quyết vấn đề trên dùng mô hình toán học, tức là các phần tử của hệ thống điều khiển đều được đặc trưng bằng mô hình toán của các phần tử sẽ cho mô hình toán của toàn bộ hệ thống. Có thể xác định đặc tính ổn định của hệ thống qua mô hình toán của hệ thống với việc sử dụng lý thuyết ổn định trong toán học. Tổng hợp hệ thống: Tổng hợp hệ thống là xác định thông số và cấu trúc của thiết bị điều khiển. Giải bài toán này, thực ra là thiết kế hệ thống điều khiển. Trong quá trình tổng hợp này thường kèm theo bài toán phân tích. Đối với các hệ thống điều khiển tối ưu và thích nghi, nhiệm vụ tổng hợp thiết bị điều khiển giữ vai trò rất quan trọng. Trong các hệ thống đó, muốn tổng hợp được hệ thống phải xác định Algorit điều khiển tức là xác định luật điều khiển Đ(t). Hệ thống điều khiển yêu cầu chất lượng cao thì việc tổng hợp càng trở nên phức tạp. Trong một số trường hợp cần đơn giản hoá một số yêu cầu và tìm phương pháp tổng hợp thích hợp để thực hiện. 1.2- Các mô hình diễn tả hệ thống điều khiển. Để tiện việc nghiên cứu về các vấn đề điều khiển cần sử dụng các sơ đồ (mô hình) diễn tả các thành phần của hệ thống sao cho rõ ràng mọi mối quan hệ bên trong và ngoài hệ thống để dễ dàng phân tích, thiết kế và đánh giá hệ thống. Thực tế sử dụng các mô hình sau là phổ biến và thuận tiện: 1) Hệ thống các phương trình vi phân 2) Sơ đồ khối. 3) Graph tín hiệu. 4) Hàm truyền đạt 5) Không gian trạng thái (Sơ đồ khối và Graph tín hiệu là cách biểu diễn bằng đồ hoạ để diễn tả một hệ thống vật lý hoặc một hệ phương trình toán đặc trưng cho các phần tử của hệ thống - Diễn tả một cách trực quan hơn). * Về mặt lý thuyết mỗi hệ thống điều khiển đều có thể diễn tả bằng các phương trình toán. Giải các phương trình này và nghiệm của chúng sẽ diễn tả trạng thái của hệ thống. Tuy nhiên việc giải phương trình thường khó tìm nghiệm (có trường hợp không tìm được) lúc đó cần đặt các giả thiết để đơn giản hoá nhằm dẫn tới các phương trình vi phân tuyến tính thường – Hệ điều khiển tuyến tính liên tục. * Phần lớn kỹ thuật điều khiển hiện đại, là sự phát triển của các mô hình toán học cho các hiện tượng vật lý. Sau đó dựa vào các mô hình toán học để nghiên cứu các tính chất của hệ thống điều khiển. 1.2.1. Phương trình vi phân Các hệ thống vật lý (hoặc các quá trình) cần được diễn tả chính xác mọi quan hệ giữa những đại lượng biến động bên trong của chúng. Từ đó ta dễ dàng nghiên cứu được các hiện tượng diễn biến của hệ thống; các định luật cơ bản của vật lý có thể giúp ta giải quyết vấn đề đó. Các quan hệ của các đại lượng cơ bản nói chung có thể biểu diễn bằng các phương trình vi phân ( gọi là mô hình toán của hệ thống). Ví dụ: Phương trình của định luật II Newton F = m.a Trong phương trình đại số giá trị các đại lượng không thay đổi theo thời gian, vì thế nó chỉ diễn tả trạng thái ổn định của hệ. Nhưng trong thực tế hệ không tĩnh. Đầu ra thường biến động đối với các thay đổi của đầu vào, thêm vào đó tác động của nhiễu cũng thay đổi theo thời gian, nên hệ không ổn định tức là đầu ra dao động. Vì thế cần phải phân tích hệ trong các điều kiện động lực hoặc gọi là trong trạng thái quá độ, lúc này các biến số không cố định mà thay đổi theo thời gian. Phương trình vi phân mô tả hệ ở trạng thái động lực không chỉ chứa bản thân các biến số mà còn chứa tốc độ thay đổi hoặc gọi là đạo hàm của các biến số đó. * Các nội dung cơ bản của phương trình vi phân: Phương trình dạng: an. n n dt yd + an-1. 1n 1n dt yd − − + ... + a1. dt dy + a0. y = x(t) (1.1) x(t) và y(t) là các biến phụ thuộc, t là biến độc lập. * Các tính chất của phương trình vi phân: Mọi hệ là tuyến tính nếu quan hệ vào- ra của nó có thể biểu thị bằng phương trình vi phân tuyến tính: ∑∑ = = i i i n i i i i dt xdb dt yda .. 0 Hoặc một hệ là tuyến tính nếu quan hệ vào ra của nó có thể biểu thị bằng tích phân: y(t) = ∫∞ ∞− τττ dxtW )(),(¦ Trong đó W(t,τ ) là hàm thể hiện các tính chất bên trong của hệ, y(t) là đầu ra và x(t) là đầu vào. Hàm 2 biến W(t,τ ) là hàm trọng lượng của hệ. - Đáp ứng y(t) của một hệ tuyến tính do nhiều đầu vào x1(t), x2(t), ...., xn(t) tác động đồng thời lên hệ bằng tổng các đáp ứng của mỗi đầu vào tác động riêng biệt (nguyên lý chồng chất) y(t) = ∑ = n i i ty 0 )( Ví dụ: Phương trình vi phân thuần nhất: A. dt tdyB dt tyd )(.)(2 2 + + C.y(t) = 0 Có hai nghiệm y1(t), y2(t). theo nguyên lý chồng chất thì y1(t) + y2(t) cũng là một nghiệm của phương trình đó. - Toán tử vi phân và phương trình đặc trưng: Xét phương trình vi phân tuyến tính hệ số hằng cấp n an n n dt yd + an-1. 1 1 − − n n dt yd + ... + a1. dt dy + a0. y = x(t) Gọi toán tử vi phân D = dt d , D n = n n dt d Phương trình trên có thể viết thành: D n y + an-1 D 1−n y + ... + a1Dy + a0y = x (D n + an-1 D 1−n + ... + a1D + a0 )y = x (1.2) Đa thức D n + an-1 D 1−n + ... + a1D + a0 gọi là đa thức đặc trưng. Phương trình D n + an-1 D 1−n + ... + a1D + a0 = 0 là phương trình đặc trưng. Nghiệm của phương trình đặc trưng rất có ý nghĩa khi xét tính ổn định của hệ thống. 1.2.2- Sơ đồ khối. * Sơ đồ khối được biểu thị bằng các khối liên kết với nhau để diễn tả mối quan hệ đầu vào và đầu ra của một hệ thống vật lý. * Sơ đồ khối thuận tiện để diễn tả mối quan hệ giữa các phần tử của hệ thống điều khiển. Ví dụ: a) b) c) Vào A Phần tử G Ra B G1 A G2 B C x d dt y = dt dx Hình 1-2 * Các khối có thể là một thiết bị hoặc dụng cụ và có thể là một hàm (chức năng) xảy ra trong hệ thống. Khối: Ký hiệu thuật toán phải thực hiện đầu vào để tạo đầu ra. Đường nối: Đường nối giữa các khối biểu thị đại lượng hoặc biến số trong hệ thống. Mũi tên: Chỉ tiêu của dòng thông tin hoặc tín hiệu “Các khối nối tiếp nhau thì đầu ra của khối trước là đầu vào của khối sau” Điểm tụ: Biểu hiện thuật toán cộng hoặc trừ ký hiệu bằng một vòng tròn đầu ra của điểm tụ là tổng đại số của các đầu vào. Hình 1-3 * Điểm tán: Cùng một tín hiệu hoặc một biến số phân ra nhiều nhánh tại điểm đó gọi là điểm tán, tức là tại đó đầu ra áp lên nhiều khối khác “ký hiệu là một nốt tròn đen”. Hình 1-4 Cấu trúc sơ đồ khối của hệ thống điều khiển kín Hình 1-5 Hình (1-5) diễn tả một hệ thống điều khiển kín bằng sơ đồ khối. Các khối mô tả các phần tử trong hệ được nối với nhau theo quan hệ bên trong của hệ thống. * Các biến số của hệ: (1) Giá trị vào V: tín hiệu ngoài áp vào hệ. (2) Tín hiệu vào chuẩn R: rút từ giá trị vào V là tín hiệu ngoài hệ áp lên hệ điều khiển như một lệnh xác định cấp cho đối tượng. R biểu thị cho một đầu vào lý tưởng dùng làm chuẩn để so sánh với tín hiệu phản hồi B. x + - y (x-y) x x x x C C C E G1 G2 M C GV V R + H B - u x + + y (x+y) x + + y (x+y-u) - u (3) Biến số điều khiển M (tín hiệu điều chỉnh): là đại lượng hoặc trạng thái mà phần tử điều khiển G1 áp lên phần từ (đối tượng) điều khiển G2 (quá trình được điều khiển). (4) Biến số ra C (tín hiệu ra): là đại lượng hoặc trạng thái của đối tượng (hoặc quá trình) đã được điều khiển. (5) Tín hiệu phản hồi B: là một hàm của tín hiệu ra C được cộng đại số với vào chuẩn R để được tín hiệu tác động E. (6) Tín hiệu tác động E (cũng gọi là sai lệch hoặc tác động điều khiển) là tổng đại số (thường là trừ) giữa đầu vào là R với phần tử B là tín hiệu áp lên phần tử điều khiển. (7) Nhiễu u: là tín hiệu vào không mong muốn ảnh hưởng tới tín hiệu ra C. Có thể vào đối tượng theo M hoặc một điểm trung gian nào đó (mong muốn đáp ứng của hệ đối với nhiễu là nhỏ nhất). * Các phần tử của hệ: (1) Phần tử vào chuẩn GV: chuyển đổi giá trị vào V thành tín hiệu vào chuẩn R (thường là một thiết bị chuyển đổi). (2) Phần tử điều khiển G1: là thành phần tác động đối với tín hiệu E tạo ra tín hiệu điều khiển M áp lên đối tượng điều khiển G2 (hoặc quá trình). (3) Đối tượng điều khiển G2 là vật thể, thiết bị, quá trình mà bộ phận hoặc trạng thái của nó được điều khiển. (4) Phần tử phản hồi H: là thành phần để xác định quan hệ (hàm) giữa tín hiệu phản hồi B và tín hiệu ra C đã được điều khiển (đo hoặc cảm thụ trị số ra C để chuyển thành tín hiệu ra B (phản hồi). (5) Kích thích: là các tín hiệu vào từ bên ngoài ảnh hưởng tới tín hiệu ra C. Ví dụ tín hiệu vào chuẩn R và nhiều u là các kích thích. (6) Phản hồi âm: điểm tụ là một phép trừ E = R - B (7) Phản hồi dương: ở điểm tụ là phép cộng: E = R + B (Điều khiển kín gồm hai tuyến: Tuyến thuận truyền tín hiệu từ tác động E đến tín hiệu ra C. Các phần tử trên tuyến thuận ký hiệu G (G1 , G2, ...) tuyến phản hồi truyền từ tín hiệu ra C đến phản hồi B các phần tử ký hiệu là H (H1 , H2 , ...). 1.2.3. Hàm truyền đạt: Hàm truyền đạt của hệ thống. * Hàm truyền đạt của hệ thống đối với hệ thống điều khiển liên tục một đầu vào và một đầu ra được định nghĩa: - Là tỷ số của biến đổi Laplace của đầu ra với biến đổi Laplace của đầu vào với giả thiết toàn bộ các điều kiện đầu đồng nhất bằng không (điều kiện dừng). G(s) = o1 1n 1n n o1 1m 1m m m aS.a...S.aS bsb...S.bSb ++++ ++++ − − − − (1.3) Đối với hệ thống vật lý thực các chỉ số trong hàm truyền n ≥ m. * Trong lĩnh vực thời gian gián đoạn (điều khiển rời rạc) việc biến đổi Z đóng vai trò của biến đổi Laplace: Hàm truyền có dạng sau: G(z) = o1 1n 1n n o1 1m 1m m m az.a...z.az bzb...z.bzb ++++ ++++ − − − − (1.4) * Đối với hệ thống nhiều đầu vào nhiều đầu ra với r đầu vào, p đầu ra, các hàm truyền là các phần tử của ma trận cấp p×r phần tử , với chỉ số i của phần tử thứ i của đầu vào, chỉ số thứ j của phần tử thứ j đầu ra. G11(s) G12(s) ..... G1r(s) G21(s) G22(s) ..... G2r(s) G(s) = ..... ..... Gji(s) ..... (1.5) ..... ..... ..... ..... GP1(s) ..... ..... GPr(s) Ở đây: Gji(s) = )s(u )s(Y i j ; các đầu vào khác ui(s) đều coi là bằng không. (Nguyên lý độc lập tác dụng). * Một cách tương tự với hệ thống điều khiển gián đoạn ta có hàm truyền của hệ thống nhiều đầu vào nhiều đầu ra. G11(z) G12(z) ..... G1r(z) p×r G21(z) G22(z) ..... G2r(z) G(z) = ..... ..... Gji(z) ..... (1.6) ..... ..... ..... ..... GP1(z) ..... ..... GPr(z) Ở đây: s - số phức - biến Laplace. z = eS.T - biến của phép biến đổi z. 1.2.4. Không gian trạng thái Khi phân tích và thiết kế hệ thống điều khiển tuyến tính thường sử dụng một trong hai hình thức sau: + Đối với lĩnh vực thời gian sử dụng hàm trạng thái . + Trong lĩnh vực tần số dùng hàm truyền đạt. Như ở trên, ta xét hệ phương trình vi phân, sai phân đạo hàm đến bậc n (hệ thống bậc n) ; n thực chất là trạng thái của các biến. Các trạng thái của biến được mô tả như là vectơ x. Các phương trình trạng thái được mô tả dưới dạng sau (hệ thống tuyến tính). . x (t) = A.x(t) + B.u(t) ; x(o) = xo y(t) = C.x(t) + D. u(t) (1.7) và x(k+1) = A. x(k) + B.u(k) ; x(o) = xo y(k) = C.x(k) + D. u(k) (1.8) Ở đây: A, B, C, D là các ma trận hệ số hằng có kích thước. An×n , Bn×r , CP×n , DP×r Các hệ phương trình viết dạng (1-11); (1-12) các phương trình trạng thái của hệ thống điều khiển. * Không gian trạng thái: Một hệ thống có r tín hiệu vào u1(t), u2(t), u3(t) ... ur(t) m tín hiệu ra: y1(t), y2(t), y3(t).... ym(t) Xác định n biến trạng thái: x1(t), x2(t)..... xn(t) Vậy hệ thống được mô tả bởi phương trình không gian trạng thái như sau: =(t)x 1 . f1(x1, x2,..., xn; u1, u2,..., ur; t) . . . =)(. txn fn(x1, x2,..., xn; u1, u2,..., ur; t) Đại lượng ra: y1(t) = g1(x1, x2,..., xn; u1, u2,..., ur; t) . . . ym(t) = gm(x1, x2,..., xn; u1, u2,..., ur; t) + + D(t) C(t) A(t) dt u(t) y(t) B(t) x(t) x(t) . + + =(t)x. ⎥⎥ ⎥⎥ ⎥⎥ ⎥⎥ ⎥ ⎦ ⎤ ⎢⎢ ⎢⎢ ⎢⎢ ⎢⎢ ⎢ ⎣ ⎡ (t)x . . . (t)x (t)x n . 2 . 1 . f(x, u, t) = ⎥⎥ ⎥⎥ ⎦ ⎤ ⎢⎢ ⎢⎢ ⎣ ⎡ t);u,...,u,u;x,...,x,(xf ... t);u,...,u,u;x,...,x,(xf t);u,...,u,u;x,...,x,(xf r21n21n r21n212 r21n211 (1.9) Phương trình trạng thái: =)(. tx f(x, u, t) y(t) = g(x, u, t) Hoặc dưới dạng ma trận: =)(. tx A(t). x(t) + B(t). u(t) y(t) = C(t). x(t) + D(t). u(t) Sơ đồ khối: Hình 1-6 1.3. Mô tả toán học của các phần tử điều khiển a. Phần tử di động thẳng: Tác dụng vào lò xo có chiều dài L0 để lò xo di động một lượng X thì cần một lực: PL = k .X (k: là độ cứng lò xo hay là hằng số lò xo) k = X PL Δ Δ O L P L P L L0 L X K = PLX H×nh 1-7. §¦ êng ®Æc tÝnh H×nh 1-8. S¬ ®å khèi 1/k R=XPV= PL Đối với lò xo thông thường tín hiệu vào là lực PV = PL, tín hiệu ra là lượng di động R = X. Vậy mô hình toán đặc trưng và sơ đồ khối biểu diễn chức năng như hình 1-8 b. Bộ giảm chấn bằng không khí hoặc bằng dầu ép: Để di động piston với vận tốc V, cần tác dụng lực PV có giá trị: PV = C.V= C. dt dR áp dụng toán tử Laplace: s = dt d PV = C.V= C. dt dR = C.s.R Lực PV coi là tín hiệu vào Tín hiệu ra: Lượng di động R. Từ các yếu tố trên thành lập sơ đồ khối thể hiện mô hình toán của bộ giảm chấn. c. Trọng khối Theo định luật II Newton tổng các lực P ở bên ngoài tác dụng vào một trọng khối sẽ có biểu thức: ∑P = M.A = M. 22dtRd Dùng toán tử Laplace: s = dt d nên ∑ P = M. S2.R R = ∑ P . 2. 1 SM Sơ đồ khối thể hiện mô hình toán như sau: d. Phần tử quay Định luật II Newton: Đối với chuyển động quay gia tốc góc của vật thể quay tỷ lệ thuận với tổng mô men tác dụng lên nó. Dạng toán học của định luật: P V R H×nh 1-9 PV R1/C.s H×nh 1-10 H×nh 1-11 1/M.S2 RP 2 2 dt d ϕ = θ M∑ ⇒ M d dt ∑= .2 2 ϕθ Trong đó: ϕ là góc quay θ là momen quán tính của vật thể M là momen bên ngoài tác dụng vào vật thể. Momen bên ngoài được tạo ra từ động cơ, do tải trọng tác dụng lò xo hoặc giảm chấn. Xét một đĩa quay trong chất lỏng và nối với một bánh đà như hình vẽ: -Phân tích để xây dựng mô hình toán: Quay đĩa được phải tác dụng một momen xoắn Mx, trục quay đi một góc là j tạo mo men của lò xo: M1 = kx. j (1.10) Trục có đường kính D, chiều dài l, hệ số lò xo xoắn là: kx = l GD 32 .. 4π (G: Mô đun đàn hồi) Momen cần thiết để thắng lực ma sát của chất lỏng: Mm = C.w = C. dt dϕ = C. p. j (1.11) w: là vận tốc góc C: hệ số ma sát của chất lỏng Nếu quay đĩa với momen xoắn Mx (momen xoắn của trục lò xo) và momen ma sát sẽ ngăn cản sự quay của đĩa do đó có thể viết thành: ∑ M = Mx – M1 – Mm = 2 2 . dt d ϕθ = q. s2. j Thay các trị số (1.10) và (1.11) ta có: Mx = q. s2. j + kx. j + C. s. j = (q. s2 + kx + C.s). j Từ phương trình trên ta có sơ đồ khối của hệ thống như hình vẽ. e. Các phần tử điện Các phần tử cơ bản của các mạch điện H×nh 1-12 M1 Mm Mx ω ϕ H×nh 1-13 1(θ.S2+ C.S + k x) ϕMx y = Rx =V C1 A.P b) uR = R. I ⇒ I = R 1 .uR uL= L. dt dI = LP. I ⇒ dt dI = p. I = dt d .I uC= C 1 . ∫ dtI. = PC 1 .I f.Các phần tử thuỷ khí Xét phần tử dầu ép: -Nếu van trượt được đẩy lên phía trên , dầu có áp suất P0 sẽ vào buồng trên của xi lanh 3 và dầu của buồng dưới sẽ qua van trượt về bể dầu. - Nếu van trượt được đưa xuống phía dưới , dầu sẽ qua buồng dưới của xilanh 1 và dầu ở buồng trên sẽ chảy về bể dầu. Với hiệu áp không đổi được hình thành ở cửa van, tức là tỷ lệ thuận với lượng di động x. Gọi q là lượng dầu chảy vào xilanh, ta có: q = C1.x q đồng thời cũng là sự thay đổi thể tích của xilanh: q = A.Py (A là diện tích bề mặt của xilanh) ⇒ A.Py = C1.x − + u R uL + − uC + − R L C 1 R u R I Iu L Lp 1 IuC Cp 1 H×nh 1-14 x 1 2 3 y P 0 M H×nh 1-15 ⇒ y = PA C . 1 .x Từ phương trình trên tín hiệu vào là x ( lượng di động của xilanh 1) và tín hiệu ra y lượng di động của xilanh 2. g.Phần tử phi tuyến Ta xét một phần tử phi tuyến và trên cơ sở đó tiến hành tuyến tính hoá mô hình toán học đặc trưng cho chức năng của cơ cấu. Xét cơ cấu nâng vuông góc bằng cơ khí: Thanh nâng vuông góc tại điểm A (a + b = 900) và có thể chuyển động cưỡng bức trong rãnh thẳng đứng. Một nhánh của thanh nâng có thể trượt trên con trượt ở điểm B , con trượt này di động cưỡng bức theo phương ngang. Nhánh kia của thanh nâng có thể di động trong bạc của khớp nối cố định ở điểm C. - Phân tích: Tam giác AOB luôn đồng dạng tam giác AOC nên: K X X Y = ⇒ K XY 2 = ( K = const) Nếu tín hiệu vào là X, thì vị trí của điểm B là tín hiệu ra Y tỷ lệ với bình phương của X. Còn tín hiệu vào là Y và tín hiệu ra là X sẽ tỷ lệ với căn bậc hai của Y: X = YK. Để viết phương trình toán và xây dựng mô hình toán học ta cần tuyến tính hoá các phương trình phi tuyến trên. Phương pháp như sau. 1.4- Phân loại hệ thống điều khiển. * Việc phân loại hệ thống điều khiển (Controller System) có rất nhiều hình thức tuỳ theo góc độ nhìn nhận đánh giá: phân loại theo tín hiệu vào, theo các lớp phương trình vi phân mô tả quá trình động lực học của hệ thống. Theo số vòng kín trong hệ, v.v... Tuy nhiên đây chỉ là tương đối. Xét về tính chất làm việc và nội dung cơ bản của điều khiển thì hệ thống điều khiển có 2 loại làm cơ sở trong phân tích tính năng (Phân biệt tác động vào hệ và đáp ứng ra): Hệ thống kín Hệ thống hở. H×nh 1-16 α β O Y K X B C *Theo đặc điểm mô tả toán học thì có các hệ thống sau: Hệ thống liên tục Hệ thống gián đoạn Hệ thống tuyến tính Hệ thống phi tuyến Hệ thống tuyến tính hoá * Theo dạng năng lượng tiêu thụ: Hệ thống điều khiển bằng điện Hệ thống điều khiển bằng dầu Hệ thống điều khiển bằng khí ép .... 1.4.1. Các hệ thống điều khiển hở và hệ thống kín a. Hệ thống điều khiển hở (Open- Loop Control Systems) *Khái niệm: Hệ thống điều khiển hở là hệ thống mà tác động điều khiển độc lập với đầu ra (Hoặc đầu ra không được đo và không được phản hồi so với đầu vào) Ví dụ: Quá trình hoạt động của máy giặt hoàn toàn tự động mà chúng ta chỉ cần tác động trước khi máy hoạt động là đóng điện và nhấn công tắc sau khi máy hoàn thành công việc thì chúng ta lấy sản phẩm ra. Trong máy có diễn ra các quá trình như sau: quá trình làm ướt quần áo (Soaking), quá trình giặt (Washing), quá trình vắt khô (Rinsing) đều làm việc với một thời gian tổng chuẩn (time basic) Và các quá trình này không được đo kết quả (Tức là không được kiểm tra là đã làm sạch quần áo hay chưa) Sơ đồ khối của hệ thống (Control System in Washing Machine) t = ts + tW + tR = const Từ ví dụ trên ta thấy hệ thống điều khiển hở có dáp ứng ra không so sánh đáp ứng vào. Mỗi tác động vào có trạng thái (hoạt động) ổn định, kết quả của hệ thống có độ chính xác phụ thuộc hệ thống chia độ (hệ thống đo). Trong quá trình có nhiễu, hệ thống không thực hiện nhiệm vụ yêu cầu. * Đặc tính của hệ thống điều khiển hở: - Độ chính xác của hệ quyết định bởi điều chỉnh (căn) và có duy trì độ chính xác đó được lâu hay không. H×nh 1-17 Soaking Washing Rinsing Turn on Finish Cleanliness - Nhạy cảm với các biến đổi xung quanh như: nhiệt độ, dao động, xung lực, điện thế, phụ tải... - Đáp ứng chậm khi tín hiệu vào thay đổi. * Ưu điểm: - Đơn giản - Giá thành thấp (Độ chính xác vừa phải) - Vấn đề mất ổn định không nghiêm trọng. b. Hệ thống điều khiển kín Khái niệm: Hệ thống điều khiển kín là hệ thống mà tác động điều khiển phụ thuộc đáp ứng ra. còn gọi là hệ thống điều khiển phản hồi. E: Sai lệch điều khiển E = R – B R: Tín hiệu vào B: Tín hiệu phản hồi. Trong hệ thống điều khiển kín sai lệch điều khiển là sự chênh lệch giữa tín hiệu vào và tín hiệu phản hồi. Quá trình điều khiển nhằm giảm sai lệch và đáp ứng ra đạt giá trị mong muốn. Ví dụ: Hệ thống điều khiển nhiệt độ trong lò là một hệ thống điều khiển kín. Nhiệt độ trong lò điện được đo bởi nhiệt kế ( là thiết bị Analog(tương tự)) Nhiệt độ dưới dạng tín hiệu tương tự được biến đổi thành tín hiệu nhiệt độ dạng số bởi bộ A/D. Tín hiệu nhiệt độ được chuyển về máy tính trung tâm qua Interface. và nhiệt độ được so sánh với tín hiệu nhiệt độ mà chương trình của máy tính đã lập, nếu có bất kỳ sai số nào (discrepancy: sự chênh lệch, sự khác nhau) thì máy tính trung tâm có tín hiệu qua Interface và tín hiệu này được khuếch đại nhờ thiết bị Amplifier và tác động lên Relay làm cho nhiệt độ trong lò tăng hay giảm tuỳ theo yêu cầu của chương trình đã lập. ER + - G1 H H×nh 1-18 G2 C B Lß ®iÖn (E.Furnace) A/D Converter Interface Relay Amplifier Interface Computer Programming input H×nh 1-19. Ví dụ 2: Để điều khiển một bình nước sao cho mực nước trong bình luôn là hằng số không đổi thì độ cao cột nước trong bình sẽ là một trong những thông số kỹ thuật cần quan tâm của hệ thống. Giá trị về độ cao cột nước tại thời điểm t được đo cảm biến và được biểu diễn thành một đại lượng điện áp dưới dạng hàm số phụ thuộc thời gian u(t) có đơn vị Volt. Đại lượng vật lý ở đây là điện áp đã được sử dụng để truyền tải hàm thời gian u(t) mang thông tin về độ cao cột nước. ( Phần mô hình toán học) * Đặc tính của hệ thống điều khiển kín( hệ thống phản hồi) Đặc trưng của hệ thống điều khiển kín là phản hồi. - Nâng cao độ chính xác có khả năng tạo lại đầu ra - Tốc độ đáp ứng nhanh - Độ chính xác phụ thuộc các điều kiện làm việc - Giảm tính chất phi tuyến và nhiễu - Giảm độ nhạy cảm của tỷ số đầu ra và đầu vào đối với sự thay đổi tính chất của hệ. - Tăng bề rộng dải tần (dãy tần số của đầu vào) - Có khuynh hướng dao động hoặc không ổn định. - Điều khiển mềm . 1.4.2.- Các hệ thống điều khiển liên tục và gián đoạn. Các hệ thống thực được mô tả ở trạng thái tĩnh hoặc động lực học. Các hệ thống tĩnh thường được diễn tả bởi hệ thống các phương trình đại số. Trong điều khiển kỹ thuật các hệ thống tĩnh không diễn tả đầy đủ trạng thái của hệ thống. Vì vậy người ta dùng các phương trình vi phân/sai phân mô tả trạng thái động lực học của hệ thống (được biết như là các hệ thống với các tham số cục bộ hoặc tập trung) hoặc các phương trình vi phân đạo hàm riêng (như là các hệ thống có các tham số phân tán). Trong nội dung giáo trình ta nghiên cứu các hệ thống được mô tả bởi hệ các phương trình vi phân/sai phân tuyến tính, nghĩa là các tham số của hệ thống độc lập tuyến tính. Ví dụ hệ thống động lực học được mô tả dưới dạng các phương trình vi phân/sai phân vô hướng: x& (t) = fc(x(t)) , x(to) = xo (1.12) x(k +1) = fd (x(k)) , x (ko) = xo (1.13) Ở đây: t : biến thời gian liên tục. k : biến thời gian gián đoạn. Chỉ số e: (continuous- Time) - thời gian liên tục. d: (discrete - Time) - thời gian gián đoạn. Nếu hệ thống chịu tác động của ngoại lực, hay các tác động vật lý khác. Ta nói nó chịu tải động điều khiển và phương trình vi phân/sai phân mô tả trạng thái động lực của hệ thống. x& (t) = fc (x(t), u(t)) ; x(to) = xo (1.14) x(k+1) = fd (x(k), u(k)) ; x(ko) = xo (1.15) Ở đây: u(t) ; u(k) đóng vai trò biến điều khiển. Với mục đích của điều khiển ta thay đổi biến điều khiển nhận được các đáp ứng của hệ thống kỹ thuật theo yêu cầu như vậy, nhìn chung vấn đề chính của điều khiển có thể mô hình hoá theo dạng sau: tìm biến điều khiển bằng cách giải hệ thống phương trình vi phân đặc trưng của hệ. Nếu các hệ phương trình vi phân (1.12) ÷ (1.15) là tuyến tính ta gọi hệ thống là tuyến tính. Nếu là phi tuyến ta gọi là hệ thống phi tuyến. Việc nghiên cứu hệ thống phi tuyến tương đối khó. Trong thực tế, người ta tìm cách tuyến tính hoá. Trong phạm vi giáo trình này, chúng ta chỉ nghiên cứu hệ thống điều khiển tuyến tính. 1.5- Tuyến tính hoá hệ thống phi tuyến. Trong thực tế không có một hệ thống vật lý nào có thể mô tả tuyệt đối chính xác bằng phương trình vi phân hệ số hằng tuy nhiên nhiều hệ phi tuyến có thể xấp xỉ hoặc coi như tuyến tính trong từng đoạn làm việc. Có nhiều phương pháp được áp dụng cho việc tuyến tính hoá hệ thống phi tuyến. Phương pháp trung bình gần điểm làm việc. Phương pháp tuyến tính hoá điều hoà và phương pháp sai lệch nhỏ. 1.5.1- Phương pháp trung bình gần điểm làm việc. Đây là phương pháp đơn giản được dùng trong thiết kế các hệ thống khi đặc tính trên không thể xấp xỉ hoá được bằng các hàm giải tích. Phương pháp này áp dụng cho các hệ có những phần tử chỉ phi tuyến ở trạng thái tĩnh, quan hệ giữa đầu ra y với đầu vào u là ở trạng thái xác lập (ổn định). Giả thiết trong đoạn: - uM < u < um đặc tính phi tuyến có thể xấp xỉ hoá bằng đường thẳng. Trong đó: y = K . u ; k = m m u y = tgα ; α là độ dốc. 1.5.2- Phương pháp tuyến tính hoá điều hoà. Phương pháp này được dùng khi hệ có một phần tử tuyến tính nối sau một phần tử phi tuyến làm việc ở chế độ tự dao động. Các tín hiệu trong hệ là làm tuần hoàn theo thời gian. Phương pháp này dựa trên cơ sở khai triển hàm sóng thành chuỗi hàm dạng sin (chuỗi Fonricr) điều hoà có tần số là ω, 2ω, 3ω, ... có biên độ và góc pha khác nhau. Giả thiết các hàm điều hoà bậc cao khác (2ω, 3ω, ...) có biên độ nhỏ bỏ qua chỉ giữ lại thành phần điều hoà bậc nhất (ω) (giả thiết lọc) nghĩa là: Hình 1-20 Trong đó: u(t) = Um . sin (ωt + ψ) y(t) = Ym1 . sin (ωt + ϕ) Trong đó Um = Ym1 và ϕ - ψ = π được gọi là điều kiện cân bằng điều hoà. 1.5.3- Phương pháp sai lệch nhỏ. Theo phương pháp này việc tuyến tính hoá được thực hiện bằng cách khai triển hàm phi tuyến thành chuỗi Taylor tại vùng lân cận điểm ổn định (tương ứng với chế độ xác lập). Chỉ khảo sát các sai lệch bậc nhất trong chuỗi đó. Sai lệch so với trạng thái ổn định càng nhỏ thì việc đánh giá các quá trình của phần tử phi tuyến có sai số càng bé sau khi biến đổi tuyến tính. a) Hệ thống (bậc nhất) phi tuyến. x& (t) = f(x(t) , u(t) ) (1.16) Giả thiết rằng hệ thống làm việc ở trạng thái xác lập với quĩ đạo xn(t) khi nó được điều khiển bởi tín hiệu vào un(t). Chúng ta gọi xn(t) và un(t) là quĩ đạo danh nghĩa và đầu vào danh nghĩa theo phương trình (1.16) ta có: x& n(t) = f(xn(t) , un(t) ) (1.17) Bây giờ ta giả thiết rằng thay đổi của hệ phi tuyến (1.16) lân cận quĩ đạo danh định một lượng nhỏ (vô cùng bé). x(t) = xn(t) + Δx(t) (1.18) Lượng biến đ...m lại từ bước 1 đến bước 5 đối với mỗi tín hiệu vào. 2.3. Graph tín hiệu và qui tắc Mason 2.3.1. Graph tín hiệu Các hệ thống điều khiển còn được mô tả bằng mô hình toán là Graph tín hiệu. Graph tín hiệu thể hiện bằng đồ thị sự truyền tín hiệu trong hệ thống, nhưng dễ dàng hơn các dạng mô hình toán khác. Xét phương trình đơn giản: Xi = Aij. Xj Các biến Xi, Xj : là hàm thời gian, hàm biến phức hoặc hằng số, hoặc là hằng số. NótNót Nh¸nh Xj XiAij Ain XiX2 X1 Xn Ai2 Ai1 Aij là một toán tử ánh xạ Xj vào trong Xi nên Aij gọi là hàm truyền ( hàm truyền đạt). Khi Xi, Xj các hàm của biến Laplace S ( biến phức). Mỗi biến số trong Graph được Mỗi biến số trong Graph được kí hiệu bằng một nút mỗi hàm chuyển được ký hiệu bằng một nhánh, các nhánh đều có hướng ký hiệu bằng mũi tên diễn tả dòng tín hiệu. Hình 2.35 * Quy tắc hội tụ ( cộng vào): Tổng các tín hiệu đi vào một nút bằng giá trị các nút đó. Tổng quát: Xi = ∑ = n 1j jij .XA Hình 2.36 *Quy tắc phân kỳ ( chuyển ra): Giá trị của một nút có thể chuyển ra từng nhánh rời khỏi nút đó. Nếu ta có: Xi = Aik ; i = 1,2,..., n. Thì Graph như hình vẽ: AjkXk X2 X1 Xn A2k A1k Xn Ank Hình 2.37 * Quy tắc nhân: Nhiều nhánh nối tiếp nhau có thể thay bằng một nhánh có hàm chuyển bằng tích các hàm chuyển của các nhánh đó. A21 Xn-1X2X1 Xn An(n-1) X1 Xn A21.A21...An(n-1) = A43 X4X1 X2 X3 A21 A33 A32 A23 A23 A32A21 X3X2X1 = X1 X2 X3 A21 A32 A23 1 X4 X3=X4 Xn = A21. A32. A43... An(n-1) .X1 Hình 2.38 * Các thành phần trong Graph tín hiệu: Cho một Graph tín hiệu như hình vẽ sau Hình 2.39 - Một tuyến: Là một trình tự nối tiếp, đơn hướng của các nhánh, trong đó không có nút nào bị xuyên qua quá một lần. X1 đến X2 đến X3 đến X4 X2 đến X3 và trở về X2 X1 đến X2 đến X4. - Nút vào: Là một nút chỉ có các nhánh đi khỏi nó ( X1). - Nút ra: là một nút chỉ có các nhánh đi tới nó ( X4) Có thể thêm một nút giả với hàm chuyển bằng 1 để thoả mãn định nghĩa này. Hình 2.40 - Tuyến thuận: là tuyến đi từ nút vào đến nút ra ( bằng bất cứ đường nào) X1 đến X2 đến X3 đến X4; X1 đến X2 đến X4 - Tuyến phản hồi: là tuyến xuất phát và kết thúc tại cùng một nút. X2 đến X3 đến X2 - Tuyến đơn: Là tuyến phản hồi chỉ có một nhánh. - Hai tuyến ( hoặc hai vòng kín) gọi là không chạm nhau nếu chúng không có nút chung. - Hàm truyền của tuyến hoặc của vòng kín bằng tích hàm truyền của các nhánh nằm trong tuyến hoặc vòng kín đó. Tuyến thuận X1 đến X2 đến X3 đến X4 có gia lượng A21. A32. A43 Tuyến phản hồi: X2 đến X3 đến X2 có gia lượng A32. A23 x2 x1 x3x3 dt 2 2 d -1 dt d 1 x1 A32 A23 A21 A31 A33 x4A43 A42 x1 x2 x3 x4 A21 A23 A31 A32 A42 A43 A33 x2 x3 G1 G2 G3 G4 G5 G7 G6 -H1 -H2 u y1 1 Các ví dụ: Ví dụ 1: Dựng Graph tín hiệu cho hệ thống được mô tả bởi phương trình vi phân sau: x3 = 112 2 2 x dt dx dt xd −+ Từ phương trình ta thấy có 3 biến số x1, x2, x3 nên cần có 3 nút ( không kể nút giả). Các toán tử trong phương trình là dt d và 2 2 dt d Viết lại phương trình trên: x3 = 12 2 x dt d dt d −+ )()( 12 xx Sơ đồ Graph tín hiệu: Hình 2.41 Ví dụ 2: Dựng Graph tín hiệu cho nhóm phương trình xét đồng thời sau: x2 = A21. x1 + A23. x3 x3 = A31. x1 + A32. x2 + A33.x3 x4 = A42. x2 + A43. x3 Nhận xét: Phương trình trên có 4 biến số x1, x2, x3, x4 ta có sơ đồ Graph tín hiệu sau Hình 2.42 2.3.2. Quy tắc Mason Từ sơ đồ Graph tín hiệu có thể rút gọn sơ đồ và tìm hàm truyền đạt của cả hệ thống. Để tìm hiểu về quy tắc Mason ta có ví dụ minh họa sau: Hình 2.43 Bước 1: Xác định tất cả những tuyến thẳng Pk có thể có của hệ thống. Đó là những đường nối liền nhau không chứa đường phản hồi đi từ điểm nút nguồn u(t) tới điểm nút đích y(t) và Pk có giá trị bằng tích các giá trị các đường nối có trong Pk. Hệ trên có 3 tuyến thẳng: P1 = G1. 1. G2. G7 P2 = G1. 1. G6. G4. G5 P3 = G1. G2. G3. G4. G5 Bước 2: Xác định tất cả những vòng lặp Lk có thể có của hệ thống. Đó là những đường nối liền nhau tạo thành một vòng kín. Hệ trên có 4 vòng lặp: L1 = -1. G4. H1 L2 = -1. G2. G3. G4. G5. H2 L3 = -1. G6. G4. G5. H2 L4 = - 1. G2. G7. H2 Bước 3: Tính ....L.LL.LLL1Δ nm,l, nml ji, ji k k +−+−= ∑∑∑ (2.3.2.1) Trong đó: Li, Lj là những cặp hai vòng lặp không trùng nhau ( không có chung một nhánh nào) Ll, Lm, Ln là bộ 3 vòng lặp không trùng nhau,... Hệ trên chỉ có 2 vòng lặp L1, L2 là không trùng nhau ( không có đoạn nào giống nhau). ....L.LL.LLL1 nm,l, nml ji, ji k k +−+−= ∑∑∑Δ = 1 – ( L1 + L2 + L3 + L4) + L1. L4 = = 1 + G4. H1 + G2. G3. G4. G5. H2 + G6. G4. G5. H2 + G2. G7. H2 Bước 4: Xác định Δ k từ Δ bằng cách trong công thức (2.3.2.1) ta bỏ đi tất cả những vòng lặp có đoạn nối chung với Pk . Tức là: Δ 1= 1 – L1 = 1 + G4. H1 ( tất cả các vòng lặp đều không có đoạn nối chung với P1) Δ 2 = 1 ( tất cả các vòng lặp đều có đoạn nối chung với P2 ( có G1) Δ 3 = 1 ( Vòng lặp có đoạn chung với P3 ) Bước 5: Xác định hàm truyền đạt G(s) theo công thức Mason: G(s) = ).Δ(P Δ 1 k kk∑ Vậy G(s) = Δ 1 . ( P1. Δ 1 + P2. Δ 2 + P3. Δ 3) = - u G2 y G3G1 - H2 H1 1u G3G2G1 -H2 H1 -1 y h1 h2 A1 A2q u(t) y(t) r1 r2 γ,p1 γ,p2 = H .G .G H .G .G .G H .G .G .G .G H .G 1 G .G .G .G .G G .G .G 1. .G)HG.(1G .G 1. .G 27225462543214 54321546114721 ++++ +++ Ví dụ 1: Cho hệ thống có sơ đồ khối như sau, sơ đồ Graph tín hiệu tương đương như hình vẽ Hình 2.44 Hệ chỉ có một tuyến thẳng đó là: P1 = G1. G2. G3 Hệ có 3 vòng lặp từng đôi một có đoạn nối chung: L1 = G1. G2. H1 L2 = -G2. G3. H2 L3 = -G1. G2. G3 Vậy, ....L.LL.LLL1 nm,l, nml ji, ji k k +−+−= ∑∑∑Δ = 1 – ( L1 + L2 + L3) = = 1 - G1. G2. H1 +G2. G3. H2 + G1. G2. G3 Do tất cả các vòng lặp cũng đều có tuyến thẳng P1 nên Δ 1= 1 Hàm truyền của hệ thống là: G(s) = ).(P1 k kk∑ ΔΔ = Δ1 . P1. Δ 1 = 321232121 321 G .G .G H .G .G H .G .G - 1 G .G .G ++ Ví dụ 2: Xét một hệ thống gồm 2 bình chứa chất lỏng như sau Hình 2.45 A1s 1 γ 1 r1 1 A2s γ r2 1u(t) y(t)h2 p2qp1h1 1 A1s γ γr11 1A2s r21 -1 -1 -1 u(t) y(t) Chất lỏng được bơm vào bình thứ nhất với lưu lượng u(t). Nếu chất lỏng trong bình thứ nhất có độ cao h1, áp suất p1, hệ số chuyển đổi áp suất, lưu lượng r1, hệ số áp suất, độ cao g. lưu lượng chảy sang bình thứ hai là q và h2, p2, r2 là độ cao, áp suất, hệ số chuyển đổi áp suất, lưu lượng của chất lỏng trong bình thứ 2. Theo các định luật vật lý, giữa những thông số kỹ thuật đó có quan hệ: A1. dt dh1 = u(t) – q q = . r 1 1 (p1 – p2) A2. dt dh2 = q – y(t) y(t) = . r 1 2 p2 ( áp suất tại đầu ra được xem như bằng 0) p1 = γ.h1 p2 = γ.h2 Trong đó y(t) là lưu lượng chất lỏng chảy ra khỏi bình thứ 2. Từ những hiểu biết lý thuyết ban đầu đó của hệ thống ta có sơ đồ khối và sơ đểu Graph mô tả tín hiệu mô tả hệ thống. Hình 2.46 Từ sơ đồ trên ta thấy hệ chỉ có một tuyến thẳng: P1 = 2 2121 2 .s.A.A.rr γ Hệ có 3 vòng lặp: L1 = - .s.Ar 11 γ L2 = - .s.Ar 21 γ L3 = - .s.Ar 22 γ Trong đó có 2 vòng lặp L1 và L2 không có nhánh nào chung. Nên HÖ thèng Ku(t) y(t)K ....L.LL.LLL1 nm,l, nml ji, ji k k +−+−= ∑∑∑Δ = 1 – ( L1 + L2 + L3) + L1.L3 = 1 + ( .s.Ar 11 γ + .s.Ar 21 γ + .s.Ar 22 γ ) + .s.Ar 11 γ . .s.Ar 22 γ = 2 2121 221211 2 2112 .s.A.A.rr .Ar.Ar.Ar.s.A.r.Ar 2)..( γγ ++++ s Vì cả 3 vòng lặp trên đều có nhánh nối chung với P1 nên 11 =Δ Vậy hàm truyền đạt: G(s) = Δ .ΔP 11 = 2 2121 2 .s.A.A.rr γ . 2)..( γγ ++++ s22121122112 2 2121 .Ar.Ar.Ar.s.A.r.Ar .s.A.A.rr = = 2)..( γγ γ ++++ s22121122112 2 .Ar.Ar.Ar.s.A.r.Ar 2.4. Các hệ thống lấy mẫu dữ liệu Như đã biết, hệ thống liên tục là hệ có các biến số vào và ra được truyền đi và biến đổi liên tục theo thời gian, có thể quan sát vào bất cứ thời điểm nào. Nhưng trong điều khiển còn có nhiều hệ thống mà các biến số chỉ được đưa vào và xử lý gián đoạn, nó cho đáp ứng tại các thời điểm gián đoạn đó là các hệ thống rời rạc mà các tín hiệu truyền đi không liên tục. Có các dạng hệ thống gián đoạn: - Các hệ thống lấy mẫu gián đoạn từ các hệ liên tục, biến đổi tín hiệu liên tục thành gián đoạn gọi là lượng tử hoá - Các hệ thống làm việc theo chu kỳ - Các hệ thống có cấu trúcc chu kỳ Hệ rời rạc, gián đoạn có những ưu điểm: - Làm việc ít tốn năng lượng, có tính kinh tế - Có thể điều khiển nhiều kênh đồng thời, chống nhiễu tốt - Truyền và giữ tin được lâu - Về lý thuyết không cần phép tính tích phân và vi phân nên đơn giản hơn - Có nhiều tính chất giống như hệ liên tục - Mô hình toán là các phương trình lặp ( phục hồi lại) * Mô hình toán của hệ thống rời rạc Xét hệ xung lấy mẫu gián đoạn: Hình 2.47 0 y(kT) t T 2T 3T 4T 5T 5T4T3T2TT 1 2 3 4 554321 k y(k) 0 f) Chuçi rêi r¹c 0 u(k) k u(1) u(2) u(3) u(4) u(5) y(1) y(2) y(3) y(4) y(5) e) Chuçi rêi r¹c t u(kT) 0 τ TÝn hiÖu vμo b) c) TÝn hiÖu ra τ Đóng và mở bộ ngắt K theo chu kỳ để mạch của nó không liên tục được nữa; ta sẽ được các xung gián đoạn liên tiếp nhau tạo thành một chuỗi tín hiệu xung. Mỗi xung kéo dài một thời gian t . Giả sử thao tác bộ ngắt K sao cho t càng nhỏ ( t → 0) với một chu kỳ lấy mẫu cố định T thì các xung càng thu hẹp lại và ta chọn tỷ lệ thời gian sao cho chu kỳ lấy mẫu T = 1, tức là u(kT) = u(k) y(kT) = y(k) k = 0, 1, 2, 3,... là các số nguyên Phương trình lặp đại số có dạng sau: any(k+n) + an-1y(k + n -1) +... + a1y(k + 1) + a0 y(k) = bm u( k+m) + bm-1u(k + m -1) + ... + b1 u(k + 1) + b0 u(k) Trong phương trình trên không có vi phân, cũng không có tích phân gọi là phương trình lặp lại để diễn tả hệ rời rạc ( lấy mẫu) tương đương với phương trình vi phân của hệ liên tục. k: là biến độc lập với các giá trị 0, 1, 2, 3,... u(k): là một chuỗi rời rạc mô tả tín hiệu vào y(k): là một chuỗi rời rạc khác mô tả tín hiệu ra. Hình 2.48 * Toán tử gián đoạn: Hệ thống gián đoạn cũng quy định một vài toán tử với hàm cần tìm. - Toán tử cộng thêm 1: E (k) = k + 1 f(k+1) f(k) k k+1 f(k) k f(k) 0 E[f(k)] = f( k+ 1) E f(k) f(k+1) E[f(k)] = f(k+1) k+mk E m n E k+m+n Tính chất của toán tử E: Tính lặp lại: En(k) = k + n Tính nghịch đảo: E-n(k) = k – n Tính gộp: Em. En = Em+n E[Cf(k)] = C.f(k+1) = C. E[f(k)] ; C là hằng số E[f1(k) + f2(k)] = f1( k+1) + f2(k + 1) = E[f1(k)] + E[f2(k)] E[C1.f1(k) + C2.f2(k) + ... + Cn. fn(k)] = ∑∑ == =+ n 1i ii n 1i ii (k)].E[fC1)(k.fC Hàm truyền đạt: Phương trình lặp: any(k+n) + an-1y(k + n -1) +... + a1y(k + 1) + a0 y(k) = = bm u( k+m) + bm-1u(k + m -1) + ... + b1 u(k + 1) + b0 u(k) ⇔ an.En[y(k)] + an-1.En-1[ y(k)] +... + a1. E[y(k)] + a0 y(k) = bm. Em[u(k)] + + bm-1. Em-1[u(k)] + ... + b1. E[u(k)] + b0. u(k) Ta có : D(E).y(k) = N(E) u(k) Trong đó: D(E) = an.En + an-1.En-1 +... + a1. E + a0 N(E) = bm. Em + bm-1. Em-1 + ... + b1. E + b0 Hàm truyền của hệ thống: H(E) = D(E) N(E) = 0 1n 1n n n 0 1m 1m m m a...EaEa b...Eb.Eb +++ +++ − − − − *Toán tử sai phân Δ [f(k)] = f(k+1) – f(k) Hình 2.49 Các tính chất: 1. Δ [Cf(k)] = Cf(k+1) – C f(k) = C Δ f(k) 2. Δ [f1(k) + f2(k)] = [ f1(k+1) + f2(k+1)] – [ f1(k) + f2(k)] = Δ f1(k) - Δ f2(k) 0 f(t) t T 2T 3T 4T 5T t y(kT) 0 f(T) f(2T) f(3T) f(5T) 0 y(kT) t 5T4T3T2TT f(0) f(kT)f(t) y(kT) 3.Δ[C1f1(k) + C2f2(k) + ... + Cnfn(k)] = ∑∑∑ === Δ=−+ n 1i ii n 1i ii n 1i ii (k)f.C (k).fC1)(k.fC 4.Δ[f(k)] = f(k+1) – f(k) = E f(k) – f(k) = ( E – 1)f(k) ; Δ = E – 1 5.Δ[f(k)] = f(k+1) – f(k) 6.Δ2[f(k)] = Δ[Δf(k)] = Δf(k+1) - Δf(k) 7.Δn[f(k)] = Δ[Δn-1f(k)] = Δn-1f(k+1) - Δn-1f(k) Δn[f(k)] = (E-1)n [f(k)] = [ En - ∑ = −−−− −+−=+−+−++−+−+=+−+ ++−=+−++−+ n 0r r n rr n rr n r 1nnrnr n r2n1n r)nf(kC1)(...r)nf(kC1)(....1)nnf(kn)f(k ...[f(k)]C1)( ...[f(k)]E 1! n [f(k)]E...].f(k)E.C1)(....E 2! 1)n(n .E 1! n Trong đó: r)!r!.(n n! C rn −= 8.Δmn[f(k)] = Δm.Δn[f(k)] = Δn. Δm[f(k)] Còn dùng toán tử sai phân ngược: ∇f(k) = f(k) – f( k- 1); ∇ = 1 – E-1 * Biến đổi Z: Trong các hệ tuyến tính liên tục ta đã dùng biến đổi Laplace; các hệ này có tính nhân quả ( tích phân một phía từ 0 đến ∞ ); Với định nghĩa biến đổi Laplace của hàm f(t) là: L {f(t)} = ∫ ∞ + − 0 ).( dtstetf Giữa biến đổi Laplace của hệ tuyến tính liên tục và biến đổi Z của hệ tuyến tính rời rạc có mối liên quan chặt chẽ. Phép nhân 2 tín hiệu có thể thực hiện được bằng một sự biến điệu (modulation). Các xung lấy mẫu gián đoạn từ một tín hiệu liên tục f(t) có thể xem như kết quả của một sự biến điệu của một chuỗi rời rạc Y(kT) theo biên độ của tín hiệu liên tục và có cùng chu kỳ với chu kỳ lấy mẫu gián đoạn. y(k) 1 2 3 4 5 k 0 1 d(k) d(k-1) d(k-2) d(k-3) d(k-4) d(k-5) Hình 2.50 Ta có : f(kT) = f(t). Y(kT) Chuỗi rời rạc Y(kT) có thể là chuỗi Kronecker hoặc chuỗi Dirac. Chuỗi Kronecker nhân quả ( 1 phía) là một chuỗi xung có biên độ bằng 1, tác động tại các giá trị bằng 0 hoặc nguyên dương của k ( k ≥ 0): Y(kT) = ∑∞ = − 0j j)Td(k , với T = 1 ta có chuỗi như hình vẽ sau: Hình 2.51 Y(k) = ∑∞ = − 0j j)d(k Trong đó: d( k-j) = ⎩⎨ ⎧ ≠ = j knÕu; 0 j k nÕu; 1 với j là số nguyên, j ≥ 0 Ta có các xung lấy mẫu gián đoạn của f(t): f(kT) = f(t). Y(kT) – f(t) ∫∑ ∞ −∞ = 0 st 0 skT- dtf(t).ef(kT).e Đặt esT = Z với T = 1 nên Z = es Vậy biến đổi Z của chuỗi f(k). Y(k): Z[f(k).Y(k)] = F(Z) = ∑∞ =0k k-f(k).Z Với k = 0,1, 2, 3,... Biến đổi Z là một thuật toán qua đó một chuỗi đưa vào f(k) sẽ cho ra một chuỗi vô tận f(k).Z-k. Tức là, nếu f(k) = [ f(0), f(1), f(2), ..., f(k),...] Z[f(k)] = f(0), f(1).Z-1, f(2). Z-2,.... Biến đổi Z của một chuỗi f(k) chỉ tồn tại nếu chuỗi biến đổi Z đó tuyệt đối hội tụ: [Z] > Rc : Chuỗi hội tụ [Z] < Rc : Chuỗi phân kỳ. Rc là bán kính của vòng tròn hội tụ ( Tâm của vòng tròn là gốc của mặt phẳng phức. Các điểm ở ngoài vòng tròn diễn tả tính hội tụ của chuỗi F(k), các điểm trong vòng tròn diễn tả tính không hội tụ, còn các điểm nằm trên đường tròn là điểm đặc biệt phải xét riêng. 2.5. Hàm truyền đạt của hệ thống rời rạc a) Đáp ứng xung của hệ rời rạc 0 d(k-k0) k k0 §iÒu kiÖn ®Çu b»ng 0 k0 k h(k,k0) 0 k HÖ thèng Cũng như ở hệ liên tục, đáp ứng xung của hệ rời rạc, vô hướng ( 1 biến) là đáp ứng của hệ đó với xung Dirac ( Kronecker). Hệ lúc đầu ở trạng thái dừng, điều kiện đầu bằng 0 x(0) = 0; tín hiệu xung Dirac áp lên hệ vào thời điểm k0T và đáp ứng được xét ở thời điểm kT ( k > k0). Đáp ứng xung là một chuỗi [ h(kT, k0T)] Giả sử T = 1: Hình 2.52 Vì vậy với một hệ rời rạc, vô hướng , nhân quả, tuyến tính, hệ số cố định ta có: Tín hiệu vào là xung d( k- j ) thì tín hiệu đầu ra là đáp ứng xung h ( k- j). Tín hiệu vào là chuỗi u(k) = ∑∞ = − 0j j)d(k (j). u thì tín hiệu ra là chuỗi y(k) = ∑∞ = − 0j j) h(ku(j). Do đó: y(k) = ∑∞ = − 0j j) h(ku(j). với i, j = 0, 1, 2, 3 ,... y(k) = u(k). h(k) Nếu u(k) = d(k) là tín hiệu xung đơn vị ( Dirac hoặc Kronecker) thì đầu ra đáp ứng xung là h(k). “Như vậy đáp ứng của một hệ rời rạc, vô hướng , nhân quả, tuyến tính, bất biến là tích chập của tín hiệu vào u(k) với đáp ứng xung h(k) của hệ đó”. b) Hàm truyền của hệ rời rạc y(k) = u(k). h(k) Ta có biến đổi Z là: Z[y(k)] = Z[h(k)]. Z[u(k)] Y(Z) = H(Z). U(Z) H(Z) là hàm truyền của hệ rời rạc, đó là biến đổi Z của đáp ứng xung h(k) của hệ. H(Z) = U(Z) Y(Z) 2.6. Ứng dụng MatLab Nhập mô hình của hệ thống điều khiển trong MatLab: G(S) = o1 1n 1n n o1 1m 1m m m aS.a...S.aS bsb...S.bSb ++++ ++++ − − − − = )( )( Sden Snum >> num = [ bm bm-1 . . . b1 b0 ] >> den = [ 1 an-1 . . . a1 a0 ] >> sys = tf (num, den) Kết quả: Tranfer function )( )( Sden Snum Hoặc: >> S = tf(’s’) >> G(s) = )( )( Sden Snum Mô hình điểm không - điểm cực: >> [z, p, k] = residue (num, den) >> z = zero (sys) >> [p,z] = pzmap (sys) ( Hiển thị đồ thị cực - không) >> p = pole (sys) Tìm nghiệm mẫu số của hàm truyền đạt: >> c = [ 1 an-1 . . . a1 a0 ] >> p = roots (c) Đồ thị đáp ứng của hệ thống điều khiển: >> impulse (num, den,t) >> Step (num, den, t) >> Lsim (num,den,u, t) Chuyển từ hàm truyền đạt sang phương trình trạng thái: >> [A, B, C, D] = tf2ss (num, den) Và ngược lại: >> [num, den] = ss2tf ( A,B, C, D) Chuyển sang mô hình hệ thống gián đoạn: >> sysd = c2d (sys , Ts) >> sysc = d2c (sysd) Đáp ứng của hệ thống gián đoạn: >> dimpulse (num, den) >> Dstep(num, den) >> dlsim(num, den) Chương 3 PHƯƠNG TRÌNH TRẠNG THÁI * Đặt vấn đề: - Các hệ thống tuyến tính liên tục được mô tả bởi hệ n phương trình vi phân cấp một mô tả n trạng thái của hệ thống mô hình toán hệ thống viết dưới dạng ma trận. x& (t) = A. x(t) + B. u(t) ; xo = x(o) (3-1) và y(t) = C. x(t) + D. u(t) (3-2) Ở đây: x ∈ ℜn, u ∈ ℜr, y ∈ ℜP tương ứng là các vectơ trạng thái, các đầu vào, các đầu ra. Ma trận hệ số An×n mô tả các mối liên hệ bên trong hệ thống. Các ma trận Bn×r , CP×n , DP×r , đặc trưng cho mối liên hệ với bên ngoài của hệ thống. Nếu không có đường dẫn trực tiếp giữa các đầu vào với đầu ra thì DP×r là ma trận zero. * Mô hình không gian trạng thái của hệ thống điều khiển gián đoạn (tuyến tính) là các phương trình sai phân. x(k+1) = Ad . x(k) + Bd.u(k) , x(o) = xo (3-3) y(k) = Cd x(k) + Ddu(k) (3-4) 3.1- Các mô hình không gian trạng thái. Mô hình không gian trạng thái của hệ thống động lực học liên tục đều có thể diễn tả hệ thống trong lĩnh vực thời gian bằng các phương trình vi phân hoặc hàm truyền dưới bốn dạng (form) sau: - Dạng điều khiển (không gian pha). (Controller canonical form). - Dạng quan sát (không gian quan sát). (observer canonical form). - Dạng modal (không gian modal). (Modal canonical form). - Dạng Jordan (không gian Jordan). 3.2- Mô hình không gian trạng thái và các phương trình vi phân. Hệ thống động lực học cấp n được mô tả bằng phương trình vi phân cấp n. n n dt )t(yd + an-1 1n 1n dt )t(yd − − +... + a1 dt )t(dy + aoy(t) = = bn n n dt )t(ud + bn-1 1n 1n dt )t(ud − − +... + b1 dt )t(du + bou(t) (3-5) Ta giả thiết các điều kiện đầu của hệ thống y(o-) , dt ody )( − , ... , 1n 1n dt )o(yd − − đồng thời bằng không, ta tiến hành biến đổi phương trình vi phân cấp n thành hệ n phương trình vi phân cấp 1. + Xét phương trình vi phân cấp n sau: n n dt )t(yd + an-1 1n 1n dt )t(yd − − + ...+ a1 dt )t(dy + aoy(t) = u(t) (3-6) Đổi biến theo: x1(t) = y(t) x2(t) = dt )t(dy x3(t) = 2 2 dt )t(yd (3-7) . . . . . . . . xn(t) = 1n 1n dt )t(yd − − Tiến hành lấy đạo hàm hai vế các phương trình (3-7). dt )t(dx1 = x& 1 = dt )t(dy = x2(t) dt )t(dx 2 = x& 2 = 2 2 dt )t(yd = x3(t) (3-8) . . . . . . . . )t(d )t(dx n = x& n = n n dt )t(yd = - ao(y(t) - a1 dt )t(dy - ... - an-1 1n 1n dt )t(yd − − + u(t) = - aox1(t) - a1x2(t) - ... - an-1 xn(t) + u(t) Vậy không gian trạng thái (3-8) viết dưới dạng ma trận x& 1 0 1 0 L L 0 x1(t) 0 x& 2 0 0 1 L L 0 x2(t) 0 x& 3 M M O O O M M M M = M M M O O 0 × M + M × u(t) x& n-1 0 0 L L 0 1 xn-1(t) 0 x& n -ao -a1 L L L -an-1 xn(t) 1 (3-9) Đầu ra được viết theo (3-7). y(t) = [1 0 0 ... ... 0] × [x1(t) x2(t) ... xn(t)]T (3-10) (3-9) và (3-10) được gọi là dạng chính tắc của không gian pha. + Đối với hệ thống được mô tả bởi phương trình (3-5) ta có: y(t) = [(bo - aobn) (b1-a1bn) ... ... (bn-1 - an-1)] × [x1(t) x2(t) ... xn(t)]T + bn u(t) (3-11) Với: bn = 0 ta có: y(t) = [bo b1 ... ... bn-1] × [x1(t) x2(t) ... ... xn(t)]T (3-12) 3.3- Xác định các biến trạng thái từ hàm truyền. Phần này giới thiệu các kỹ thuật hình thành mô hình không gian trạng thái từ hàm truyền của hệ thống thường được áp dụng trong thực tế. Đó là kỹ thuật chương trình trực tiếp và kỹ thuật chương trình song song. Để đơn giản ta xét với hệ thống một đầu vào một đầu ra. 3.3.1- Mô phỏng HT theo dạng điều khiển chính tắc. Kỹ thuật này được sử dụng thuận lợi khi hàm truyền của thiết bị dạng đa thức không phân tích ra thừa số được. )s(U )s(Y = o1 1n 1n n o1 1n 1n n n aSa...SaS bSb...SbSb ++++ ++++ − − − − (3-13) Ở đây ta sử dụng biến phụ V(s). (Tính điều khiển được của hệ thống là với một tác động vào liệu có chuyển được trạng thái của hệ từ thời điểm đầu to đến thời điểm cuối trong khoảng thời gian hữu hạn không?). )s(V )s(Y = bnSn + bn-1Sn-1+ ... + b1S + bo (3-14a) )( )( sU sV = o1 1n 1n n aSa...SaS 1 ++++ −− (3-14b) Sơ đồ khối mô tả hệ thống có sử dụng biến phụ V(s). Hình 3.1 Phương trình (3-14a) được viết lại như sau: Y(s) = bnSnV(s) + bn-1Sn-1V(s) + ... + b1S.V(s) + boV(s) (3-15) Điều này chỉ ra rằng y(t) là sự chồng chất của V(t) và các đạo hàm của nó vì ta có thể trình bày (3-14a, b) dưới dạng phương trình vi phân khi điều kiện đầu đồng nhất bằng không bằng cách thay: S ≡ dt d ; Si ≡ i i dt d ; V(s) ≡ v(t), V(s)/U(s) Y(s)/V(s) U(s) V(s) Y(s) Xây dựng mô hình không gian trạng thái của hệ thống từ các hàm truyền bằng cách sử dụng sơ đồ mô phỏng rất thuận tiện. Trong các trường hợp hệ thống liên tục sơ đồ mô phỏng các máy tính tương tự giải các phương trình vi phân mô tả các hệ thống động lực học sử dụng các bộ tích phân, bộ cộng bộ trừ và nhân được thực hiện như là bộ khuếch đại thuật toán. Số khối tích phân phụ thuộc vào cấp của phương trình vi phân. + Sơ đồ mô phỏng (3-14a, b) như sau: Sử dụng kỹ thuật chương trình trực tiếp: đặt n khối tích phân nối tiếp với đầu vào tương ứng là V(n)(t) , v(n-1)(t) , ..., V(1)(t) , v(t). Áp dụng (3-15) xác định y(t) bằng cách nhân đầu vào vi(t) với các hệ số bi và cộng bằng bộ cộng ∑ . + Từ (3-14b) ta có: v(n)(t) = u(t) - an-1 v(n-1)(t) - ... - a1v(1)(t) - aov(t) (3-16) Các phép trừ mô phỏng bằng mối liên hệ ngược trên sơ đồ ta có: Hình 3-2: Sơ đồ mô phỏng kỹ thuật chương trình trực tiếp (dạng điều khiển chính tắc). Theo hình 3-2 ta có mô hình không gian trạng thái của hệ thống dạng điều khiển chính tắc. 0 1 0 L L 0 0 0 0 1 L L 0 0 x& (t) = M M O O O M × x(t) + 0 × u(t) M M M O O 1 M -ao -a1 -a2 L L -an-1 1 ∑ 1/S K 1/S 1/S bo ∑ U(t) x& n V(n-1) xn x& 2 x2 x& 1 x1 y(t) b1 b2 bn -an-1 V(n) -a1 -ao v(1) v(0) (3-17) Và y(t) = [(bo-aobn) (b1 - a1bn) ... (bn-1 - an-1 bn)] x(t) + u(t) . bn (3-18) ** (Để chuyển từ hàm truyền sang dạng không gian trạng thái trong Matlab sử dụng hàm tf2ss). Ví dụ: Cho hàm truyền. G(s) = S811S1213185.97S463S996.0S 1908065.90S576S331.0S65.1 23456 234 +++++ ++−− a) Vẽ sơ đồ mô phỏng hệ thống (dạng điều khiển). b) Dựng mô hình không gian trạng thái của hệ thống. 3.3.2- Mô phỏng HT theo dạng quan sát chính tắc. Cùng với dạng điều khiển, dạng quan sát chính tắc là quan hệ quan trọng đối với lý thuyết điều khiển hiện đại. * Quan sát được của một hệ thống là với các toạ độ đo được ở biến ra y(t) của hệ thống liệu ta có thể khôi phục được các vectơ trạng thái x (t) trong thời gian hữu hạn không? Không gian trạng thái của hệ thống và dạng quan sát chính tắc của nó được xác định có cấu trúc rất đơn giản. Xuất phát từ hàm truyền (3-13) ta có: Y(s)(Sn + an-1Sn-1 + ...+ a1S + ao) = U(s) (bnSn + bn-1Sn-1 + ...+ b1S + bo) (3-18) Y(s) = - nS 1 (an-1.Sn-1 +...+ a1S + ao) Y(s) + nS 1 . U(s) (bnSn + + bn-1Sn-1 + ...+ b1S + bo) (3-19) Khai triển ra ta có: Y(s) = - an-1 . S 1 Y(s) - an-2 2S 1 Y(s) - ... - a1 1nS 1 − Y(s) - ao nS 1 Y(s) + + bnU(s) + bn-1 S 1 U(s) + ...+ b1 1nS 1 − U(s) + bo nS 1 U(s) Mối quan hệ (3-20) được thể hiện trên sơ đồ mô phỏng qua n tầng tích phân. Nhãn của các tín hiệu đặt quá tầng tích phân ví dụ một khối S 1 chỉ một tầng tích phân. Tín hiệu an-2 y(t) và bn-2U(t) chỉ vượt qua hai tầng tích phân, aoy(t) và bo u(t) vượt qua n tầng tích phân. Hình 3-3: Sơ đồ khối mô phỏng dạng quan sát chính tắc. + Các biến trạng thái như là đầu ra của các khối tích phân quan hệ đầu ra với các biến trạng thái theo sơ đồ trên ta có: Y(t) = xn(t) + bnu(t) (3-21) x& 1(t) = - aoy(t) + bou(t) = - aoxn(t) + (bo - aobn) u(t) x& 2(t) = - a1y(t) + b1u(t) + x1 = x1(t) - a1xn(t) + (b1 - a1bn) u(t) x& 3(t) = - a2y(t) + b2u(t) + x2 = x2(t) - a2xn(t) + (b2 - a2bn) u(t) x& n(t) = - a1-1y(t) + bn-1u(t) + xn-1 = = xn-1(t) - an-1xn(t) + (bn-1 - an-1bn) u(t) (3-22) Từ (3-21) và (3-22) ta dễ dàng viết dưới dạng ma trận của dạng quan sát chính tắc: 0 0 L L -ao bo - aobn 1 0 L L -a1 b1 - a1bn x& (t) = 0 1 L L -a2 x(t) + b2 - a2bn u(t) M 0 1 L L . 0 0 0 L -an-2 . 0 0 0 1 -an-1 bn-1 - an-1bn (3-23) và y(t) = [ 0 0 ... ... 0 1] x(t) + bn u(t) (3-24) Ví dụ: G(s) = S11,8S12131S8,94S463S996,0S 19080S6,90S576S331,0S65,1 23456 234 +++++ ++−− Hãy viét dạng quan sát chính tắc dưới dạng ma trận. 3.3.3- Kỹ thuật mô phỏng chương trình song song. Đối với kỹ thuật này ta phân ra làm hai trường hợp: đa thức mẫu có nghiệm thực riêng biệt và có nghiệm lặp. + bo - ao + b1 - a1 + bn-1 - an-1 1/S 1/S 1/S L 1/S + bn U(s) y(s) x& 1 x& n-1 x2 xn-1 x& n xn a) Đa thức mẫu có hàm truyền, có nghiệm riêng biệt. Dạng không gian trạng thái này thuận tiện cho các ứng dụng kiểu này bắt nguồn từ việc khai triển hàm truyền thành tổng các phân thức. Một cách tổng quát m < n thì: )s(U )s(Y = )pS)...(pS)(pS( )s(P n21 m +++ = 1 1 pS r + + 2 2 pS r + + ... + n n pS r + + k (3-25) Ở đây p1 , p2 , ..., pn các nghiệm riêng biệt (các cực) của đa thức mẫu của hàm truyền. - Sơ đồ khối mô phỏng dạng này như sau: (Dạng modal chính tắc) Hình 3-4: Sơ đồ khối mô phỏng kỹ thuật lập trình song song. Mô hình không gian trạng thái theo sơ đồ khối này như sau: -p1 0 L L 0 1 0 -p2 L L 0 1 x& (t) = 0 1 L L L x(t) + M u(t) (3-26) L L L L 0 M 0 0 0 0 -pn 1 y(t) = [ k1 k2 ... ... kn ] x(t) (3-27) b) Đa thức mẫu có nghiệm lặp. ∑ y(t) + 1/S r2 x & 2 x2 - p2 + 1/S r1 x& 1 x1 - p1 + 1/S rn x& n xn - pn M u(t) Khi hàm truyền có cực thực lặp. Giả thiết cực p1 lặp r lần. )s(U )s(Y = )pS)...(pS()pS( )s(N n1r r 1 +++ + Dạng khai triển của nó là: )s(U )s(Y = 1 11 pS k + + 21 12 )pS( k + +...+ r1 r1 )pS( k + + 1r 1r pS k + + + +...+ n n pS k + Hình 3-5: Sơ đồ mô phỏng dạng Jordan chính tắc. 3.3.4- Các mô hình của hệ thống gián đoạn. (tương tự trong sơ đồ chỉ thay khối S 1 ≡ Z-1 ). 3.4. Xác định hàm đáp ứng từ phương trình trạng thái 3.4.1. Hệ thống điều khiển liên tục Phương trình trạng thái của hệ theo (3.1) và (3.2) Nghiệm của (3.1): x(t) = eAt. x(0) + ∫t At duBe 0 )(.. ττ (3-28) y(t) = C. eAt. x(0) + C. ∫t At duBe 0 )(.. ττ + D. u(t) (3-29) y(t) = yqđ(t) + yôđ(t) + 1/S kr+1 x& r+1 xr+1 - p1 M u(t) + 1/S kn - pn + 1/S k1r+1 x& r xr - p1 K + 1/S - p1 x& 2 x2 + 1/S - p1 x& 1 k1r ∑ y(t)x1 k1r+1 Đáp ứng quá độ: Là đáp ứng của hệ thống không phụ thuộc vào kích thích u(t) mà do các điều kiện đầu của hệ ( trạng thái ban đầu). Gọi là dao động tự do của hệ thống. Đáp ứng ổn định: Đáp ứng phụ thuộc vào u(t). Đặc trưng cho quá trình cưỡng bức của u(t) làm cho hệ thống ổn định. 3.4.2. Hệ thống điều khiển gián đoạn Phương trình trạng thái được biểu diễn ở (3-3) và (3-4) Nghiệm của phương trình (3-3): x(k) = Ak-k0. x(k0) + ∑− = −−1 1 0 )(.. n kj jk juBA (3-30) y(k) = C. Ak-k0. x(k0) + C.∑− = −−1 1 0 )(.. n kj jk juBA + D.u(k) (3-31) 3.4.3. Các phương pháp tìm đáp ứng Tìm ma trận trạng thái: eAt - Toán tử Laplace: Áp dụng công thức: eAt = l-1 [ (SI - A)-1] - Phương pháp Sylvester: Dựa vào trị riêng của ma trận A: Tìm trị riêng bằng cách tìm nghiệm của phương trình sau det (λI - A) = 0, giải phương trình được các nghiệm: λ1, λ2, ..., λn. Ta có: eAt = α0(t).Ι + α1(t).Α + α2(t).Α2+ ... + αn(t).An Trong đó: Các hệ số α0(t), α1(t), α2(t), ..., αn(t) xác định từ hệ phương trình sau α0(t) + α1(t).λ1 + α2(t).λ 21 + ... + αn(t).λ 11 −n = eλ 1 t α0(t) + α1(t).λ2 + α2(t).λ 22 + ... + αn(t).λ 12−n = eλ 2 t . . . α0(t) + α1(t).λn + α2(t).λ 2n + ... + αn(t).λ 1−nn = eλ n t Chương 4 ỔN ĐỊNH CỦA HỆ THỐNG TUYẾN TÍNH 4.1- Khái niệm chung. Các chương II và III đã trình bày mô tả toán học của hệ thống của hệ thống điều khiển truyền động. Chương này sẽ sử dụng các tư liệu của các chương đã trình bày trước đây để giải quyết nhiệm vụ đầu tiên khi phân tích hệ thống điều khiển tự động là xác định tính ổn định của nó. Thực ra việc nói một hệ thống ổn định là nói đến một số đại lượng nào đó được điều khiển ổn định. Một hệ thống thường biểu diễn bằng phương trình vi phân tổng quát: n n dt )t(yd + an-1 1n 1n dt )t(yd − − + ... + a1 dt )t(dy + aoy(t) = = bm m m dt )t(ud + bm-1 1m 1m dt )t(yd − − + ... + b1 dt )t(du + bou(t) (4-1) Hoặc phương trình sai phân: y(k+n) + an-1y(k+n-1) + ... + a1y(k+1) + aoy(k) = = bmu(k+n) + bm-1u(k+m-1) + ...+ b1u(k+1) + bou(k) (4-2) Sẽ bao gồm hai quá trình: Quá trình xác lập và quá trình quá độ. Đặc trưng bằng nghiệm: y(t) = yo(t) + yqđ(t) (4-3) Hoặc y(k) = yo(k) + yqđ(k) (4-4) Trong đó: - yo là nghiệm riêng của (4-1) hoặc (4-2) đặc trưng cho quá trình xác lập. - yqđ là nghiệm tổng quát của (4-1) hoặc (4-2) khi không có vế phải đặc trưng cho quá trình quá độ. Quá trình xác lập là một quá trình ổn định vấn đề chỉ còn xét quá trình quá độ yqđ. 4.2- Khái niệm ổn định và các định nghĩa chính. Đối với