龙格库塔法matlab编程
发布网友
发布时间：2022-05-29 21:40
共1个回答
热心网友
时间：2024-10-05 02:46
用matlab的四、五阶龙格库塔方法（ode45）可以求得其数值解。设定其微分方程组的初始条件X(0)=0，S(0)=0。解决方法：
1、常微分方程组的自定义函数odefun（t，y），其主要内容
dy(1)=3*y(2)/(60+y(2))*y(1)-0.16*y(1);
dy(2)=4*250-4*y(2)-5*y(2)/(60+y(2))*y(1)-0.06*y(1);
2、t、y1、y2的求解
[t,y]=ode45(@odefun,[0 10],y0)
3、计算结果及t与S、X的关系图，S与X的关系图
追问能不能把完整程序发一下
