一种基于混沌免疫遗传机制的UUV智能体行为学习与演化模型的制作方法

本发明属于水下无人系统建模与仿真技术领域，具体涉及一种基于混沌免疫遗传机制的uuv智能体行为学习与演化模型。

背景技术：

军用uuv作为一种武器装备、水下环境和作战任务等作战要素相互耦合、相互影响的复杂水下无人作战系统，具有续航力大，隐蔽性好、风险性低，可回收等特点。该无人系统要完成搜索侦查、远程攻击、潜艇防御和组网探测等作战使命中担负重要角色，其就必须要具备强大的智能性，具体表现在uuv能够完成探测、识别、通信等功能，并且实现自主行为反应、规划和学习，确保水下作业任务的高质量和高效率完成。

如若uuv拥有高分辨率且完备的行为规则库，则其可以根据环境和任务特征信息进行行为反应、规划和学习。但是，目前uuv智能化程度和导航控制计算能力有限，形成的知识库和规则库规模较小，只能根据专家知识和仿真结果获得部分规则，距离uuv智能性建模所要求高分辨率且完备的行为规则库还有距离。如何形成完备的行为规则库，行为学习与演化方法成为uuv智能性建模的研究重点。

行为学习与演化方法主要分为两大类：机器学习和智能优化，其中机器学习通常采用神经网络、专家系统等人工智能领域的建模工具进行行为的学习，但是存在计算复杂度高、收敛慢的缺点，而且严重依赖于初始参数假设。智能优化则利用粒子群算法、遗传算法、模拟退后算法等进化算法对行为规则进行迭代寻优，具有实现容易、精度高、收敛快的特点，不但可以完成规则学习，而且能够满足行为规则的自适应演化，但是这类算法的全局搜索与局部搜索的兼容性较差。

免疫进化算法作为一种新兴的智能优化方法被提出用来实现并行分布式自适应系统的学习与演化，该算法具有抗原提取、自动识别、学习记忆等功能。目前在此基础上设计多种免疫算法，大多数方法是将遗传算法进行形式改变，或是结合粒子群优化算法等形成组合型免疫算法。

尽管这些免疫算法能够实现问题求解，但是在全局搜索空间中的搜索效率和准确性还有待提高。因此，希望开发出具有快速性和准确性显著的一种行为规则学习与演化方法，实现uuv水下无人系统自主行为的反应、规划与学习。

技术实现要素：

为了解决现有技术中存在的上述问题，本发明提供了一种基于混沌免疫遗传机制的uuv智能体行为学习与演化模型，该模型功能强大、知识完备，能够快速、准确进行规则的匹配与更新，实现uuv水下无人系统自主行为的反应、规划与学习。

本发明要解决的技术问题通过以下技术方案实现：

包括以下步骤：

步骤1：抗原识别：载入待求解问题及约束条件作为抗原ag；

步骤2：疫苗提取：利用专家和先验知识及依据抗原ag的特征作为疫苗信息hj，采用二进制编码方式形成基因组，构建规模大小为n1的疫苗群ah：

其中：t为矩阵转置符号；

步骤3：初始化抗体种群：k＝0时，记忆种群am的规模大小为n2，其中k为迭代次数，种群规模为n的初始化抗体种群ab＝[ab1,ab2,…,abi,…,abn]^t,i＝1,2,…,n，n2＜n，初始化抗体种群ab根据记忆种群am的具体情况由疫苗群ah和随机种群ar组成，或由疫苗群ah、记忆种群am和随机种群ar组成，或由记忆种群am和随机种群ar组成；其中随机种群ar采用混沌机制产生；具体混沌机制通过设定初始值ab0，ab0是启动logistic映射模型的初始值，采用一维logistic映射模型产生n/2规模的随机种群ar，具体混沌模型和随机种群ar分别为：

rl＝4×rl-1×(1-rl-1),l＝1,2,…,n/2(2)

ar＝[r1,r2,…,rl,…,rn/2]^t,l＝1,2,…,n/2(3)

其中，rl代表初始值ab0迭代的中间变量，直至形成n/2个随机种群ar；

