24小时热门版块排行榜    

查看: 2363  |  回复: 15
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

_alarrn90

无虫 (小有名气)

[求助] 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















回复此楼

» 猜你喜欢

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

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

_alarrn90

无虫 (小有名气)

多谢,我试一下。
3楼2011-12-02 09:47:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

_alarrn90

无虫 (小有名气)

多谢,但是我感觉加不加end没什么影响吧。去掉了end计算出来的结果也不太对,你还看出其他什么问题没有?
5楼2011-12-03 11:26:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

_alarrn90

无虫 (小有名气)

确实是氧气浓度的初值有问题,我忘记乘上环境压力,但是修改后的图趋势没什么变化。扩散的影响在这里可以忽略。[PH]的初值确实是20.问一下您对照下面微分方程的话,程序本身有没有错误?因为计算出来的图的趋势与文献明显不同。另外我想问一下初值以及步长的设定对图的变化趋势影响大不大,或者还有没有其他影响图形趋势的因素?非常感谢。
7楼2011-12-10 10:43:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

_alarrn90

无虫 (小有名气)

这是模拟出的图和文献中图的对比





9楼2011-12-10 11:09:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

_alarrn90

无虫 (小有名气)

那您觉得怎么调整可能会好一些?
10楼2011-12-10 11:14:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

_alarrn90

无虫 (小有名气)

这个确实应该是时间的问题,运行的时候出现警告:Warning: Failure at t=1.847468e-002.  Unable to meet integration tolerances without reducing the step size below the smallest value allowed (5.551115e-017) at time t.
> In ode45 at 355,我在论坛搜了一下好像没有统一的解决办法。
动力学常数的单位都是s,这个问题您有什么建议?


12楼2011-12-10 15:04:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

_alarrn90

无虫 (小有名气)

具体微分方程组及相关参数在http://g.zhubajie.com/urllink.php?id=11829823oelq8zutbei57rpi 这里,麻烦了。
13楼2011-12-10 15:07:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 _alarrn90 的主题更新
信息提示
请填处理意见