L thuyt KT chuyn thi c Ngi vit

Lý thuyết ĐKTĐ chuyện thi cử Người viết: Bùi Trung Hiếu Ngành Điều khiển tự động Khoa Điện-Điện tử Trường ĐHBK tp Hồ Chí Minh

Lời thưa: Như đã biết, với Matlab, công việc học tập môn ĐKTĐ trở nên rất đơn giản và thú vị. Tuy nhiên, để đối phó với kì thi, dù bạn là một người học rất tốt lý thuyết nhưng không chú trọng đến cách làm bài vẫn có thể bị điểm thấp. p Đã một lần bị như thế, tôi đành phải bỏ ra một khoảng thời gian để có thể thích nghi với công việc tất nhiên của SV: thi cử! Trong bài này, tôi trình bày với các bạn 2 bài toán rất cơ bản của lý thuyết ĐKTĐ. v Vẽ biểu đồ Bode. v Thiết kế một khâu rời rạc. Tất nhiên, chúng sẽ được trình bày để giải với Caculator, tôi sử dụng FX 570 MS. p

Vẽ giản đồ Bode với sự trợ giúp của FX 570 MS Với Matlab, công việc này rất đơn giản dùng dòng lệnh: bode(hàm_truyền) với hàm truyền đã được khai báo dưới dạng: v Hàm_truyền=tf(tử_số, mẫu_số) v Hàm_truyền=zpk(zero, cực, độ_lợi) …các thông số phụ q Tuy nhiên, để đối phó với kì thi, bạn phải vẽ được biểu đồ Bode dùng Caculator, lý thuyết trong sách ĐKTĐ đã hướng dẫn các bạn một cách rất chi tiết, tôi chỉ nêu cách các bạn dùng Caculator để tính ra các kết quả chú ý: p

Vẽ giản đồ Bode biên độ với sự trợ giúp của FX 570 MS Lý thuyết được trình bày chi tiết trong sách LT ĐKTĐ nhà xuất bản ĐHQG tp Hồ Chí Minh trang 112 -113. p v Bước 1: Xác định tất cả các tần số gãy chúng theo trật tự tăng dần v p và xếp Bước 2: Dùng FX 570 MS ở Mode 2 (CMPLX) nhập hàm truyền cần khảo sát, chú ý thay ω bằng A*i. Việc khảo sát sẽ cần các tần số gãy, ta chỉ đơn giản thay chúng để kiểm soát việc vẽ đúng hay sai. Xét ví dụ sau để làm rõ điều đó:

Vẽ giản đồ Bode biên độ với sự trợ giúp của FX 570 MS (ví dụ 1) p Vẽ biểu đồ Bode biên độ gần đúng của hệ thống có hàm truyền: (đề thi học kì 2 năm 2004 -2005) v Bước 1: Các tần số gãy (lưu ý hàm mũ không ảnh hưởng tới Bode biên độ) Nhận xét: Hàm truyền có một khâu tích phân lý tưởng

Vẽ giản đồ Bode biên độ với sự trợ giúp của FX 570 MS (ví dụ 1) v Bước 2: Dùng Caculator (FX 570 MS): Nhấn Mode→ 2(CMPLX) tức Complex v Nhập số liệu như sau: (Tính bằng đơn vị d. B) 20 log. Abs(100(Ai+10)÷(Ai+1)÷(Ai+100)) (Tôi lấy A là biến số trong ví dụ trên) Việc tiếp theo, bạn nhấn phím Calc, sau đó thay A bằng các giá trị Calc: A? 0. 01→’=’→ 60 Calc: A? 1→’=’→ 17 Calc: A? 10→’=’→-17 Calc: A? 100→’=’→-42 Calc: A? 1000→’=’→-80

Vẽ giản đồ Bode biên độ với sự trợ giúp của FX 570 MS (ví dụ 1) v Nối các điểm trên lại, bạn sẽ có được bản đồ Bode biên độ cần vẽ. v Kết quả như sau:

Màu xanh: Bode dùng Matlab Màu đỏ : Vẽ xấp xỉ các giá trị

