2 views (last 30 days)

Show older comments

clc;

Cd=xlsread('trial',1,'b2:b6');

rho = xlsread('trial',1,'a1:a6');

S=3.2365;

m=1000;

gama=10;

H=input('enter altitude in meter (m)');

V=input('enter the velocity (m/s)');

dt=1;

t=0:dt:10^3;

n=(10^3/dt+1);

Ux=zeros(1,n);

Ux = V*cos(gama);

Uy = V*sin(gama);

for i=1:n

D=0.5.*rho*V^2*S.*Cd;

Ux(1,i+1)=Ux(1,i)-(dt)*(D/m*cos(gama));

Uy(1,i+1)=Uy(1,i)+(dt)*(g-(D/m*sin(gama)));

V(1,i+1)=((Ux(1,i+1))^2+(Uy(1,i+1)^2)^0.5;

i=i+1;

end

James Tursa
on 6 Nov 2019

Lasse Jakobsen
on 6 Nov 2019

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!