把你上面的过程整理成matlab代码

好的，以下是一个完整的 MATLAB 代码示例，用于求解非线性振动微分方程：

% 定义常数
m = 1; % 质量
c = 0.1; % 阻尼系数
k = 1; % 弹性系数

% 定义外力函数
F = @(t) sin(t);

% 定义微分方程组
odefun = @(t,y) [y(2); (F(t) - c*y(2) - k*sin(y(1)))/m];

% 定义初始条件
y0 = [0; 0];

% 求解微分方程组
[t,y] = ode45(odefun, [0 10], y0);

% 绘制位移-时间图像
plot(t,y(:,1))
xlabel('时间 t')
ylabel('位移 y')
title('非线性振动')

在这个示例中，我们定义了常数 m、c、k，并定义了外力函数 F(t)。然后，我们定义了微分方程组 odefun，以及初始条件 y0。最后，我们使用 ode45 函数求解微分方程组，并绘制出非线性振动的位移-时间图像。