movements of the descent vehicle in the atmosphere of Earth
5 次查看(过去 30 天)
显示 更早的评论
I need help writting a Matlab script file, In the matlab I need to write the equations of motion of the space capsule that are as follows:
and I need the graph of time (x axis) vs speed (y axis).
3 个评论
Ameer Hamza
2020-5-22
编辑:Ameer Hamza
2020-5-22
Is it necessary to write your own RK-4 code, or can you use ode45?
采纳的回答
James Tursa
2020-5-22
You have four differential equations involving four variables: V, theta, H, and L. You have to solve these simultaneously since they depend on each other. You can't just pick the one you want and try to solve it by itself like this:
f = @(t,v) v*sind(theta);
You need an f that takes in a 4-element state vector containing V, theta, H, and L, and produces a 4-element derivative vector containing dV/dt, dtheta/dt, dH/dt, and dL/dt. Then that is used in your RK4 looping code.
15 个评论
James Tursa
2020-5-26
Please show us your current code. You should have two separate files. One script file that sets constants and executes the RK4 code, and another derivative function file.
Also, I did not notice the large stepsize z that you have. That simply won't work. You need the t vector to have a much smaller stepsize. E.g.,
z = 0.5;
t = 0:z:900;
And H should not be a vector that you initialize. H should just be a single scalar number that is the initial value. Your RK4 code will be calculating the other H values dynamically as part of the simulation. E.g.,
H = 100000; % initial height (m)
James Tursa
2020-5-27
You are going to get garbage out if you continue to use a huge stepsize. You need to drop the stepsize as I suggested. E.g.,
z = 0.5;
t = 0:z:900;
更多回答(0 个)
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!