MATLAB代码求解二阶常微分方程

以下是MATLAB代码实现二阶常微分方程的求解：

% 定义初始条件和参数
y0 = 0; % 初始位置
v0 = 1; % 初始速度
tspan = [0 10]; % 时间区间
omega = 2; % 角频率
f = @(t,y) [y(2); -omega^2*y(1)]; % 定义常微分方程

% 求解常微分方程
[t,y] = ode45(f, tspan, [y0 v0]);

% 绘制解
plot(t, y(:,1), '-');
xlabel('时间');
ylabel('位移');
title('二阶常微分方程的解');

其中，`y0` 和 `v0` 分别表示初始位置和初始速度，`tspan` 表示时间区间，`omega` 表示角频率，`f` 表示常微分方程。`ode45` 函数用于求解常微分方程，返回时间和解的数组 `t` 和 `y`。最后，使用 `plot` 函数绘制解。