Scribd
Upload
20 views6 pages

Determining The Trajectory of Projectile Motion in Gravity

Uploaded by

Diệp Sương
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
20 views6 pages

Determining The Trajectory of Projectile Motion in Gravity

Uploaded by

Diệp Sương
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 6

Determining the trajectory of projectile motion in gravity

with drag force


Group Information:
Group No: 3 Class: CC21 Completion day: December 12, 2022
1. Hoang Thi Diep Suong Student ID: 2153771
2. Pham Hong Thai Student ID: 2153793
3. Tran Minh Duy Student ID: 2153257
4. Nguyen Huy Tai Student ID: 2110513

1.Introduction
Motion, in physics, is defined as the change in position of an object with respect to its
surroundings in a given interval of time.

Projectile motion is the motion of an object projected into the air with an initial velocity v 0 with
horizontal angle α (throw angle) and moving under the influence of gravity.

Projectile motion with linear drag is the motion of an object thrown or projected into the air,
subject is affected by gravity and linear drag force (such as air resistance, wind, …). The object
is called a projectile, and its path is called its trajectory. The motion of falling objects, is a two-
dimensional type of projectile motion in which there are horizontal and vertical movement. This
project considers two-dimensional projectile motion, such as that of the cannon ball or other
object with air resistance (figure 1) [1].

Figure 1: Projectile motion with drag force [1]


2.Theory
Motion equation of projectile motion in gravity with drag force is represented as follows
m ⃗a =m⃗g −h ⃗v
- Projecting on x-axis (from the left hand): ma=−h v x (1)
'' 2 h '
o Applying the differential equation: (1)  x ( t ) d t + x ( t ) dt=0
m
with x 0=0 and x 'o=v 0 cosα (since v 0 x =v 0 cosα )
- Protecting on y-axis (from the bottom up): ma=mg−h v y (2)
'' h '
o Applying the differential equation: (2)  y ( t )+ g + y ( t )=0
m
with 0y =0 '
v =v
and y o =v 0 sinα (since 0 y 0 sinα )

Using the symbolic calculation of MATLAB, there is a mentioned equations and initial
conditions to calculate the trajectory, and plot the trajectory depending on the angle.

3. MATLAB Code and Explanation

Enter the values m, h, v0, , t (flight time)


clear; clc;
syms x(t) y(t);

m = input ('Enter its weight (kg): ');


v0 = input ('Enter initial velocity (m/s): ');
a1 = input ('Enter the angle (degree): ');
a = a1/180 * pi;
h = input ('Enter the coefficient of drag force (kg/s): ');
t0 = input ('Enter the fight time (s): ');
g = 9.8;

Set up the differential equations with x (t) and y (t). Use symbolic commands to solve the system
of equations.
Dx = diff(x,t);
Dy = diff(y,t);
ode1 = diff(x,2) == -(h/m)*Dx;
ode2 = diff(y,2) == -g -(h/m)*Dy;
ic = [Dx(0) == v0*cos(a); Dy(0) == v0*sin(a); x(0) == 0; y(0) ==
0];

h = dsolve(ode1,ode2,ic);
disp ('The motion equation of particular is: ');
disp ('x= '); disp(h.x)
disp ('y= '); disp(h.y)

Plot the trajectory depending on the angle  (15, 30, 45, 60, 750).
fplot(h.x,h.y,[0,t0]);
xlabel('x');
ylabel('y');

4. Results and discussion


4.1 Plot the trajectory depending on the angle α =15 °

Figure 4.1: Trajectory with the angle α =15 ° [2]

4.2 Plot the trajectory depending on the angle α =30 °


Figure 4.2: Trajectory with the angle α =30 ° [3]

4.3 Plot the trajectory depending on the angle α =45°

Figure 4.3: Trajectory with the angle α =45° [4]


4.4 Plot the trajectory depending on the angle α =60 °

Figure 4.4: Trajectory with the angle α =60 ° [5]

4.5 Plot the trajectory depending on the angle α =75 °


Figure 4.5: Trajectory with the angle α =75 ° [6]

Discussion:
Above results exactly match the results manually calculated. With MATLAB calculation, we can
replace appropriately many others values of quantities to study other special cases.

5. Conclusion
The project has completed the solution of the projectile motion problem using MATLAB
symbolic calculation. With this tool we can solve more complex motion situations that cannot be
solved by the analytical method.

References:
[1] PhD. Paola La Rocca & Prof. Francesco Riggi (Professor), Projectile motion with a drag
force: Were the Midievals right after all?, 2009. DOI: 10.1088/0031-9120/44/4/009
https://www.researchgate.net/publication/
228919307_Projectile_motion_with_a_drag_force_Were_the_Medievals_right_after_all
[2] OER Service, Projectile motion
https://courses.lumenlearning.com/suny-physics/chapter/3-4-projectile-motion/

You might also like