给出Duffing振子matlab求解的代码

以下是用MATLAB求解Duffing振子运动方程的代码：

% Duffing振子的运动方程为：mx'' + cx' + kx + ax^3 = F(t)
% 其中，m表示质量，c表示阻尼系数，k表示弹性系数，a表示非线性系数，F(t)表示外力
% 下面的代码使用MATLAB的ode45求解Duffing振子的运动方程

% 定义Duffing振子的参数
m = 1;
c = 0.5;
k = 1;
a = 0.1;
F = 1;

% 定义Duffing振子的初始状态
x0 = 0.5;
v0 = 0;

% 定义时间范围和时间步长
tspan = [0 100];
dt = 0.1;

% 定义运动方程
f = @(t, y) [y(2); (1/m) * (F - c*y(2) - k*y(1) - a*y(1)^3)];

% 使用ode45求解运动方程
[t, y] = ode45(f, tspan, [x0 v0]);

% 绘制Duffing振子的位移-时间图像
plot(t, y(:,1));
xlabel('时间 t');
ylabel('位移 x');
title('Duffing振子位移-时间图像');

以上代码可以绘制Duffing振子的位移-时间图像。如果需要绘制其它图像，可以根据需要进行修改。