生成duffing振子matlab

生成Duffing振子的MATLAB代码如下：

% 定义常数
alpha = -1;
beta = 1;
delta = 0.3;
omega = 1;

% 定义函数
f = @(t, x) [x(2); -delta*x(2)-alpha*x(1)-beta*x(1)^3+omega*cos(t)];

% 定义时间范围和初始条件
tspan = [0, 1000];
x0 = [0.1; 0];

% 求解微分方程
[t, x] = ode45(f, tspan, x0);

% 绘制相图
plot(x(:, 1), x(:, 2));
xlabel('x');
ylabel('dx/dt');
title('Duffing Oscillator Phase Portrait');

在这个代码中，我们首先定义了Duffing振子的常数，然后定义了微分方程的右侧函数f。接着，我们指定了时间范围和初始条件，使用ode45求解微分方程，得到时间和位置的数组t和x。最后，我们绘制了相图，展示了Duffing振子的相轨迹。