| 查看: 2520 | 回复: 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       |
回复此楼» 猜你喜欢
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组招收2026硕士研究生1名
已经有0人回复
-
山西大同大学物理学硕士研究生招收调剂生
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有153人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
［调剂信息］211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
一志愿物理所凝聚态物理352求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有2人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
-
070200求调剂，一志愿某211，288分
已经有18人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用matlab求解一偏微分方程，求助！
已经有1人回复
- 流体力学一维定常不可压缩流体ns方程，已知压差变化，怎么求解速度u关于x的常微分方程
已经有0人回复
- matlab书籍
已经有103人回复
- 求大神帮忙解题啊！！！
已经有4人回复
- Matlab求解微分方程里面有参数带入怎么求解，求救！！！
已经有0人回复
- 求助CFD软件可以自定义模拟算法
已经有0人回复
- 【紧急求助】 带有微分平方项的ODE求解
已经有11人回复
- 在控制器设计中关于大家如何给定？如果给定是方波怎么给？
已经有2人回复
- 请教高手：PDE工具箱能解两个偏微分方程相互迭代求解吗
已经有2人回复
- 【求助】求解形如y''+a*y'^2+by'+c=0 二阶微分方程
已经有14人回复
- 【资料】一些数学实验,有程序源代码
已经有7人回复
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
5楼2011-12-03 11:26:16
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 计算强帖: 12
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
6楼2011-12-09 23:39:40
7楼2011-12-10 10:43:54
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 计算强帖: 12
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
8楼2011-12-10 11:08:07
9楼2011-12-10 11:09:40
10楼2011-12-10 11:14:33
