代码只是分配一堆变量,然后退出。它永远不会调用任何函数的文件。为什么不呢?这就是为什么没有得到策划——你永远不会称为函数做任何阴谋。
亲爱的,我有一个问题,封闭的代码。我不能得到任何结果for循环的for循环的结果仍然0 +策划这些结果。任何帮助请。
2视图(30天)
显示旧的评论
clc
清晰的所有
关闭
抽搐
K = 273.15;%转换常数摄氏开尔文
Tev = 10 + K;
Tcond = 40 + K;
% Tgen = 200 + K;
解决= 1;
dT =解决* 2;
Mz = 1;%的沸石质量公斤。Bed1和Bed2等于质量
Cz = 0.920;%的比热沸石(kJ /公斤K)
连续波= 4.187;%比热的水(kJ /公斤K)
简历= 1.85;%比热的水蒸气(kJ /公斤K)
Pcond = (Tcond-K);
Pev = (Tev-K);
k = 100;
%最大吸附物/吸附剂比发现通过蒸发器
%的温度函数AdsW (Tev)和吸附平衡方程。
运行= 40;
Qreg = 0(运行,1);
qad = 0(运行,1);
Qsa = 0(运行,1);
Qsve = 0(运行,1);
量化宽松政策= 0(运行,1);
Qa = 0(运行,1);
Qa1 = 0(运行,1);
Qic1 = 0(运行,1);
Qd1 = 0(运行,1);
Qtot = 0(运行,1);
警察= 0(运行,2);
Ta1 = 0(运行,1);
T1new = 0(运行,1);
Thot = 0(运行,1);
数= 0(运行,1);
是从= 0(运行,1);
WWmin = 0(运行,1);
sayac = 0(运行,1);
Temp = 0(运行,1);
色= 0(运行,2);
Ta2new = 0(运行,1);
P_hot = 0(运行,1);
P_cool = 0(运行,1);
阈值= 0.01;
regtemp = 59;%初始再生温度
函数AdsW = AdsW (Tgen Pcond)
为i = 1:长度(运行)
Tgen = regtemp + K + i;
Wmin = AdsW1_TP (Tgen Pcond);
%假设吸附温度等于Tcond (Ta2)
Wmax = AdsW1_TP (Tcond, Pev);
一号= AdsW3_PW (Pcond Wmax);
Ta2 = AdsW3_PW (Pev Wmax);
Ta1 (i) = AdsW3_PW (Pev Wmin);
WWmin (i) = Wmin;
T1new (i) = Ta2;%床1被认为是凉爽的
是从(i) = (Tgen-Ta1(我))/解决;
qd = 0;
Qsd = 0;
W2 = Wmin;
W1 = Wmax;
Ta2new (i) = Ta2;%用于建模大规模复苏
%质量恢复过程
为T = Tgen:解决:(Ta1 +解决)
delW = (Cz + W2 * Cw) *解决/ AdsH (W2);
Wnew = W2-delW;
如果((abs ((W2-Wnew) / W2)) >阈值)
delW = (Cz + (W2-delW / 2) * Cw) *解决/ AdsH (W2-delW / 2);
Wnew = W2-delW;
结束
P_hot (i) = AdsW2_TW (T, Wnew);
= (AdsH (W1 + delW / 2) * delW) / (Cz + (W1 + delW / 2) * Cw);
T_cool = Ta2new (i) +一个;
Ta2new (i) = T_cool;
P_cool (i) = AdsW2_TW (T_cool W1 + delW / 2);
W1 = W1 + delW / 2;
W2 = W2-delW / 2;
如果(P_cool(我)> P_hot(我))
打破
结束
结束
T1new (i) = AdsW3_PW (P_cool(我),W1);
色彩(i) = T1new(我);
一号= 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) <一号)%定义等排的冷却
%等排的冷却过程对BED2 P-T-w图(3 - > 4)
Qic = Mz * (Cz +连续波* Wmin) *(解决);
Qic1 (i) = Qic1 (i) + Qic;
%等排的加热过程BED1 P-T-w图(1 - > 2)
T1new (i) = T1new (i) + Qic / (Mz * (Cz +连续波* Wmax));
%数(i) = count (i) + 1;
elseif(T > T1new (i) & & T1new (i) >一号)
%(2 - > 3)吸收过程P-T-w图
W_before = AdsW1_TP (T, Pev);
W_after = AdsW1_TP ((Pev T-delT));
%的平均热解吸的每一步
H_avg = (AdsH (W_before) + AdsH (W_after)) / 2;
% relased由于热解吸
qad系统(i) = Mz * H_avg * (W_after-W_before);
%释放热量的合理的冷却
Qsa (i) = Mz * (Cz +连续波* (W_after-W_before)) *解决;
%的热量吸收由于加热蒸汽
%蒸发吸附温度
Qsve (i) = Mz * (W_after-W_before) * *(简历(W_after) / (0.5 * (b…
(日志(Pev)——(W_after))) + (b (W_before) /…
(日志(Pev)——(W_before)))) tev);
Qa (i) = qad系统(我)+ Qsa (i) - Qsve(我);
Qa1 (i) = Qa1(我)+ qad系统(我)+ Qsa (i) - Qsve(我);
为Tvar = T1new(我):(dT / k): (T1new(我)+ dT)
如果(Qa (i) > Qd1(我))
W_before = AdsW1_TP (Tvar Pcond);
W_after = AdsW1_TP ((Pcond Tvar + dT / k));
%的平均热解吸的每一步
H_avg = (AdsH (W_before) + AdsH (W_after)) / 2;
%热解吸
qd = qd + Mz * H_avg * (W_before-W_after);
%显热
Qsd = Qsd + Mz * (Cz +连续波* (W_before-W_after)) * (dT / k);
Qd1 (i) = qd + Qsd;
sayac (i) = sayac(我)+ 1;
其他的
打破
结束
结束
T1new (i) = T1new(我)+ (sayac (i) * (dT / k));
其他的
打破%恢复过程完成
结束
临时(i) = T;
结束
Qreg (i) = Qic1(我)+ Qa1(我);
% intermiitent周期(必要的热不热回收)
Tgen = regtemp + K + i;
% Wmax = AdsW1_TP (Tcond, Pev);
Wmin = AdsW1_TP (Tgen Pcond);
一号= % AdsW3_PW (Pcond Wmax);
Qih = 0;
qd = 0;
Qsd = 0;
%的热量,必须提供吸附剂等排的加热
Qih = Qih + Mz * (Cz +连续波* Wmax) * (Tg1-Tint(我));
为一号:T =解决:(Tgen-delT)
W_before = AdsW1_TP (T, Pcond);
W_after = AdsW1_TP ((T +解决),Pcond);
%的平均热解吸的每一步
H_avg = (AdsH (W_before) + AdsH (W_after)) / 2;
%热解吸
qd = qd + Mz * H_avg * (W_before-W_after);
Qsd = Qsd + Mz * (Cz +连续波* (W_before-W_after)) *解决;
结束
从Cacciola % & Restuccia“六角设计的影响…
LH2O = 3172 - 2.4425 * Tev;
量化宽松政策(i) = Mz * (W1-W2) * (LH2O - Cw * (Tcond-Tev));
Qd = Qd + Qsd;
Qtot (i) = Qd + Qih - Qreg(我);
警察(我,1)= Tgen-K;
警察(我,2)=量化宽松(i) / Qtot(我);
Qd = 0;
结束
toc
t = toc / 60% #好< NOPTS >
我= 1:运行;
情节(警察(我,1),警察(我,2),“- - -”);
结束
