| 查看: 2035 | 回复: 4 | ||
xiaofenger01金虫 (初入文坛)
|
[求助]
如何该这段MATLAB程序（把一个常量改成变量） 已有1人参与
|
|
原题目是这样的： 在反应器中进行液相反应制备产物B，反应可在180~260℃的温度范围内进行，反应物X大量过剩，而C, D和E为副产物。各反应均为一级动力学关系：r＝－kC，式中已知参数：k01＝5.78052×1010，k02＝3.92317×1012，k03＝1.64254×104，k04＝6.264×108，Ea1=124670，Ea2=150386，Ea3=77954，Ea4=111528。初始浓度：CA=1kmol/m3，其余物质浓度为0。已知是产物B收率最大的最优反应温度为224.6℃ 试计算在最优反应温度（224.6℃）下各组分浓度随时间的动态变化（这一段程序如下），如果我想看从180℃升温到（每分钟10℃）224.6℃，并且在224.6℃恒温10分钟再升温到260℃（每分钟10℃）这段过程中的各组分浓度随时间的动态变化应该怎么改？（就是把程序里的T从常量改成变量的意思） function Cha5demo4 T = 224.6 + 273.15; % 1℃ = 1+273.15K R = 8.31434; k0 = [5.78052E+10 3.92317E+12 1.64254E+4 6.264E+8]; Ea = [124670 150386 77954 111528]; C0 = [1 0 0 0 0]; % Initial concentration:C0(i), kmol/m^3 tspan = [0 1e4]; opt=odeset('reltol',1e-4,'outputfcn','odephas2','outputsel',[1;4]) [t,C] = ode45(@MassEquations, tspan, C0,opt,k0,Ea,R,T) plot(t,C(:,1),'r-',t,C(:,2),'k:',t,C(:,3),'b-.',t,C(:,4),'k--'); xlabel('Time (s)'); ylabel('Concentration (kmol/m^3)'); legend('A','B','C','D') % ------------------------------------------------------------------ function dCdt = MassEquations(t,C,k0,Ea,R,T) k = k0.*exp(-Ea/(R*T)); k(5) = 2.16667E-04; % Reaction rate constants, 1/s rA = -(k(1)+k(2))*C(1); rB = k(1)*C(1)-k(3)*C(2);rC = k(2)*C(1)-k(4)*C(3); rD = k(3)*C(2)-k(5)*C(4); rE = k(4)*C(3)+k(5)*C(4); % Reaction rates, kmoles/m3 s dCdt = [rA; rB; rC; rD; rE]; % Mass balances |
回复此楼» 猜你喜欢
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
物理学I论文润色/翻译怎么收费?
已经有284人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有23人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何求这个matlab程序画出的legend曲线的法矢量？
已经有6人回复
- 求MATLAB批量读取文件代码，并赋值指定的变量名
已经有4人回复
- 投稿需提交Matlab代码，不知道要怎么写代码解释等信息！
已经有9人回复
- matlab中断循环
已经有9人回复
- 如何用Matlab拟合大量散点的边界曲线
已经有5人回复
- matlab如何编程由因变量求解自变量
已经有5人回复
- 使用VC编写一个数据采集程序遇到一个问题，求指教
已经有5人回复
- 算法程序用matlab，C语言，软件平台，硬件平台运行耗时之间的关系
已经有9人回复
- 请问这段小matlab代码怎么修改？
已经有10人回复
- 最小二乘法的matlab 实现
已经有4人回复
- 怎样把这段代码变成matlab代码
已经有5人回复
- 如何改变matlab默认的变量最大数值
已经有3人回复
- 如何知道一个MATLAB程序在计算中占的内存是多大？谢谢！
已经有3人回复
- 帮忙详解示范一下MATLAB中ginput函数的应用吧
已经有3人回复
- Matlab数据如何连续写入Excel？
已经有5人回复
- 求助：资源优化问题的遗传算法matlab编程
已经有5人回复
- 请matlab高手看一下这个问题怎样编程才好，附程序代码，谢谢！
已经有9人回复
- matlab语句问题
已经有5人回复
- matlab txt文件数据提取转换
已经有10人回复
- 大家爱帮我看看这段MATLAB程序为啥出错了
已经有5人回复
- matlab生成随机行走网格 模拟材料生长
已经有26人回复
- 100BB急求matlab编程中的问题解决，如何修改这个程序呢？
已经有4人回复
- 看不懂matlab程序该怎么办？
已经有7人回复
 本帖仅楼主可见 |
2楼2015-04-03 16:50:25
感谢参与，应助指数 +1
|
本帖内容被屏蔽 |
3楼2015-04-03 21:18:53
xiaofenger01
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 709.2
- 帖子: 44
- 在线: 92.7小时
- 虫号: 2318750
- 注册: 2013-03-05
- 性别: MM
- 专业: 化学反应工程
4楼2015-04-03 23:02:27
xiaofenger01
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 709.2
- 帖子: 44
- 在线: 92.7小时
- 虫号: 2318750
- 注册: 2013-03-05
- 性别: MM
- 专业: 化学反应工程
5楼2015-04-03 23:02:37