步骤4：适应度计算：当k>1，计算初始化抗体种群ab中每一个抗体abi适应度fitnessk(abi),i＝1,2,…,n，选择适应度最高maxfitnessk(abi)的抗体进入抗体记忆库，记忆库规模更新为n2+1；

步骤5：疫苗接种：在初始化抗体种群ab中自适应选择规模为h(h<n)的抗体进行疫苗接种，接种规模随着迭代次数k自适应变化为：

其中，α和β为可调参数，e代表指数；

步骤6：结束判断：结束条件设置为是否满足最大迭代次数m，如果k满足k>＝m，停止，输出结果最优行为规则rule；如果k满足k<m，执行步骤7；

步骤7：抗体选择：计算抗体之间的亲和力以及抗体浓度，根据抗体期望繁殖率采用轮盘赌的方式进行抗体选择，每个抗体被选择的概率与其期望繁殖概率成正比，形成种群规模为n的新种群a_n；

抗体与抗体之间的亲和力反映了抗体之间的相似性，抗体浓度被提出衡量抗体之间的相似性，其主要计算相似抗体在抗体种群中所占的比例：

其中，concentrationk(abi)表示第k代抗体abi的浓度，si,j表示抗体abi和抗体abj的相似度，i(j)表示除抗体abi之外其余抗体的数量，ki,j表示抗体abi和抗体abj编码相同的位数，l表示抗体编码长度，τ表示相似度阈值，如果ki,j/l大于等于τ，则si,j＝1，如果ki,j/l小于τ，则si,j＝0；

基于抗体浓度，抗体abi的期望繁殖率计算为：

式中，fitnessk(abi)表示第k代抗体abi的适应度，λ为多样性评价系数；

根据抗体的期望繁殖率计算进行选择操作，采用轮盘赌的方式按照每次赌盘指针所落概率范围确定所选抗体，直到选择n个抗体种群a_n为止，具体为随机产生浮点数π∈(0,1]作为赌盘指针，如果pk(abi-1)＜π≤pk(abi)，则选择抗体abi，其中pk(abi)表示抗体abi期望繁殖率；

步骤8：抗体交叉：对种群a_n按照自适应交叉概率进行单点抗体交叉，形成种群规模为n的新种群b_n；

抗体依据自适应交叉概率进行单点交叉操作，形成新种群b_n，其中第k代交叉概率pk^c采用自适应机制进行计算：

式中，表示抗体种群a_n中的最大适应值，表示抗体种群的平均适应值，fitnessk表示要交叉的两个抗体abi和abj中较大的适应度，i,j＝1,2,…,n,k1和k2为预设常数；

步骤9：抗体变异：对种群b_n采用基于多项式与高斯混合的变异操作，形成种群规模为n的新种群c_n；

多项式变异具有良好的局部逃逸，其形式为和xi为抗体abi变异前后的反向译码数值，xmax和xmin分别为种群b_n反向译码的上限和下限，δk表示第k代的变异控制参数，其计算为：

其中，ak为满足[0,1]均匀分布的随机数，η表示变异分布指数；

高斯变异是指用符合均值为μ，方差为σ²的正态分布的一个随机数来替代原有抗体的基因值，具备较强的局部搜索能力，采用形式为xi^*＝xi+δk，和xi为抗体abi变异前后的反向译码数值，δk表示第k代的变异控制参数，其计算为：

δk＝0.1×n(0,1)(10)

其中，n(0,1)表示0均值1标准差的一维正态分布随机数；

为了结合多项式变异和高斯变异的优势，根据参数θk来实现两种变异操作的转换，其计算为：

其中，b为高斯变异的预设概率，c为预设参数，k为当前迭代次数，m为最大迭代次数，fitnessk(abi)为第k代抗体种群b_n中抗体abi(i＝1,2,…,n)的适应值，表示第k代抗体种群b_n中的最大适应值,表示第k代抗体种群b_n中的最小适应值。如果θk＜0.5，采用高斯变异，反之采用多项式变异，生成规模为n的新种群c_n；

