báo cáo chuyên đề Ứng dụng MATLAB/Simulink trong điều khiển

TRƯỜNG ĐẠI HỌC ĐIỆN LỰC KHOA ĐIỀU KHIỂN & TỰ ĐỘNG HĨA BÁO CÁO CHUN ĐỀ NGÀNH: CƠNG NGHỆ KTĐK&TĐH CHUYÊN NGÀNH: TDH&DKTBDCN HỌC PHẦN: Ứng dụng MATLAB/Simulink điều khiển Giảng viên hướng dẫn: Nguyễn Ngọc Khốt Nhóm sinh viên/ sinh viên thực hiện: Nhóm 1: Võ Nam An – MSV: 20810410072 Hồ Tuấn Anh – MSV: 20810410076 Nguyễn Huy Việt Anh – MSV: 20810410060 Nguyễn Tùng Anh – MSV: 20810410069 STT: N = Chu Tam Thức – MSV: 20810430312 Lớp: D15TDH&DKTBDCN1 HÀ NỘI, 9/2022 MỤC LỤC {Để dòng trống đây} {Font: Time New Roman; thường; cỡ chữ: 13; dãn dòng: 1,2; lề: justified} Trang Chương 1: Ứng dụng MATLAB giải tích Bài 1: a) b) c) Bài 2: a) b) c) 10 12 12 Chương : Khảo sát hệ thống điều khiển dùng MATLAB Bài : a) 14 b) 21 28 Bài : KẾT LUẬN TÀI LIỆU THAM KHẢO ĐỀ BÀI BÁO CÁO CHUYÊN ĐỀ Chương 1: Ứng dụng MATLAB giải tích Bài 1: Giải phương trình, hệ phương trình nhiều cách Bài 2: (Nhóm 1->4) Cho hàm số: có đồ thị (C) hệ tọa độ thực Oxy Hãy viết chương trình m-file (đặt tên TT_TEN_Exam_Cau_1.m) sử dụng phần mềm MATLAB để: (VD: TT_24_Khoi_exam_cau_1.m) a) Vẽ đồ thị (C) miền Dx = [-10 - N/10; 10 + N/10] b) Tìm giao điểm đồ thị (C) vừa vẽ với trục tọa độ phương pháp đại số giải phương trình phương pháp đồ thị c) Vẽ đồ thị hàm số y = g(x) = -5x + N/2 (d) hệ trục tọa độ tìm giao điểm hai đồ thị (d) (C) (nếu có) Chương 2: Khảo sát hệ thống điều khiển dùng MATLAB 2.1 Nhiệm vụ riêng cho nhóm Bài 3: Nhóm 2, 4, 6, 8: Khảo sát đặc tính động học hệ thống cho hàm truyền đạt dạng zpk sau: Yêu cầu: - Viết m-file nhập hàm truyền trên; - Xác định điểm cực, điểm khơng Xác định tính ổn định hệ thống - Vẽ đồ thị Bode, Nyquist - Vẽ đặc tính độ (đáp ứng bước nhảy), đặc tính xung Dirac (hàm trọng lượng).Xác định độ điều chỉnh, thời gian xác lập, sai số xác lập, thời gian tăng - Xác định tính điều khiển được, quan sát hệ thống - Vẽ mơ hình Simulink tương ứng với hàm truyền hệ hở phản hồi âm đơn vị Vẽ đặc tính đầu y(t) với đầu vào hàm 1(t), ramp, sin, pulse 2.2 Nhiệm vụ chung cho tất nhóm Bài 4: Hãy viết m-file để thực sơ đồ điều khiển hệ LTI sau Vẽ đáp ứng hàm độ, hàm trọng lượng, đồ thị Bode, Nyquist xác định tính ổn định hệ thống Chương 1: Ứng dụng MATLAB giải tích Bài 1: Giải phương trình, hệ phương trình nhiều cách ( N = ) a) Cách 1: ● Editor : syms x eqn = 5*x^3-5*x+2*cos(x)==0 S = solve(eqn,x) ● Command Window : >> Bai1 S= 0.45231269708657210661962396478477 Cách 2: ● Editor : syms x S = solve('5*x^3-5*x+2*cos(x)==0','x') ● Command Window : >> Bai1 S= 0.45231269708657210661962396478477 b) Cách 1: ● Editor : syms x y pt1 = x^2 - 4*y^2 + 9*x ==2 pt2 = 8*x^2 + y^2 - 9*y == 10 [x,y] = solve (pt1,pt2,x,y) ● Command Window : >> Bai1 pt1 = x^2 + 9*x - 4*y^2 == pt2 = 8*x^2 + y^2 - 9*y == 10 x= (11*root(z^4 - (6*z^3)/11 + (125*z^2)/121 - (60*z)/11 - 716/121, z, 1)^2)/24 - root(z^4 - (6*z^3)/11 + (125*z^2)/121 - (60*z)/11 - 716/121, z, 1)/8 + 1/12 (11*root(z^4 - (6*z^3)/11 + (125*z^2)/121 - (60*z)/11 - 716/121, z, 2)^2)/24 - root(z^4 - (6*z^3)/11 + (125*z^2)/121 - (60*z)/11 - 716/121, z, 2)/8 + 1/12 (11*root(z^4 - (6*z^3)/11 + (125*z^2)/121 - (60*z)/11 - 716/121, z, 3)^2)/24 - root(z^4 - (6*z^3)/11 + (125*z^2)/121 - (60*z)/11 - 716/121, z, 3)/8 + 1/12 (11*root(z^4 - (6*z^3)/11 + (125*z^2)/121 - (60*z)/11 - 716/121, z, 4)^2)/24 - root(z^4 - (6*z^3)/11 + (125*z^2)/121 - (60*z)/11 - 716/121, z, 4)/8 + 1/12 y= root(z^4 - (6*z^3)/11 + (125*z^2)/121 - (60*z)/11 - 716/121, z, 1) root(z^4 - (6*z^3)/11 + (125*z^2)/121 - (60*z)/11 - 716/121, z, 2) root(z^4 - (6*z^3)/11 + (125*z^2)/121 - (60*z)/11 - 716/121, z, 3) root(z^4 - (6*z^3)/11 + (125*z^2)/121 - (60*z)/11 - 716/121, z, 4) >> Cách 2: ● Editor : syms x y a = input('Nhap a= ') b = input('Nhap b= ') c = input('Nhap c= ') d = input('Nhap d= ') e = input('Nhap e= ') f = input('Nhap f= ') g = input('Nhap g= ') h = input('Nhap h= ') pt1=a*x^2-b*y^2+c*x==d pt2=e*x^2+f*y^2-g*y==h [x,y] = solve(pt1,pt2,x,y) ● Command Window : Nhap a= a= Nhap b= -2 b= -2 Nhap c= c= Nhap d= d= Nhap e= e= Nhap f= f= Nhap g= -9 g= -9 Nhap h= 10 h= 10 pt1 = x^2 + 9*x + 2*y^2 == pt2 = 8*x^2 + y^2 + 9*y == 10 x= root(z^3 - z^2/5 + (56*z)/25 + 716/25, z, 1)/8 - (5*root(z^3 - z^2/5 + (56*z)/25 + 716/25, z, 1)^2)/24 + 1/12 root(z^3 - z^2/5 + (56*z)/25 + 716/25, z, 2)/8 - (5*root(z^3 - z^2/5 + (56*z)/25 + 716/25, z, 2)^2)/24 + 1/12 root(z^3 - z^2/5 + (56*z)/25 + 716/25, z, 3)/8 - (5*root(z^3 - z^2/5 + (56*z)/25 + 716/25, z, 3)^2)/24 + 1/12 y= root(z^3 - z^2/5 + (56*z)/25 + 716/25, z, 1) root(z^3 - z^2/5 + (56*z)/25 + 716/25, z, 2) root(z^3 - z^2/5 + (56*z)/25 + 716/25, z, 3) >> >> Bai1 c) Cách 1: Ma trận ● Editor : A = [3 -2 12;4 -6 16;3 -1 -4] B = [13;17;0] X=A\B ● Command Window : >> Bai1 A= ... 1: Ứng dụng MATLAB giải tích Bài 1: a) b) c) Bài 2: a) b) c) 10 12 12 Chương : Khảo sát hệ thống điều khiển dùng MATLAB Bài : a) 14 b) 21 28 Bài : KẾT LUẬN TÀI LIỆU THAM KHẢO ĐỀ BÀI BÁO CÁO... dùng MATLAB Bài : a) 14 b) 21 28 Bài : KẾT LUẬN TÀI LIỆU THAM KHẢO ĐỀ BÀI BÁO CÁO CHUYÊN ĐỀ Chương 1: Ứng dụng MATLAB giải tích Bài 1: Giải phương trình, hệ phương trình nhiều cách Bài 2: (Nhóm... - Vẽ đặc tính độ (đáp ứng bước nhảy), đặc tính xung Dirac (hàm trọng lượng).Xác định độ điều chỉnh, thời gian xác lập, sai số xác lập, thời gian tăng - Xác định tính điều khiển được, quan sát