来自毕文青的问题
求matlab解一阶微分方程组,并画相图.dx/dt=y;dy/dt=-2.6*y^3+0.6*y-0.1*x
求matlab解一阶微分方程组,并画相图.
dx/dt=y;
dy/dt=-2.6*y^3+0.6*y-0.1*x
4回答
2020-12-29 04:42
求matlab解一阶微分方程组,并画相图.dx/dt=y;dy/dt=-2.6*y^3+0.6*y-0.1*x
求matlab解一阶微分方程组,并画相图.
dx/dt=y;
dy/dt=-2.6*y^3+0.6*y-0.1*x
你没给初值,我随便设一个[11]
新建一个脚本文件,输入:
clearall
clc
f=@(t,x)([x(2);-2.6*x(2)^3+0.6*x(2)-0.1*x(1)]);
[t,x]=ode45(f,[01],[11]);
plot(t,x(:,1),t,x(:,2))
xlabel('t')
legend('x','y')
初值[00],可不可以横轴x竖轴Y
可以clearallclcf=@(t,x)([x(2);-2.6*x(2)^3+0.6*x(2)-0.1*x(1)]);[t,x]=ode45(f,[00.1],[00]);plot(x(:,1),x(:,2))xlabel('x'),ylabel('y')初值是[00]的话计算结果全为0,检查一下初值
搞了半天还是初值问题,t取大一点就是极限环了。。