24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

_alarrn90

无虫 (小有名气)

应助: 0 (幼儿园)
金币: -28
红花: 1
帖子: 291
在线: 118.5小时
虫号: 1016060
注册: 2010-05-11
性别: GG
专业: 高分子材料的加工与成型

[求助] 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人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
物理学I论文润色/翻译怎么收费? 已经有210人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有23人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

用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人回复

1楼 2011-12-01 18:53:32

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

_alarrn90

无虫 (小有名气)

应助: 0 (幼儿园)
金币: -28
红花: 1
帖子: 291
在线: 118.5小时
虫号: 1016060
注册: 2010-05-11
性别: GG
专业: 高分子材料的加工与成型

确实是氧气浓度的初值有问题，我忘记乘上环境压力，但是修改后的图趋势没什么变化。扩散的影响在这里可以忽略。[PH]的初值确实是20.问一下您对照下面微分方程的话，程序本身有没有错误？因为计算出来的图的趋势与文献明显不同。另外我想问一下初值以及步长的设定对图的变化趋势影响大不大，或者还有没有其他影响图形趋势的因素？非常感谢。