基于Matlab的水箱液位级联控制动态系统模型与整定

水箱液位级联控制的动态系统模型，这是我用Matlab搭过最典型的一套过程控制仿真。项目本身不算复杂：一个水箱、一个调节阀、一台流量计，外环测液位，内环测流量，两个PID串起来，目标就是让液位顶住进水扰动。真正动手写代码之后才发现，单回路改级联不是简单加一个PID，差在对象动态、测量变送器、执行机构以及积分饱和这一堆细节上。这篇文章把模型怎么建、代码怎么写、参数怎么整定完整过一遍，适合正在学过程控制、做Matlab课程设计，或者刚接触串级控制想找个能跑的参考模型的人。

1. 项目思路：为什么液位控制要选级联方案

1.1 单回路液位控制的局限性

只用一个液位控制器直接去控调节阀，也就是最传统的单回路结构，从理论课上讲好像没什么问题：液位高了就关阀，液位低了就开阀。但实际水箱系统不是这么理想。

水箱的物料平衡方程是 A·dh/dt = qin − qout，其中 A 是截面积，h 是液位。这个方程在控制眼里就是一个积分环节，输入侧流量差一点，液位就会一直累积下去，没有自平衡能力。更麻烦的是，调节阀和执行机构本身有惯性，流量计也有测量滞后，进水扰动又直接叠加在对象入口。单回路里液位控制器要同时克服阀的非线性、流量的波动和进水扰动的冲击，而且由于被控对象是积分特性，控制器的增益和积分时间必须放在一个很窄的区间里。增益大了容易等幅振荡，增益小了扰动来了半天压不住。

在实际装置上，进水流量并不稳定，管道压力波动会让同一阀位下的实际流量发生变化。单回路液位控制器看到的是“液位偏了”，然后去调阀，但阀调了之后流量变化还要经过执行机构、管道、流量计等一串滞后，等液位真正有反应已经过了好几秒。这段时间里扰动可能又变了，于是控制器来回追，液位就在设定值附近慢慢晃。

1.2 级联控制的结构与角色分工

级联控制的思路是把“流量”这个中间变量单独拎出来，做成一个快速的内回路。整体结构可以写成下面这条链路：

液位设定值 → 外环液位PI → 流量设定值 → 内环流量PI → 阀位指令 → 执行机构动态 → 实际流量 → 水箱积分 → 液位

流量测量值再反馈回内环。外环只负责“液位该对应多大流量”这件事，内环负责“让实际流量尽快跟上这个设定值”。两层PID分工清晰，外环不需要再关心阀的非线性、管道压差波动这些乱七八糟的问题，因为内环已经把流量控制得比较听话了。

这个结构对水箱特别合适，因为流量是一个变化快、可测量、可快速控制的中间信号。把快回路放里面，慢回路放外面，系统带宽拉开，控制性能比单回路好一个档次。用一句通俗的话讲：外环是管理层，定指标；内环是执行层，抓落实。指标定得慢点没关系，但执行层一定要反应快。

1.3 为什么用Matlab代码完成动态模型

水箱液位级联控制的动态系统模型，用理论公式推算也能分析，但真要研究扰动抑制、参数整定、积分饱和这些工程问题时，手算基本算不动。Matlab的好处在于能把非线性、限幅、测量滞后、执行机构惯性都放进同一个仿真里，而且参数改起来非常直观。

我选择用纯m脚本而不是Simulink来搭这套模型，原因后面会详细说。核心是，只有自己在代码里把每个状态变量写出来，才会真正理解“液位状态、流量状态、测量状态、积分状态”各自是怎么变化的。这也是这篇模型代码最大的价值：它不是黑盒，每一步都能拆开看。

2. 动态系统模型的数学化过程

2.1 从物料平衡方程到液位百分比模型

首先要统一单位，否则写代码的时候会乱成一团。我习惯把液位、流量、阀位都按百分比处理。设最大液位量程为 Hmax，阀门全开时最大体积流量为 Qmax，那么液位百分比 h% = h/Hmax×100，流量百分比 q% = q/Qmax×100。

把前面物料平衡方程两边同时除以 A·Hmax，再换算成百分比：

dh%/dt = Qmax/(A·Hmax) × (qin% − qout%)

