| 查看: 432 | 回复: 1 | ||
| 本帖产生 1 个 博学EPI ，点击这里进行查看 | ||
zj34331197木虫 (正式写手)
木虫
|
[求助]
微分方程求解
|
|
|
微分方程求解： 方程如下： da/dT=20*exp(-70/T)*(1-a)^1.03 采用MATLAB如何求取a~T曲线。附上MATLAB的M文件代码和命令窗口代码。 感谢各位大神了。 |
回复此楼» 猜你喜欢
085701求调剂
已经有8人回复
-
321求调剂
已经有11人回复
-
26考研调剂0710 0860
已经有12人回复
-
270化工求调剂
已经有9人回复
-
一志愿0817化学工程与技术，求调剂
已经有23人回复
-
309求调剂
已经有7人回复
-
一志愿北京化工大学，初试成绩350求调剂
已经有7人回复
-
081700化学工程与技术 一志愿中海洋 323 求调剂学校
已经有13人回复
-
求调剂：085600材料与化工，考材科基，总分319
已经有24人回复
-
土木304求调剂
已经有4人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
zj34331197(金币+100, 博学EPI+1): ★★★很有帮助 用前部分代码可以获得a T曲线，后边代码没有用。 2012-03-11 11:16:10
zj34331197(金币+100, 博学EPI+1): ★★★很有帮助 用前部分代码可以获得a T曲线，后边代码没有用。 2012-03-11 11:16:10
|
function dy=g0(x,y) data=20 dy=data*exp(-70/T)*(1-y)^1.03; 在命令窗口中输入如下代码： [x,num_y]=ode45('g0',[0,300],0);y=abs(num_y); plot(x, y, 'k'); %获得a~T曲线; kk=100*diff(y)./diff(x);x(1,  =[];plot(x, kk, 'k*'); %获得da/dT~T曲线。由上式计算出T0-T1范围内转化率随温度的变化曲线。 动力学方程数值解，初值选取对求得da/dT~T曲线的最高点的影响。 |
2楼2012-03-11 10:59:02