步骤10：抗体早熟判断：计算种群c_n平均适应值进行早熟判断，如果早熟，进行正弦混沌变异的多样性维持，生成规模为n的新种群d_n，返回步骤7；否则，将种群c_n直接作为下一代抗体群，执行步骤11；

早熟判断基于种群平均适应值fitnessarg进行判断，具体为：

式中，为第k代抗体种群c_n的平均适应值，ω1和ω2为很小的预设正参数；

如果不满足公式(12)和公式(13)，则说明抗体没有发生早熟现象；如果同时满足公式(12)和公式(13)，则说明抗体种群发生早熟现象，采用基于正弦的混沌变异操作，维持种群多样性，具体形成种群规模为n的新种群d_n，计算如下：

其中，xi表示抗体abi的二进制反向解码后的量值，需要指出的是通过混沌变异生成的新抗体xi+1需要进行正向二进制编码；

步骤11：更新迭代次数k←k+1,返回步骤4。

进一步的，所述步骤3中的初始化抗体种群ab通过以下步骤来确定：

步骤3.1：如果抗体记忆种群am为空，即n2＝0，则初始化抗体种群ab由两部分组成：(1)n/2规模的种群根据疫苗群ah随机组合产生(2)n/2规模采用混沌机制由随机种群ar产生，即初始化抗体种群ab＝ah∪ar；

步骤3.2：如果抗体记忆种群am不为空，并且规模n2＜n/2，初始化抗体种群ab由三部分组成：(1)采用混沌机制随机产生的规模为n/2的随机种群ar(2)记忆库提供的规模为n2的抗体记忆种群am(3)由疫苗库产生的规模为n/2-n2的疫苗群ah，即初始化抗体种群ab＝ah∪ar∪am；

步骤3.3：如果抗体记忆种群am不为空，并且规模n2＞n/2，初始化抗体种群ab由两部分组成：(1)采用混沌机制随机产生的规模为n/2的随机种群ar(2)记忆库提供的规模为n/2的抗体记忆种群am，即初始化抗体种群ab＝ar∪am。

进一步的，所述步骤5中，比较接种前父代抗体abi与抗原的适应度fitnessk(abi)和接种后子代抗体abi^s与抗原的适应度

步骤5.1：如果适应度提高，则在初始化抗体种群ab中用接种过的子代抗体替代接种前的父代抗体，即则

步骤5.2：如果适应度降低，表示抗体种群在退化，则遗弃子代抗体，保留父代抗体，即则abi维持不变。

与现有技术相比，本发明的有益效果：

本发明模型结合基本遗传算法全局搜索能力和免疫与混沌机制局部搜索能力的优势，通过不断调整和优化问题解的搜索空间，从而促进行为规则的快速学习与演化。

具体的，通过抗体适应度计算，不断的对抗体记忆库规模进行更新；依次设计基于轮盘赌的选择算子、基于自适应调整的交叉算子和基于高斯与多项式混合的变异算子实现抗体种群的多样性，并进行早熟抑制，进而实现抗体种群的更新与迭代。

附图说明

图1是基于混沌免疫遗传机制的uuv智能体行为学习与演化模型示意图。

图2是uuv规避障碍示意图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚明了，下面通过附图中示出的具体实施例来描述本发明。但是应该理解，这些描述只是示例性的，而并非限制本发明的范围。

实施例1：

为了能够更好的理解本发明的技术方案，作了以下更具体的说明：

如图1，一种基于混沌免疫遗传机制的uuv智能体行为学习与演化模型，包括以下步骤：

步骤1：抗原识别：载入待求解问题及约束条件作为抗原ag；

以uuv在10000m×10000m范围的水下空间实现自主规避水下障碍为求解问题，将作战环境、目标状态及自身状态等作为约束条件，实现抗原的显示表达为ag。

步骤2：疫苗提取：利用专家和先验知识并依据抗原ag的特征作为疫苗信息hj，其采用二进制编码方式形成基因组，构建规模大小为n1的疫苗群ah：

其中：t为矩阵转置符号；

利用专家和先验知识依据问题及约束条件的特征信息如作战环境的洋流流速velocity、目标状态的位置信息position、自身状态的生命状态state和行为动作rule作为疫苗信息hj，采用四元组表示：

