参考书籍:《自动控制原理》(第七版).胡寿松主编.
《自动控制原理PDF版下载》
典型环节分为:最小相位环节、非最小相位环节;
最小相位环节:
非最小相位环节:
开环传递函数的典型环节分解可将开环系统表示为若干个典型环节的串联形式:
G
(
s
)
H
(
s
)
=
∏
i
=
1
N
G
i
(
s
)
G(s)H(s)=\prod_{i=1}^NG_i(s)
G(s)H(s)=i=1∏NGi(s)
设典型环节的频率特性为:
G
i
(
j
ω
)
=
A
i
(
ω
)
e
j
φ
i
(
ω
)
G_i({\rm j}\omega)=A_i(\omega){\rm e}^{{\rm j}\varphi_i(\omega)}
Gi(jω)=Ai(ω)ejφi(ω)
则系统开环频率特性为:
G
(
j
ω
)
H
(
j
ω
)
=
[
∏
i
=
1
N
A
i
(
ω
)
]
e
j
[
∑
i
=
1
N
φ
i
(
ω
)
]
G({\rm j}\omega)H({\rm j}\omega)=\left[\prod_{i=1}^NA_i(\omega)\right]{\rm e}^{{\rm j}\left[\sum_{i=1}^N\varphi_i(\omega)\right]}
G(jω)H(jω)=[i=1∏NAi(ω)]ej[∑i=1Nφi(ω)]
系统开环幅频特性和开环相频特性为:
A
(
ω
)
=
∏
i
=
1
N
A
i
(
ω
)
,
φ
(
ω
)
=
∑
i
=
1
N
φ
i
(
ω
)
A(\omega)=\prod_{i=1}^NA_i(\omega),\varphi(\omega)=\sum_{i=1}^N\varphi_i(\omega)
A(ω)=i=1∏NAi(ω),φ(ω)=i=1∑Nφi(ω)
系统开环对数幅频特性为:
L
(
ω
)
=
20
lg
A
(
ω
)
=
∑
i
=
1
N
20
lg
A
i
(
ω
)
=
∑
i
=
1
N
L
i
(
ω
)
L(\omega)=20\lg{A(\omega)}=\sum_{i=1}^N20\lg{A_i(\omega)}=\sum_{i=1}^NL_i(\omega)
L(ω)=20lgA(ω)=i=1∑N20lgAi(ω)=i=1∑NLi(ω)


