系统阻尼比的定义与计算

阻尼比的定义与计算

阻尼比 $\zeta$ 是描述系统振荡衰减程度的重要参数。在控制系统、振动分析等领域，阻尼比通常通过以下公式计算：

$$\zeta =-\frac{\text{实部}(\lambda )}{\sqrt{\text{实部}(\lambda {)}^2+\text{虚部}(\lambda {)}^2}}$$

其中：

• $\lambda =\sigma +j\omega$ 是系统的特征值，
  • $\sigma =\text{实部}(\lambda )$：特征值的实部，对应系统的衰减率；
  • $\omega =\text{虚部}(\lambda )$：特征值的虚部，对应系统的振荡角频率。

---

推导背景

对于二阶系统（例如机械振动、RLC电路），标准状态方程为：
$$\ddot{x}+2\zeta {\omega }_n\dot{x}+{\omega }_n^{2}x=0$$
其特征方程为：
$${\lambda }^2+2\zeta {\omega }_n\lambda +{\omega }_n^2=0$$
特征值解为：
$$\lambda =-\zeta {\omega }_n\pm j{\omega }_n\sqrt{1-{\zeta }^2}$$

由此可得：

• 实部 $\sigma =-\zeta {\omega }_n$
• 虚部 $\omega ={\omega }_n\sqrt{1-{\zeta }^2}$

通过上述关系，可以计算阻尼比 $\zeta$。

---

具体步骤

给定特征值计算阻尼比

1. 确定特征值：
   • 从系统矩阵（如状态矩阵 $A$）中求解特征值 $\lambda$。
2. 提取实部和虚部：
   • $\sigma =\text{实部}(\lambda )$
   • $\omega =\text{虚部}(\lambda )$
3. 代入公式：
   $$\zeta =-\frac{\sigma }{\sqrt{{\sigma }^2+{\omega }^2}}$$

---

阻尼比的物理意义

1. $\zeta =0$：无阻尼，系统振荡永不衰减（理想条件）。
2. $0<\zeta <1$：欠阻尼，系统存在振荡，但振幅逐渐衰减。
3. $\zeta =1$：临界阻尼，系统快速返回平衡，无振荡。
4. $\zeta >1$：过阻尼，系统缓慢返回平衡，无振荡。

---

MATLAB 计算示例

以下是用 MATLAB 计算特征值和阻尼比的代码示例：

% 定义状态矩阵 A
A = [...]; % 替换为实际矩阵

% 计算特征值
eig_values = eig(A);

% 初始化阻尼比数组
zeta = zeros(size(eig_values));

% 计算每个特征值对应的阻尼比
for i = 1:length(eig_values)
    sigma = real(eig_values(i)); % 实部
    omega = imag(eig_values(i)); % 虚部
    zeta(i) = -sigma / sqrt(sigma^2 + omega^2); % 阻尼比公式
end

% 输出结果
disp(zeta);

---

总结

阻尼比 $\zeta$ 是通过特征值的实部和虚部计算得出的，公式为
$$\zeta =-\frac{\sigma }{\sqrt{{\sigma }^2+{\omega }^2}}$$
它是系统性能的重要指标，用于判断系统的振荡和稳定特性。