Matlab R2014a 一阶系统响应：5种极点位置对阶跃/冲激响应的影响对比

在自动控制系统的分析与设计中，理解系统极点位置对动态响应的影响至关重要。一阶系统作为最基本的动态系统模型，其响应特性直接反映了极点分布与系统性能的内在联系。本文将基于Matlab R2014a环境，通过5种典型极点配置案例，系统分析极点位于s平面不同区域时对阶跃响应和冲激响应的决定性影响。

1. 一阶系统建模基础

一阶系统的传递函数标准形式为：

G(s) = K / (τs + 1)

其中，时间常数τ决定了系统响应速度，极点位置为s = -1/τ。通过调整分母系数，我们可以实现极点在s平面的不同位置分布：

% 创建传递函数的通用命令格式 num = K; % 分子系数 den = [τ 1]; % 分母系数（从高次到低次） sys = tf(num, den);

关键参数说明：

• K：系统增益，影响稳态输出幅值
• τ：时间常数，决定响应速度（τ越小响应越快）
• 极点位置：s = -1/τ，直接影响系统稳定性和动态特性

提示：在Matlab中，tf()函数创建的模型对象可直接用于step()和impulse()函数进行时域响应分析。

2. 左半平面极点分析

2.1 稳定快响应（τ=0.5）

num = 1; den = [0.5 1]; % 极点s=-2 sys1 = tf(num, den); step(sys1); hold on; impulse(sys1);

响应特性：

• 阶跃响应：快速上升（调节时间约3τ=1.5秒），无超调
• 冲激响应：快速衰减（时间常数0.5秒）
• 稳定性：绝对稳定（极点位于负实轴）

2.2 稳定慢响应（τ=2）

den = [2 1]; % 极点s=-0.5 sys2 = tf(num, den); figure; step(sys2);

特性变化：

• 响应速度显著降低（调节时间增至6秒）
• 冲激响应衰减更缓慢
• 系统保持稳定但动态性能下降

3. 右半平面极点分析

3.1 不稳定系统（τ=-1）

den = [-1 1]; % 极点s=1 sys3 = tf(num, den); figure; subplot(1,2,1); step(sys3); title('阶跃响应'); subplot(1,2,2); impulse(sys3); title('冲激响应');

不稳定现象：

• 阶跃响应呈指数发散
• 冲激响应同样发散
• 系统完全不稳定

工程意义：

• 右半平面极点导致系统无法稳定工作
• 实际系统中必须通过控制消除右半平面极点

4. 原点极点分析

4.1 临界稳定系统

den = [1 0]; % 极点s=0 sys4 = tf(num, den); figure; step(sys4);

特殊响应：

• 阶跃响应为斜坡函数（无界输出）
• 冲激响应为阶跃函数
• 系统处于临界稳定状态

应用限制：

• 实际工程中应避免
• 可作为积分器使用（需配合其他稳定环节）

5. 极点位置对比实验

5.1 综合对比脚本

% 定义五种极点配置 systems = { tf(1, [0.2 1]), % 快速稳定 tf(1, [1 1]), % 标准稳定 tf(1, [5 1]), % 慢速稳定 tf(1, [-1 1]), % 不稳定 tf(1, [1 0]) % 临界稳定 }; % 绘制阶跃响应对比 figure; hold on; for i = 1:length(systems) step(systems{i}); end legend('τ=0.2','τ=1','τ=5','τ=-1','积分器'); title('不同极点位置的阶跃响应对比'); % 绘制冲激响应对比 figure; hold on; for i = 1:length(systems) impulse(systems{i}); end legend('τ=0.2','τ=1','τ=5','τ=-1','积分器'); title('不同极点位置的冲激响应对比');

5.2 关键发现

通过对比实验可得出以下规律：

1. 极点实部符号决定稳定性：

   • 实部为负 → 稳定
   • 实部为正 → 不稳定
   • 实部为零 → 临界稳定

2. 极点实部大小影响响应速度：

   • 绝对值越大 → 响应越快
   • 绝对值越小 → 响应越慢

3. 极点虚部影响（二阶系统）：

   • 一阶系统无虚部，响应无振荡

6. 工程应用指导

基于上述分析，在实际工程中处理一阶系统时应：

1. 稳定性设计：

   • 确保极点位于左半平面
   • 避免右半平面极点（会导致系统不稳定）

2. 动态性能调整：

   • 通过调整时间常数τ改变响应速度
   • 快速响应需要更大的|Re(pole)|

3. Matlab实现技巧：

   • 使用pole(sys)直接查看极点位置
   • stepinfo()函数获取详细响应参数

% 获取阶跃响应性能指标 info = stepinfo(sys1); disp(info);

输出结果示例：

RiseTime: 0.2197 SettlingTime: 1.4765 Overshoot: 0 Peak: 1