MATLAB中如何使用传递函数s表示一阶系统并绘制频率响应曲线？

Question

在MATLAB中如何正确表示一阶系统的传递函数并绘制其频率响应曲线？ 假设一阶系统传递函数为 \( G(s) = \frac{1}{\tau s + 1} \)，其中 \( \tau \) 是时间常数。如何使用 `tf` 函数定义该传递函数，并通过 `bode` 或 `freqs` 函数绘制频率响应曲线？需要注意哪些参数设置以确保结果准确，例如频率范围的选择和单位转换（如从弧度到赫兹）？此外，如何解读幅频特性和相频特性曲线中的关键点，比如截止频率和相位裕度？

Answer 1

1. MATLAB中一阶系统传递函数的表示

在MATLAB中，使用`tf`函数可以方便地定义传递函数。对于一阶系统传递函数 \( G(s) = \frac{1}{\tau s + 1} \)，我们可以通过指定分子和分母系数来定义。

        tau = 1; % 时间常数
        num = [1]; % 分子系数
        den = [tau 1]; % 分母系数
        G = tf(num, den); % 定义传递函数
    

上述代码中，`num`表示分子系数向量，`den`表示分母系数向量。通过`tf`函数，我们可以将这些系数转换为传递函数对象。

2. 频率响应曲线的绘制

频率响应曲线可以通过`bode`或`freqs`函数绘制。以下是两种方法的实现：

• 使用`bode`函数：直接生成幅频特性和相频特性曲线。
• 使用`freqs`函数：适合离散频率点分析，需要手动设置频率范围。

以下代码展示了如何使用`bode`函数绘制频率响应：

        bode(G); % 绘制Bode图
        grid on; % 添加网格
    

如果需要自定义频率范围，可以使用以下代码：

        w = logspace(-2, 2, 1000); % 设置频率范围（从0.01到100 rad/s）
        bode(G, w); % 使用指定频率范围绘图
    

3. 参数设置与单位转换

为了确保结果准确，需注意以下参数设置：

参数	说明
频率范围	通常选择对数分布的频率范围，例如`logspace(-2, 2, 1000)`。
单位转换	Bode图默认以弧度/秒为单位。若需转换为赫兹，可使用`bodeplot`并调整选项。

示例代码如下：

        opts = bodeoptions; % 创建Bode图选项
        opts.FreqUnits = 'Hz'; % 设置频率单位为赫兹
        bodeplot(G, w, opts); % 使用选项绘图
    

4. 关键点解读

在幅频特性和相频特性曲线中，关键点包括：

1. 截止频率：幅频特性下降到-3dB时对应的频率。
2. 相位裕度：当增益为1时，相位距离-180°的角度差。

以下是获取截止频率和相位裕度的代码：

        [Gm, Pm, Wcg, Wcp] = margin(G); % 获取增益裕度、相位裕度及对应频率
        fprintf('截止频率: %.2f rad/s\n', Wcg);
        fprintf('相位裕度: %.2f degrees\n', Pm);
    

5. 分析流程图

以下是整个分析过程的流程图：

graph TD;
            A[定义时间常数] --> B[创建传递函数];
            B --> C[设置频率范围];
            C --> D[绘制Bode图];
            D --> E[获取关键点信息];