此微分方程组可用ode45（）函数求得其数值解。
首先，建立自定义函数
function dy = rigid(t,y)
dy = zeros(3,1);
dy(1) = y(2) * y(3);
dy(2) = -y(1) * y(3);
dy(3) = -0.51 * y(1) * y(2);
end
其二，用ode45（）函数求出x、y值
[x,y] = ode45(@rigid,[0 12],[0 1 1]);
其三，用plot（）函数绘出x—y，x—dy/dx，x—d²y/dx²曲线图
plot(x,y(:,1),'-',x,y(:,2),'-.',x,y(:,3),'.')
运行结果
内容全部来源于网络收集，如有侵权，请联系网站删除！
QQ：24596024