第1个回答 2019-12-08
这说明此微分方程没有显示解,应用数值方法求解,例如ode45等函数
就一个例子(选自matlab帮助文档)
解入下微分方程组
1,建立方程组函数
function
dy
=
rigid(t,y)dy
=
zeros(3,1);
%
a
column
vectordy(1)
=
y(2)
*
y(3);dy(2)
=
-y(1)
*
y(3);dy(3)
=
-0.51
*
y(1)
*
y(2);
2,求解并绘图
[T,Y]
=
ode45(@rigid,[0
12],[0
1
1]);plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.')
结果