% rocket - Compute motion of a rock

% rocket - Compute motion of a rocket using the Verlet method
clear all; % Clear memory
help rocket; % Print header

%@ Initialize variables (e.g., mass, burn rate, exhaust velocity)
y = 0; % Initial position (m)
mass = 2.5e6; % Initial mass of the rocket (kg)
dm_dt = -1.6e4; % Burn rate (kg/s)
gravity = 9.81; % Gravitational acceleration (m/s^2)
ExhaustU = input('Enter exhaust velocity (m/s): ');
FinalT = input('Enter final time (s): ');
dt = 0.5; % Time increment (s)
NSteps = FinalT/dt; % Number of steps

%@ Do the first step using the Euler method
a = -dm_dt*ExhaustU/mass - gravity; % Acceleration
yNew = y + dt*a; % Compute the new height (Euler)
yOld = y; % Current position becomes old position
y = yNew; % New position becomes current position
mass = mass + dt*dm_dt; % Compute new mass
tPlot(1) = dt; % Record time
yPlot(1) = y; % Record height

%@ Do the rest of the steps using the Verlet method
for iStep=2:NSteps

%@ Compute the new height using Verlet
a = -dm_dt*ExhaustU/mass - gravity; % Acceleration
yNew = 2*y - yOld + dt^2*a;
yOld = y; % Current position becomes old position
y = yNew; % New position becomes current position

%@ Compute the new mass
mass = mass + dt*dm_dt;

%@ Record time and height for plotting
tPlot(iStep) = iStep*dt; % Record time
yPlot(iStep) = y; % Record height

end

%@ Plot the position versus time
clf; figure(gcf); % Clear figure; bring figure window forward
plot(tPlot,yPlot); % Plot height versus time
xlabel('Time (s)'); ylabel('Height (m)');

Kết quả (Việt) 1:

tên lửa % - tính toán chuyển động của một tên lửa sử dụng phương pháp Verletrõ ràng tất cả; Bộ nhớ rõ ràng %giúp tên lửa; % In tiêu đề% @ Khởi tạo biến (ví dụ, khối lượng, ghi tốc độ, thải vận tốc)y = 0; Vị trí ban đầu % (m)khối lượng = 2.5e6; % Khối lượng ban đầu của tên lửa (kg)dm_dt = - 1.6e4; Tốc độ ghi % (kg/s)lực hấp dẫn = 9,81; % Tăng tốc hấp dẫn (m/s ^ 2)ExhaustU = đầu vào ('Enter ống xả tốc độ (m/s): ');FinalT = đầu vào (' nhập cuối cùng thời gian (s): ');DT = 0,5; % Thời gian tăng (s)NSteps = FinalT/dt; % Số bước% @ Làm bước đầu tiên bằng cách sử dụng phương pháp Euler a =-dm_dt * ExhaustU/trọng lượng - trọng lực; % Tăng tốcyNew = y + dt * a; % Tính toán tầm cao mới (Euler)yOld = y; % Vị trí hiện tại sẽ trở thành vị trí cũy = yNew; % Vị trí mới trở thành vị trí hiện tạikhối lượng = khối lượng + dt * dm_dt; % Tính mới hàng loạttPlot(1) = dt; Thời gian kỷ lục %yPlot(1) = y; Chiều cao kỷ lục %% @ Làm phần còn lại của các bước bằng cách sử dụng phương pháp Verletcho iStep = 2:NSteps % @ Tính toán tầm cao mới bằng cách sử dụng Verlet a =-dm_dt * ExhaustU/trọng lượng - trọng lực; % Tăng tốc yNew = 2 * y - yOld + dt ^ 2 * a; yOld = y; % Vị trí hiện tại sẽ trở thành vị trí cũ y = yNew; % Vị trí mới trở thành vị trí hiện tại % @ Tính toán khối lượng mới khối lượng = khối lượng + dt * dm_dt; % @ Thời gian kỷ lục và chiều cao cho âm mưu tPlot(iStep) = iStep * dt; Thời gian kỷ lục % yPlot(iStep) = y; Chiều cao kỷ lục % kết thúc% @ Âm mưu vị so với thời gian CLF; figure(GCF); % Con số rõ ràng; mang hình cửa sổ về phía trướcPlot(tPlot,yPlot); Chiều cao âm mưu % so với thời gianxlabel ('thời gian (s)'); ylabel ('chiều cao (m)');