频域传递函数构建与稳定性裕度分析

一、传递函数的推导

1.1 状态空间的拉普拉斯变换

由状态空间方程：

$$\dot{\mathbf{x}}=\mathbf{Ax}+\mathbf{B}\delta$$

取拉普拉斯变换（零初始条件）：

$$s\mathbf{X}(s)=\mathbf{AX}(s)+\mathbf{B}\Delta (s)$$

整理得：

$$(s\mathbf{I}-\mathbf{A})\mathbf{X}(s)=\mathbf{B}\Delta (s)$$

1.2 传递函数矩阵

解得：

$$\mathbf{X}(s)=(s\mathbf{I}-\mathbf{A}{)}^{-1}\mathbf{B}\Delta (s)$$

定义传递函数矩阵：

$$\mathbf{G}(s)=(s\mathbf{I}-\mathbf{A}{)}^{-1}\mathbf{B}=\left[\begin{array}{c}{G}_{v_y}(s)\\ G_r(s)\end{array}\right]$$

1.3 逆矩阵展开

对于 $2\times 2$ 系统：

$$(s\mathbf{I}-\mathbf{A}{)}^{-1}=\frac{1}{\det (s\mathbf{I}-\mathbf{A})}\left[\begin{array}{cc}s-A_{22}& -A_{12}\\ -A_{21}& s-A_{11}\end{array}\right]$$

其中特征多项式：

$$\det (s\mathbf{I}-\mathbf{A})=s^2-(A_{11}+A_{22})s+(A_{11}{A}_{22}-A_{12}{A}_{21})$$

1.4 横摆角速度传递函数

$$G_r(s)=\frac{B_{2}s+(A_{12}{B}_2-A_{22}{B}_1)}{s^2-(A_{11}+A_{22})s+(A_{11}{A}_{22}-A_{12}{A}_{21})}$$

代入状态空间矩阵元素，经整理得：

$$G_r(s)=\frac{b_{1}s+b_0}{s^2+2\zeta {\omega }_{n}s+{\omega }_n^2}$$

其中：

$$\begin{array}{rl}{b}_1& =\frac{C_{\alpha f}{L}_f}{I_z}\\ b_0& =\frac{C_{\alpha f}(C_{\alpha f}{L}_f^2+C_{\alpha r}{L}_r^2)}{m{I}_z{v}_x^2}-\frac{C_{\alpha f}{L}_f(C_{\alpha f}+C_{\alpha r})}{m{I}_z{v}_x}\end{array}$$

1.5 侧向加速度传递函数

侧向加速度 $a_y={\dot{v}}_y+v_{x}r$ 的传递函数：

$$G_{a_y}(s)=s{G}_{v_y}(s)+v_x{G}_r(s)$$

---

二、频率响应特性

2.1 频率响应的定义

令 $s=j\omega$，得到频率响应：

$$G_r(j\omega )=\frac{b_{1}j\omega +b_0}{-{\omega }^2+2\zeta {\omega }_{n}j\omega +{\omega }_n^2}$$

2.2 幅频特性

$$|G_r(j\omega )|=\frac{\sqrt{b_0^2+(b_1\omega {)}^2}}{\sqrt{({\omega }_n^2-{\omega }^2{)}^2+(2\zeta {\omega }_n\omega {)}^2}}$$

2.3 相频特性

$$\angle G_r(j\omega )=\arctan \left(\frac{b_1\omega }{b_0}\right)-\arctan \left(\frac{2\zeta {\omega }_n\omega }{{\omega }_n^2-{\omega }^2}\right)$$

2.4 典型频率响应曲线

graph TB
    subgraph 幅频特性
        A[低频区] --> |"ω << ω_n"| B[稳态增益 G_r(0)]
        C[共振区] --> |"ω ≈ ω_n"| D["峰值 Mp = 1/(2ζ√(1-ζ²))"]
        E[高频区] --> |"ω >> ω_n"| F["-40 dB/dec 衰减"]
    end
    
    subgraph 相频特性
        G[低频] --> |"ω → 0"| H["相位 → 0°"]
        I[共振] --> |"ω = ω_n"| J["相位 = -90°"]
        K[高频] --> |"ω → ∞"| L["相位 → -180°"]
    end

