¿Matlab y el descenso un paracaidista?
Publicado por Simona (15 intervenciones) el 07/08/2013 09:08:42
Hola es que tengo un error, debo graficar la caída de un paracaidista con respecto al tiempo me da una velocidad constante y la masa
debo grafica velocidad con respecto al tiempo después de abierto el paracaidas
No me grafica de esta manera :S
function [detect,stopint,direction]=detener(t,x)
% Determina el instante cuando pasa por la altura cero y detiene la integración
detect=x(1); % Detecta cuando llega al suelo
stopint=1; % Detiene la integración
direction=0; % No tiene interés
end
m=78;
A=input('área del paraciaidas, A: ');
%condiciones iniciales
t0=zeros(1,2);
t0(1)=input('t inicial (s), t0: ');
t0(2)=0; %parte del reposo
k=1.29*A*0.8/2;
lambda=7482.2;
f=@(t,x)[x(2); -9.8+(k*x(2)*x(2)/m)*exp(-x(1)/lambda)];
tspan=[0 inf]; %inf es infinito positivo
opts=odeset('events',@detener);
[t,x,te,xe]=ode45(f,tspan,t0,opts);
%plot(x0(1)-x(:,1),-x(:,2),'r')
plot(t0(1)-t(:,1),-t(:,2),'r')
text(t0(1),-xe(2),num2str(-xe(2)))
grid on
xlabel('x')
ylabel('v');
title('caída de un paracaidista')
Gracias por la ayuda :)
debo grafica velocidad con respecto al tiempo después de abierto el paracaidas
No me grafica de esta manera :S
function [detect,stopint,direction]=detener(t,x)
% Determina el instante cuando pasa por la altura cero y detiene la integración
detect=x(1); % Detecta cuando llega al suelo
stopint=1; % Detiene la integración
direction=0; % No tiene interés
end
m=78;
A=input('área del paraciaidas, A: ');
%condiciones iniciales
t0=zeros(1,2);
t0(1)=input('t inicial (s), t0: ');
t0(2)=0; %parte del reposo
k=1.29*A*0.8/2;
lambda=7482.2;
f=@(t,x)[x(2); -9.8+(k*x(2)*x(2)/m)*exp(-x(1)/lambda)];
tspan=[0 inf]; %inf es infinito positivo
opts=odeset('events',@detener);
[t,x,te,xe]=ode45(f,tspan,t0,opts);
%plot(x0(1)-x(:,1),-x(:,2),'r')
plot(t0(1)-t(:,1),-t(:,2),'r')
text(t0(1),-xe(2),num2str(-xe(2)))
grid on
xlabel('x')
ylabel('v');
title('caída de un paracaidista')
Gracias por la ayuda :)
Valora esta pregunta


0