Numerical Method Terminal Velocity
59 views (last 30 days)
Show older comments
Parachutist of mass(m) = 68.1 kg
drag coefficient(c) = 12.5 kg/s
g= gravitational acceleration = 9.81
ti=0
vi=0
ti+1 - ti = 0.1
here is the eqn
v(ti+1)=v(ti)+(g-(c/m)*v(ti))*(ti+1 - ti)
I just want to plot within such a point that v(ti+1)-v(t)<0.001
thanks
3 Comments
Torsten
on 23 Feb 2017
If you set the initial condition for v to be zero, the solution is v=0 for all times. I guess this is not what you want.
Best wishes
Torsten.
Answers (2)
Jan
on 23 Feb 2017
Edited: Jan
on 23 Feb 2017
The terminal velocity is reached, when there is no further acceleration. This means that g-(c/m)*v(ti) must be 0.0 and you can calculate the result without any iterations or rough limits.
If you really want to calculate this by a loop:
v(1) = 0;
ti = 1;
tStep = 0.1;
vStep = inf; % Arbitrary large value to allow entering the loop
while vStep > 0.001
... increase ti by 1 (not by 0.1)
... calculate new speed and store it in v(ti) using tStep (not ti)
... calculate the step in the velocity vStep
end
0 Comments
See Also
Categories
Find more on Matrix Indexing 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!