利用MATLAB软件设计一套程序，能够用于一阶测试系统响应特性和二阶测试系统响应特性的仿真计算，实现带入任意时间常数、固有频率、阻尼比、和被测频率，能够计算振幅比和相位角并画出图形

MATLAB是一个强大的数值计算和可视化工具，非常适合用来设计这样的信号处理程序。对于一阶和二阶系统的响应特性，你可以创建函数来分别模拟这两个系统的动态行为。以下是简要步骤：

1. **一阶系统**:
- 设计一个`first_order_response`函数，接受输入参数：时间常数τ (time constant)、初始条件（如零状态响应的初始值）以及被测频率ω0。
- 计算响应函数，通常是指数衰减形式 (`y = A * exp(-t/τ) * cos(ω0*t + Φ)`)，其中振幅A、相位Φ依赖于初始条件。
- 计算振幅比（幅值响应的最大值除以稳态值）和相位角。
- 使用`plot`命令绘制时间域内的波形。

2. **二阶系统**:
- 创建`second_order_response`函数，除了τ、ω0外，还需要阻尼比ζ (damping ratio)。
- 根据二阶微分方程的解，可以采用典型形式 `y = e^(-ζ*ωn*t) * (Acos(ω_d*t + Φ) + Bsin(ω_d*t))`，其中ωn是自然频率，ω_d = ω0 - ζ*ωn。
- 振幅比和相位角的计算稍微复杂些，需要额外的数学处理。
- 绘制Bode图（幅频响应和相频响应），这通常涉及`bode`函数。

**编写示例**:

function [amp_ratio, phase_angle] = first_order_response(τ, initial_condition, ω0)
    % ... 具体计算过程 ...
end

function [response, bode_data] = second_order_response(τ, ζ, ω0, initial_conditions)
    % ... 算法包括计算振幅比和相位角，以及绘制Bode图 ...
end

% 调用函数
[amp_ratio1, phase_angle1] = first_order_response(1, 1, 1);
[response2, bode_data2] = second_order_response(0.5, 0.7, 1, [1; 0]);

% 绘图
figure;
plot(time_vector, response2); % 时间域波形
figure;
bode(bode_data2);            % Bode图