1.一种基于改进粒子群算法的鲁棒泊位岸桥联合分配方法,其特征在于:具体步骤如下:S1:种群初始化;具体包括种群规模、迭代次数、粒子位置和速度的上下界、惯性权重和学习因子;
S2:计算最早靠泊时间tbj(ES):对于船舶j,计算tbj(ES)公式如下:其中sijk表示船舶j在泊位i第k个靠泊则为1,否则为0,taj表示船舶j的来港时间,tdj′表示在泊位i上船舶j′的离港时间;
S3:计算最晚靠泊时间tbj(LS):对于船舶j,若存在与其在同一泊位靠泊且邻接在其后靠泊的船舶为j′,则令t=tbj′(LS),若不存在这样的船舶则令t=tdj,计算tbj(LS)的公式如下:其中twj表示船舶j的工作时间;
S4:插入缓冲区;具体步骤如下:
S41:更新权重系数wj;权重系数代表船舶的服务优先权,对于船舶j,若存在使得所求tbj(ES)与tbj(LS)方案中有时间和空间有冲突的船舶j',则wj设为1,否则设为0;
S42:计算累计权重;计算对某艘船舶j在同一泊位其先前靠港的船舶的累积权重αj和在其后靠港的船舶的累积权重βj;
定义集合Fa(j)记录与船舶j在同一泊位靠泊且在其前靠泊的船舶,定义集合Fb(j)记录与船舶j在同一泊位靠泊且在其后靠泊的船舶,定义W作为所有船舶的总权重,则计算αj和βj的公式如下:S43:最终得到的插入缓冲区后的船舶靠泊时间tbj由tbj(ES)、tbj(LS)、αj、βj求得:S5:粒子自适应变异;设计如下自适应变异策略对传统PSO算法进行改进,令自适应变异概率P的更新公式如下:其中,tn为当前迭代次数,TN为总的迭代次数;对于粒子i,随机生成一个0到1间的随机数,如果该随机数大于此次迭代的变异概率,则对该粒子进行变异操作;
S6:计算样本的目标函数值;即船舶延迟离港时间,公式如下:
S7:更新粒子的个体最优值与全局最优值;对于每个粒子,将此次迭代所求的适应度值与其最优适应度值pbest做比较,如果小于pbest则更新,否则pbest不变;再将pbest与全局最优值gbest作比较,若小于gbest,则令pbest替换gbest;
S8:更新粒子的速度和位置,对第i个粒子的速度更新公式为:
vij(tn+1)=wvij(tn)+c1r1j(pij‑xij(tn))+c2r2j(gj‑xij(tn)),j∈V位置更新公式为:xij(tn+1)=xij(tn)+vij(tn+1),j∈V其中,vij(tn+1)表示第tn+1次迭代粒子i的第j维速度值,w为权重系数,c1和c2为学习因子,r1j与r2j为基于Logistic混沌序列的混沌变量,pij为粒子i的最优解的第j维位置值,gj为全局最优解的第j维位置值,xij(tn)表示第tn次迭代时粒子i的第j维位置值;对于更新速度和位置后的粒子,要对其进行边界处理,以让其编码满足约束;
S9:令tn=tn+1,判断是否满足达到最大迭代次数TN的终止条件,是则终止迭代输出最优解,否则返回第S2步继续迭代。
2.根据权利要求1所述的一种基于改进粒子群算法的鲁棒泊位岸桥联合分配方法,其特征在于:所述S5中的变异策略设计有两种,分别为交换和逆序,具体如下,交换:随机选择两列交换;逆序:从样本中随机选择两列,将包含所选两列及其之间的部分进行逆序排列。