XÁC ĐỊNH QUĨ ĐẠO CỦA CHẤT ĐIỂM KHI CÓ PHƯƠNG TRÌNH VẬN TỐC

ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA  BÁO CÁO BÀI TẬP LỚN MÔN VẬT LÝ 1 ĐỀ TÀI 2 XÁC ĐỊNH QUĨ ĐẠO CỦA CHẤT ĐIỂM KHI CÓ PHƯƠNG TRÌNH VẬN TỐC GVHD Ths Nguyễn Ngọc Quỳnh Phan Ngọc.

ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA - - BÁO CÁO BÀI TẬP LỚN MÔN VẬT LÝ ĐỀ TÀI 2: XÁC ĐỊNH QUĨ ĐẠO CỦA CHẤT ĐIỂM KHI CÓ PHƯƠNG TRÌNH VẬN TỐC GVHD: Ths Nguyễn Ngọc Quỳnh Phan Ngọc Khương Cát Lớp: L26 Nhóm số: MỤC LỤC TĨM TẮT BÁO CÁO I CHƯƠNG 1: MỞ ĐẦU 1.1 Yêu cầu 1.2 Điều kiện 1.3 Nhiệm vụ CHƯƠNG 2: CƠ SỞ LÝ THUYẾT 2.1 Chuyển động chất điểm không gian hai chiều 2.2 Vectơ vận tốc 2.3 Áp dụng CHƯƠNG 3: MATLAB 3.1 Tổng quan Matlab 3.2 Các hàm Matlab sử dụng toán 3.3 Giải toán Matlab 3.3.1 Đoạn code hoàn chỉnh 3.3.2 Giải thích đoạn code 3.4 Sơ đồ khối CHƯƠNG 4: KẾT QUẢ VÀ KẾT LUẬN 4.1 Kết trường hợp 4.1.1 Trường hợp 4.1.2 Trường hợp 10 4.2 Kết luận 10 TÀI LIỆU THAM KHẢO Error! Bookmark not defined DANH MỤC HÌNH ẢNH Hình 4.1 Hình chụp kết từ hình Command Window Matlab .9 Hình 4.2 Quĩ đạo chất điểm theo thời gian Hình 4.3 Kết chụp từ hình command window 10 Hình 4.4 Quĩ đạo chất điểm theo thời gian 10 I TÓM TẮT BÁO CÁO Bài báo cáo nhằm mục đích vận dụng sở lí thuyết học liên quan đến quỹ đạo, vận tốc, gia tốc vật, từ xác định vẽ chuyển động vật đồng thời tính tốn thơng số liên quan phần mềm hỗ trợ MATLAB Phần nội dung trình bày báo cáo viết CODE MATLAB 1 CHƯƠNG 1: MỞ ĐẦU 1.1 Yêu cầu - Vận tốc chất điểm mặt phẳng Oxy xác định phương trình: v  a cos  bt  i  cx j - Cho trước giá trị a,b,c xác định quỹ đạo vật vẽ quỹ đạo 1.2 Điều kiện - Sinh viên cần có kiến thức lập trình MATLAB - Tìm hiểu lệnh MATLAB liên quan symbolic đồ họa 1.3 Nhiệm vụ Xây dựng chương trình MATLAB: a) Nhập giá trị ban đầu (những đại lượng đề cho) b) Thiết lập phương trình tương ứng Sử dụng lệnh symbolic để giải hệ phương trình Từ đưa phương trình chuyển động vât kết luận quỹ đạo c) Vẽ hình quỹ đạo vật theo thời gian 2 CHƯƠNG 2: CƠ SỞ LÝ THUYẾT 2.1 Chuyển động chất điểm không gian hai chiều Gắn vào hệ trục tọa độ Descartes với hai trục Ox, Oy với hai vecto đơn vị, ta bán kính vectơ vẽ từ gốc hệ tọa độ Khi chất điểm chuyển động, vecto vị trí 𝑟⃗ chuyển động với thay đổi tọa độ hai phương trình x(t) y(t) theo thời gian t Để xác định quỹ đạo ta thực cách khử tham số t từ phương trình chuyển động, từ ta phương trình quỹ đạo 2.2 Vectơ vận tốc Vận tốc tức thời: 𝑑𝑟⃗ 𝑑𝑥 𝑑𝑦 = 𝑖⃗ + 𝑗⃗ = 𝑣𝑥 𝑖⃗ + 𝑣𝑦 𝑗⃗ 𝑑𝑡 𝑑𝑡 𝑑𝑡 Ta phương trình vận tốc theo hai phương x y: 𝑣⃗ = dx dy 𝑣𝑥 = dt 𝑣𝑦 = dt Vectơ vận tốc 𝑣⃗ đạo hàm vectơ vị trí theo thời gian, có gốc đặt điểm chuyển động, phương tiếp tuyến với quĩ đạo điểm đó, chiều chiều chuyển động có độ lớn v 2.3 Áp dụng Đề cho ta phương trình vận tốc theo hai phương x y, nhờ vào kiến thức học ta biết vectơ vận tốc 𝑣⃗ đạo hàm vectơ vị trí theo thời gian, ta suy nguyên hàm vectơ vận tốc 𝑣⃗ vectơ vị trí, từ ta xác định phương trình chuyển động vật khử bỏ tham số t để quĩ đạo vật theo thời gian t Đề bài: v  a cos  bt  i  cx j a dx 𝑣𝑥 = dt = acos(bt) → x(t) = ∫ 𝑣𝑥 dx = sin bt + C1 b ac dy 𝑣𝑦 = = cx = sin bt + cC1 dt b { a x(t) = sin bt + C1 b →{ 𝑎𝑐 𝑦(𝑡) = ∫ 𝑣𝑦 = − cos 𝑏𝑡 + 𝑐𝐶1𝑡 + 𝐶2 𝑏 Quĩ đạo chất điểm phụ thuộc vào điều kiện ban đầu chất điểm: * Trường hợp 1: Nếu C1 = 0, C2 = tương ứng với x(0) = 0, y(0) = ta được: a x(t) = sin bt b { 𝑎𝑐 𝑦(𝑡) = − cos 𝑏𝑡 𝑏 → x2 y2 + =1 a −ac ( ) ( 2) b b Khử tham số t ta quĩ đạo vật có dạng hình elip * Trường hợp 2: Nếu C1≠ C2 ≠ 0, ta quĩ đạo chất điểm đường cong không xác định 4 CHƯƠNG 3: MATLAB 3.1 Tổng quan Matlab - Matlab (viết tắt Matrix Laboratory) ngơn ngữ lập trình bậc cao bốn hệ, mơi trường để tính tốn số học, trực quan lập trình - Cơng cụ cho phép thao tác với ma trận, vẽ biểu đồ với hàm số liệu, thực thuật toán, tạo giao diện người dùng, phân tích liệu, phát triển thuật tốn, tạo kiểu mẫu ứng dụng - Phần mềm cơng cụ hỗ trợ đắc lực việc tính tốn, vẽ hình, vẽ biểu đồ thơng dụng thực thi phương pháp tính tốn 3.2 Các hàm Matlab sử dụng toán Lệnh Cú pháp Function function phuongtrinhquydao Ý nghĩa Tạo hàm mới, tên tập tin phuongtrinhquydao Clear all Clear all Xoá tất biến khỏi vùng làm việc Clc clc Xoá văn khỏi command window Close all close all Đóng chương trình cũ Syms syms x Khai báo biến x biến ký hiệu 5 Input x=input(‘nhập giá trị biến’) Hiện thị dấu nhắc lệnh chờ đầu vào Đạo hàm x theo t Diff diff(x,t) Cond cond = x(0) == C1 Điều kiện ban đầu x cond = x(0) == C2 y thời điểm ban đầu Dsolve dsolve(Vx(t),cond) Nguyên hàm hàm x với điều kiện ban đầu Disp disp([‘chuỗi ký tự’,char(x)]) Hiển thị nội dung mảng chuỗi Char Char(x) Chuyển đổi dạng hàm symbolic dạng chuỗi ký tự Fplot Fplot(x,y) Vẽ đồ thị đường cong theo hàm x y theo t phạm vi mặc định Title title(‘tên đồ thị’) Label xlabel(‘tên’) Thêm nhãn cho trục x yabel(‘tên’) Thêm nhãn cho trục y Grid on Hiển thị ô ly phần đồ Grid Tên đồ thị thị 6 3.3 Giải toán Matlab 3.3.1 Đoạn code hoàn chỉnh function phuongtrinhquydao Clear all; clc; close all; % khai bao bien syms x(t) y(t) a b c % nhap cac gia tri input a = input('nhap gia tri a = '); b = input('nhap gia tri b = '); c = input('nhap gia tri c = '); Vx(t) = input('nhap Vx(t) = '); C1 = input('nhap dieu kien tai thoi diem ban dau cua x(0) = '); C2 = input('nhap dieu kien tai thoi diem ban dau cua y(0) = '); % chay code Vx(t) = diff(x,t) == Vx(t) cond = x(0) == C1; x(t) = dsolve(Vx(t),cond); Vy(t) = input('nhap Vy(t) = '); Vy(t) = diff(y,t)== Vy(t); cond = y(0) == C2; y(t) = dsolve(Vy(t),cond); disp(['x(t) =',char(x)]) disp(['y(t) =',char(y)]) % Ve quy dao cua vat theo thoi gian t fplot(x,y); title('Phuong trinh quy dao cua vat theo thoi gian') xlabel('truc Ox'); ylabel('truc Oy'); grid on; end 3.3.2 Giải thích đoạn code Ta giải tốn theo bước: Bước 1:Tạo hàm với tên tập tin phuongtrinhquydao function phuongtrinhquydao Bước 2: Nhập giá trị a,b,c, C1 C2 thời điểm ban đầu a = input('nhap gia tri a = '); b = input('nhap gia tri b = '); c = input('nhap gia tri c = '); C1 = input('nhap dieu kien tai thoi diem ban dau cua x tai thoi diem ban dau = '); C2 = input(''nhap dieu kien tai thoi diem ban dau cua y tai thoi diem ban dau = '); Bước 3: Nhập phương trình vận tốc theo x Vx(t) = input('nhap Vx(t) = '); Bước 4:Nhập Vx(t) đạo hàm x(t) theo t, ta tính x(t) nguyên hàm Vx(t) với điều kiện ban đầu x(0) = C1 Vx(t) = diff(x,t) == Vx(t); cond = x(0) == C1; x(t) = dsolve(Vx(t),cond); Bước 5: Nhập Vx(t) đạo hàm x(t) theo t, ta tính y(t) nguyên hàm Vy(t) với điều kiện ban đầu y(0) = C2 Vy(t) = input('nhap Vy(t) = '); Vy(t) = diff(y,t)== Vy(t); cond = y(0) == C2; y(t) = dsolve(Vy(t),cond); Bước 6: Đưa kết cần tìm hình: disp(['x(t) =', char(x)]) disp(['y(t) =', char(y)]) Bước 7: Vẽ quỹ đạo vật đồ thị hình: fplot(x,y); title('Phuong trinh quy dao cua vat theo thoi gian') xlabel('truc Ox'); ylabel('truc Oy'); grid on; 3.4 Sơ đồ khối Nhập hệ số a, b, c, điều kiện ban đầu C1 C2 Nhập phương trình vận tốc Vx(t) Vy(t) Tìm phương trình chuyển động x(t) y(t) Vẽ quĩ đạo vật từ x(t) y(t) 9 CHƯƠNG 4: KẾT QUẢ VÀ KẾT LUẬN 4.1 Kết 4.1.1 Trường hợp 1: C1 = 0, C2 = Hình 4.1 Hình chụp kết từ hình Command Window Matlab Hình 4.1 Kết Ta nhập giá trị a = 1, b = 2, c = 3, phương trình Vx(t) = a*cos(b*t) phương trình Vy(t) = c*x Ta xác định phương trình x(t) y(t) vật, từ ta vẽ quỹ đạo vật Hình 4.2 Quĩ đạo chất điểm theo thời gian 10 4.1.2 Trường hợp 2: C1≠ C2 ≠ Hình 4.3 Kết chụp từ hình command window Hình 4.4 Quĩ đạo chất điểm theo thời gian 4.2 Kết luận Đề tài hỗ trợ xác định quỹ đạo vật giúp ta vẽ quỹ đạo vật biết phương trình vận tốc vật Với phương pháp sử dụng phần mềm Matlab giúp thuận tiện dễ dàng việc 11 giải tốn tương tự mà khơng thể giải tay giúp có nhìn trực quan TÀI LIỆU THAM KHẢO https://www.mathworks.com/help/matlab/ref/diff.html https://www.youtube.com/watch?v=EjskX-MSRIQ A L Garcia and C Penland, MATLAB Projects for Scientists and Engineers, Prentice Hall, Upper Saddle River, NJ, 1996 http://www.algarcia.org/fishbane/fishbane.html 4.Sách “Bài tập vật lý A1”- Đại học Bách kho TP HCM 12 ... - Vận tốc chất điểm mặt phẳng Oxy xác định phương trình: v  a cos  bt  i  cx j - Cho trước giá trị a,b,c xác định quỹ đạo vật vẽ quỹ đạo 1.2 Điều kiện - Sinh viên cần có kiến thức lập trình. .. ta phương trình quỹ đạo 2.2 Vectơ vận tốc Vận tốc tức thời: