基于改进差分进化机理的动力定位节能推力分配方法组成比例

本发明专利技术提供一种基于改进差分进化机理的动力定位节能推力分配方法，其包括以下步骤：步骤S1：确定动力定位功率需求；步骤S2：建立动力定位功率分配数学模型；步骤S3：针对动力定位中的功率分配问题初始化解；步骤S4：根据突变常量λ来选择突变方式；步骤S5：进行突变操作；步骤S6：进行交叉操作；步骤S7：评价所有解的适应值，并记录最优解；步骤S8：判断是否满足终止条件，是则执行步骤S9，否则返回步骤S4；步骤S9：停止循环，输出最优解，根据最优解进行功率分配。通过本发明专利技术，将动力定位船舶所需推力合理分配到各推进器，在保证船舶安全的基础上，减少船舶能耗，减低船舶运营成本，降低船舶污染。

【技术实现步骤摘要】
基于改进差分进化机理的动力定位节能推力分配方法
本专利技术涉及船舶自动化领域，尤其涉及一种基于改进差分进化机理的动力定位节能推力分配方法。
技术介绍
船舶动力定位系统(DPS)的主要功能是保持船舶在期望的位置或驱动船舶沿着指定的轨迹航行。它主要根据船舶的当前的位置状态与期望值之间的偏差量和外部环境干扰，计算出船舶实时所需要的总的推力和力矩，然后再根据推力分配逻辑方法计算出船舶装备的各个推进器的所需产生的推力及其角度，并将控制指令传送到推进器，推进器根据该指令产生要求的推力完成船舶定位。因此，推力分配是DPS中重要的一环，它将控制器与推力系统联系起来形成一个整体。推进器通常情况下不少于5个，但并不是越多越好，还要必须考虑到两个推进器之间和推进器与船体之间的互相作用等因素，同时还要兼顾船舶所必需的操纵性。此时整个推进系统便有很多种能够满足定位要求的推力大小及方向的组合，而推力分配逻辑的主要功能就是找出一个最优秀的组合，并将其转化为控制指令实时的分配给各个推进器。推力分配优化问题是一个多约束优化问题，它要求推进器实时高效合理地产生船舶所需要的推力，并满足在保持一定程度的操纵性的同时尽量使能耗最小的要求。优秀的分配方法不仅能够提高船舶定位精度，而且还拥有降低能量消耗，减少机械磨损和噪声等功能。动力定位船舶一般装备的是全回转推力器，此类推进器及其推进系统一般拥要有以下限制条件：(1)受驱动机功率限制，每个推进器产生的推力都有一个最大值，推进器无法产生比最大值更大的推力。(2)推进器推力及其角度的变化速率有限制，无法短时间变化很大。(3)全回转推进器可绕其主轴360°转动，但又因为临近推进器之间的互相作用，存在推力禁区，因此需要对各推进器的转动角度设定一个变化范围。此外，还要考虑奇异结构的出现。对动力定位船舶来说，当其装备的所有推进器都拥有一样的推进方位角时，奇异结构就会发生。此时，船舶操纵性能就会大幅降低甚至失去，又由于推进器旋转速度的限制，整个系统将会较长的时间来恢复，这将有可能在相关的工程作业中产成严重的后果。因此，在船舶动力定位推力分配的研究与应用中，应当避免奇异结构的发生。因此为了保持能耗与操作性之间的平衡，提出一种折衷的方案，即根据实际的需求设置调节系数，来平衡船舶对操作性和能耗的需求。现有技术中，主要采用以下几种技术方案：专利《采用遗产算法的船舶动力定位系统推力分配方法》，申请号201110253236.2，采用遗传算法方式来进行推力的分配过程。专利《动力定位系统的推力分配方法》，申请号201310003584.3，采用低通滤波和序列二次规划法两种优化算法分布式计算推力器方位角和推力大小。专利《采用动态禁止角的动力定位推力分配装置及其分配方法》，申请号201210356674.6，采用二次规划算法来进行推力的分配过程。目前船舶动力定位推力分配的主要方法是忽略推力分配的能耗特性，难以找到推力分配的最优解。
技术实现思路
本专利技术的目的是根据所构建的动力定位推力分配数学模型，通过改进的差分进化机制，有效地找到分配各推进器推力和方向指令的最优解，在保证推进器合力满足控制要求的基础上，降低船舶推进器能耗。为实现上述目的，本专利技术采用以下技术方案：一种基于改进差分进化机理的动力定位节能推力分配方法，其包括以下步骤：步骤S1：确定动力定位功率需求；步骤S2：建立动力定位功率分配数学模型；步骤S3：针对动力定位中的功率分配问题初始化解；步骤S4：根据突变常量λ来选择突变方式；步骤S5：进行突变操作；步骤S6：进行交叉操作；步骤S7：评价所有解的适应值，并记录最优解；步骤S8：判断是否满足终止条件，是则执行步骤S9，否则返回步骤S4；步骤S9：停止循环，输出最优解，根据最优解进行功率分配。在本专利技术一实施例中，步骤S2包括以下具体步骤：步骤S21：关于推力分配问题的非线性优化数学模型如下：约束调节如下：s＝τ-B(α)FFmin≤F≤Fmax△Fmin≤F-F0≤△Fmaxαmin≤α≤αmax△αmin≤α-α0≤△αmax其中，第一项中W是总能耗，P为权值系数，Fi为第i个推进器的推力，ki为计算参数；第二项sTQs为惩罚项，s是广义推力误差向量；权值矩阵Q为对角正定矩阵，其应取大值，用以保证误差s趋近于零；第三项(α-α0)TΩ(α-α0)是用来约束推进角度的变化速度，其中α为本时刻的推进器的方位角，α0为前一时刻的推进器的方位角，权值矩阵Ω>0用来调节优化目标；第四项是用来避免奇异结构的，其中lxn和lyn分别是第n个推进器的X方向坐标和Y方向坐标；X方向为从船中到船首，Y方向为从船中到船右舷；如果推进系统奇异或趋近于奇异，即det(B(α)B′(α))等于零或近似于零，则第四项的值会很大，相当于惩罚函数；式中ε>0，δ大于零，δ是调节系数，用来平衡船舶的能耗与操纵性，δ取值越大操纵性越好，相应的能耗也会有所增加，δ取值越小，则情况正好相反；步骤S22：在约束条件中，τ＝(τx,τy,τM)T是期望得到的力和力矩，F是推进系统的推力矩阵，B(α)F则是推进系统实际产生的合推力及合力矩，由B(α)计算公式计算出实际与期望之间的广义推力误差向量；Fmax和Fmin分别表示推进器推力的最大值和最小值，限制了推进器的推力范围；△Fmax和△Fmin分别表示推进器推力在单位时间内变化的上下限，B(α)计算公式规定了推力变化率的范围；相应的，αmax和αmin为推进器旋转角度的范围，△αmax和△αmin是两时刻间推进器角度变化幅值的上下限。在本专利技术一实施例中，步骤S4包括以下步骤：步骤S41：在改进差分进化中，按照如下公式将人工蜂群算法加入差分进化算法的方式作为局部搜索策略；式中，xi为选中的进行邻域搜索的个体，xk为xi相邻的个体，x′i为邻域搜索后的得到的粒子，即x′i是xi通过和其相邻的随机个体xk相互比较后的修正解，是[-1,1]中随机生成的值；k∈{1,2,3,...,SN}其中SN是种群的数量；步骤S42：按照如下公式加入粒子群算法的社会认知部分来优化差分进化算法：x′i＝xr1+F(xr2-xr3)+ψ(xgbest-xr1)；其中ψ是范围[0,1]中随机生成的数值，xgbest是全局最优解；xr1，xr2和xr3是随机选取的变量且r1≠r2≠r3；F是用于控制在突变过程中的摄动总量和提高收敛速度的缩放因子，其取值范围为[0,1]；新解x′i由三个部分合成产生：第一部分被选取的目标个体x′i；第二部分是由随机选择的父代个体的差分产生的向量；最后一部分是根据由被选取的目标个体和当前整个种群中全局最优解的差值来进行差分操作生成的解向量；步骤S43：按照以下步骤将这步骤S41和步骤S42改进方式进行组合：如果rand<λ则利用否则按照x′i＝xr1+F(xr2-xr3)+ψ(xgbest-xr1)进行突变操作；rand为产生随机数的一个随机函数；其中λ是范围在[0,1]的控制变量。进一步的，按照以下公式来调节λ值来选择寻优策略：其中t是当前循环次数，T是选定的最大循环次数。在本专利技术一实施例中步骤S5包括以下步骤：分别利用公式ABC改进策略或PSO改进策略进行突变操作；本文档来自技高网...

