你没给初值,我随便设一个[1 1]
新建一个脚本文件,输入：
clear all
clc
f=@(t,x)([x(2);-2.6*x(2)^3+0.6*x(2)-0.1*x(1)]);
[t,x]=ode45(f,[0 1],[1 1]);
plot(t,x(:,1),t,x(:,2))
xlabel('t')
legend('x','y')初值[0 0]，可不可以横轴x竖轴Y可以 clear allclcf=@(t,x)([x(2);-2.6*x(2)^3+0.6*x(2)-0.1*x(1)]);[t,x]=ode45(f,[0 0.1],[0 0]);plot(x(:,1),x(:,2))xlabel('x'),ylabel('y') 初值是[0 0]的话计算结果全为0，检查一下初值搞了半天还是初值问题，t取大一点就是极限环了。。