Need help writing a simple pseudo code for the following script

dt = 0.1;
omega1 = 10;
domega_dt = 1/1000;
R = 0.11;
S = pi*R^2;
rho = 1.225;
m = 0.43;
t=0:dt:2;
u=zeros(length(t));
v=zeros(length(t));
x=zeros(length(t));
y=zeros(length(t));
u(1) = 25;
for i = 2:length(t)
omega = omega1 - domega_dt* t(i);
vspin = omega*R;
speed = sqrt(u(i-1)^2+v(i-1)^2);
CL = 1/(2+(speed/vspin));
CD = 0.55 + 1/(22.5 + 4.2*(speed/vspin)^2.5)^0.4;
Lift = 0.5*rho*S*speed^2*CL;
Drag = 0.5*rho*S*speed^2*CD;
u(i) = u(i-1) + (dt/m)*(-Drag * u(i-1) - Lift*v(i-1))/max(speed,1e-16);
v(i) = v(i-1) + (dt/m)*(-Drag * v(i-1) + Lift*u(i-1))/max(speed,1e-16);
x(i) = x(i-1) + dt*(u(i-1)+u(i))/2;
y(i) = y(i-1) + dt*(v(i-1)+v(i))/2;
end

2 Commenti

Why do you need pseudo code if you already have the code?
Consider going through MATLAB Onramp. It will give you the basics you need to understand this code, and you can complete it while you wait for an answer.
Start with the first line. Think out what it does, and then think bout what it does in context of what the overall code should be used for. After all, it is you who apparently needs to use this code. So you probably have some clue what it SHOULD do. But just a list of lines of undocumented code is often worth no more than something produced by a chimp pecking away at a typewriter.

Accedi per commentare.

Risposte (0)

Categorie

Scopri di più su Programming in Centro assistenza e File Exchange

Richiesto:

il 21 Mar 2021

Commentato:

il 21 Mar 2021

Community Treasure Hunt

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

Start Hunting!

Translated by