这个式子的含义很清楚：进出流量差一个百分比，液位会以固定速率变化。用我代码里的参数 A=0.5 m²、Hmax=1.0 m、Qmax=5e-3 m³/s，系数 Qmax/(A·Hmax)=0.01，也就是说流量差1%，液位每秒变化0.01%。水箱从空到满如果一直开最大流量，大概需要100秒，这个时间尺度决定了外环控制器不能像内环那样激进。

积分对象写成传递函数就是典型的 1/s 形式，这也是为什么液位回路整定要特别谨慎。仿真里如果直接用物理单位，比如 m³/h、Pa、kg/h，很容易出现系数差很多个数量级的问题；归一化成百分比之后，控制器增益直观得多。

2.2 执行机构与流量变送器的动态不可省略

真正的水箱系统里，调节阀不是瞬间开到指定开度的。阀芯、气动执行器或者电动执行器都有惯性，我把它简化成一阶惯性环节：

dqout%/dt = (mv% − qout%)/Tv

mv% 是控制器给出的阀位指令，qout% 是实际出流，Tv 是执行机构时间常数。这个惯性在单回路里会造成相位滞后，在级联系统里则直接影响内环的控制性能。如果把这个惯性设成0，内环会快得离谱，外环参数整定结果完全不可信。

流量计同样不能当理想传感器。电磁流量计、涡街流量计都有响应滞后，尤其在小流量时还可能带噪声。模型里用一阶测量动态：

dqm%/dt = (qout% − qm%)/Tt

Tt 是流量变送器时间常数。这个环节在实际项目中经常被忽略，但忽略之后内环反馈信号等于直接取真值，控制器会“看得太准”，整定出来的增益往往偏大，投到现场就振荡。

2.3 控制器的正反作用与积分饱和

级联系统里有两种作用方向，必须分清楚。

内环是正作用。流量设定值 rq 增大，要求阀位 mv 也要增大，所以误差取 e2 = rq − qm，控制器输出直接加。

外环是反作用。阀开大、流量增大，水箱液位反而是下降的。也就是说，从阀位到液位的通道是负增益。如果外环还用 e = SP − PV，那液位一旦偏高，控制器反而把流量设定值减小，液位会越偏越高，这是典型的正反馈。正确处理是取 e1 = h − SP，液位偏高时误差为正，输出流量设定值增大，出流增大，液位回落。

控制器本身我采用位置式PI，输出限幅：

v = Kp×e + Ki×int，且 0 ≤ v ≤ 100

位置式PI必然要处理积分饱和。我对每个回路都加了条件积分，只有控制器输出没有卡在限幅边界时，或者已经卡住但误差方向有助于退出饱和时，才更新积分项。这段逻辑在代码里看起来不起眼，但省掉它，外环一旦被扰动顶到满量程，积分状态可能冲得很高，扰动消失后液位要过很久才能回来，这就是典型的“恢复太慢”问题。

3. Matlab代码实现与仿真流程

3.1 为什么我选了纯m脚本而不是Simulink

Simulink搭级联控制非常快，双击PID模块就能设参数，模型结构也直观。但我在教学和调试中更推荐纯m脚本，尤其是第一次接触动态系统模型的人。

原因有几个。第一，m脚本里每个状态变量都显式写在积分循环里，不容易出现“模型能跑但不知道里面在算什么”的情况。第二，Simulink里外环PID模块默认是正作用还是反作用容易搞混，经常出现液位设定值越给越高但液位往下掉的奇怪现象，脚本里一行代码就能看到信号方向。第三，采样控制、抗积分饱和等逻辑在Simulink里要么需要额外搭一堆模块，要么容易产生代数环，而脚本里用 if 条件就能实现。

当然Simulink也不是没用，后面想加非线性阀特性、变速泵、多水箱互连，用Simulink改模型结构更省事。我做了一个简单对比：

如果你只是想验证级联控制效果，纯脚本完全够用。下面这段代码就是可以直接复制到Matlab里跑的完整模型。

3.2 水箱液位-流量级联控制完整代码

