| 查看: 2272 | 回复: 15 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
matlab解微分方程组
|
|||
|
哪位高手指点一下matlab解微分方程组！ 图片是微分方程组。 下面是我根据网上查到的解法写的程序，能运行，但是计算出来的图像完全不对。 主函数： tspan=[0:0.01:2];%求解时间为2秒也可以 tspan=[0 2]表示自动时间步长 x0=[0;0;2.4*10^(-4);20;1.15*10^(-7);0];%初值 [t,x]=ode45(@ff,tspan,x0);%用龙格库塔法求解 figure plot(t,x(:,1)) figure plot(t,x(:,2)) figure plot(t,x(:,3)) figure plot(t,x(:,4)) figure plot(t,x(:,5)) figure plot(t,x(:,6)) M文件 %微分方程 function dx=ff(t,x) dx=zeros(6,1); %x(1)=X1 %x(2)=X2 %x(3)=X3 %x(4)=X4 %x(5)=X5 %x(6)=X6 k1u=1.5*10^(-7) k1b=1.0*10^(-4) k2=1.3*10^(7) k3=7*10^(-1) k4=1.0*10^(10) k5=2.6*10^(9) k6=3*10^(3) Xc=0.5 MMT=0 dx(1)=2*k1u*x(3)+k1b*x(3)^2-k2*x(1)*x(5)+k3*x(2)*x(4)-2*k4*x(1)^2-k5*x(1)*x(2); dx(2)=k1b*x(3)^2+k2*x(1)*x(5)-k3*x(2)*x(4)-k5*x(1)*x(2)-2*k6*x(2); dx(3)=-k1u*x(3)-2*k1b*x(3)^2+k3*x(2)*x(4); dx(4)=-k1u*x(3)-k1b*x(3)^2-k3*x(2)*x(4); dx(5)=-k2*x(1)*x(5)+k6*x(2)^2; dx(6)=(k1u*x(3)+k1b*x(3)^2)*(1-Xc-MMT)*1.1111; end       |
回复此楼» 猜你喜欢
求助VISSIM破解版软件
已经有0人回复
-
求2026年在台湾举行的物理和材料领域国际学术会议信息
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有61人回复
-
求国际会议网站
已经有1人回复
-
求取一些关于纳米材料和纳米技术相关的英文PPT。
已经有0人回复
-
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有17人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有0人回复
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 计算强帖: 12
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
11楼2011-12-10 11:19:56
buxiudehun
木虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 1710.4
- 红花: 1
- 帖子: 60
- 在线: 151.9小时
- 虫号: 634805
- 注册: 2008-10-23
- 专业: 计算数学与科学工程计算
2楼2011-12-02 09:22:51
3楼2011-12-02 09:47:19
zt_chem
木虫 (正式写手)
- 应助: 12 (小学生)
- 金币: 2147.8
- 红花: 6
- 帖子: 872
- 在线: 981.6小时
- 虫号: 286920
- 注册: 2006-10-15
- 性别: GG
- 专业: 催化化学
4楼2011-12-03 10:26:51