时域响应 Time Response

时域响应是评判系统的第一现场：直接看输出随时间如何演化。一阶、二阶系统提供标尺，上升时间、超调量、整定时间、稳态误差四项指标把”好不好”量化成可检验的数，而全部动态品质由闭环极点在 s 平面的位置唯一决定。

一阶系统

传递函数 Transfer Function最简形态：

$G(s)=\frac{1}{Ts+1}⟹y(t)=1-e^{-t/T}$

一切由时间常数 $T$ 决定：$t=T$ 时输出达 63.2%，$t=3T$ 达 95%，$t=4T$ 达 98%。响应单调无超调，$t_r\approx 2.2T$（10%→90%），$t_s\approx 3T$（5% 误差带）。想让系统快 $n$ 倍，唯一办法是把极点 $-1/T$ 向左推 $n$ 倍。

二阶系统

标准形式：

$G(s)=\frac{{\omega }_n^2}{s^2+2\zeta {\omega }_{n}s+{\omega }_n^2}$

特征方程 $s^2+2\zeta {\omega }_{n}s+{\omega }_n^2=0$ 就是常微分方程的特征方程问题，两个参数 $\zeta$（阻尼比）与 ${\omega }_n$（自然频率）刻画全部行为：

阻尼比	极点	响应形态
$\zeta =0$	纯虚数 $\pm j{\omega }_n$	等幅振荡
$0<\zeta <1$	共轭复根 $-\zeta {\omega }_n\pm j{\omega }_d$	欠阻尼：衰减振荡
$\zeta =1$	二重实根 $-{\omega }_n$	临界阻尼：无超调最快
$\zeta >1$	两个负实根	过阻尼：无超调但迟缓

欠阻尼（工程主流）阶跃响应，其中 ${\omega }_d={\omega }_n\sqrt{1-{\zeta }^2}$：

$y(t)=1-\frac{e^{-\zeta {\omega }_{n}t}}{\sqrt{1-{\zeta }^2}}\sin \left({\omega }_{d}t+\phi \right),\phi =\arctan \frac{\sqrt{1-{\zeta }^2}}{\zeta }$

动态性能指标

指标	公式	由谁决定
上升时间（0→100%）	$t_r=\frac{\pi -\phi }{{\omega }_d}$	$\zeta$ 与 ${\omega }_n$ 共同
峰值时间	$t_p=\frac{\pi }{{\omega }_d}$	振荡频率
超调量	$\sigma \%=e^{-\pi \zeta /\sqrt{1-{\zeta }^2}}\times 100\%$	只由 $\zeta$ 决定
整定时间（进入 ±2% 带）	$t_s\approx \frac{4}{\zeta {\omega }_n}$	只由极点实部决定

关键读法：$\sigma \%$ 是 $\zeta$ 的单调减函数——$\zeta =0.5$ 超调 16.3%，$\zeta =0.707$ 超调 4.3%，后者即经典设计目标；$t_s$ 由衰减包络 $e^{-\zeta {\omega }_{n}t}$ 决定，与振荡频率几乎无关。

import numpy as np
import matplotlib.pyplot as plt
from scipy import signal
 
wn = 1.0  # 固定自然频率，只变阻尼比
t = np.linspace(0, 12, 600)
for zeta in [0.2, 0.5, 0.707, 1.0, 2.0]:
    G = signal.TransferFunction([wn**2], [1, 2*zeta*wn, wn**2])
    _, y = G.step(T=t)
    plt.plot(t, y, label=f"zeta={zeta}")
plt.axhline(1.0, ls="--", lw=0.8)  # 设定值参考线
plt.legend(); plt.xlabel("t / s"); plt.title("二阶系统阶跃响应")
plt.show()

极点几何：s 平面即品质地图

复根 $s=-\zeta {\omega }_n\pm j{\omega }_d$ 的几何坐标直接编码动态品质：

• 实部 $-\zeta {\omega }_n$：衰减速度，决定 $t_s$；越左收敛越快
• 虚部 ${\omega }_d$：振荡频率，决定 $t_p$ 与振铃周期
• 模长 ${\omega }_n$ 与负实轴夹角 $\theta$（$\cos \theta =\zeta$）：从原点出发的射线是等阻尼比线，同心圆是等 ${\omega }_n$ 线

设计指标因此翻译成 s 平面上的禁区：超调上限 $\Rightarrow$ 等阻尼比射线左侧；整定时间上限 $\Rightarrow$ 垂直线 $\sigma \le -4/t_s$ 以左。极点落在左半平面才收敛，这是系统稳定性 Stability的核心判据；全平面所有极点离虚轴的距离还给出收敛快慢的排序。

高阶系统降阶：距虚轴最近、附近无零点的共轭极点对主导整个响应（主导极点）；当其余极点实部比它远 5 倍以上，高阶系统近似为该极点对的二阶系统。互相靠近的零极点对（偶极子）则可对消忽略。这一近似是整个频域、时域设计方法通向现实的桥梁。

稳态误差与系统型别

动态指标管过渡过程，稳态误差 $e_{ss}$ 管终值精度。单位负反馈下由终值定理：

$e_{ss}={\lim }_{s\to 0}\frac{sR(s)}{1+G(s)}$

开环传函含积分环节的个数（系统型别）决定能无差跟踪的输入阶次：

型别	阶跃	斜坡	抛物线
Type 0	$\frac{1}{1+K_p}$	$\infty$	$\infty$
Type I	$0$	$\frac{1}{K_v}$	$\infty$
Type II	$0$	$0$	$\frac{1}{K_a}$

静态误差系数 $K_p,K_v,K_a$ 是 $G(s)$ 在原点的极限：积分环节数给无差度，开环增益给误差大小——这就是 PID 控制 积分项能清稳态误差的传函解释。

阻尼直觉的跨界迁移

• 动画与物理引擎：spring-damper 模型 $\ddot{x}+2\zeta {\omega }_n\dot{x}+{\omega }_n^{2}x=f(t)$ 就是二阶系统。iOS/CSS spring 动画的 stiffness/damping 参数即 ${\omega }_n^2$ 与 $2\zeta {\omega }_n$；UI 默认临界阻尼（最快无过冲），游戏摄像机跟拍故意用欠阻尼制造”呼吸感”
• 优化器：梯度下降加动量后连续时间极限是欠阻尼二阶系统——动量压低有效阻尼，loss 曲线振荡 = $\zeta$ 偏小的读数；学习率扮演增益角色，增益过大沿根轨迹 Root Locus把闭环极点推向虚轴甚至右半平面，对应训练震荡直至发散
• 调度与 RL：整定时间对应收敛的样本效率，超调对应越过安全边界的风险，性能指标表可直接改写成 reward 设计的工程语言

一句话收拢：控制器设计（PID 整定、根轨迹、极点配置）无论工具多花哨，本质都是在 s 平面上把闭环主导极点搬进指标允许的区域——快与稳的张力永远由 $\zeta$ 与 ${\omega }_n$ 的取舍来仲裁。