Bài tập Matlab căn bản Nguyễn Hoài Sơn 1. Sử dụng lệnh điều kiện if viết hàm để tính giá trị hàm g(x) như sau: 2. Viết một hàm(functon) sắp xếp lại giá trị của một mảng từ nhỏ dến lớn (theo cột) và ngược lại và xuất ra giá trị trung bình của mảng input A=5 9 6 3 5 4 Output B=5 3 4 5 6 9 và C=9 6 5 4 3 5 Ví dụ : 3. Một Rocket được phóng lên từ mặt đất và được mô tả bởi hệ phương trình sau: 0.4 1.8 0 52 y x x ⎧ ⎪ = ⎨ ⎪ ⎩ ≤ ≤ a) Tạo một Script file để vẽ quỹ đạo của Rocket trên. b) Hiệu chỉnh script file trên thành function với hai thông số nhập vào là a,b và phương trình quỹ đạo của Rocket theo phương y la y a bx = 4. Ưng suất pháp và ứng suất tiếp trên mặt cắt xiên của thanh bị kéo nén theo góc alpha theo công thức như sau: Ứng suất pháp trên mặt cắt xiên z z z cos cos 2 2 σ α σ α ′ = = Ưng suất tiếp trên mặt cắt xiên sin 2 z sin 2 N F 2 2 z z NF σ τ α α ′ = = Yêu cầu : Viết 1 function tính ứng suất pháp và ứng suất tiếp trên mặt cắt xiên. Kiểm tra tham số đưa vào hàm, Đối với 3 tham số(Nz,F,Alpha) và hai tham số đưa vào(σ z ,alpha) dùng biến nargin và 1 function đổi độ sang rad Viết 1 script file vẽ đồ thị biến thiên của ứng suất pháp và ứng suất tiếp trên mặt cắt xiên như hình sau.
Bài tập Matlab Nguyễn Hoài Sơn Sử dụng lệnh điều kiện if viết hàm để tính giá trị hàm g(x) sau: Viết hàm(functon) xếp lại giá trị mảng từ nhỏ dến lớn (theo cột) ngược lại xuất giá trị trung bình mảng input A=[5 -5 4] Ví dụ : Output B=[-5 9] vaø C=[9 -5] Một Rocket phóng lên từ mặt đất mô tả hệ phương trình sau: ⎧⎪ y = 0.4 1.8 x ⎨ ⎪⎩ ≤ x ≤ 52 a) Tạo Script file để vẽ quỹ đạo Rocket b) Hiệu chỉnh script file thành function với hai thông số nhập vào a,b phương trình quỹ đạo Rocket theo phương y la ø y = a bx ng suất pháp ứng suất tiếp mặt cắt xiên bị kéo nén theo góc alpha theo công thức sau: N Ứng suất pháp mặt cắt xiên σ z′ = z cos α = σ z cos α F Nz σ ng suất tiếp mặt cắt xiên τ z′ = sin 2α = z sin 2α 2F Yêu cầu : Viết function tính ứng suất pháp ứng suất tiếp mặt cắt xiên Kiểm tra tham số đưa vào hàm, Đối với tham số(Nz,F,Alpha) hai tham số đưa vào( σ z ,alpha) dùng biến nargin function đổi độ sang rad Viết script file vẽ đồ thị biến thiên ứng suất pháp ứng suất tiếp mặt cắt xiên hình sau Giải hệ phương trình phi tuyến sau : Bài tập Matlab Nguyễn Hoài Sơn Dùng x=A\b Một người chạy từ dốc có độ nghiên alpha=20o để tránh đá Nếu người đến vị trí x=0 trước viên đá thoát Hãy viết chương trình mô chuyển động người đá HD: Phương trình chuyển động viên đá Trong : m khối lượng viên đá v& Vận tốc đá B hệ số ma sát G gia tốc trọng trường Vận tốc vị trí đá , Hai pt cho người dt : bước thời gian v(i) : vận tốc thời điểm t=to+t*i v(i-1) : vận tốc thời điểm t=to+t*(i-1) Giá trị: B1=10, m1=1000, T=10 (thời gian mô phỏng) B2=0.2, m2=100, gia tốc =1.5 , v_max=5 m/s Xây dựng ma trận 16x16 sau: Bài tập Matlab Với I= Nguyễn Hoài Sơn Viết function I=trap(f,a,b,n) dùng tính tích phân Trong : h=(a-b)/n xi = a+i*h Kiểm tra kết hàm với hàm f=cos(x)+sin(x) cho a=0,b=pi/3 Kiểm tra kết quy luật Simpson I= 10 Dựa vào lưu đồ sau xây dựng script file dùng lệnh input, sprinft 11 Sử dụng vòng lặp for lệnh để xây dựng script file dựa lưu đồ sau: Bài tập Matlab 12 Nguyễn Hoài Sơn Sử dụng vòng lặp while xây dựng script file dựa lưu đồ sau: Bài tập Matlab Nguyễn Hoài Sơn 13 14 15 Chuyển đổi script file(8,9,10) thành hàm (function) Vẽ hình tròn Vẽ mặt (surface ) sau: 16 Vẽ mặt dựa phương trình tham số sau : Với u,v=0 :2*pi HD : Tính u,v sau sử dụng lệnh surf(X,Y,Z) 14 Vẽ mặt contour sau : Bài tập Matlab Nguyễn Hoài Sơn Với 17 sử dụng hàm contour Dưới lệu thu từ thí nghiệm mà mô hình toán học có công thức : Trong xm,ym giá trị đo δ y sai số giá trị xm Sử dụng hàm errorbar để thể kết 16 Xây dựng structure sau : Bài tập Matlab Nguyễn Hoài Sơn BÀI T P MATLAB C N B N Vi t ch ng trình tính ti n n cho m t gia đình bi t r ng n u s Kwh n s d ng: - < 100 giá Kwh 500đ - 100 ≤ Kwh < 200 giá Kwh 700đ - 200 ≤ Kwh < 500 giá Kwh 1000đ - Kwh > 500 giá Kwh 1200đ ng th i: + N u Kwh < in ra: khơng h p l + N u Kwh s d ng ≥ 500 sau dịng in t ng s ti n có ghi kèm theo câu: H n ch s d ng n Vi t ch ng trình gi i bi n lu n b t ph ng trình : ax > 3x + Vi t ch ng trình tính doanh thu trung bình DTTB n m c a Công ty X theo doanh thu c a quý: Q1, Q2, Q3 Q4, đ n v tính Tri u đ ng N u DTTB < 100 in ra: Phá s n, t 100 đ n nh h n 200 in ra: Doanh thu T 200 đ n nh h n 300 in ra: Doanh thu trung bình, t 300 đ n nh h n 500 in ra: Doanh thu khá, t 500 tr lên in ra: Doanh thu t t Vi t ch ng trình tính n! Vi t ch ng trình tính c n b c c a s : - L y xn = x/2 - 1⎛ x⎞ S d ng công th c: xn+1 = ⎜⎜ xn + ⎟⎟ 2⎝ xn ⎠ - Ki m tra u ki n: Vi t ch xn +1 − xn 0 OUTPUT N INPUT R Area = π * R INPUT R Circ = * π * R If True Area > 20.0 OUTPUT R, Area, Circ False N=N+1 Stop Bài tập Matlab Nguyễn Hoài Sơn ng trình tính di n tích hình thang theo đáy l n a, đáy nh b chi u cao Vi t ch h, có xét đ n s t n t i ho c không t n t i c a hình thang ( s d ng tốn t &) 10 Vi t ch ng trình th c hi n s l a ch n: tính di n tích hình vng theo c nh a ho c tính di n tích hình trịn theo đ ng kính d, có xét đ n s t n t i hay không t n t i c a hình vng hình trịn 11 Vi t ch ng trình h ng d n mua vé máy bay ( th c hi n vi c l a ch n): mu n t TPHCM Hà N i: đ n qu y s 1, TPHCM N ng: đ n qu y s 2, TPHCM Nha Trang: đ n qu y s 3, TPHCM V ng Tàu: đ n qu y s 4, TPHCM L t: đ n qu y s 5, TPHCM Côn 12 Vi t ph ng trình gi i bi n lu n ph o: đ n qu y s ng trình: ⎧a1 x + b1 y + c1 = ⎨ ⎩a2 x + b2 y + c2 = Theo tham s a1 , b1 , c1 , a2 , b2 , c2 CÁC PHÉP TÍNH TRÊN MA TR N VÀ M NG 13 Cho: 5⎤ ⎡2 ⎢ A = ⎢8 12 20 ⎥⎥ ⎢⎣4 12 ⎥⎦ a Tính B = A + 12 b Tính C = 4*A c Gõi D ma tr n chuy n trí c a C Xác đ nh D 14 Cho: ⎡10 ⎤ ⎡2 ⎤ ⎢ ⎥ A = ⎢8 ⎥ B = ⎢⎢8 12 ⎥⎥ ⎣⎢4 14 20 ⎥⎦ ⎣⎢6 10 ⎦⎥ a Tính ma tr n t ng C c a A B b Tính ma tr n tích D = A*B c Xác đ nh ma tr n ngh ch đ o E c a D 15 Cho ma tr n A véct hàng B: Bài tập Matlab Nguyễn Hoài Sơn ⎡1 ⎤ A = ⎢⎢4 6⎥⎥ B = [2 ⎢⎣7 ⎥⎦ Vi t ch 8] ng trình th c hi n liên ti p công vi c sau: a Nh p A B b T o ma tr n C c a ma tr n A b ng cách l y s h ng th 1, 2, c Tính ma tr n t ng D c a B C d Tính ma tr n E ma tr n chuy n trí c a B e Tính F = D*E Vi t ch ng trình cho k t qu đ c hi n th hình 16 Cho ⎡2 12 15 ⎤ A = ⎢⎢24 27 ⎥⎥ ⎢⎣18 20 ⎥⎦ ⎡4 ⎤ B = ⎢⎢12 ⎥⎥ ⎢⎣0 10 ⎥⎦ Hãy th c hi n phép chia A./B cho bi t INF nói lên u ? 17 Vi t ch ng trình tính giá tr y = sin(x), giá tr c a x ≤ x ≤ π b c nh y π / 18 Hãy s d ng ones, eye ph ng pháp xây d ng ma tr n t kh i đ t o nên ma tr n sau: ⎡2 ⎢2 ⎢ ⎢2 A= ⎢ ⎢4 ⎢0 ⎢ ⎣⎢0 2 0⎤ 2 ⎥⎥ 2 0 4⎥ ⎥ 0 2 2⎥ 2 2⎥ ⎥ 2 2⎦⎥ 19 Hãy s d ng ones, eye ph sau: ⎡1 ⎢1 C =⎢ ⎢0 ⎢ ⎣0 0⎤ 0⎥⎥ 1⎥ ⎥ 1⎦ ng pháp xây d ng t kh i đ t o nên ma tr n Bài tập Matlab Nguyễn Hoài Sôn Bi n đ i C thành: ⎡1 ⎢1 D=⎢ ⎢0 ⎢ ⎣0 20 Gi i ph 0⎤ ⎡1 0⎤ 0⎥⎥ E = ⎢⎢1 0⎥⎥ 1⎥ ⎢⎣0 1 ⎥⎦ ⎥ 1⎦ ng trình sau b ng ph ng pháp đ nh th c: ⎧15 x +13 y = 440 ⎨ ⎩24 x + y = 428 21 Gi i h ph ng trình sau b ng ph ng pháp ngh ch đ o ma tr n: ⎧4 x1 + 3x2 + x3 + x4 = 884 ⎪5 x + x + x + 10 x = 1248 ⎪ ⎨ = 544 ⎪12 x1 + x2 + x3 ⎪8 x + x + 12 x + 20 x = 1908 ⎩ 22 Gi i h ph ng trình sau b ng ph ng pháp chia trái ma tr n: ⎧5 x1 + x2 + x3 = 82 ⎪ ⎨4 x1 + x2 + 10 x3 = 96 ⎪8 x + 3x + x = 52 ⎩ 23 Liệt kê phần tử với điều kiện x= -0.4326 -1.6656 0.1253 0.2877 -1.1465 1.1909 1.1892 ind = (x >= 1) | (x < -0.2) ind = 1100111 y = x(ind) y= -0.4326 -1.6656 -1.1465 1.1909 1.1892 a = [1 1] a= 11341 ind = (a == 1) ind = 11001 b = a(ind) b= 111 24 Tìm phần tử ma trận vectơ cho biết vị trí chúng %tim dia chi voi a = ind = find(a == 1) ind = Bài tập Matlab 125 25 xếp từ nhỏ đến lớn function [b, j] = descsort(a) [b ,j] = sort(-a); b = -b; % main program a = [pi –10 35 0.15]; [b, j] = descsort(a) ket qua: b= 35.0000 3.1416 0.1500 -10.0000 j= 26 tạo hàm tính giá trị hàm “inline” feval f = inline('sqrt(x.^2+y.^2)','x','y') f= Inline function: f(x,y) = sqrt(x.^2+y.^2) % tinh f: f(3,4) ans = A = [1 2;3 4] A= 12 34 B = ones(2) B= 11 11 C = f(A, B) C= 1.4142 2.2361 3.1623 4.1231 Nguyễn Hoài Sơn Bài tập Matlab Nguyễn Hoài Sơn BÀI T P MATLAB C N B N PH N I N Baøi 1: cho mạch điện hình vẽ, tìm điện nút V1, V2, V3 Tại nút 1: Tại nút 2: Tại nút 3: Ta có hệ: Bài 2: cho mạch điện hình vẽ, tìm điện nút V1, V2, V3, V4 Bài tập Matlab Tại nút 1: Tại nút 2: Nút 2, 3: Tại nút 4: Cuối ta có hệ: Bài 3: Xác định giá trị dòng điện chạy qua điện trở Ta có: Nguyễn Hoài Sơn Bài tập Matlab Vòng 1: Vòng 2: Vòng 3: Ta có hệ phương trình Bài 4: Tìm điện nút mạch điện Phân tích mạch điện ta có hệ phương trình: Nguyễn Hoài Sơn Bài tập Matlab Nguyễn Hoài Sơn Thế số ta được: Bài 5: cho mạch điện hình vẽ: Từ mạch điện ta có phương trình: Cho C = 10 μ F sử dụng Matlab để vẽ điện qua tụ điện với giá trị điện trở R = k Ω, R = 10 k Ω, R = 0.1 k Ω Bài 6: Cho mạch điện RLC hình vẽ Neáu di ( ) L = 10 H , R = 400 Ω, C = 100 μ F , vs ( t ) = 0, i ( ) = A, = 15 A / s dt Tìm i ( t ) Cho Bài 7: Cho mạch điện hình vẽ Nếu thời điểm ban đầu t = công tắc mở Tìm điện v ( t ) Với L = 1/ 32 H , R = 10 Ω, C = 50 μ F , I s = A Bài tập Matlab Cho Nguyễn Hoài Sơn ... Nguyễn Hoài Sơn Bài tập Matlab Nguyễn Hoài Sơn BÀI T P MATLAB C N B N PH N I N Bài 1: cho mạch điện hình vẽ, tìm điện nút V1, V2, V3 Tại nút 1: Tại nút 2: Tại nút 3: Ta có hệ: Bài 2: cho mạch... nút V1, V2, V3, V4 Bài tập Matlab Tại nút 1: Tại nút 2: Nút 2, 3: Tại nút 4: Cuối ta có hệ: Bài 3: Xác định giá trị dòng điện chạy qua điện trở Ta có: Nguyễn Hoài Sơn Bài tập Matlab Vòng 1: Vòng... lệnh để xây dựng script file dựa lưu đồ sau: Bài tập Matlab 12 Nguyễn Hoài Sơn Sử dụng vòng lặp while xây dựng script file dựa lưu đồ sau: Bài tập Matlab Nguyễn Hoài Sơn 13 14 15 Chuyển đổi script