24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3019  |  回复: 0

七名镇南

铁虫 (初入文坛)

[求助] matlab求解隐式常微分方程组问题

我想利用matlab ode15i来求解图1中这个常微分方程组,建立相应的Matlab函数:
function fun=c7eqdae(t,y,dy)
fun=[dy(1)-sin(y(3));
dy(2)+cos(y(4));
dy(3)-dy(4);
1/(19-32*y(1))-4/3*pi*y(2)*(2*sin(y(4))+1)*(y(1)+2/3*y(2)*cos(y(4))*sin(y(3)))];
将M矩阵输入给MATLAB工作空间,并在命令窗口中给出如下命令:
clear all;
y0=[0.3;0.027;pi/2;pi/2];dy0=[1;0;0;0]
res=ode15i('c7eqdae',[0.6 0],y0,dy0);
plot(res.x,res.y)

代入matlab中,提示我的t步长过小,无法求解。求大神帮忙解答一下,问题在哪?谢谢了

matlab求解隐式常微分方程组问题
1.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 七名镇南 的主题更新
信息提示
请填处理意见