clc
clear all
close all
y0=[0 0];
alpha=1;
beta=-1;
epsilon=1
delta=0.2;
omega=1;
gamma1=0.1;
option=odeset('RelTol',1e6);
[t,y]=ode45(@(t,y)duffing_02(t,y,alpha,beta,delta,epsilon,omega,gamma1),[0 600],[1 1],option);
比如我想输出上面程序t=[100,600]的部分,谢谢!!!
用上面代码得到的t和y,后续处理:
ind=find(t>=100&t<=600);即为满足题主所要求时间范围的数据,后续怎样画图由题主自行决定了,例如画相平面图很可能就是
plot(y(:,1),y(:,2))追问谢谢大神,分全给你了!