use of ode45 for projectile trajectory including drag
3 views (last 30 days)
Show older comments
equation of motion
dvx/dt = -k*v*vx
dvy/dt = -k*v*vy - g
where k = drag coefficient
v = sqrt((vx)^2 + (vy)^2)
after changing it into space variable it becomes
dvx/dx = -k*v,
dvy/dx = (-k*v*vy - g)/vx,
dt/dx = 1/vx,
dy/dx = vy/vx,
form these equation i have to find for fixed range :- time taken ,vx,vy,y-displacement at particular x-displacement using odesolver if i use ode45 it gives error because of only two variable in right side of each equation (vx and vy). how wil i solve this equation using ode45 or any other ode solver
please reply soon
thanks
0 Comments
Accepted Answer
Grzegorz Knor
on 5 Sep 2011
To solve this system, create a function ode_func containing the equations:
function dy = ode_func(x,y)
k = 1;
g = 9.8;
v = sqrt(y(1)^2+y(2)^2);
dy = zeros(4,1);
dy(1) = -k*v;
dy(2) = (-k*v*y(2)-g)/y(1);
dy(3) = 1/y(1);
dy(4) = y(2)/y(1);
And solve with ode45:
[X,Y] = ode45(@ode_func,...
0 Comments
More Answers (0)
See Also
Categories
Find more on Ordinary Differential Equations in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!