你沒給初值,我隨便設(shè)一個(gè)[1 1]
新建一個(gè)腳本文件,輸入：
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]的話計(jì)算結(jié)果全為0，檢查一下初值搞了半天還是初值問(wèn)題，t取大一點(diǎn)就是極限環(huán)了。。