5、MATLAB数值计算与数据分析

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
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值