实战宝典 原创大赛 仪器问答 仪友会 总帖：8042824今日发帖：25 我的社区 签到
快速导航 采购仪器 提醒
数据处理
版主:
入住本版
专家:
居民列表

主题：【原创】MATLAB在化工领域的应用 举例

浏览0 回复4 电梯直达
省部重点实验室
头像
结帖率：
100%
关注：0 |粉丝：0
新手级： 新兵
发表于：2012/03/10 23:51:06 楼主 管理 分享 倒序浏览 只看楼主 回复 私聊
在可视为全混流的流化床反应器中迸行原料油催化裂化反应

            A（原料油，g） B（汽油，g）十C（裂化气，g）

反应器进料中含80％原料油和20％惰性组分。在反应过程中催化剂会因积炭而失活。裂化反应可近似看作原料油浓度的一级反应：

               

为您推荐
近期热榜
热门活动
您可能想找: 气相色谱仪(GC) 询底价
选参数 看心得 找厂商 查方案
专属顾问快速对接
立即提交
相关阅读
省部重点实验室
gl19860312
结帖率：
100%
关注：0 |粉丝：0
新手级： 新兵
2012/3/10 23:52:22 沙发 管理 分享 倒序浏览 只看楼主 回复 私聊
完整的MATLAB程序如下所列：

function decayinCSTR



global ct0 yA0 tao k kd

ct0=1.0；

yA0＝0.8；

tao＝0.02；

k＝45；

kd＝9；

c0＝［0.8 1］；

t0＝0；t1＝0.5；

[t,c]＝ode45（@decay，［t0 t1］,c0）；

nt＝length(c(:,1));

for i=1:nt

  x(i)=1-(1+yA0)/(1+c(i,1)/ct0)*c(i,1)/ct0/yA0;

end

plot(t,c(:,1),’b-‘,t,c(:2),’b--‘,t,x,’b*’)

xlabel(’time(h)’)

ylabel(‘cA,a,x’)

legend(‘cA’,’a’.’x’)



function dcdt=decay(t,c)



global ct0 yA0 tao k kd

dcdt(1)=ct0*yA0/tao-((1+yA0)/(1+c(1)/ct0/yA0)+c(2)*tao*k)/tao*c(1);

dcdt(2)=-kd*c(1)*c(2);

dcdt=[dcdt(1);dcdt(2)];







运行上述程序得反应器出口原料油浓度、转化率和反应器中催化剂活性随

时间的变化如下图所示。
赞贴
0
收藏
0
拍砖
0
举报
展清风
cztshou
结帖率：
100%
关注：0 |粉丝：0
新手级： 新兵
2012/7/18 20:14:38 板凳 管理 分享 倒序浏览 只看楼主 回复 私聊
.......想学学minitab.......没看懂.......
赞贴
0
收藏
0
拍砖
0
举报
zsj201204
zsj201204
结帖率：
100%
关注：0 |粉丝：0
新手级： 新兵
2012/9/28 10:47:13 马扎 管理 分享 倒序浏览 只看楼主 回复 私聊
现在已经看不懂了，以前还学过的
赞贴
0
收藏
0
拍砖
0
举报
yang198000301
yang198000301
结帖率：
100%
关注：0 |粉丝：0
新手级： 新兵
2012/9/29 15:31:41 地毯 管理 分享 倒序浏览 只看楼主 回复 私聊
非常好，谢谢
赞贴
0
收藏
0
拍砖
0
举报
手机版： MATLAB在化工领域的应用 举例
猜你喜欢最新推荐热门推荐更多推荐
头像
高级回复 发表评论
品牌合作伙伴
岛津
日立科学仪器
珀金埃尔默仪器（上海）有限公司（PerkinElmer）
日本电子株式会社
丹纳赫
安捷伦
赛默飞世尔科技
普析通用
欧波同
天美
天瑞仪器
德国耶拿
海能技术
马尔文帕纳科
磐诺科技
上海仪电科仪
梅特勒托利多
聚光科技
莱伯泰科
盛瀚
多宁生物
丹东百特
科哲
卓立汉光
屹尧科技
华谱科仪
宝德仪器
优莱博
HORIBA
布鲁克核磁