基于改进WOA预测PID的电缆线径控制.pdf

Journal of Sensor Technology and Application 传感器技术与应用传感器技术与应用,2023,11(3),247-256 Published Online May 2023 in Hans.https:/www.hanspub.org/journal/jsta https:/doi.org/10.12677/jsta.2023.113027 文章引用文章引用:程元清,郭金妹.基于改进 WOA 预测 PID 的电缆线径控制J.传感器技术与应用,2023,11(3):247-256.DOI:10.12677/jsta.2023.113027 基于改进基于改进WOA预测预测PID的电缆线径控制的电缆线径控制 程元清程元清1，郭金妹，郭金妹2 1江西理工大学电气工程与自动化学院，江西 赣州 2江西应用技术职业学院机械与电子工程学院，江西 赣州 收稿日期：2023年4月13日；录用日期：2023年5月17日；发布日期：2023年5月25日 摘摘 要要 针对电线电缆生产推挤系统自身具有的时滞性针对电线电缆生产推挤系统自身具有的时滞性、非线性和大扰动等特点导致生产推挤过程线径厚度控制非线性和大扰动等特点导致生产推挤过程线径厚度控制较难和建模难度较高的问题，提出一种改进鲸鱼算法较难和建模难度较高的问题，提出一种改进鲸鱼算法(IWOA)、动态控制矩阵、动态控制矩阵(DMC)与传统与传统PID相结合的相结合的算法对推挤生产系统进行优化的策略。这种策略集合了算法对推挤生产系统进行优化的策略。这种策略集合了PID算法的较少的稳态误差优点、算法的较少的稳态误差优点、DMC预测克服预测克服时滞问题的优点和时滞问题的优点和IWOA自动寻优整定参数的优点于一身。通过在自动寻优整定参数的优点于一身。通过在MATLAB平台上对推挤系统进行建模平台上对推挤系统进行建模和仿真，对比传统和仿真，对比传统PID算法，未改进的算法，未改进的WOA算法的效果。实验仿真表明，具备预测能力的算法的效果。实验仿真表明，具备预测能力的IWOA控制策控制策略可以使系统达到较小的超调量、较快的响应速度和较强抗扰性等优点，具有一定的实略可以使系统达到较小的超调量、较快的响应速度和较强抗扰性等优点，具有一定的实际生产应用价值。际生产应用价值。关键词关键词 线径系统线径系统，动态矩阵预测动态矩阵预测，改进鲸鱼算法改进鲸鱼算法，PID Cable Diameter Control Based on Improved WOA Predictive PID Yuanqing Cheng1,Jinmei Guo2 1School of Electrical Engineering and Automation,Jiangxi University of Science and Technology,Ganzhou Jiangxi 2College of Mechanical and Electronic Engineering,Jiangxi College of Applied Technology,Ganzhou Jiangxi Received:Apr.13th,2023;accepted:May 17th,2023;published:May 25th,2023 Abstract For the current stagnation,non-linearity and large disturbances of the production and push sys-tem of electric wire and cable production lead to the problem of difficult control and high model-ing of the production push process.IWOA,dynamic control matrix(DMC)and traditional PID algo-rithms optimize the push production system.This strategy integrates the advantages of less steady-state errors in the PID algorithm,the advantages of DMC prediction to overcome the problem of 程元清，郭金妹 DOI:10.12677/jsta.2023.113027 248 传感器技术与应用 time stagnation,and the advantages of IWOA to automatically find out the parameters.By model-ing and simulation on the push system on the MATLAB platform,comparing the traditional PID algorithm and the effect of the unrestrained WOA algorithm,experimental simulation shows that IWOA control strategies with predictive ability can enable the system to achieve a small over-adjust-ment,a fast response speed,and strong antipity,and have certain actual production and applica-tion value.Keywords Wire Diameter System,Dynamic Matrix Prediction,Improve Whale Algorithm,PID Copyright 2023 by author(s)and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License(CC BY 4.0).http:/creativecommons.org/licenses/by/4.0/1.引言引言 推挤绝缘是 PTFE 电线电缆生产制造之中非常重要的一项工序。与之对应的绝缘推挤机是完成推挤动作的最主要的机器设备。推挤动作要求达到的工艺效果是在线芯圆周紧密的挤出覆盖一层熔融的 PTFE材质外壳，要求其厚度要均匀，不得超出规定的误差，对于厚度的精度有较高的要求，同时不得有外壳断裂导致内芯线裸露。在过去线径生产工艺常常用到的控制方法是 PID 控制，因为其具有效果较为稳定，方法简单等优点而被大量应用到该工艺之中，但不足之处是其参数的选取比较困难，需要人有十分丰富的经验才能快速得到较为满意的三个参数。在参数整定结束后，参数的值保持恒定不变，比较难达到最好的控制效果。在真实的生产条件下，由于 PTFE 绝缘电缆生产系统具有时变性、非线性、大滞后性的特点 1。所以其控制效果未达到电缆生产的高精度、抗干扰性能强等要求。故需要寻找新的控制方法来优化上述问题，以满足实际电缆生产的需求。21 世纪以来，智能控制科学发展突飞猛进，涌现出非常多智能控制算法。故许多专家学者将智能控制算法引入线径控制之中，以此来帮助解决线径生产中的问题。针对线径生产系统中的大滞后问题，Hussein H.2将动态矩阵预测的方法引入系统中。针对控制系统 PID 参数整定问题。Li H.3等人将蚁群优化算法引入系统，并结合卡尔曼滤波算法来增强系统整体的抗干扰能力。但会导致整体的运算量剧增，对传递函数也有影响，同时对设备的要求较高，会使控制结果出现一定的偏差。上述对于线径推挤系统的改进方法都在一定程度上增强了的调节能力，但仍然存在一些不足，例如控制精度有待提升，调节时间稍长，易陷入局部最优。本文引入改进的鲸鱼算法(Improved Whale Opti-mization Algorithm,IWOA)和动态矩阵控制(Dynamic Matrix Control,DMC)结合的方法。这种方法的优势在于利用 IWOA 对 PID 的三个参数进行自动整定，解决 PID 调参问题。同时用 DMC 4建立系统模型的预测控制算法，解决系统本身的大滞后问题，具有一定的预判性，使系统可以对干扰做出及时准确的调整。较好的解决了线径生产过程中的控制精度差、非线性和大滞后等问题，具有一定的实际生产价值。1.1.线径生产推挤系统线径生产推挤系统 1.1.1.线径控制系统结构线径控制系统结构 在理想条件下，电缆线径的生产制作可理解为将 PTFT 材质熔融之后包裹在一根铜芯线上。完成这一动作的系统按照从头到尾排序主要由放线设备、牵引设备、加热设备，挤塑设备、冷却设备、测量反Open AccessOpen Access程元清，郭金妹 DOI:10.12677/jsta.2023.113027 249 传感器技术与应用 馈设备、牵引设备和质量反馈设备组成。其系统的整体结构示意图如图 1 所示。Figure 1.Structure of wire diameter extrusion system 图图 1.线径推挤系统结构 金属芯线先经由放线装置和矫正装置，推挤电机作用于螺杆，芯线被螺杆以一定的推挤速度进入推挤机。同时 PTFE 原料经由供料管道被吸入供料部分，并熔融至可以被推挤的状态。芯线和熔融态的 PTFE相遇后经推挤处理后一起被拉出挤塑机头部。然后经过一小段间隔，之后被测径仪检测半径是否合格。之后电缆进过烧结炉和电火花机将质量情况的数据上传，最后电缆由收线装置一起整装待出厂。1.1.2.对象模型建立对象模型建立 出于对测径仪的保护要求，机尾和测径仪之间要间隔一定的距离，防止烧坏测径仪。这导致反馈给系统的数据有一定的延迟，带来了一定的滞后性。为方便生产，通常将推挤速度设置为恒定。于是系统的滞后时间 可用公式求得为 tractionLV=(1)其中 L 为机尾到测径仪的长度，Vtraction为系统的推挤速度，即生产速度。其中滞后的环节可以看做是一阶惯性加纯滞后环节，螺杆调速系统为直流闭环调速，模型可等效为二阶惯性环节 5。故其数学模型可以用如下近似表达：()()()12e11sKG sT sT s=+(2)其中，K 为比例系数，T1、T2为惯性时间常数，为时滞时间常数 6。根据相关研究成果 1 2 3，并且对现场数据收集分析拟合得出较为符合本系统的传递函数为：()()()0.880.61e0.841 1.461sG sss=+(3)2.基于基于 IWOA 预测预测 PID 控制系统设计控制系统设计 2.1.系统运作原理系统运作原理 在实际线径制造过程中，如果两节胚料的链接处连接不牢固或者两节之间存在空气等干扰都会影响程元清，郭金妹 DOI:10.12677/jsta.2023.113027 250 传感器技术与应用 线径的推挤效果。从现场的线径波动数据分析，其差距缓缓减小，然后剧增，振幅较大。此类扰动可等效为一干扰信号。根据分析建立 IWOA 预测 PID 线径控制系统，如图 2 所示。Figure 2.IWOA predictive PID control principle 图图 2.IWOA 预测 PID 控制原理 系统的控制部分由控制和预测两部分组合共同作用。控制部分的算法是 IWOA 与 PID 结合，其原理是将系统控制的目标线径 yr与实际线径 y 的差值 er和差值变化率 ec作为整个系统的输入量，之后按照适应度作用最后整定出最佳的 Kp、Ki和 Kd三个参数。预测部分的算法是将系统各参数拟合到一个预测模型之中，在矫正环节之后输出下一时刻的预测线径。然后 u 是推挤机的控制量(电压)，作用于控制对象推挤机。最终使得推挤机的控制效果达到最优，来达到最佳的线径控制的效果。2.2.鲸鱼算法及其改进鲸鱼算法及其改进 2.2.1.鲸鱼算法鲸鱼算法 鲸鱼优化算法 7的灵感来源与白鲸在大海中的群体猎食行为，目标鱼群代表的是一个解，每个鲸鱼的位置可看做是一个可行解。算法的主体精髓包含三种策略，一种是围猎，群体向目标游动。另一种是围绕猎物呈圆形游动并吐出气泡驱逐猎物。最后一种是基于本身的位置向周围寻找猎物。其动作原理如下：1)围猎。在目标鱼群位置为止的情况下，选取鲸群中当前最优的位置为目标位置，剩下的鲸鱼向当前目标位置聚集，更新迭代，直至找出最优解。此一系列动作可用以下数学模型解释：()()()()122X tXtA DDC XtX tAa raCr+=(4)模型中的 t 表示为目前的迭代次数，X(t)和 X*(t)分别为第 t 代和截止 t 为止的最优解。D 为 X(t)和 X*(t)的空间距离的绝对值。A、C 为矩阵系数，两者由随机向量0,1r和线性衰减系数 a 决定。a 的表达式为：Max22tat=(5)其中 tMax是算法的最大迭代次数。2)气泡攻击。鲸鱼呈螺旋状游动并吐出气泡，气泡因为浮力上升并组成空气墙，攻击鱼群。此阶段的动作可用如下数学模型表示：()()()1ecos 2blX tDlXt+=+(6)程元清，郭金妹 DOI:10.12677/jsta.2023.113027 251 传感器技术与应用 其中 b 为螺旋系数，0,1l为向量。3)搜寻目标。鲸鱼在全局范围内寻找目标鱼群，此阶段动作可用如下模型表示：()()()()randrand1DC XtX tX tXtA D=+=(7)模型中 Xrand为鲸群中随机一个体。在算法的初期，1A 时算法执行全局搜索，在1A 时，算法执行局部搜索，便于快速找到最优解。局部搜索时为实现包围和缩圈的同时运行，采取用相同的 p 执行包围，用不同的概率 1 p 执行螺旋缩圈，可用如下数学模型表示：()()()()0.51ecos 20.5blXtA DpX tDlXtp+=+(8)2.2.2.鲸鱼算法的改进鲸鱼算法的改进 由于初代鲸鱼算法的一些特性，比如其在前期执行的是全局搜索，且收敛速度较为迅速，后期执行的是局部搜索，这会导致其陷入局部最优解的困境。这也是诸多生物智能算法的通病。另外鲸鱼在搜寻目标阶段的策略较为单一古板，搜寻效果偏向于“随机”，相应的运算量比较繁琐且效率有待提升。1)针对鲸鱼算法之所以容易陷入局部最优解的问题，是因为其 a 的分布呈线性，而原始算法的 a 是在初始时就呈线性下降，不利于全局搜索，而在后期变化较慢，不利于全局搜索。应尽量多的使 a 停留在“全局搜索”的区域内，故本文引入新的函数(9)代替原始鲸鱼算法的函数。()2Max2 1ta tt=(9)函数(9)在算法执行的前期变化较为缓慢，有利于增强前期俯瞰全局的能力，大大增强其全局搜索性。而在算法执行后期其迅速收敛直至变为 0，使其局部搜索能力大幅提升。2)针对搜寻效率不高且搜索不均衡的问题，现引入 Levy 飞行策略 8。Levy 被论证为是一种较为高效的仿生随机游走搜索的策略。其执行动作大部分为短步搜索，一小部分为长步跨越穿插在其中。这种长步跨越使得 Levy 自身具有鲸鱼算法那样精确搜索特性的同时，也具有一定的探索性。有利于算法在迭代时在区域性最优解附近长步探索更优解，避免其陷入局部最优解的误区。为方便在数学模型上复现此策略，本文所使用的 Levy 飞行数学表示如下 9：()()1sin12Levy ss+(10)其中 s 为步长，s 的表达式如下：1us=(11)其中()20,uuN，()20,uvN，的取值为 1.5，u和 v的取值为如下：()112 1sin2,1122uv+=+(12)程元清，郭金妹 DOI:10.12677/jsta.2023.113027 252 传感器技术与应用 将新的 Levy 飞行策略优化至原始鲸鱼算法中，新的搜寻策略为：()()()()()rand1X tX ts x txt+=(13)即把式(13)换做式(7)。2.3.适应度函数适应度函数 为了使由算法得出的 PID 参数的调节效果能有一个衡量的标准，引入适应度函数的概念。适应度可用 PID 优化的评价指标中的误差绝对值对时间的积分的性能指标，并且加入了输入的平方项来抑制控制能量的增加 10。适应度函数可用下式表示。()()()20dJt e tutt=+(14)其中，式中：e(t)为输入与输出间的误差值，t，w 为权值常数，u(t)为控制量。考虑到超调量对于评价指标的影响，将超调量也加入评价指标之中，引入惩罚措施来避免产生超调 11。此时的适应度函数式如式(15)所示。()()()()()()200,df e tJt e tute tt=+(15)2.4.DMC 预测控制预测控制 由于线径制造时芯线从机头被推挤完毕后与测径仪的反馈之间有一定的间隔，从芯线出机头到被测径仪测量反馈至系统的这段时间为，由推挤速度和机头与测径仪之间的间距l决定，的表达式为wl v=，vw这里的是指线径生产的平均速度。由于系统之后时间的存在导致了系统的控制不是实时进行的，所以需要进行预测控制。本文采用 DMC 预测的方法来预测从机头产出的线径数据，在原来时滞的基础上提前 个时间将预测结果输入系统，大幅提高线径控制系统的性能。其算法的本质是从阶跃响应来调整增量控制，之后以预测模型来对当前时刻的和之前的偏差值为数据进行滚动优化，最后将最合适的结果输入至系统。2.4.1.预测模型预测模型 模型的建立是以被测对象的输入和输出的数据进行采集为基础的。第一步，对被测对象进行单位响应的采样，其值用()iaa iT=，1,2,3,i=?表示。第二步依据得到的采样数据建立模型的预测向量，用T1iNaaa=?来表示，其中 N 为建模过程中时域的长度。第三步，采取控制变量法，除了其余的条件不变，只在 K 时间点对对象施加一增量控制()u k，采集之后由对象反馈的 N 个输出的数据。那么在时间点之后对于对象的输出预测值可用如下描述：()()()0mykyka u k=+(16)()()()T1,mmmykykykp=+?(17)()()(),1,u ku k ku kMk=+?(18)其中，是 K 时间点起系统加入()u k之后的预测输出。表示系统未加()u k时的输出值；p 表示为系统的优化时域长度；M 表示控制时域长度；N 表示建模时域长度；(k+1,k)表示在 k 时对 k+1 的预测。根据对推挤机对象模型的单位阶跃响应进行参考，拟使用采样周期 T=0.5，N=20，M=1，P=15。2.4.2.在线校正在线校正 实际生产环境下，会存在干扰、建模不准确等因素的存在，得到的预测输出常常与真实的输出有一程元清，郭金妹 DOI:10.12677/jsta.2023.113027 253 传感器技术与应用 定的偏差。故必须对预测值进行在线校正以保证预测的拟合性，具体做法是在线校正。即在 k 的下一时刻把对象的实际输出值与(14)得出的预测输出值 ym(k)进行差值处理。()()()111,me ky kykk+=+(19)接下来用由加权系数组成的 N 维校正向量 h 以滚动优化的方式实施校正更新：()()()111pmykykhe k+=+(20)其中，T12,Nhh hh=?，这里的12130.9hhh=?是由本系统决定的。p(k+1)为校正后的预测输出值，在相位偏移过后其自身以下一时刻的初始预测值的身份进入系统之中。p(k+1)可用向量表示为：()()021ppyks yk+=+(21)其中 s 是位移矩阵，可表示为：0100001000010001s=?(22)3.系统仿真与结果分析系统仿真与结果分析 3.1.试验仿真数据试验仿真数据 此次仿真的实验选择在 MATLAB 的 Simulink 内实施。数据的选取方式是采用到赣州某电缆生产工厂的实际现场采集取得。实验的目的是比较 IWOA-DMC-PID 控制算法、WOA-DMC 控制算法与传统 PID算法的控制效果。仿真实验开始时，先对 PID 三个参数进行初始化，令 kp=ki=0.7，kd=2.1。传递函数的 T1=0.84，T2=0.46，时间滞后 =0.9。设定鲸鱼的最大迭代次数为 100，种群大小设置为 30。3.2.寻优迭代分析寻优迭代分析 Figure 3.Relationship between fitness value and iteration number 图图 3.适应值与迭代次数的关系 由图 3 可知，WOA 的迭代次数为 53 次左右的时候迭代完成，并且适应值在 101.66 左右，而 IWOA的迭代次数为 63 次数左右，并且相较于 WOA 的适应度值小，这说明了 WOA 算法虽然迭代次数少但是陷入了局部最优情况，没有得出全局最优的最佳一组参数，而 IWOA 则有效避免了这种局部最优解的情况，准确的找到了全局最优解。程元清，郭金妹 DOI:10.12677/jsta.2023.113027 254 传感器技术与应用 3.3.脉冲响应分析脉冲响应分析 为了比较三种控制方法的控制效果，设置线径值为 2 mm。当系统未受外界扰动时分别对三种控制方式进行仿真，三种方式的控制效果如图 4 所示。由图中曲线可以得到，采用传统的 PID 算法时的系统超调量 约为 27.5%，达到运行平稳所需的调节时间 ts约为 30 秒，调节过程中线径波动较大，响应较为迟钝，调节时间较为漫长。采用WOA-PID的结合算法，系统调节的效果有所提升，系统超调量约为15.4%，调节时间 ts约为 21 秒，整体过程的控制线径波动比传统 PID 调节效果稳定，但依然有很大提升空间。而采用 IWOA-DMC-PID 复合控制算法时，系统的超调量 约为 0，调节时间 ts为 10 秒，但系统的响应时间稍长。Figure 4.Three types of pulse response curves 图图 4.三种方式的脉冲响应曲线 3.4.抗扰动性能分析抗扰动性能分析 在三种控制方式达到稳定运行状态后，在 30 秒对系统施加一个单位为 1 的阶跃扰动，其三种控制方式的仿真曲线如图 5 所示：Figure 5.Three ways to add disturbance response curve 图图 5.三种方式加扰动响应曲线 由图 5 中的相应曲线可知，普通 PID 在扰动出现后可做出最快反应，但是其控制线径的波动也十分程元清，郭金妹 DOI:10.12677/jsta.2023.113027 255 传感器技术与应用 剧烈，到达稳定的调节时间为 23 秒。采用 WOA-PID 结合算法反应时间较传统 PID 的稍慢，但是线径的波动情况有所改善，但线径值仍有波动。采用 IWOA-DMC-PID 复合算法的反应时间与前两者相当，但是在线径波动和调节时间方面都有极大的改善，线径值程平滑改变的趋势，几乎无波动。这说明了三种方式中 IWOA-DMC-PID 具有良好的抗干扰性能。3.5.鲁棒性测试分析鲁棒性测试分析 由于在实际的工厂生产条件下，由于电磁干扰等因素的存在，会导致数学模型发生变化，而模型的变化会影响控制效果甚至对系统失去控制。所以要求整个系统具有较强的鲁棒性，即使系统参数或者系统特性因为外因发生变化时，系统也能有维持现状，也能保持较好控制效果的性能。本文将数学模型的增益与时滞各自增大 40%后以三种控制方式进行仿真，得到的响应结果如图 6 所示。Figure 6.System response diagram after parameter changes 图图 6.参数改变后的系统响应图 由图 6 可知，在数学模型改变后，普通 PID 的调节效果整体呈现出较大幅度的震荡，且调节时间最长，约为27 秒。WOA-PID的调节效果较传统PID有一定的提升，调节时间约为11 秒。最后IWOA-DMC-PID的调节效果对比前两者有大幅度提升，系统震荡最小，几乎无震荡，且调节速度很快，只需 7 秒即可达到稳定。这表明三种方式中 IWOA-DMC-PID 的鲁棒性最佳。4.结论结论 本文针对传统 PID 参数调节困难，控制效果不佳，智能算法与 PID 结合寻优易陷入局部最优的问题，本文引入鲸鱼算法与 DMC 与 PID 结合并应用于线径推挤系统。通过对鲸鱼算法进行改进，并且改进之后在 Matlab 中编程并用 SIMULINK 对 PID，WOA-PID 与 IWOA-DMC-PID 三种控制方式进行仿真，结果表明，采用 IWOA-DMC-PID 控制效果比其他两种控制方式的效果好。线径挤出系统能够将挤出的线径稳定控制在 3.5 mm 左右，并且对于外界的干扰有较好的快速响应的能力，具有较强的鲁棒性，可以满足实际线径挤出的生产要求。参考文献参考文献 1 周克良,董世镇,张自建.基于改进磷虾群预测 PID 的电缆线径控制J.传感器与微系统,2022,41(4):96-99,103.程元清，郭金妹 DOI:10.12677/jsta.2023.113027 256 传感器技术与应用 2 Hussein,H.(2014)Feedforward,Hot and Cold Control System for Controlling Plastic Insulation Thickness in Manu-facturing of Cable Products.Arabian Journal for Science&Engineering,39,1399-1408.https:/doi.org/10.1007/s13369-013-0723-x 3 Li,H.and Yu,Q.(2016)The Wire Beltline Diameter ACO-KF-PID Control Research.2016 Prognostics and System Health Management Conference,Chengdu,19-21 October 2016,1-6.https:/doi.org/10.1109/PHM.2016.7819758 4 袁宇鹏,胡杨,李军,等.基于动态矩阵控制的自主移动机器人模型预测控制方法研究J.重庆邮电大学学报(自然科学版),2018,30(4):537-543.5 周克良,胡梁眉,洪智慧.模糊分数阶 PID 在线径控制系统中的应用J.制造业自动化,2019,41(7):63-69.6 周克良,过振宇.基于模糊自适应 PID 的推挤机线径预测控制J.计算机仿真,2019,36(3):302-306.7 Mirjalili,S.and Lewis,A.(2016)The Whale Optimization Algorithm.Advances in Engineering Software,95,51-67.https:/doi.org/10.1016/j.advengsoft.2016.01.008 8 孙运.基于改进鲸鱼算法的逆变器控制及应用研究D:硕士学位论文.南昌:南昌大学,2021.9 王曦璐.改进鲸鱼算法在函数优化和片上网络中的应用D:硕士学位论文.西安:西安电子科技大学,2019.10 赵广元,王超.基于人群搜索算法的直流电机优化设计J.舰船电子工程,2019,39(10):242-246.11 汪谦松,夏链,韩江.改进人群搜索算法优化的直线电机模糊 PID 控制J.合肥工业大学学报(自然科学版),2022,45(4):458-463,474.