hj＝<velocity,position,state,rule>(1.1)

步骤3：初始化抗体种群：k＝0时，记忆种群am的规模为n2，其中k为迭代次数，算法初始化从k＝0开始；种群规模为n的初始化抗体群ab＝[ab1,ab2,…,abi,…,abn]^t,i＝1,2,…,n，n2＜n，初始化抗体种群ab根据记忆种群am的具体情况由疫苗群ah和随机种群ar组成，或由疫苗群ah、记忆种群am和随机种群ar组成，或由记忆种群am和随机种群ar组成；其中随机种群ar采用混沌机制产生。具体混沌机制通过设定初始值ab0，采用一维logistic映射模型产生n/2规模的随机种群ar，具体混沌模型和随机种群ar分别为：

rl＝4×rl-1×(1-rl-1),l＝1,2,…,n/2(2)

ar＝[r1,r2,…,rl,…,rn/2]^t,l＝1,2,…,n/2(3)

其中，rl代表初始值ab0迭代的中间变量，直至形成n/2个随机种群ar；

步骤3.1：如果抗体记忆库为am为空，即n2＝0，则抗体群由两部分组成：(1)n/2规模的种群根据疫苗群ah随机组合产生，(2)n/2规模采用混沌机制由随机种群ar产生，即抗体群ab＝ah∪ar；

步骤3.2：如果抗体记忆库am不为空，并且规模n2＜n/2，抗体群由三部分组成：(1)采用混沌机制随机产生的规模为n/2的随机种群ar，(2)记忆库提供的规模为n2的记忆种群am，(3)由疫苗库产生的规模为n/2-n2的疫苗群ah，即抗体群ab＝ah∪ar∪am；

步骤3.3：如果抗体记忆库am不为空，并且规模n2＞n/2，抗体群由两部分组成：(1)采用混沌机制随机产生的规模为n/2的随机种群ar，(2)记忆库提供的规模为n/2的记忆种群am，即抗体群ab＝ar∪am；

步骤4：适应度计算：当k>1，计算抗体种群ab中每一个抗体abi适应度fitnessk(abi),i＝1,2,…,n，选择适应度最高maxfitnessk(abi)的抗体进入抗体记忆库，记忆库规模更新为n2+1；

步骤5：疫苗接种：在抗体群ab中自适应选择规模为h(h<n)的抗体进行疫苗接种，接种规模随着迭代次数k自适应变化为：

其中，α和β为可调参数，e代表指数；

进一步，比较接种前父代抗体abi与抗原的适应度fitnessk(abi)和接种后子代抗体与抗原的适应度

步骤5.1：如果适应度提高，则在抗体群中用接种过的子代抗体替代接种前的父代抗体，即则

步骤5.2：如果适应度降低，表示抗体种群在退化，则遗弃子代抗体，保留父代抗体，即则abi维持不变；

步骤6：结束判断：结束条件设置为是否满足最大迭代次数m，如果k满足k>＝m，停止，输出结果最优行为规则rule；如果k满足k<m，执行步骤7；

步骤7：抗体选择：计算抗体之间的亲和力以及抗体浓度，根据抗体期望繁殖率采用轮盘赌的方式进行抗体选择，每个抗体被选择的概率与其期望繁殖概率成正比，形成种群规模为n的新种群a_n；

抗体与抗体之间的亲和力反映了抗体之间的相似性，抗体浓度被提出衡量抗体之间的相似性，其主要计算相似抗体在抗体种群中所占的比例：

其中，concentrationk(abi)表示第k代抗体abi的浓度，si,j表示抗体abi和抗体abj的相似度，i(j)表示除抗体abi之外其余抗体的数量，ki,j表示抗体abi和抗体abj编码相同的位数，l表示抗体编码长度，τ表示相似度阈值，如果ki,j/l大于等于τ，则si,j＝1，如果ki,j/l小于τ，则si,j＝0。

基于抗体浓度，抗体abi的期望繁殖率计算为：

式中，fitnessk(abi)表示第k代抗体abi的适应度，λ为多样性评价系数。

