mathcad 传递函数

在 Mathcad 中处理传递函数（Transfer Function），通常涉及到以下几个关键步骤：

定义拉普拉斯变量 s: 传递函数是复变量 s（拉普拉斯算子）的函数。你需要告诉 Mathcad s 是一个复数变量。
- 方法： 输入 s:j 或者 s:1j。这表示将 s 定义为虚数单位 i 或 j（Mathcad 默认使用 i 或 j）。这样 Mathcad 就知道 s 代表复数频率变量。
表达传递函数： 传递函数通常是一个有理函数（分子分母都是 s 的多项式）。Mathcad 有几种方式来表示它：
- 直接分式形式： 这是最常见的形式。
```
G(s) := (分子多项式) / (分母多项式)
```
  - 例子 (一阶系统)：
```
G1(s) := 1 / (0.5 * s + 1)   // 时间常数形式
```
  - 例子 (二阶系统)：
```
G2(s) := (s + 1) / (s^2 + 5 * s + 6) // 零点在 s=-1, 极点在 s=-2, s=-3
```
- 使用系数向量 (polyroots / polyval 辅助)： 如果你想用系数向量 (aₙ, aₙ₋₁, ..., a₀) 表示多项式 aₙsⁿ + aₙ₋₁sⁿ⁻¹ + ... + a₀，可以：
```
num_coeffs := [1, 1]  // 对应于 [s, 常数项] -> s + 1
den_coeffs := [1, 5, 6]  // 对应于 [s², s, 常数项] -> s² + 5s + 6
G2(s) := polyval(num_coeffs, s) / polyval(den_coeffs, s)
```
  - 注意： polyval(v, x) 函数计算多项式 v₀xⁿ + v₁xⁿ⁻¹ + ... + vₙ。确保系数向量是按 s 的降幂顺序排列。polyroots(v) 可用于求根（极点/零点）。
计算和评估：
- 定义一个具体的复频率点 (例如 ω := 2 rad/s)，可以计算传递函数在该点的复数值：
```
ω := 2
point_value := G2(j * ω) // s = jω，j 是虚数单位
```
  point_value 将是一个复数，包含幅值和相位信息。
- 创建频率向量 ω_vec，并计算对应的频率响应向量：
```
ω_vec := 0.1, 0.2..100  // 定义频率范围
H_vec := G2(j * ω_vec)  // 计算每个ω点的复数值
```
分析和可视化： 计算得到复数值 H_vec 后，就可以生成工程上常用的图：
- 伯德图 (Bode Plot):
  - 幅频特性：20 * log10(|H_vec|) vs. log(ω_vec) 或 ω_vec (横轴常用对数坐标)
  - 相频特性：arg(H_vec) * (180/π) (度) vs. log(ω_vec) 或 ω_vec (横轴常用对数坐标)
- 奈奎斯特图 (Nyquist Plot):
  - Im(H_vec) vs. Re(H_vec) (虚部 vs. 实部)
- 阶跃响应 (Step Response):
  - 需要计算反拉普拉斯变换。Mathcad 有 invlaplace 变换（在符号计算区域使用），或者使用控制系统模块（如有）模拟。
- 零极点图 (Pole-Zero Plot):
  - 使用 polyroots(den_coeffs) 求极点位置 p。
  - 使用 polyroots(num_coeffs) 求零点位置 z。
  - 在复平面上绘制 p 和 z (Re(p), Im(p) 和 Re(z), Im(z))。

重要提示：

符号 vs. 数值计算： 上述分式定义 (G(s):= ...) 主要适用于数值计算（比如计算特定频率点的响应或绘制响应曲线）。对于符号运算（如拉普拉斯反变换），你需要使用 Mathcad 的符号计算区域（通常以 -> 或 Ctrl + . 启动）和 invlaplace 关键字。
控制系统模块： 较新版本的 Mathcad Prime 可能包含专门的控制系统模块(需要加载项)，该模块提供了创建传递函数对象（如 tf）、直接计算阶跃/脉冲响应、绘制伯德图/奈奎斯特图等更专业的函数。如果可用，这是更推荐的方法。
多项式顺序： 使用向量系数时 (polyval)，务必注意系数顺序是按 s 的降幂排列（最高次幂系数在前）。

简单示例模板：

// 定义拉普拉斯变量 s
s := j

// 定义传递函数 G(s) = (s + 1) / (s^2 + 5s + 6)
G(s) := (s + 1) / (s^2 + 5 * s + 6)

// 定义计算频率点 (单位：rad/s)
ω_point := 1
response_at_point := G(j * ω_point)
amplitude_at_point := |response_at_point|  // 幅值
phase_at_point := arg(response_at_point) * (180 / π) // 相位(度)

// 生成频率向量绘制伯德图
ω_min := 0.1
ω_max := 100
N := 1000
ω_log := logspace(ω_min, ω_max, N) // 生成对数间隔频率点 (如果 Mathcad 有此函数，否则需手动构造)
// 或者线性步进： ω_vec := ω_min, ω_min + (ω_max-ω_min)/(N-1) .. ω_max
H_vec := G(j * ω_log)

// 计算幅度 (dB) 和相位 (度)
Mag_dB := 20 * log10(|H_vec|)
Phase_deg := arg(H_vec) * (180 / π)

// 绘制伯德图 - 幅频图
// 在图上放置 X-Y 绘图区域，设置 X 轴为 log(ω_log) 或直接选对数坐标，Y 轴为 Mag_dB

// 绘制伯德图 - 相频图
// 在另一个图上放置 X-Y 绘图区域，设置 X 轴为 log(ω_log) 或直接选对数坐标，Y 轴为 Phase_deg

// (可选) 计算并绘制阶跃响应（使用符号运算或控制系统模块更简单）

希望这个中文解释能帮你在 Mathcad 中有效处理传递函数！关键在于正确定义 s 和用分式形式表达分子分母的多项式。