---

三、稳定性裕度分析

3.1 稳定性的频域判据

线性系统稳定的充要条件：

• 所有极点位于复平面左半平面
• 等价于：特征多项式所有系数为正（二阶系统）

3.2 特征多项式的系数

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

稳定性条件：

$$\begin{array}{rl}2\zeta {\omega }_n& >0\\ {\omega }_n^2& >0\end{array}$$

3.3 不足转向梯度与稳定性

由 ${\omega }_n^2$ 的表达式：

$${\omega }_n^2=\frac{L^2{C}_{\alpha f}{C}_{\alpha r}}{m{I}_z{v}_x^2}(1+K{v}_x^2)$$

转向特性	$K$ 符号	${\omega }_n^2$	稳定性
不足转向	$K>0$	恒为正	所有车速稳定
中性转向	$K=0$	恒为正	临界稳定
过多转向	$K<0$	$v_x<v_{\text{crit}}$ 时为正	条件稳定

3.4 阻尼比与稳定性裕度

阻尼比 $\zeta$ 反映系统的相对稳定性：

$\zeta$ 范围	稳定性裕度	响应特性
$\zeta >0.707$	充足	无共振峰，响应平稳
$0.5<\zeta <0.707$	良好	小超调，快速响应
$0.3<\zeta <0.5$	不足	明显振荡
$\zeta <0.3$	临界	强烈振荡
$\zeta <0$	不稳定	发散

3.5 相位裕度与增益裕度

对于开环传递函数 $G_{\text{open}}(s)$：

相位裕度（Phase Margin）：

$$\text{PM}=180^{\circ }+\angle G_{\text{open}}(j{\omega }_c)$$

其中 ${\omega }_c$ 为增益穿越频率（$|G_{\text{open}}(j{\omega }_c)|=1$）。

增益裕度（Gain Margin）：

$$\text{GM}=\frac{1}{|G_{\text{open}}(j{\omega }_{180})|}$$

其中 ${\omega }_{180}$ 为相位穿越频率（$\angle G_{\text{open}}(j{\omega }_{180})=-180^{\circ }$）。

稳定裕度要求

典型设计要求：

• 相位裕度：$\text{PM}>45^{\circ }$
• 增益裕度：$\text{GM}>6\mathrm{dB}$

---

四、参数敏感性分析

4.1 车速对稳定性的影响

参数	低速 ($v_x\to 0$)	高速 ($v_x\to \infty$)
${\omega }_n$	$\to \infty$	$\to 0$（过多转向）
$\zeta$	$\to \infty$（过阻尼）	减小
$G_r(0)$	$\to 0$	$\to 0$（不足转向）

4.2 侧偏刚度的影响

前轮刚度 $C_{\alpha f}$ 增大：

$$\begin{array}{rl}{\omega }_n& \uparrow \\ \zeta & \uparrow \\ K& \downarrow \text{（趋向过多转向）}\end{array}$$

后轮刚度 $C_{\alpha r}$ 增大：

$$\begin{array}{rl}{\omega }_n& \uparrow \\ \zeta & \uparrow \\ K& \uparrow \text{（趋向不足转向）}\end{array}$$

4.3 质心位置的影响

质心前移（$L_f$ 减小，$L_r$ 增大）：

$$\begin{array}{rl}K& \uparrow \text{（更不足转向）}\\ {\omega }_n& \text{变化复杂}\end{array}$$

质心后移（$L_f$ 增大，$L_r$ 减小）：

$$\begin{array}{rl}K& \downarrow \text{（更过多转向）}\\ {\omega }_n& \text{变化复杂}\end{array}$$

