系统分析与校正方法——时域法

一、概述

时域法是一种直接在时间域中对系统进行分析和校正的方法。

• 优点：可以提供系统时间响应的全部信息，直观、准确。
• 缺点：研究系统参数改变引起系统性能指标变化的趋势，及对系统进行校正设计时，时域法不是非常方便。

系统的时域性能指标：

1. 动态性能指标：延迟时间，上升时间，峰值时间，调节时间，超调量。
2. 稳态性能指标：稳态误差。

二、一阶系统的时间响应及动态性能

系统的传递函数为：
$$\Phi (s)=\frac{1}{Ts+1}$$

系统单位阶跃响应的拉普拉斯变换为
$$C(s)=\Phi (s)\cdot R(s)=\frac{1}{Ts+1}\cdot \frac{1}{s}=\frac{1}{s}-\frac{1}{s+1/T}$$

单位阶跃响应
$$h(t)={\mathcal{L}}^{-1}[C(s)]=1-e^{-\frac{t}{T}}$$

令
$$h(t_s)=1-e^{-\frac{t_s}{T}}=0.95$$

解得调节时间为
$$t_s=3T$$

三、二阶系统的时间响应及动态性能

系统传递函数标准形式为：
$$\Phi (s)=\frac{{\omega }_n^2}{s^2+2\zeta {\omega }_{n}s+{\omega }_n^2}(首1型)=\frac{1}{T^2{s}^2+2T\zeta s+1}(尾1型)$$

其中，$T=\sqrt{\frac{T_0}{K}},{\omega }_n=\frac{1}{T}=\sqrt{\frac{K}{T_0}},\zeta =\frac{1}{2}\sqrt{\frac{1}{K{T}_0}}$。$\zeta ,{\omega }_n$ 分别称为系统的阻尼比和无阻尼自然频率。

二阶系统闭环特征方程为
$$D(s)=s^2+2\zeta {\omega }_{n}s+{\omega }_n^2=0$$

其特征根为
$${\lambda }_{1,2}=-\zeta {\omega }_n\pm {\omega }_n\sqrt{{\zeta }^2-1}$$

当 $\zeta >1$ 系统为过阻尼，当 $\zeta =1$ 系统为临界阻尼，当 $0<\zeta <1$ 系统为欠阻尼，当 $\zeta =0$ 系统为零阻尼。

3.1 过阻尼二阶系统动态性能指标计算

设过阻尼二阶系统的极点为
$${\lambda }_1=-\frac{1}{T_1}=-(\zeta -\sqrt{{\zeta }^2-1}){\omega }_n{\lambda }_2=-\frac{1}{T_2}=-(\zeta +\sqrt{{\zeta }^2-1}){\omega }_n(T_1>T_2)$$

系统单位阶跃响应的拉普拉斯变换
$$C(s)=\Phi (s)R(s)=\frac{{\omega }_n^2}{(s+1/T_1)(s+1/T_2)}\frac{1}{s}$$

进行拉普拉斯反变换，得出系统单位阶跃响应
$$h(t)=1+\frac{e^{-\frac{t}{T_1}}}{\frac{T_2}{T_1}-1}+\frac{e^{-\frac{t}{T_2}}}{\frac{T_1}{T_2}-1}(t\ge 0)$$

当 $T_1/T_2$ （或 $\zeta$）很大时，特征根 ${\lambda }_2=-1/T_2$ 比 ${\lambda }_1=-1/T_1$ 远离虚轴，模态 $e^{-t/T_2}$ 很快衰减为零，系统调节时间主要由 ${\lambda }_1=-1/T_1$ 对应的模态 $e^{-t/T_1}$ 决定。此时可将过阻尼二阶系统近似看作由 ${\lambda }_1$ 确定的一阶系统，估算其动态性能指标。

3.2 欠阻尼二阶系统动态性能指标计算

系统单位阶跃响应为
$$h(t)=1-\frac{e^{-\zeta {\omega }_{n}t}}{\sqrt{1-{\zeta }^2}}\sin (\sqrt{1-{\zeta }^2}{\omega }_{n}t+\arctan \frac{\sqrt{1-{\zeta }^2}}{\zeta })$$

• 峰值时间 $t_p=\frac{\pi }{\sqrt{1-{\zeta }^2}{\omega }_n}$
• 超调量 $\sigma \%=e^{-\zeta \pi \sqrt{1-{\zeta }^2}}\times 100\%$
• 调节时间 $t_s\approx \frac{3.5}{\zeta {\omega }_n}(0.3<\zeta <0.8)$

