Reto Motor de Simulación del Movimiento
Programa principal:
for i=1:3
Xo = 0;
Yo = 1000;
dt = 0.5;
cd = 0.5;
rho_aire = 1.2;
rho_rock = 4500;
gravedad = -9.81;
Vo = rand()*100;
angulo = rand()*180;
radio = rand()*10;
are = pi*radio^2;
volumen = 4/3*pi*radio^3;
masa = volumen*rho_rock;
b = 0.5*are*cd*rho_aire;
vox = Vo*cosd(angulo);
voy = Vo*sind(angulo);
ax = -b*vox^2/masa;
ay = gravedad-(voy/abs(voy)*b*voy^2/masa);
x_1 = Xo-vox*dt;
y_1 = Yo-voy*dt+0.5*gravedad*dt^2;
eje_x = -320:1:320;
eje_y = -0.01*eje_x.^2 + Yo;
area(eje_x, eje_y, "FaceColor", "[0, 0.4470, 0.7410]") ;
y = 0;
lista_x = (Xo);
lista_y = (Yo);
while y >= 0
x = 2*Xo-x_1+0.5*ax*dt^2; %ECUACIONES EULER
y = 2*Yo-y_1+0.5*ay*dt^2; %ECUACIONES EULER
lista_x(end+1) = x
lista_y(end+1) = y
x_1 = Xo;
y_1 = Yo;
Xo = x;
Yo = y;
vox = (Xo-x_1)/dt;
voy = (Yo-y_1)/dt;
ax = -b*vox^2/masa;
ay = gravedad- (voy/abs(voy)*b*voy^2/masa) ;
end
lista_x = lista_x(1:end-1);
lista_y = lista_y(1:end-1);
hold on
title('Simulacion del volcan')
comet(lista_x,lista_y)
end
1