---

五、波特图与奈奎斯特图

5.1 波特图的绘制

波特图包含幅频特性和相频特性两条曲线：

幅频特性波特图：

• 横轴：${\log }_{10}(\omega )$
• 纵轴：$20{\log }_{10}|G_r(j\omega )|(\mathrm{dB})$

相频特性波特图：

• 横轴：${\log }_{10}(\omega )$
• 纵轴：$\angle G_r(j\omega ){(}^{\circ })$

5.2 奈奎斯特图的绘制

奈奎斯特图为复平面上的轨迹：

$$\text{Re}[G_r(j\omega )]\text{vs.}\text{Im}[G_r(j\omega )]$$

参数为频率 $\omega$。

5.3 稳定性判据的奈奎斯特表述

奈奎斯特稳定判据：

开环传递函数 $G_{\text{open}}(s)$ 的奈奎斯特图：

• 不包围 $(-1,j0)$ 点：系统稳定
• 包围 $(-1,j0)$ 点：系统不稳定
• 穿越 $(-1,j0)$ 点：临界稳定

---

六、数值示例

6.1 车辆参数

参数	符号	数值
整车质量	$m$	$1500\mathrm{kg}$
轴距	$L$	$2.7\mathrm{m}$
质心到前轴	$L_f$	$1.1\mathrm{m}$
质心到后轴	$L_r$	$1.6\mathrm{m}$
前轮侧偏刚度	$C_{\alpha f}$	$60000N/rad$
后轮侧偏刚度	$C_{\alpha r}$	$55000N/rad$
横摆转动惯量	$I_z$	$2500\mathrm{kg}\cdot {\mathrm{m}}^2$

6.2 不同车速下的频域参数

车速	${\omega }_n(rad/s)$	$\zeta$	$G_r(0)({\mathrm{s}}^{-1})$
$40km/h$	5.8	0.82	0.52
$80km/h$	4.2	0.65	0.92
$120km/h$	3.1	0.55	1.15

6.3 共振峰分析

在 $v_x=80km/h$ 时：

$$\begin{array}{rl}{\omega }_n& =4.2rad/s\approx 0.67\mathrm{Hz}\\ \zeta & =0.65\\ M_p& =\frac{1}{2\zeta \sqrt{1-{\zeta }^2}}\approx 1.02\end{array}$$

共振峰仅为稳态增益的 2%，系统响应平稳。

6.4 相位裕度计算

对于单位反馈系统，等效开环传递函数：

$$G_{\text{open}}(s)=\frac{G_r(s)}{1-G_r(s)}$$

在 $v_x=80km/h$ 时：

• 相位裕度：$\text{PM}\approx 65^{\circ }$
• 增益裕度：$\text{GM}\approx 12\mathrm{dB}$

---

七、控制器设计启示

7.1 电子稳定程序（ESP）的设计目标

ESP 通过制动干预调节横摆角速度，等效增加阻尼：

$${\zeta }_{\text{eff}}=\zeta +\Delta {\zeta }_{\text{ESP}}$$

典型设计：

• 目标阻尼比：${\zeta }_{\text{eff}}>0.7$
• 响应延迟：$<30\mathrm{ms}$

7.2 主动后轮转向的稳定性增强

主动后轮转向可等效调节不足转向梯度：

$$K_{\text{eff}}=K+\Delta K_{\text{ARS}}$$

高速时可通过同向偏转后轮减小 $K$，提高响应性。

7.3 扭矩矢量分配的横摆控制

独立驱动电机可通过差动扭矩产生附加横摆力矩：

$$M_{z,\text{ext}}=\frac{t_w}{2}(F_{x,\text{right}}-F_{x,\text{left}})$$

等效增加系统阻尼。

---

八、相关内容

• 01-阶跃与脉冲转向输入的时域响应解析
• 02-特征车速与临界车速的理论计算与物理意义
• 02-状态空间矩阵构建与系统特征参数分析