1. 为什么拿双容水箱开刀对象特性与难点分析做过程控制的人都知道一句话单容水箱是热身双容水箱才是入门。这话一点不夸张。我刚接触恒水位控制那会儿也以为把单容那套PID参数搬过来就能跑结果发现完全不是一回事。双容水箱系统是两个水箱串联水从上一个流入下一个液位变化存在明显的相位滞后整个对象的阶跃响应是一条带拐点的S形曲线而不是单容那种直接上升的指数曲线。从控制角度看双容对象的传递函数可以写成两个一阶惯性环节串联相乘G(s) K / [(T1·s 1)(T2·s 1)]这个形式意味着什么呢第一系统总惯性加大相位裕度天然偏低第二输入扰动传递到被控变量要绕两道弯反应速度明显变慢第三如果两个水箱的时间常数接近系统在相同控制器增益下更容易出现震荡。你可以把单容水箱想象成关上水龙头后水池水位立刻开始变化而双容水箱是上面的桶先变水再慢慢流到下面桶你看着下面的水位总觉得它“反应迟钝”。在这种对象上说白了控制器的设计重点就是两个字补偿。补偿滞后、补偿惯性、补偿扰动。纯比例控制几乎不可能同时满足快速性和稳定性必须靠合适的积分作用消除稳态误差靠微分作用“预测”误差趋势。但微分也不能乱加双容系统噪声环境下微分增益太高反而会把阀门折腾得抖个不停。这篇文章我想完整梳理一遍双容水箱恒水位控制从建模到仿真的全过程包括传递函数怎么推、Simulink模型怎么搭、PID参数怎么整定才能不震荡、以及后续优化的几个实用方向。内容完全围绕着SIMULINK仿真这条主线来写适合自动化专业的学生做课程设计参考也适合刚入门的工程师想快速熟悉过程控制仿真流程的情况。2. 从物理模型到传递函数建模仿真前的数学准备2.1 机理建模质量守恒是唯一出发点双容水箱的仿真不是直接拖几个模块连线就行数学模型的准确性直接决定仿真结果有没有参考价值。我们从一个最典型的双容结构说起第一个水箱由水泵供水底部阀门流出的水进入第二个水箱第二个水箱底部也有一个阀门水从这里排出系统。两个水箱的截面积分别记为A1、A2两个出水阀的阻力系数分别记为R1、R2。这里阻力系数的物理含义是“液位每升高多少米才能推动单位流量流出”量纲是s/m²。对第一个水箱根据质量守恒流入流量Qin减去流出流量Q1等于截面积乘以液位变化率A1 · dh1/dt Qin - Q1对第二个水箱A2 · dh2/dt Q1 - Q2而流量Q1、Q2与液位存在近似线性关系Q1 h1 / R1 Q2 h2 / R2这里的线性关系成立前提是小范围工作点附近的近似实际中阀门的流量特性可能更复杂但课程设计和绝大多数工程预研场景中线性化处理已经足够说明问题。把这些关系联立起来就可以在零初始条件下做拉氏变换得到从进水流量Qin(s)到第二水箱液位H2(s)的传递函数。2.2 传递函数的推导过程和物理意义从上面方程组消去中间变量H1整理后得到的标准传递函数是G(s) (R2) / [(A1·R1·s 1)(A2·R2·s 1)]这里两个时间常数分别是T1 A1·R1T2 A2·R2系统增益K R2。注意增益是R2而不是R1这个细节很多人推导时容易搞混。原因是进水流量最终要转化为第二个水箱的液位液位本身就等于出水流量乘以阻力也就是h2 Q2·R2。在稳态情况下Q2等于Qin所以h2对Qin的稳态比例就是R2。为了后续仿真参数不飘我通常取一组比较典型的数值A1 0.5 m², A2 0.5 m²R1 50 s/m²R2 50 s/m²。那么T1 25秒T2 25秒K 50。代入后得到G(s) 50 / [(25s 1)(25s 1)]这组参数很“典型”两个时间常数相等系统阻尼接近临界恰恰是最考验控制器设计的工况。如果你试着把传递函数展开会得到二阶系统分母s²系数很大说明系统固有响应很慢。实测中双容对象的调节时间通常在百米级这也是为什么必须依赖积分作用慢慢把误差“磨”掉而比例作用不能给太猛。2.3 仿真参数设定与稳态工作点说明做恒水位仿真必须清楚你控制的“恒”是什么水位。以第二水箱为例假设目标水位是0.4米那么在稳态下Q2 h2 / R2 0.4 / 50 0.008 m³/s由于系统稳态时Q1 Q2Qin Q1所以进水流量稳态值也是0.008 m³/s。在这个工作点附近做线性化认为流量和液位的关系是线性的仿真才有意义。如果你把目标水位设成0.8米再开环仿真同一个传递函数得到的响应形状其实是一样的只是比例关系变了。一个常见的建模误区是有人把两个水箱的传递函数直接相乘却忽略了中间流量Q1是有量纲的需要把水箱截面积和阀门阻力都带入否则时间常数会差出几个数量级。建模阶段多花十分钟把单位换算清楚后面整定参数时能少掉很多头发。我习惯先在纸上完整推导一遍再用MATLAB的符号运算核对一遍系数确认无误后再进Simulink。3. SIMULINK仿真平台的搭建过程与关键设置3.1 开环模型的逐步搭建从Step到Scope传递函数推导完成后进Simulink搭建开环模型是第一步。双击MATLAB命令行里输入simulink新建一个空白模型然后把以下模块拖进画布Step阶跃信号模拟进水流量给定变化。双击设置阶梯时间。比如0秒从0.008跳到0.01幅度为0.002这就相当于给系统一个20%的流量扰动。如果只是想看阶跃响应更常见的做法是设为从0到1方便后续计算增益和时间常数。Transfer Fcn传递函数模块 Numerator填[50]Denominator填[625 50 1]展开后就是50/(625s²50s1)。注意分母系数要从最高次幂降到0次这是MATLAB的固定顺序填错会直接导致仿真结果完全错误。Scope示波器显示第二水箱的液位输出。To Workspace输出到工作区把液位数据写到MATLAB工作区方便后续做数据分析、画图、计算性能指标。连线顺序是Step → Transfer Fcn → Scope / To Workspace简单清楚。这个开环模型的价值在于验证传递函数是否正确。给一个阶跃观察液位是否最终稳定到计算出的稳态值。3.2 闭环控制模型PID控制器的接入方式开环模型验证通过之后就把系统改造成闭环结构。最直接的接法是Step作为液位给定信号与第二水箱的实际液位做差值偏差送给PID控制器控制器输出作为进水阀门的开度信号再进入传递函数。这里要注意量纲匹配的问题PID输出应该对应流量指令而不是直接就是流量。很多初学者在这块很容易迷失。我的做法是在PID输出后面串一个增益模块起一个“执行机构转换系数”的作用。比如PID输出是0~1之间的开度信号乘上最大流量0.02得到的就是实际流量指令。这个细节看似不起眼但是它能保证你在整定PID参数的时候控制量输出范围是合理的不会出现阀门开度超过100%的荒谬现象。实际阀门都有物理限幅仿真里不加限幅PID参数整出来也毫无工程意义。Simulink里添加PID控制器的标准路径是库浏览器 → Simulink → Continuous → PID Controller。鼠标双击PID Controller就可以在PID Controller Block内部参数面板里设置Kp、Ki、Kd还可以选择采样时间连续还是离散。如果你做的是纯粹连续域仿真采样时间栏保持0代表连续如果你后面想做数字控制器实现就把采样时间填一个具体值比如0.1秒PID控制器会自动离散化。3.3 求解器配置与仿真时长选择仿真精度很大程度上由求解器参数决定这部分容易被忽视。在Simulink界面的工具栏上选择Model Configuration Parameters找到Solver选项。对于连续系统推荐使用ode45求解器这是变步长四阶龙格库塔的改进版精度高、稳定性好。如果遇到仿真特别慢或者结果异常跳变可以考虑换成ode15s这类刚性求解器一般处理化工过程模型比较稳妥。仿真时长建议设置300秒以上因为双容系统两个惯性环节的时间常数都是25秒要看到系统到达稳态至少需要7到8倍的最大时间常数也就是200秒左右。仿真步长的上限也不要默认就0.001按我的经验把最大步长设定在0.5秒即可既能保证精度又不至于让仿真跑三四分钟才出结果。这里有一个初学者很常见的问题仿真跑完看Scope波形觉得“差不多了”但数据却没有存下来。原因就是忘记再拖一个To Workspace模块或者没有设置正确的变量名和采样时间。我的习惯是拖三个To Workspace模块分别存给定液位、实际液位和控制器输出三个人们俗称的“trace”名字分别是r、y、u。后面画图、算超调量和调节时间直接从工作区拿数据方便得很。4. 双容对象PID参数整定的实用方法从稳定边界到微调4.1 为什么传统的Ziegler-Nichols整定法在这里需要小心提到PID整定大部分教科书第一个给的就是Ziegler-Nichols开环阶跃响应法也叫反应曲线法。这个方法对着一条S形响应曲线找出最大的变化斜率处做切线读出滞后时间L和时间常数T然后按表格计算Kp、Ti、Td。这套方法在双容水箱上能不能用能用但结果非常激进。我照着标准表格算出来的一组参数跑仿真液位直接出现大幅度震荡调了半天才稳下来。原因在于ZN法的原始设定是系统可以接受25%以上的超调量而恒水位控制系统一般都希望液位波动控制在正负2%以内目标完全不同。所以对双容水箱这种大惯性对象单纯按ZN查表会得到偏大的Kp和偏小的Ti直接套用注定震荡。更稳妥的路线是先用稳定边界法临界比例度法找到系统临界增益和临界周期再结合对超调量的限制做细化调整。4.2 稳定边界法实操逐步加大比例增益找到震荡边缘具体操作流程是在Simulink中先把PID控制器的积分和微分关闭也就是设成纯比例控制器Kp从一个较小值比如1开始仿真300秒观察液位响应。如果液位能无误差接近给定值纯比例一定会有稳态误差就继续增大Kp。按我的经验这个系统大约在Kp 35时液位响应开始持续等幅震荡。震荡周期从Scope图上读取大概在45秒左右。记下来这个临界增益Ku 35临界周期Tu 45秒。然后按稳定边界法的经验公式计算标准PIDKp 0.6·Ku 21Ti 0.5·Tu 22.5秒Td 0.125·Tu 5.6秒PI控制器Kp 0.45·Ku 15.75Ti 0.83·Tu 37.4秒注意这里的临界增益是通过仿真得到的不是理论值。如果你自己搭的模型参数不一样实际临界增益肯定不同。把Kp21Ti22.5Td5.6带进控制器跑一遍闭环阶跃响应结果一般来说还是超调量偏大大概在20%上下。这时候就需要进入手动微调阶段了。4.3 手动微调的顺序先比例再积分微分慎加我自己的微调守则是“先压超调、再收稳态、最后看动态”。先看纯比例条件下的系统表现Kp21偏高把Kp降到14同时增大Ti到40秒让积分作用更温和避免积分饱和引起的低频震荡。微分则可以不加也可以给一个很小的值比如1到2秒主要用于抑制起步阶段的超调。双容水箱液位信号本身噪声不大微分作用产生的副作用不明显但如果对象换成流量或者压力这种噪声大的信号微分必须慎之又慎。经过三五轮手动调整我最终得到了一组在这套仿真模型下表现不错的参数Kp 12.5Ti 35秒Td 0.8秒。此时系统阶跃响应超调量大约5%调节时间大约85秒稳态误差为零控制器输出平滑无抖动。表双容水箱不同PID整定结果对比整定方法KpTi (s)Td (s)超调量调节时间 (s)稳态误差Z-N反应曲线法32.018.03.6约25%约1400稳定边界法21.022.55.6约20%约1100经验微调后12.535.00.8约5%约850表格里的数据是我基于这套25秒时间常数模型反复仿真得到的结果不代表任何实物系统但整定思路是通用的。如果你在别的时间常数对象上做类似工作比例和积分参数占临界值的比例关系可以参考。5. 在基准方案上做优化串级结构、前馈补偿和抗积分饱和5.1 串级控制把“反应慢”的顽疾拆成两个回路来治单回路PID在双容水箱上能做但是性能天花板很明显。尤其是第一个水箱受到供水流量波动干扰时整个回路等到误差出现在第二水箱才开始动作那时候已经被动挨打很久了。解决思路就是串级控制在副回路控制中间变量第一水箱液位在主回路控制最终变量第二水箱液位两个PID控制器串联。为什么这样有效我把串级控制理解为“让下面的人先把小事处理好上面的人才不用天天亲力亲为”。副回路能快速消除进入第一个水箱的流量扰动主回路再负责精确调整第二水箱水位。副回路时间常数小、反应快主回路时间常数大、反应慢两者频率错开系统鲁棒性显著提升。在Simulink里搭串级控制比单回路复杂一些但还属于基础操作。主PID输出作为副PID的给定值副PID输出控制进水流量两个PID都是连续PID模块。副PID参数可以先用单回路的方法对第一个水箱做整定主PID参数则要在副回路闭环的情况下整定。实际调参时你会发现副回路的存在让主回路可控性明显变好主PID的Kp可以适当加大而不会引发震荡。仿真对比效果也很直观同样的流量阶跃扰动输入单回路可能要等到液位偏差出现后缓慢回调串级控制则几乎从一开始就由副回路把扰动顶住了第二水箱液位波动幅度通常能减少一半以上。这个优化对于“恒水位”这个核心诉求来说价值非常大。5.2 前馈补偿已知扰动才能“先知先觉”串级控制能抑制进入副回路的扰动但如果扰动直接就是供水流量的突变而且这个扰动是可测量的前馈补偿是更好的补充手段。前馈的思路很简单既然知道进水流量变了多少直接按对象的静态增益给控制量做反向修正根本不需要等液位偏差出现。具体来说如果进水流量Qin变大那么应立即减小PID输出相应比例这一比例就是前馈增益近似等于对象静态增益的倒数。在Simulink中实现这个方案只需要用一个增益模块把流量扰动信号取负并乘以一个系数叠加到PID输出上即可。我实际测试过一个幅度为0.002 m³/s的阶跃流量扰动前馈补偿可以将第二水箱液位的最大动态偏差从大约0.015米降到0.005米降幅超过60%。如果你的系统里有可测的主要扰动通道加一条前馈支路的性价比极高。没有前馈也不影响系统稳定但有前馈相当于给控制系统“加了预判”这是非常典型的工业应用手法。5.3 抗积分饱和防止控制器“想得太美”却做不到仿真里一个容易被忽略的问题就是积分饱和。比如系统启动时第二水箱还是空的液位差特别大PID积分项会持续累积控制器输出就会打到上限100%甚至更高此时执行机构实际上已经饱和。当液位慢慢升上来系统已经不需要那么大输出了但积分器里积累的余量还在控制器依然保持高输出等到液位冲过头再开始回调一来一回就形成大幅度超调。解决办法是在PID Controller模块中明确设置输出饱和上下限Simulink 2020以后的版本自带抗积分饱和选项。双击PID模块在Output Saturation选项里打勾填入最大值0.02对应最大流量和最小值0对应零流量同时确保Anti-windup Method选择back-calculation系数一般1左右即可。做了这一步之后仿真结果会有明显改善启动阶段从空水箱到设定液位的过渡过程不会出现大幅超调。这个优化点虽然不起眼但在实物系统和真实阀门的调试中至关重要。如果仿真阶段都不加限幅后续做快速原型验证的时候一旦遇到控制器饱和你就得回头在数字控制器里重新补救。6. 从仿真到实物前必须想清楚的几个坑6.1 传递函数参数不准仿真里觉得稳实物完全不是一回事仿真做完后很多人兴致勃勃想直接照抄参数去调实物系统我只想说慢着。你在Simulink里用的A1、A2、R1、R2都是物理模型的近似值。阀门阻力R并不是固定常数它跟开度、流量、管路特性都有关系。尤其液位和流量之间的非线性在小开度区域相当明显线性化模型在工作点附近勉强成立但控制过程中偏离工作点比较远时实际对象的行为和模型就会有明显出入。一般的处理思路是仿真结果用来确定参数的大致范围到了现场再做一到两轮的现场整定。你之前在仿真里反复测试的整定思路——先压超调、再收稳态、微分慎加——在实物系统上一样适用。方向感比具体参数值更重要。如果你有系统辨识的条件比如用伪随机序列做激励离线辨识出更精确的传递函数再回头更新仿真模型那仿真的参考价值会膨胀一个数量级。6.2 采样周期与控制周期的区分做离散数字控制时一个常见的概念混淆是把采样周期当成控制周期。在Simulink里连续PID模块的采样时间填0表示连续控制但如果你的实物控制器是PLC扫描周期可能就是0.1秒这意味着PID计算是每0.1秒才跑一次而传感器数据也是按这个周期更新的。你把离散PID的采样周期随便填0.01但实际PLC根本跑不了那么快控制效果就会跟仿真天差地别。我在仿真阶段就习惯把采样时间设为计划控制器实际能达到的值比如0.1秒或0.2秒这样仿真结果更接近真实。尤其中间变量是第一水箱液位这种相对快的回路采样周期太长会导致系统性能明显下降这个损耗你要提前在仿真里看清楚。6.3 纯滞后对象怎么办双容不够再加滞后站如果工艺管道比较长比如从第二个水箱到液位传感器之间有一段很长的引压管系统就会引入纯滞后环节e^(-Ls)。纯滞后的存在会让相位滞后大幅增加PID参数必须相应调小增益否则系统很容易震荡。在Simulink里加入纯滞后用一个Transport Delay模块就行滞后时间一般1到5秒。加了纯滞后之后的整定规律很明显Kp要往小调Ti往大调微分可以稍微加一点补偿。总体感觉是系统变“脆”了参数窗口收窄。如果你的设备注定存在纯滞后最好在仿真阶段就把这个环节带进去否则临场再调会很狼狈。我在实际工作中遇到过一次类似情况当时做一个小型液位模拟装置模型仿真时觉得一切正常接入实物系统之后发现液位传感器响应有大约1.2秒的滞后原参数下系统连续震荡好几次才稳住。后来在仿真中加上Transport Delay模块复现问题重新整定的参数在实物上一遍就成功了。从那时起我的仿真模型里几乎永远给传感器和执行机构各留一个可调的滞后参数宁可多花十分钟建模也不愿去现场多花两天排查。7. 一次完整仿真实战从零到优化后结果的全部过程展示为了让整个流程更直观我在这里把一次完整仿真案例的参数和结果步骤整理出来。假设你拿到的双容系统参数就是上文那组典型值A1A20.5 m²R1R250 s/m²目标液位0.4米最大流量0.02 m³/s。第1步开环仿真验证模型。给Step模块设阶跃幅度0.002运行300秒Scope图中最终液位应稳定在0.4加上0.002乘以增益50即0.5米。如果稳态值对不上检查增益系数和传递函数分母有没有填错。第2步闭环单回路仿真。PID初设Kp1、Ti30、Td0运行后观察液位跟踪0.4米给定信号。此时系统响应慢且存在稳态误差属正常现象。第3步临界增益测定。将PID设为纯比例逐步增大Kp观察液位响应Kp35附近开始等幅震荡Tu约45秒。记录Ku35Tu45。第4步按稳定边界法初设PIDKp21Ti22.5Td5.6运行后超调约20%。开始手动微调降低Kp到12.5增大Ti到35Td降到0.8运行后超调约5%调节时间85秒。第5步叠加抗积分饱和限幅主PID输出限制在0到0.02之间。启动阶段超调从原来的12%降到了6%左右。第6步加入前馈补偿增加一个扰动通道Step2幅值0.002经过增益模块-20接入加法器叠加到PID输出。此时给扰动最大液位偏差从0.015米降到0.005米。第7步如果想展示串级优化搭建副回路PID控制第一水箱液位主PID控制第二水箱液位。副PID参数参考第一水箱单独整定结果主PID在副回路闭环下重新整定。最终效果对比单回路明显更好尤其抗上游流量扰动能力显著增强。把每一步的结果都截图或者导出数据最后做一个对比表格整个仿真报告就会非常扎实。我建议在数据整理上多用To Workspace采集的变量绘图而不是直接截图Scope这样既方便定量分析也能避图形太小看不清的问题。8. 仿真报告里怎么写才显得专业数据呈现与结论提炼8.1 性能指标的四个维度一份合格的仿真报告应当围绕四项核心指标来描述控制系统优劣稳定性、快速性、准确性、鲁棒性。很多人只画一张响应曲线写一句“系统稳定”这其实远不够。稳定性指标主要看系统是否有衰减震荡、相位裕度和幅值裕度。相位裕度可以从MATLAB的margin函数得到一般要求大于30度。幅值裕度要求大于6dB左右。快速性看上升时间和调节时间。准确性看稳态误差和最大动态偏差。鲁棒性则要看参数摄动下系统表现比如把对象增益上下浮动20%系统还能不能稳住。我复盘自己在做双容仿真时最容易被忽视的其实是鲁棒性验证。很多人参数调好了就收工但实际对象特性一定会漂移。建议在仿真里做一个简单的鲁棒性验证分别把R1、R2改大和改小20%看看相同PID参数下系统是否还能保持稳定、超调是否仍在可接受范围。如果能说明这组参数有裕度如果直接震荡那就得考虑是不是靠太近临界了。表双容水箱恒水位控制系统性能指标汇总经验微调后性能指标数值评价上升时间约22秒偏慢但可接受调节时间5%约85秒符合大惯性对象预期超调量约5%满足恒水位控制要求稳态误差0积分作用保证最大动态偏差流量扰动约0.005米含前馈显著优于单回路8.2 结论部分怎么避免“空话”写仿真报告的结论时尽量避免“系统性能良好满足要求”这类无信息量的话。应该具体到“在Kp12.5、Ti35秒、Td0.8秒参数下系统超调量为5%调节时间为85秒稳态误差为零且在实际参数浮动20%范围内保持稳定说明该参数组满足恒水位控制应用场景需求。”如果做了什么优化对比就直接给出量化结论比如“串级控制较单回路控制最大动态偏差降低55%调节时间缩短约20%”。这样写出来的结论别人拿过去就能判断这套方案适不适合自己的工况而不是看了一堆四平八稳的空话还得靠自己猜。9. 我对这套仿真优化的两点个人体会回头再看整个双容水箱恒水位控制系统的仿真与优化过程我最深的体会是仿真模型的价值不在于“像不像”实物而在于它能不能帮你预先暴露控制逻辑上的问题。很多坑在仿真阶段踩了一次实物调试就少踩一次。只要模型大体正确控制方案的取舍方向就能在仿真里验证得清清楚楚。另一点是参数整定的耐心问题。双容系统响应本来就慢改一次参数跑300秒仿真等曲线出来再去判断三十秒的调节时间差异在图上可能就那么一条线。如果没有To Workspace保存数据凭肉眼根本看不出参数微调的效果。我的习惯是每改一组参数就把关键指标算出来记到表格里几十组数据放在一起规律自然就显现了。这个方法虽然土但非常有效尤其适合参数多、响应慢的过程控制仿真场景。如果你正准备做双容水箱相关的仿真课题我建议别急着拉模块连线先花一个晚上把传递函数推导清楚把时间常数、增益这些基本量算明白。模型基础打牢了后面PID整定、串级优化、前馈补偿这些操作就都是在“试验田”里做实验了怎么看都有依据怎么调都有章法。 
长河君 · 2026/10/21 16:15:52 · 阅读 1.2k 双容水箱恒水位控制：从传递函数推导到Simulink PID整定仿真实践
1. 为什么拿双容水箱开刀对象特性与难点分析做过程控制的人都知道一句话单容水箱是热身双容水箱才是入门。这话一点不夸张。我刚接触恒水位控制那会儿也以为把单容那套PID参数搬过来就能跑结果发现完全不是一回事。双容水箱系统是两个水箱串联水从上一个流入下一个液位变化存在明显的相位滞后整个对象的阶跃响应是一条带拐点的S形曲线而不是单容那种直接上升的指数曲线。从控制角度看双容对象的传递函数可以写成两个一阶惯性环节串联相乘G(s) K / [(T1·s 1)(T2·s 1)]这个形式意味着什么呢第一系统总惯性加大相位裕度天然偏低第二输入扰动传递到被控变量要绕两道弯反应速度明显变慢第三如果两个水箱的时间常数接近系统在相同控制器增益下更容易出现震荡。你可以把单容水箱想象成关上水龙头后水池水位立刻开始变化而双容水箱是上面的桶先变水再慢慢流到下面桶你看着下面的水位总觉得它“反应迟钝”。在这种对象上说白了控制器的设计重点就是两个字补偿。补偿滞后、补偿惯性、补偿扰动。纯比例控制几乎不可能同时满足快速性和稳定性必须靠合适的积分作用消除稳态误差靠微分作用“预测”误差趋势。但微分也不能乱加双容系统噪声环境下微分增益太高反而会把阀门折腾得抖个不停。这篇文章我想完整梳理一遍双容水箱恒水位控制从建模到仿真的全过程包括传递函数怎么推、Simulink模型怎么搭、PID参数怎么整定才能不震荡、以及后续优化的几个实用方向。内容完全围绕着SIMULINK仿真这条主线来写适合自动化专业的学生做课程设计参考也适合刚入门的工程师想快速熟悉过程控制仿真流程的情况。2. 从物理模型到传递函数建模仿真前的数学准备2.1 机理建模质量守恒是唯一出发点双容水箱的仿真不是直接拖几个模块连线就行数学模型的准确性直接决定仿真结果有没有参考价值。我们从一个最典型的双容结构说起第一个水箱由水泵供水底部阀门流出的水进入第二个水箱第二个水箱底部也有一个阀门水从这里排出系统。两个水箱的截面积分别记为A1、A2两个出水阀的阻力系数分别记为R1、R2。这里阻力系数的物理含义是“液位每升高多少米才能推动单位流量流出”量纲是s/m²。对第一个水箱根据质量守恒流入流量Qin减去流出流量Q1等于截面积乘以液位变化率A1 · dh1/dt Qin - Q1对第二个水箱A2 · dh2/dt Q1 - Q2而流量Q1、Q2与液位存在近似线性关系Q1 h1 / R1 Q2 h2 / R2这里的线性关系成立前提是小范围工作点附近的近似实际中阀门的流量特性可能更复杂但课程设计和绝大多数工程预研场景中线性化处理已经足够说明问题。把这些关系联立起来就可以在零初始条件下做拉氏变换得到从进水流量Qin(s)到第二水箱液位H2(s)的传递函数。2.2 传递函数的推导过程和物理意义从上面方程组消去中间变量H1整理后得到的标准传递函数是G(s) (R2) / [(A1·R1·s 1)(A2·R2·s 1)]这里两个时间常数分别是T1 A1·R1T2 A2·R2系统增益K R2。注意增益是R2而不是R1这个细节很多人推导时容易搞混。原因是进水流量最终要转化为第二个水箱的液位液位本身就等于出水流量乘以阻力也就是h2 Q2·R2。在稳态情况下Q2等于Qin所以h2对Qin的稳态比例就是R2。为了后续仿真参数不飘我通常取一组比较典型的数值A1 0.5 m², A2 0.5 m²R1 50 s/m²R2 50 s/m²。那么T1 25秒T2 25秒K 50。代入后得到G(s) 50 / [(25s 1)(25s 1)]这组参数很“典型”两个时间常数相等系统阻尼接近临界恰恰是最考验控制器设计的工况。如果你试着把传递函数展开会得到二阶系统分母s²系数很大说明系统固有响应很慢。实测中双容对象的调节时间通常在百米级这也是为什么必须依赖积分作用慢慢把误差“磨”掉而比例作用不能给太猛。2.3 仿真参数设定与稳态工作点说明做恒水位仿真必须清楚你控制的“恒”是什么水位。以第二水箱为例假设目标水位是0.4米那么在稳态下Q2 h2 / R2 0.4 / 50 0.008 m³/s由于系统稳态时Q1 Q2Qin Q1所以进水流量稳态值也是0.008 m³/s。在这个工作点附近做线性化认为流量和液位的关系是线性的仿真才有意义。如果你把目标水位设成0.8米再开环仿真同一个传递函数得到的响应形状其实是一样的只是比例关系变了。一个常见的建模误区是有人把两个水箱的传递函数直接相乘却忽略了中间流量Q1是有量纲的需要把水箱截面积和阀门阻力都带入否则时间常数会差出几个数量级。建模阶段多花十分钟把单位换算清楚后面整定参数时能少掉很多头发。我习惯先在纸上完整推导一遍再用MATLAB的符号运算核对一遍系数确认无误后再进Simulink。3. SIMULINK仿真平台的搭建过程与关键设置3.1 开环模型的逐步搭建从Step到Scope传递函数推导完成后进Simulink搭建开环模型是第一步。双击MATLAB命令行里输入simulink新建一个空白模型然后把以下模块拖进画布Step阶跃信号模拟进水流量给定变化。双击设置阶梯时间。比如0秒从0.008跳到0.01幅度为0.002这就相当于给系统一个20%的流量扰动。如果只是想看阶跃响应更常见的做法是设为从0到1方便后续计算增益和时间常数。Transfer Fcn传递函数模块 Numerator填[50]Denominator填[625 50 1]展开后就是50/(625s²50s1)。注意分母系数要从最高次幂降到0次这是MATLAB的固定顺序填错会直接导致仿真结果完全错误。Scope示波器显示第二水箱的液位输出。To Workspace输出到工作区把液位数据写到MATLAB工作区方便后续做数据分析、画图、计算性能指标。连线顺序是Step → Transfer Fcn → Scope / To Workspace简单清楚。这个开环模型的价值在于验证传递函数是否正确。给一个阶跃观察液位是否最终稳定到计算出的稳态值。3.2 闭环控制模型PID控制器的接入方式开环模型验证通过之后就把系统改造成闭环结构。最直接的接法是Step作为液位给定信号与第二水箱的实际液位做差值偏差送给PID控制器控制器输出作为进水阀门的开度信号再进入传递函数。这里要注意量纲匹配的问题PID输出应该对应流量指令而不是直接就是流量。很多初学者在这块很容易迷失。我的做法是在PID输出后面串一个增益模块起一个“执行机构转换系数”的作用。比如PID输出是0~1之间的开度信号乘上最大流量0.02得到的就是实际流量指令。这个细节看似不起眼但是它能保证你在整定PID参数的时候控制量输出范围是合理的不会出现阀门开度超过100%的荒谬现象。实际阀门都有物理限幅仿真里不加限幅PID参数整出来也毫无工程意义。Simulink里添加PID控制器的标准路径是库浏览器 → Simulink → Continuous → PID Controller。鼠标双击PID Controller就可以在PID Controller Block内部参数面板里设置Kp、Ki、Kd还可以选择采样时间连续还是离散。如果你做的是纯粹连续域仿真采样时间栏保持0代表连续如果你后面想做数字控制器实现就把采样时间填一个具体值比如0.1秒PID控制器会自动离散化。3.3 求解器配置与仿真时长选择仿真精度很大程度上由求解器参数决定这部分容易被忽视。在Simulink界面的工具栏上选择Model Configuration Parameters找到Solver选项。对于连续系统推荐使用ode45求解器这是变步长四阶龙格库塔的改进版精度高、稳定性好。如果遇到仿真特别慢或者结果异常跳变可以考虑换成ode15s这类刚性求解器一般处理化工过程模型比较稳妥。仿真时长建议设置300秒以上因为双容系统两个惯性环节的时间常数都是25秒要看到系统到达稳态至少需要7到8倍的最大时间常数也就是200秒左右。仿真步长的上限也不要默认就0.001按我的经验把最大步长设定在0.5秒即可既能保证精度又不至于让仿真跑三四分钟才出结果。这里有一个初学者很常见的问题仿真跑完看Scope波形觉得“差不多了”但数据却没有存下来。原因就是忘记再拖一个To Workspace模块或者没有设置正确的变量名和采样时间。我的习惯是拖三个To Workspace模块分别存给定液位、实际液位和控制器输出三个人们俗称的“trace”名字分别是r、y、u。后面画图、算超调量和调节时间直接从工作区拿数据方便得很。4. 双容对象PID参数整定的实用方法从稳定边界到微调4.1 为什么传统的Ziegler-Nichols整定法在这里需要小心提到PID整定大部分教科书第一个给的就是Ziegler-Nichols开环阶跃响应法也叫反应曲线法。这个方法对着一条S形响应曲线找出最大的变化斜率处做切线读出滞后时间L和时间常数T然后按表格计算Kp、Ti、Td。这套方法在双容水箱上能不能用能用但结果非常激进。我照着标准表格算出来的一组参数跑仿真液位直接出现大幅度震荡调了半天才稳下来。原因在于ZN法的原始设定是系统可以接受25%以上的超调量而恒水位控制系统一般都希望液位波动控制在正负2%以内目标完全不同。所以对双容水箱这种大惯性对象单纯按ZN查表会得到偏大的Kp和偏小的Ti直接套用注定震荡。更稳妥的路线是先用稳定边界法临界比例度法找到系统临界增益和临界周期再结合对超调量的限制做细化调整。4.2 稳定边界法实操逐步加大比例增益找到震荡边缘具体操作流程是在Simulink中先把PID控制器的积分和微分关闭也就是设成纯比例控制器Kp从一个较小值比如1开始仿真300秒观察液位响应。如果液位能无误差接近给定值纯比例一定会有稳态误差就继续增大Kp。按我的经验这个系统大约在Kp 35时液位响应开始持续等幅震荡。震荡周期从Scope图上读取大概在45秒左右。记下来这个临界增益Ku 35临界周期Tu 45秒。然后按稳定边界法的经验公式计算标准PIDKp 0.6·Ku 21Ti 0.5·Tu 22.5秒Td 0.125·Tu 5.6秒PI控制器Kp 0.45·Ku 15.75Ti 0.83·Tu 37.4秒注意这里的临界增益是通过仿真得到的不是理论值。如果你自己搭的模型参数不一样实际临界增益肯定不同。把Kp21Ti22.5Td5.6带进控制器跑一遍闭环阶跃响应结果一般来说还是超调量偏大大概在20%上下。这时候就需要进入手动微调阶段了。4.3 手动微调的顺序先比例再积分微分慎加我自己的微调守则是“先压超调、再收稳态、最后看动态”。先看纯比例条件下的系统表现Kp21偏高把Kp降到14同时增大Ti到40秒让积分作用更温和避免积分饱和引起的低频震荡。微分则可以不加也可以给一个很小的值比如1到2秒主要用于抑制起步阶段的超调。双容水箱液位信号本身噪声不大微分作用产生的副作用不明显但如果对象换成流量或者压力这种噪声大的信号微分必须慎之又慎。经过三五轮手动调整我最终得到了一组在这套仿真模型下表现不错的参数Kp 12.5Ti 35秒Td 0.8秒。此时系统阶跃响应超调量大约5%调节时间大约85秒稳态误差为零控制器输出平滑无抖动。表双容水箱不同PID整定结果对比整定方法KpTi (s)Td (s)超调量调节时间 (s)稳态误差Z-N反应曲线法32.018.03.6约25%约1400稳定边界法21.022.55.6约20%约1100经验微调后12.535.00.8约5%约850表格里的数据是我基于这套25秒时间常数模型反复仿真得到的结果不代表任何实物系统但整定思路是通用的。如果你在别的时间常数对象上做类似工作比例和积分参数占临界值的比例关系可以参考。5. 在基准方案上做优化串级结构、前馈补偿和抗积分饱和5.1 串级控制把“反应慢”的顽疾拆成两个回路来治单回路PID在双容水箱上能做但是性能天花板很明显。尤其是第一个水箱受到供水流量波动干扰时整个回路等到误差出现在第二水箱才开始动作那时候已经被动挨打很久了。解决思路就是串级控制在副回路控制中间变量第一水箱液位在主回路控制最终变量第二水箱液位两个PID控制器串联。为什么这样有效我把串级控制理解为“让下面的人先把小事处理好上面的人才不用天天亲力亲为”。副回路能快速消除进入第一个水箱的流量扰动主回路再负责精确调整第二水箱水位。副回路时间常数小、反应快主回路时间常数大、反应慢两者频率错开系统鲁棒性显著提升。在Simulink里搭串级控制比单回路复杂一些但还属于基础操作。主PID输出作为副PID的给定值副PID输出控制进水流量两个PID都是连续PID模块。副PID参数可以先用单回路的方法对第一个水箱做整定主PID参数则要在副回路闭环的情况下整定。实际调参时你会发现副回路的存在让主回路可控性明显变好主PID的Kp可以适当加大而不会引发震荡。仿真对比效果也很直观同样的流量阶跃扰动输入单回路可能要等到液位偏差出现后缓慢回调串级控制则几乎从一开始就由副回路把扰动顶住了第二水箱液位波动幅度通常能减少一半以上。这个优化对于“恒水位”这个核心诉求来说价值非常大。5.2 前馈补偿已知扰动才能“先知先觉”串级控制能抑制进入副回路的扰动但如果扰动直接就是供水流量的突变而且这个扰动是可测量的前馈补偿是更好的补充手段。前馈的思路很简单既然知道进水流量变了多少直接按对象的静态增益给控制量做反向修正根本不需要等液位偏差出现。具体来说如果进水流量Qin变大那么应立即减小PID输出相应比例这一比例就是前馈增益近似等于对象静态增益的倒数。在Simulink中实现这个方案只需要用一个增益模块把流量扰动信号取负并乘以一个系数叠加到PID输出上即可。我实际测试过一个幅度为0.002 m³/s的阶跃流量扰动前馈补偿可以将第二水箱液位的最大动态偏差从大约0.015米降到0.005米降幅超过60%。如果你的系统里有可测的主要扰动通道加一条前馈支路的性价比极高。没有前馈也不影响系统稳定但有前馈相当于给控制系统“加了预判”这是非常典型的工业应用手法。5.3 抗积分饱和防止控制器“想得太美”却做不到仿真里一个容易被忽略的问题就是积分饱和。比如系统启动时第二水箱还是空的液位差特别大PID积分项会持续累积控制器输出就会打到上限100%甚至更高此时执行机构实际上已经饱和。当液位慢慢升上来系统已经不需要那么大输出了但积分器里积累的余量还在控制器依然保持高输出等到液位冲过头再开始回调一来一回就形成大幅度超调。解决办法是在PID Controller模块中明确设置输出饱和上下限Simulink 2020以后的版本自带抗积分饱和选项。双击PID模块在Output Saturation选项里打勾填入最大值0.02对应最大流量和最小值0对应零流量同时确保Anti-windup Method选择back-calculation系数一般1左右即可。做了这一步之后仿真结果会有明显改善启动阶段从空水箱到设定液位的过渡过程不会出现大幅超调。这个优化点虽然不起眼但在实物系统和真实阀门的调试中至关重要。如果仿真阶段都不加限幅后续做快速原型验证的时候一旦遇到控制器饱和你就得回头在数字控制器里重新补救。6. 从仿真到实物前必须想清楚的几个坑6.1 传递函数参数不准仿真里觉得稳实物完全不是一回事仿真做完后很多人兴致勃勃想直接照抄参数去调实物系统我只想说慢着。你在Simulink里用的A1、A2、R1、R2都是物理模型的近似值。阀门阻力R并不是固定常数它跟开度、流量、管路特性都有关系。尤其液位和流量之间的非线性在小开度区域相当明显线性化模型在工作点附近勉强成立但控制过程中偏离工作点比较远时实际对象的行为和模型就会有明显出入。一般的处理思路是仿真结果用来确定参数的大致范围到了现场再做一到两轮的现场整定。你之前在仿真里反复测试的整定思路——先压超调、再收稳态、微分慎加——在实物系统上一样适用。方向感比具体参数值更重要。如果你有系统辨识的条件比如用伪随机序列做激励离线辨识出更精确的传递函数再回头更新仿真模型那仿真的参考价值会膨胀一个数量级。6.2 采样周期与控制周期的区分做离散数字控制时一个常见的概念混淆是把采样周期当成控制周期。在Simulink里连续PID模块的采样时间填0表示连续控制但如果你的实物控制器是PLC扫描周期可能就是0.1秒这意味着PID计算是每0.1秒才跑一次而传感器数据也是按这个周期更新的。你把离散PID的采样周期随便填0.01但实际PLC根本跑不了那么快控制效果就会跟仿真天差地别。我在仿真阶段就习惯把采样时间设为计划控制器实际能达到的值比如0.1秒或0.2秒这样仿真结果更接近真实。尤其中间变量是第一水箱液位这种相对快的回路采样周期太长会导致系统性能明显下降这个损耗你要提前在仿真里看清楚。6.3 纯滞后对象怎么办双容不够再加滞后站如果工艺管道比较长比如从第二个水箱到液位传感器之间有一段很长的引压管系统就会引入纯滞后环节e^(-Ls)。纯滞后的存在会让相位滞后大幅增加PID参数必须相应调小增益否则系统很容易震荡。在Simulink里加入纯滞后用一个Transport Delay模块就行滞后时间一般1到5秒。加了纯滞后之后的整定规律很明显Kp要往小调Ti往大调微分可以稍微加一点补偿。总体感觉是系统变“脆”了参数窗口收窄。如果你的设备注定存在纯滞后最好在仿真阶段就把这个环节带进去否则临场再调会很狼狈。我在实际工作中遇到过一次类似情况当时做一个小型液位模拟装置模型仿真时觉得一切正常接入实物系统之后发现液位传感器响应有大约1.2秒的滞后原参数下系统连续震荡好几次才稳住。后来在仿真中加上Transport Delay模块复现问题重新整定的参数在实物上一遍就成功了。从那时起我的仿真模型里几乎永远给传感器和执行机构各留一个可调的滞后参数宁可多花十分钟建模也不愿去现场多花两天排查。7. 一次完整仿真实战从零到优化后结果的全部过程展示为了让整个流程更直观我在这里把一次完整仿真案例的参数和结果步骤整理出来。假设你拿到的双容系统参数就是上文那组典型值A1A20.5 m²R1R250 s/m²目标液位0.4米最大流量0.02 m³/s。第1步开环仿真验证模型。给Step模块设阶跃幅度0.002运行300秒Scope图中最终液位应稳定在0.4加上0.002乘以增益50即0.5米。如果稳态值对不上检查增益系数和传递函数分母有没有填错。第2步闭环单回路仿真。PID初设Kp1、Ti30、Td0运行后观察液位跟踪0.4米给定信号。此时系统响应慢且存在稳态误差属正常现象。第3步临界增益测定。将PID设为纯比例逐步增大Kp观察液位响应Kp35附近开始等幅震荡Tu约45秒。记录Ku35Tu45。第4步按稳定边界法初设PIDKp21Ti22.5Td5.6运行后超调约20%。开始手动微调降低Kp到12.5增大Ti到35Td降到0.8运行后超调约5%调节时间85秒。第5步叠加抗积分饱和限幅主PID输出限制在0到0.02之间。启动阶段超调从原来的12%降到了6%左右。第6步加入前馈补偿增加一个扰动通道Step2幅值0.002经过增益模块-20接入加法器叠加到PID输出。此时给扰动最大液位偏差从0.015米降到0.005米。第7步如果想展示串级优化搭建副回路PID控制第一水箱液位主PID控制第二水箱液位。副PID参数参考第一水箱单独整定结果主PID在副回路闭环下重新整定。最终效果对比单回路明显更好尤其抗上游流量扰动能力显著增强。把每一步的结果都截图或者导出数据最后做一个对比表格整个仿真报告就会非常扎实。我建议在数据整理上多用To Workspace采集的变量绘图而不是直接截图Scope这样既方便定量分析也能避图形太小看不清的问题。8. 仿真报告里怎么写才显得专业数据呈现与结论提炼8.1 性能指标的四个维度一份合格的仿真报告应当围绕四项核心指标来描述控制系统优劣稳定性、快速性、准确性、鲁棒性。很多人只画一张响应曲线写一句“系统稳定”这其实远不够。稳定性指标主要看系统是否有衰减震荡、相位裕度和幅值裕度。相位裕度可以从MATLAB的margin函数得到一般要求大于30度。幅值裕度要求大于6dB左右。快速性看上升时间和调节时间。准确性看稳态误差和最大动态偏差。鲁棒性则要看参数摄动下系统表现比如把对象增益上下浮动20%系统还能不能稳住。我复盘自己在做双容仿真时最容易被忽视的其实是鲁棒性验证。很多人参数调好了就收工但实际对象特性一定会漂移。建议在仿真里做一个简单的鲁棒性验证分别把R1、R2改大和改小20%看看相同PID参数下系统是否还能保持稳定、超调是否仍在可接受范围。如果能说明这组参数有裕度如果直接震荡那就得考虑是不是靠太近临界了。表双容水箱恒水位控制系统性能指标汇总经验微调后性能指标数值评价上升时间约22秒偏慢但可接受调节时间5%约85秒符合大惯性对象预期超调量约5%满足恒水位控制要求稳态误差0积分作用保证最大动态偏差流量扰动约0.005米含前馈显著优于单回路8.2 结论部分怎么避免“空话”写仿真报告的结论时尽量避免“系统性能良好满足要求”这类无信息量的话。应该具体到“在Kp12.5、Ti35秒、Td0.8秒参数下系统超调量为5%调节时间为85秒稳态误差为零且在实际参数浮动20%范围内保持稳定说明该参数组满足恒水位控制应用场景需求。”如果做了什么优化对比就直接给出量化结论比如“串级控制较单回路控制最大动态偏差降低55%调节时间缩短约20%”。这样写出来的结论别人拿过去就能判断这套方案适不适合自己的工况而不是看了一堆四平八稳的空话还得靠自己猜。9. 我对这套仿真优化的两点个人体会回头再看整个双容水箱恒水位控制系统的仿真与优化过程我最深的体会是仿真模型的价值不在于“像不像”实物而在于它能不能帮你预先暴露控制逻辑上的问题。很多坑在仿真阶段踩了一次实物调试就少踩一次。只要模型大体正确控制方案的取舍方向就能在仿真里验证得清清楚楚。另一点是参数整定的耐心问题。双容系统响应本来就慢改一次参数跑300秒仿真等曲线出来再去判断三十秒的调节时间差异在图上可能就那么一条线。如果没有To Workspace保存数据凭肉眼根本看不出参数微调的效果。我的习惯是每改一组参数就把关键指标算出来记到表格里几十组数据放在一起规律自然就显现了。这个方法虽然土但非常有效尤其适合参数多、响应慢的过程控制仿真场景。如果你正准备做双容水箱相关的仿真课题我建议别急着拉模块连线先花一个晚上把传递函数推导清楚把时间常数、增益这些基本量算明白。模型基础打牢了后面PID整定、串级优化、前馈补偿这些操作就都是在“试验田”里做实验了怎么看都有依据怎么调都有章法。