典型环节频率特性曲线的若干特点:
非最小相位与对应的最小相位环节
最小相位的比例环节
G
(
s
)
=
K
(
K
>
0
)
G(s)=K(K>0)
G(s)=K(K>0),简称比例环节,幅频和相频特性如下:
A
(
ω
)
=
K
,
φ
(
ω
)
=
0
°
A(\omega)=K,\varphi(\omega)=0°
A(ω)=K,φ(ω)=0°
非最小相位的比例环节
G
(
s
)
=
−
K
(
K
>
0
)
G(s)=-K(K>0)
G(s)=−K(K>0),幅频和相频特性如下:
A
(
ω
)
=
K
,
φ
(
ω
)
=
−
180
°
A(\omega)=K,\varphi(\omega)=-180°
A(ω)=K,φ(ω)=−180°
最小相位的惯性环节
G
(
s
)
=
1
T
s
+
1
,
(
T
>
0
)
G(s)=\displaystyle\frac{1}{Ts+1},(T>0)
G(s)=Ts+11,(T>0),幅频和相频特性为:
A
(
ω
)
=
1
1
+
T
2
ω
2
,
φ
(
ω
)
=
−
arctan
T
ω
A(\omega)=\frac{1}{\sqrt{1+T^2\omega^2}},\varphi(\omega)=-\arctan{T\omega}
A(ω)=1+T2ω21,φ(ω)=−arctanTω
非最小相位的惯性环节
G
(
s
)
=
1
−
T
s
+
1
,
(
T
>
0
)
G(s)=\displaystyle\frac{1}{-Ts+1},(T>0)
G(s)=−Ts+11,(T>0),幅频和相频特性为:
A
(
ω
)
=
1
1
+
T
2
ω
2
,
φ
(
ω
)
=
arctan
T
ω
A(\omega)=\frac{1}{\sqrt{1+T^2\omega^2}},\varphi(\omega)=\arctan{T\omega}
A(ω)=1+T2ω21,φ(ω)=arctanTω
传递函数互为倒数的典型环节
最小相位典型环节中,积分环节和微分环节、惯性环节和一阶微分环节、振荡环节和二阶微分环节的传递函数互为倒数,有如下关系:
G
1
(
s
)
=
1
G
2
(
s
)
G_1(s)=\frac{1}{G_2(s)}
G1(s)=G2(s)1
设
G
1
(
j
ω
)
=
A
1
(
ω
)
e
j
φ
1
(
ω
)
G_1(j\omega)=A_1(\omega){\rm e}^{{\rm j}\varphi_1(\omega)}
G1(jω)=A1(ω)ejφ1(ω),则有:
{
φ
2
(
ω
)
=
−
φ
1
(
ω
)
L
2
(
ω
)
=
20
lg
A
2
(
ω
)
=
20
lg
1
A
1
(
ω
)
=
−
L
1
(
ω
)
传递函数互为倒数的典型环节,对数幅频曲线关于
0
d
B
{\rm 0dB}
0dB线对称,对数相频曲线关于
0
°
0°
0°线对称;此结论在非最小相位环节中亦适用;
振荡环节和二阶微分环节
振荡环节的传递函数为:
G
(
s
)
=
1
(
s
/
ω
n
)
2
+
2
ζ
(
s
/
ω
n
)
+
1
;
ω
n
>
0
,
0
<
ζ
<
1
G(s)=\frac{1}{(s/\omega_n)^2+2\zeta(s/\omega_n)+1};\omega_n>0,0<\zeta<1
G(s)=(s/ωn)2+2ζ(s/ωn)+11;ωn>0,0<ζ<1
振荡环节的频率特性:
A
(
ω
)
=
1
(
1
−
ω
2
ω
n
2
)
2
+
4
ζ
2
ω
2
ω
n
2
φ
(
ω
)
=
−
arctan
2
ζ
ω
ω
n
1
−
ω
2
ω
n
2
=
{
−
arctan
2
ζ
ω
ω
n
1
−
ω
2
ω
n
2
,
ω
≤
ω
n
−
(
180
−
arctan
2
ζ
ω
ω
n
ω
2
ω
n
2
−
1
)
,
ω
>
ω
n
φ ( 0 ) = 0 ° , φ ( ∞ ) = − 180 ° \varphi(0)=0°,\varphi(\infty)=-180° φ(0)=0°,φ(∞)=−180°,相频特性曲线从 0 ° 0° 0°单调减至-180°;当 ω = ω n \omega=\omega_n ω=ωn时, φ ( ω n ) = − 90 ° \varphi(\omega_n)=-90° φ(ωn)=−90°, A ( ω n ) = 1 2 ζ , A(\omega_n)=\displaystyle\frac{1}{2\zeta}, A(ωn)=2ζ1,振荡环节与虚轴的交点为 − j 1 2 ζ -{\rm j}\displaystyle\frac{1}{2\zeta} −j2ζ1;
A
(
0
)
=
1
,
A
(
∞
)
=
0
A(0)=1,A(\infty)=0
A(0)=1,A(∞)=0,求
A
(
ω
)
A(\omega)
A(ω)的极值,
d
A
(
ω
)
d
ω
=
−
[
−
2
ω
ω
n
2
(
1
−
ω
2
ω
n
2
)
+
4
ζ
2
ω
ω
n
2
]
[
(
1
−
ω
2
ω
n
2
)
2
+
4
ζ
2
ω
2
ω
n
2
]
3
2
=
0
\frac{{\rm d}A(\omega)}{{\rm d}\omega}=\frac{-\left[-\displaystyle\frac{2\omega}{\omega_n^2}\left(1-\frac{\omega^2}{\omega_n^2}\right)+4\zeta^2\frac{\omega}{\omega_n^2}\right]}{\left[\left(1-\displaystyle\frac{\omega^2}{\omega_n^2}\right)^2+4\zeta^2\displaystyle\frac{\omega^2}{\omega_n^2}\right]^{\frac{3}{2}}}=0
dωdA(ω)=[(1−ωn2ω2)2+4ζ2ωn2ω2]23−[−ωn22ω(1−ωn2ω2)+4ζ2ωn2ω]=0
解得谐振频率:
ω
r
=
ω
n
1
−
2
ζ
2
,
0
<
ζ
≤
2
/
2
\omega_r=\omega_n\sqrt{1-2\zeta^2},0<\zeta≤\sqrt{2}/2
ωr=ωn1−2ζ2,0<ζ≤2/2
谐振峰值:
M
r
=
A
(
ω
r
)
=
1
2
ζ
1
−
ζ
2
,
0
<
ζ
≤
2
/
2
M_r=A(\omega_r)=\frac{1}{2\zeta\sqrt{1-\zeta^2}},0<\zeta≤\sqrt{2}/2
Mr=A(ωr)=2ζ1−ζ21,0<ζ≤2/2
当
0
<
ζ
<
2
2
0<\zeta<\displaystyle\frac{\sqrt{2}}{2}
0<ζ<22时,
d
M
r
d
ζ
=
−
(
1
−
2
ζ
2
)
ζ
2
(
1
−
ζ
2
)
3
2
<
0
\frac{{\rm d}M_r}{{\rm d}\zeta}=\frac{-(1-2\zeta^2)}{\zeta^2(1-\zeta^2)^{\frac{3}{2}}}<0
dζdMr=ζ2(1−ζ2)23−(1−2ζ2)<0
ω
r
,
M
r
\omega_r,M_r
ωr,Mr均为阻尼比
ζ
\zeta
ζ的减函数
(
0
<
ζ
≤
2
2
)
(0<\zeta≤\displaystyle\frac{\sqrt{2}}{2})
(0<ζ≤22);当
0
<
ζ
<
2
2
0<\zeta<\displaystyle\frac{\sqrt{2}}{2}
0<ζ<22时,且
ω
∈
(
0
,
ω
r
)
\omega\in(0,\omega_r)
ω∈(0,ωr)时,
A
(
ω
)
A(\omega)
A(ω)单调增;
ω
∈
(
ω
r
,
∞
)
\omega\in(\omega_r,\infty)
ω∈(ωr,∞)时,
A
(
ω
)
A(\omega)
A(ω)单调减;当
2
2
<
ζ
<
1
\displaystyle\frac{\sqrt{2}}{2}<\zeta<1
22<ζ<1时,
A
(
ω
)
A(\omega)
A(ω)单调减;
二阶微分环节的传递函数为振荡环节传递函数的倒数,按对称性可得二阶微分环节的对数频率曲线,有:
{
A
(
0
)
=
1
φ
(
0
)
=
0
°
,
{
A
(
ω
n
)
=
2
ζ
φ
(
ω
n
)
=
90
°
,
{
A
(
∞
)
=
∞
φ
(
∞
)
=
180
°
当阻尼比
2
2
<
ζ
<
1
\displaystyle\frac{\sqrt{2}}{2}<\zeta<1
22<ζ<1时,
A
(
ω
)
A(\omega)
A(ω)从
1
1
1单调增至
∞
\infty
∞;当阻尼比
0
<
ζ
<
2
2
0<\zeta<\displaystyle\frac{\sqrt{2}}{2}
0<ζ<22,且
ω
∈
(
0
,
ω
r
)
\omega\in(0,\omega_r)
ω∈(0,ωr)时,
A
(
ω
)
A(\omega)
A(ω)从
1
1
1单调减至
{
A
(
ω
r
)
=
2
ζ
1
−
ζ
2
<
1
ω
r
=
ω
1
−
2
ζ
2
在
ω
∈
(
ω
r
,
∞
)
\omega\in(\omega_r,\infty)
ω∈(ωr,∞)时,
A
(
ω
)
A(\omega)
A(ω)单调增;
对数幅频渐近特性曲线
为简化惯性环节、一阶微分环节、振荡环节和二阶微分环节的对数幅频曲线的作图,常用低频和高频渐近线近似表示对数幅频曲线,称为对数幅频渐近特性曲线;
惯性环节的对数幅频渐近特性为:
L
a
(
ω
)
=
{
0
,
ω
<
1
T
−
20
lg
ω
T
,
ω
>
1
T
L_a(\omega)=

低频部分是零分贝线,高频部分是斜率为 − 20 d B / d e c {\rm -20dB/dec} −20dB/dec的直线,两条直线交于 ω = 1 T \omega=\displaystyle\frac{1}{T} ω=T1处,称频率 1 T \displaystyle\frac{1}{T} T1为惯性环节的交接频率;用渐近特性近似表示对数幅频特性存在误差: Δ L ( ω ) = L ( ω ) − L a ( ω ) \Delta{L(\omega)=L(\omega)-L_a(\omega)} ΔL(ω)=L(ω)−La(ω),在交接频率处误差最大,约为 − 3 d B -3{\rm dB} −3dB;一阶微分环节和非最小相位一阶微分环节与惯性环节的对数幅频渐近特性曲线以 0 d B 0{\rm dB} 0dB线互为镜像;
振荡环节的对数幅频特性为:
L
(
ω
)
=
−
20
lg
(
1
−
ω
2
ω
n
2
)
2
+
4
ζ
2
ω
2
ω
n
2
L(\omega)=-20\lg\sqrt{\left(1-\frac{\omega^2}{\omega_n^2}\right)^2+4\zeta^2\frac{\omega^2}{\omega_n^2}}
L(ω)=−20lg(1−ωn2ω2)2+4ζ2ωn2ω2
当
ω
<
<
ω
n
\omega<<\omega_n
ω<<ωn时,
L
(
ω
)
≈
0
L(\omega)≈0
L(ω)≈0,低频渐近线为
0
d
B
0{\rm dB}
0dB线;当
ω
>
>
ω
n
\omega>>\omega_n
ω>>ωn时,
L
(
ω
)
=
−
40
lg
ω
ω
n
L(\omega)=-40\lg\displaystyle\frac{\omega}{\omega_n}
L(ω)=−40lgωnω,高频渐近线为过
(
ω
n
,
0
)
(\omega_n,0)
(ωn,0)点,斜率为
−
40
d
B
/
d
e
c
{\rm -40dB/dec}
−40dB/dec的直线,振荡环节的交接频率为
ω
n
\omega_n
ωn,对数幅频渐近特性为:
L
a
(
ω
)
=
{
0
,
ω
<
ω
n
−
40
lg
ω
ω
n
,
ω
>
ω
n
L_a(\omega)=
非最小相位振荡环节与振荡环节对数幅频渐近特性曲线相同,二阶微分环节和非最小相位二阶微分环节与振荡环节的对数幅频渐近特性曲线关于
0
d
B
{\rm 0dB}
0dB线对称;
半对数坐标系中的直线方程为:
k
=
L
a
(
ω
2
)
−
L
a
(
ω
1
)
lg
ω
2
−
lg
ω
1
k=\frac{L_a(\omega_2)-L_a(\omega_1)}{\lg\omega_2-\lg\omega_1}
k=lgω2−lgω1La(ω2)−La(ω1)
其中:
[
ω
1
,
lg
(
ω
1
)
]
、
[
ω
2
,
lg
(
ω
2
)
]
[\omega_1,\lg(\omega_1)]、[\omega_2,\lg(\omega_2)]
[ω1,lg(ω1)]、[ω2,lg(ω2)]为直线上的两点,
k
(
d
B
/
d
e
c
)
k{\rm (dB/dec)}
k(dB/dec)为直线斜率;
绘制概略开环幅相特性曲线的方法:
开环幅相特性曲线的起点 ( ω = 0 + ) (\omega=0_+) (ω=0+)和终点 ( ω = ∞ ) (\omega=\infty) (ω=∞);
开环幅相特性曲线与实轴的交点。设
ω
=
ω
x
\omega=\omega_x
ω=ωx时,
G
(
j
ω
x
)
H
(
j
ω
x
)
G({\rm j}\omega_x)H({\rm j}\omega_x)
G(jωx)H(jωx)的虚部为:
I
m
[
G
(
j
ω
x
)
H
(
j
ω
x
)
]
=
0
或
φ
(
ω
x
)
=
∠
[
G
(
j
ω
x
)
H
(
j
ω
x
)
]
=
k
π
,
k
=
0
,
±
1
,
±
2
,
…
{\rm Im}[G({\rm j}\omega_x)H({\rm j}\omega_x)]=0或\varphi(\omega_x)=\angle[G({\rm j}\omega_x)H({\rm j}\omega_x)]=k\pi,k=0,±1,±2,\dots
Im[G(jωx)H(jωx)]=0或φ(ωx)=∠[G(jωx)H(jωx)]=kπ,k=0,±1,±2,…
其中:
ω
x
\omega_x
ωx称为穿越频率;
开环频率特性曲线与实轴交点坐标值为:
R
e
[
G
(
j
ω
x
)
H
(
j
ω
x
)
]
=
G
(
j
ω
x
)
H
(
j
ω
x
)
{\rm Re}[G({\rm j}\omega_x)H({\rm j}\omega_x)]=G({\rm j}\omega_x)H({\rm j}\omega_x)
Re[G(jωx)H(jωx)]=G(jωx)H(jωx)
开环幅相特性曲线的变化范围(象限、单调性);
实例分析:
E
x
a
m
p
l
e
1
{\rm Example1}
Example1: 设系统开环传递函数为:
G
(
s
)
H
(
s
)
=
K
s
(
T
1
s
+
1
)
(
T
2
s
+
1
)
;
K
,
T
1
,
T
2
>
0
G(s)H(s)=\frac{K}{s(T_1s+1)(T_2s+1)};K,T_1,T_2>0
G(s)H(s)=s(T1s+1)(T2s+1)K;K,T1,T2>0
绘制系统概略开环幅相特性曲线。
解:
系统开环频率特性:
G
(
j
ω
)
H
(
j
ω
)
=
K
(
1
−
j
T
1
ω
)
(
1
−
j
T
2
ω
)
(
−
j
)
ω
(
1
+
T
1
2
ω
2
)
(
1
+
T
2
2
ω
2
)
=
K
[
−
(
T
1
+
T
2
)
ω
+
j
(
−
1
+
T
1
T
2
ω
2
)
]
ω
(
1
+
T
1
2
ω
2
)
(
1
+
T
2
2
ω
2
)
G({\rm j}\omega)H({\rm j}\omega)=\frac{K(1-{\rm j}T_1\omega)(1-{\rm j}T_2\omega)(-{\rm j})}{\omega(1+T_1^2\omega^2)(1+T_2^2\omega^2)}=\frac{K\left[-(T_1+T_2)\omega+{\rm j}(-1+T_1T_2\omega^2)\right]}{\omega(1+T_1^2\omega^2)(1+T_2^2\omega^2)}
G(jω)H(jω)=ω(1+T12ω2)(1+T22ω2)K(1−jT1ω)(1−jT2ω)(−j)=ω(1+T12ω2)(1+T22ω2)K[−(T1+T2)ω+j(−1+T1T2ω2)]
幅值变化:
A
(
0
+
)
=
∞
,
A
(
∞
)
=
0
A(0_+)=\infty,A(\infty)=0
A(0+)=∞,A(∞)=0;
相角变化:
∠
(
1
j
ω
)
:
−
90
°
~
−
90
°
;
∠
(
1
1
+
j
T
1
ω
)
:
0
°
~
−
90
°
∠
(
1
1
+
j
T
2
ω
)
:
0
°
~
−
90
°
,
∠
K
:
0
°
~
0
°
⇒
φ
(
ω
)
:
−
90
°
~
−
270
°
起点处:
R
e
[
G
(
j
0
+
)
H
(
j
0
+
)
]
=
−
K
(
T
1
+
T
2
)
,
I
m
[
G
(
j
0
+
)
H
(
j
0
+
)
]
=
−
∞
{\rm Re}[G({\rm j}0_+)H({\rm j}0_+)]=-K(T_1+T_2),{\rm Im}[G({\rm j}0_+)H({\rm j}0_+)]=-\infty
Re[G(j0+)H(j0+)]=−K(T1+T2),Im[G(j0+)H(j0+)]=−∞
与实轴交点:
I
m
[
G
(
j
ω
)
H
(
j
ω
)
]
=
0
,
⇒
ω
x
=
1
T
1
T
2
{\rm Im}[G({\rm j}\omega)H({\rm j}\omega)]=0,\Rightarrow\omega_x=\frac{1}{T_1T_2}
Im[G(jω)H(jω)]=0,⇒ωx=T1T21
G ( j ω x ) H ( j ω x ) = R e [ G ( j ω x ) H ( j ω x ) ] = − K T 1 T 2 T 1 + T 2 G({\rm j}\omega_x)H({\rm j}\omega_x)={\rm Re}[G({\rm j}\omega_x)H({\rm j}\omega_x)]=-\frac{KT_1T_2}{T_1+T_2} G(jωx)H(jωx)=Re[G(jωx)H(jωx)]=−T1+T2KT1T2

绘制概略开环幅相特性曲线规律小结:
开环幅相特性曲线的起点,取决于比例环节 K K K和系统积分或微分环节的个数 ν \nu ν(系统型别);
开环幅相特性曲线的终点,取决于开环传递函数分子、分母多项式中最小相位环节和非最小相位环节的阶次和;
设系统开环传递函数的分子、分母多项式的阶次分别为
m
m
m和
n
n
n,记除
K
K
K外,分子多项式中最小相位环节的阶次和为
m
1
m_1
m1,非最小相位环节的阶次和为
m
2
m_2
m2,分母多项式中最小相位环节的阶次和为
n
1
n_1
n1,非最小相位环节的阶次和为
n
2
n_2
n2,则有:
m
=
m
1
+
m
2
,
n
=
n
1
+
n
2
m=m_1+m_2,n=n_1+n_2
m=m1+m2,n=n1+n2
φ
(
∞
)
=
{
[
(
m
1
−
m
2
)
−
(
n
1
−
n
2
)
]
×
90
°
,
K
>
0
[
(
m
1
−
m
2
)
−
(
n
1
−
n
2
)
]
×
90
°
−
180
°
,
K
<
0
\varphi(\infty)=
当开环系统为最小相位系统时,若:
n
=
m
,
G
(
j
∞
)
H
(
j
∞
)
=
K
∗
n
>
m
,
G
(
j
∞
)
H
(
j
∞
)
=
0
∠
[
(
n
−
m
)
×
(
−
90
°
)
]
其中:
K
∗
K^*
K∗为系统开环根轨迹增益;
若开环系统存在等幅振荡环节,重数
l
l
l为正整数,即开环传递函数具有如下形式:
G
(
s
)
H
(
s
)
=
1
(
s
2
ω
n
2
+
1
)
l
G
1
(
s
)
H
1
(
s
)
G(s)H(s)=\frac{1}{(\displaystyle\frac{s^2}{\omega_n^2}+1)^l}G_1(s)H_1(s)
G(s)H(s)=(ωn2s2+1)l1G1(s)H1(s)
G
1
(
s
)
H
1
(
s
)
G_1(s)H_1(s)
G1(s)H1(s)不含
±
j
ω
n
±{\rm j}\omega_n
±jωn的极点,则当
ω
\omega
ω趋于
ω
n
\omega_n
ωn时,
A
(
ω
)
A(\omega)
A(ω)趋于无穷,而:
φ
(
ω
n
−
)
≈
φ
1
(
ω
n
)
=
∠
[
G
1
(
j
ω
n
)
H
1
(
j
ω
n
)
]
φ
(
ω
n
+
)
≈
φ
1
(
ω
n
)
−
l
×
180
°
即
φ
(
ω
)
\varphi(\omega)
φ(ω)在
ω
=
ω
n
\omega=\omega_n
ω=ωn附近,相角突变
−
l
×
180
°
-l\times180°
−l×180°;
系统开环对数幅频渐近特性:
L
a
(
ω
)
=
∑
i
=
1
N
L
a
i
(
ω
)
L_a(\omega)=\sum_{i=1}^NL_{a_i}(\omega)
La(ω)=i=1∑NLai(ω)
对于任意的开环传递函数,按典型环节分解,将组成系统的典型环节分为三部分:
记 ω m i n \omega_{\rm min} ωmin为最小交接频率,称 ω < ω m i n \omega<\omega_{\rm min} ω<ωmin的频率范围为低频段;
开环对数幅频渐近特性曲线绘制步骤:
开环传递函数典型环节分解;
确定一阶环节、二阶环节的交接频率,将各交接频率标注在半对数坐标图的 ω \omega ω轴上;
绘制低频段渐近特性线:在 ω < ω m i n \omega<\omega_{\rm min} ω<ωmin频段内,开环系统幅频渐近特性的斜率取决于 K ω ν \displaystyle\frac{K}{\omega^{\nu}} ωνK,因而直线斜率为 − 20 -20 −20 ν \nu ν d B / d e c {\rm dB/dec} dB/dec;获得低频渐近线,需要确定直线上的一点,方法如下:
过点 ( ω 0 , L a ( ω 0 ) ) (\omega_0,L_a(\omega_0)) (ω0,La(ω0))在 ω < ω m i n \omega<\omega_{\rm min} ω<ωmin范围内可作斜率为 − 20 ν d B / d e c -20\nu{\rm dB/dec} −20νdB/dec的直线;
作 ω ≥ ω m i n \omega≥\omega_{\rm min} ω≥ωmin频段渐近特性线;在每个交接频率点处,斜率发生变化,变化规律取决于该交接频率对应的典型环节的种类,当系统的多个环节具有相同交接频率时,该交接频率点处斜率的变化应为各个环节对应的斜率变化值的代数和;
交接频率点处斜率的变化表:

实例分析:
E
x
a
m
p
l
e
2
:
{\rm Example2:}
Example2: 已知系统开环传递函数为:
G
(
s
)
H
(
s
)
=
2000
s
−
4000
s
2
(
s
+
1
)
(
s
2
+
10
s
+
400
)
G(s)H(s)=\frac{2000s-4000}{s^2(s+1)(s^2+10s+400)}
G(s)H(s)=s2(s+1)(s2+10s+400)2000s−4000
绘制系统开环对数幅频渐近特性曲线。
解:
开环传递函数的典型环节分解形式为:
G
(
s
)
H
(
s
)
=
−
10
(
1
−
s
2
)
s
2
(
s
+
1
)
(
s
2
2
0
2
+
1
2
⋅
s
20
+
1
)
G(s)H(s)=\frac{-10(1-\displaystyle\frac{s}{2})}{s^2(s+1)(\displaystyle\frac{s^2}{20^2}+\frac{1}{2}·\frac{s}{20}+1)}
G(s)H(s)=s2(s+1)(202s2+21⋅20s+1)−10(1−2s)
开环系统由
6
6
6个典型环节串联而成:非最小相位比例环节、两个积分环节、非最小相位一阶微分环节、惯性环节和振荡环节;
S T E P 1 {\rm STEP1} STEP1:确定各交接频率 ω i , i = 1 , 2 , 3 \omega_i,i=1,2,3 ωi,i=1,2,3及斜率变化值;
非最小相位一阶微分环节: ω 2 = 2 \omega_2=2 ω2=2,斜率增加 20 d B / d e c {\rm 20dB/dec} 20dB/dec;
惯性环节: ω 1 = 1 \omega_1=1 ω1=1,斜率减小 20 d B / d e c {\rm 20dB/dec} 20dB/dec;
振荡环节: ω 3 = 20 \omega_3=20 ω3=20,斜率减小 40 d B / d e c \rm {40dB/dec} 40dB/dec;
S T E P 2 {\rm STEP2} STEP2:绘制低频段渐近特性曲线;
因为 ν = 2 \nu=2 ν=2,低频渐近线斜率 k = − 40 d B / d e c k=-40{\rm dB/dec} k=−40dB/dec,直线上一点 ( ω 0 , L a ( ω 0 ) ) = ( 1 , 20 d B ) (\omega_0,L_a(\omega_0))=(1,20{\rm dB}) (ω0,La(ω0))=(1,20dB);
S
T
E
P
3
{\rm STEP3}
STEP3:绘制中高频段渐近特性曲线;
ω
m
i
n
≤
ω
<
ω
2
,
k
=
−
60
d
B
/
d
e
c
ω
2
≤
ω
<
ω
3
,
k
=
−
40
d
B
/
d
e
c
ω
≥
ω
3
,
k
=
−
80
d
B
/
d
e
c

输出量经恒定延时后不失真地复现输入量变化的环节称为延迟环节;含有延迟环节的系统称为延迟系统;
延迟环节的输入输出时域表达式为:
c
(
t
)
=
1
(
t
−
τ
)
r
(
t
−
τ
)
c(t)=1(t-\tau)r(t-\tau)
c(t)=1(t−τ)r(t−τ)
其中:
τ
\tau
τ为延迟时间;
延迟环节的传递函数为:
G
(
s
)
=
C
(
s
)
R
(
s
)
=
e
−
τ
s
G(s)=\frac{C(s)}{R(s)}={\rm e}^{-\tau{s}}
G(s)=R(s)C(s)=e−τs
延迟环节的频率特性为:
G
(
j
ω
)
=
e
−
j
τ
ω
=
1
⋅
∠
(
−
57.3
τ
ω
)
G({\rm j}\omega)={\rm e}^{-{\rm j}\tau\omega}=1·\angle(-57.3\tau\omega)
G(jω)=e−jτω=1⋅∠(−57.3τω)
频率响应实验

传递函数确定
从低频段起,将实验所得的对数幅频特性曲线用斜率为 0 d B / d e c , ± 20 d B / d e c , ± 40 d B / d e c , ⋯ ⋯ {\rm 0dB/dec,±20dB/dec,±40dB/dec},\cdots\cdots 0dB/dec,±20dB/dec,±40dB/dec,⋯⋯直线分段近似,获得对数幅频渐近特性曲线;
由对数幅频渐近特性曲线确定最小相位条件下系统的传递函数,是对数幅频渐近特性曲线绘制的逆问题;
实例分析:
E x a m p l e 3 : {\rm Example3:} Example3: 下图由频率响应实验获得的某最小相位系统的开环对数幅频特性曲线和对数幅频渐近特性曲线,确定系统的传递函数。

解:
S T E P 1 : {\rm STEP1:} STEP1:确定系统积分或微分环节的个数。
由图可知,对数幅频渐近特性曲线低频渐近线的斜率为 + 20 d B / d e c {\rm +20dB/dec} +20dB/dec,因此 ν = − 1 \nu=-1 ν=−1,即系统含有一个微分环节;
S T E P 2 : {\rm STEP2:} STEP2:确定系统传递函数结构形式。
ω = ω 1 \omega=\omega_1 ω=ω1处,斜率变化 − 20 d B / d e c {\rm -20dB/dec} −20dB/dec,对应惯性环节;
ω = ω 2 \omega=\omega_2 ω=ω2处,斜率变化 − 40 d B / d e c {\rm -40dB/dec} −40dB/dec,可以对应振荡环节,也可以是重惯性环节,由图可知,在 ω 2 \omega_2 ω2附近存在谐振现象,因此为振荡环节;
系统传递函数形式为:
G
(
s
)
=
K
s
(
1
+
s
ω
1
)
(
s
2
ω
2
2
+
2
ζ
s
ω
2
+
1
)
G(s)=\frac{Ks}{\left(1+\displaystyle\frac{s}{\omega_1}\right)\left(\displaystyle\frac{s^2}{\omega_2^2}+2\displaystyle\zeta\frac{s}{\omega_2}+1\right)}
G(s)=(1+ω1s)(ω22s2+2ζω2s+1)Ks
其中参数
ω
1
、
ω
2
、
ζ
、
K
\omega_1、\omega_2、\zeta、K
ω1、ω2、ζ、K待定;
S T E P 3 : {\rm STEP3:} STEP3:由给定条件确定传递函数参数。
低频渐近线方程为:
L
a
(
ω
)
=
20
lg
K
ω
ν
=
20
lg
K
−
20
ν
lg
ω
L_a(\omega)=20\lg\frac{K}{\omega^{\nu}}=20\lg{K}-20\nu\lg\omega
La(ω)=20lgωνK=20lgK−20νlgω
由给定点
(
ω
,
L
a
(
ω
)
)
=
(
1
,
0
)
(\omega,L_a(\omega))=(1,0)
(ω,La(ω))=(1,0)及
ν
=
−
1
\nu=-1
ν=−1可得:
K
=
1
K=1
K=1.
根据直线方程式:
L
a
(
ω
a
)
−
L
a
(
ω
b
)
=
k
(
lg
ω
a
−
lg
ω
b
)
L_a(\omega_a)-L_a(\omega_b)=k(\lg\omega_a-\lg\omega_b)
La(ωa)−La(ωb)=k(lgωa−lgωb)
由给定点:
ω
a
=
1
,
L
a
(
ω
a
)
=
0
,
ω
b
=
ω
1
,
L
a
(
ω
b
)
=
12
,
k
=
20
\omega_a=1,L_a(\omega_a)=0,\omega_b=\omega_1,L_a(\omega_b)=12,k=20
ωa=1,La(ωa)=0,ωb=ω1,La(ωb)=12,k=20
可得:
ω
1
=
1
0
12
20
=
3.98
\omega_1=10^{\frac{12}{20}}=3.98
ω1=102012=3.98
由给顶点:
ω
a
=
100
,
L
a
(
ω
a
)
=
0
,
ω
b
=
ω
2
,
L
a
(
ω
b
)
=
12
,
k
=
−
40
\omega_a=100,L_a(\omega_a)=0,\omega_b=\omega_2,L_a(\omega_b)=12,k=-40
ωa=100,La(ωa)=0,ωb=ω2,La(ωb)=12,k=−40
可得:
ω
2
=
1
0
(
−
12
40
+
lg
100
)
=
50.1
\omega_2=10^{\left(-\frac{12}{40}+\lg100\right)}=50.1
ω2=10(−4012+lg100)=50.1
在谐振频率
ω
r
\omega_r
ωr处,振荡环节的谐振峰值为:
20
lg
M
r
=
20
lg
1
2
ζ
1
−
ζ
2
20\lg{M_r}=20\lg\frac{1}{2\zeta\sqrt{1-\zeta^2}}
20lgMr=20lg2ζ1−ζ21
根据:
20
lg
M
r
=
20
−
12
=
8
(
d
B
)
20\lg{M_r}=20-12=8(dB)
20lgMr=20−12=8(dB),因此:
M
r
=
2.512
M_r=2.512
Mr=2.512,则有:
ζ
4
−
ζ
2
+
0.04
=
0
⇒
ζ
1
=
0.204
,
ζ
2
=
0.979
\zeta^4-\zeta^2+0.04=0\Rightarrow\zeta_1=0.204,\zeta_2=0.979
ζ4−ζ2+0.04=0⇒ζ1=0.204,ζ2=0.979
因为
0
<
ζ
<
0.707
0<\zeta<0.707
0<ζ<0.707才存在谐振峰值,因此
ζ
=
0.204
\zeta=0.204
ζ=0.204;
综上,系统的传递函数为:
G
(
s
)
=
s
(
s
3.98
+
1
)
(
s
2
50.
1
2
+
0.408
⋅
s
50.1
)
+
1
G(s)=\frac{s}{\left(\displaystyle\frac{s}{3.98}+1\right)\left(\displaystyle\frac{s^2}{50.1^2}+0.408·\frac{s}{50.1}\right)+1}
G(s)=(3.98s+1)(50.12s2+0.408⋅50.1s)+1s