时域输入信号python

第三章 线性系统的时域分析法

接下来的三个章节将分别介绍线性系统的时域分析法、根轨迹分析法、频域分析法，这三个部分的结构十分相似，每个部分基本都由三个板块构成：分析系统的稳定性（判稳）、系统瞬态性能和稳态性能（性能）、改善性能的方法（设计）。

本章讨论时域分析，就是控制系统在一定的输入信号作用下，根据系统输出量的时域表达式，分析系统的稳定性、瞬态和稳态性能。时域分析的特点就是直观准确，表达式是关于时间$t$的函数，所以也称作系统的时间响应。

一、典型的输入信号

名称	信号定义	拉氏变换
脉冲信号	$\delta(t)=\bigg\{ \begin{matrix} 0,\ t\not= 0 \\ \infty,\ t=0 \end{matrix}$$\int_{-\infty}^{+\infty}\delta(t)dt=1$	$L[\delta(t)]=1$
阶跃信号	$r_s(t)=\bigg\{ \begin{matrix} 0,\ t< 0 \\ A·1(t),\ t\geq0 \end{matrix}$	$L[1(t)]=\frac{1}{s},\ A=1$
斜坡信号（速度阶跃）	$r_v(t)=\bigg\{ \begin{matrix} 0,\ t< 0 \\ At,\ t\geq0 \end{matrix}$	$L[t]=\frac{1}{s^2},\ A=1$
抛物线信号（加速度阶跃）	$r_a(t)=\bigg\{ \begin{matrix} 0,\ t< 0 \\ \frac 1 2 At^2,\ t\geq0 \end{matrix}$	$L[\frac{1}{2}t^2]=\frac{1}{s^3},\ A=1$
正弦信号	$r(t)=\bigg\{ \begin{matrix} 0,\ t< 0 \\ Asin\omega t,\ t\geq0 \end{matrix}$	$L[sin\omega t]=\frac{\omega}{s^2+\omega^2},\ A=1$

二、瞬态响应、稳态响应与性能指标

2.1 瞬态响应

又称瞬态过程或过渡过程，是系统在典型输入信号作用下，输出量从初始状态到最终状态的响应过程。瞬态响应过程曲线可能表现为衰减振荡、等幅振荡或发散等形式，一个稳定运行的控制系统的瞬态过程必须是衰减的。

瞬态性能指标

1. 上升时间$t_r$：单位阶跃响应第一次达到稳态值的时间。上升时间越短，响应速度越快。
2. 延迟时间$t_d$：单位阶跃响应第一次达到稳态值50%的时间。
3. 峰值时间$t_p$：单位阶跃响应到达第一个峰值的时间
4. 最大超调量$\delta \%$：单位阶跃响应的最大值$y_{max}=y(t_p)$与稳态值$y(\infty)$之差与稳态值的之比的百分数，即：$\delta \%=\frac{y(t_p)-y(\infty)}{y(\infty)}\times 100\%$单调上升的阶跃响应没有超调量和峰值时间。
5. 调整时间$t_s$：当$y(t)$与$y(\infty)$的误差绝对值小于等于规定允许值，且以后不再超过次值所需要的最短时间。即：$|y(t)-y(\infty)|\leq y(\infty)\times \Delta \%,\ t\geq t_s$
6. 振荡次数$N$：从开始到调节时间内，单位阶跃响应穿越稳态值次数的一半。

2.2 稳态响应

又称稳态过程，是系统在典型输入信号作用下，当时间趋近于无穷大时，系统的输出响应状态。工程上把瞬态响应在调节时间以后的响应过程视为稳态过程。

稳态性能指标

主要是稳态误差，详见下方介绍。

三、典型一阶系统

3.1 数学模型

一阶系统微分方程为：
$T\frac{dy(t)}{dt}+y(t)=r(t)$
零初始条件下，一阶系统传递函数为：
$\Phi(s)=\frac{Y(s)}{R(s)}=\frac{1}{Ts+1}$

3.2 一阶系统对典型输入信号的响应（上一个是下一个的导数关系）

输入信号	输出响应
$N$

3.3 一阶系统的瞬态性能指标（针对单位阶跃响应）

1. 延迟时间：$t_d \approx 0.693T$
2. 上升时间：$t_r\approx 2.197T$
3. 调整时间：$t_s=-Tln\Delta\%\approx\bigg\{\begin{matrix}4T,\ \Delta=2\\ 3T,\ \Delta=5 \end{matrix}$

3.4 改善一阶系统瞬态性能（减小时间常数）

方法一：通过负反馈减小时间常数

$\Phi(s)=\frac{\frac{1}{\alpha+1}}{\frac{T}{\alpha+1}s+1}$

方法二：通过增加开环放大系数减小时间常数

$\Phi(s)=\frac{1}{\frac{1}{\alpha K}s+1}=\frac{1}{\frac{T}{\alpha }s+1}$

四、典型二阶系统的瞬态性能

4.1 数学模型

二阶系统微分方程为：
$T^2\frac{d^2d(t)}{dt^2}+2\zeta T\frac{dy(t)}{dt}+y(t)=r(t),\ t\geq0$
零初始条件下，二阶系统传递函数为：
$\Phi(s)=\frac{Y(s)}{R(s)}=\frac{1}{T^2s^2+2\zeta Ts+1}=\frac{\omega_n^2}{s^2+2\zeta \omega_n s+\omega_n^2}$

4.2 二阶系统的单位阶跃响应

根据传递函数分母多项式（特征方程）解的情况分类：

阻尼系数	特征方程根	特征方程根的位置	单位阶跃响应的形式
无阻尼		虚轴上一对共轭虚根	等幅周期振荡
欠阻尼		s左半平面一对共轭复根	衰减振荡
临界阻尼		负实轴上一对重根	单调上升
过阻尼		负实轴上两个互异根	单调上升

4.3 二阶系统的瞬态性能指标

欠阻尼典型二阶系统的瞬态性能指标

1. 上升时间：$t_r=\frac{\pi-arccos\zeta}{\omega_n\sqrt{1-\zeta^2}}=\frac{\pi-\beta(阻尼角)}{\omega_d(阻尼振荡频率)}$
2. 峰值时间：$t_p=\frac{\pi}{\omega_n\sqrt{1-\zeta^2}}=\frac{\pi}{\omega_d}$
3. 超调量：$\delta\%=e^{-\frac{\zeta\pi}{\sqrt{1-\zeta^2}}}\times100\%$
4. 调整时间：$t_s\approx\bigg\{\begin{matrix}\frac{4}{\zeta\omega_n},\ \Delta=2\\ \frac{3}{\zeta\omega_n},\ \Delta=5\end{matrix}$
5. 振荡次数：$N=\frac{t_s}{2\pi/\omega_d}$

过阻尼（临界阻尼）二阶系统的瞬态性能指标

1. 上升时间：$t_r=\frac{1+1.5\zeta+\zeta^2}{\omega_n}$
2. 调整时间：$t_s\approx\bigg\{\begin{matrix}8.4/\omega_n, \Delta=2\\ 6.6/\omega_n, \Delta=5\end{matrix}$

4.4 二阶系统瞬态性能改善

比例微分校正

传递函数为：$\Phi(s)=\frac{\omega_n^2(k_p+k_d s)}{s^2+(2\zeta\omega_n+\omega_n^2k_d)s+\omega_n^2k_p}=\frac 1 z ·\frac{\omega_{kd}^2(s+z)}{s^2+2\zeta_{kd}\omega_{kd}s+\omega_{kd}^2}$

其中，$z=\frac{k_p}{k_d},\ \omega_{kd}=\sqrt{k_p\omega_n},\ \zeta_{kd}=(\zeta+\frac{\omega_nk_d}{2})/\sqrt{k_p}$

速度反馈校正

传递函数为：$\Phi(s)=\frac{\omega_n^2}{s^2+(2\zeta\omega_n+\tau\omega_n^2)s+\omega_n^2}$

五、高阶系统的时域分析

跳过讨论，直接上结论：

1. 若某极点远离原点，其相应的瞬态响应分量的系数很小
2. 若某极点接近一个零点，而远离其他零极点，则相应的瞬态响应分量的系数也很小
3. 若某极点远离零点而又接近原点或其他极点，则相应的瞬态响应分量的系数比较大

主导极点

满足以下情况的极点称为主导极点：

1. 离虚轴最近且周围没有零点
2. 其他极点与虚轴的距离比该极点与虚轴的距离的五倍还要远

保留主导极点，略去其他极点，可以化简系统。但是要注意简化系统的稳态值要与原系统的稳态值一致

六、线性控制系统的稳定性分析

6.1 线性控制系统渐进稳定的充分必要条件

系统的所有特征根必须位于s平面的左半开平面。

6.2 代数稳定性判据（只讨论劳斯判据）

假设线性控制系统的特征方程为：$a_ns^n+a_{n-1}s^{n-1}+\dots+a_1s+a_0=0$
劳斯阵列定义为：$\begin{array}{c|llll} {s^n}&{a_n}&{a_{n-2}}&{a_{n-4}}&{\dots}\\{s^{n-1}}&{a_{n-1}}&{a_{n-3}}&{a_{n-5}}&{\dots}\\{s^{n-2}}&{b_1}&{b_2}&{b_3}&{\dots}\\{s^{n-3}}&{c_1}&{c_2}&{c_3}&{\dots}\\ {\vdots}&{\vdots}\\ {s^0}&{x_1}\end{array}$
其中，$b_i=\frac{-1}{a_{n-1}}\begin{vmatrix}a_n&a_{n-2i}\\ a_{n-1} &a_{n-1-2i} \end{vmatrix},\ i=1,2,\dots$
下方各行都按照这种方式生成相应位置的元素。
系统特征方程具有正实部根的数目与劳斯阵列第一列元素符号变化的次数相等。因此线性系统稳定的充要条件是：劳斯阵第一列元素没有符号变化。

特殊情况：

1. 某一行的第一列元素为0，用一个小正数$\epsilon$代替，根据此数继续计算，若它与其上面或下面元素符号相反，则记一次符号变化。
2. 某一行元素全为0，说明系统的特征方程存在着大小相等而径向位置相反的根，至少存在下述几种特征根之一：存在大小相等、符号相反的一对实根；或共轭虚根；或对称于虚根的两对共轭复根。此时可以用全零行的上一行元素构造一个辅助方程，并将该辅助方程对复变量s求导，用求导后方程的系数取代全零行元素，继续构建劳斯阵。辅助方程的根一般就是共轭复根。

七、线性控制系统的稳态性能分析

7.1 系统型别

系统的开环传函可表示为：$G_k(s)=\frac{K}{s^{\nu}}·\frac{\prod\limits_{i=1}^{m_1}(\tau_is+1)\prod\limits_{k=1}^{m_2}(\tau_k^2s^2+2\zeta_k\tau_ks+1)}{\prod\limits_{j=1}^{n_1}(T_js+1)\prod\limits_{l=1}^{n_2}(T_l^2s^2+2\zeta_lT_ls+1)}=\frac{K}{s^{\nu}}·G_0(s)$
满足：$m_1+2m_2=m,\ v+n_1+2n_2=n$
将$\nu$的取值，也就是积分环节的个数定义为系统型别（0型、Ⅰ型、Ⅱ型等等）。

7.2 系统稳态误差

假设被研究的控制系统是稳定的，得到控制系统的给定稳定误差为：
$e_{ssr}=\lim\limits_{s\to 0}sE(s)=\lim\limits_{s\to 0}\frac{sR(s)}{1+G_k(s)}=\lim\limits_{s\to 0}\frac{sR(s)}{1+\frac{K}{s^{\nu}}·G_0(s)}$

典型参考输入下系统的给定稳定误差和静态误差系数

7.3 减小或消除稳态误差的措施

（1）比例积分控制

（2）复合控制

顺馈控制系统

前馈控制系统