MATLAB控制系统仿真：从建模到分析实战

1. MATLAB控制系统仿真实战指南

作为一名控制工程师，我每天的工作都离不开MATLAB这个强大的工具。它就像我的数字实验室，从传递函数分析到系统响应测试，几乎所有控制系统的仿真需求都能在这里找到解决方案。今天我就带大家深入几个典型场景，分享如何用MATLAB快速实现控制系统的建模与分析。

2. 连续系统建模与分析

2.1 传递函数构建与阶跃响应

让我们从一个简单的二阶系统开始。在MATLAB中构建传递函数就像搭积木一样简单：

matlab复制sys_s = tf([1], [1 2 5]);
step(sys_s);
grid on

这段代码定义了一个传递函数为1/(s²+2s+5)的系统。这里有几个关键细节需要注意：

• tf函数的第一个参数是分子多项式系数，第二个是分母多项式系数
• 系数排列顺序必须是从高阶到低阶（s²在前，s在中间，常数项最后）
• 如果顺序写反了，比如写成[5 2 1]，系统特性将完全改变

提示：新手常犯的错误就是系数顺序弄反。建议在输入多项式时，先在纸上写出标准形式，再按降幂顺序提取系数。

2.2 零极点图与系统稳定性

判断系统稳定性最直观的方法就是观察零极点分布：

matlab复制pzmap(sys_s)
title('连续系统零极点分布')

运行后会显示极点在复平面上的位置。对于这个系统，你会看到两个共轭极点位于左半平面，这说明系统是稳定的。如果极点出现在右半平面，系统将会发散。

3. 微分方程求解实战

3.1 弹簧阻尼系统建模

实际工程中，很多系统都是用微分方程描述的。比如这个经典的弹簧-质量-阻尼系统：

matlab复制m = 2;  % 质量(kg)
c = 1;  % 阻尼系数(N·s/m)
k = 5;  % 弹簧刚度(N/m)

ode_fun = @(t,x) [x(2); (1/m)*(-c*x(2)-k*x(1))];
[t, x] = ode45(ode_fun, [0 10], [0; 1]);  % 初始条件：位置0m，速度1m/s
plot(t, x(:,1));  % 绘制位移响应
xlabel('时间(s)'); ylabel('位移(m)');

这里使用了状态空间法，将二阶微分