% ========================================================= % 水箱液位-流量级联控制动态系统仿真 % 被控对象：单水箱液位 % 副变量：出口流量 % 结构：外环液位PI -> 流量设定值 -> 内环流量PI -> 阀位 % ========================================================= clear; clc; close all; %% 1. 物理参数 A = 0.5; % 水箱截面积 m^2 Hmax = 1.0; % 液位量程 m Qmax = 5.0e-3; % 阀门全开时的最大体积流量 m^3/s Tv = 2.0; % 阀门执行机构时间常数 s Tt = 0.5; % 流量变送器时间常数 s %% 2. 控制器参数 % 先整定内环，再整定外环，下面这组参数是我在仿真里反复试出来的 Kp1 = 0.8; Ki1 = 0.02; % 外环液位PI，输出为流量设定值，反作用 Kp2 = 1.5; Ki2 = 0.5; % 内环流量PI，输出为阀位，正作用 %% 3. 仿真设置 Ts = 0.02; % 控制周期，也是固定步长 Tfinal = 400; % 总仿真时间 t = 0:Ts:Tfinal; N = length(t); qin0 = 40; % 初始进水量，百分比 qinStep = 50; % 扰动后的进水量，百分比 tStep = 200; % 加入进水阶跃扰动的时刻 %% 4. 初始化状态，保证初始稳态 % 初始液位50%，进水量等于出水量40%，所有变量都平衡 h = 50; % 液位百分比 qout = 40; % 实际出口流量百分比 qm = 40; % 流量测量值百分比 mv = 40; % 阀门开度百分比 int1 = 40 / Ki1; % 外环积分初值，使PI输出稳态为40% int2 = 40 / Ki2; % 内环积分初值，使PI输出稳态为40% %% 5. 日志数组 hLog = zeros(N,1); qoutLog = zeros(N,1); qmLog = zeros(N,1); mvLog = zeros(N,1); rqLog = zeros(N,1); qinLog = zeros(N,1); hLog(1) = h; qoutLog(1) = qout; qmLog(1) = qm; mvLog(1) = mv; rqLog(1) = 40; qinLog(1) = qin0; %% 6. 主仿真循环 % 固定步长欧拉法 + 采样PI控制器，模拟PLC周期扫描 for k = 1:N-1 % 进水扰动，200秒后进水量从40%阶跃到50% qin = qin0; if t(k) >= tStep qin = qinStep; end % ---------- 外环：液位PI，反作用 ---------- e1 = h - 50; % 液位高则误差为正，需要加大出流 v1 = Kp1 * e1 + Ki1 * int1; v1 = max(0, min(100, v1)); % 流量设定值限幅 % 条件积分：未饱和时积分；饱和后只允许向退出饱和的方向积分 s1 = (v1 > 0 && v1 < 100) || (v1 >= 100 && e1 < 0) || (v1 <= 0 && e1 > 0); if s1 int1 = int1 + e1 * Ts; end rq = v1; % ---------- 内环：流量PI，正作用 ---------- e2 = rq - qm; % 流量设定值与测量值之差 v2 = Kp2 * e2 + Ki2 * int2; v2 = max(0, min(100, v2)); % 阀位限幅 s2 = (v2 > 0 && v2 < 100) || (v2 >= 100 && e2 < 0) || (v2 <= 0 && e2 > 0); if s2 int2 = int2 + e2 * Ts; end mv = v2; % ---------- 被控对象动态 ---------- qout = qout + (mv - qout) / Tv * Ts; % 执行机构惯性 qm = qm + (qout - qm) / Tt * Ts; % 流量测量动态 h = h + Qmax / (A*Hmax) * (qin - qout) * Ts; % 水箱物料平衡 % ---------- 记录 ---------- hLog(k+1) = h; qoutLog(k+1) = qout; qmLog(k+1) = qm; mvLog(k+1) = mv; rqLog(k+1) = rq; qinLog(k+1) = qin; end %% 7. 绘制结果 figure('Name','水箱液位级联控制响应'); subplot(3,1,1); plot(t, hLog, 'LineWidth', 1.5); hold on; plot(t, ones(size(t))*50, 'k--', 'LineWidth', 1); grid on; ylim([45 55]); title('水箱液位响应'); ylabel('液位 %'); xlabel('t/s'); legend('h','设定值','Location','best'); subplot(3,1,2); plot(t, rqLog, 'LineWidth', 1.5); hold on; plot(t, qmLog, 'LineWidth', 1.5); grid on; title('流量内环设定值与测量值'); ylabel('流量 %'); xlabel('t/s'); legend('rq','qm','Location','best'); subplot(3,1,3); plot(t, mvLog, 'LineWidth', 1.5); hold on; plot(t, qoutLog, 'LineWidth', 1.5); plot(t, qinLog, 'k--', 'LineWidth', 1); grid on; title('阀位、实际出流与进水扰动'); ylabel('%'); xlabel('t/s'); legend('mv','qout','qin','Location','best'); %% 8. 性能指标统计 idx = t >= tStep; ise = sum((hLog(idx) - 50).^2) * Ts; peak = max(hLog(idx)); fprintf('扰动后ISE = %.3f，液位峰值 = %.2f%%\n', ise, peak);

