具有不稳定极点的传递函数的波特图

信息处理 matlab 预兆
2022-02-15 23:14:49

假设我有一个 LTI 系统,由以下传递函数

G(s)=(s+1)s(s101).
只要我按照 Matlab 代码运行

s = tf('s');
G = (s + 1)/(s*(s/10 - 1));
bode(G);

我收到以下波特图

在此处输入图像描述

量级部分符合我的预期(手绘草图)。至于相位部分,Matlab 给我的和我的预期之间存在差异。也就是说,我预计直到第一个断点ω=1rads1由于原点处的极点,相位应该是90在我看来,相移为180,但我不明白为什么。

2个回答

如果你弹出一个非常小的s,你会得到

G(s)=(s+1)s(s101)1s=jω.

那确实有一个+90或-270度的相位。

这与杆位无关,但它是一个简单的标志翻转:

G(s)=(s+1)s(1s10)

传递函数相位是分子相位减去分子相位

分子的相位可以表示为 atan(ω)

对于分母,它有点复杂。您需要将分母扩展为s210s

将分母中的替换的正值,这会产生一个具有负实值和负虚值的复数,您不能简单地使用 Atan,因为函数仅限于但是,您可以使用 atan 并将添加到分母阶段。s=jωω210jωω±π2π

因此分母相位可以表示为π+atan(10ωω2)π+atan(10ω)

最后,传递函数是分子相位减去分子相位,它给你

θ(ω)=Atan(ω)(π+(Atan(10ω))

对于的小值, ωatan(ω)0atan(10ω)π2

因此,对于的小值,ωθ(ω)3π2

编辑:谢谢马特。L 指出我原来的错误。