如何以编程方式找到伯德图中给定相位的幅度和频率？

Question

如果我有植物，可以说

Gp（s）= 1 /（s + 1）

我可以找到相位裕度

使用MATLAB命令

Gp = tf([1],[1 1]);
[G P] = margin(Gp);

我的问题是，如果我想知道特定增益超频中的相位超频怎么办。 我如何找到它而不希望预示情节？

通常，我可以通过命令bode（Gp）找到它，然后将鼠标移到我想知道其相位裕度的特定增益上。

对于我之前的示例，-20相过频时的增益过频为0.363。

我如何将其编写为不在伯德图中查找的命令？

提前致谢

Answer 1

您似乎误解了“ 增益”和“ 相位余量”的实际含义，而不是解释它的地方。 我认为您真正想要的是一种无需单击即可评估波特图的方法。 例如，您想知道-20 相点的大小和频率 。

让我们看一下这三种情况：

---

情况1：您知道频率，并且正在搜索幅度和相位

最简单的情况：

w = 0.363;                 % specify given frequency
[mag,phase] = bode(Gp,w)   % output of according magnitude and phase

收益：

mag =

    0.9400


phase =

  -19.9509

---

情况2：您想知道某个相位的幅度和频率

p = -20;
[mag,phase,wout] = bode(Gp);

mag_p = interp1( squeeze(phase), squeeze(mag), p)
w_p   = interp1( squeeze(phase), wout, p)

收益：

mag_p =

    0.9394

w_p =

    0.3642

---

情况3：您想知道一定幅度的相位和频率

m = 0.9394;
[mag,phase,wout] = bode(Gp);

phase_m = interp1( squeeze(mag), squeeze(phase), m)
w_m     = interp1( squeeze(mag), wout, m)

收益：

phase_m =

  -19.9998


w_m =

    0.3642

---

squeeze命令是必需的，因为1x1x...输出一个1x1x...矩阵用于相位和幅度。 您还可以使用interp1不同插值方法。