MATLAB小白也能懂的LTI系统时域分析：从零输入响应到阶跃响应全攻略

MATLAB零基础玩转LTI系统时域分析：从微分方程到响应曲线实战指南

刚接触信号与系统课程时，看到那些复杂的微分方程和响应曲线总让人望而生畏。但别担心，今天我们就用MATLAB这把"瑞士军刀"，带你轻松拆解LTI（线性时不变）系统的时域分析难题。不同于教科书上的理论推导，我们将通过可复用的代码模板和可视化技巧，让你在30分钟内掌握零输入响应、零状态响应、冲激响应和阶跃响应的完整分析方法。

1. 环境准备与基础概念速成

1.1 快速搭建MATLAB分析环境

工欲善其事，必先利其器。在开始前，请确保你的MATLAB已经安装了Symbolic Math Toolbox（符号数学工具箱），这是进行符号运算的关键。验证方法很简单：

>> ver symbolic

如果看到类似下面的输出，说明工具箱已就绪：

Symbolic Math Toolbox Version 8.1 (R2018b)

常见问题排查：

• 报错"未定义函数'dsolve'"：说明缺少符号工具箱，需要从MATLAB附加功能管理器安装
• heaviside函数报错：较旧版本可能需要改用stepfun函数，或更新到新版MATLAB

1.2 LTI系统时域分析核心概念

用最直白的语言理解几个关键术语：

• 零输入响应：系统初始状态单独作用产生的响应（没有外部输入）
• 零状态响应：仅由外部输入引起的响应（初始状态为零）
• 完全响应= 零输入响应 + 零状态响应
• 冲激响应：系统对单位冲激信号δ(t)的零状态响应
• 阶跃响应：系统对单位阶跃信号u(t)的零状态响应

提示：可以把LTI系统想象成一个弹簧振子，零输入响应就像你拉开弹簧后松手的振动，而零状态响应则是外力持续作用的结果。

2. 微分方程求解实战：从符号计算到数值解

2.1 dsolve函数解决零输入响应

让我们从一个典型二阶系统开始：

D²y + 3Dy + 2y = Dx + 3x

初始条件：y(0)=1, Dy(0)=2

零输入响应求解代码：

eq = 'D2y + 3*Dy + 2*y = 0'; cond = 'y(0)=1, Dy(0)=2'; y_zi = dsolve(eq, cond); pretty(simplify(y_zi))

运行后会得到解析解：

-2 t y_zi = e (4 e - 3)

代码解读：

• D2y表示y的二阶导数，Dy为一阶导数
• simplify()用于简化表达式
• pretty()让输出更易读

2.2 零状态响应求解技巧

假设输入信号x(t)=e^(-3t)u(t)，求零状态响应：

eq_input = 'D2y + 3*Dy + 2*y = Dx + 3*x'; eq_x = 'x = exp(-3*t)*heaviside(t)'; cond_zero = 'y(-0.01)=0, Dy(-0.01)=0'; % 零状态条件 y_zs = dsolve(eq_input, eq_x, cond_zero);

注意：这里使用heaviside(t)表示单位阶跃函数u(t)，在旧版本中可能需要替换为stepfun(t,0)

2.3 全响应可视化完整流程

将零输入和零状态响应叠加，并绘制曲线：

t = 0:0.01:5; y_zi_num = eval(vectorize(y_zi)); % 符号表达式转数值 y_zs_num = eval(vectorize(y_zs.y)); % 提取解的结构体字段 figure; plot(t, y_zi_num, 'r--', 'LineWidth', 2); hold on; plot(t, y_zs_num, 'b-.', 'LineWidth', 2); plot(t, y_zi_num + y_zs_num, 'k-', 'LineWidth', 2); legend('零输入响应', '零状态响应', '完全响应'); xlabel('时间(s)'); ylabel('幅值'); grid on; title('系统响应分量分解');

3. 系统响应分析进阶：冲激与阶跃响应

3.1 使用tf函数创建系统模型

对于数值解法，首先需要建立系统传递函数模型。以上述微分方程为例：

D²y + 3Dy + 2y = Dx + 3x

转换为传递函数形式：

H(s) = (s + 3)/(s² + 3s + 2)

MATLAB实现：

num = [1 3]; % 分子系数 den = [1 3 2]; % 分母系数 sys = tf(num, den);

3.2 一键生成专业响应曲线

冲激响应与阶跃响应对比分析：

t = 0:0.01:4; figure; subplot(2,1,1); impulse(sys, t); grid on; title('冲激响应 h(t)'); subplot(2,1,2); step(sys, t); grid on; title('阶跃响应 g(t)');

参数调整技巧：

• 时间向量t的步长影响曲线平滑度，建议0.01-0.001
• 添加grid on让图形更专业
• 使用subplot创建多图对比

3.3 响应特性参数提取

从阶跃响应曲线获取关键指标：

[g, t] = step(sys); stepinfo = stepinfo(g, t); disp(stepinfo);

输出包含：

• RiseTime（上升时间）
• SettlingTime（稳定时间）
• Overshoot（超调量）
• 等10余个动态特性参数

4. 工程实践中的常见问题解决方案

4.1 高频报错及解决方法

问题1：heaviside函数未定义

% 解决方案A：改用stepfun（需自定义） function y = stepfun(t,t0) y = (t >= t0); end % 解决方案B：更新MATLAB到较新版本

问题2：dsolve返回空解

• 检查微分方程书写格式，特别是导数表示（D2y不是D²y）
• 确认初始条件与方程阶次匹配（n阶方程需要n个初始条件）

问题3：曲线显示异常

% 调整图形显示范围 axis([t_start t_end y_min y_max]);

4.2 高阶系统分析模板

对于三阶系统示例：

D³y + 4D²y + 5Dy + 2y = 2Dx + x

代码模板：

% 符号解法 eq_high = 'D3y + 4*D2y + 5*Dy + 2*y = 2*Dx + x'; cond_high = 'y(0)=0, Dy(0)=0, D2y(0)=0'; y_high = dsolve(eq_high, 'x=heaviside(t)', cond_high); % 数值解法 sys_high = tf([2 1], [1 4 5 2]); step(sys_high);

4.3 自定义输入信号分析

处理任意输入信号f(t)的响应：

t = 0:0.01:10; u = sin(2*pi*0.5*t).*(t>=0); % 0.5Hz正弦信号 lsim(sys, u, t); % 模拟任意输入响应 xlabel('时间(s)'); ylabel('幅值'); title('自定义输入信号响应');

5. 分析结果应用与扩展思路

5.1 系统特性判断技巧

通过响应曲线快速判断系统特性：

5.2 频域与时域分析联动

结合频域分析更全面理解系统：

figure; bode(sys); % 伯德图分析 grid on;

关联技巧：

• 时域响应速度 ↔ 频域带宽
• 时域超调量 ↔ 频域谐振峰值
• 时域稳态误差 ↔ 频域低频增益

5.3 自动化分析脚本开发

创建可复用的分析函数：

function analyze_LTI(num, den, tspan) sys = tf(num, den); figure; subplot(2,2,1); impulse(sys, tspan); subplot(2,2,2); step(sys, tspan); subplot(2,2,3); bode(sys); subplot(2,2,4); pzmap(sys); % 零极点图 end

调用示例：

analyze_LTI([1 3], [1 3 2], 0:0.01:4);