亲爱的所有人，我有一个问题与代码附上。我不能得到for循环的任何结果，因为for循环的结果仍然是零，再加上绘制这些结果。请帮忙。

4ビュ(過去30日間)

穆罕默德萨因 2018 年 12 月 26 日

0
リンク

直接链接到这个问题

https://jp.mathworks.com/matlabcentral/answers/437264-dear-all-i-am-having-a-problem-with-the-code-enclosed-i-can-not-get-any-results-of-the-for-loop-as

編集済み: 图像分析 2019 年 3 月 29 日

                         clc
                        
                         清晰的所有
                        
                         关闭
                        
                         抽搐
                        
                         K = 273.15;%转换常数摄氏度到开尔文
                        
                         Tev = 10+K;
                        
                         Tcond = 40+K;
                        
                         %Tgen = 200+K;
                        
                         delT = .1;
                        
                         dT = delta *2;
                        
                         Mz = 1;沸石质量%，单位为kg。床1和床2质量相等
                        
                         Cz = 0.920;沸石比热% [kJ/kg K]
                        
                         Cw = 4.187;水比热% [kJ/kg K]
                        
                         Cv = 1.85;水蒸气比热% [kJ/kg K]
                        
                         Pcond = (Tcond-K);
                        
                         Pev = (Tev-K);
                        
                         k = 100;
                        
                         最大吸附剂/吸附剂比由蒸发器确定
                        
                         %温度(Tev)与吸附平衡方程函数AdsW。
                        
                         Run = 40;
                        
                         Qreg = 0(运行，1);
                        
                         Qads = 0(运行，1);
                        
                         Qsa = 0(运行，1);
                        
                         Qsve = 0(运行，1);
                        
                         Qe = 0(运行，1);
                        
                         Qa = 0(运行，1);
                        
                         q1 = 0(运行，1);
                        
                         Qic1 = 0(运行，1);
                        
                         Qd1 = 0(运行，1);
                        
                         Qtot = 0(运行，1);
                        
                         COP = 0(运行，2);
                        
                         Ta1 = 0(运行，1);
                        
                         T1new = 0 (run,1);
                        
                         Thot = 0 (run,1);
                        
                         计数=零(运行，1);
                        
                         Count2 = 0(运行，1);
                        
                         WWmin = 0(运行，1);
                        
                         Sayac = 0(运行，1);
                        
                         Temp = 0 (run,1);
                        
                         Tint = 0(运行，2);
                        
                         Ta2new = 0 (run,1);
                        
                         P_hot = 0 (run,1);
                        
                         P_cool = 0 (run,1);
                        
                         阈值= 0.01;
                        
                         Regtemp = 59;初始再生温度%
                        
                         函数AdsW = AdsW(Tgen,Pcond)
                        
                         为I = 1:长度(运行)
                        
                         Tgen = regtemp+K+i;
                        
                         Wmin = AdsW1_TP(Tgen,Pcond);
                        
                         假设吸附温度(Ta2)等于Tcond
                        
                         Wmax = AdsW1_TP(Tcond,Pev);
                        
                         Tg1 = AdsW3_PW(Pcond,Wmax);
                        
                         Ta2 = AdsW3_PW(Pev,Wmax);
                        
                         Ta1(i) = AdsW3_PW(Pev,Wmin);
                        
                         WWmin(i) = Wmin;
                        
                         T1new(i) = Ta2;%床1被认为是凉爽的
                        
                         count2(i) = (Tgen-Ta1(i))/delT;
                        
                         Qdes = 0;
                        
                         Qsd = 0;
                        
                         W2 = Wmin;
                        
                         W1 = Wmax;
                        
                         Ta2new(i) = Ta2;%用于建模质量回收
                        
                         质量回收过程
                        
                         为T = Tgen:解决:(Ta1 +解决)
                        
                         delW = (Cz + W2*Cw)*delT / AdsH(W2);
                        
                         Wnew = W2-delW;
                        
                         如果((abs((W2- wnew)/W2)) > threshold)
                        
                         delW = (Cz + (W2-delW / 2) * Cw) *解决/ AdsH (W2-delW / 2);
                        
                         Wnew = W2-delW;
                        
                         结束
                        
                         P_hot(i) = AdsW2_TW(T,Wnew);
                        
                         A = (AdsH(W1+delW/2)*delW)/(Cz + (W1+delW/2)*Cw);
                        
                         T_cool = Ta2new(i) + A;
                        
                         Ta2new(i) = T_cool;
                        
                         P_cool(i) = AdsW2_TW(T_cool,W1+delW/2);
                        
                         W1 = W1+delW/2;
                        
                         W2 = W2- delw /2;
                        
                         如果(P_cool(i) > P_hot(i))
                        
                         打破
                        
                         结束
                        
                         结束
                        
                         T1new(i) = AdsW3_PW(P_cool(i)，W1);
                        
                         Tint(i) = T1new(i);
                        
                         Tg1 = AdsW3_PW(Pcond,W1);
                        
                         Tgen = AdsW3_PW(P_hot(i)，W2);
                        
                         Wmin = W2;
                        
                         Wmax = W1;
                        
                         %热回收过程
                        
                         为T = Tgen:解决:(Ta2 +解决)
                        
                         Thot(i) = T;
                        
                         Qic = 0;
                        
                         如果(T > T1new(i) && T1new(i) < Tg1)%定义等空间冷却
                        
                         P-T-w图中BED2(3->4)的等等冷却过程
                        
                         Qic = Mz*(Cz + Cw*Wmin)*(delT);
                        
                         Qic1(i) = Qic1(i) + Qic;
                        
                         P-T-w图中BED1(1->2)的等等加热过程
                        
                         T1new(i) = T1new(i) + Qic / (Mz*(Cz + Cw*Wmax));
                        
                         % count(i) = count(i)+1;
                        
                         elseif(T > T1new(i) && T1new(i) > Tg1)
                        
                         P-T-w图中%吸收过程(2->3)
                        
                         W_before = AdsW1_TP(T,Pev);
                        
                         W_after = AdsW1_TP((t - delta)，Pev);
                        
                         每一步平均解吸热%
                        
                         H_avg = (AdsH(W_before) + AdsH(W_after))/2;
                        
                         %由于解吸而释放的热量
                        
                         Qads(i) = Mz*H_avg*(W_after-W_before);
                        
                         %由于合理冷却而释放的热量
                        
                         Qsa(i) = Mz*(Cz + Cw*(W_after-W_before))*delT;
                        
                         %由于加热蒸汽而吸收的热量
                        
                         %蒸发到吸附温度
                        
                         Qsve (i) = Mz * (W_after-W_before) * *(简历(W_after) / (0.5 * (b.．.
                        
                         (日志(Pev)——(W_after))) + (b (W_before) /.．.
                        
                         (日志(Pev)——(W_before)))) tev);
                        
                         Qa(i) = Qads(i) + Qsa(i) - Qsve(i);
                        
                         q1 (i) = q1 (i) + Qads(i) + Qsa(i) - Qsve(i);
                        
                         为Tvar = T1new(我):(dT / k): (T1new(我)+ dT)
                        
                         如果(Qa(i) > Qd1(i))
                        
                         W_before = AdsW1_TP(Tvar,Pcond);
                        
                         W_after = AdsW1_TP((Tvar+dT/k)，Pcond);
                        
                         每一步平均解吸热%
                        
                         H_avg = (AdsH(W_before) + AdsH(W_after))/2;
                        
                         %解吸热
                        
                         Qdes = Qdes + Mz*H_avg*(W_before-W_after);
                        
                         %显热
                        
                         Qsd = Qsd + Mz * (Cz +连续波* (W_before-W_after)) * (dT / k);
                        
                         Qd1(i) = Qdes + Qsd;
                        
                         Sayac (i) = Sayac (i)+1;
                        
                         其他的
                        
                         打破
                        
                         结束
                        
                         结束
                        
                         T1new(i) = T1new(i)+ (sayac(i)*(dT/k));
                        
                         其他的
                        
                         打破回收过程完成百分比
                        
                         结束
                        
                         温度(i) = T;
                        
                         结束
                        
                         Qreg(i) = Qic1(i) + q1 (i);
                        
                         %间歇循环(必需的热量，没有热量回收)
                        
                         Tgen = regtemp+K+i;
                        
                         %Wmax = AdsW1_TP(Tcond, Pev);
                        
                         Wmin = AdsW1_TP(Tgen,Pcond);
                        
                         %Tg1 = AdsW3_PW(Pcond,Wmax);
                        
                         Qih = 0;
                        
                         Qdes = 0;
                        
                         Qsd = 0;
                        
                         %必须供给吸附剂进行等空间加热的热量
                        
                         Qih = Qih + Mz*(Cz + Cw*Wmax)*(Tg1-Tint(i));
                        
                         为一号:T =解决:(Tgen-delT)
                        
                         W_before = AdsW1_TP(T,Pcond);
                        
                         W_after = AdsW1_TP((T+delT)，Pcond);
                        
                         每一步平均解吸热%
                        
                         H_avg = (AdsH(W_before) + AdsH(W_after))/2;
                        
                         %解吸热
                        
                         Qdes = Qdes + Mz*H_avg*(W_before-W_after);
                        
                         Qsd = Qsd + Mz*(Cz + Cw*(W_before-W_after))*delT;
                        
                         结束
                        
                         %来自Cacciola & Restuccia“HEX设计的影响…
                        
                         LH2O = 3172-2.4425*Tev;
                        
                         Qe(i) = Mz*(W1-W2)*(LH2O - Cw*(Tcond-Tev));
                        
                         Qd = Qdes + Qsd;
                        
                         Qtot(i) = Qd + Qih - Qreg(i);
                        
                         COP(i,1) = Tgen-K;
                        
                         COP(i,2) = Qe(i)/Qtot(i);
                        
                         Qd = 0;
                        
                         结束
                        
                         toc
                        
                         T = toc/60% #好< NOPTS >
                        
                         I = 1:run;
                        
                         情节(警察(我,1),警察(我,2),“- - -”）;
                        
                         结束