计算二阶系统幅值裕度与相角裕度的代码

以下是计算二阶系统幅值裕度与相角裕度的Python代码：

import numpy as np
import control.matlab as cm

# 定义二阶系统的传递函数
num = [1]
den = [1, 0.1, 1]
sys = cm.tf(num, den)

# 计算幅值裕度和相角裕度
gm, pm, wg, wp = cm.margin(sys)

print("幅值裕度：", gm)
print("相角裕度：", pm)

其中，`num` 和 `den` 分别表示系统的分子多项式和分母多项式，`sys` 表示系统的传递函数。`gm` 表示系统的幅值裕度，`pm` 表示系统的相角裕度，`wg` 表示系统的幅值裕度交叉频率，`wp` 表示系统的相角裕度交叉频率。