- Bài tập MATLAB ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 1 Bài 1 : chương trình bảng cửu chương % lap bang cuu chuong tu 2 den 9 clear con = 'y. - end SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 2 case'6' for i = 1:10 fact = 6. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 3 fact = fact*i end % end of program Bài 3 : dung lenh for tinh tong day so s=1+3+5+..+n % dung lenh for tinh tong day so s=1+3+5+..+n n = input('positive integer. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 4 operator = input('Enter operator s. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 5 disp(c). - end SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 13 i = i + 1. - lop(i) end % end of program Ví Dụ Như Sau : enter STT :1 enter name :Hung enter diem :9 continue thanh vien khac y/n :y enter STT :2 enter name :Binh enter diem :8 continue thanh vien khac y/n :y enter STT :3 enter name :Dat enter diem :8 continue thanh vien khac y/n :n ans = STT: '1' name: 'Hung' diem: 9 hang: 'Gioi' ans = STT: '2' name: 'Binh' diem: 8 hang: 'Kha' ans = STT: '3' name: 'Dat' diem: 8 hang: 'Kha' SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 14 Bài 20 : vẽ đồ thị sin(x) va cos(x) trên cùng 1 hệ trục tọa độ %Begining of program x = -2*pi:pi/100:5*pi. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 15 Bài 21 : Vẽ sin(x) và cos(x) trên 2 hệ tọa độ xy nhưng cùng trong 1 khung %Begining of program x = -2*pi:pi/100:2*pi. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 16 Bài 22 : Vẽ sin(x. - legend ('sinc(x)',-1) subplot(2,2,4) plot(x,y4,'linewidth',1.0,'linestyle','.','color','g') grid SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 17 xlabel('x') ylabel('y') title('y4 = 1 - sin(x. - subplot(2,2,2) plot(x,y2,'linewidth',1.0,'linestyle','.','color','b') grid xlabel('x') ylabel('y') SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 18 title('y2 = cos(x. - A = pi A e+000 SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 59 So sánh kết quả khi lần lượt thực hiện các dòng lệnh trên cho thấy : kết quả bằng nhau. - 2/2*3 ans = 3 SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 60 b. - X= 0 0 4 SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 61. - 2+3i 7-4i 9+5i ] z i i i SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 62 - m ng 2 chiều: Thực hiện các dịng lệnh sau lưu lại kết qu. - ans i 7.2000 SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 63 VI .Th c hành một s hàm với mảng trên cửa sổ l nh. - L = trace(z) L i SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 64. - truy cập phần tử thứ 1 đến 3 của hàng 2 và hàng 3 % của ma trận A y SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 65 c. - D SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 66 a. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 67. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 68 BÀI T P TH C HÀNH S 2 Th c hành l p trình với m-file Th c hành các l nh trên cửa sổ soạn thảo m-file. - End of program SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 69 Hay ta cĩ thể gõ lần lượt các lênh sau trước d u nhắc của cửa số commandWindow. - giai pt: x^2 + 2*x -3 = 0 x1 = SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 70 -3 x2 = 1 b. - End of program Save lại với tên là quadratic2.m , trong cửa sổ command window gõ lệnh: SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 71. - End of program Chạy chương trình trên cho ta kết qu như hình dưới dây: SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 72 c. - chunhat(a,b) c = total(a,b) function c = total(a,b) c =sqrt( a.^2 + b.^2) %End of function SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 74 Chạy chương trình và cho ta kết qu sau. - ans = Dai Hoc DL Cong Nghe Sai Gon SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 75 Tiếp tục th c hi n. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 76 - Tạo mảng cấu trúc. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 77 BÀI T P TH C HÀNH S 3 Xử lý đồ họa I. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 78 2. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 79 chương trình thực hiện yên cầu trên. - grid on SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 80. - End of program SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 81 Chương trình vẽ 2 đồ thi y = sin(x) và z = x.*sin(x) trên 2 khung trên cùng một cửa sổ Figue: %Begining of program x = 0:pi/100:2*pi. - lech cho phep hien thi ten va mau do thi ve % End of program SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 82 Chương trình vẽ 4 đồ thi trên 4 khung trên cùng một cửa sổ Figue: %Begining of program x = -2*pi:pi/100:2*pi. - End of program SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 83 4. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 84 II. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 85. - End of program Chạy chương trình cho ta kết qu : Cách 2: đánh tr c tiếp các l nh trong cửa sổ Command Window: SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 86. - Với giao diện đề nghị như sau: SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 87 Trước tiên ta tạo giao diện guide như hình dưới đây: Để có được hình như trên ta cần tạo các đối tượng sau: Edit Text ( số lượng 3 ) STT Fontsize String Tag Color SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 88 1 12 0 Hesoa White 2 12 0 Hesob White 3 12 0 Hesoc White Static Text ( số lượng 8) STT Fontsize String Tag color 2 1 12 Giai pt: ax + bx +c Text1 Tùy ý 2 12 a Text2 Tùy ý 3 12 b Text3 Tùy ý 4 12 c Text4 Tùy ý 5 12 Bỏ trống Ketquax1 White 6 12 Bỏ trống Ketquax2 white 7 12 X1 = Text7 Tùy ý 8 12 X2 = Text8 Tùy ý Push Button ( số lượng 1) STT Fontsize String Tag Color 1 12 Giai PT Giai_pt Tùy ý Sau khi tạo song các đối tượng trên, ta nên save nó với tên ( Giaipt.fig. - set(handles.ketquax1,'string',x(1)) set(handles.ketquax2,'string',x(2)) SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 89. - end SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 90 if nargout [varargout{1:nargout. - set(handles.ketquax1,'string',x(1)) set(handles.ketquax2,'string',x(2)) SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 91 %En of proram. - bcos(x) Text1 Tùy ý 2 12 He so a Text2 Tùy ý 3 12 He so b Text3 Tùy ý SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 92 Axes ( số lượng 1) STT Xgrid Ygrid Tag NextPlot 1 on on DoThi replacechilden Push Button ( số lượng 1) STT Fontsize String Tag Color 1 12 BIEU DIEN DoThi Tùy ý 2 12 EXIT Close Tùy ý Cũng tương tự như bài 1, Sau khi tạo song các đối tượng trên, ta nên save nó với tên ( DotThi.fig. - End of proram Sau khi soạn song đoạn chương trình đó, ta save và chạy chương trình và nhập lần lượt các hệ số a và b vào, nhấn vào nút BIEU DIEN cho ta kết quả như hình dưới : SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 93 SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 94 Khi nhấn vào nút EXIT lập tức cửa sổ giao diện Guide sẽ đóng lại Chương trình M-file tạo ra giao diện trên là. - end % End initialization code - DO NOT EDIT SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 95. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 96 plot(x,y,'b. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 97 BÀI T P TH C HÀNH S 4 Simulink ứng dụng trong các ngành kỹ thu t I. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 98 Lưu mơ hình lên đĩa với lệnh file-save từ cửa sổ model với phần đuơi mặc định là .mdl. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 99 Xu t phát của việc xây dựng các mơ hình hệ thống từ các phương trình vi phân tương ứng là các khối tích phân (Integrator). - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 100 Đến đây việc xây dựng mơ hình xe t i với ngõ vào u và ngõ ra v coi như hồn thành. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 101 Khi ta nhập các giá trị và chạy mơ phịng được kết qu là hình dưới. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 102 IV. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 103 Sau khi l y ra các khối, ta thay đổi tên và thơng số cho khối cũng như đường kết nối giữa các khối như hình. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 104 Mơ hình Simulink hồn chỉnh của hệ điều khiển vị trí motor DC Gán các giá trị vào cho các tham số và chạy mơ phỏng cho ta kết qua như hình dưới. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 105 Gán các giá trị vào cho các tham số và chạy mơ phỏng cho ta kết qua như hình dưới. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 106 1. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 107 BÀI THÍ NGHIỆM 1 ỨNG DỤNG MATLAB KHẢO SÁT HỆ THỐNG DÙNG GIẢN ĐỒ BODE Câu .1. - nhập hàm truyền G2 Transfer function: SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 108 s. - 2 s^5 + 20 s^4 + 70 s^3 + 124 s^2 + 120 s SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 109. - grid on Kết quả ta có hình như sau: SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 110 SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 111 b. - Hệ thống trên có ổn định không, giải thích. - HD: Hệ thống trên có ổn định. - Vẽ đáp ứng quá độ của hệ thống trên với đầu vào hàm nấc đơn vị trong khoảng thời gian t = 0 ÷ 10s để minh họa kết luận ở câu c HD: Để vẽ đáp úng của hệ kín ta nhập các lệnh sau: SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 112. - title('ve dap ung qua do') Kết quả có được: ta thấy hệ thống ổn định. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 113 Câu .2.b. - Dựa vào biểu đồ Bode, tìm tần số cắt biên, độ dự trữ pha, tần số cắt pha, độ dự trữ biên của hệ thống HD: Kết quả cần tìm như hình vẽ: SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 114 Tần số cắt biên wc = 0.249 dB, tần số cắt pha = -179 0C Độ dự trữ biên GM. - HD: Hệ thống trên không ổn định. - Gk = feedback(G,1) Transfer function: 400 SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 115. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 116 BÀI THÍ NGHIỆM 2 ỨNG DỤNG MATLAB KHẢO SÁT HỆ THỐNG BẰNG PHƯƠNG PHÁP QUỸ ĐẠO NGHIỆM SỐ Câu .1. - Vẽ QĐNS của hệ thống. - Lúc này, giá trị K sẽ hiển thị lên như trên hình vẽ sau: Kgh = Ka = 42.9 SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 117 b. - Ke = 185 SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 118 Câu .2. - Thí nghiệm: Với hệ thống như ở phần III.1 : a. - G = tf(TS,MS) Transfer function: SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 119 1. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 120 HD. - 25% SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 121 c. - title('ve dap ung qua do he kin voi K = 185') Kết quả có được: Ta thấy Txl sec SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 122 d. - tạo hà truyền hệ kín Gk2 SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 123 Transfer function: 185. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 124 BÀI THÍ NGHIỆM 3 ỨNG DỤNG SIMULINK MÔ PHỎNG VÀ ĐÁNH GIÁ CHẤT LƯỢNG HỆ THỐNG Câu .1. - Khảo sát mô hình hệ thống điều khiển nhiệt độ: Câu.1.a. - Giá trị của hàm nấc từ 0 1 tương ứng công suất cung cấp 0% 100% SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 125 Transfer Fcn - Transfer Fcn1 : mô hình lò nhiệt tuyến tính hóa. - Mô phỏng và vẽ quá trình quá độ của hệ thống trên. - Cửa sổ Parameters hiện ra, nhấp chuột vào trang Data history và tiến hành cài đặt các thông số như hình bên dưới: SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 126 Tiến hành chạy mô phỏng lại để tín hiệu lưu vào biến ScopeData. - grid on %ke luoi SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 127 Lúc này cửa sổ Figure hiện ra với hình vẽ giống như hình vẽ ở cửa sổ Scope. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 128 Câu .2. - Thí nghiệm: Xây dựng mô hình hệ thống điều khiển nhiệt độ ON-OFF như sau: Trong đó. - Khảo sát quá trình quá độ của hệ thống với các giá trị của khâu Relay theo bảng sau: SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 129 SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 130 Với Switch on = +1 , Switch off = -1 ta có kết quả sau: SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 131 Với Switch on = +5 , Switch off = -5 : SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 132 Với Switch on = +10 , Switch off = -10, ta có kết quả sau : SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 133 Với Switch on = +20 , Switch off = -20, ta có kết quả sau : b. - Tính sai số ngõ ra so với tín hiệu đặt và thời gian đóng ngắt ứng với các trường hợp của khâu Relay ở câu a theo bảng sau: SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 134 ■ Hướng dẫn: Khi điều khiển ON-OFF , ngõ ra của hệ thống có dạng dao động quanh giá trị đặt, sai số của nó được đánh giá qua biên độ của sai lệch nhiệt độ: Δ e = phan hoi - Đặt khi hệ thống có dao động ổn định. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 135 SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 136 Câu 3. - Thí nghiệm: Xây dựng mô hình hệ thống điều khiển nhiệt độ PID như sau: Trong đó. - Chú ý, SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 137 giá trị K đã cho trước ở mô hình hàm truyền lò nhiệt K = 300. - 100 ( tượng trưng nhiệt độ đặt 100oC): Sau đó ta nhập các thông số cho bộ điều khiển PID như sau: SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 138 Tiếp theo nhập khâu bão hòa Saturation có giới hạn là upper limit = 1, lower limit = 0 Chạy mô phỏng ta được kết quả SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 139 c. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 140 Đối với phương pháp điều khiển PID nếu giá trị của K càng lớn thì sai số sẽ giảm nhưng ngõ ra se xuất hiện dao động và độ vọt lố. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 141 BÀI THÍ NGHIỆM 5 HỆ THỐNG ĐIỀU KHIỂN NHIỆT ĐỘ Phần 1. - do ta sử dụng c m biến cĩ thơng số là (K) Temperature Unit :đơn vị đo nhiệt độ là oC SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 142 Sp upper limit :nhiệt độ giới hạn trên 100 oC Sp lower limit :nhiệt độ giới hạn dưới 0 oC Select the control system :hệ thống điều khiển ON/OFF Standard or heating/cooling : kiểu điều khiển Standard Direct/reverse operation : chọn kiểu điều khiển reverse Alarm 1 ÷ 3: chọn 0. - n chọn stop và n hệ thống tạm dừng . - Ta có kết quả ở Bảng 1 T(min SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 143 T 0C T(min T 0C T(min T 0C Bảng 1 : Điều khiển ON – OFF QTQĐ nhiệt độ lò Từ kết quả ghi theo bảng 1 ta vẽ QTQĐ nhiệt độ lò như hình trên. - Sai số thực tế là ± 4 0 C Phần 2 .Điều khiển tuyến tính : SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 144 Mục đích : Khảo sát hệ thống điều khiển tuyến tính nhiệt độ với hiệu chỉnh PID Thí nghiệm : Mở cửa lò để hạ nhiệt độ xuống dưới 40oC . - T(min T C T(min T 0C T(min T 0C Bảng 2 : Điều khiển tuyến tính (PID) SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 145 QTQĐ nhiệt độ lò Kết quả ghi theo bảng 2 và vẽ QTQĐ nhiệt độ lò như hình trên. - SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3 ĐH CƠNG NGH SÀI GỊN * BÀI TẬP MATLAB Trang 146 SVTH : VƯƠNG VĂN HÙNG * CLASS : ĐĐT307.3