Đây là báo cáo của môn học Matlab hoàn chỉnh, được soạn thảo và thực hành theo đúng quy chuẩn các trường đại học. Các em sinh viên tải về và tham khảo các trình bày, số liệu cũng như kết quả để phục vụ tốt nhất cho môn học. Chúc các em đạt kết quả tốt trong học tập
MỤC LỤC CÁC BÀI THỰC HÀNH-THÍ NGHIỆM BÀI 1: VẼ ĐỒ THỊ TRONG MATLAB BÀI : MA TRẬN TRONG MATLAB BÀI : GIẢI PHƯƠNG TRÌNH TRONG MATLAB .9 BÀI : SIMULINK 12 T r a n g | 27 BÀI 1: VẼ ĐỒ THỊ TRONG MATLAB I Cơ sở lý thuyết Các lệnh thông dụng đồ thị không gian 2D Plot(x,y) : Vẽ đồ thị theo tọa độ x ( Tên biến ), y (Tên hàm) Plot( tên biến , tên hàm 1, tên biến 2, tên hàm 2…) Xlabel(‘…’) : Đặt tên cho trục x Ylabel(‘…’) : Đặt tên cho trục y Title(‘…’) : Đặt tên cho đồ thị Grid: Xuất lưới Hold on : Bảo vệ đồ vẽ Kiểu đường , dấu màu đồ thị : Ký hiệu Màu Ký hiệu Kiểu nét vẽ, đường b Xanh da trời - Nét liền g Xanh : Đường chấm r Đỏ - Đường gạch chấm c Xanh xám Đường gạch gạch m Đỏ tím O Đường o y Vàng X Đường x k Đen + Đường dấu + w Trắng * Đường hình * T r a n g | 27 II Nội dung thực hành : Bài : Vẽ đồ thị hàm y = x2 + sinx + khoảng x = (0, π) Câu lệnh : x=0:pi/20:pi; y=x.*x+sin(x)+5; plot(x,y) xlabel('Truc X') ylabel('Truc Y') title('Bai1a') grid on T r a n g | 27 Bài : Vẽ đồ thị biểu diễn mối quan hệ momen độ trượt động không đồng roto lồng sóc theo biểu thức : Trong : Mth=18N.m sth = 0.45 , n0= 2000v/ph, n=n0(1-s) Đồ thị có màu đỏ, đường liền nét, đánh dấu điểm chọn dấu o Câu lệnh : Mth=18; sth=0.45; n0=2000; s=0:0.01:1; n=n0.*(1-s); M=2.*Mth./((sth./s)+(s./sth)); plot(M,n,'r-o') xlabel('Momen M') ylabel('Do truot n') title('Nguyen Khanh Du') grid on T r a n g | 27 Bài : Vẽ đồ hàm số y = sin(1/x) , với x lấy khoảng phù hợp Câu lệnh : x=0:0.01:1; y=sin(1./x); plot(x,y) grid on T r a n g | 27 Bài : Vẽ đồ thị hàm y = sin(x) với đồ thị màu lam, đường liền nét đánh dấu điểm chọn dấu * , trục x thay đổi từ tới 2π , bước nhảy π/8 Câu lệnh x=0:pi/8:2*pi y=sin(x) plot(x,y,'b-*') xlabel('Truc X') ylabel('Truc Y') title('Do thi y=sin(x)') T r a n g | 27 T r a n g | 27 BÀI : MA TRẬN TRONG MATLAB I Lý thuyết Khái niệm: Trong Matlab liệu để đưa vào xử lý dạng ma trận Ma trận A có n hàng , m cột gọi ma trận cỡ n x m Được kí hiệu An x m Phần tử aij ma trận An x m phần tử nằm hàng thứ i , cột j Ma trận đơn ma trận có hàng cột Ma trận hàng ( x m) số liệu bố trí hàng a11 a12 ….a1n Ma trận cột ( n x ) bố trí cột Xử lí ma trận a Tạo Véc tơ từ ma trận b Gọi phần tử ma trận c Gọi ma trận từ ma trận lớn Các phép toán ma trận a b c d e II Phép chuyển vị Phép cộng - trừ ma trận (+,-) Phép nhân, chia ma trận (*, /) Phép quay ma trận Phép đảo ma trận Nội dung thực hành Bài a Nhập ma trận A : A=[16 13;5 10 11 8;9 12;4 15 14 1] b Kích thước ma trận A: size(A) ans = 4 c Lấy dòng ma trận A : B=A(1,:) B= 16 13 d Tạo ma trận C hai dòng cuối ma trận A : C=A(3:4,:) B= T r a n g | 27 15 12 14 e Tính tổng phần tử cột A Tongcot=sum(A(:,1:2:3:4)) Tongcot = 34 34 34 34 f Tính tổng phần tử dòng A Tongdong=sum(A(1:2:3:4,:)) Tongdong = 34 34 34 34 Bài a Nhập Ma trận A: A=[1 2; 4] A= Nhập Ma trận B: B=[14 13; 12 11] B= 14 13 T r a n g | 27 12 11 b Gọi phần tử ma trận A : *A11=A(1,1) A11 = *A12=A(1,2) A12 = *A21=A(2,2) A21 = *A22=A(2,2) A22 = c Gọi ma trận từ ma trận A: N=A(1:2,2) N= T r a n g 10 | 27 BÀI : GIẢI PHƯƠNG TRÌNH TRONG MATLAB I Lý thuyết Phương trình bậc hai thường hiểu phương trình đại số bậc hai ẩn Dạng tổng quát phương trình bậc hai ẩn a ≠ 0, số a, b c số (thực phức) gọi hệ số: a hệ số x², b hệ số x c số hay số hạng tự Công thức sau cho phép tính nghiệm phương trình bậc hai theo giá trị a, b c, tạm thời coi chúng số thực (xem thêm phần tổng quát hóa đây) với a khác Công thức sau: Các cấu trúc thường để giải phương trình : Cấu trúc if – else – and if biểu thức điều kiện khối lệnh thực điều kiện elseif biểu thức điều kiện khối lệnh thực điều kiện else khối lệnh thực điều kiện sai end Cấu trúc switch case switch biểu thức điều kiện case {giá trị thử} khối lệnh case {giá trị thử 2, giá trị thử 3, giá trị thử 4} khối lệnh otherwise T r a n g 13 | 27 khối lệnh end Lập trình dạng Function file Việc xây dựng hàm thực tương tự scripfile Tuy nhiên, hàm ta cần quan tâm đến tham số truyền vào cho hàm kết trả sau thực II Hàm có đối số vào : function tenham( giá trị vào) Hàm có đối số ra: function [out1, out2,…]= tenham Hàm có đối số vào/ra : function [x,y,…]=tenham(a,b,c,…) Nội dung thực hành Bài Giải phương trình bậc : ax2 + bx + c = Lập trình dạng Function file : function gptb2(a,b,c); if nargin u3 = tan(u1) – u2 – u1 Hệ thống Simulink : T r a n g 21 | 27 Khối scope Khối XY graph Xây dựng Simulink mơ hình động điện chiều có hệ phương trình đơn giản hóa : Me = ia.ϕ U = Ra.ia + ωϕ + Ta dia/dt T r a n g 22 | 27 Khai triển rút biến M e Ia Me=ia Φ Ia=[(-0.4.ia - + 1) Hệ thống Simulink : Khối scope : T r a n g 23 | 27 Xây dựng Simulink mơ hình động khơng đồng có hệ phương trình đơn giản : Ud = r.id – Ψq Uq = r.iq + Ψd = + ΨQ.s.ωST = - - ΨD.s.ωST Ψd = x.id + iD Ψq = x.iq + iQ ΨD = µ.x.iq + iD ΨQ = µ.x.id + iQ TM= Me - Mc Me = Ψd.iq – Ψq.id ω.r = – s Biết : Ud = 0;Uq = 1;r = 0.066;Tr = 0.253 ωST = 314; x = 2.12; µ = 0.91; TM = 1.5;Mc = 0.3 Biến đổi rút ẩn q = 0.066.id d = 1-0.066.iq D = [iD(-3.95) + (1 - r) Q.314] Q = [iQ(-3.95) - (1 - r) Q.314] id = ( d – iD).0,47 iq = ( q – iQ).0,47 iD = D – id.1,9292 iQ = Q – iq.1,9292 r = ( Me = d.id - q.id T r a n g 24 | 27 Hệ thống Simulink : T r a n g 25 | 27 Khối Scope biểu diễn thuộc tính ω.r , Me , id , iq Khối X,Y biểu diễn quan hệ Me ω.r T r a n g 26 | 27 T r a n g 27 | 27 ... Khối Mux: Tập hợp nhiều đầu vào thành đầu (nằm thư viện Signal routing ) Khối hàm Sin: Tính hàm sin đầu vào (nằm thư viện Math operations ) Nội dung thực hành Mô tả hệ thống Simulink , khảo... Y') title('Do thi y=sin(x)') T r a n g | 27 T r a n g | 27 BÀI : MA TRẬN TRONG MATLAB I Lý thuyết Khái niệm: Trong Matlab liệu để đưa vào xử lý dạng ma trận Ma trận A có n hàng , m cột gọi ma... trừ ma trận (+,-) Phép nhân, chia ma trận (*, /) Phép quay ma trận Phép đảo ma trận Nội dung thực hành Bài a Nhập ma trận A : A=[16 13;5 10 11 8;9 12;4 15 14 1] b Kích thước ma trận A: size(A)