Scribd
Upload
10 views6 pages

General Physic 1 - Matlab Project - Group 5

This report details a MATLAB project focused on determining the trajectory of projectile motion under the influence of gravity and drag force. It includes theoretical background, the formulation of differential equations, MATLAB code for solving these equations, and results showing the trajectory for various launch angles. The project demonstrates the capability of MATLAB to analyze complex motion scenarios that are challenging to solve analytically.

Uploaded by

khoa.ngo2005
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
10 views6 pages

General Physic 1 - Matlab Project - Group 5

This report details a MATLAB project focused on determining the trajectory of projectile motion under the influence of gravity and drag force. It includes theoretical background, the formulation of differential equations, MATLAB code for solving these equations, and results showing the trajectory for various launch angles. The project demonstrates the capability of MATLAB to analyze complex motion scenarios that are challenging to solve analytically.

Uploaded by

khoa.ngo2005
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

Matlab project 3 [Date]

Ho Chi Minh City University of Technology


🙞···☼···🙜

REPORT
GENERAL PHYSIC 1
……………………………………………………………………………………
…………………………………………………………………………………….

Lecturer: …Nguyễn Xuân Thanh Trâm ……….


Class : CC12
Group 5

Student Student code


Ngô Minh Khoa 2353567
Lê Tấn Minh Khoa 2352563
Nguyễn Hữu Nguyên Khoa 2352574
Huỳnh Minh Khôi 2352608
Nguyễn Hữu Minh Khôi 2352614
Nguyễn Lê Chí Kiên 2352643
Matlab project 3 [Date]

Project 3:

Determining the trajectory of projectile motion in gravity

with drag force

1.Introduction.
What Is Projectile Motion?
Projectile motion is the motion of an object thrown (projected) into the air when,
after the initial force that launches the object, air resistance is negligible and the
only other force that object experiences is the force of gravity. The object is
called a projectile, and its path is called its trajectory. Air resistance is a
frictional force that slows its motion and can significantly alter the trajectory of
the motion. Due to the difficulty in calculation, only situations in which the
deviation from projectile motion is negligible and air resistance can be ignored
are considered in introductory physics.

projectile motion projectile motion with drag force

2. Theory
The motion equation of projectile motion can be described as follows:

With is gravity acceleration vector and h is the coefficient of drag. Expressing


by component in the Cartesian coordinates reference with Ox – the ground
surface and Oy – the height, we can obtain following differential equations:
Matlab project 3 [Date]

And:

With the initial conditions (at t = 0) as follows:

Using the symbolic calculation of MATLAB, we can solve mentioned


differential equations to obtain the motion equation of the projectile motion,
from which the necessary features of the motion can be derived.

3. MATLAB Code and Explanation


Formulate the differential equations of the motion including the initial
conditions:
% Define symbolic variables
syms x(t) y(t)
syms m v0 a g t h

% Define the differential equations


vx = diff(x,t);
vy = diff(y,t);
eqns_x = diff(x,2) == -h/m*vx;
eqns_y = diff(y,2) == -g - h/m*vy;

Solve the differential equations by the command dsolve of symbolic calculation:


% Solve the system of equations
cond = [x(0) == 0, y(0) == 0, vx(0) == v0*cos(a), vy(0) ==
v0*sin(a)];
sol = dsolve(eqns_x, eqns_y, cond);

Represent the results in the Command Window:


% Motion equations
disp('Motion equation for x:'); disp(sol.x);
disp('Motion equation for y:'); disp(sol.y);
Matlab project 3 [Date]

Enter the values:


% Define the values
g = 9.81;
m = input('Enter the mass m[kg]: ');
h = input('Enter the coefficient of drag: ');
v0 = input('Enter the initial velocity v0[m/s]: ');

Plot the trajectory for the launch angle a = 15°, 30°, 45°, 60°, 75° in the figure:
% Plot the trajectory for different angles
figure
hold on
for a = [15 30 45 60 75]*pi/180
% Calculate the trajectory
x = (m*v0*cos(a))/h - (m*v0*exp(-(h*t)/m)*cos(a))/h;
y = (m*(g*m + h*v0*sin(a)))/h^2 - (m*(g*t + (exp(-(h*t)/m)*(g*m
+ h*v0*sin(a)))/h))/h;
plot(x, y)
end
hold off
xlabel('Distance (m)')
ylabel('Height (m)')
legend('15°', '30°', '45°', '60°', '75°')
title('Projectile motion with drag for different launch angles')

4. Results and discussion


We assume that
*m=1 kg ,
*Vo= 1 m/s,
*h= 0.09 m
Matlab project 3 [Date]

-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
As we can see , the project has deciphered the projectile motion in gravity with
drag force by using the MATLAB symbolic calculation. With this tool we can
solve more complicated motion situations that are beyond analytic methode .

6. References:
University Physics with Modern Physics (14th Edition) 14th Edition .
Fundamentals of Physics by Resnick, Halliday, Walker
Matlab project 3 [Date]

You might also like