Is this formula correct for simulate 3D Water rocket

Question

Le Long 2021-6-28

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/866730-is-this-formula-correct-for-simulate-3d-water-rocket

评论： Alan Stevens 2021-6-29

在 MATLAB Online 中打开

I want to simulate 3D Water Rocket only use Acceleration data collected from Sensor GY86

But after I run

The results gave me have a lot of line and isn't correct

I don't know if i use this Calculation is correct or enough or i am missing some steps in formula

I have tried this

function Water_rocket_simulator
clc
close all
clear all
load('Acce and Gyro xyz')
%% CONSTANTS
g = 9.81;
%% INPUT DATA
x = 0;
y = 0;
z = 0;
v0 = 30;
alpha = 90;
t = 0;
dt = 0.01;
%% FIGURE
figure('name','Water Rocket','color','white','numbertitle','off');
hold on
fig_rocket = plot3(x,y,z,'ro','MarkerSize',10,'markerfacecolor','r');
ht = title(sprintf('t = %0.2f s',t));
xlabel('X axis Latitude');
ylabel('Y axis Longtitude');
zlabel('Z axis heights');
grid on
axis equal
axis([1 300 1 300 1 300]);
%% CALCULATION
alpha = alpha/180*pi;
vx = v0*cos(alpha);
vy = v0*cos(alpha);
vz = v0*sin(alpha);
while z>-0.01
t = t+dt;
ax = AccelerationX(:,1)/1000;
ay = AccelerationY(:,1)/1000;
az = AccelerationZ(:,1)/1000;
vx = vx+ax*dt;
vy = vy+ay*dt;
vz = vz+az*dt;
x = x+vx*dt+0.5*ax*dt.^2;
y = y+vy*dt+0.5*ay*dt.^2;
z = z+vz*dt+0.5*az*dt.^2;
plot3(x,y,z,'o','markersize',1,'color','k');
xlabel('X axis Latitude');
ylabel('Y axis Longtitude');
zlabel('Z axis Heights');
grid on
set(fig_rocket,'xdata',x,'ydata',y,'zdata',z);
set(ht,'string',sprintf('t = %0.2f s',t));
pause(0.002);
end
end

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Alan Stevens 2021-6-29

在 MATLAB Online 中打开

Perhaps you need something like

for i = 2:numel(AccelerationX)
    t(i) = t(i-1) + dt;
    ax(i) = Acceleration(i,1)/1000;
    % similarly for ay and az
    vx(i) = vx(i-1)+ax(i)*dt;
   % ...etc
   x(i) = x(i-1) + vx(i)*dt +0.5*ax(i)*dt^2;
   
  % ...etc     
end