| 查看: 2864 | 回复: 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       |
回复此楼» 猜你喜欢
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
物理学I论文润色/翻译怎么收费?
已经有242人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有11人回复
-
散金币祈福
已经有94人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用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人回复
7楼2011-12-10 10:43:54
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
20