3.3 代码里的几个关键细节

这段代码里最容易看漏的地方是初始积分值。我直接让 int1=40/Ki1、int2=40/Ki2，而不是从0开始。原因是，如果积分初值设0，那么初始时刻外环误差为0，输出也为0，阀门会被内环直接关掉，水箱只能靠进水慢慢涨上去，仿真一开始就会有一个很长的过渡过程，看不出级联控制本身的特性。设好积分初值，系统从一开始就处于稳态，200秒以后加的进水阶跃扰动才能被单独观察。

正反作用我也在代码里做了特别处理。外环误差 e1 = h − 50，这个符号和很多教科书上的“设定值减测量值”不一样。因为外环输出是流量设定值，流量越大液位越低，所以必须是反作用。你要是想保持 e1 = 50 − h 的写法，就得把 Kp1 和 Ki1 同时取负号，效果一样，但更容易把自己绕晕。

限幅之后的条件积分逻辑是：

• 输出没有卡在0%或100%时，正常积分；
• 输出已经顶到100%，但误差方向是让控制量减小的，允许积分，因为这是退出饱和；
• 输出已经压到0%，但误差方向是让控制量增加的，允许积分，同样是为退出饱和。

很多初学代码会简单写成“只要输出饱和就把积分清零”，这样看起来合理，其实会让控制器在饱和边界附近非常生硬。上面这种条件积分虽然简单，但用在这套模型里已经足够。

3.4 仿真结果怎么看

运行完代码后，三个子图分别对应液位、流量内环、阀位和扰动。正常情况下，第200秒之前所有曲线都是平的，说明初始稳态设置正确。第200秒进水从40%跳到50%，水面开始缓慢上升，外环液位控制器检测到偏差后，逐渐提高流量设定值，内环流量控制器马上跟踪，阀门开大，实际出流随之增大。

由于外环是PI控制，扰动结束后液位会慢慢回到50%，不会留下余差。内环的 qm 和 rq 两条曲线几乎贴在一起，这说明内环完成得很快。如果内环跟踪不上，液位响应就会出现明显振荡。

性能指标方面，ISE是在扰动后液位偏差平方的累积，峰值是扰动后液位最高点。这两项指标都能用来对比不同参数的好坏。我自己跑这组参数时，液位峰值大约在51%上下，恢复时间在100秒左右，属于比较均衡的水平。如果想要更快恢复，可以适当加大 Ki1，但峰值可能会略微上升。

4. 内环外环整定顺序与实测记录

4.1 为什么必须先从内环开始

级联控制整定的一条铁律是“先内环，后外环”。内环如果把流量控制得又快又稳，外环眼中的对象就近似成一个“给定多少流量就出多少流量”的理想环节，外环只需要按水箱积分特性来整定。反过来，如果内环本身就晃晃悠悠，外环再怎么调都别扭，甚至会出现两个回路互相激励的振荡。

在这套模型里，内环的对象是阀加流量计，时间常数Tv和Tt加起来大约2.5秒，属于快回路。外环的水箱积分时间大约100秒，属于慢回路。两者时间尺度拉开后，级联结构才能发挥优势。如果内环太慢，或者外环太快，两者带宽接近，就会出现频率耦合。

整定步骤我建议这样走：

1. 外环暂时断开，把流量设定值 rq 手动固定在一个工作点，比如40%。
2. 只调内环 Kp2 和 Ki2，给流量设定值加一个5%~10%的阶跃，观察 qm 能否快速无超调地跟上。
3. 内环满意后，恢复外环，把流量设定值交给外环输出。
4. 给进水量加扰动，观察液位恢复情况，再调整外环 Kp1 和 Ki1。

4.2 一组可行的参数整定过程

我在仿真里试过的整定记录大致是这样的。

首先设 Kp2=1.5、Ki2=0.5，内环流量阶跃基本在5秒内稳定，没有明显超调。这个参数相当于流量偏差1%时，阀位变化1.5%，再靠积分把稳态误差消掉。

