| 查看: 405 | 回复: 1 | |||
| 本帖产生 1 个 博学EPI ，点击这里进行查看 | |||
zj34331197木虫 (正式写手)
木虫
|
[求助]
微分方程求解
|
||
|
微分方程求解： 方程如下： da/dT=20*exp(-70/T)*(1-a)^1.03 采用MATLAB如何求取a~T曲线。附上MATLAB的M文件代码和命令窗口代码。 感谢各位大神了。 |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有4人回复
-
基金申报
已经有5人回复
-
基金委咋了？2026年的指南还没有出来？
已经有7人回复
-
国自然申请面上模板最新2026版出了吗？
已经有17人回复
-
纳米粒子粒径的测量
已经有8人回复
-
疑惑？
已经有5人回复
-
计算机、0854电子信息（085401-058412）调剂
已经有5人回复
-
Materials Today Chemistry审稿周期
已经有5人回复
-
溴的反应液脱色
已经有7人回复
-
推荐一本书
已经有12人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
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