MATLAB数值计算与数据分析
1. 乙烷 - 蒸汽裂解反应模拟
下面是一个用于模拟乙烷 - 蒸汽裂解反应的MATLAB脚本 ethanrxn.m :
% ethanrxn.m: ethane-steam cracking reaction
x0 = [0.001 0.001 0.001 0.993 1 0.0001 5.992 1 0.001 10 10 10]';
T = 1000;
[x, fval] = fsolve(@rxnfun, x0,[],T);
comp = {'CH4','C2H4','C2H2','CO2','CO','O2','H2','H2O','C2H6'};
lamda = {'lambda1','lambda2','lambda3'};
fprintf('\n i\tComp. \t Initial Val.\t\tFinal val.\n');
for k = 1:length(x)-3
fprintf('%d\t%s \t%12.9f\t%12.9f\n',k,comp{k},x0(k),x(k));
end
fprintf('\n i\tLambda\tInitial Val.\tFinal val.\n');
for i = k+1:length(x)
fprintf('%g\t%s\t\t%4.1f\t%15.9f\n',i,lamda{i-length(comp)},x0(i),x(i));
end
运行该脚本后,得到如下结果:
Equation solved.
fsolve com
超级会员免费看
订阅专栏 解锁全文
2万+

被折叠的 条评论
为什么被折叠?



