请给出一个MATLAB示例，要求用ode45求解单自由度的振动微分方程，在每个频率下会得到该频率对应的位移响应向量，如何将这个频率下对应的位移响应向量得到一个二维点，包含这个频率与该频率下的振幅，最后在一定频率范围内得到一个二维的幅频曲线 请给出具体的MATLAB代码

好的，下面是一个MATLAB的示例代码：

% 定义微分方程
f = @(t, y) [y(2); -y(1)-0.2*y(2)];

% 定义频率范围
frequencies = linspace(0, 5, 51);

% 定义时间范围
tspan = [0, 20];

% 初始化振幅数组
amplitudes = zeros(size(frequencies));

% 求解微分方程
for i = 1:length(frequencies)
    % 定义频率对应的初始条件
    y0 = [1; 0];
    w = 2*pi*frequencies(i);
    y0(2) = w*y0(1);
    
    % 求解微分方程
    [t, y] = ode45(f, tspan, y0);
    
    % 计算振幅
    amplitudes(i) = max(abs(y(:, 1)));
end

% 绘制幅频曲线
figure;
plot(frequencies, amplitudes);
xlabel('Frequency');
ylabel('Amplitude');
title('Amplitude-Frequency Response');

在代码中，我们先定义了微分方程，然后定义了频率范围和时间范围。接着，我们使用一个循环，在每个频率下求解微分方程，然后计算振幅。最后，我们将频率和振幅作为二维点，绘制出幅频曲线。

这段代码的输出结果应该是一个幅频曲线，横坐标是频率，纵坐标是振幅。