i want to write a code for lumped vortex model for unsteady flow,

for t=1:10;
hm=0.5*sin((t));
alp=(0.5)*cos(t);
a=[1/(2*pi*(0.75-0.25)); 1];
b=[1/(2*pi*(0.75-10)) ; 1];
A=[a b];
B=[sin(alp)+M;0];
X=linsolve(A,B);
end
At t=1, M=0
For t>1, M depends on the values of X(2) at previous timesteps,
for example:
at t=2 M=(X(2) at t=1)/(2*pi*(1+1))
at t=3 M=(X(2) at t=2)/(2*pi*(1+1)) + (X(2) at t=1)/(2*pi*(1+2))
at t=4 M=(X(2) at t=3)/(2*pi*(1+1)) + (X(2) at t=2)/(2*pi*(1+2))+
(X(2) at t=1)/(2*pi*(1+3))
where X(2) is the 2nd element of the solutions from X=linsolve(A,B);
how can i do that by using a loop

Answers (0)

Categories

Asked:

on 3 Nov 2013

Edited:

on 3 Nov 2013

Community Treasure Hunt

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

Start Hunting!