CN109885061A - 一种基于改进nsga-ⅱ的动力定位多目标优化方法 - Google Patents
一种基于改进nsga-ⅱ的动力定位多目标优化方法 Download PDFInfo
- Publication number
- CN109885061A CN109885061A CN201910194185.7A CN201910194185A CN109885061A CN 109885061 A CN109885061 A CN 109885061A CN 201910194185 A CN201910194185 A CN 201910194185A CN 109885061 A CN109885061 A CN 109885061A
- Authority
- CN
- China
- Prior art keywords
- ship
- dynamic
- dynamic positioning
- solution
- nsga
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 61
- 238000005457 optimization Methods 0.000 claims abstract description 45
- 238000005265 energy consumption Methods 0.000 claims abstract description 32
- 230000000694 effects Effects 0.000 claims abstract description 20
- 239000011159 matrix material Substances 0.000 claims description 13
- 230000008569 process Effects 0.000 claims description 8
- 238000013016 damping Methods 0.000 claims description 6
- 239000012530 fluid Substances 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 230000002194 synthesizing effect Effects 0.000 claims description 2
- 230000007613 environmental effect Effects 0.000 abstract description 7
- 238000004422 calculation algorithm Methods 0.000 description 33
- 238000004364 calculation method Methods 0.000 description 5
- 230000002068 genetic effect Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 239000000446 fuel Substances 0.000 description 4
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 4
- 230000009471 action Effects 0.000 description 2
- 238000004873 anchoring Methods 0.000 description 2
- 238000004880 explosion Methods 0.000 description 2
- 239000003345 natural gas Substances 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Landscapes
- Feedback Control In General (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明属于船舶动力定位控制领域,具体涉及一种基于改进NSGA‑Ⅱ的动力定位多目标优化方法,包括以下步骤:建立动力定位船的三自由度动力学模型和运动学模型;利用动态面控制方法设计船舶动力定位控制律;建立动力定位船的能源消耗模型;选取优化目标函数,利用改进的NSGA‑Ⅱ调节动态面控制器的控制参数,对控制器效果及船舶能源消耗进行多目标问题优化。本发明考虑了动力定位改变自身位置及姿态过程中的能源消耗问题,利用改进的NSGA‑Ⅱ整定控制器参数,减小推进器的能耗,增强环保性和经济性。
Description
技术领域
本发明属于船舶动力定位控制领域,具体涉及一种基于改进NSGA-Ⅱ的动力定位多目标 优化方法。
背景技术
随着人们对于资源需求的不断增加,陆地石油及天然气资源日渐枯竭,海洋资源成为了 探索和开发的重点。动力定位作为深海开发必不可少的关键性技术受到了重视。动力定位, 是指利用船舶自身推进器产生的推力,抵消外部环境如:风、浪、流产生的干扰,维持船舶 的位置和艏向,使得船舶能够在指定位置范围内工作。与传统的锚泊系统相比,动力定位技 术具有不受水深限制,操作灵活,定位精度高的优点。但是,由于推进器在调整船舶位置和 艏向的过程中会消耗燃料,存在经济性和环保性问题。
本发明采用改进NSGA-Ⅱ的智能控制算法,对于动力定位动态面控制器的控制效果以及 船舶的燃料消耗两方面问题进行优化。目前,国内外文献未见利用本方法的动力定位多目标 问题优化方法。
发明内容
本发明的目的在于提供一种基于NSGA-Ⅱ的动力定位多目标优化方法,它解决了控制器 控制效果和船舶能源消耗的多目标优化问题。
一种基于改进NSGA-Ⅱ的动力定位多目标优化方法,包括以下步骤:
(1)建立动力定位船的三自由度动力学模型和运动学模型;
(2)利用动态面控制方法设计船舶动力定位控制律;
(3)建立动力定位船的能源消耗模型;
(4)选取优化目标函数,利用改进的NSGA-Ⅱ调节动态面控制器的控制参数,对控制器 效果及船舶能源消耗进行多目标问题优化。
所述建立动力定位船的三自由度动力学模型和运动学模型,包括:
动力定位船的三自由度动力学模型和运动学模型为以下形式:
其中,η=[x,y,ψ]T表示在大地坐标系下船的位置(x,y)和艏向ψ;υ=[u,v,r]T表示船 体坐标系下的纵荡、横荡方向上的线速度u,v和艏摇方向上的角速度r;M为系统惯性矩阵; D(υ)为水动力阻尼系数矩阵;τ为船舶推进系统产生的控制输入量;w为外界环境干扰产生 的力和力矩;
其中,m为船舶质量,Iz为船舶的转动惯量;xG为重心G在船体坐标系下的位置;为船舶在纵荡、横荡、艏摇方向上产生的附加质量;Xu,Yv,Yr,Nv,Nr为船舶三个方向上的线性阻尼值。
所述利用动态面控制方法设计船舶动力定位控制律,包括:
建立第一个船舶位置姿态误差动态面:
s1=η-ηd
其中:ηd为期望的船舶位置姿态,对s1求导得:
则υ的虚拟反馈控制率为:
φ1=-R-1(ψ)k1s1
其中,k1∈R3×3为正定的对角矩阵;
为φ1设计一个一阶低通滤波器:
其中,Td∈R3×3为滤波器的时间常数矩阵,Xd为滤波器的输出值;
建立第二个船舶速度误差动态面:
s2=υ-Xd
综合上述公式及船舶模型,得到:
船舶动力定位的动态面控制器的控制律设计如下:
其中,k1,k2,Td为动态面控制器控制参数;
将设计的控制律τ代入动力定位船的三自由度动力学和运动学模型,得到:
由此,得出每一时刻动力定位船的状态η(k1,k2,Td,t)、υ(k1,k2,Td,t)。
所述建立动力定位船的能源消耗模型,包括:
船体的螺旋桨推进力为:
FT=(1-tp0)ρn2Dp 4KT(J)
其中,tp0为推力减额系数;ρ表示流体密度,n为螺旋桨转速;Dp为螺旋桨直径;KT(J) 为无因次推力系数且满足下式:
KT(J)=K0+K1J+K2J2
其中,是进速比,V=[u,v]T是船速,w是螺旋桨伴流系数,K0,K1,K2是根据 螺旋桨的无因次推力性能曲线得到的无因次常系数;
船舶的推进系统的推力为:
FT=AV2+BVn+Cn2
其中,A,B,C为有因次系数,且满足
A=(1-tp0)(1-wp)2ρD2 pK2
B=(1-tp0)(1-wp)2ρD3 pK1
C=(1-tp0)ρD4 pK0
实际的螺旋桨推进器的功率消耗为:
其中,KQ,KT分别表示螺旋桨推进器的推力系数和转矩系数。
所述选取优化目标函数,利用改进的NSGA-Ⅱ调节动态面控制器的控制参数,对控制器 效果及船舶能源消耗进行多目标问题优化,包括:
选取的两个优化目标函数如下:
f1(X)=mine=min{η(X)-ηd}
其中:X=(x1,x2,x3)=(k1,k2,Td);
改进的NSGA-II新生成解的接受概率为:
其中,m为优化的目标数,Ri为第i个目标的范围;
依据接受概率的局部搜索过程包括:选择要进行局部搜索的个体x,选择单位长度随机 向量ν、指定的步长λ遍历周围的解x″=x′+λ·v,x′=x,判断x′和x″的支配关系;
若x″支配x′,则x″解较优,令x′=x″,λ=λ/2,缩小步长,进一步在x″附近搜索更优解;
若二者不存在支配关系,属于同一层级,则把x″作为非支配解加入外部种群;
若x′支配x″,需要依据概率接受x″,计算方式如下:
若任意生成的rand(·)∈(0,1),使得rand(·)≤prob,则依据概率选择接受解x″,缩小步长 λ=λ/2;反之,则丢弃这个解,继续下一步搜索,直到满足所规定的总步长,即将指定解的 周围搜索完毕,局部搜索结束。
本发明的有益效果在于:
1、本发明应用改进NSGA-Ⅱ方法对动态面控制器参数进行调整,解决了控制器参数多靠人工 经验选取,易造成适应性差、控制效果不佳的问题,提高了优化性能,缩短整定时间,增强对控制 目标的适应性。
2、本发明考虑了动力定位改变自身位置及姿态过程中的能源消耗问题,利用改进的NSGA-Ⅱ 整定控制器参数,减小推进器的能耗,增强环保性和经济性。
3、本发明不同于传统的多目标优化问题解决方法,即选取权重将多目标问题转化为单目标问 题解决。提出了利用改进的NSGA-Ⅱ的方法,直接优化多个目标,给出Pareto最优解,并可结合 实际情况,选出最适应需求的解。
4.本发明提出了改进的NSGA-Ⅱ方法,增加外部种群,记录存储非支配解,增加了种群的丰富 性;引入依据接受概率的局部搜索方法,在一定程度上解决了算法陷入局部最优的问题,提高了 算法的搜索能力和性能。
附图说明
图1是船体坐标系和大地固定坐标系;
图2是NSGA-Ⅱ的算法流程图;
具体实施方式
下面结合附图对本发明做进一步描述。
随着人们对于资源需求的不断增加,陆地石油及天然气资源日渐枯竭,海洋资源成为了 探索和开发的重点。动力定位作为深海开发必不可少的关键性技术受到了重视。动力定位, 是指利用船舶自身推进器产生的推力,抵消外部环境如:风、浪、流产生的干扰,维持船舶 的位置和艏向,使得船舶能够在指定位置范围内工作。与传统的锚泊系统相比,动力定位技 术具有不受水深限制,操作灵活,定位精度高的优点。但是,由于推进器在调整船舶位置和 艏向的过程中会消耗燃料,存在经济性和环保性问题。
本发明采用改进NSGA-Ⅱ的智能控制算法,对于动力定位动态面控制器的控制效果以及 船舶的燃料消耗两方面问题进行优化。目前,国内外文献未见利用本方法的动力定位多目标 问题优化方法。
1、一种基于改进NSGA-Ⅱ的动力定位多目标优化方法,该方法针对动力定位控制器控制 效果和系统的能源消耗两个目标问题,采用了改进的NSGA-Ⅱ的算法调整动力定位动态面控 制器参数,优化控制效果,同时降低能源消耗。首先
建立动力定位船的三自由度动力学模型和运动学模型,
利用动态面控制方法设计船舶动力定位控制律,得到每一时刻下船舶位置姿态与期望的 位置姿态的误差;
其次建立动力定位螺旋桨推进系统能源消耗模型,计算在不同时刻的船速下推进系统的 能源消耗功率;
最后选取船舶位置姿态误差和推进系统功率消耗作为优化目标函数,两个优化目标函数 的自变量是动态面控制器的控制参数。引入带有外部种群和依据接受概率的局部搜索的改进 NSGA-Ⅱ算法,针对上述的两个目标函数进行多目标优化,即对动态面控制器的参数进行整 定,能够得到Pareto最优解,根据实际情况,选择相应参数,即在可接受的误差范围内,选 取使得能源消耗小的解。
选取能够解决反步控制微分爆炸问题的动态面控制方法,来设计动力定位控制器。动力 定位的动态面控制律设计为:
式中:k1,k2,Td为动态面控制器控制参数。
将设计的控制律τ代入动力定位船的三自由度动力学和运动学模型,可以得到
根据上述公式,可以得出每一时刻动力定位船的状态η(k1,k2,Td,t)、υ(k1,k2,Td,t),并且 它们与动态面控制器参数有关。
选取船舶位置姿态误差和推进系统能源消耗作为优化目标函数,利用NSGA-Ⅱ的方法对 动态面控制器参数进行整定,可以得到Pareto最优解,并根据实际情况,选取所需的解。
选取的两个优化目标函数如下:
f1(X)=mine=min{η(X)-ηd}
式中:X=(x1,x2,x3)=(k1,k2,Td)
在NSGA-Ⅱ的基础上,提出了增加种群丰富性和多样性的方法。在算法中设立两个种群, 外部种群和内部种群,内部种群在算法的作用下进行选择、交叉、变异,外部种群用来记录、 存储通过快速非支配排序得到的非支配解。当内部种群趋于收敛时,由外部种群中的个体随 机生成来对原始种群进行更新,即二次初始化。同时,对外部种群的规模进行限制,在外部 种群个体数达到最大时,用拥挤度算子来优先删除稠密区域的解。这样增加了种群的多样性, 提高了算法的搜索能力和性能。
针对NSGA-Ⅱ易陷入局部最优解的问题,在遗传算法的选择操作中引入了依据接受概率 的局部搜索方法。提高了NSGA-II的局部搜索能力,增加最优解的收敛速度和精度。下面给 出新生成解的接受概率。
式中,m为优化的目标数,在这里取m=2,Ri为第i个目标的范围。
依据接受概率的局部搜索过程具体描述如下:选择要进行局部搜索的个体x,选择单位 长度随机向量ν、指定的步长λ遍历周围的解x″=x′+λ·v,x′=x,判断x′和x″的支配关系。
若x″支配x′,则x″解较优,令x′=x″,λ=λ/2,缩小步长,进一步在x″附近搜索更优解;
若二者不存在支配关系,属于同一层级,则把x″作为非支配解加入外部种群;
若x′支配x″,需要依据概率接受x″,计算方式如下:
若任意生成的rand(·)∈(0,1),使得rand(·)≤prob,则依据概率选择接受解x″,缩小步长 λ=λ/2;反之,则丢弃这个解,继续下一步搜索,直到满足所规定的总步长,即将指定解的 周围搜索完毕,局部搜索结束。
本发明的目的按以下步骤实现:
1.建立动力定位船的三自由度动力学模型和运动学模型
按照图1建立大地坐标系NOE和船体坐标系xbO′yb。
动力定位船的三自由度动力学模型和运动学模型为以下形式:
式中,η=[x,y,ψ]T表示在大地坐标系下船的位置(x,y)和艏向ψ;υ=[u,v,r]T表示船体 坐标系下的纵荡、横荡方向上的线速度u,v和艏摇方向上的角速度r;M为系统惯性矩阵; D(υ)为水动力阻尼系数矩阵;τ为船舶推进系统产生的控制输入量;w为外界环境干扰产生 的力和力矩。
式中,m为船舶质量,Iz为船舶的转动惯量;xG为重心G在船体坐标系下的位置;为船舶在纵荡、横荡、艏摇方向上产生的附加质量;Xu,Yv,Yr,Nv,Nr为船舶三个方向上的线性阻尼值。
2.利用动态面控制方法设计船舶动力定位控制律
动态面控制方法是基于反步控制发展起来的,反步控制在对虚拟控制量求导的过程中容 易引起微分爆炸的情况,动态面控制引入一阶积分滤波器来计算虚拟控制量的导数,消除了 微分项的膨胀,减少了计算的复杂性。这里为动力定位设计动态面控制器。
(1)建立第一个船舶位置姿态误差动态面
s1=η-ηd
式中:ηd为期望的船舶位置姿态,对s1求导可得:
则υ的虚拟反馈控制率可以设计为:
φ1=-R-1(ψ)k1s1
式中:k1∈R3×3为正定的对角矩阵。为φ1设计一个一阶低通滤波器。
式中:Td∈R3×3为滤波器的时间常数矩阵。Xd为滤波器的输出值。
(2)建立第二个船舶速度误差动态面
s2=υ-Xd
综合上述公式及船舶模型,可以得到:
由此,船舶动力定位的动态面控制器的控制律设计如下:
式中:k2∈R3×3为正定的对角矩阵。
将上面的公式代入计算,可得关于控制参数k1,k2,Td的控制律:
τ=D(υ)υ-w-k2(υ-Xd)+MTd -1(-R-1(ψ)k1(η-ηd)-Xd)
将设计的控制律τ代入动力定位船的三自由度动力学和运动学模型,可以得到
根据上述公式,可以得出每一时刻动力定位船的状态η(k1,k2,Td,t)、υ(k1,k2,Td,t),并且 它们与动态面控制器参数有关。
3.建立动力定位船的能源消耗模型
计算船体的螺旋桨推进力为:
FT=(1-tp0)ρn2Dp 4KT(J)
式中,tp0为推力减额系数;ρ表示流体密度,n为螺旋桨转速;Dp为螺旋桨直径;KT(J) 为无因次推力系数且满足
KT(J)=K0+K1J+K2J2
式中:是进速比,V=[u,v]T是船速,w是螺旋桨伴流系数,K0,K1,K2是根据 螺旋桨的无因次推力性能曲线得到的无因次常系数。由此,船舶的推进系统的推力可以写成:
FT=AV2+BVn+Cn2
式中:A,B,C为有因次系数,且满足
A=(1-tp0)(1-wp)2ρD2 pK2
B=(1-tp0)(1-wp)2ρD3 pK1
C=(1-tp0)ρD4 pK0
实际的螺旋桨推进器的功率消耗为:
式中:KQ,KT分别表示螺旋桨推进器的推力系数和转矩系数。
4.选取优化目标函数,利用改进的NSGA-Ⅱ调节动态面控制器参数,对控制器效果及船 舶能源消耗进行多目标问题优化
本发明欲达到的效果是使得控制器的控制效果精确的同时保证系统能源消耗低,但由于 控制效果精确就要控制器不停地调整船舶的位置及姿态来抵御环境力及外界干扰,这就使系 统的能源消耗增加,这是两个相互矛盾的多目标优化问题。改进的NSGA-Ⅱ算法是用来解决 这样的多目标优化问题,能够得到Pareto最优解,并根据需求选取Pareto前沿面上的点作为 所求解,即可以在误差可接受的范围内选取能源消耗小的方案。
(1)选取优化目标函数
本发明利用改进的NSGA-Ⅱ算法调整动态面控制器的参数,使得控制效果精确,控制误 差减小,同时,使得推进器能源消耗低。因此,需要有两个值来量化、评测这两个目标,这 里选取位置姿态误差和推进器功率消耗。
位置姿态误差
评价控制器的控制效果,将动态面控制器控制后输出的船舶位置姿态与期望的船舶位置 姿态相比较,用二者之间的误差来描述控制器的控制效果。选取船舶的位置姿态误差作为第 一个目标函数:
f1(X)=mine=min{η(X)-ηd}
推进器功率消耗
选取推进器功率消耗来评价推进系统的能源消耗,为降低系统所消耗的能源,应使得推 进器功率消耗最小,选取的第二个目标函数如下:
式中:X=(x1,x2,x3)=(k1,k2,Td)
(2)改进的NSGA-Ⅱ算法
本文选取的用来解决多目标优化问题的NSGA-Ⅱ算法,是在NSGA的基础上提出的,它 比NSGA算法更加优越:它采用了快速非支配排序算法,解决了NSGA一直被人诟病的计算复杂性问题;利用计算较为简便的拥挤度和拥挤度比较算子作为快速非支配排序分层后同一 层级中比较的标准,代替了需要指定的共享半径,使准Pareto域中的个体能均匀分布到整个 Pareto域上去,保持了种群的多样性;引入了精英策略,扩大了采样空间,防止最佳个体的丢 失,提高了算法的运算速度和鲁棒性。
本发明在NSGA-Ⅱ算法的基础上进行了改进,提出了增加外部种群进行二次初始化以及 依据接受概率的局部搜索方法。通过这两处改进,增加了种群的多样性,提高了算法的搜索 能力,解决了NSGA-Ⅱ算法易陷入局部最优的问题,提高了算法的性能。
NSGA-Ⅱ算法的流程图如附图2。具体执行步骤如下:
步骤一:初始化种群
根据所选的优化目标函数和约束条件,随机产生规模为N的初始种群,快速非支配排序 后,增加外部种群用来记录存储排序后的非支配解,以备稍后用来进行二次初始化。通过遗 传算法的选择、交叉、变异三个基本操作得到内部种群的第一代子代种群。该算法的优越处 在于在选择算子执行前,根据个体之间的支配关系进行分层,可以加快后续操作的过程。同 时增加了外部种群记录非支配解,以待内部种群收敛时进行二次初始化,对内部种群进行更 新,增加了种群多样性,防止陷入局部最优解。
快速非支配排序
初始化种群为N,根据两个优化目标函数,能够得到个体之间的支配关系,计算每个个 体i的两个参数np和Sp,即种群中支配个体i的个体数和个体i支配种群中其他个体的个体数。 算法的主要步骤是:找到种群中所有np=0的个体,并保存在当前集合F中;对于当前集合F 中的每个个体,其所支配的个体集合为Si,遍历Si中的每个个体l,执行nl=nl-1,如果nl=0 则将个体l保存在集合H中;记F中得到的个体为第一个非支配层的个体,并以H作为当前 集合,重复上述操作,直至整个种群被分层
外部种群
本发明在基础的算法上进行了改进,增加了一组外部种群,用来记录存储通过快速非支 配排序后得到的非支配解,内部种群在遗传算子的作用下进行选择、遗传、变异,在其趋于 收敛时,由外部种群中的个体随机生成来对内部种群进行更新,即二次初始化。同时,对外 部种群的规模进行限制,在外部种群个体数达到最大时,用拥挤度算子来优先删除稠密区域 的解。这样增加了种群多样性,提高了算法的搜索能力和性能。
步骤二:从第二代开始,将父代种群与子代种群合并,进行快速非支配排序,同时对每 个非支配层中的个体进行拥挤度计算,根据非支配关系以及个体的拥挤度选取合适的个体组 成新的父代种群;
拥挤度计算及拥挤度比较算子
计算拥挤度nd,对于每个目标函数,基于该目标函数对种群进行排序,令边界的两个个 体的拥挤度为无穷,计算nd=nd+(fm(i+1)-fm(i-1)),n=2,3,…,N-1。通过快速非支配排序 及拥挤度计算后,种群的每个个体都有非支配序rank和拥挤度nd两个属性,可以区分种群中 任意两个个体的支配和非支配关系。定义拥挤度比较算子≥n,个体的优劣比较依据为:i≥n j, 即个体i优于个体j,当且仅当irank<jrank或irank=jrank且id>jd
依据接受概率的局部搜索方法
在利用遗传算法进行选择操作的时候,为了提高NSGA-II的局部搜索能力,增加最优解 的收敛速度和精度,这里引入了一种依据概率的局部搜索方法。下面给出新生成解的接受概 率。
式中,m为优化的目标数,在这里取m=2,Ri为第i个目标的范围。
依据接受概率的局部搜索过程具体描述如下:选择要进行局部搜索的个体x,选择单位 长度随机向量ν、指定的步长λ遍历周围的解x″=x′+λ·v,x′=x,判断x′和x″的支配关系。
若x″支配x′,则x″解较优,令x′=x″,λ=λ/2,缩小步长,进一步在x″附近搜索更优解;
若二者不存在支配关系,属于同一层级,则把x″作为非支配解加入外部种群;
若x′支配x″,需要依据概率接受x″,计算方式如下:
若任意生成的rand(·)∈(0,1),使得rand(·)≤prob,则依据概率选择接受解x″,缩小步长 λ=λ/2;反之,则丢弃这个解,继续进行下一步搜索,直到满足所规定的总步长,即将指定 解的周围搜索完毕,局部搜索结束。通过这种依据接受概率的局部搜索方式加快了最优解的 收敛速度,增加了最优解的精度,避免算法陷入局部最优解。
步骤三:通过遗传算法的基本操作产生新的子代种群:依此类推,直到满足程序结束的条件, 输出最优解。
本发明的目的是这样实现的:
1.建立动力定位船的三自由度动力学模型和运动学模型
2.利用动态面控制方法设计船舶动力定位控制律
3.建立动力定位船的能源消耗模型
4.选取优化目标函数,利用改进的NSGA-Ⅱ调节动态面控制器的控制参数,对控制器效 果及船舶能源消耗进行多目标问题优化。
1、本发明应用改进NSGA-Ⅱ方法对动态面控制器参数进行调整,解决了控制器参数多靠人工 经验选取,易造成适应性差、控制效果不佳的问题,提高了优化性能,缩短整定时间,增强对控制 目标的适应性。
2、本发明考虑了动力定位改变自身位置及姿态过程中的能源消耗问题,利用改进的NSGA-Ⅱ 整定控制器参数,减小推进器的能耗,增强环保性和经济性。
3、本发明不同于传统的多目标优化问题解决方法,即选取权重将多目标问题转化为单目标问 题解决。提出了利用改进的NSGA-Ⅱ的方法,直接优化多个目标,给出Pareto最优解,并可结合 实际情况,选出最适应需求的解。
4.本发明提出了改进的NSGA-Ⅱ方法,增加外部种群,记录存储非支配解,增加了种群的丰富 性;引入依据接受概率的局部搜索方法,在一定程度上解决了算法陷入局部最优的问题,提高了 算法的搜索能力和性能。
Claims (5)
1.一种基于改进NSGA-Ⅱ的动力定位多目标优化方法,其特征在于,包括以下步骤:
(1)建立动力定位船的三自由度动力学模型和运动学模型;
(2)利用动态面控制方法设计船舶动力定位控制律;
(3)建立动力定位船的能源消耗模型;
(4)选取优化目标函数,利用改进的NSGA-Ⅱ调节动态面控制器的控制参数,对控制器效果及船舶能源消耗进行多目标问题优化。
2.根据权利要求1所述的一种基于改进NSGA-Ⅱ的动力定位多目标优化方法,其特征在于,所述建立动力定位船的三自由度动力学模型和运动学模型,包括:
动力定位船的三自由度动力学模型和运动学模型为以下形式:
其中,η=[x,y,ψ]T表示在大地坐标系下船的位置(x,y)和艏向ψ;υ=[u,v,r]T表示船体坐标系下的纵荡、横荡方向上的线速度u,v和艏摇方向上的角速度r;M为系统惯性矩阵;D(υ)为水动力阻尼系数矩阵;τ为船舶推进系统产生的控制输入量;w为外界环境干扰产生的力和力矩;
其中,m为船舶质量,Iz为船舶的转动惯量;xG为重心G在船体坐标系下的位置;为船舶在纵荡、横荡、艏摇方向上产生的附加质量;Xu,Yv,Yr,Nv,Nr为船舶三个方向上的线性阻尼值。
3.根据权利要求1所述的一种基于改进NSGA-Ⅱ的动力定位多目标优化方法,其特征在于,所述利用动态面控制方法设计船舶动力定位控制律,包括:
建立第一个船舶位置姿态误差动态面:
s1=η-ηd
其中:ηd为期望的船舶位置姿态,对s1求导得:
则υ的虚拟反馈控制率为:
φ1=-R-1(ψ)k1s1
其中,k1∈R3×3为正定的对角矩阵;
为φ1设计一个一阶低通滤波器:
其中,Td∈R3×3为滤波器的时间常数矩阵,Xd为滤波器的输出值;
建立第二个船舶速度误差动态面:
s2=υ-Xd
综合上述公式及船舶模型,得到:
船舶动力定位的动态面控制器的控制律设计如下:
其中,k1,k2,Td为动态面控制器控制参数;
将设计的控制律τ代入动力定位船的三自由度动力学和运动学模型,得到:
由此,得出每一时刻动力定位船的状态η(k1,k2,Td,t)、υ(k1,k2,Td,t)。
4.根据权利要求1所述的一种基于改进NSGA-Ⅱ的动力定位多目标优化方法,其特征在于,所述建立动力定位船的能源消耗模型,包括:
船体的螺旋桨推进力为:
其中,为推力减额系数;ρ表示流体密度,n为螺旋桨转速;Dp为螺旋桨直径;KT(J) 为无因次推力系数且满足下式:
KT(J)=K0+K1J+K2J2
其中,是进速比,V=[u,v]T是船速,w是螺旋桨伴流系数,K0,K1,K2是根据螺旋桨的无因次推力性能曲线得到的无因次常系数;
船舶的推进系统的推力为:
FT=AV2+BVn+Cn2
其中,A,B,C为有因次系数,且满足
实际的螺旋桨推进器的功率消耗为:
其中,KQ,KT分别表示螺旋桨推进器的推力系数和转矩系数。
5.根据权利要求1所述的一种基于改进NSGA-Ⅱ的动力定位多目标优化方法,其特征在于,所述选取优化目标函数,利用改进的NSGA-Ⅱ调节动态面控制器的控制参数,对控制器效果及船舶能源消耗进行多目标问题优化,包括:
选取的两个优化目标函数如下:
f1(X)=mine=min{η(X)-ηd}
其中:X=(x1,x2,x3)=(k1,k2,Td);
改进的NSGA-II新生成解的接受概率为:
其中,m为优化的目标数,Ri为第i个目标的范围;
依据接受概率的局部搜索过程包括:选择要进行局部搜索的个体x,选择单位长度随机向量ν、指定的步长λ遍历周围的解x″=x′+λ·v,x′=x,判断x′和x″的支配关系;
若x″支配x′,则x″解较优,令x′=x″,λ=λ/2,缩小步长,进一步在x″附近搜索更优解;
若二者不存在支配关系,属于同一层级,则把x″作为非支配解加入外部种群;
若x′支配x″,需要依据概率接受x″,计算方式如下:
若任意生成的rand(·)∈(0,1),使得rand(·)≤prob,则依据概率选择接受解x″,缩小步长λ=λ/2;反之,则丢弃这个解,继续下一步搜索,直到满足所规定的总步长,即将指定解的周围搜索完毕,局部搜索结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910194185.7A CN109885061B (zh) | 2019-03-14 | 2019-03-14 | 一种基于改进nsga-ⅱ的动力定位多目标优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910194185.7A CN109885061B (zh) | 2019-03-14 | 2019-03-14 | 一种基于改进nsga-ⅱ的动力定位多目标优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109885061A true CN109885061A (zh) | 2019-06-14 |
CN109885061B CN109885061B (zh) | 2021-11-23 |
Family
ID=66932324
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910194185.7A Active CN109885061B (zh) | 2019-03-14 | 2019-03-14 | 一种基于改进nsga-ⅱ的动力定位多目标优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109885061B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110794843A (zh) * | 2019-11-15 | 2020-02-14 | 山东交通学院 | 基于观测器的非线性船舶时滞动力定位船鲁棒镇定系统 |
CN110909950A (zh) * | 2019-11-29 | 2020-03-24 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种采用非劣排序算法进行燃油消耗优化方法及装置 |
CN111930123A (zh) * | 2020-08-13 | 2020-11-13 | 上海船舶研究设计院(中国船舶工业集团公司第六0四研究院) | 多目标综合优化决策方法、装置及电子设备 |
CN114019799A (zh) * | 2021-11-06 | 2022-02-08 | 易如(山东)智能科技有限公司 | 基于多目标遗传算法优化拖拉机燃油经济及动力性的方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102708406A (zh) * | 2012-05-10 | 2012-10-03 | 湖北省电力公司 | 一种基于多目标遗传算法的调度图优化方法 |
CN103580020A (zh) * | 2013-03-07 | 2014-02-12 | 长沙理工大学 | 一种基于NSGA-II和Look-ahead的含风电场电力系统多目标动态优化调度方法 |
EP2799944A1 (en) * | 2013-04-29 | 2014-11-05 | Siemens Aktiengesellschaft | A method for controlling a power generation plant |
CN104156944A (zh) * | 2014-07-15 | 2014-11-19 | 西安电子科技大学 | 基于nsga-ⅱ进化算法的变化检测方法 |
US20150032306A1 (en) * | 2012-03-22 | 2015-01-29 | Kongsberg Maritime As | Dynamic load compensation |
CN105911867A (zh) * | 2016-06-16 | 2016-08-31 | 哈尔滨工程大学 | 基于nsga-ii算法的船舶推力分配方法 |
CN105929825A (zh) * | 2016-05-16 | 2016-09-07 | 哈尔滨工程大学 | 一种基于神经动态模型的船舶动力定位反步控制方法 |
CN106055830A (zh) * | 2016-06-20 | 2016-10-26 | 吉林大学 | 基于动态规划的phev控制门限参数优化方法 |
CN107330131A (zh) * | 2016-04-29 | 2017-11-07 | 南京理工大学 | 机械零部件结构尺寸参数及其尺寸公差的区间优化方法 |
CN108319138A (zh) * | 2018-01-29 | 2018-07-24 | 哈尔滨工程大学 | 一种欠驱动无人艇的滑模-反步双回路轨迹跟踪控制方法 |
CN108490770A (zh) * | 2018-02-28 | 2018-09-04 | 哈尔滨工程大学 | 一种基于混合算法的船舶动力定位系统推力分配方法 |
-
2019
- 2019-03-14 CN CN201910194185.7A patent/CN109885061B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150032306A1 (en) * | 2012-03-22 | 2015-01-29 | Kongsberg Maritime As | Dynamic load compensation |
CN102708406A (zh) * | 2012-05-10 | 2012-10-03 | 湖北省电力公司 | 一种基于多目标遗传算法的调度图优化方法 |
CN103580020A (zh) * | 2013-03-07 | 2014-02-12 | 长沙理工大学 | 一种基于NSGA-II和Look-ahead的含风电场电力系统多目标动态优化调度方法 |
EP2799944A1 (en) * | 2013-04-29 | 2014-11-05 | Siemens Aktiengesellschaft | A method for controlling a power generation plant |
CN104156944A (zh) * | 2014-07-15 | 2014-11-19 | 西安电子科技大学 | 基于nsga-ⅱ进化算法的变化检测方法 |
CN107330131A (zh) * | 2016-04-29 | 2017-11-07 | 南京理工大学 | 机械零部件结构尺寸参数及其尺寸公差的区间优化方法 |
CN105929825A (zh) * | 2016-05-16 | 2016-09-07 | 哈尔滨工程大学 | 一种基于神经动态模型的船舶动力定位反步控制方法 |
CN105911867A (zh) * | 2016-06-16 | 2016-08-31 | 哈尔滨工程大学 | 基于nsga-ii算法的船舶推力分配方法 |
CN106055830A (zh) * | 2016-06-20 | 2016-10-26 | 吉林大学 | 基于动态规划的phev控制门限参数优化方法 |
CN108319138A (zh) * | 2018-01-29 | 2018-07-24 | 哈尔滨工程大学 | 一种欠驱动无人艇的滑模-反步双回路轨迹跟踪控制方法 |
CN108490770A (zh) * | 2018-02-28 | 2018-09-04 | 哈尔滨工程大学 | 一种基于混合算法的船舶动力定位系统推力分配方法 |
Non-Patent Citations (2)
Title |
---|
SOUMYADIP SENGUPTA,等: "Energy-efficient differentiated coverage of dynamic objects using an improved evolutionary multi-objective optimization algorithm with fuzzy-dominance", 《IEEE CEC》 * |
王奎民,等: "基于主尺度的深海采矿船多目标优化", 《船舶标准化工程师》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110794843A (zh) * | 2019-11-15 | 2020-02-14 | 山东交通学院 | 基于观测器的非线性船舶时滞动力定位船鲁棒镇定系统 |
CN110794843B (zh) * | 2019-11-15 | 2022-11-01 | 山东交通学院 | 基于观测器的非线性船舶时滞动力定位船鲁棒镇定系统 |
CN110909950A (zh) * | 2019-11-29 | 2020-03-24 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种采用非劣排序算法进行燃油消耗优化方法及装置 |
CN111930123A (zh) * | 2020-08-13 | 2020-11-13 | 上海船舶研究设计院(中国船舶工业集团公司第六0四研究院) | 多目标综合优化决策方法、装置及电子设备 |
CN114019799A (zh) * | 2021-11-06 | 2022-02-08 | 易如(山东)智能科技有限公司 | 基于多目标遗传算法优化拖拉机燃油经济及动力性的方法 |
CN114019799B (zh) * | 2021-11-06 | 2023-10-13 | 易如(山东)智能科技有限公司 | 基于多目标遗传算法优化拖拉机燃油经济及动力性的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109885061B (zh) | 2021-11-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109885061B (zh) | 一种基于改进nsga-ⅱ的动力定位多目标优化方法 | |
CN109901403A (zh) | 一种自主水下机器人神经网络s面控制方法 | |
CN110363343B (zh) | 一种混合自适应的水电站群智能优化调度方法及系统 | |
CN107203687B (zh) | 吸收塔脱硫过程多目标协同智能优化控制方法 | |
CN108460548B (zh) | 一种基于改进灰狼算法的多金属露天矿生产作业计划编制方法 | |
CN109931943B (zh) | 无人船舶全局路径规划方法及电子设备 | |
CN102411305A (zh) | 单旋翼直升机/涡轴发动机综合抗扰控制系统设计方法 | |
CN112132469B (zh) | 一种基于多种群协作粒子群算法的水库群调度方法和系统 | |
CN106934459A (zh) | 一种基于种群进化过程的自适应遗传算法 | |
CN108572550A (zh) | 一种基于机器学习的在线实时推力分配方法 | |
CN114442640B (zh) | 一种水面无人艇轨迹跟踪控制方法 | |
CN104252132A (zh) | 基于自适应遗传算法的行星际轨道控制优化方法 | |
CN115655279A (zh) | 一种基于改进鲸鱼算法的海上无人救援飞艇路径规划方法 | |
CN109901622B (zh) | 一种基于机理模型的自主水下机器人预测s面控制方法 | |
CN116300406A (zh) | 基于动态种群代价诱导粒子群优化的船舶模糊pid控制方法 | |
CN116108982A (zh) | 一种水库群多目标调度合作搜索方法及系统 | |
CN116542882B (zh) | 光伏发电平滑方法、系统和存储介质 | |
CN117452827A (zh) | 一种欠驱动无人艇轨迹跟踪控制方法 | |
Gong et al. | A path planning method based on improved particle swarm optimization algorithm | |
CN112464367B (zh) | 一种自主水下航行器的外形和结构两层设计优化方法 | |
Ding et al. | Thrust allocation of dynamic positioning based on improved differential evolution algorithm | |
CN114879481A (zh) | 一种抗水动力干扰的船舶动力定位鲁棒h∞控制方法 | |
CN109754108A (zh) | 基于变加速度系数混沌粒子群的机组经济负荷分配方法 | |
CN113238486A (zh) | 一种自适应的多机水电站调速器参数调控方法 | |
CN116929367A (zh) | 一种基于神经网络算法多步规划的推力分配方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |