what is wrong with this while loop?

2018 9 6
1 个回答
更新时间：2018 9 6
6 次查看（30 天）

0 个投票

function [Gj,Lj,Xij,Yij,Tj]=SR(N,Tj,Fl,Fg,Zij,Qj,TF,Pj)
%inicio
Fj=Fl+Fg;
[Gj,Lj]=MetodoSR(N,Fg,Fl);
Tj0=0;
k=0;
while (((Tj-Tj0)/Tj)^2)>=1e-10 && k<=100
k=k+1;
for k=1:N
Tj0=Tj(N-1);
Kij=ctteK1(Pj,Tj,N);
%Conversion de C a F
Xij=MatrizxijT(N,Fg,Fl,Gj,Kij,Zij);
[Lj,Gj,Xij,Yij]=Flujos(N,Lj,Xij,Fg,Fl,Kij);
Tj=Tj*1.8+32;
TF=Tj*1.8+32;
%Conversion de kmol a lbmol
[MatrizTH]=TH(N,Tj,Xij,Yij,Lj,Gj);
[HFlj,HFgj]=THFeed(N,TF,Zij);
[Tj]=Balance(N,Tj,Yij,Xij,Lj,Gj,Fj,HFlj,HFgj,Qj,MatrizTH);
end
end

2 个评论
显示 无 隐藏 无

madhan ravi
madhan ravi 2018-9-6
编辑：madhan ravi 2018-9-6
provide all the datas for (N,Tj,Fl,Fg,Zij,Qj,TF,Pj) and all the functions and paste it here
Giulliana Marin
Giulliana Marin 2018-9-6
There are way too many functions. The code used to run normally until I added the loop, meaning I added this:
Tj0=0;
k=0;
while (((Tj-Tj0)/Tj)^2)>=1e-10 && k<=100
k=k+1;
for k=1:N
Tj0=Tj(N-1);
...
end
end
Also, Tj is a vector.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Walter Roberson
Walter Roberson 2018-9-6

0 个投票

k=k+1;
for k=1:N
The second of those two lines overwrites the variable 'k' referred to in the first of the lines. After the for loop, k will be left at the last value from the for loop, which is to say N . The while loop will be testing that value against 100. If N < 100, then k will never exceed 100.

6 个评论
显示 4更早的评论 隐藏 4更早的评论

Giulliana Marin
Giulliana Marin 2018-9-6
What should I do to fix it? delete the for? I just want the loop to run while (((Tj-Tj0)/Tj)^2)>=1e-10 && k<=100
Walter Roberson
Walter Roberson 2018-9-6
We do not know if you want Tj0 to be repeatedly calculated for each k value. But yes, probably you should be deleting the "for k" loop.
Giulliana Marin
Giulliana Marin 2018-9-6
Well, the main variable is Tj, so the loop should continue running until the difference between Tj and the previously calculated Tj (I wrote it as Tj0) is less than 1e-10. Actually, I made a mistake because I had written >= and it's <=. It still doesn't work though
Walter Roberson
Walter Roberson 2018-9-6
for k = 1 : 100
Kij=ctteK1(Pj,Tj,N);
%Conversion de C a F
Xij=MatrizxijT(N,Fg,Fl,Gj,Kij,Zij);
[Lj,Gj,Xij,Yij]=Flujos(N,Lj,Xij,Fg,Fl,Kij);
Tj=Tj*1.8+32;
TF=Tj*1.8+32;
%Conversion de kmol a lbmol
[MatrizTH]=TH(N,Tj,Xij,Yij,Lj,Gj);
[HFlj,HFgj]=THFeed(N,TF,Zij);
[Tj]=Balance(N, Tj, Yij, Xij, Lj, Gj, Fj, HFlj, HFgj, Qj, MatrizTH);
if abs((Tj-Tj0)/Tj) <= 1e-10; break; end
Tj0 = Tj;
end
Giulliana Marin
Giulliana Marin 2018-9-6
I get this: Undefined function or variable 'Tj0'
Walter Roberson
Walter Roberson 2018-9-6
That should not be the case, since you had
Tj0=0;
before your loop.

请先登录，再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Loops and Conditional Statements 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by