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

回复此楼

» 猜你喜欢

求助VISSIM破解版软件已经有0人回复
求2026年在台湾举行的物理和材料领域国际学术会议信息已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有138人回复
求国际会议网站已经有1人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有19人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！已经有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的回帖

dbb627

荣誉版主 (著名写手)

专家经验: +4
计算强帖: 12
应助: 289 (大学生)
贵宾: 0.589
金币: 24640.4
散金: 551
红花: 61
沙发: 1
帖子: 1246
在线: 1794.8小时
虫号: 149791
注册: 2005-12-29
性别: GG
专业: 污染控制化学
管辖: 计算模拟

【答案】应助回帖

引用回帖:

7楼: Originally posted by _alarrn90 at 2011-12-10 10:43:54:
确实是氧气浓度的初值有问题，我忘记乘上环境压力，但是修改后的图趋势没什么变化。扩散的影响在这里可以忽略。[PH]的初值确实是20.问一下您对照下面微分方程的话，程序本身有没有错误？因为计算出来的图的趋势与 ...

程序本身没看出来啥问题，初值和步长对图形趋势很有影响的，这个步长取太大，有些极值，拐点就跳过了