Bài tập Matlab căn bản

Chia sẻ: Trần Công Khánh | Ngày: | Loại File: PDF | Số trang:0

Thêm vào BST

Báo xấu

615
lượt xem 183
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

Tài liệu về bài tập Matlab căn bản...

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: Bài tập Matlab căn bản

Baøi taäp Matlab caên baûn Nguyeãn Hoaøi Sôn 1. Söû duïng leänh ñieàu kieän if vieát haøm ñeå tính giaù trò haøm g(x) nhö sau: 2. Vieát moät haøm(functon) saép xeáp laïi giaù trò cuûa moät maûng töø nhoû deán lôùn (theo coät) vaø ngöôïc laïi vaø xuaát ra giaù trò trung bình cuûa maûng Ví duï : input A=[5 9 6 3 -5 4] Output B=[-5 3 4 5 6 9] vaø C=[9 6 5 4 3 -5] 3. Moät Rocket ñöôïc phoùng leân töø maët ñaát vaø ñöôïc moâ taû bôûi heä phöông trình sau: ⎧ ⎪ y = 0.4 1.8 x ⎨ ⎪ 0 ≤ x ≤ 52 ⎩ a) Taïo moät Script file ñeå veõ quyõ ñaïo cuûa Rocket treân. b) Hieäu chænh script file treân thaønh function vôùi hai thoâng soá nhaäp vaøo laø a,b vaø phöông trình quyõ ñaïo cuûa Rocket theo phöông y la ø y = a bx 4. Öùng suaát phaùp vaø öùng suaát tieáp treân maët caét xieân cuûa thanh bò keùo neùn theo goùc alpha theo coâng thöùc nhö sau: N ÖÙng suaát phaùp treân maët caét xieân σ z′ = z cos 2 α = σ z cos 2 α F σ Nz Öùng suaát tieáp treân maët caét xieân τ z′ = sin 2α = z sin 2α 2F 2 Yeâu caàu : Vieát 1 function tính öùng suaát phaùp vaø öùng suaát tieáp treân maët caét xieân. Kieåm tra tham soá ñöa vaøo haøm, Ñoái vôùi 3 tham soá(Nz,F,Alpha) vaø hai tham soá ñöa vaøo( σ z ,alpha) duøng bieán nargin vaø 1 function ñoåi ñoä sang rad Vieát 1 script file veõ ñoà thò bieán thieân cuûa öùng suaát phaùp vaø öùng suaát tieáp treân maët caét xieân nhö hình sau. 5. Giaûi heä phöông trình phi tuyeán sau :
Baøi taäp Matlab caên baûn Nguyeãn Hoaøi Sôn Duøng x=A\b Moät ngöôøi ñang chaïy töø treân moät doác coù ñoä nghieân alpha=20o ñeå traùnh moät 6. hoøn ñaù Neáu ngöôøi naøy ñeán vò trí x=0 tröôùc vieân ñaù thì thoaùt ñöôïc. Haõy vieát chöông trình moâ phoûng chuyeån ñoäng cuûa ngöôøi vaø hoøn ñaù. HD: Phöông trình chuyeån ñoäng cuûa vieân ñaù Trong ñoù : m khoái löôïng vieân ñaù v Vaän toác hoøn ñaù & B heä soá ma saùt G gia toác troïng tröôøng Vaän toác vaø vò trí hoøn ñaù , Hai pt naøy cuõng ñuùng cho ngöôøi dt : böôùc thôøi gian v(i) : vaän toác taïi thôøi ñieåm t=to+t*i v(i-1) : vaän toác taïi thôøi ñieåm t=to+t*(i-1) Giaù trò: B1=10, m1=1000, T=10 (thôøi gian moâ phoûng) B2=0.2, m2=100, gia toác =1.5 , v_max=5 m/s 7. Xaây döïng ma traän 16x16 sau:
Baøi taäp Matlab caên baûn Nguyeãn Hoaøi Sôn 8. Vieát function I=trap(f,a,b,n) duøng tính tích phaân. Vôùi I= Trong ñoù : h=(a-b)/n vaø xi = a+i*h Kieåm tra keát quaû cuûa haøm treân vôùi haøm f=cos(x)+sin(x) cho a=0,b=pi/3 9. Kieåm tra keát quaû baøi treân baèng quy luaät Simpson I= 10. Döïa vaøo löu ñoà sau haõy xaây döïng moät script file duøng caùc leänh input, sprinft 11. Söû duïng voøng laëp for vaø caùc leänh treân ñeå xaây döïng script file döïa treân löu ñoà sau:
Baøi taäp Matlab caên baûn Nguyeãn Hoaøi Sôn 12. Söû duïng voøng laëp while xaây döïng script file döïa treân löu ñoà sau:
Baøi taäp Matlab caên baûn Nguyeãn Hoaøi Sôn 13. Chuyeån ñoåi caùc script file(8,9,10) treân thaønh haøm (function) 14. Veõ hình troøn 15. Veõ caùc maët (surface ) sau: 16. Veõ maët döïa treân phöông trình tham soá sau : Vôùi u,v=0 :2*pi HD : Tính u,v sau ñoù söû duïng leänh surf(X,Y,Z) 14. Veõ maët contour sau :
Baøi taäp Matlab caên baûn Nguyeãn Hoaøi Sôn Vôùi söû duïng haøm contour 17. Döôùi ñaây laø baûn döõ leäu thu ñöôïc töø moät thí nghieäm maø moâ hình toaùn hoïc coù coâng thöùc : Trong ñoù xm,ym laø giaù trò ño ñöôïc coøn δ y laø sai soá taïi giaù trò xm Söû duïng haøm errorbar ñeå theå hieän keát quaû treân. 16. Xaây döïng structure sau :
Baøi taäp Matlab caên baûn Nguyeãn Hoaøi Sôn BÀI TẬP MATLAB CĂN BẢN 1. Viết chương trình tính tiền điện cho một gia đình biết rằng nếu số Kwh điện sử dụng: - < 100 thì giá 1 Kwh là 500đ. 100 ≤ Kwh < 200 thì giá 1 Kwh là 700đ. - 200 ≤ Kwh < 500 thì giá 1 Kwh là 1000đ. - - Kwh > 500 thì giá 1 Kwh là 1200đ. Đồng thời: + Nếu Kwh < 0 thì in ra: không hợp lệ. + Nếu Kwh sử dụng ≥ 500 thì sau dòng in tổng số tiền có ghi kèm theo câu: Hạn chế sử dụng điện. 2. Viết chương trình giải và biện luận bất phương trình : ax > 3x + 6 3. Viết chương trình tính doanh thu trung bình DTTB trong 1 năm của Công ty X theo doanh thu của 4 quý: Q1, Q2, Q3 và Q4, đơn vị tính là Triệu đồng. Nếu DTTB < 100 thì in ra: Phá sản, từ 100 đến nhỏ hơn 200 thì in ra: Doanh thu kém. Từ 200 đến nhỏ hơn 300 thì in ra: Doanh thu trung bình, từ 300 đến nhỏ hơn 500 thì in ra: Doanh thu khá, từ 500 trở lên thì in ra: Doanh thu tốt. 4. Viết chương trình tính n! 5. Viết chương trình tính căn bậc 2 của 1 số: - Lấy xn = x/2 1⎛ x⎞ Sử dụng công thức: xn+1 = ⎜ xn + ⎟ - ⎜ xn ⎟ 2⎝ ⎠ xn +1 − xn
Baøi taäp Matlab caên baûn Nguyeãn Hoaøi Sôn xn +1 − xn 20.0 R, Area, Circ False N=N+1 Stop
Baøi taäp Matlab caên baûn Nguyeãn Hoaøi Sôn 8. Sử dụng vòng lặp while xây dựng script file dựa trên lưu đồ sau: Start N=0 R=0 Area = 0 Circ = 0 OUTPUT While R>0 N INPUT R Area = π * R 2 Circ = 2 * π * R INPUT R If True OUTPUT Area > 20.0 R, Area, Circ False N=N+1 Stop
Baøi taäp Matlab caên baûn Nguyeãn Hoaøi Sôn 9. Viết chương trình tính diện tích hình thang theo đáy lớn a, đáy nhỏ b và chiều cao 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 toá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 vuông 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 vuông và 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 đi từ TPHCM đi Hà Nội: đến quầy số 1, TPHCM đi Đà Nẵng: đến quầy số 2, TPHCM đi Nha Trang: đến quầy số 3, TPHCM đi Vũng Tàu: đến quầy số 4, TPHCM đi Đà Lạt: đến quầy số 5, TPHCM đi Côn Đảo: đến quầy số 6. 12. Viết phương trình giải và biện luận phương trình: ⎧a1 x + b1 y + c1 = 0 ⎨ ⎩a2 x + b2 y + c2 = 0 Theo các tham số a1 , b1 , c1 , a2 , b2 , c2 CÁC PHÉP TÍNH TRÊN MA TRẬN VÀ MẢNG 13. Cho: ⎡2 4 5⎤ ⎢8 12 20 ⎥ A=⎢ ⎥ ⎢4 2 12 ⎥ ⎣ ⎦ a. Tính B = A + 12 b. Tính C = 4*A c. Gõi D là ma trận chuyển trí của C. Xác định D 14. Cho: ⎡2 4 5 ⎤ ⎡10 4 6 ⎤ ⎢8 2 4 ⎥ và B = ⎢8 2 12 ⎥ A=⎢ ⎥ ⎢ ⎥ ⎢ ⎣6 8 10 ⎥ ⎢4 14 20 ⎥ ⎦ ⎣ ⎦ a. Tính ma trận tổng C của A và 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à véctơ hàng B:
Baøi taäp Matlab caên baûn Nguyeãn Hoaøi Sôn ⎡1 2 3 ⎤ A = ⎢4 5 6⎥ B = [2 4 6 8] ⎢ ⎥ ⎢7 8 9 ⎥ ⎣ ⎦ Viết chương trình thực hiện liên tiếp các công việc sau: a. Nhập A và B b. Tạo ma trận con C của ma trận A bằng cách lấy số hạng thứ 1, 2, 3 và 4 c. Tính ma trận tổng D của B và C d. Tính ma trận E là ma trận chuyển trí của B e. Tính F = D*E Viết chương trình sao cho các kết quả được hiển thị trên màn hình. 16. Cho ⎡2 12 15 ⎤ ⎡4 6 5 ⎤ A = ⎢24 6 27 ⎥ B = ⎢12 0 9 ⎥ và ⎢ ⎥ ⎢ ⎥ ⎢18 4 20 ⎥ ⎢0 2 10 ⎥ ⎣ ⎦ ⎣ ⎦ Hãy thực hiện phép chia A./B và hãy cho biết INF nói lên điều gì ? 17. Viết chương trình tính các giá trị y = sin(x), trong đó giá trị của x 0 ≤ x ≤ π và bước nhảy là π / 4 . 18. Hãy sử dụng ones, eye và phương pháp xây dựng ma trận từ các khối để tạo nên ma trận sau: ⎡2 2 2 4 0 0⎤ ⎢2 2 2 0 4 0⎥ ⎢ ⎥ ⎢2 2 2 0 0 4⎥ A= ⎢ ⎥ ⎢4 0 0 2 2 2⎥ ⎢0 4 0 2 2 2⎥ ⎢ ⎥ ⎢0 0 4 2 2 2⎥ ⎣ ⎦ 19. Hãy sử dụng ones, eye và phương pháp xây dựng từ các khối để tạo nên ma trận sau: ⎡1 1 0 0⎤ ⎢1 1 0 0⎥ C =⎢ ⎥ ⎢0 0 1 1⎥ ⎢ ⎥ ⎣0 0 1 1⎦
Baøi taäp Matlab caên baûn Nguyeãn Hoaøi Sôn Biến đổi C thành: ⎡1 0 0⎤ ⎡1 0 0⎤ ⎢1 0 0⎥ ⎥ và E = ⎢1 0 0⎥ D=⎢ ⎢ ⎥ ⎢0 1 1⎥ ⎢0 1 1 ⎥ ⎢ ⎥ ⎣ ⎦ ⎣0 1 1⎦ 20. Giải phương trình sau bằng phương pháp định thức: ⎧15 x +13 y = 440 ⎨ ⎩24 x + 7 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 + 5 x3 + 8 x4 = 884 ⎪5 x + 8 x + 6 x + 10 x = 1248 ⎪1 2 3 4 ⎨ ⎪12 x1 + 6 x2 + x3 = 544 ⎪8 x + x + 12 x + 20 x = 1908 ⎩1 2 3 4 22. Giải hệ phương trình sau bằng phương pháp chia trái ma trận: ⎧5 x1 + 6 x2 + 8 x3 = 82 ⎪ ⎨4 x1 + 7 x2 + 10 x3 = 96 ⎪8 x + 3x + 4 x = 52 ⎩1 2 3 23. Lieät keâ caùc phaàn töû vôùi ñieàu kieä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 3 4 1] a= 11341 ind = (a == 1) ind = 11001 b = a(ind) b= 111 24. Tìm phaàn töû ma traän vaø vectô vaø cho bieát vò trí cuûa chuùng %tim dia chi voi a = 1. ind = find(a == 1) ind =
Baøi taäp Matlab caên baûn Nguyeãn Hoaøi Sôn 125 25. saép xeáp töø nhoû ñeá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= 3 1 4 2 26. taïo haøm vaø tính giaù trò haøm baèng “inline” vaø 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 = 5 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
Baøi taäp Matlab caên baûn Nguyeãn Hoaøi Sôn BÀI TẬP MATLAB CĂN BẢN PHẦN ĐIỆN Baøi 1: cho maïch ñieän nhö hình veõ, tìm ñieän theá nuùt V1, V2, V3 Taïi nuùt 1: Taïi nuùt 2: Taïi nuùt 3: Ta coù heä: Baøi 2: cho maïch ñieän nhö hình veõ, tìm ñieän theá nuùt V1, V2, V3, V4
Baøi taäp Matlab caên baûn Nguyeãn Hoaøi Sôn Taïi nuùt 1: Taïi nuùt 2: Nuùt 2, 3: Taïi nuùt 4: Cuoái cuøng ta coù heä: Baøi 3: Xaùc ñònh giaù trò doøng ñieän chaïy qua ñieän trôû Ta coù:
Baøi taäp Matlab caên baûn Nguyeãn Hoaøi Sôn Voøng 1: Voøng 2: Voøng 3: Ta coù heä phöông trình Baøi 4: Tìm ñieän theá nuùt trong maïch ñieän Phaân tích maïch ñieän ta coù heä phöông trình:
Baøi taäp Matlab caên baûn Nguyeãn Hoaøi Sôn Theá soá ta ñöôïc: Baøi 5: cho maïch ñieän nhö hình veõ: Töø maïch ñieän ta coù phöông trình: Cho C = 10 μ F söû duïng Matlab ñeå veõ ñieän theá ñi qua tuï ñieän vôùi giaù trò cuûa ñieän trôû R = 1 k Ω, R = 10 k Ω, R = 0.1 k Ω Baøi 6: Cho maïch ñieän RLC nhö hình veõ. Neáu di ( 0 ) L = 10 H , R = 400 Ω, C = 100 μ F , vs ( t ) = 0, i ( 0 ) = 4 A, = 15 A / s dt Tìm i ( t ) Cho Baøi 7: Cho maïch ñieän nhö hình veõ. Neáu taïi thôøi ñieåm ban ñaàu t = 0 coâng taéc môû. Tìm ñieän theá v ( t ) . Vôùi L = 1 / 32 H , R = 10 Ω, C = 50 μ F , I s = 2 A
Baøi taäp Matlab caên baûn Nguyeãn Hoaøi Sôn Cho