简介：这份文档面向具备基础控制理论与编程能力的高年级本科生、研究生及工程技术人员，围绕自动控制原理实验教学，系统讲解控制系统的建模、分析与校正方法，帮助读者在MATLAB与MWORKS环境下完成从数学模型建立到性能验证的完整实践。内容涵盖六个实验：控制系统数学模型的MWORKS与MATLAB实现、时域与频域分析、典型环节动态性能与稳定性分析、频率特性测量、系统串联校正，以及EF-Servo集成与滤波，每个实验均配有实验目的、内容、步骤与思考题，便于对照操作并加深理解。资源包为1个docx文档，约1.43MB，结构完整、条理清晰，适合作为课程实验指导书或自学练习材料。目前已有677人学习，读者可借此熟悉传递函数、状态空间与Simulink建模，掌握频率法校正与二阶系统时域、频域指标分析，提升动手实践与系统分析能力。
1. 从一张“跑不通”的阶跃响应图说起
做自动控制原理实验，最让人头疼的不是推导传递函数，而是模型建好了、参数填对了，阶跃响应曲线却跟教材对不上。我带过几轮课程设计，几乎每届都有同学卡在同一个地方：MATLAB 里step()出来的曲线和 MWORKS 里搭的模型结果差一截，或者二阶系统阻尼比一改，超调量直接飞到天上。这套实验方案就是围绕六个递进实验展开的——从数学模型建立、时域频域分析、典型环节动态性能、频率特性测量，到串联校正和 EF-Servo 集成滤波，用 MATLAB 与 MWORKS 双工具对照，把控制理论从纸面公式落到可复现的仿真曲线上。适合正在做控制类课程设计、实验报告或毕设前期验证的读者，照着这套结构走，能少踩不少坑。
2. 数学模型建立：TF、ZPK、SS 三种描述与相互转换
控制系统实验的第一步永远是把物理系统翻译成数学对象。MATLAB 和 MWORKS 都支持传递函数（TF）、零极点增益（ZPK）、状态空间（SS）三种 LTI 模型，但两者的函数命名和调用习惯有差异，这是后面所有实验的地基。
2.1 传递函数模型的建立与系数向量约定
传递函数模型的核心是分子分母多项式系数向量。连续系统用降幂排列，缺项必须补零，这是最常见的翻车点。比如 $G(s)=\frac{4}{s(s+1)^2(s+3)}$，分母展开是 $s^4+5s^3+7s^2+3s$，系数向量必须写成[1 5 7 3 0]，末尾那个 0 代表常数项为零，漏掉它系统阶次就错了。
MATLAB 里直接：
% 实验一：建立传递函数模型并做 TF-SS-TF 往返转换 num = 4; den = [1 5 7 3 0]; % s^4+5s^3+7s^2+3s，常数项补0 Gtf = tf(num, den); % 建立TF模型 Gss = ss(Gtf); % 转换为状态空间模型 Gtf_back = tf(Gss); % 再转回TF，验证一致性MWORKS 的 TyControlSystems 包语法接近 Julia 风格，建立方式为G = tf(num, den)，但要注意在脚本开头using TyControlSystems，否则tf函数不识别。我一般会让学生把两种工具各跑一遍，对比Gtf_back和原始Gtf的系数——理论上完全一致，实际因为浮点运算会出现 $10^{-15}$ 量级的尾数误差，这是正常的，不用慌。
2.2 零极点增益模型与状态空间的转换验证
ZPK 模型适合快速判断系统稳定性：极点全在左半平面则稳定，有右半平面极点则发散。实验一要求把 $G(s)=\frac{6}{(s+1)(s+2)(s+3)}$ 建成 ZPK 模型，无零点时零点向量留空z = []，这是 MWORKS 和 MATLAB 都一致的约定。
% ZPK模型建立与转换 z = []; % 无零点 p = [-1 -2 -3]; % 极点 k = 6; % 增益 Gzpk = zpk(z, p, k); Gtf_from_zpk = tf(Gzpk); % 转回TF验证状态空间模型不唯一，ss()转换出来的 A、B、C、D 矩阵取决于算法实现，所以实验里要求“转回 TF 再比较”而不是“比较 SS 矩阵”。这一点很多同学不理解，以为转换出问题，其实是正常现象。判断转换是否正确的标准只有一个：转回传递函数后分子分母系数是否一致。
2.3 子系统串联、并联与反馈连接的实现
实验一第三部分要求对两个子系统 $G_1(s)=\frac{s+1}{s^2+4s+4}$ 和 $G_2(s)=\frac{1}{s+1}$ 做串联、并联、反馈连接。MATLAB 里串联用G = G2 * G1（注意顺序，乘法不交换），并联用G = G1 + G2，反馈用feedback(G1, G2)，默认负反馈。
% 子系统连接 G1 = tf([1 1], [1 4 4]); G2 = tf(1, [1 1]); G_series = G2 * G1; % 串联：G2在前，G1在后 G_parallel = G1 + G2; % 并联 G_feedback = feedback(G1, G2); % 负反馈连接这里有个容易忽略的细节：串联顺序影响结果，G1*G2和G2*G1在 SISO 系统里数学上等价，但状态空间实现不同，如果后续要做状态反馈设计，顺序不能随意调。反馈连接默认负反馈，正反馈要加参数feedback(G1, G2, +1)，实验报告里如果没写清楚正负反馈，答辩时容易被追问。
3. 时域与频域分析：从阶跃响应到奈奎斯特判据
模型建好之后，下一步是看系统“动起来”什么样。时域看阶跃响应、斜坡响应，频域看 Bode 图和 Nyquist 图，两条线互为印证，这是控制实验的核心分析手段。
3.1 典型环节的动态性能对比与参数影响
实验二要求对比例、积分、微分、惯性、延迟五个典型环节分别做阶跃响应，改变参数观察曲线变化。比例环节 $G(s)=K$ 输出直接是输入的 K 倍，无动态过程；积分环节 $G(s)=K/s$ 输出线性上升，K 越大斜率越陡；微分环节 $G(s)=Ts$ 对阶跃输入输出一个冲激，实际仿真里因为步长限制会显示成一个窄脉冲，这是数值仿真的正常现象，不是模型错了。
惯性环节 $G(s)=\frac{K}{Ts+1}$ 的时间常数 T 决定响应快慢，T 越大上升越缓。延迟环节 $G(s)=e^{-\tau s}$ 在 Simulink 里用 Transport Delay 模块实现，MWORKS 里对应延迟模块，参数 $\tau$ 直接填延迟时间。我一般让学生把五个环节的曲线画在同一张图上对比，一眼就能看出谁有超调、谁有稳态误差、谁响应最快。
% 惯性环节不同时间常数的阶跃响应对比 t = 0:0.01:10; for T = [0.1 1 10] G = tf(1, [T 1]); [y, t_out] = step(G, t); plot(t_out, y); hold on; end legend('T=0.1','T=1','T=10');这段代码的关键是step(G, t)指定时间向量，不指定的话 MATLAB 自动选时间范围，不同 T 值下时间轴不一致，曲线没法直接对比。实验报告里要求“记录相应输出曲线并说明参数影响”，曲线时间轴统一是基本要求。
3.2 二阶系统阻尼比与无阻尼振荡频率的独立影响
二阶系统 $G(s)=\frac{\omega_n^2}{s^2+2\zeta\omega_n s+\omega_n^2}$ 是控制实验的绝对主角。实验二要求固定 $\omega_n=5$，分别取 $\zeta=0, 0.5, 0.707, 1, 1.5$，观察超调量、上升时间、调节时间的变化。
$\zeta=0$ 时等幅振荡，超调量 100%；$\zeta=0.707$ 是工程上常用的“最佳阻尼比”，超调约 4.3%；$\zeta \geq 1$ 时无超调但响应变慢。这个规律教材上都有，但亲手跑一遍曲线，印象完全不一样。反过来固定 $\zeta=0.5$，改变 $\omega_n=0.5, 1, 3, 5$，会发现 $\omega_n$ 只影响响应速度，不影响超调量——超调量只由 $\zeta$ 决定，这是二阶系统的核心结论，答辩必问。
% 二阶系统：固定wn，扫描zeta wn = 5; t = 0:0.01:5; for zeta = [0 0.5 0.707 1 1.5] num = wn^2; den = [1 2*zeta*wn wn^2]; G = tf(num, den); step(G, t); hold on; end legend('zeta=0','0.5','0.707','1','1.5');性能指标提取可以用stepinfo()函数，直接返回 RiseTime、Overshoot、SettlingTime，比手动在图上点准确得多。MWORKS 里对应有类似函数，如果版本不支持，就手动按定义算：超调量 = (峰值 - 稳态值)/稳态值 × 100%。
3.3 状态空间模型的响应求解与稳定性判定
实验二还要求对给定状态空间模型分别用step()和lsim()求单位阶跃和单位斜坡响应。step()专用于阶跃输入，lsim()可以喂任意输入信号，包括斜坡、正弦、指数。
% 状态空间模型响应求解 A = [0 1 0; 0 0 1; -6 -10 -4]; B = [0; 0; 1]; C = [6 0 0]; D = 0; G = ss(A, B, C, D); t = 0:0.01:10; u = t; % 单位斜坡输入 [y, t_out, x] = lsim(G, u, t); plot(t_out, y); % 输出响应稳定性判定用eig(A)看特征值实部是否全负，或者pole(G)看极点是否全在左半平面。实验里要求两种方法都做并比较结果，实际上pole(G)内部就是算 A 的特征值，结果一致是必然的，但亲手验证一遍能加深对“极点决定稳定性”的理解。
频域部分，实验二要求对开环 $G(s)=\frac{2}{s^3+9s^2+14s}$ 画 Nyquist 图并用奈奎斯特判据判断闭环稳定性。nyquist(G)直接出图，判据是 Z = P - N，P 是开环右半平面极点数（这里 P=0），N 是 Nyquist 曲线包围 -1 点的圈数，Z=0 则闭环稳定。同时用pole(feedback(G,1))算闭环极点验证，两种方法结论必须一致，不一致说明 Nyquist 曲线包围方向数错了。
4. 典型环节电路模拟与频率特性测量：仿真与实测的对照
实验三和实验四把纯软件仿真往硬件电路上靠，要求设计模拟电路、算传递函数、做软件仿真，再和虚拟实验软件或实测数据对比。这部分是整份方案里最“工程”的环节，也是实验报告最容易写空的地方。
4.1 运放模拟电路的传递函数推导与参数计算
比例、积分、PI、PD、PID、惯性环节都可以用运放加电阻电容搭出来。以积分环节为例，输入电阻 $R_0=100\text{k}\Omega$，反馈电容 $C=1\mu\text{F}$，传递函数 $G(s)=\frac{1}{R_0Cs}=\frac{1}{0.1s}$，时间常数 $T=R_0C=0.1\text{s}$。实验里要求“求取传递函数并在图中标注 $R_1$、$R_0$”，就是让你把电路参数和传递函数系数对应起来。
PID 环节的模拟电路是近似的，理想 PID 的微分项在物理电路里无法实现纯微分（噪声放大无穷大），所以实际电路加了一个小时间常数滤波。实验文档里给出的参数 $R_0=200\text{k}\Omega$、$R_1=100\text{k}\Omega$、$R_2=10\text{k}\Omega$、$R_3=1\text{k}\Omega$、$C_1=1\mu\text{F}$、$C_2=10\mu\text{F}$，对应 $K_P=\frac{R_1}{R_0}=0.5$，$T_i=R_0C_1=0.2\text{s}$，$T_d=\frac{R_1R_2}{R_0}C_2$ 需要按实际电路拓扑算。这部分我一般要求学生先手算再仿真，手算和仿真对不上就回去查电路图，十有八九是电阻编号看串了。
4.2 二阶与三阶系统模拟电路的稳定性边界
实验三第二部分给出典型二阶系统模拟电路，通过电位器 $R_x$ 改变阻尼比，观察过阻尼、临界阻尼、欠阻尼三种响应。三阶系统更狠，开环 $G(s)H(s)=\frac{K}{s(0.1s+1)(0.5s+1)}$，特征方程 $s^3+12s^2+20s+20K=0$，劳斯判据算出稳定范围 $0<K<12$，临界稳定 $K=12$，不稳定 $K>12$。$K=500/R_x$（$R_x$ 单位 kΩ），所以 $R_x$ 的调节范围直接对应稳定性边界。
% 三阶系统稳定性验证：扫描K值看阶跃响应 for K = [5 12 20] num = K; den = [1 12 20 20*K]; % 闭环特征多项式 G = tf(num, den); step(G, 0:0.01:20); hold on; end legend('K=5 稳定','K=12 临界','K=20 不稳定');K=20 时曲线发散，仿真时间一长数值会溢出，这是正常的，实验报告里截取前几秒的发散趋势就够了。这里有个血泪经验：仿真步长不能太大，否则临界稳定那条曲线可能因为数值误差看起来像发散或衰减，步长取 0.01s 比较稳妥。
4.3 频率特性曲线的测量与传递函数反求
实验四要求用实验方法测一阶惯性环节和二阶系统的频率特性，再根据曲线反求传递函数。一阶惯性环节 $G(s)=\frac{K}{Ts+1}$ 的 Nyquist 图是个半圆，半圆直径对应低频增益 K，曲线上某频率点 $\omega_k$ 对应的相角 $\varphi_k$ 满足 $T=\frac{\tan\varphi_k}{\omega_k}$，这是反求时间常数的依据。
软件仿真用bode(G)和nyquist(G)直接出图，和实测数据画在同一张图上对比。实测数据点少、有噪声，曲线不会像仿真那样光滑，实验报告里要说明“实测与仿真的偏差来源”，一般归因于运放非理想特性、电阻电容容差、测量读数误差。我一般要求学生至少测 8 到 10 个频率点，低频段密一点，高频段稀一点，这样 Bode 图的转折频率附近才能看清。
5. 串联校正与 EF-Servo 集成滤波：从分析到设计
前四个实验是“看系统”，实验五和实验六是“改系统”。串联校正用频率法设计超前或滞后网络，EF-Servo 集成滤波则是把控制算法落到具体执行机构上，这部分最接近实际工程。
5.1 频率法串联校正的设计步骤与参数整定
实验五要求对给定开环系统用频率法设计串联校正装置。步骤是：先画未校正系统的 Bode 图，看截止频率和相角裕度；根据指标要求确定需要补偿的相角；选择超前校正网络 $G_c(s)=\frac{\alpha Ts+1}{Ts+1}$，$\alpha>1$；计算 $\alpha$ 和 T；验证校正后系统的相角裕度和截止频率是否达标。
% 串联超前校正设计示例 s = tf('s'); G = 2 / (s*(s+1)*(s+2)); % 未校正系统 margin(G); % 查看原始裕度 % 假设需要提升相角裕度到50度，设计超前网络 alpha = 5; % 超前比 T = 0.1; % 时间常数 Gc = (alpha*T*s + 1) / (T*s + 1); G_open = Gc * G; margin(G_open); % 验证校正后裕度参数整定没有唯一解，$\alpha$ 取大了相角补偿多但高频增益也大，噪声敏感；T 取小了转折频率右移，可能来不及补偿。我一般让学生先按公式算一版，再手动微调，直到margin()显示的相角裕度落在 45° 到 60° 之间，截止频率满足指标。实验报告里要把设计过程和验证结果都写上，只写最终参数不写推导，答辩时会被问住。
5.2 EF-Servo 集成与滤波环节的实现
实验六涉及 EF-Servo 集成，本质是把控制器输出经过滤波后送给伺服执行机构。滤波环节通常是一阶低通 $G_f(s)=\frac{1}{T_f s+1}$，$T_f$ 根据执行机构带宽选，太大响应慢，太小滤不掉噪声。MATLAB 里把滤波器和控制器串联，再和伺服模型闭环连接，仿真看整体响应。
这部分实验文档给的信息相对少，我一般按“控制器 + 滤波器 + 执行机构”三段式搭 Simulink 模型，滤波器时间常数从执行机构手册里查，查不到就按截止频率的 1/5 到 1/10 估。仿真时重点看两点：阶跃响应有没有高频抖振（滤波器没起作用），稳态误差是否达标（控制器增益够不够）。
6. 避坑与排查：实验报告里最容易翻车的五个点
6.1 传递函数系数缺项未补零导致阶次错误
现象：tf(num, den)建出来的模型阶次比预期低，step()曲线形状完全不对。原因：分母多项式中间某次项系数为零时，系数向量里漏写了那个 0。比如 $s^3+2s$ 写成[1 2]就变成了一阶系统。解决：展开多项式后逐项对照，缺项补零，用length(den)-1检查阶次是否等于最高次幂。
6.2 仿真步长过大导致临界稳定系统误判
现象：劳斯判据算出临界稳定，仿真曲线却发散或快速衰减。原因：step()或sim()的步长太大，数值积分误差累积。解决：把仿真步长降到 0.01s 以下，或者用step(G, t)显式指定密集时间向量，临界稳定系统至少跑 20 秒以上看趋势。
6.3 MWORKS 与 MATLAB 函数名混淆
现象：在 MWORKS 脚本里写feedback()报错，或者 MATLAB 里写 Julia 风格语法不识别。原因：两个工具函数库不同，MWORKS 的 TyControlSystems 包函数名和 MATLAB 接近但不完全一样。解决：脚本开头明确using TyControlSystems，不确定的函数先?函数名查帮助，别凭记忆硬写。
6.4 反馈连接正负号搞反导致闭环不稳定
现象：开环稳定系统接成闭环后阶跃响应发散。原因：feedback(G1, G2)默认负反馈，如果实际电路是正反馈接法，仿真和实测对不上。解决：确认电路反馈极性，正反馈加+1参数，实验报告里注明反馈类型。
6.5 频率特性实测点数不足导致反求参数偏差大
现象：根据实测 Bode 图反求的 T 和 K 与电路标称值差很多。原因：测量频率点太少，转折频率附近没有数据点。解决：转折频率 $\omega=1/T$ 附近至少测 3 个点，低频段测增益，高频段测衰减斜率，总共不少于 8 个点。
7. 验收与答辩：怎么证明你的实验真的跑通了
验收这套实验，我一般不看报告写得多厚，就看三样东西：曲线能不能复现、参数改了结果跟不跟着变、异常现象能不能解释清楚。答辩时导师最爱问的是“你这个超调量为什么是 16.3% 而不是教材上的 16.3%”——数字对上了说明你真跑了，对不上但能解释误差来源也算过关，最怕的是支支吾吾说不清曲线怎么来的。
具体验收标准我通常按这张表卡：
| 验收项 | 合格标准 | 检查方法 |
|---|---|---|
| 模型建立 | TF/ZPK/SS 互转后系数一致 | 转回 TF 对比 num/den |
| 时域分析 | 性能指标与理论值偏差 < 5% | stepinfo()读数对比公式 |
| 频域分析 | Nyquist 判据与极点法结论一致 | 两种方法各跑一遍 |
| 电路仿真 | 手算传递函数与仿真曲线吻合 | 改变 R、C 参数看响应变化 |
| 校正设计 | 校正后相角裕度 45°~60° | margin()验证 |
答辩技巧上，我建议每个实验准备一张“对比图”：仿真曲线和实测（或理论）曲线画在一起，偏差用阴影标出来，讲的时候先说结论再说原因。比如二阶系统 $\zeta=0.707$ 时仿真超调 4.3%，实测 5.1%，偏差来自运放带宽限制和电阻容差，这样讲既展示了工作量，又体现了分析能力。
还有一个容易被忽略的点：实验报告里的图必须带坐标轴标签、单位、图例，曲线颜色能区分。我见过太多报告图截得模模糊糊，答辩时导师问“这条线是哪个参数”都答不上来。从那以后我每次让学生交报告前，都强制走一遍“关掉 MATLAB 只看图能不能读懂”的检查——图能自己说话，报告就成功了一半。希望这套结构能帮到你，把六个实验串成一条完整的分析链，而不是六份孤立的作业。
本文还有配套的精品资源，点击获取