然后恢复外环，初始选 Kp1=0.3，Ki1=0.01，进水量扰动后液位能回到50%，但很慢，差不多要到350秒才稳定。为了加快恢复，我把 Ki1 加到0.02，Kp1 加到0.8。此时扰动后液位峰值有所上升，但恢复时间缩短到100秒左右。再往上加 Kp1 到1.5，液位开始出现连续小幅振荡，说明外环增益已经太接近稳定边界。

这套试凑过程在m脚本里修改参数非常快，每次改完直接运行，看曲线变化。做课程设计或者工程验证时，先跑出几组数据，再画在一张图上对比，比凭感觉调要有说服力得多。

4.3 级联控制和单回路的直观差异

我特意用同一套水箱参数搭过一个单回路模型，也就是只用液位PI直接控制阀位，外环输出直接当 mv，不经过流量内环。对比下来有两个非常明显的现象。

第一，单回路在同样进水扰动下，液位峰值明显更高。因为阀位变化之后，实际流量还要经过执行机构惯性慢慢变，液位已经被扰动推出去好远了，控制器才开始发挥作用。级联系统里内环可以提前把流量顶上去，液位偏差自然更小。

第二，单回路参数稍微激进一点，液位就振荡。原因就是前面说的，液位是积分对象，加上阀门滞后和测量滞后之后，相位裕度很低。级联系统由于内环把流量整定成接近一阶跟随环节，外环面对的对象特性改善了很多，参数区间更宽。

如果只看最终稳态，两者都能把液位拉回50%，区别全在动态过程。做项目汇报时，把这两条液位曲线放到同一张图上，效果比讲十页原理都直接。

5. 常见问题与排查技巧实录

5.1 仿真发散或者严重振荡时先查什么

这套模型的仿真发散，绝大多数情况不是控制器设计问题，而是数值问题或者符号问题。

最容易被忽略的是仿真步长。固定步长欧拉法对一阶惯性环节的稳定有限制，步长必须远小于时间常数。我代码里 Tv=2、Tt=0.5，Ts=0.02，已经余量很大。如果把 Ts 改成0.5甚至1，内环就会出现数值振荡，看起来像控制器没调好，实际上是积分步长太大。

步长没问题的话，下一步查符号。外环反作用如果搞成 e1=50−h，Kp1又是正数，那液位越高输出越小，水箱会越积越高直到超过量程。这种发散是单调向上，不是振荡。内环符号错误则表现为流量设定值和测量值越差越远。

还有一种常见情况是控制器参数太激进，输出频繁顶在0%或100%。此时虽然条件积分已经做了，但系统仍会出现极限环振荡。判断方法是看阀位曲线，是不是长期贴在边界上反复弹跳。如果是，先把 Kp1 和 Kp2 往下降50%再试。

5.2 积分饱和处理的一点经验

积分饱和在纯脚本仿真里其实已经比现场好处理了，因为你随时能看 int1 和 int2 的变化。我这边踩过的坑是，一开始为了简化，把饱和判断写成“只要输出限幅就让积分停止”，结果扰动结束后液位恢复得非常慢。

举个例子，外环输出顶到100%，积分停止。这时候误差仍然为正，但输出已经最大，所以停止积分不会让误差继续增大。问题是扰动一撤，液位开始回落，误差变成负的，如果还是“饱和就停止积分”，积分状态会一直停在很高的位置，控制器迟迟无法退出来。所以正确的逻辑必须允许“向退出饱和方向积分”的情况。

条件积分虽然简单，但真正常用的抗积分饱和方案还有积分分离和back-calculation。积分分离是偏差大时直接切掉积分，偏差小时恢复。back-calculation是把实际限幅输出和控制器计算输出的差值反馈回去，像额外加了一个小的负反馈。固定步长仿真里条件积分够了，如果要移植到实时系统，建议用back-calculation，行为会更平滑。

5.3 常见问题速查表

这套速查表是我在实际调模型时总结出来的，大部分问题都集中在“单位不统一、符号搞反、初值不匹配”三个根源上。如果代码跑出来的曲线不对，不要急着改 PID 参数，先按这三类排查。

我再多说一句，做水箱液位级联模型，真正难的不是会写PID，而是能不能把物理过程翻译成动态方程，再把动态方程落实到可运行的代码里。代码里的执行机构时间常数、流量计时间常数、水箱截面积这些参数，都不是摆设，它们决定了整个系统的时间尺度。以后接到真实水箱项目，先把这几个时间常数实测出来，跑一遍这套模型，再决定参数，会省非常多的事。