| 查看: 489 | 回复: 1 | |||
| 本帖产生 1 个 博学EPI ,点击这里进行查看 | |||
zj34331197木虫 (正式写手)
木虫
|
[求助]
微分方程求解
|
||
|
微分方程求解: 方程如下: da/dT=20*exp(-70/T)*(1-a)^1.03 采用MATLAB如何求取a~T曲线。附上MATLAB的M文件代码和命令窗口代码。 感谢各位大神了。 |
» 猜你喜欢
售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急
已经有8人回复
售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
申博发邮件
已经有4人回复
售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急
已经有3人回复
售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
上海工程技术大学激光智能制造课题组招收博士研究生
已经有6人回复
售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急
已经有3人回复
售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急
已经有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









回复此楼
=[];plot(x, kk, 'k*'); %获得da/dT~T曲线。