| 查看: 2356 | 回复: 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 |
回复此楼» 猜你喜欢
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有124人回复
-
物理学I论文润色/翻译怎么收费?
已经有267人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有5人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有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
- 在线: 93.6小时
- 虫号: 2318750
- 注册: 2013-03-05
- 性别: MM
- 专业: 化学反应工程
4楼2015-04-03 23:02:27
xiaofenger01
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 709.2
- 帖子: 44
- 在线: 93.6小时
- 虫号: 2318750
- 注册: 2013-03-05
- 性别: MM
- 专业: 化学反应工程
5楼2015-04-03 23:02:37
200