根据抗体的期望繁殖率计算进行选择操作，采用轮盘赌的方式按照每次赌盘指针所落概率范围确定所选抗体，直到选择n个抗体种群a_n为止，具体为随机产生浮点数π∈(0,1]作为赌盘指针，如果pk(abi-1)＜π≤pk(abi)，则选择抗体abi，其中pk(abi)表示抗体abi期望繁殖率。

步骤8：抗体交叉：对种群a_n按照自适应交叉概率进行单点抗体交叉，形成种群规模为n的新种群b_n；

抗体依据自适应交叉概率进行单点交叉操作，形成新种群b_n，其中第k代交叉概率采用自适应机制进行计算：

式中，表示抗体种群a_n中的最大适应值，表示抗体种群的平均适应值，fitnessk表示要交叉的两个抗体abi和abj中较大的适应度，i,j＝1,2,…,n,k1和k2为预设常数。

步骤9：抗体变异：对种群b_n采用基于多项式与高斯混合的变异操作，形成种群规模为n的新种群c_n；

多项式变异具有良好的局部逃逸，形式为和xi为抗体abi变异前后的反向译码数值，xmax和xmin分别为抗体群b_n反向译码的上限和下限，δk表示第k代的变异控制参数，其计算为：

其中，ak为满足[0,1]均匀分布的随机数，η表示变异分布指数。

高斯变异是指用符合均值为μ，方差为σ²的正态分布的一个随机数来替代原有抗体的基因值，具备较强局部搜索能力，采用形式为xi^*＝xi+δk，xi^*和xi为抗体abi变异前后的反向译码数值，δk表示第k代的变异控制参数，其计算为：

δk＝0.1×n(0,1)(10)

其中，n(0,1)表示0均值1标准差的一维正态分布随机数。

为了结合多项式变异和高斯变异的优势，根据参数θk来实现两种变异操作的转换，其计算为：

其中，b为高斯变异的预设概率，c为预设参数，k为当前迭代次数，m为最大迭代次数，fitnessk(abi)为第k代抗体种群b_n中抗体abi(i＝1,2,…,n)的适应值，表示第k代抗体种群b_n中的最大适应值,表示第k代抗体种群b_n中的最小适应值。如果θk＜0.5，采用高斯变异，反之采用多项式变异，生成规模为n的新种群c_n。

步骤10：抗体早熟判断：计算种群c_n平均适应值进行早熟判断，如果早熟，进行正弦混沌变异的多样性维持，生成规模为n的新种群d_n，返回步骤7；否则，将种群c_n直接作为下一代抗体群，执行步骤11；

早熟判断基于种群平均适应值fitnessarg进行判断，具体为：

式中，为第k代抗体种群c_n的平均适应值，ω1和ω2为很小的预设正参数。

如果同时不满足公式(12)和公式(13)，则说明抗体没有发生早熟现象；如果同时满足公式(12)和公式(13)，则说明抗体种群发生早熟现象，采用基于正弦的混沌变异操作，维持种群多样性，具体形成种群规模为n的新种群d_n，计算如下：

其中，xi表示抗体abi的二进制反向解码后的量值，需要指出的是通过混沌变异生成的新抗体xi+1需要进行正向二进制编码。

步骤11：更新迭代次数k←k+1,返回步骤4。

根据本发明提出的基于混沌免疫遗传机制的uuv智能体行为学习与演化模型，形成最优行为规则库，uuv依据最优行为规则实现了规避障碍如图2所示：图2表示以uuv水下规避障碍为实例进行智能体行为学习与演化，在10000m×10000m范围的水下空间，uuv抽象为二维运动质点，初始航向，速度2.5m/s，探测范围1000m，从起始位置(0,0)自主规避水下障碍运动到目标位置(10000m,10000m)，其中水下静态障碍区域为：障碍#1：半径(800m),位置(4200m,4500m)，障碍#2：半径(800m),位置(8000m,7700m)。

总的来说，本发明提出的基于混沌免疫遗传机制的uuv智能体行为学习与演化模型可以实现uuv行为的自主学习和演化等智能特性。

以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护范围。