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!