clear all Clcr=设定值;f=@(t,x)[-8/3*x(1)+x(2)*x(3);-10*x(2)+10*x(3);-x(1)*x(2)+r*x(2)-x(3)];t_final=100;x0=[0;0;1e-10][t,x]=ode45(f,[0,t_final],x0)figure(2)plot3(x(:,1),x(:,2),x(:,3));axis([10 42 -20 20 -20 20])
来源:学生作业帮助网 编辑:作业帮 时间:2024/12/02 05:08:58
clear all Clcr=设定值;f=@(t,x)[-8/3*x(1)+x(2)*x(3);-10*x(2)+10*x(3);-x(1)*x(2)+r*x(2)-x(3)];t_final=100;x0=[0;0;1e-10][t,x]=ode45(f,[0,t_final],x0)figure(2)plot3(x(:,1),x(:,2),x(:,3));axis([10 42 -20 20 -20 20])
clear all
Clc
r=设定值;
f=@(t,x)[-8/3*x(1)+x(2)*x(3);-10*x(2)+10*x(3);-x(1)*x(2)+r*x(2)-x(3)];
t_final=100;
x0=[0;0;1e-10]
[t,x]=ode45(f,[0,t_final],x0)
figure(2)
plot3(x(:,1),x(:,2),x(:,3));
axis([10 42 -20 20 -20 20])
clear all Clcr=设定值;f=@(t,x)[-8/3*x(1)+x(2)*x(3);-10*x(2)+10*x(3);-x(1)*x(2)+r*x(2)-x(3)];t_final=100;x0=[0;0;1e-10][t,x]=ode45(f,[0,t_final],x0)figure(2)plot3(x(:,1),x(:,2),x(:,3));axis([10 42 -20 20 -20 20])
1.清除工作空间中的所有变量
2.清除指令窗中显示内容
3.设置参数r
4.f为函数句柄
5.设置终止时间为100
6.函数的初值
7.求微分方程组的数值解
8.图像2
9.画三维图
10.设置三个坐标轴的显示范围