典型欠阻尼二阶系统动态性能、系统参数及极点分布之间的关系
$$\{\begin{array}{ll}\sigma \%& =e^{-\zeta \pi \sqrt{1-{\zeta }^2}}\times 100\%\\ t_s& \approx \frac{3.5}{\zeta {\omega }_n}\\ {\lambda }_{1,2}& =\sigma \pm j{\omega }_d=-\zeta {\omega }_n\pm j\sqrt{1-{\zeta }^2}{\omega }_n\end{array}$$

1. 当 ${\omega }_n$ 固定，$\zeta$ 增加（$\beta$ 减小），系统超调量减小；同时由于极点远离虚轴，调节时间减小。
2. 当 $\zeta$ 固定，${\omega }_n$ 增加，系统超调量不变；由于极点远离虚轴，调节时间减小。
3. 当 $T_0$ 固定，$K$ 增加，阻尼比变小，超调量增加。
   
   综上所述，要获得满意的系统动态性能，应该适当选择参数，使二阶系统的闭环极点位于 $\beta =45^{\circ }$ 线附近，使系统具有合适的超调量，并根据情况尽量使其远离虚轴，以提高系统的快速性。

3.3 改善二阶系统动态性能的措施

采用测速反馈和**比例加微分（PD）**控制方式，可以有效改善二阶系统的动态性能。

系统结构图$(K_t=0.216)$	图$(a)$原系统	图$(b)$测速反馈	图$(c)$比例+微分
开环传递函数	$G_{(a)}(s)=\frac{10}{s(s+1)}$	$G_{(b)}(s)=\frac{10}{s(s+1+10K_t)}$	$G_{(c)}(s)=\frac{10(K_{t}s+1)}{s(s+1)}$
开环增益	$K_{(a)}=10$	$K_{(b)}=\frac{10}{1+10K_t}$	$K_{(c)}=10$
闭环传递函数	${\Phi }_{(a)}(s)=\frac{10}{s^2+s+10}$	${\Phi }_{(b)}(s)=\frac{10}{s^2+(1+10K_t)s+10}$	${\Phi }_{(c)}(s)=\frac{10(K_{t}s+1)}{s^2+(1+10K_t)s+10}$
$\zeta$	0.158	0.5	0.5
${\omega }_n$	3.16	3.16	3.16
零点	—	—	$-4.63$
极点	$-0.5\pm j3.12$	$-1.58\pm j2.74$	$-1.58\pm j2.74$
$t_p$	1.01	1.15	1.05
$\sigma \%$	60%	16.3%	23%
$t_s$	7	2.2	2.1

由于引入了测速反馈和PD控制，$(b)$ 和 $(c)$ 所示系统的闭环极点较 $(a)$ 所示系统闭环极点远离虚轴（相应调节时间小）；且 $\beta$ 减小（对应阻尼比较大，超调量较小）。因而动态性能优于图 $(a)$ 所示系统。

比较 $(b)$ 和 $(c)$ 所示两系统，它们闭环传递函数的分母相同，只是后者较前者多一个闭环零点。附件闭环零点不会影响闭环极点，因而不会影响单位阶跃响应中的各模态。但它会改变单位阶跃响应中各模态的加权系数，由此影响系统的动态性能。

四、高阶系统的阶跃响应及动态性能

对稳定的闭环系统，远离虚轴的极点对应的模态因为收敛速度快，只影响阶跃响应的起始段，而距虚轴近的极点对应的模态衰减缓慢，系统动态性能主要取决于这些极点对应地响应分量。

距离虚轴最近而且附近又没有零点的极点对系统的动态性能起主导作用，称相应极点为主导极点。

一般规定，若某极点的实部大于主导极点实部的 5~6 倍以上时，则可以忽略相应分量的影响；若两相邻零、极点间的距离比它们本身的模值小一个数量级时，则称该零、极点为“偶极子”，其作用近似抵消，可以忽略相应分量的影响。

在绝大多数实际系统的闭环零、极点中，可以选留最靠近虚轴的一个或几个极点作为主导极点，略去比主导极点距虚轴远 5 倍以上的闭环零、极点，以及不十分接近虚轴的相互靠得很近（该零、极点距虚轴距离是其相互之间距离的10倍以上）的偶极子，忽略其对系统动态性能的影响。