[发明专利]一种基于多目标人工蜂群算法的生物多序列对比方法

权利要求书

1.一种基于多目标人工蜂群算法的生物多序列对比方法，其特征在于：所述方法包括以下步骤：

步骤1：设置多序列对比方法的初始参数，包括预先设定种群规模M、最大迭代次数G_max、蜜源开采限制次数Limit以及当前迭代次数iter＝0；

步骤2：对蜜源进行初始化及个体编码，并随机产生M个引领蜂；

设有N条序列参与比对，各条序列的长度依次为l₁，l₂，Λ，l_N，并随机产生l_max和1.2l_max之间的随机数L作为比对后序列的长度；其中l_max＝max(l₁,l₂,Λ,l_N)；

根据预设的序列比对的数学模型，将每条序列需要插入的空位数分别设为L-l_i，并随机生成N个1到L-l_i随机排列的一维矩阵α_i后采用二维矩阵进行个体编码，得到每一个蜜源的编码为γ＝[a₁,a₂,Λ,a_N]；其中，一维矩阵中存储的是各条序列所需插入空位的位置，在输入比对结果时将蜜源编码转化成对应的比对并以二维矩阵的形式存储；1≤i≤N；

确定随机产生的M个引领蜂及其对应的序列长度和序列；

步骤3：确定种群当前最优蜜源位置E_best；

对每一个随机产生的引领蜂均根据下述公式(4)进行计算，计算出目标函数f₁(S_i)和f₂(S_i)；其中，

minimize F(S′)=(f1(S′),f2(S′))Tf1(S′)=-Score(S′)f2(S′)=g(S′)S′∈Ω---(4)]]>

式(4)中，g(S')＝n₁×g_open+n₂×g_extend表示对齐序列组S'的空位罚分，n₁，n₂分别为空位开放的数量和空位扩展的数量；Ω表示多序列比对空间；

令F₁＝(f₁(S₁),f₁(S₂),Λ,f₁(S_L)以及F₂＝(f₂(S₁),f₂(S₂),Λ,f₂(S_L)，则计算出F＝(min(F₁),min(F₂))，并记录种群当前最优蜜源位置E_best；

步骤4：通过种群当前最优蜜源位置E_best对每一个随机产生的引领蜂进行单点交叉操作产生新的引领蜂V_i；

步骤5：按步骤3中的公式(4)计算每一个随机产生的引领蜂的适应度值F(E_i)及其对应新引领蜂的适应度值F(V_i)并进行比较；如果F(V_i)<F(E_i)，则用新引领蜂V_i的序列比对替代及其对应的随机产生的引领蜂E_i的序列比对；否则继续保持随机产生的引领蜂E_i的原始序列比对不变；

步骤6：对步骤5中所得的引领蜂的适应度值从小到大进行排序，再根据锦标赛选择策略式(7)，计算排序后的选择概率P_i(t)；

Pi(t)=ci(t)Σi=1Nci(t).---(7)]]>

其中c_i为每个个体的得分；

步骤7：跟随蜂根据步骤6中所得的选择概率P_i(t)选择蜜源，当一个引领蜂被选择作为蜜源时，跟随蜂在引领蜂邻域位置通过变异操作产生一个新的蜜源位置，引入一个新的变异操作，步骤如下：

步骤7.1：移动空位组：对各条序列随机选取一个空位组，对该空位组向左或向右移动一个残基位；

步骤7.2：合并空位组：先随机选择一条序列，接着在此序列中随机选择的一组空位，并将它和最近的空位组合并；

步骤7.3：划分空位组：先随机选择一条序列，接着在此序列中随机选择的一组空位，并将此空位组划分成两个大小相等的空位组；

步骤7.4：删除全为空位的列，产生新的蜜源位置；

步骤8：按步骤3中的公式(4)，计算每一个所选引领蜂的适应度值及其对应产生的新蜜源位置的适应度值F(V_i^/)并进行比较，如果F(V^/_i)＜F(E^/_i)，则用新蜜源位置的序列比对替代及其对应所选引领蜂的序列比对，且设置开采度标志向量trial(i)＝0；否则保持所选引领蜂的原始序列比对不变，且设置开采度标志向量trial(i)＝trial(i)+1；

步骤9：判断trial(i)是否大于Limit，若trial(i)＞Limit，则第i个引领蜂放弃当前蜜源而成为侦察蜂，侦察蜂按照如下局部搜索策略对种群进行搜索产生一个新的蜜源替代蜜源，

步骤9.1：随机选取序列中某个块，块大小为序列长度的[5％-25％]之间；

步骤9.2：删除被选中块的所有空位；

步骤9.3：对被选中块利用Kalign2比对算法重新产生比对结果；

步骤9.4：用Kalign2比对算法产生的新的比对结果替代原序列中的选中块，即得到一个新的蜜源；

步骤10：按步骤3中的公式(4)计算步骤9中所得各个蜜源的适应度值，并保存前M个最优蜜源作为下一次迭代的引领蜂；

步骤11：更新迭代次数iter＝iter+1，判断是否达到最大迭代次数G_max，如满足，则输出最优解，否则返回步骤4。

2.根据权利要求1所述的一种基于多目标人工蜂群算法的生物多序列对比方法，其特征在于：预先设定种群规模M取100、最大迭代次数G_max取5000以及蜜源开采限制次数Limit取25。