zunia.vn

Tuyển sinh 2024 dành cho Gen-Z

zunia.vn

» Công Nghệ Thông Tin

» Kỹ thuật lập trình

Matlab - Tin học ứng dụng - Bài tập

Chia sẻ: Nguyen Nhi | Ngày: | Loại File: PDF | Số trang:6

Báo xấu

268
lượt xem 63
download

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

Tài liệu tham khảo bài giảng Matlab - Tin học ứng dụng ( Nguyễn Bá Tuyên & Nguyễn Quang Chiến ) dành cho sinh viên chuyên ngành kỹ thuật biển gồm 7 chương và phần bài tập - Bài tập

Chủ đề:

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

Lưu

Nội dung Text: Matlab - Tin học ứng dụng - Bài tập

Bài giả g Tin Ứ g dụ g 2 – Matlab 7 n n n 8. PHẦ BÀI TẬ N P Bài tậ số : p1 Tính toán sóng theo số iệ sau: lu o T = 8 s, H = 2 m, α= 30 . o o •G嫕 gi trịho 3 bi ế T, Ho , al ph a0. c n •T嫕 Co. h •T嫕 L h o •Mở1 fil e d isp erse.m b ằg cách gõ>> edit disperse n •Vớ fi l e mớ mởhãy gõ vào nộ du n g sau: fu n ct io n [L] = d isp er se(h, Lo) i i i L = Lo; err = Inf; while err > 0.1; Ltry = Lo * tanh(2*pi*h / L); err = abs(Ltry - L); L = Ltry; end; •Ta đ thiếlậ đ ợ hàm tính L từcác gi á t rị và Lo cho trư c t p ưc h ớ ã •Hãy thiếl ậ mộ vec-t ơ chứ đ sâu kh ác nh au : 3 m, 3 .2 m, . .., 4 m. Tín h L tạ tp t aộ i h các đ s âu cho b ởh . ộ i •T嫕 C tư ng ứg vớ L h ơn i •Tính sin_alpha theo công thức sinα sinα* C / Co(chú ý đn vịộ ) ơ đ! = o •T嫕 al p ha (d 嫕 lệh asi n) h gn •T嫕 Kr = sq rt (cosα/ cos α h ) o •T嫕 Ksh t heo c嫕g t hứ K sh = h c •Tính H = Ho × Ksh × Kr •Ghi lạc嫕gi i trịủ H ca - 61 -
Bài giả g Tin Ứ g dụ g 2 – Matlab 7 n n n •Viếkếqu ảa f il e ketq ua. tx t tt r Bài tậ số p2 Quan hệ iữ vậ chu yể b ùn cát S và lư tố u có q uan hệ ạg: S = aub . Đ x ác ga n n uc dn ể đ h các hệ ố và b ngưi t a t iế hành thí nghi ệ và thu đ ợ kếqu ả au : ị n sa ờ n m ưc t s u (m/s) 0.2 0.35 0.57 0.68 0.81 0.96 1.12 S 0.0002 0.0018 0.0159 0.0282 0.0609 0.1258 0.2858 (m3 /s.m) 1) Hãy nhậ các gi á tr ị và S vào hai vec -t ơ sau đ tính X và Y là logarit tư ng p ơ u , ó ứg củ u và S . n a 2) Vẽ ồh ịi ể củ X và Y. đ t đm a 3) Hai hệ số a và b đợ x ác đ h bằg cách d ự v ào quan hệ ưc ịn n a : log S = log a + b log u , hay Y = b X + log a. Sửdụg h àm polyval đ tí nh b và log a trong phư ng trình trên như n ể ơ sau: p = polyval(X, Y, 1) Phầ tửđu củ vec-t ơ chính là b; phầ tửthứ2 củ p bằg l og a. n ầa n a n p Hãy tính a. 4) Vẽ ồ h ịư n g t hẳg : Y = b X + log a lên cùng hệrụ v ớ các đ m chấ ở đ t đờ n tc i iể m câu 2). Bài tậ s ố p3 1) Mộ ch uỗ sối ệ đ đc vậ tố dòng chả đợ p hát si nh bở: t i luoạ n c y ưc i t = 0:0.5:48; N là dãy số g ẫ n hi ên có chiề d ài bằg l en gt h( t) nu u n u = 0.4 + 0.12 sin(2π / 24) + 0.05 N t 2) Tính ứg suấti ế tạven bờ bi ế n t pi , t 0.5 C = 60 m /s τ ρ u2 / C2 =g Thiếl ậ su bp l ot 3 hàng × 1 cộ Hai vùng trên vẽ ai đ th ị ~ t, τ t. tp t. h ồu ~ 3) Nế bờcó t hành ph ầ đt sét vớ ứ g su ất i ế tớ hạ τ= 0.65 N/m2, hãy ch ỉ u nấ in t p i nc ra xem có bao nhiêu thờ đ m x uấhi ệ τ τ. i iể t n >c - 62 -
Bài giả g Tin Ứ g dụ g 2 – Matlab 7 n n n 4) Công thứ x ói l ở ư n g b ờ ư c xác đ h b ở c đờ đợ ị n i: Hãy tính khoảg cách x ói l ở . n E E = 10 (τ– τ n ế τ τ u >c c) trư n g h ợ còn lạ ờ p i =0 Vẽ ~ t lên vùng đ th ịư i cùng tạ bở su bpl o t. E ồ dớ oi Bài tậ số : p4 Cho véctơ = [10 20 30] X và Y = [10 40 30]. 1. Tính diệ tí ch t am gi ác t ạ b ở 3 đ h có tọ đ X( 1)Y (1 ), X (2 )Y( 2), n o i ỉ n aộ X(3)Y(3). 2. Tổg q uát h ơ tính diệ t ích đ gi ác n cạh (t rên mặp hẳg 2 ch iề) vớ tọ đ n n a n t n u i aộ n, xi và yi củ đ h th ứ ch o bở phầ t ử h ứ củ véct ơX và Y tư ng ứ g (d o đ X aỉ n i i ntia ơn ó và Y là 2 véctơ cùng chiề d ài = n). u có LỜI GIẢI Bài tậ số p 1: Hư ng dẫ: t ạ mộ script file có tên baitap1.m, soạ t hả nộ d un g d ưi đy , v à ớ no t n oi ớâ ghi lạ Sau đ chạ fi le n ày, ta sẽhu đ ợ kếqu ả i. ó y t ưc t . % baitap1 % Inputs: T, H0, alpha0 % Outputs: C0, L0, L %------------------------------------------------------------- clear all; g=9.81; T=8; H0=2; alpha0=30; %------------------------------------------------------------- L0=g*T^2/(2*pi) C0 = L0/T; h=[3:0.2:4] % n=length(h); % for i = 1:n - 63 -
Bài giả g Tin Ứ g dụ g 2 – Matlab 7 n n n % L(i)=disperse(h(i), L0); % end; L = disperse(h,L0); % Wave length C=L/T % Wave celerity sin_alpha=sin(alpha0)/C0*C alpha = asin(sin_alpha) % arcsin Kr=sqrt(cos(alpha0)./cos(alpha)) % Refraction coefficient k=2*pi./L % Wave number Ksh=sqrt(1./tanh(k.*h)/(1+2*k.*h/sinh(2*k.*h))) % Shoaling coefficient H=H0*Ksh.*Kr % Wave height Bài tậ số p 2: Hư ng dẫ: t ạ mộ scri pt fi l e có tên baitap2.m , soạ t hả nộ d un g d ưi đy , v à ớ no t n oi ớâ ghi lạ Sau đ chạ fi le n ày, ta sẽhu đ ợ kếqu ả i. ó y t ưc t . (Bài này làm theo các bư c n hưtrên lớ: t í nh lo g( u), l og (S ), vẽlên trụ t ọ đ ớ p caộ thư ng (khô ng vẽên trụ lo ga vì ở â y ta đ tính giá trịủ l o ga rồ , và xác đ h ờ l c đ ca i) ịn ã các hệ ố , b củ đ ờg t hẳg mộcá ch g ầ đng t rên đ th ị sa a ưn n t nú ồ ). Bài tậ số p 3: Hư ng dẫ: t ạ mộ scri pt fi l e có tên baitap3.m , soạ t hả nộ d un g d ưi đy , v à ớ no t n oi ớâ ghi lạ Sau đ chạ fi le n ày, ta sẽhu đ ợ kếqu ả i. ó y t ưc t . clear all; t=0:0.5:48 N=rand(1,length(t)); u=0.4+0.12*sin(2*pi*t/24)+0.05*N plot(N); hold on; plot(u,'+'); C=60; Rho = 1000; g=9.81; To = Rho*g/C^2*u.^2 plot(t,To) hold on grid on Toc = 0.65 10 * (To - Toc) .* (To > Toc) + 0 .* (To
Bài giả g Tin Ứ g dụ g 2 – Matlab 7 n n n function Sdagiac = dientich(x,y); % Tinh dien tich da giac n = length (x); m = length (y); if m ~=n 'Error: 2 vecto X va Y can co chieu dai bang nhau!' else Sdagiac = 0; for i = 1:(n-1) Sdagiac = Sdagiac + 0.5*(y(i)+y(i+1))*(x(i+1)-x(i)); end Sdagiac = Sdagiac + 0.5*(y(n)+y(1))*(x(1)-x(n)); end 2. Vừ rồ t a đ t ạ ra mộ hàm mớ t ro n g M at l ab , vớ tên là hàm dientich. Sử ai ão t i i dụg h àm này ta có thể ễ àng tính đợ d iệ tí ch đ gi ác nó i chu ng và tam giác n dd ưc n a nói riêng bằg cách gọ hàm từ ử sổ h ậ l ệh n i ca n p n  >> dientich(X,Y) Lư ý r ằg hàm tính di ệ tí ch ngầ đ h r ằg 2 v éc tơX, Y ở ây đ ch ứ số iệ n n mị n n đã a lu u cho trư c vềọ đ các đ h củ đ gi ác th eo nh ư ớ ta ộ ỉ aa n quy ưc. H ơ n ữ, khi nh ậ t ọ ớ na pa đ ta đnh số ác đ h từ1 đn n t heo chi ề k im đn g hồ ộá c ỉ n ế u ồ . * Đ t í nh d iệ t ích tam gi ác S v ớ ch i ề dài 3 cạh ch o t rư c l à a, b, c, b ạ có t h ể ể n i u n ớ n thử huật oán kh ác n hưsau: tt S  p ( p a)( p  )( p  ) b c vớ p=(a+b+c)/2 là mộ nử ch u v i. i ta - 65 -
Bài giả g Tin Ứ g dụ g 2 – Matlab 7 n n n 9. TÀI LI Ệ THAM KHẢ : U O 1. Matlab R14 - Helps & Demos, The Mathworks, Inc., 2004 2. David F. Griffths, An introduction to Matlab, Department of Mathematics- The University of Dundee – Sweden, 1996 3. John M. Stockie, A Whirlwind Tour of MATLAB for Students of CS 3113, Department of Mathematics and Statistics, University of New Brunswick - Canada, 2003 4. Bill Mason, Introduction to Matlab, Northeastern University - College of Computer and Information Science – USA, 2003. 5. Tobin A. Driscoll, Crash course in MATLAB , Department of Mathematical Sciences - University of Delaware – USA, 2006 6. Nguyễ Hoàng Hả& N guy ễ ViệAn h, Lậ trình Matlab và ứ g dụg , n i n t p n n NXB Khoa họ và Kỹ hu ậ– Hà Nộ 20 05 c tt i, 7. Nguyễ Phùng Quang, Matlab & Simulink dành cho k ỹ ư i ề kh i ể t ự n s đu n đn g, NXB Khoa h ọ & Kỹh uậ– Hà N ộ 2 00 6 ộ c tt i, 8. Nguyễ Ph ư ng Thả, Programming in Matlab, handouts, 2007 n ơ o 9. Đ Thủ Lợ, Matlab version 7.0 cơ ả, Khó a ti n họ nâ ng cao cho cán yi bn c H. bộ i ảg dạ củ d ự n 9 5 b ộ N&PTNT, handouts, 2007 gn ya á N 10. Knight A., Basics of MatLab® and beyond, CRC Press, 2000 11. Timothy A. Davis & Kermit Sigmon, Matlab ® Primer, Chapman & Hall/CRC 12. R. J. Braun, Beginning Matlab Exercises, Department of Mathematical Sciences – University of Delaware – USA. - 66 -

CÓ THỂ BẠN MUỐN DOWNLOAD

THÔNG TIN

TRỢ GIÚP

HỖ TRỢ KHÁCH HÀNG

Theo dõi chúng tôi

Chịu trách nhiệm nội dung:

Nguyễn Công Hà - Giám đốc Công ty TNHH TÀI LIỆU TRỰC TUYẾN VI NA

LIÊN HỆ

Địa chỉ: P402, 54A Nơ Trang Long, Phường 14, Q.Bình Thạnh, TP.HCM

Hotline: 093 303 0098

Email: support@tailieu.vn

Giấy phép Mạng Xã Hội số: 670/GP-BTTTT cấp ngày 30/11/2015 Copyright © 2022-2032 TaiLieu.VN. All rights reserved.