【技术保护点】
1.一种基于改进差分进化机理的动力定位节能推力分配方法，其特征在于：包括以下步骤：步骤S1：确定动力定位功率需求；步骤S2：建立动力定位功率分配数学模型；步骤S3：针对动力定位中的功率分配问题初始化解；步骤S4：根据突变常量λ来选择突变方式；步骤S5：进行突变操作；步骤S6：进行交叉操作；步骤S7：评价所有解的适应值，并记录最优解；步骤S8：判断是否满足终止条件，是则执行步骤S9，否则返回步骤S4；步骤S9：停止循环，输出最优解，根据最优解进行功率分配。

【技术特征摘要】
1.一种基于改进差分进化机理的动力定位节能推力分配方法，其特征在于：包括以下步骤：步骤S1：确定动力定位功率需求；步骤S2：建立动力定位功率分配数学模型；步骤S3：针对动力定位中的功率分配问题初始化解；步骤S4：根据突变常量λ来选择突变方式；步骤S5：进行突变操作；步骤S6：进行交叉操作；步骤S7：评价所有解的适应值，并记录最优解；步骤S8：判断是否满足终止条件，是则执行步骤S9，否则返回步骤S4；步骤S9：停止循环，输出最优解，根据最优解进行功率分配。2.根据权利要求1所述的基于改进差分进化机理的动力定位节能推力分配方法，其特征在于：步骤S2包括以下具体步骤：步骤S21：关于推力分配问题的非线性优化数学模型如下：约束调节如下：s＝τ-B(α)FFmin≤F≤Fmax△Fmin≤F-F0≤△Fmaxαmin≤α≤αmax△αmin≤α-α0≤△αmax其中，第一项中W是总能耗，P为权值系数，Fi为第i个推进器的推力，ki为计算参数；第二项sTQs为惩罚项，s是广义推力误差向量；权值矩阵Q为对角正定矩阵，其应取大值，用以保证误差s趋近于零；第三项(α-α0)TΩ(α-α0)是用来约束推进角度的变化速度，其中α为本时刻的推进器的方位角，α0为前一时刻的推进器的方位角，权值矩阵Ω>0用来调节优化目标；第四项是用来避免奇异结构的，其中lxn和lyn分别是第n个推进器的X方向坐标和Y方向坐标；X方向为从船中到船首，Y方向为从船中到船右舷；如果推进系统奇异或趋近于奇异，即det(B(α)B′(α))等于零或近似于零，则第四项的值会很大，相当于惩罚函数；式中ε>0，δ大于零，δ是调节系数，用来平衡船舶的能耗与操纵性，δ取值越大操纵性越好，相应的能耗也会有所增加，δ取值越小，则情况正好相反；步骤S22：在约束条件中，τ＝(τx,τy,τM)T是期望得到的力和力矩，F是推进系统的推力矩阵，B(α)F则是推进系统实际产生的合推力及合力矩，由B(α)计算公式计算出实际与期望之间的广义推力误差向量；Fmax和Fmin分别表示推进器推力的最大值和最小值，限制了推进器的推力范围；△Fmax和△Fmin分别表示推进器推力在单位时间内变化的上下限，B(α)计算公式规定了推力变化率的范围；相应的，αmax和αmin为推进器旋转角度的范围，△αmax和△αmin是两时刻间推进器角度变化幅值的上下限。3.根据权利要求1所述的基于改进差分进化机理的动力定位节能推力分配方法，其特征在于：步...

【专利技术属性】
技术研发人员：吴德烽，赵珂，顾佳栋，
申请(专利权)人：集美大学，
类型：发明
国别省市：福建,35