báo cáo bài tập lớn vật lý 1 đề tài xác định quỹ đạo của vật

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA TP HCM - BÁO CÁO BÀI TẬP LỚN VẬT LÝ Đề Tài 21 Xác Định Qũy Đạo Của Vật Giảng viên tập : Nguyễn Ngọc Thiện Giảng viên lý thuyết:Phạm Thị Hải Miền LỚP BÀI TẬP: L29-A TÊN NHĨM: Nhóm 02 21/4/2021 Sinh viên thực hiện: Họ Và Tên Mã Số Sinh Viên Lê Phương Các 2110833 Phạm Trường Giang 2111105 Phạm Nguyễn Thành Danh 2110072 Hồ Lệ Giang 2111093 Lê Khánh Duy 2113009 LỚP BÀI TẬP: L29-A TÊN NHĨM: Nhóm 02 MỤC LỤC GIỚI THIỆU ĐỀ BÀI YÊU CẦU ĐỀ BÀI 2.1 INPUT .4 2.2 OUTPUT CƠ SỞ LÝ THUYẾT: 4 ĐOẠN CODE MATLAB VÀ DIỄN GIẢI CHI TIẾT 5 KẾT QUẢ .7 GIỚI THIỆU ĐỀ BÀI Với đề tài xác định quỹ đạo vật mặt phẳng nhóm chúng em tìm hiểu chuyển động vật để đưa sở lí thuyết tốn quỹ đạo dùng cơng thức học để đưa thông số chuyển động, thiết lập phương trình chuyển động sử dụng phần mềm Matlab để giải số toán đo kết hình vẽ chuyển động vật YÊU CẦU ĐỀ BÀI Sử dụng Matlab để giải toán sau: 2.1 INPUT “Vận tốc chất điểm chuyển động mặt phẳng Oxy xác định biểu thức ⃗v =acos ( bt ) ⃗i + cx ⃗j Cho trước giá trị a, b c, 2.2 OUTPUT Xác định quỹ đạo vật vẽ quỹ đạo đó?” CƠ SỞ LÝ THUYẾT: Chọn gốc tọa độ O, trục Ox hướng sang phải, trục Oy hướng lên Ta có phương trình vận tốc ban đầu vật: ⃗v =acos ( bt ) ⃗i + cx ⃗j Với a,b,c giá trị cho trước Chiếu vật lên trục Ox, Oy: cos ( bt ) {Vx=aVy =cx Để tìm phương trình chuyển động theo phương y, ta cần tìm phương trình chuyển động theo phương x: Ta có: x = x o+∫ v x dt ⇒ x o+ a sin ⁡(bt) b Ta tìm phương trình chuyển động theo phương y y = y o + ∫ v y dt ⇒ y = y o + c x ot - ac cos ⁡(bt ) b2 Phương trình chuyển động vật : { a x=x o + sin ( bt ) b ac y= y o +c x o t− cos ( bt ) b Vậy với vị trí ban đầu khác vật có quỹ đạo khác nhau, nên quỹ đạo vật phụ thuộc vào vị trí ban đầu ĐOẠN CODE MATLAB VÀ DIỄN GIẢI CHI TIẾT %Khai bao cac bien syms x y Vx Vy t n %Nhap he so a b c a=input('Nhap a: '); b=input('Nhap b: '); c=input("Nhap c: "); %Tinh van toc theo phuong x Vx=a*cos(b*t); disp('Phuong trinh Vx la: ') disp(Vx) %Chon cac truong hop co the disp('Hay lua chon vi tri ban dau cua chat diem: ') disp('1: Cx=Cy=0') disp('2: Cx=0 & Cy!=0') disp('3: Cx!=0 & Cy=0') disp('4: Cx!=0 & Cy!=0') n=input('Nhap n: '); switch n case %Tinh toan phuong trinh chuyen dong theo phuong x va y disp('Phuong trinh x la: ') x=int(Vx,t); disp(x) Vy=c*x; disp('Phuong trinh Vy la: ') disp(Vy) disp('Phuong trinh y la: ') %Neu b=2 thi nguyen ham sin(2*x) cho ket qua: sin(x)^2 dung lenh %intergal thi khac biet ve hang so C khong dang ke Do phai xet truong %hop rieng cho b=2 %Source: https://www.mathworks.com/matlabcentral/answers/419375-why-doi-get-sin-t-2-when-i-int-sin-2-t if b==2 y = (-c*a)/(b*b)*cos(b*t); else y=int(Vy,t); end disp(y) case %Tinh toan phuong trinh chuyen dong theo phuong x va y disp('Phuong trinh x la: ') x=int(Vx,t); disp(x) Vy=c*x; disp('Phuong trinh Vy la: ') disp(Vy) syms Cy Cy=input('Nhap vi tri ban dau cua chat diem theo phuong y: '); disp('Phuong trinh y la: ') %Neu b=2 thi nguyen ham sin(2*x) cho ket qua: sin(x)^2 dung lenh %intergal thi khac biet ve hang so C khong dang ke Do phai xet truong %hop rieng cho b=2 %Source: https://www.mathworks.com/matlabcentral/answers/419375-why-doi-get-sin-t-2-when-i-int-sin-2-t if b==2 y = (-c*a)/(b*b)*cos(b*t) +Cy; else y=int(Vy,t)+Cy; end disp(y) case %Tinh toan phuong trinh chuyen dong theo phuong x va y syms Cx Cx=input('Nhap vi tri ban dau cua chat diem theo phuong x: '); disp('Phuong trinh x la: ') x=int(Vx,t)+Cx; disp(x) Vy=c*x; disp('Phuong trinh Vy la: ') disp(Vy) disp('Phuong trinh y la: ') %Neu b=2 thi nguyen ham sin(2*x) cho ket qua: sin(x)^2 dung lenh %intergal thi khac biet ve hang so C khong dang ke Do phai xet truong %hop rieng cho b=2 %Source: https://www.mathworks.com/matlabcentral/answers/419375-why-doi-get-sin-t-2-when-i-int-sin-2-t if b==2 y = (-c*a)/(b*b)*cos(b*t)+c*Cx*t; else y=int(Vy,t); end disp(y) case %Tinh toan phuong trinh chuyen dong theo phuong x va y syms Cx Cx=input('Nhap vi tri ban dau cua chat diem theo phuong x: '); disp('Phuong trinh x la: ') x=int(Vx,t)+Cx; disp(x) Vy=c*x; disp('Phuong trinh Vy la: ') disp(Vy) syms Cy Cy=input('Nhap vi tri ban dau cua chat diem theo phuong y: '); disp('Phuong trinh y la: ') %Neu b=2 thi nguyen ham sin(2*x) cho ket qua: sin(x)^2 dung lenh %intergal thi khac biet ve hang so C khong dang ke Do phai xet truong %hop rieng cho b=2 %Source: https://www.mathworks.com/matlabcentral/answers/419375-why-doi-get-sin-t-2-when-i-int-sin-2-t if b==2 y = (-c*a)/(b*b)*cos(b*t)+c*Cx*t+Cy; else y=int(Vy,t)+Cy; end disp(y) end %Ve thi fplot(x,y) grid on title('Quy Dao Chuyen Dong') xlabel('x(t)') ylabel('y(t)') KẾT QUẢ Chạy code trường trình với hệ số a=1 b=2 c=3 //Nhập hệ số //Trường hợp 1: Các vị trí ban đầu //Trường hợp 2: Vị trí ban đầu theo phương y khác Xét Cy=5 //Trường hợp 3: Vị trí ban đầu theo phương x khác Xét Cx=6 //Trường hợp 4: Vị trí ban đầu theo phương x y khác Xét Cx=6 Cy=5 10 11