BÁO cáo bài tập lớn xác ĐỊNH QUỸ đạo của vật GIẢI QUYẾT bài TOÁN BẰNG MATLAB

    ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA KHOA KHOA HỌC ỨNG DỤNG BÁO CÁO BÀI TẬP LỚN XÁC ĐỊNH QUỸ ĐẠO CỦA VẬT Giáo viên giảng dạy: TRƯƠNG VĂN MINH NHÓM LỚP L13, NHÓM 1: PHAN TRƯỜNG AN ĐẶNG THÁI DƯƠNG LẠI NGUYỄN DUY NGUYỄN TRI BÌNH Tp HCM, 01/2021 II   ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA KHOA KHOA HỌC ỨNG DỤNG BÁO CÁO BÀI TẬP LỚN XÁC ĐỊNH QUỸ ĐẠO CỦA VẬT Nhóm 1: Thành viên PHAN TRƯỜNG AN ĐẶNG THÁI DƯƠNG LẠI NGUYỄN DUY NGUYỄN TRI BÌNH MSSV 2010826 2011031 2010999 2010936 II   LỜI CÁM ƠN Trong suốt q trình thực tiểu luận nói trên, nhóm nhận nhiều quan tâm ủng hộ, giúp đỡ tận tình thầy cơ, anh chị em bè bạn  Ngồi ra, nhóm xin gửi lời tri ân chân thành đến thầy Trương Văn Minh, giảng viên hướng dẫn cho đề tài matlab Nhờ có thầy hết lịng bảo mà nhóm hồn thành báo cáo nghiên cứu tiến độ giải tốt vướng mắc gặp phải Sự hướng dẫn thầy kim nam cho hành động nhóm phát huy tối đa mối quan hệ hỗ trợ thầy trị mơi trường giáo dục Lời cuối, xin lần gửi lời biết ơn sâu sắc đến cá nhân, thầy cô dành thời gian dẫn cho nhóm Đây niềm tin, nguồn động lực to lớn để nhóm đạt kết III   TÓM TẮT BÀI BÁO CÁO Có thể nói việc xác định quỹ đạo vật không gian vấn đề vô quan trọng nhiều cơng trình nghiên cứu vấn đề thực tiễn ngày Quỹ đạo chuyển động vật đường vạch vật chuyển động Khi chuyện chuyện động, vật thể chịu ảnh hưởng lực tác động không Ở đề tài này, nghiên cứu vị trí chuyển động vật mặt phẳng Oxy cho vecto bán kính kính Với đại lượng tọa độ x, y cho hàm lượng giác từ vecto bán kính Các đại lượng tọa độ biến đổi theo thời gian tạo nên thay đổi vị trí MATLAB công cụ đặc biệt giúp xác định quỹ đạo chuyển động vật cách đơn giản xác Từ dễ dàng đưa kết luận phương trình, quỹ đạo yếu tố liên quan QUỸ ĐẠO CHUYỂN ĐỘNG CỦA VẬT trường hợp đặc biệt IV   DANH MỤC CÁC HÌNH ẢNH Hình 0.1 Logo trường I Hình 1.2 Hệ tọa độ Đề-các Hình 1.3 Vectơ vị trí trong hệ tọa độ 22 Hình 1.4 Đồ thị hàm sin Hình 1.5 Đồ thị hàm cos cos 66 Hình 3.1 Code toán matlab 11 Hình 3.2 Kết chương trình cho trường hợp 12 Hình 3.3 Kết chương trình cho trường hợp 12 Hình 3.4 Kết chương trình cho trường hợp 13 Hình 3.5 Kết chương trình cho trường hợp 14 V   MỤC LỤC LỜI CẢM ƠN IV IV TÓM TẮT BÁO CÁO DANH MỤC CÁC HÌNH ẢNH V MỤC LỤC .VI THUYẾT .1 CHƯƠNG CƠ SỞ LÝ THUYẾT 1.1 Khái niệm mở đầu 1.1.1 Hệ tọa độ 1.1.2 Chuyển động, động, không gian, thời gian, hệ quy chiếu, chất điểm điểm .11 1.1.3 Vectơ vị trí 1.1.4 Phương trình chuyển động 1.1.5 Quỹ đạo phương phương trình quỹ đạo đạo 33 1.2 Chuyển động chất điểm 1.2.1 Chuyển động thẳng 1.2.2 Chuyển động elip .4 1.2.3 Chuyển động tròn 55 1.2.4 Chuyển động parabol 1.2.5 Chuyển động động hyperbol hyperbol .55 1.3 Hàm lượng giác giác 1.3.1 Hàm lượng giác sin .6 1.3.2 Hàm lượng giác cos 1.3.3 Một số công thức lượng giác .6 TOÁN 77 CHƯƠNG PHÂN TÍCH BÀI TỐN 2.1 Yêu cầu toán toán 77 2.2 Input Output VI   2.2.1.Input VI   2.2.2 Output 2.3 Cách làm 2.2.2 Cách làm thông thường thường 77 2.3.2 Sử dụng công cụ matlab matlab 88 CHƯƠNG GIẢI QUYẾT BÀI TOÁN BẰNG MATLAB 3.1 Giới thiệu thiệu Matlab lệnh sử dụng dụng 3.1.1 Giới thiệu Matlab Matlab .99 3.1.2 Các lệnh Matlab sử dụng 3.1.2.1 SYMS SYMS 3.1.2.2 IF-ELSE 3.1.2.3 EZPLOT EZPLOT 3.1.2.4 CÁC LỆNH NHẬP XUẤT .10 3.1.2.5 LỆNH CLC .10 3.1.2.6 LỆNH CLEAR 10 3.1.2.7 LỆNH GRID 10 3.2 Giải toán 10 3.3 Ví dụ .11 KẾT 15 15 CHƯƠNG TỔNG KẾT TÀI LIỆU THAM KHẢO 16 .?? PHỤ LỤC (nếu có) VII    CHƯƠNG CƠ SỞ LÝ THUYẾT 1.1 Khái niệm mở đầu 1.1.1 Hệ tọa độ Vì chuyển động xảy khơng gian thời gian nên để mô tả chuyển động trước tiên ta phải tìm cách định vị vật vị vật không gian Muốn vậy, ta  phải đưa thêm vào hệ qui chiếu hệ tọa độ. Trong độ. Trong vật lý người ta sử dụng nhiều hệ tọa độ khác Ở ta dùng hệ tọa độ Descartes độ  Descartes *Hệ tọa độ Đề-các Đề -các (Descar (D escartes) tes) : Hệ tọa độ Đề-các gồm ba trục Ox, Oy, Oz tương ứng vng góc với đơi một, chúng tạo thành tam diện thuận Điểm O gọi gốc tọa độ Vị trí điểm M hồn tồn xác định bán kính vectơ , hay tập hợp  ba số (x,y,z) (x,y,z) x,y,z hình hình chiếu ccủa điểm mút mút M vect vectơơ = O lên trục tương ứng Ox, Oy, Oy, Oz gọi là ba tọa độ của độ của điểm M hệ tọa độ Đề - các. Nguồn: Internet   Hình 1.2 Hệ tọa độ Đề-các Đề-các  Nếu gọi i ,  j , k  là  là vectơ đơn vị hướng theo trục Ox, Oy, Oz ta viết : ⃗ ⃗ ⃗   r = x i +  y  j + z k  ⃗ x, y, z llàà hình hình chiếu ⃗ ⃗   ⃗ bbaa trục Ox, Oy, Oy, Oz gọi gọi là thành phần số 1.1.2 Chuyển động; không gian; thời gian; hệ quy chiếu; chất điểm -Chuyển động vật thay đổi liên tục vị trí vật theo thời gian   -Hệ vật quy ước đứng yên dùng để xác định vị trí vật khác chuyển động gọi hệ quy chiếu Người ta thường gắn điểm gốc hệ trục tọa độ vào hệ quy chiếu, hệ trục tọa độ gọi hệ quy chiếu -Để xác định thời gian chuyển động vật, người ta gắn vào hệ quy chiếu đồng hồ, vật chuyển động vị trí thay đổi theo thời gian -Nếu vật có kích thước nhỏ, khơng đáng kể so với khoảng cách chúng kích thước vật khác mà ta xem xét, kích thước vật bỏ qua, ta có khái niệm chất điểm Tập hợp chất điểm gọi hệ chất điểm 1.1.3 Vectơ vị trí:   Để xác xác định định vị ttrí rí của một chất chất điểm điểm M tron trongg không gian, người ta thường gắn vào hệ quy chiếu hệ trụ tọa độ, hệ tọa độ thường dùng hệ tọa độ Descart Descartes es với ba ba trụ Ox, Oy Oz vng vng góc với đôi một, hợp thành tam diện thuận Vị trí điểm M hồn tồn xác định ta xác định thành phần x, y, z vectơ vị trí  Nguồn: Slide giảng thầy  Nguyễn Minh Châu  Hình 1.3 Vectơ vị trí trong  hệ tọa độ  ⃗OM =r  x ( x , y , z ) ¿  gọi bbánán kính vectơ được vẽ từ gốc hệ tọa độ đến vị ttrí rí chất chất ⃗ điểm M) *Vectơ vị trí chất điểm: r = x i +  y  j + z k  ⃗ ⃗ ⃗   ⃗   1.1.4 Phương trình chuyển động Để xác định chuyển động chất điểm cần biết vị trí chất điểm thời điểm khác Nói cách khác, cần biết phụ thuộc theo thời  gian  của bán kính vectơ  gian chất điểm : {  x = f 1( t ) r =r ( t )=  y = f 2 ( t )  z = f 3 ( t ) ⃗ ⃗ Phương trình biểu diễn vị trí chất điểm theo thời gian gọi phương trình chuyển động chất điểm Trong hệ tọa độ Đề-các, phương trình chuyển động chất điểm hệ gồm ba  phương trình trình : x = x(t);y = y(t) ;z = z(t) Ví dụ sau phương trình chuyển động chất điểm hệ tọa độ Đề-các : {   x = x cosωt   y = y cos  ( ωt + φ )  z =0 1.1.5 Quỹ đạo phương trình quỹ đạo: Khi chuyển động, vị trí chất điểm thời điểm khác vạch không gian đường cong liên tục gọi quĩ đạo của đạo của chuyển động Phương trình mơ tả đường cong quĩ đạo gọi là phương là phương trình quĩ đạo Trong hệ tọa độ Đề-các  phương trình trình quĩ đạo đạo có dạng :  f (x,y,z) (x,y,z) = C đó f  đó f  là  là hàm tọa độ x, y, z C số Về nguyên tắc, ta biết phương trình chuyển động khử tham số t  ta  ta tìm mối liên hệ tọa độ x, y, z tức tìm phương trình quĩ đạo   Vì vậy, đơi người ta cịn gọi phương trình chuyển động là phương là phương trình quĩ đạo cho dạng tham số Có dạng khử t phương trình chuyển động chất điểm: Dạng 1: Phương pháp Dạng 2: sin & cos theo t, áp dụng số công thức lượng giác   1.2 Chuyển động chất điểm - Bằng cách khử tham số phương trình chuyển động, ta tìm phương trình quỹ đạo từ đối chiếu với phương trình quỹ đạo chuyển động thẳng, elip, tròn, parabol, hyperol học 1.2.1 Chuyển động thẳng  - Quỹ đạo đường thằng - Phương trình quỹ đạo đường thẳng   y=ax +b 1.2.2 Chuyển động elip - Quỹ đạo hình elip - Trong tốn học, học, hình elíp đường cong phẳng phẳng xung  xung quanh hai tiêu điểm, điểm, cho với điểm đường cong, tổng khoảng cách đến hai tiêu điểm số - Phương trình tắc elip:  x a  y b + =1 ( ∀ a , b ≠ , 1)    1.2.3 Chuyển động tròn -Quỹ đạo đường trịn bán kính R Trong hìn hình học phẳn phẳngg, đường ờng trò rònn (hoặ vòng òng tr tròòn) tập hợp tất điểm trên điểm trên mặt phẳng, phẳng, cách điểm cho trước khoảng cách nào Điểm cho trước gọi tâm đường tròn, khoảng cho trước gọi cách nào là bán kính của kính của đường trịn Đường trịn định nghĩa hình elíp đặc elíp đặc biệt với hai tiêu điểm trùng tâm sai bằng sai bằng Đường tròn tâm O bán kính R ký hiệu (O;R) Phương trình tắc đường trịn: 2   x  + y = A 1.2.4 Chuyển động parabol - Quỹ đạo parabol - Phương trình quỹ đạo:   y = ax2 + bx + c (∀ a ≠ ¿ 1.2.5 Chuyển động hyperbol  - Quỹ đạo hyperbol - Phương trình tắc hyperbol: 1.3 Hàm lượng giác học,, hàm tuần hoàn là hoàn là hàm số số lặp  lặp lại giá trị *  Trong toán học khoảng đặn hay chu kỳ Ví dụ quan trọng hàm tuần hồn đó hồn đó hàm lượng giác, giác, mà lặp lại khoảng 2π radian radian Hàm tuần hoàn sử dụng     thường xuyên để miêu tả dao động, động, sóng, sóng, tượng khác thể tính tuần hồn hồn 1.3.1 Hàm lượng giác sin: - Hàm sin góc định nghĩa tam giác vng tỷ lệ cạnh vng góc đối diện chia cho cạnh huyền Đồ thị gọi đồ thị hàm sin Hàm sin định nghĩa khoảng từ -∞ đến ∞ có giá trị từ -1 đến - Hàm sin có tính chất tuần hồn  Nguồn: Internet   Hình 1.4 Đồ thị hàm sin 1.3.2 Hàm lượng giác cos: - Hàm cos góc định nghĩa tam giác vng tỷ lệ cạnh kề vng góc chia cho cạnh huyền Đồ thị gọi đồ thị hàm cos Hàm cos định nghĩa khoảng từ ∞ đến ∞ có giá trị từ -1 đến - Hàm cos có tính chất tuần hồn.  hồn.   Nguồn: Internet  1.3.3 Một số công thức lượng giác: sin cos a + cos b=1 π  π  cos − α  = α − = sinα   ( ) ( ) ( + )=− cos  Hình 1.5 Đồ thị hàm cos π  α  sinα  cos ( α + π )=−cos α  cos ( α + π  ) )=cos α  cos ( a + b ) =cos a∗cos b −sin a∗¿ sin b ¿   CHƯƠNG 2.1 PHÂN TÍCH BÀI TỐN u cầu toán Sử dụng Matlab để giải toán sau: “Vị trí chất điểm chuyển động mặt phẳng Oxy xác định vectơ bán kính vật?” 2.2 Cho trước giá trị x0, y 0 và φ, xác định quỹ đạo Input output Các giá trị x0, y0 và φ 2.2.1 Input  Các giá trị x0, y0 và φ 2.2.2 Output  - Thiết lập phương trình tương ứng:  x = x 0∗ cos( t )  y = y 0∗cos ( t + φ ) - Đưa phương trình chuyển động kết luận quỹ đạo: Quỹ đạo vật đoạn thẳng, elip, đường tròn tùy thuộc vào giá trị input - Vẽ hình quỹ đạo: Sử dụng matlab để vẽ   2.3 Cách làm 2.3.1 Cách làm thông thường  - Để có hình vẽ quỹ đạo vật Chúng ta phải thiết lập phương trình tương ứng từ giá trị x 0, y0 và φ - 0 Xét các Nếu TH xđặc biệt với giá trị x , y  vàthẳng φ: nằm Oy chạy từ -  y0 = ta có quỹ đạo đoạn  đến + y0 Tương tự với y0 = Với giá trị đặc biệt φ ta dễ dàng rút phương trình quỹ đạo y = f(x), từ dễ dàng có hình vẽ quỹ đạo - Với TH mà giá trị x 0, y0 và φ không dễ dàng để có phương trình quỹ đạo việc vẽ hình vơ khó khăn   2.3.2 Sử dụng cơng cụ matlab Chúng ta dễ dàng có hình vẽ quỹ đạo với cơng cụ tính tốn MATLAB Sẽ trình bày chương sau    CHƯƠNG GIỚI THIỆU CHUNG 3.1 Giới thiệu Matlab lệnh sử dụng: 3.1.1 Giới thiệu Matlab: MATLAB là phần mềm cung cấp mơi trường tính tốn số và MATLAB là số và lập trình, cơng ty MathWorks  MathWorks th thiế iếtt kế kế MATLAB MATLAB   cho phép tính tốn số với ma trận, số hay biểu đồ thông tin, thực thuật toán vẽ đồ thị  thị hàm số hay - Matlab vừa ngơn ngữ lập trình vừa phần mềm ứng dụng tính tốn hiệu - 3.1.2 Các lệnh Matlab sử dụng: 3.1.2.1 SYMS: -   Lệnh biến đổi số, biến hay đối tượng thành kiểu Symbolic VD: >> syms x y t; 3.1.2.2 IF-ELSE:   - Lệnh giúp thực thi đoạn lệnh thỏa điều kiện cho trước   - Cú pháp: if (Biểu thức điều kiện 1) Lệnh thực thi ĐK elseif (Biểu thức điều kiện 2) Lệnh thực thi ĐK End 3.1.2.3 EZPLOT: - Lệnh Lệnh  vẽ đồ thị hàm theo biến khoảng muốn - Cú pháp:  Vẽ hàm y=f(x): ezplot (x, f, [khoảng x]); VD: syms x; ezplot (x^2+1, [0,2])  Vẽ hàm tham số: ezplo t(x(t), y(t), [t1,t2]) VD: syms t; x = t; y = t^2; ezplot(x ,y, [0,2])   3.1.2.4 CÁC LỆNH NHẬP XUẤT: - Lệnh nhập VD: Nhập x >>x=input(‘Nhap x=’) - Lệnh xuất: VD: Xuất hình dịng “Hello World” >>disp(‘Hello World’) 3.1.2.5 LỆNH CLC: - Là lệnh xóa cửa sổ lệnh - Cú pháp: clc 3.1.2.6 LỆNH CLEAR: - Là lệnh để xóa đề mục nhớ - Cú pháp: clear  3.1.2.7 LỆNH GRID: - Là lệnh để tạo lưới tọa độ - Cú pháp:  grid on (để hiển thị lưới tọa độ)  grid off (dể tắt lưới tọa độ) 3.2 Giải toán: Sau code để giải toán matlab: clc; % Xoa cua so lenh clear; % Xoa tat ca cac bien luu lai bo nho syms x y t ; % Khai bao cac bien thi x0 = input('Nhap input('Nhap x0: ''); ); % Nhap x0 y0 = input('Nhap input('Nhap y0: ''); ); % Nhap y0  phi = input('Nhap input('Nhap goc phi: '); '); % Nhap goc phi disp('Phuong disp( 'Phuong trinh chuyen dong cua vat la: '')) % Xuat dong chu man hinh   x = x0 * cos(5*t) y = y0 * cos(5*t + phi) % Dua phuong trinh x,y theo t   ezplot(x, y); % Ve thi x,y theo t grid on  on  % Hien thi luoi toa 10     if  (x  (x == y) && (x == 0) 'Vat KHONG chuyen dong '); disp('Vat disp( '); elseif  (0  (0 == mod(phi/pi,1)) || (x == 0) || (y == 0)   disp('Quy disp('Quy dao cua vat la doan thang '); ' ); elseif  (0  (0 == mod(0.5*((2*phi/pi)-1),1)) && (x0 == y0) 'Quy dao cua vat la duong tron '); disp('Quy disp( ' ); else disp( else disp('Quy 'Quy dao cua vat la mot duong elip ''); );   % Xet dieu kien de kiem tra hinh dang quy dao cua vat   end end   % Ket thuc chuong trinh 3.2.Ví dụ: -  Nhập code code vào matlab: matlab:  Hình 3.1 Code toán toán matlab 11    - Cho chạy với tham số: x0=4, y0=3, φ= π/3 chương trình chạy xuất hình phương trình x,y theo t hình dạng quỹ đạo chất điểm đồ thị hình elip:  Hình 3.2 Kết chương trình cho trường hợp   - Với trường hợp khác x0=4, y0=4, φ= π/2 quỹ đạo hình trịn 12  Hình 3.3 Kết chương trình cho trường trường hợp   - Với trường hợp x0=7, y0=2, φ= π quỹ đạo đoạn thẳng:  Hình 3.4 Kết chương trình cho trường hợp - Với trường hợp đặc biệt x0=0, y0=0, φ= π/6 chẳng hạn vật khơng chuyển động nên đồ thị điểm:  Hình 3.5 Kết của chương trình cho trường trường hợp 13   CHƯƠNG TỔNG KẾT Bài báo cáo hoàn thành xong việc giải toán từ việc đặt vấn đề chung, nêu kiến thức bản, yêu cầu vấn đề đặt ra, cách thức giải vấn đề Qua việc tìm hiểu thực hành, chúng em rút vấn đề trọng tâm thể qua báo cáo Qua nắm vững cách sử dụng phần mềm Matlab để giải toán, nắm vững kiến thức cách giải vấn đề khác đời sống Bằng việc giải toán từ vấn đề chung đến vấn đề cốt lỗi, phần mềm Matlab thể nhiều ưu điểm giúp việc giúp chúng em tiếp cận, khảo sát giải toán trở nên sinh động, trực quan lý thú 14   TÀI LIỆU THAM KHẢO [1] Bộ môn Vật lý Ứng dụng: Vật lý đại cương A1 (Tài liệu lưu hành nội bộ), Trường ĐH Bách Khoa – ĐH Quốc gia TP HCM, 2009 [2] Tài liệu hướng dẫn sử dụng phần mềm Matlab : http://hnue.edu.vn/Portals/0/TeachingSubject/tungpk/0e9a673a-ae76-45e3-99e9 bb923c75e577Thuc bb923c75 e577Thuc-hanh-MatLab hanh-MatLab.pdf .pdf https://uet.vnu.edu.vn/~tantd/dieukhien/Chuong%201%20Matlab%20co%20ban.p Df [3] Khái quát chuyển động học - http://www4.hcmut edu.vn/~huynhqlinh/ edu.vn/~huynhqlinh/   VLDC1/index.html [4] Chương : Động học chất điểm - Slide giảng thầy Nguyễn Minh Châu [5] Hàm lượng giác - http:// Wikipedia.com 15 ... thay đổi vị trí MATLAB công cụ đặc biệt giúp xác định quỹ đạo chuyển động vật cách đơn giản xác Từ dễ dàng đưa kết luận phương trình, quỹ đạo yếu tố liên quan QUỸ ĐẠO CHUYỂN ĐỘNG CỦA VẬT trường hợp... to lớn để nhóm đạt kết III   TĨM TẮT BÀI BÁO CÁO Có thể nói việc xác định quỹ đạo vật không gian vấn đề vô quan trọng nhiều cơng trình nghiên cứu vấn đề thực tiễn ngày Quỹ đạo chuyển động vật. .. HỌC QUỐC GIA TP HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA KHOA KHOA HỌC ỨNG DỤNG BÁO CÁO BÀI TẬP LỚN XÁC ĐỊNH QUỸ ĐẠO CỦA VẬT Nhóm 1: Thành viên PHAN TRƯỜNG AN ĐẶNG THÁI DƯƠNG LẠI NGUYỄN DUY NGUYỄN