用matlab举例一个模拟谐振系统运动的振荡矩阵，并画图

下面是一个简单的谐振系统的振荡矩阵模拟和画图的示例代码：

% 定义系统参数
m = 1;      % 质量
k = 1;      % 劲度系数
omega0 = sqrt(k/m);     % 自然角频率
tspan = [0, 10];        % 模拟时间范围
y0 = [0.5; 0];          % 初始条件

% 定义系统方程
f = @(t, y) [y(2); -omega0^2*y(1)];

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

% 画图
figure;
plot(t, y(:, 1), 'b', 'LineWidth', 2);
xlabel('时间');
ylabel('位移');
title('谐振系统的位移随时间的变化');

运行上面的代码，可以得到如下的图像：


该图展示了谐振系统的位移随时间的变化，可以看到系统在自然角频率下做简谐振动。