Vẽ giản đồ Bode pha với sự trợ giúp của FX 570 MS q v Trên, tôi đã trình bày cách vẽ Bode biên độ, còn cách vẽ Bode pha, sẵn đây tôi cũng xin dẫn ra: Trong giáo trình ĐKTĐ, ta không thấy hướng dẫn cách vẽ Bode pha, một mặt vì các bước tiến hành của nó hơi rắc rối, mà kết quả cho cũng không thật chính xác (Không đảm bảo sai số như Bode biên độ(<=3 d. B)). Cách vẫn hay thường làm với SV là thế các giá trị của tần số và lấy giá trị tương ứng. Thường, SV lấy Arc. Tan hàm thích hợp, với lý thuyết từ trang 112 Sđd. Thật ra, ta có thể dùng chức năng Argument của FX 570 MS để tính toán.

Vẽ giản đồ Bode pha với sự trợ giúp của FX 570 MS q v Cơ sở lý thuyết: Sử dụng FX 570 MS, Mode→ 2 (CMPLX). Nhập hàm số dạng:

Vẽ giản đồ Bode pha với sự trợ giúp của FX 570 MS q Lưu ý: Phải phân rã các tích số của hàm truyền thành các tổng Arg. Chuyển Arg của hàm mũ về dạng: v v q v Xét ví dụ sau để rõ hơn: Vẽ biểu đồ Bode pha của hàm truyền cho ở Ví dụ 1

Vẽ giản đồ Bode pha với sự trợ giúp của FX 570 MS (Ví dụ 2) q Nhập hàm khảo sát như sau: (đơn vị: độ) arg(Ai+10)-0. 5 A*180/pi-arg(Ai)-arg(Ai+100) v Công việc tiếp theo: i Calc Calc Calc i i i i A? A? A? 0. 01 0. 5 1 2 10 20 50 100 → -90. 8 → -98. 05 → -128 → -158 → -200 → -421 → -697 → -1580 → -3000

Màu xanh: Bode dùng Matlab Màu đỏ : Vẽ xấp xỉ các giá trị

Vẽ giản đồ Bode pha với sự trợ giúp của FX 570 MS q v v v Các sai lầm thường gặp: Khi làm việc với hàm mũ, không đưa về dạng –Tω mà vẫn để ở dạng arg(exp(-Tjω)), bạn lưu ý là hàm mũ ở miền số phức có chu kì tuần hoàn là 2Π. Khi làm việc với hàm mũ quên chuyển đơn vị từ rad sang deg. Không phân rã các tích của hàm truyền thành các tổng Arg khi gặp một số hàm bất thường gây sai kết quả. q Việc tính toán khi thế các giá trị rất nhanh, không mất nhiều thời gian, chỉ khoảng 3’-5’ để hoàn thành bài Bode pha.

Kết luận về phương pháp vẽ Bode: Như đã phân tích trên, việc vẽ Bode không phải là công việc quá khó khăn, tuy nhiên, nó là công việc nặng về toán học, và nếu bạn đã rất thành thục với cách vẽ bằng tay, dùng lý thuyết trang 112 -113(Sđd), khi bạn xác định các khâu của hàm truyền, có thể bạn sẽ không mất thời gian khi vẽ Bode biên độ(khoảng 3’), tuy nhiên, với bạn không cần xác định các khâu của hàm truyền, vẫn có thể vẽ được mà không có trục trặc gì, bạn nên nhớ rằng cách vẽ mà yêu cầu đề chỉ là vẽ gần đúng. Thật ra, tôi có thể vẽ chính xác đến 99% giản đồ Bode nếu dùng FX 570 MS. Lúc đó, tôi đang làm lại cái công việc của một chiếc máy, và tôi chợt nghĩ, làm tốt công việc của một cái máy thì có gì phải tự hào

Khảo sát hệ rời rạc dùng PP Kg trạng thái: Đây cũng là một dạng bài tập rất thường hay gặp trong các kì thi, tôi cũng thử dùng FX 570 MS để giải: Tôi lấy ví dụ trong sách giáo trình để đơn giản: (trang 280 -281) r(t) +- e(t) T e. R(t) C(t) e(k. T) a. Thành lập hệ phương trình trạng thái mô tả hệ thống trên. b. Tính đáp ứng của hệ đối với tín hiệu vào là hàm nấc đơn vị (nhân quả)

