Scribd
Upload
5 views2 pages

Codebtl

The document outlines a simulation of projectile motion, where initial velocity and launch angle are input to calculate the trajectory of an object under gravity. It includes equations for horizontal and vertical motion, as well as calculations for the radius of curvature at specific time points. The results are visualized in a plot and specific values for the radius of curvature at landing and at a given time are computed and displayed.

Uploaded by

36-Tuyết Trinh
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
5 views2 pages

Codebtl

The document outlines a simulation of projectile motion, where initial velocity and launch angle are input to calculate the trajectory of an object under gravity. It includes equations for horizontal and vertical motion, as well as calculations for the radius of curvature at specific time points. The results are visualized in a plot and specific values for the radius of curvature at landing and at a given time are computed and displayed.

Uploaded by

36-Tuyết Trinh
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd

%%nhap cac gia tri dau vao

v0=input(‘nhap gia tri v0=’)


nhap gia tri v0=15
alpha=input(‘nhap gia tri goc alpha=’)
nhap gia tri goc alpha=pi/6
alpha=0.5236
g=9.8
sysms t %%bien thoi gian t
%%thiet lap phuong trinh chuyen dong
x=v0*cos(alpha)*t
y=v0*sin(alpha)*t-1/2*g*t^2
%%tinh ban kinh quy dao
vx=diff(x,t)
vy=diff(y,t)
v=sqrt(vx^2+vy^2)
a=g
a=9.8
at=diff(v,t)
an=sqrt(a^2-at^2)
R=v^2/an

%%ve do thi
t_end=10
fplot(x,y,[0 t_end])
axis equal
axis ([0 25 0 5])
xlabel(‘Ox’)
ylabel(‘Oy’)
title(‘quy dao nem xien cua vat’)
%%thu tinh cac gia tri ban kinh quy dao tai
thoi diem bat ki
t0=0
%%ban kinh quy dao tai vi tri cham dat
Rt0=subs(R,t,t0)
fprintf(‘ban kinh quy dao tai vi tri cham=’)
disp(Rt0)
Rt0=(750*7203^(1/2))/2401
t1=2
Rt1=subs(R,t,t1)
Fprintf(‘ban kinh quy dao tai thoi diem t1=’)
Disp(Rt1)
Rt1=(7879*31516^(1/2)*1620675^(1/2))/40516875

You might also like