我想利用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 |