I have problem to get the Qb. Here is my code:
a=0.993; for j=1:length(Q) B=0.03; if (0<=B && B<=1) Qb(j)=((1-B)*a*Qb(j-1)+(1-a)*B*Q(j))/(1-a*B); elseif (B == 0) Qb(j)= a*Qb(j-1); else Qb(j)= Q(j)-Qs(j); end
end
Please help me to figure out.
Thanks!
Helen
Best Answer