| 查看: 444 | 回复: 1 | |||
| 本帖产生 1 个 博学EPI ，点击这里进行查看 | |||
zj34331197木虫 (正式写手)
木虫
|
[求助]
微分方程求解
|
||
|
微分方程求解： 方程如下： da/dT=20*exp(-70/T)*(1-a)^1.03 采用MATLAB如何求取a~T曲线。附上MATLAB的M文件代码和命令窗口代码。 感谢各位大神了。 |
回复此楼» 猜你喜欢
一篇MDPI论文改变了学习工作和生活
已经有4人回复
-
一个化合物的合成路线：CAS:367929-02-0 名称：8β-乙烯基雌二醇
已经有4人回复
-
太白金星有点烦
已经有3人回复
-
中国地质大学（北京）博士招生补录，数理学院材料科学与工程专业和材料与化工专业
已经有3人回复
-
河北省自然基金
已经有8人回复
-
收到国自然专家邀请后几年才会有本子送过来评
已经有3人回复
-
考博
已经有3人回复
-
有没有快的中文核心比较快录用的，纳米材料光催化
已经有4人回复
-
本人42，博士刚毕业，现在找不到工作，怎么办？:(
已经有21人回复
-
有人投过CCC中国控制会议吗？
已经有3人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
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