matlab数值解微分方程组问题请用 ode45 解下列微分方程式:dx1/dt = -x2(x3-1) dx2/dt = x1(x3-3) dx3/dt = x1x2 假设启始条件为 [1 sqrt(6) 0],时间由 0 到 200 秒.请画出在三度相位空间的曲线图.
来源:学生作业帮助网 编辑:作业帮 时间:2024/11/15 22:35:51
matlab数值解微分方程组问题请用 ode45 解下列微分方程式:dx1/dt = -x2(x3-1) dx2/dt = x1(x3-3) dx3/dt = x1x2 假设启始条件为 [1 sqrt(6) 0],时间由 0 到 200 秒.请画出在三度相位空间的曲线图.
matlab数值解微分方程组问题
请用 ode45 解下列微分方程式:dx1/dt = -x2(x3-1)
dx2/dt = x1(x3-3)
dx3/dt = x1x2
假设启始条件为 [1 sqrt(6) 0],时间由 0 到 200 秒.
请画出在三度相位空间的曲线图.
matlab数值解微分方程组问题请用 ode45 解下列微分方程式:dx1/dt = -x2(x3-1) dx2/dt = x1(x3-3) dx3/dt = x1x2 假设启始条件为 [1 sqrt(6) 0],时间由 0 到 200 秒.请画出在三度相位空间的曲线图.
在Matlab下输入:edit,然后将下面两行百分号之间的内容,复制进去,保存
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function y=zhidao_rk4_tangwei0206302(t,x)
x1=x(1);
x2=x(2);
x3=x(3);
y=[-x2*(x3-1); x1*(x3-3); x1*x2];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
在Matlab下面输入:
t_end=200;
x0=[1;sqrt(6);0];
[t,x]=ode45('zhidao_rk4_tangwei0206302',[0,t_end],x0);
plot3(x(:,1),x(:,2),x(:,3));
xlabel('x1');
ylabel('x2');
zlabel('x3');