THỰC HÀNH lập TRÌNH TRÊN MATLAB

PHẦN I: THỰC HÀNH LẬP TRÌNH TRÊN MATLAB Bai 1_VẼ ĐỒ THỊ TRONG MATLAB Lí thuyết: - Đồ thị không gian 2D: • Lệnh “plot” plot(x,y) với y=f(x) – hàm số cần vẽ x – vectoer miền giá trị hàm y • Các tham số plot Plot(x,y,’linestye_marker_color’) Color b G R C M Y K w Linestyle : - -O X + * Marker S D V ^ < > P H Xanh da trời Xanh Đỏ Xanh xám Đỏ tím Vàng Đen Trắng Nét liền Đường chấm Đường gạch chấm Đường gạch gạch Đường o Đường x Đường dấu + Đường dấu * Vuông Diamond Triangle(down) Triangle(up) Triangle(left) Triangle(right) Pentagram Hexagram • • • • • - Lệnh “hold on” Khi muốn vẽ thêm đô thị đồ thị hiên có ta sử dụng lệnh “hold on” Để tắt chế độ ta dùng lệnh “hold off” Lệnh “ grid on “ Một lệnh dùng để tùy chỉnh kểu trục tọa độ,giúp nhìn đồ thị trực quan Lệnh “axis” axis=([xmin xmax ymin ymax]) Dùng để giới hạn đồ thị theo mong muốn Chú thích đồ thị Các lệnh “ xlabel,ylabel ” để đặt tên cho trục đồ thị xlabel(‘truc x’) ylabel(‘truc y’) Lệnh “ title(‘ ’) “ để đặt tên cho đồ thị Lệnh “ linspace “ để xác định giới hạn trục Ox,Oy x=linspace(0,pi,30) Đồ thị không gian 3D: • Lệnh “plot” plot(x,y,z) • Các lệnh khác tương tự đồ thị 2D Bài làm: Bài 1-1: Vẽ đồ thị hàm sau: y = x2 + sinx + khoảng (0,π) Bài Làm Nhập : x=0:0.1:pi; y=x.^2+sin(x)+5; plot(x,y) xlabel('Truc x') ylabel('Truc y') grid on title('Tran Thanh Duong') Nhấn Enter Bài 1-2 Cho biết mối quan hệ mô men độ trượt động không đồng rô to lồng sóc tuân theo biểu thức: M = M th S th S + S S th - Lập trình Matlab vẽ đồ thị đặc tính động với đồ thị màu đỏ, đường liền nét đánh dấu điểm chọn dấu o, biết Mth=18N.m, Sth = 0.45, n0 = 2000v/ph Đặt tên đồ thị theo họ tên Bài Làm Nhập : Mth=18; Sth=0.45; n0=2000; s=0:0.1:1; M=2*Mth./(Sth./s+s./Sth); n=n0*(1-s); plot(M,n,'r-o') xlabel('Momen') ylabel('Tocdo') grid on title('Tran Thanh Duong') Nhấn Enter Bài 1-3: Vẽ đồ thị hàm số sin(1/x), với x lấy khoảng phù hợp Bài Làm Nhập : x=0:0.2:2*pi; y=sin(1./x); plot(x,y) xlabel('Truc x') ylabel('Truc y') grid on title('Tran Thanh Duong') Nhán Enter Bài 1-4: 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 thay đổi π/8 Bài Làm Nhập : x=0:pi./8:2*pi; y=sin(x); grid on plot(x,y,'c-*') title('Tran Thanh Duong') Nhán Enter Bai 2_MA TRAN TRONG MATLAB Lí thuyết: - - - Khái niệm: 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ỡ m x n ,kí hiệu Am x n Phần tử aij phần tử nằm hàng thứ i ,cột j Cách để nhập ma trận • Liệt kê trực tiếp: >> A=[1 3;4 6;7 9] • Thông qua lệnh nhập: >>input(‘ nhap gia tri cho ma tran A’) nhap gia tri cho ma tran A=[1 3;4 6;7 9] Tìm kích thước ma trận >>size(_ten_matran_) Trích xuất phần tử ma trận Vidu: A=[16 13;5 10 11 8;9 12;4 15 14 1]; B=A(1,:) >> B = 16 13 - Các phép toán với ma trận • Tính tổng phần tử hàng, cột A=[16 13;5 %Tinh tong cac B=[sum(A(:,1)) %Tinh tong cac B=[sum(A(1,:)) >> B = 34 >> 34 34 34 34 34 34 B = 34 • 10 11 8;9 12;4 15 14 1]; phan tu tren cac cot cua A sum(A(:,2)) sum(A(:,3)) sum(A(:,4))] phan tu tren cac dong cua A sum(A(2,:)) sum(A(3,:)) sum(A(4,:))] Tìm ma trận chuyển vị A=[5 6;7 1;3 4] %Tim ma tran chuyen vi cua ma tran tren M=A' >> M = • 9 Phép cộng, trừ, nhân, chia A=[5 6;7 1;3 4] B=[15 62 13;23 31 34;41 90 36] %Phép c?ng tr? nhân chia% C=A+B D=A-B E=A*B F=A/B >> C = 20 30 44 >> 19 35 40 -59 -22 -81 -7 -33 -32 D = -10 -16 -38 >> 65 40 99 E = 390 353 416 >> 943 803 825 383 433 489 F = -0.1978 -0.4728 0.1506 0.1142 -0.3786 0.1008 0.1303 0.5561 -0.0385 Phép quay ma trận • %Thuc hien phep Quay va Dao Z=rot90(A,2) >> Z = 9 3 Bài làm: Bài 2-1: Nhập vào ma trận: A=[16 13; 10 11 8; 12; 15 14 1] - Tìm kích thước ma trận A Lấy dòng ma trận A Tạo ma trận B dòng cuối A Tính tổng phần tử cột A (gợi ý: tính tổng phần tử cột 1: sum(A(:,1))) - Tính tổng phần tử dòng A Bài Làm Nhập : A=[16 13;5 10 11 8;9 12;4 15 14 1]; %Tim kich thuoc cua ma tran A size(A) %Lay dong dau tien cua ma tran A B=A(1,:) %Tao ma tran bang dong cuoi cung cua A B=A(2:3,:) %Tinh tong cac phan tu tren cac cot cua A B=[sum(A(:,1)) sum(A(:,2)) sum(A(:,3)) sum(A(:,4))] %Tinh tong cac phan tu tren cac dong cua A B=[sum(A(1,:)) sum(A(2,:)) sum(A(3,:)) sum(A(4,:))] Nhấn Enter ans = 4 16 13 10 11 12 34 34 34 34 34 34 34 34 B = B = B = B = Bài 2-2: Nhập vào ma trận A, B bất kỳ, thực yêu cầu sau: - Gọi phần tử ma trận A Gọi ma trận từ ma trận A B Tìm ma trận chuyển vị ma trận Thực phép nhân, chia, cộng, trừ ma trận Thực phép quay đảo ma trận Bài Làm Nhập : A=[5 6;7 1;3 4] B=[15 62 13;23 31 34;41 90 36] %Goi cac phan tu ma tran A C=A(1:2) %Goi mot ma tran tu ma tran A hoac B D=B(1:3,2:2) %Tim ma tran chuyen vi cua ma tran tren M=A' N=B' %Thuc hien cac phep tinh Cong,Tru,Nhan,Chia cac ma tran tren Y=A+B U=A-B I=A*B O=A./B %Thuc hien phep Quay va Dao Z=rot90(A,2) X=fliplr(A) Nhán Enter C = D = 62 31 90 M = 9 15 62 13 23 31 34 41 90 36 20 30 44 65 40 99 19 35 40 -10 -16 -38 -59 -22 -81 -7 -33 -32 943 803 825 383 433 489 N = Y = U = I = 390 353 416 O = 0.3333 0.3043 0.0484 0.2903 0.4615 0.0294 10 0.0732 0.1000 0.1111 Z = 9 3 9 X = Bai 3_GIẢI PHƯƠNG TRÌNH TRONG MATLAB Lí thuyết: - Kiến thức biện luận & giải phương trình bậc hai Cho phương trình bậc có dạng ax2 + bx + c = (+) Biện luận : • TH1: a =  phương trình trở dạng phuong trình bậc : bx + c = Suy PT nghiệm x = -c/b • TH2 : a ≠ phương trình PT bậc Ta tính delta: ∆ = b2 -4*a*c Suy PT có khả năng: KN1: PT vô nghiệm ∆ < KN2: PT có nghiệm kép∆ = KN3: PT có nghiệm phân biệt∆ > x1= (- b – sqrt(∆))/2*a x2= (- b + sqrt(∆))/2*a - Định luật Ohm cho toàn mạch : 11 - I=U/R Định luật Kirhoff 1,2: Bài làm: Bài 3-1: Biện luận số nghiệm phương trình a – bx + c=0 Bài Làm Nhập : function [x1,x2]=FileThuNghiem(a,b,c) if nargin Constant => kích kéo sang tab untitle.Làm tương tự với khối lại • B4 : Nối khối theo sơ đồ cấu trúc Nối đầu mũi tên khối vào đầu mũi tên vào khối kia.Khi cần nối nhiều nối chập nối từ đầu mũi tên vào khối vào đường nối cần thiết thả Ngoài thay đổi kích thước khối,di chuyển khối, copy khối • • • B5 : Mở khối kích đúp vào khối điều chỉnh số liệu cần thiết B6 ; Mô Kích vào Simulation => Start kích trực tiếp vào biểu tượng Start/Pause Simulation B7 : Thay đổi thông số trình mô cần Bài làm: Bài 1: Cho hệ thống mô tả bởi: u1 = (u sin u1 + 5)u + u1 = (u sin u1 + 5)u +    du  = u1 + − u  u = (u1 + − u ) s  dt  u + u + u1 = tan u1 ⇒ u = −u − u1 + tan u1 Mô phỏng: 17 18 Bài 2: Xây dựng Simulink mô hình động điện chiều có hệ phương trình đơn giản hóa: 1  dω  ω = ( M e − M c )  TM s TM dt = M e − M c   M e = ia φ M e = ia φ   1 di ia = (U − ia Ra − ω.φ ) U = Ra ia + ω.φ + Ta a Ta s dt ⇒   Mô phỏng: 19 Bài 3: Xây dựng Simulink mô hình động không đồng có hệ phương trình đơn giản hóa viết hệ trục (d,q) sau: 20 U d = r.id − ψ q U q = r.iq + ψ d dψ D i i = − D + ψ Q s.ϖ ST = − D + (1 − ω r ).314.ψ Q dt Tr Tr dψ Q dt =− iQ Tr − ψ D s.ϖ ST = − iQ Tr − (1 − ω r ).314.ψ D ψ d = x.id + i D ψ q = x.i q + iQ ψ D = µ x.i d + i D ψ Q = µ.x.i q + iQ dω r = Me − Mc dt M e = ψ d iq − ψ q id TM ⇒ ψ q = −U d + r.id ψ d = U q − r.i q  i 1 ψ D = − D + (1 − ω r ).314.ψ Q   Tr  s  iQ 1 ψ Q = − − (1 − ω r ).314.ψ D   Tr  s i d = (ψ d − i D ) x i q = (ψ q − iQ ) x i D = ψ D − µ x.i d iQ = ψ Q − µ x.iq ω r = ( M e − M c ) 1 TM s M e = ψ d iq − ψ q i d Mô phỏng: 21 22 23 24 [...]... Z = 4 1 6 9 9 3 3 7 5 6 1 4 3 9 9 5 7 3 X = Bai 3_GIẢI PHƯƠNG TRÌNH TRONG MATLAB Lí thuyết: - Kiến thức về biện luận & giải phương trình bậc hai Cho phương trình bậc 2 có dạng ax2 + bx + c = 0 (+) Biện luận : • TH1: a = 0  phương trình trở về dạng phuong trình bậc nhất : bx + c = 0 Suy ra PT nghiệm duy nhất x = -c/b • TH2 : a ≠ 0 phương trình là PT bậc 2 Ta tính delta: ∆ = b2 -4*a*c Suy ra PT có 2... 6.2500 khong hop le PHẦN II: THỰC HÀNH TRÊN SIMULINK Li thuyết: - - Các khối sử dụng trong bài thực hành • Khối nguồn không đổi: Constant • Khối hiển thị : Scope,XYgraph • Các khối trong thư viện hàm liên tục : Integrator,Derivative • Các khối trong thư viện Math Operator : Gain,Product,Sum,Trigonometric Function • Khối chính trong thư viện Signal routing : Mux Các bước thực hiện Mô phỏng hệ thống bằng... I=input('nhap I:'); Ur=input('nhap Ur:'); f=input('nhap f:'); R=Ur/I Uc=sqrt(U^2-Ur^2); Zc=Uc/I; C=1/(2*pi*f*Zc) Chạy trên common window nhap U:220 nhap I:2 nhap Ur:110 13 nhap f:50 R = 55 C = 3.3414e-005 Bài Về Nhà: Bài 1:Giải hệ phương trình Bằng cách chuyển về dạng ma trận và thực hiện nó trên Matlab Bài Làm Nhập : A=[4 -2 -10;2 10 -12;-4 -6 16] B=[-10;32;-16] inv(A)*B Nhấn Enter A = 4 2 -4 -2 10 -6 -10... Kirhoff 1,2: Bài làm: Bài 3-1: Biện luận số nghiệm của phương trình a – bx + c=0 Bài Làm Nhập : function [x1,x2]=FileThuNghiem(a,b,c) if nargin