Matlab（R2014a版本）一阶、二阶系统响应：从极点分布到系统稳定性实战解析

1. 一阶系统响应实战解析

在Matlab R2014a中分析一阶系统响应，就像用放大镜观察水滴落入平静水面的过程。我们先从最简单的传递函数开始：

num = 1; den = [1 1]; sys = tf(num, den);

这个传递函数表示系统在s域的数学模型，分母中的"1"对应s的一次项系数，"1"是常数项。我习惯把传递函数想象成弹簧阻尼系统——分母的常数项相当于阻尼系数，一次项系数相当于弹簧刚度。

1.1 阶跃响应特性分析

执行step(sys)命令后，你会看到一条平滑上升的曲线。这个响应有三个关键特征：

• 无振荡：曲线像爬坡一样单调上升
• 快速稳定：约3秒后基本达到稳态值
• 时间常数：曲线达到63.2%稳态值的时间正好是1秒（与分母常数项对应）

改变分母系数时，响应速度会明显变化。试试这个对比实验：

den_list = {[1 1], [1 2], [1 4], [1 6], [1 8]}; hold on for i = 1:length(den_list) sys = tf(1, den_list{i}); step(sys) end legend('1','2','4','6','8')

你会发现随着分母一次项系数增大，系统响应变慢，就像在糖浆中移动的物体阻力越来越大。

1.2 极点位置对稳定性的影响

极点位置是理解系统稳定性的钥匙。在Matlab中可以通过pole(sys)查看极点位置：

1. 左半平面极点（如den=[1 1]）：

   • 极点：s=-1
   • 响应稳定衰减
   • 实际物理系统大多属于这种情况

2. 右半平面极点（如den=[1 -1]）：

   • 极点：s=1
   • 响应呈指数发散
   • 系统不稳定，像倒立的摆

3. 原点极点（如den=[1 0]）：

   • 极点：s=0
   • 响应无限增长
   • 典型例子是积分器

我曾在项目中遇到过右半平面极点的情况，当时系统输出不断增大直到饱和，后来通过调整控制器参数才解决。

1.3 脉冲响应对比分析

脉冲响应能揭示系统的本质特性。Matlab提供两种方法：

方法一：人工构造脉冲信号

imp = [1; zeros(99,1)]; % 生成脉冲信号 t = 0:0.1:9.9; y = filter(num, den, imp); plot(t, y)

方法二：使用impulse函数

impulse(sys)

对于稳定系统（极点左半平面），脉冲响应会衰减到零；不稳定系统则呈指数增长。建议同时观察两种响应，能更全面理解系统特性。

2. 二阶系统深度解析

二阶系统就像钟摆运动，比一阶系统复杂得多。典型二阶传递函数形式为：

wn = 2; % 自然频率 zeta = 0.5; % 阻尼比 sys = tf(wn^2, [1 2*zeta*wn wn^2]);

2.1 极点分布与响应形态

二阶系统的极点位置决定了四种典型响应：

1. 过阻尼（ζ>1）：

   • 两个实极点
   • 响应缓慢无超调
   • 示例：den=[1 5 4]

2. 临界阻尼（ζ=1）：

   • 重实极点
   • 最快无超调响应
   • 示例：den=[1 4 4]

3. 欠阻尼（0<ζ<1）：

   • 共轭复极点
   • 振荡衰减
   • 示例：den=[1 2 4]

4. 无阻尼（ζ=0）：

   • 纯虚极点
   • 持续振荡
   • 示例：den=[1 0 4]

通过下面代码可以直观比较不同阻尼比的影响：

zeta_list = [0.2, 0.5, 0.7, 1, 1.5]; hold on for zeta = zeta_list sys = tf(wn^2, [1 2*zeta*wn wn^2]); step(sys) end

2.2 性能指标量化分析

二阶系统的动态性能可以用几个关键指标衡量：

1. 上升时间：响应从10%到90%稳态值的时间

   • 与wn和ζ都相关

2. 峰值时间：达到第一个峰值的时间

   • 计算公式：π/(wn√(1-ζ²))

3. 超调量：最大超出稳态值的百分比

   • 计算公式：100*e^(-ζπ/√(1-ζ²))

4. 调节时间：进入±5%误差带的时间

   • 近似公式：3/(ζwn)

在Matlab中可以通过右键点击响应曲线，选择"Characteristics"查看这些指标。我曾用这些指标优化过伺服系统的响应速度。

2.3 复极点特殊情况

当极点同时具有实部和虚部时，响应呈现振荡衰减特性：

sys = tf([3], [1 2 10]); % 极点：-1±3i step(sys)

这种响应在电机控制中很常见。实部决定衰减速度，虚部决定振荡频率。通过调整这两个参数，可以平衡系统的快速性和稳定性。

3. 稳定性判定方法

3.1 极点分布判据

最直观的稳定性判据就是观察极点位置：

• 所有极点左半平面：系统稳定
• 任一极点右半平面：系统不稳定
• 虚轴极点：临界稳定（持续振荡）

在Matlab中快速检查极点：

p = pole(sys); if all(real(p) < 0) disp('系统稳定'); else disp('系统不稳定'); end

3.2 时域响应观察法

即使不计算极点，通过观察响应曲线也能判断稳定性：

• 收敛到稳态值：稳定
• 发散或持续振荡：不稳定
• 等幅振荡：临界稳定

这种方法在调试实际系统时特别实用，我经常先用step()函数快速评估系统大体性能。

4. 综合应用案例

4.1 温度控制系统分析

假设有一个温度控制系统，其传递函数为：

sys = tf([5], [1 3 6]); step(sys);

1. 计算极点：

   p = pole(sys) % 得到-1.5±1.9365i

2. 判断稳定性：极点实部为负，系统稳定
3. 性能评估：欠阻尼响应，超调量约16%

4.2 机械振动系统调试

对于机械系统：

m = 1; % 质量 b = 4; % 阻尼 k = 5; % 刚度 sys = tf([1], [m b k]);

通过调整b值可以改变阻尼特性。当b=2√(mk)≈4.472时，系统达到临界阻尼，响应最快且无超调。

5. 常见问题排查

在实际使用中容易遇到几个典型问题：

1. 响应曲线异常：

   • 检查传递函数系数是否正确
   • 确认采样时间设置合理

2. 数值不稳定：

   • 避免极高阶系统
   • 尝试不同的求解器

3. 结果与理论不符：

   • 检查单位是否一致
   • 确认模型简化合理

记得有次我的仿真结果总是发散，后来发现是传递函数分子分母系数顺序写反了。这种低级错误往往最难发现，建议把常用代码封装成函数。