Khảo sát hệ rời rạc dùng PP Kg trạng thái: Bước 1: Một cách máy móc, ta tìm được các ma trận A, B, D mô tả hệ liên tục: C(s) E (s) R 1. Chú ý rằng hệ phương trình mô tả hệ liên tục có dạng:

Khảo sát hệ rời rạc dùng PP Kg trạng thái: Bước 2: Tính ma trận quá độ. 1. Cách tìm ma trận nghịch đảo:

Khảo sát hệ rời rạc dùng PP Kg trạng thái: Với Aij là det một ma trận vuông cấp n-1 được tính bằng cách: Bỏ đi hàng i, cột j

Khảo sát hệ rời rạc dùng PP Kg trạng thái: Tính 1. Cách tính Laplace ngược: Giả sử phương trình cần biến đổi có dạng tử số/mấu số. Với bậc tử số bé hơn bậc mẫu số. 2. Ta chia ra 3 trường hợp sau: 1. Nghiệm mẫu số là nghiệm đơn. 2. Nghiệm mẫu số là nghiệm thực kép. 3. Nghiệm mẫu số là nghiệm phức liên hợp.

Khảo sát hệ rời rạc dùng PP Kg trạng thái: 1. Nghiệm mẫu số là nghiệm thực đơn. .

Khảo sát hệ rời rạc dùng PP Kg trạng thái: 2. Nghiệm mẫu số là nghiệm thực kép:

Khảo sát hệ rời rạc dùng PP Kg trạng thái: 3. Nghiệm mẫu số là nghiệm phức liên hợp: s 1, 2= -α±jβ

Khảo sát hệ rời rạc dùng PP Kg trạng thái: 1. Tôi lấy một ví dụ đơn giản để làm rõ vấn đề này:

Khảo sát hệ rời rạc dùng PP Kg trạng thái: 1. Trở lại ví dụ của bài toán rời rạc, bây giờ ta tính các ma trận Ad, Bd, Dd : 1. Bước 3: Một cách máy móc, ta thế Ad = Ф(T) 1. Với lưu ý rằng, ta có thể dùng chức năng tích phân của FX 570 MS 1. Bước 4: Hệ pt trạng thái mô tả hệ thống rời rạc với tín hiệu vào r(k. T)

Khảo sát hệ rời rạc dùng PP Kg trạng thái: 1. Dùng chức năng ma trận của FX 570 MS:

Khảo sát hệ rời rạc dùng PP Kg trạng thái: 1. Ta đã tìm được phương trình trạng thái: Nhập ma trận A là ma trận vừa tìm được. Gán ma trận C=(0 0)T sau đó gán ma trận Ans là ma trận C. (rõ ràng, có thể dùng tới ma trận cấp 3, tức đạo hàm cấp 3) 1. Ta cứ nhấn nút ‘=‘đến chừng nào k=n theo yêu cầu bài toán. 2. Đến đây, ta đã giải xong bài toán khảo sát hệ rời rạc dùng pp kg trạng thái, tuy các bước tính đều có thể dùng Caculator, nhưng bạn phải thực hành vài lần thì mới mong có thể làm không sai sót. 1. Và tôi xin nhắc lại rằng, đây chỉ là biện pháp đối phó với thi cử, bạn cần phải có sự tìm hiểu thích hợp khi dùng Matlab. Hiện tại, tôi sử dụng chương trình matlab 7. 0 và Mathematica 5. 0, các chương trình trên phục vụ rất tốt công việc tính toán.

Nếu bạn có hứng thú trao đổi với tôi về phương pháp dùng Caculator để giải các bài tập trong thi cử các môn ở trường ĐH xin gởi mail về địa chỉ hộp thư: buitrunghieu@khvt. com, rất vui lòng để trao đổi kinh nghiệm sử dụng với bạn.