该系列博客主要讲述Matlab软件在自动控制方面的应用,如无自动控制理论基础,请先学习自动控制系列博文,该系列博客不再详细讲解自动控制理论知识。
自动控制理论基础相关链接:https://blog.csdn.net/qq_39032096/category_10287468.html?spm=1001.2014.3001.5482
博客参考书籍:《MATLAB/Simulink与控制系统仿真》。
1.频域分析法理论基础
1.1 频率特性基本概念
1.1.1 频率特性概述
- 控制系统的频率特性反映的是系统对正弦输入信号的响应性能;
- 频域分析法可以通过试验直接求得频率特性来分析系统的品质,应用频率特性分析系统可以得出定性和定量的结论,并具有明显的物理含义;
- 采用频率特性作为数学模型来分析和设计系统的方法称为频率特性法,亦称频率响应法;
- 频率响应法的基本思想:把控制系统中的各个变量看成一些信号,这些信号又由很多不同频率的正弦信号合成,各个变量的运动就是系统对各个不同频率的信号的响应的总和;
1.1.2 频率特性定义
频率特性:指系统在正弦信号作用下,稳态输出与输入之比相对频率的关系特性;
频率特性函数:
G
(
j
ω
)
=
X
o
(
j
ω
)
X
i
(
j
ω
)
=
A
(
ω
)
e
j
φ
(
ω
)
G({\rm j}\omega)=\frac{X_o({\rm j}\omega)}{X_i({\rm j}\omega)}=A(\omega){\rm e}^{{\rm j}\varphi(\omega)}
G(jω)=Xi(jω)Xo(jω)=A(ω)ejφ(ω)
其中:
A
(
ω
)
=
X
o
(
ω
)
X
i
(
ω
)
A(\omega)=\displaystyle\frac{X_o(\omega)}{X_i(\omega)}
A(ω)=Xi(ω)Xo(ω)称为幅频特性,
φ
(
ω
)
=
φ
o
(
ω
)
−
φ
i
(
ω
)
\varphi(\omega)=\varphi_o(\omega)-\varphi_i(\omega)
φ(ω)=φo(ω)−φi(ω)称为相频特性;
频率特性亦表示为:
G
(
j
ω
)
=
X
o
(
j
ω
)
X
i
(
j
ω
)
=
p
(
ω
)
+
j
θ
(
ω
)
G({\rm j}\omega)=\frac{X_o({\rm j}\omega)}{X_i({\rm j}\omega)}=p(\omega)+{\rm j}\theta(\omega)
G(jω)=Xi(jω)Xo(jω)=p(ω)+jθ(ω)
其中:
p
(
ω
)
p(\omega)
p(ω)为
G
(
j
ω
)
G({\rm j}\omega)
G(jω)的实部,称为实频特性,
θ
(
ω
)
\theta(\omega)
θ(ω)为
G
(
j
ω
)
G({\rm j}\omega)
G(jω)的虚部,称为虚频特性;
有如下关系:
{
p
(
ω
)
=
A
(
ω
)
cos
φ
(
ω
)
,
θ
(
ω
)
=
A
(
ω
)
sin
φ
(
ω
)
A
(
ω
)
=
p
2
(
ω
)
+
θ
2
(
ω
)
,
φ
(
ω
)
=
arctan
θ
(
ω
)
p
(
ω
)
\begin{cases} &p(\omega)=A(\omega)\cos\varphi(\omega),&\theta(\omega)=A(\omega)\sin\varphi(\omega)\\\\ &A(\omega)=\sqrt{p^2(\omega)+\theta^2(\omega)},&\varphi(\omega)=\arctan\displaystyle\frac{\theta(\omega)}{p(\omega)} \end{cases}
⎩
⎨
⎧p(ω)=A(ω)cosφ(ω),A(ω)=p2(ω)+θ2(ω),θ(ω)=A(ω)sinφ(ω)φ(ω)=arctanp(ω)θ(ω)
当输入为非正弦周期信号,其输入可利用傅里叶级数展开成正弦波的叠加,其输出为相应的正弦波的叠加;此时系统频率特性定义为系统输出量的傅氏变换与输入量的傅氏变换之比;
1.1.3 频域分析法特点
-
适用于各环节、开环和闭环系统的性能分析。
运用奈奎斯特稳定判据,通过作图方法,可以根据系统开环频率特性分析闭环系统的稳定性及性能,不必求解系统的特征根,从而避免直接求解微分方程的困难;
-
频率特性有明确的物理意义。
很多元部件频率特性可以用实验方法确定,对于机理复杂或机理不明难以写出微分方程的元部件或系统,在实验室中采用信号发生器和精密测量仪器,可以测出其频率特性;
-
频域性能指标和时域性能指标有确定的对应关系。
对于二阶系统,频率特性与时域过渡过程性能指标有确定的对应关系,对于高阶系统,通过把系统参数和结构的变化与时域过渡过程指标联系起来,存在着近似的对应关系;
-
频率设计可兼顾动态响应和噪声抑制两方面的要求。
当系统在某些频率范围内存在严重的噪声时,应用频率分析法可以设计处能满意地抑制这些噪声的系统;
-
在校正方法中,频域分析法校正十分方便。
当系统的性能指标以幅值裕度、相位裕度和误差系数等形式给出时,采用频率分析法来分析和设计系统十分方便;
-
频率法不能全面分析非线性系统。
频率法主要应用于单输入单输出的线性定常系统的分析研究中,在多输入多输出的线性定常系统中也有应用;在非线性中只有某些局部典型的应用,不能对非线性系统进行全面的分析;
1.1.4 频域性能指标
- 谐振频率 ω r \omega_r ωr:表示幅频特性 A ( ω ) A(\omega) A(ω)出现最大值时对应的频率;
- 谐振峰值 M r M_r Mr:表示幅频特性的最大值, M r M_r Mr值大表明系统对频率的正弦信号反应强烈,即系统平稳性差,阶跃响应超调量大;
- 频带 ω b \omega_b ωb:表示幅频特性 A ( ω ) A(\omega) A(ω)的幅值衰减到初始值的 0.707 0.707 0.707倍对应的频率; ω b \omega_b ωb大表明系统复现快速变化信号的能力强,失真小,即系统快速性好,阶跃响应上升时间短,调节时间短;
- 零频 A ( 0 ) A(0) A(0):表示频率 ω = 0 \omega=0 ω=0时的幅值; A ( 0 ) A(0) A(0)表示系统阶跃响应的终值, A ( 0 ) A(0) A(0)与 1 1 1间的差反映系统的稳态精度, A ( 0 ) A(0) A(0)越接近 1 1 1,系统精度越高;
1.2 频率特性的表示方法
1.2.1 极坐标图(Nyquist图)
-
系统频率特性表示:
G ( j ω ) = A ( ω ) e j φ ( ω ) G({\rm j}\omega)=A(\omega){\rm e}^{{\rm j}\varphi(\omega)} G(jω)=A(ω)ejφ(ω) -
用向量表示某一频率 ω i \omega_i ωi下的 G ( j ω i ) G({\rm j}\omega_i) G(jωi)向量的长度 A ( ω i ) A(\omega_i) A(ωi),向量极坐标角为 φ ( ω i ) \varphi(\omega_i) φ(ωi), φ ( ω i ) \varphi(\omega_i) φ(ωi)正方向取逆时针,选极坐标与直角坐标重合,极坐标的顶点在坐标原点,如下图所示:
-
频率特性 G ( j ω ) G({\rm j}\omega) G(jω)是输入频率 ω \omega ω的复变函数,当频率 ω \omega ω由 0 → ∞ 0\to\infty 0→∞时, G ( j ω ) G({\rm j}\omega) G(jω)变化的曲线,即向量端点轨迹,称为极坐标图;
-
极坐标图中,当 ω = ω i \omega=\omega_i ω=ωi时,在实轴上的投影即为实频特性 p ( ω i ) p(\omega_i) p(ωi),在虚轴上的投影即为虚频特性;
1.2.2 对数坐标图(Bode图)
-
B o d e {\rm Bode} Bode图由对数幅频特性和对数相频特性组成,如下图所示:
-
对数幅频特性是频率特性的对数值 L ( ω ) = 20 lg A ( ω ) {\rm L(\omega)=20\lg{A(\omega)}} L(ω)=20lgA(ω)与频率 ω \omega ω的关系曲线,对数相频特性是频率特性的相角 φ ( ω ) \varphi(\omega) φ(ω)与频率 ω \omega ω的关系曲线;
-
对数幅频特性纵轴为 L ( ω ) = 20 lg A ( ω ) {\rm L(\omega)}=20\lg{A(\omega)} L(ω)=20lgA(ω),单位为 d B {\rm dB} dB(分贝),采用线性分度, A ( ω ) {\rm A(\omega)} A(ω)每增加 10 10 10倍, L ( ω ) {\rm L(\omega)} L(ω)增加 20 d B 20{\rm dB} 20dB,横坐标采用对数分度,即横轴上的 ω \omega ω取对数后为等分点;
-
对数相频特性横轴采用对数分度,纵轴采用线性分度,单位为°(度);
-
对数坐标图( B o d e {\rm Bode} Bode图)的优点:
- 将幅值相乘化为对数相加运算,极大地简化了系统频率特性的绘制工作;
- 横轴采用对数分度,缩小了比例尺,扩大了频率视野,可以在较大的频段范围内表示系统频率特性;在一张 B o d e {\rm Bode} Bode图上,既可以绘制频率特性的中、高频段,也能绘制低频段,有利于系统的分析和设计;
- 可以绘制渐近的对数幅频特性,亦可制作标准样板,绘制精确的对数频率特性;
1.2.3 对数幅相图(Nichols图)
-
对数幅相图亦称为 N i c h o l s {\rm Nichols} Nichols图,是将对数幅频特性和相频特性在角频率为参变量情况下合成为一张图;
-
对数幅相图如下所示:
-
N i c h o l s {\rm Nichols} Nichols图特点:纵轴为 L ( ω ) = 20 lg A ( ω ) {\rm L(\omega)=20\lg{A(\omega)}} L(ω)=20lgA(ω),单位为 d B {\rm dB} dB(分贝),采用线性分度;横坐标采用对数分度,单位为°(度),频率 ω \omega ω为参变量;
1.3 典型环节的频率特性
1.3.1 典型环节
典型环节分为:最小相位环节和非最小相位环节;
最小相位环节:
- 比例环节: K , ( K > 0 ) K,(K>0) K,(K>0);
- 惯性环节: 1 / ( T s + 1 ) , ( T > 0 ) 1/(Ts+1),(T>0) 1/(Ts+1),(T>0);
- 一阶微分环节: T s + 1 , ( T > 0 ) Ts+1,(T>0) Ts+1,(T>0);
- 振荡环节: 1 / ( s 2 / ω n 2 + 2 ζ s / ω n + 1 ) , ( ω n > 0 , 0 ≤ ζ < 1 ) 1/(s^2/\omega_n^2+2\zeta{s}/\omega_n+1),(\omega_n>0,0≤\zeta<1) 1/(s2/ωn2+2ζs/ωn+1),(ωn>0,0≤ζ<1);
- 二阶微分环节: s 2 / ω n 2 + 2 ζ s / ω n + 1 , ( ω n > 0 , 0 ≤ ζ < 1 ) s^2/\omega_n^2+2\zeta{s}/\omega_n+1,(\omega_n>0,0≤\zeta<1) s2/ωn2+2ζs/ωn+1,(ωn>0,0≤ζ<1);
- 积分环节: 1 / s 1/s 1/s;
- 微分环节: s s s;
非最小相位环节:
- 比例环节: K , ( K < 0 ) K,(K<0) K,(K<0);
- 惯性环节: 1 / ( − T s + 1 ) , ( T > 0 ) 1/(-Ts+1),(T>0) 1/(−Ts+1),(T>0);
- 一阶微分环节: − T s + 1 , ( T > 0 ) -Ts+1,(T>0) −Ts+1,(T>0);
- 振荡环节: 1 / ( s 2 / ω n 2 − 2 ζ s / ω n + 1 ) , ( ω n > 0 , 0 < ζ < 1 ) 1/(s^2/\omega_n^2-2\zeta{s}/\omega_n+1),(\omega_n>0,0<\zeta<1) 1/(s2/ωn2−2ζs/ωn+1),(ωn>0,0<ζ<1);
- 二阶微分环节: s 2 / ω n 2 − 2 ζ s / ω n + 1 , ( ω n > 0 , 0 < ζ < 1 ) s^2/\omega_n^2-2\zeta{s}/\omega_n+1,(\omega_n>0,0<\zeta<1) s2/ωn2−2ζs/ωn+1,(ωn>0,0<ζ<1);
开环传递函数的典型环节分解可将开环系统表示为若干个典型环节的串联形式:
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(j\omega)=A_i(\omega){\rm e}^{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(j\omega)H(j\omega)=\left[\prod_{i=1}^NA_i(\omega)\right]{\rm e}^{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(ω)
- 系统开环频率特性表现为组成开环系统的诸典型环节频率特性的合成;
- 系统开环对数频率特性,表现为诸典型环节对数频率特性叠加;
1.3.2 典型环节的频率特性
典型环节频率特性曲线的若干特点:
-
非最小相位与对应的最小相位环节
最小相位的比例环节 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ω- 最小相位惯性环节和非最小相位的惯性环节,幅频特性相同,相频特性符号相反,幅相曲线关于实轴对称;
- 最小相位惯性环节和非最小相位惯性环节,对数幅频曲线相同,对数相频曲线关于0°线对称;
- 上述两点特点,对于振荡环节和非最小相位振荡环节、一阶微分环节和非最小相位一阶微分环节、二阶微分环节和非最小相位二阶微分环节均适用;
-
传递函数互为倒数的典型环节
最小相位典型环节中,积分环节和微分环节、惯性环节和一阶微分环节、振荡环节和二阶微分环节的传递函数互为倒数,有如下关系:
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}^{j\varphi_1(\omega)} G1(jω)=A1(ω)ejφ1(ω),则有:
{ φ 2 ( ω ) = − φ 1 ( ω ) L 2 ( ω ) = 20 lg A 2 ( ω ) = 20 lg 1 A 1 ( ω ) = − L 1 ( ω ) \begin{cases} & \varphi_2(\omega)=-\varphi_1(\omega) \\\\ & L_2(\omega)=20\lg{A_2(\omega)}=20\lg\displaystyle\frac{1}{A_1(\omega)}=-L_1(\omega) \end{cases} ⎩ ⎨ ⎧φ2(ω)=−φ1(ω)L2(ω)=20lgA2(ω)=20lgA1(ω)1=−L1(ω)
传递函数互为倒数的典型环节,对数幅频曲线关于0dB线对称,对数相频曲线关于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 A(\omega)=\frac{1}{\sqrt{\begin{pmatrix}1-\displaystyle\frac{\omega^2}{\omega_n^2}\end{pmatrix}^2+4\zeta^2\displaystyle\frac{\omega^2}{\omega_n^2}}} A(ω)=(1−ωn2ω2)2+4ζ2ωn2ω21φ ( ω ) = − arctan 2 ζ ω ω n 1 − ω 2 ω n 2 = { − arctan 2 ζ ω ω n 1 − ω 2 ω n 2 , ω ≤ ω n − ( 180 − arctan 2 ζ ω ω n ω 2 ω n 2 − 1 ) , ω > ω n \varphi(\omega)=-\arctan\displaystyle\frac{2\zeta\displaystyle\frac{\omega}{\omega_n}}{1-\displaystyle\frac{\omega^2}{\omega_n^2}}= \begin{cases} &-\arctan\frac{2\zeta\displaystyle\displaystyle\frac{\omega}{\omega_n}}{1-\displaystyle\frac{\omega^2}{\omega_n^2}},\omega≤\omega_n \\\\ &-\begin{pmatrix} 180-\arctan\frac{2\zeta\displaystyle\frac{\omega}{\omega_n}}{\displaystyle\frac{\omega^2}{\omega_n^2}-1} \end{pmatrix},\omega>\omega_n \end{cases} φ(ω)=−arctan1−ωn2ω22ζωnω=⎩ ⎨ ⎧−arctan1−ωn2ω22ζωnω,ω≤ωn− 180−arctanωn2ω2−12ζωnω ,ω>ωn
φ ( 0 ) = 0 ° , φ ( ∞ ) = − 180 ° \varphi(0)=0°,\varphi(\infty)=-180° φ(0)=0°,φ(∞)=−180°,相频特性曲线从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 ° \begin{cases} &A(0)=1\\ &\varphi(0)=0° \end{cases}, \begin{cases} &A(\omega_n)=2\zeta\\ &\varphi(\omega_n)=90° \end{cases}, \begin{cases} &A(\infty)=\infty\\ &\varphi(\infty)=180° \end{cases} {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单调增至 ∞ \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单调减至
{ A ( ω r ) = 2 ζ 1 − ζ 2 < 1 ω r = ω 1 − 2 ζ 2 \begin{cases} &A(\omega_r)=2\zeta\sqrt{1-\zeta^2}<1\\ &\omega_r=\omega\sqrt{1-2\zeta^2} \end{cases} {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)= \begin{cases} 0,&\omega<\displaystyle\frac{1}{T}\\ -20\lg\omega{T},&\omega>\displaystyle\frac{1}{T} \end{cases} La(ω)=⎩ ⎨ ⎧0,−20lgωT,ω<T1ω>T1
低频部分是零分贝线,高频部分是斜率为 − 20 d B / d e c -20{\rm dB/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,低频渐近线为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)点,斜率为-40dB/dec的直线,振荡环节的交接频率为 ω n \omega_n ωn,对数幅频渐近特性为:
L a ( ω ) = { 0 , ω < ω n − 40 lg ω ω n , ω > ω n L_a(\omega)= \begin{cases} 0,&\omega<\omega_n\\ -40\lg\displaystyle\frac{\omega}{\omega_n},&\omega>\omega_n \end{cases} La(ω)=⎩ ⎨ ⎧0,−40lgωnω,ω<ωnω>ωn
非最小相位振荡环节与振荡环节对数幅频渐近特性曲线相同,二阶微分环节和非最小相位二阶微分环节与振荡环节的对数幅频渐近特性曲线关于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)为直线斜率;
1.3.3 开环幅相特性曲线
绘制概略开环幅相特性曲线的方法:
-
开环幅相特性曲线的起点 ( ω = 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) -
开环幅相特性曲线的变化范围(象限、单调性);
绘制概略开环幅相特性曲线规律小结:
-
开环幅相特性曲线的起点,取决于比例环节 K K K和系统积分或微分环节的个数 ν \nu ν(系统型别);
- ν < 0 \nu<0 ν<0,起点为原点;
- ν = 0 \nu=0 ν=0,起点为实轴上的点 K K K处, K K K为系统开环增益, K K K有正负之分;
- ν > 0 \nu>0 ν>0,设 ν = 4 k + i ( k = 0 , 1 , 2 , … ; i = 1 , 2 , 3 , 4 ) , \nu=4k+i(k=0,1,2,\dots;i=1,2,3,4), ν=4k+i(k=0,1,2,…;i=1,2,3,4),则 K > 0 K>0 K>0时为 i × ( − 90 ° ) i\times(-90°) i×(−90°)的无穷远处, K < 0 K<0 K<0时为 i × ( − 90 ° ) − 180 ° i\times(-90°)-180° i×(−90°)−180°的无穷远处;
-
开环幅相特性曲线的终点,取决于开环传递函数分子、分母多项式中最小相位环节和非最小相位环节的阶次和;
设系统开环传递函数的分子、分母多项式的阶次分别为 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)= \begin{cases} [(m_1-m_2)-(n_1-n_2)]\times90°,&K>0\\\\ [(m_1-m_2)-(n_1-n_2)]\times90°-180°,&K<0 \end{cases} φ(∞)=⎩ ⎨ ⎧[(m1−m2)−(n1−n2)]×90°,[(m1−m2)−(n1−n2)]×90°−180°,K>0K<0
当开环系统为最小相位系统时,若:
n = m , G ( j ∞ ) H ( j ∞ ) = K ∗ n > m , G ( j ∞ ) H ( j ∞ ) = 0 ∠ [ ( n − m ) × ( − 90 ° ) ] \begin{aligned} &n=m, && G({\rm j}\infty)H({\rm j}\infty)=K^*\\\\ &n>m,&& G{(\rm j}\infty)H({\rm j}\infty)=0\angle[(n-m)\times(-90°)] \end{aligned} n=m,n>m,G(j∞)H(j∞)=K∗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 ° \begin{aligned} &\varphi(\omega_{n^-})≈\varphi_1(\omega_n)=\angle[G_1(j\omega_n)H_1(j\omega_n)]\\\\ &\varphi(\omega_{n^+})≈\varphi_1(\omega_n)-l\times180° \end{aligned} φ(ωn−)≈φ1(ωn)=∠[G1(jωn)H1(jωn)]φ(ωn+)≈φ1(ωn)−l×180°
即 φ ( ω ) \varphi(\omega) φ(ω)在 ω = ω n \omega=\omega_n ω=ωn附近,相角突变 − l × 180 ° -l\times180° −l×180°;
1.3.4 开环对数频率特性曲线
系统开环对数幅频渐近特性:
L
a
(
ω
)
=
∑
i
=
1
N
L
a
i
(
ω
)
L_a(\omega)=\sum_{i=1}^NL_{a_i}(\omega)
La(ω)=i=1∑NLai(ω)
对于任意的开环传递函数,按典型环节分解,将组成系统的典型环节分为三部分:
- K s ν \displaystyle\frac{K}{s^{\nu}} sνK或 − K s ν ( K > 0 ) \displaystyle\frac{-K}{s^{\nu}}(K>0) sν−K(K>0);
- 一阶环节,包括惯性环节、一阶微分环节及对应的非最小相位环节,交接频率为 1 T \displaystyle\frac{1}{T} T1;
- 二阶环节,包括振荡环节、二阶微分环节及对应的非最小相位环节,交接频率为 ω n \omega_n ωn;
记 ω min \omega_{\min} ωmin为最小交接频率,称 ω < ω min \omega<\omega_{\min} ω<ωmin的频率范围为低频段;
开环对数幅频渐近特性曲线绘制步骤:
-
开环传递函数典型环节分解;
-
确定一阶环节、二阶环节的交接频率,将各交接频率标注在半对数坐标图的 ω \omega ω轴上;
-
绘制低频段渐近特性线:在 ω < ω min \omega<\omega_{\min} ω<ωmin频段内,开环系统幅频渐近特性的斜率取决于 K ω ν \displaystyle\frac{K}{\omega^{\nu}} ωνK,因而直线斜率为 − 20 ν d B / d e c -20\nu{\rm dB/dec} −20νdB/dec;获得低频渐近线,需要确定直线上的一点,方法如下:
- 方法1:在 ω < ω min \omega<\omega_{\min} ω<ωmin范围内,任选一点 ω 0 \omega_0 ω0,计算 L a ( ω 0 ) = 20 lg K − 20 ν lg ω 0 L_a(\omega_0)=20\lg{K}-20\nu\lg\omega_0 La(ω0)=20lgK−20νlgω0;
- 方法2:取频率为特定值 ω 0 = 1 \omega_0=1 ω0=1,则 L a ( 1 ) = 20 lg K L_a(1)=20\lg{K} La(1)=20lgK;
- 方法3:取 L a ( ω 0 ) L_a(\omega_0) La(ω0)为特殊值0,有 K ω 0 ν = 1 \displaystyle\frac{K}{\omega_0^{\nu}}=1 ω0νK=1,则 ω 0 = K 1 ν \omega_0=K^{\frac{1}{\nu}} ω0=Kν1;
过点 ( ω 0 , L a ( ω 0 ) ) (\omega_0,L_a(\omega_0)) (ω0,La(ω0))在 ω < ω min \omega<\omega_{\min} ω<ωmin范围内可作斜率为-20 ν \nu νdB/dec的直线;
-
作 ω ≥ ω min \omega≥\omega_{\min} ω≥ωmin频段渐近特性线;在每个交接频率点处,斜率发生变化,变化规律取决于该交接频率对应的典型环节的种类,当系统的多个环节具有相同交接频率时,该交接频率点处斜率的变化应为各个环节对应的斜率变化值的代数和;
交接频率点处斜率的变化表:
1.3.5 延迟环节和延迟系统
输出量经恒定延时后不失真地复现输入量变化的环节称为延迟环节;含有延迟环节的系统称为延迟系统;
延迟环节的输入输出时域表达式为:
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(j\omega)={\rm e}^{-{\rm j}\tau\omega}=1·\angle(-57.3\tau\omega)
G(jω)=e−jτω=1⋅∠(−57.3τω)
1.4 开环和闭环频率特性的性能分析
1.4.1 开环对数频率特性与时域响应的关系
-
低频段分析。
-
低频段:指 L ( ω ) = 20 lg ∣ G ( j ω ) ∣ {\rm L(\omega)=20\lg|G(j\omega)|} L(ω)=20lg∣G(jω)∣的渐近线在第一个转折频率前的频段,这一频段的特性完全由积分环节和开环放大倍数决定;
-
低频对数幅频特性为:
L d ( ω ) = 20 lg K − 20 ν lg ω L_d(\omega)=20\lg{K}-20\nu\lg\omega Ld(ω)=20lgK−20νlgω
其中: K K K为开环放大倍数, ν \nu ν为开环传递函数中积分环节个数; -
低频段的频率特性决定系统的稳态性能;
-
-
中频段分析。
- 中频段:指开环对数幅频特性曲线在开环截止频率 ω c \omega_c ωc附近( 0 d B {\rm 0dB} 0dB附近)的区段,这一频段特性集中反映了闭环系统动态响应的平稳性和快速性;
- 时域响应的动态特性主要取决于中频段的形状;
- 反映中频段形状的三个参数:开环截止频率 ω c \omega_c ωc、中频段斜率、中频段宽度;
- 开环对数幅频特性中频段斜率最好为 − 20 d B / d e c {\rm -20dB/dec} −20dB/dec,且希望长度尽可能长,以确保系统有足够的相角裕度;当中频段的斜率为 − 40 d B / d e c {\rm -40dB/dec} −40dB/dec时,中频段占据的频率范围不宜过长,否则相角裕度很小;若中频段斜率更小,系统则难以稳定;
- 截止频率 ω c \omega_c ωc越高,系统复现信号能力越强,系统快速性越好;
-
高频段分析。
- 高频段:指开环对数幅频特性曲线在中频段以后的频段,高频段的形状主要影响时域响应的起始段;
- 进行系统分析时,可将高频段进行近似处理,即用一个小惯性环节来等效多个小惯性环节,等效的小惯性环节的时间常数等于被替代的多个小惯性环节的时间常数之和;
- 系统开环对数幅频特性在高频段的幅值,直接反映系统对高频干扰信号的抑制能力;高频部分的幅值越低,系统的抗干扰能力越强;
-
良好系统的各频段形状。
- 低频段要有一定的高度和斜率;
- 中频段斜率最好为 − 20 d B / d e c {\rm -20dB/dec} −20dB/dec,且具有足够的宽度;
- 高频段采用迅速衰减的特性,以抑制不必要的高频干扰;
1.4.2 开环频率特性与时域响应的关系
典型二阶系统的频域指标和时域指标相互转换公式:
ω
c
=
1
+
4
ζ
4
−
2
ζ
2
ω
n
,
γ
(
ω
c
)
=
arctan
2
ζ
1
+
4
ζ
4
−
2
ζ
2
\begin{aligned} &\omega_c=\sqrt{\sqrt{1+4\zeta^4}-2\zeta^2}\omega_n,\gamma(\omega_c)=\arctan\displaystyle\frac{2\zeta}{\sqrt{\sqrt{1+4\zeta^4}-2\zeta^2}} \end{aligned}
ωc=1+4ζ4−2ζ2ωn,γ(ωc)=arctan1+4ζ4−2ζ22ζ
高阶系统的频域指标和时域指标相互转换公式:
σ
p
=
0.16
+
0.4
(
1
sin
γ
−
1
)
,
35
°
≤
γ
≤
90
°
t
s
=
π
ω
c
[
2
+
1.5
(
1
sin
γ
−
1
)
+
2.5
(
1
sin
γ
−
1
)
2
]
\begin{aligned} &\sigma_p=0.16+0.4\left(\displaystyle\frac{1}{\sin\gamma}-1\right),35°≤\gamma≤90°\\\\ &t_s=\displaystyle\frac{\pi}{\omega_c}\left[2+1.5\left(\displaystyle\frac{1}{\sin\gamma}-1\right)+2.5\left(\displaystyle\frac{1}{\sin\gamma}-1\right)^2\right] \end{aligned}
σp=0.16+0.4(sinγ1−1),35°≤γ≤90°ts=ωcπ[2+1.5(sinγ1−1)+2.5(sinγ1−1)2]
1.4.3 闭环频率特性与时域响应的关系
典型二阶系统闭环传递函数为:
Φ
(
s
)
=
ω
n
2
s
2
+
2
ζ
ω
n
s
+
ω
n
2
\Phi(s)=\frac{\omega_n^2}{s^2+2\zeta\omega_ns+\omega_n^2}
Φ(s)=s2+2ζωns+ωn2ωn2
可得系统闭环频率特性为:
Φ
(
j
ω
)
=
ω
n
2
(
j
ω
)
2
+
j
2
ζ
ω
n
ω
+
ω
n
2
\Phi({\rm j}\omega)=\frac{\omega_n^2}{({\rm j}\omega)^2+{\rm j}2\zeta\omega_n\omega+\omega_n^2}
Φ(jω)=(jω)2+j2ζωnω+ωn2ωn2
系统闭环幅频特性为:
M
(
ω
)
=
ω
n
2
(
ω
n
2
−
ω
2
)
2
+
(
2
ζ
ω
n
ω
)
2
M(\omega)=\frac{\omega_n^2}{\sqrt{(\omega_n^2-\omega^2)^2+(2\zeta\omega_n\omega)^2}}
M(ω)=(ωn2−ω2)2+(2ζωnω)2ωn2
系统闭环相频特性为:
φ
(
ω
)
=
−
arctan
2
ζ
ω
n
ω
ω
n
2
−
ω
2
\varphi(\omega)=-\arctan\frac{2\zeta\omega_n\omega}{\omega_n^2-\omega^2}
φ(ω)=−arctanωn2−ω22ζωnω
二阶系统的谐振峰值
M
r
M_r
Mr与时域超调量
σ
p
\sigma_p
σp间的关系为:
σ
p
=
e
−
ζ
π
/
1
−
ζ
2
×
100
%
,
M
r
=
1
2
ζ
1
−
ζ
2
\sigma_p={\rm e}^{-\zeta\pi/\sqrt{1-\zeta^2}}\times100\%,M_r=\frac{1}{2\zeta\sqrt{1-\zeta^2}}
σp=e−ζπ/1−ζ2×100%,Mr=2ζ1−ζ21
-
谐振峰值 M r M_r Mr仅与阻尼比 ζ \zeta ζ有关,超调量 σ p \sigma_p σp亦仅取决于阻尼比 ζ \zeta ζ;
-
ζ \zeta ζ越小, M r M_r Mr增加得越快,超调量 σ p \sigma_p σp亦很大,超过 40 % 40\% 40%,此时系统一般不符合瞬态响应指标要求;
-
当 0.4 < ζ < 0.707 0.4<\zeta<0.707 0.4<ζ<0.707时, M r M_r Mr与 σ p \sigma_p σp变化趋势基本一致,此时谐振峰值 M r = 1.2 ~ 1.5 M_r=1.2~1.5 Mr=1.2~1.5,超调量 σ p = 20 % ~ 30 % \sigma_p=20\%~30\% σp=20%~30%,系统响应结果较为理想;
-
当 ζ > 0.707 \zeta>0.707 ζ>0.707时,无谐振峰值, M r M_r Mr与 σ % \sigma\% σ%对应关系不再存在,因此在设计中一般 ζ \zeta ζ取值在 0.4 ~ 0.7 0.4~0.7 0.4~0.7之间;
-
二阶系统谐振频率 ω r \omega_r ωr与峰值时间 t p t_p tp间的关系:
t p ω r = π 1 − 2 ζ 2 1 − ζ 2 t_p\omega_r=\frac{\pi\sqrt{1-2\zeta^2}}{\sqrt{1-\zeta^2}} tpωr=1−ζ2π1−2ζ2 -
当 ζ \zeta ζ为常数时,谐振频率 ω r \omega_r ωr与峰值时间 t p t_p tp成反比, ω r \omega_r ωr值越大, t p t_p tp越小,系统时间响应越快;
-
二阶系统闭环截止频率 ω b \omega_b ωb与过渡过程时间 t s t_s ts间的关系:
ω b t s = 3 ~ 4 ζ 1 − 2 ζ 2 + 2 − 4 ζ 2 + 4 ζ 4 \omega_bt_s=\frac{3~4}{\zeta}\sqrt{1-2\zeta^2+\sqrt{2-4\zeta^2+4\zeta^4}} ωbts=ζ3~41−2ζ2+2−4ζ2+4ζ4 -
当阻尼比 ζ \zeta ζ给定后,闭环截止频率 ω b \omega_b ωb与过渡过程时间 t s t_s ts成反比, ω b \omega_b ωb越大,系统响应速度越快;