1.一种带监督机制的推力分配方法，其特征在于：包括如下步骤：

步骤1、初始化船舶各类参数：结构参数、状态参数、物理条件约束参数、权值参数、扇区参数；

步骤2、判别是否执行最优扇区搜寻，是则运用PSO算法进行最优扇区求解，否则保留前一时刻最优扇区解；

执行最优扇区搜寻条件为：首次执行推力分配或当前需求推力和力矩与上一时刻的推力和力矩不相等；

PSO算法进行最优扇区求解实施方法如下：

步骤2.1、初始化参数，其中粒子选取推力分配方程等式约束项的齐次线性方程组BeTe＝0的基础解系数；

Te为船体坐标系下X方向和Y方向两个相互正交的力，(lxn,lyn)为第n个推进器在船体坐标系中的位置；

步骤2.2、以能耗项和偏差项为代价函数计算粒子适应度值，并计算全局最优值，该步骤中推进器角度在禁区范围内处理方法为：计算该角度与该禁区极值的距离，靠近禁区极大值则取该禁区极大值为推进器角度值，否则取禁区极小值为推进器角度值；

步骤2.3、运用带权重系数的PSO算法进行种群进化，粒子更新；

步骤2.4、判别迭代次数是否达设定值，是执行步骤2.5，否执行步骤2.3；

步骤2.5、对所求最优推力分配解的推进器角度进行扇区判别，具体实施方法如下：首先判别角度值是否处于前一时刻选择扇区内，是则判定最优扇区取前一时刻选择扇区，否则判别角度值是否处于前一时刻之前选择扇区内，是则判定最优扇区取前一时刻之前选择扇区，否则按照扇区序数从小到大逐一判别进行最优扇区取值；

步骤3、判别推进器所处扇区是否与前一时刻选择扇区相一致，是则执行步骤4，否则执行步骤5；

步骤4、运用非扇区切换的SQP算法进行推力分配求解，然后执行步骤6，非扇区切换的SQP算法具体实施方法如下：构建二次型推力分配数学模型为：

式中H是包含能耗项、偏差项、机械磨损项和奇异项的二次项权系数矩阵，f是对应的一次项权系数向量，X是包含推力变化量、方位角变化量和松弛变量的优化变量，Aeq和beq分别为等式约束矩阵和向量，Lb和Ub为包含推力变化量、方位角变化量和松弛变量的优化变量边界；

beq＝τ‑B(αP)TP

ΔSTimin＝max(ΔTimin,Timin‑TiP)，ΔSTimax＝min(ΔTimax,Timax‑TiP)ΔSαimax＝min(Δαimax,αimax(Snow)‑αiP)，ΔSαimin＝max(Δαimin,αimin(Snow)‑αiP)其中B(αP)是前一时刻推进器配置控制矩阵，αP和TP分别为前一时刻方位角和推力，E3为三阶单位矩阵，τ为控制系统输出的推力和力矩指令，ΔSTimin和ΔTimin分别是第i个推进器当前和自身物理条件允许的推力最小变化量，ΔSTimax和ΔTimax分别是第i个推进器当前和自身物理条件允许的推力最大变化量，ΔSαimin和ΔSαimax分别是第i个推进器允许的方位角最小和最大变化量，Δαimax和Δαimin分别是第i个推进器自身物理条件允许的方位角最大和最小变化量，αimax(Snow)和αimin(Snow)为第i个推进器前一时刻选择扇区的边界，αiP和TiP分别为第i个推进器当前角度和推力，Snow为前一时刻选择扇区值；

根据构建的二次型推力分配数学模型利用SQP算法进行推力求解；

步骤5、运用扇区切换的SQP算法进行推力分配求解，具体实施方法如下：

构建二次型推力分配数学模型为：

式中HS是包含能耗项和偏差项的二次项权系数矩阵，fS是对应的一次项权系数向量，XS是包含推力变化量和松弛变量的优化变量，ASeq和bSeq分别为等式约束矩阵和向量，LSb和USb为包含推力变化量和松弛变量的优化变量边界，其中ASeq和bSeq计算方法如下：ASeq＝[B(αP+Δα)| E3]，bSeq＝τ‑B(αP+Δα)TP其中B(αP+Δα)为此刻的推进器配置控制矩阵，Δα为方位角变化量；

根据构建的二次型推力分配数学模型利用SQP算法进行推力求解；

步骤6、运用监督与切换机制进行扇区选择，该机制实施思路为：前一时刻选择扇区每次更新只有两种可能，一是保持不变，二是被最优扇区取代，而要被最优扇区取代只要满足如下三个条件中一个：(1)最优扇区合力和力矩偏差比当前扇区小于一定值，且保持一段时间；(2)最优扇区功率比当前扇区小于一定值，最优扇区合力和力矩偏差比当前扇区不大于一定值，且保持一段时间；(3)最优扇区与前一时刻之前选择扇区相等且切换时间不超过一定值；

步骤7、输出推力分配结果，更新船舶状态参数，跳至步骤2进行下一个周期推力分配计算。

2.根据权利要求1所述的一种带监督机制的推力分配方法，其特征在于：所述步骤2中，推进器角度在禁区范围内处理具体计算方法如下：设推进器在禁区内的角度为α，禁区范围为[αstrmin,αstrmax]，则α处理结果为：

3.根据权利要求1所述的一种带监督机制的推力分配方法，其特征在于：所述步骤5中，SQP算法进行推力分配求解中方位角变化量具体计算方法为：

4.根据权利要求1所述的一种带监督机制的推力分配方法，其特征在于：所述步骤6中，运用监督与切换机制进行扇区选择具体实施计算方法如下：假设前一时刻之前选择的扇区为Sago，当前最优扇区为Sbest，当前选择的扇区为Schoose，则：t是指Sbest和Snow连续不等累加的时间，tlim是指停留时间切换的最小切换时间，t0是为了防止瞬态干扰导致的切换所设置的最大返回切换时间，Pbest是全局最优消耗功率值，PSQP是扇区切换或非扇区切换的SQP算法求解所得消耗功率值，eSQP是扇区切换或非扇区切换的SQP算法求解结果偏差的平方和，ebest是全局最优求解结果偏差的平方和，eS为偏差设定值，P1和P2为设定的百分数。