一种模型参数的计算方法、装置、设备和存储介质与流程

1.本发明涉及模型参数求解技术领域，具体而言，涉及一种模型参数的计算 方法、装置、设备和存储介质。


背景技术：

2.阿奎拉鹰模型，通过模拟阿奎拉鹰的狩猎方式，设计有探索阶段和开发阶 段。在探索阶段模拟阿奎拉鹰的高空飞行搜索并垂直俯身攻击和环绕猎物飞行 并短距离滑翔攻击。在开发阶段模拟阿奎拉鹰的低空飞行并缓降攻击和地面近 距离攻击。
3.阿奎拉鹰模型，能够用于求解各种计算模型的解。但是，其容易陷入局部 最优，导致求解精度降低。在求解计算模型时，可能无法得到最优解。
4.压力容器是指盛装气体或者液体，承载一定压力的密闭设备。在各个行业 的发展过程中，必须依靠各种设备，其中压力容器在我国许多地区是应用非常 广泛且非常重要的设备。
5.压力容器的设计是一个典型的非线性混合整数规划问题。压力容器设计问 题共有4个参数需要优化，分别是容器中间部分的长度、厚度、内部半径与半 球形厚度。在先技术中的，压力容器的参数的计算方法普遍存在大量的局部最 优点，容易陷入局部最优，导致收敛精度低、收敛速度慢等问题。有鉴于此， 申请人在研究了现有的技术后特提出本技术。


技术实现要素：

6.本发明提供了一种模型参数的计算方法、装置、设备和存储介质，以改善 相关技术中的设计方法收敛精度低、收敛速度慢的问题。
7.第一方面、
8.本发明实施例提供了一种模型参数的计算方法，其至少包含步骤s1和s2。
9.s1、获取待求解的目标对象的计算模型。
10.s2、通过改进的阿奎拉鹰模型对所述计算模型进行求解，以获得所述目标 对象的设计参数。其中，所述改进的阿奎拉鹰优化模型通过logistic混沌映射计 算初始种群，且在开发阶段采用精英狩猎策略进行局部开发，以避免过早收敛。
11.在一个可选地实施例中，步骤s2至少包括步骤s21至步骤s26。
12.s21、根据所述计算模型，随机生成初始粒子种群，并初始化迭代次数和最 大迭代次数。
13.s22、根据所述初始粒子种群，通过logistic混沌映射模型，计算得到多样 化粒子种群。
14.s23、计算所述多样化粒子种群中每个粒子的适应度，并判断迭代次数是否 不大于预设迭代次数。
15.s24、当迭代次数不大于预设迭代次数时，采用大范围搜索攻击策略更新粒 子种
[0037][0038][0039][0040]
其中，x2(t 1)表示粒子的下一个更新位置，x
best
(t)代表目前已取得的最佳位 置，lf(d)代表莱维飞行函数，d代表种群维度，x
r
(t)代表任意粒子的位置，r2为 均匀分布在0到1之间的随机数，s＝0.01，u与v为0到1之间的随机数，β＝1.5。
[0041]
在一个可选地实施例中，步骤s25至少包括步骤s251至步骤s253。
[0042]
s251、获取第二随机数，并判断随机数是否小于第二预设值。
[0043]
s252、当判断到所述第二随机数小于所述第二预设值时，采用低空飞行并 缓降攻击的策略更新粒子位置。具体地，低空飞行并缓降攻击的位置更新模型 为：
[0044]
x3(t 1)＝(x
best
(t)
‑
x
m
(t))
×
α
‑
r4 ((ub
‑
lb)
×
r5 lb)
×
δ
[0045][0046]
其中，x3(t 1)表示粒子的下一个更新位置，x
best
(t)代表目前已取得的最佳位 置，x
m
(t)是当前所有粒子的平均位置，α与δ为适应参数，其值为0.1，ub是上 边界，lb为下边界，r4和r5为0到1之间的随机数。x
i
(t)为当前第i个粒子位置， n粒子的数量。
[0047]
s253、当判断到所述第二随机数不小于所述第二预设值时，采用地面近距 离攻击猎物的策略更新粒子位置。具体地，地面近距离攻击猎物的位置更新模 型为：
[0048][0049]
其中，x4(t 1)表示粒子的下一个更新位置，qf(t)是质量函数，用于平衡其 搜索策略，x
best
(t)代表目前已取得的最佳位置，g1表示阿奎拉鹰追踪猎物的移动 参数，为[
‑
1,1]之间的随机数，x(t)为当前的粒子位置，r6、r7和r8为0到1之间 的随机数，g2在空中追逐猎物的飞行斜率，从2到0线性递减，lf(d)代表莱维 飞行函数，t是当前迭代次数，t为最大迭代次数。
[0050]
在一个可选地实施例中，通过精英狩猎策略计算得到精英位置的精英位置 计算模型为：
[0051]
x
e
(t 1)＝x
e
‑
best
cd
×
(x(t)
‑
x
e
‑
normal
(t)) σ
×
(x
rand1
‑
x
rand2
)
[0052][0053]
cd＝z0 0.1
×
tan(π
×
(r
‑
0.5))
[0054]
其中，x
e
(t 1)是基于精英狩猎策略产生的新位置，x
e
‑
best
为从5个精英个体 中随机选出的个体，cd为符合柯西分布的随机数，x(t)为当前的粒子位置， x
e
‑
normal
(t)为从仓库中随机选出的15个位置，x
rand1
与x
rand2
为从剩下的15个个体 中随机选出的，exponent＝2，σ
initial
＝1σ
final
＝0，z0＝2，t为当前迭代次数，t为 最大迭代次数，r为分布在[0,1]之间的随机数。
[0055]
在一个可选地实施例中，所述预设迭代次数为最大迭代次数的三分之二。
[0056]
在一个可选地实施例中，所述目标对象为压力容器。
[0057]
第二方面、
[0058]
本发明实施例提供了一种模型参数的计算装置，其包含：
[0059]
模型获取模块，用于获取待求解的目标对象的计算模型。
[0060]
模型求解模块，用于通过改进的阿奎拉鹰模型对所述计算模型进行求解， 以获得所述目标对象的设计参数。其中，所述改进的阿奎拉鹰优化模型通过 logistic混沌映射计算初始种群，且在开发阶段采用精英狩猎策略进行局部开发， 以避免过早收敛。
[0061]
在一个可选地实施例中，所述模型求解模块，具体包括：
[0062]
初始化单元，用于根据所述计算模型，随机生成初始粒子种群，并初始化 迭代次数和最大迭代次数。
[0063]
混沌映射单元，用于根据所述初始粒子种群，通过logistic混沌映射模型， 计算得到多样化粒子种群。
[0064]
适应度单元，用于计算所述多样化粒子种群中每个粒子的适应度，并判断 迭代次数是否不大于预设迭代次数。
[0065]
大范围单元，用于当迭代次数不大于预设迭代次数时，采用大范围搜索攻 击策略更新粒子种群的位置，并将迭代次数加一，和判断是否达到最大迭代次 数。
[0066]
局部单元，用于当迭代次数大于预设迭代次数时，采用局部攻击策略更新 粒子种群的位置，以及通过精英狩猎策略计算得到精英位置，然后通过贪婪策 略获取最佳位置，并将迭代次数加一，和判断是否达到最大迭代次数。
[0067]
输出单元，用于当判断到未达到最大迭代次数时，计算更新后粒子种群中 每个粒子的适应度，并判断迭代次数是否不大于预设迭代次数。否则，输出最 大迭代次数时的最佳位置及其对应的适应度，以获得所述设计参数。
[0068]
第三方面、
[0069]
本发明实施例提供了一种模型参数的计算设备，其包括处理器、存储器， 以及存储在所述存储器内的计算机程序。所述计算机程序能够被所述处理器执 行，以实现如第一方面所说的模型参数的计算方法。
[0070]
第四方面、
[0071]
本发明实施例提供了一种计算机可读存储介质，所述计算机可读存储介质 包括存储的计算机程序，其中，在所述计算机程序运行时控制所述计算机可读 存储介质所在设备执行如第一方面所说的模型参数的计算方法。
[0072]
通过采用上述技术方案，本发明可以取得以下技术效果：
[0073]
本发明实施例通过在种群初始化阶段引入logistic map混沌映射，提高初 始化种群分布的多样性，在局部开发阶段引入精英狩猎策略，建立小范围的种 群库，利用精英个体与普通个体为导向进行局部开发，避免算法早熟收敛，能 够取得更好的计算效果。
[0074]
为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例， 并配合所附附图，作详细说明如下。
附图说明
[0075]
为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使 用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例， 因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创 造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
[0076]
图1是本发明第一实施例提供的一种模型参数的计算方法的流程示意图。
[0077]
图2是通过改进的阿奎拉鹰模型对计算模型进行求解的逻辑框图。
[0078]
图3是logistic混沌的点集图形。
[0079]
图4是压力容器结构模型。
[0080]
图5是本发明第二实施例提供的一种模型参数的计算装置的结构示意图。
具体实施方式
[0081]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清 楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是 全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造 性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0082]
为了更好的理解本发明的技术方案，下面结合附图对本发明实施例进行详 细描述。
[0083]
应当明确，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实 施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前 提下所获得的所有其它实施例，都属于本发明保护的范围。
[0084]
在本发明实施例中使用的术语是仅仅出于描述特定实施例的目的，而非旨 在限制本发明。在本发明实施例和所附权利要求书中所使用的单数形式的“一 种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。
[0085]
应当理解，本文中使用的术语“和/或”仅仅是一种描述关联对象的关联关 系，表示可以存在三种关系，例如，a和/或b，可以表示：单独存在a，同时存 在a和b，单独存在b这三种情况。另外，本文中字符“/”，一般表示前后关联 对象是一种“或”的关系。
[0086]
取决于语境，如在此所使用的词语“如果”可以被解释成为“在
……
时
”ꢀ
或“当
……
时”或“响应于确定”或“响应于检测”。类似地，取决于语境，短 语“如果确定”或“如果检测(陈述的条件或事件)”可以被解释成为“当确定 时”或“响应于确定”或“当检测(陈述的条件或事件)时”或“响应于检测 (陈述的条件或事件)”。
[0087]
实施例中提及的“第一\第二”仅仅是是区别类似的对象，不代表针对对象 的特定
排序，可以理解地，“第一\第二”在允许的情况下可以互换特定的顺序 或先后次序。应该理解“第一\第二”区分的对象在适当情况下可以互换，以使 这里描述的实施例能够以除了在这里图示或描述的那些内容以外的顺序实施。
[0088]
下面结合附图与具体实施方式对本发明作进一步详细描述：
[0089]
实施例一：
[0090]
请参阅图1至图4，本发明第一实施例提供的一种模型参数的计算方法，其 可由模型参数的计算设备来执行。特别地，由计算设备中的一个或多个处理器 来执行，以至少实现步骤s1和s2。
[0091]
s1、获取待求解的目标对象的计算模型。在本实施例中，目标对象为压力 容器。需要说明的是，压力容器的设计是一个典型的非线性混合整数规划问题。 压力容器设计问题普遍存在大量的局部最优点。因此，在设计的时候，往往无 法获取到最优解。
[0092]
如图4所示，压力容器设计问题共有4个参数需要优化，分别是容器中间 部分的长度l、厚度ts、内部半径r与半球形厚度th。具体地，压力容器的计 算模型为：
[0093][0094][0095][0096]
其中，ts为x1、th为x2、r为x3、l为x4。
[0097]
可以理解的是，所述计算设备，可以是本地的计算机、膝上型便携式计算 机，本地的服务器，或者云端的服务器，本发明对此不做具体限定。在其它实 施例中，目标对象还可以是其它研究对象，本发明对此亦不做具体限定。
[0098]
s2、通过改进的阿奎拉鹰模型对计算模型进行求解，以获得目标对象的设 计参数。其中，改进的阿奎拉鹰优化模型通过logistic混沌映射计算初始种群， 且在开发阶段采用精英狩猎策略进行局部开发，以避免过早收敛。
[0099]
阿奎拉鹰优化模型(aquila optimizer，ao)由于算法新颖、优化性能较 高而得到了大量关注，目前已被用于分类问题、神经网络训练、特征选择等问 题，并取得了不错的效果。尽管ao在解决某些问题上具有一定的优势，但算法 后期寻优时仅根据当前最优位置作为导向，容易陷入局部最优，导致算法存在 收敛精度低、收敛速度慢等问题。
[0100]
为解决上述问题，本实施例提供了一种改进的阿奎拉鹰优化模型(improved
aquila optimizer，iao)，改进的阿奎拉鹰优化模型在种群初始化阶段引入 logistic map混沌映射，提高初始化种群分布的多样性，在局部开发阶段引入 精英狩猎策略，建立小范围的种群库，利用精英个体与普通个体为导向进行局 部开发，避免算法早熟收敛。改进的阿奎拉鹰优化模型明显好于原始阿奎拉鹰 优化模型，尤其是在压力容器设计问题上达到了较优水平。
[0101]
具体地，如图2所示，在上述实施例的基础上，本发明一可选实施例中， 步骤s2至少包括步骤s21至步骤s26。
[0102]
s21、根据计算模型，随机生成初始粒子种群，并初始化迭代次数和最大迭 代次数。具体地，步骤s21至少包括步骤s211和步骤s212。
[0103]
s211、根据计算模型，获取搜索空间的上界和下界。
[0104]
s212、根据随机函数，在搜索空间内生成初始粒子种群。具体地，初始种 群生成模型为：
[0105]
x＝rand
×
(ub
‑
lb) lb
[0106]
其中，x为粒子位置，rand为均匀分布在0到1之间的随机数，ub为搜索 空间的上边界，lb为搜索空间的下边界。
[0107]
s22、根据初始粒子种群，通过logistic混沌映射模型，计算得到多样化粒 子种群。具体地，logistic混沌映射模型为：
[0108]
x(t 1)＝a
×
x(t)
×
(1
‑
x(t))
………………
(式1)
[0109]
其中，x(t 1)为更新后的粒子位置，a＝4，x(t)为当前的粒子位置。
[0110]
需要说明的是，如图3所示，logistic混沌是非线性系统的一个特征，在 数学上被定义为由简单确定性函数生成的随机数。在本实施例中，将混沌思想 融入到改进的阿奎拉鹰优化模型，其核心是利用混沌变量的随机性、遍历性和 不规律性，替换服从标准概率分布的随机变量进行优化搜索，从而搜寻整个解 空间，产生分布均匀的个体，提高初始化种群分布的多样性和初始种群质量。
[0111]
s23、计算多样化粒子种群中每个粒子的适应度，并判断迭代次数是否不大 于预设迭代次数。在本实施例中，预设迭代次数为最大迭代次数的三分之二。 在其它实施例中，预设迭代次数也可以是其它比例，或者其它预定次数，本发 明对此不做具体限定。具体地，粒子的适应度的适应度计算模型为：
[0112][0113]
其中，f(x)是惩罚函数值，f(x)是原函数值，m是惩罚因子,g
i
(x)为约束条 件，n为粒子数量。
[0114]
压力容器设计问题是一个约束优化问题，我们要将其转化为无约束优化问 题，通过惩罚函数法实现。对于种群中没有违反约束条件的个体x，由于此时所 有约束条件都满足，即g1≤0，g2≤0，g3≤0，g4≤0，则0与g1(x)二者中最 大者使0，从而可得max(0,g1(x))＝0,其中的max是求最大值函数。同理可得， max(0,g2(x))＝0,
…
,max(0,gn(x))＝0.从而满足适应度计算模型。
[0115]
对于没有违反约束条件的个体x，加入惩罚项以后的惩罚函数值f(x)与原 来的函
数值f(x)相等，则惩罚项不起任何作用；对于违反约 束条件的个体x，由于违反了约束条件，其惩罚项>0，会导致该个体的惩罚函数 值f(x)大于原来函数值f(x)，使得该个体在进化中由于惩罚函数值较差而被逐 步淘汰。
[0116]
所述的惩罚函数法是在不合格的解中加上一个惩罚项使得某些只有某一维 上违反约束的不合格解不会被放弃能够继续迭代，从而避免了种群中合格的解 越来越少的问题，同时得到质量较差的解进而使得未违反约束的解成为好解， 逐渐淘汰这些质量较差的解。
[0117]
s24、当迭代次数不大于预设迭代次数时，采用大范围搜索攻击策略更新粒 子种群的位置，并将迭代次数加一，和判断是否达到最大迭代次数。需要说明 的是，阿奎拉鹰模型具有探索阶段和开发阶段，当迭代次数不大于预设迭代次 数时处于探索阶段，当迭代次数大于预设次数时处于开发阶段。
[0118]
探索阶段阿奎拉鹰在随机范围内通过两种不同的狩猎方式实现搜索空间内 的大范围搜索猎物并攻击，两种狩猎方式的介绍如下：
[0119]
1、高空飞行搜索并垂直俯身攻击：
[0120]
阿奎拉鹰在高空飞行，通过眼睛观察位于其下方飞行的猎物，确定猎物范 围后进行垂直俯身攻击，其计算公式如下：
[0121][0122][0123]
其中，x1(t 1)表示粒子的下一个更新位置，x
best
(t)代表目前已取得的最佳位 置，x
m
(t)是当前所有粒子的平均位置，t为当前迭代次数，t为最大迭代次数， n为种群数量，r1为0到1之间的随机数，x
i
(t)为当前第i个粒子位置。
[0124]
2、环绕猎物飞行并短距离滑翔攻击：
[0125]
阿奎拉鹰通常采用此狩猎方式，逐渐下降高度并环绕猎物，采用短距离滑 翔攻击猎物，其位置更新公式如下：
[0126]
x2(t 1)＝x
best
(t)
×
lf(d) x
r
(t) (y
‑
x)
×
r2[0127][0128]
[0129][0130]
其中，x2(t 1)表示粒子的下一个更新位置，x
best
(t)代表目前已取得的最佳位 置，lf(d)代表莱维飞行函数，d代表种群维度，x
r
(t)代表任意粒子的位置，r2为 均匀分布在0到1之间的随机数，s＝0.01，u与v为0到1之间的随机数，β＝1.5。
[0131]
具体地，在步骤s24中“采用大范围搜索攻击策略更新粒子种群的位置
”ꢀ
至少包括步骤s241至步骤s243。
[0132]
s241、获取第一随机数，并判断随机数是否小于第一预设值。
[0133]
s242、当判断到第一随机数小于第一预设值时，采用高空飞行搜索并垂直 俯身攻击的策略更新粒子种群的位置。具体地，高空飞行搜索并垂直俯身攻击 的位置更新模型为：
[0134][0135][0136]
其中，x1(t 1)表示粒子的下一个更新位置，x
best
(t)代表目前已取得的最佳位 置，x
m
(t)是当前所有粒子的平均位置，t为当前迭代次数，t为最大迭代次数， n为种群数量，r1为0到1之间的随机数，x
i
(t)为当前第i个粒子位置。
[0137]
s243、当判断到第一随机数不小于第一预设值时，采用环绕猎物飞行并短 距离滑翔攻击的策略更新粒子种群的位置。具体地，环绕猎物飞行并短距离滑 翔攻击的位置更新模型为：
[0138]
x2(t 1)＝x
best
(t)
×
lf(d) x
r
(t) (y
‑
x)
×
r2……………
(式4)
[0139][0140][0141][0142]
其中，x2(t 1)表示粒子的下一个更新位置，x
best
(t)代表目前已取得的最佳位 置，
lf(d)代表莱维飞行函数，d代表种群维度，x
r
(t)代表任意粒子的位置，r2为 均匀分布在0到1之间的随机数，s＝0.01，u与v为0到1之间的随机数，β＝1.5。
[0143]
s25、当迭代次数大于预设迭代次数时，采用局部攻击策略更新粒子种群的 位置，以及通过精英狩猎策略计算得到精英位置，然后通过贪婪策略获取最佳 位置，并将迭代次数加一，和判断是否达到最大迭代次数。需要说明的是，当 迭代次数大于预设迭代次数时，阿奎拉鹰模型进入开发阶段在开发阶段，阿奎 拉鹰执行局部攻击策略。通过两种不同的狩猎方式实现搜索空间内的局部搜索 猎物并攻击，两种狩猎方式的介绍如下：
[0144]
3、低空飞行并缓降攻击
[0145]
阿奎拉鹰在探索阶段已确定猎物位置，因此，阿奎拉鹰采用低空飞行并缓 降攻击猎物，利用此方法捕食逃跑能力较低的猎物，其位置更新公式如下：
[0146]
x3(t 1)＝(x
best
(t)
‑
x
m
(t))
×
α
‑
r4 ((ub
‑
lb)
×
r5 lb)
×
δ
[0147][0148]
其中，x3(t 1)表示粒子的下一个更新位置，x
best
(t)代表目前已取得的最佳位 置，x
m
(t)是当前所有粒子的平均位置，α与δ为适应参数，其值为0.1，ub是上 边界，lb为下边界，r4和r5为0到1之间的随机数。x
i
(t)为当前第i个粒子位置， n粒子的数量。
[0149]
4、地面近距离攻击猎物
[0150]
阿奎拉鹰在低空种根据猎物逃跑轨迹追逐猎物，并在地面上对猎物发动攻 击，其位置更新公式如下：
[0151][0152]
其中，x4(t 1)表示粒子的下一个更新位置，qf(t)是质量函数，用于平衡其 搜索策略，x
best
(t)代表目前已取得的最佳位置，g1表示阿奎拉鹰追踪猎物的移动 参数，为[
‑
1,1]之间的随机数，x(t)为当前的粒子位置，r6、r7和r8为0到1之间 的随机数，g2在空中追逐猎物的飞行斜率，从2到0线性递减，lf(d)代表莱维 飞行函数，t是当前迭代次数，t为最大迭代次数。
[0153]
具体地，在步骤s25中“采用局部攻击策略更新粒子种群的位置”至少包 括步骤s251至步骤s253。
[0154]
s251、获取第二随机数，并判断随机数是否小于第二预设值。
[0155]
s252、当判断到第二随机数小于第二预设值时，采用低空飞行并缓降攻击 的策略更新粒子位置。具体地，低空飞行并缓降攻击的位置更新模型为：
[0156]
x3(t 1)＝(x
best
(t)
‑
x
m
(t))
×
α
‑
r4 ((ub
‑
lb)
×
r5 lb)
×
δ
………
(式5)
[0157][0158]
其中，x3(t 1)表示粒子的下一个更新位置，x
best
(t)代表目前已取得的最佳位 置，
x
m
(t)是当前所有粒子的平均位置，α与δ为适应参数，其值为0.1，ub是上 边界，lb为下边界，r4和r5为0到1之间的随机数。x
i
(t)为当前第i个粒子位置， n粒子的数量。
[0159]
s253、当判断到第二随机数不小于第二预设值时，采用地面近距离攻击猎 物的策略更新粒子位置。具体地，地面近距离攻击猎物的位置更新模型为：
[0160][0161]
其中，x4(t 1)表示粒子的下一个更新位置，qf(t)是质量函数，用于平衡其 搜索策略，x
best
(t)代表目前已取得的最佳位置，g1表示阿奎拉鹰追踪猎物的移动 参数，为[
‑
1,1]之间的随机数，x(t)为当前的粒子位置，r6、r7和r8为0到1之间 的随机数，g2在空中追逐猎物的飞行斜率，从2到0线性递减，lf(d)代表莱维 飞行函数，t是当前迭代次数，t为最大迭代次数。
[0162]
为了防止改进的阿奎拉鹰模型，在开发阶段陷入局部最优解，而导致过早 收敛。在本实施例中，通过精英狩猎策略，利用精英个体与普通个体进行局部 开发，避免阿奎拉鹰模型早熟收敛。具体地，通过精英狩猎策略计算得到精英 位置的精英位置计算模型为：
[0163]
x
e
(t 1)＝x
e
‑
best
cd
×
(x(t)
‑
x
e
‑
normal
(t)) σ
×
(x
rand1
‑
x
rand2
)
[0164][0165]
cd＝z0 0.1
×
tan(π
×
(r
‑
0.5))
[0166]
其中，x
e
(t 1)是基于精英狩猎策略产生的新位置，x
e
‑
best
为从5个精英个体 中随机选出的个体，cd为符合柯西分布的随机数，x(t)为当前的粒子位置， x
e
‑
normal
(t)为从仓库中随机选出的15个位置，x
rand1
与x
rand2
为从剩下的15个个体 中随机选出的，exponent＝2，σ
initial
＝1σ
final
＝0，z0＝2，t为当前迭代次数，t为 最大迭代次数，r为分布在[0,1]之间的随机数。
[0167]
需要说明的是，精英狩猎策略通过建立一个种群库，选派个体精英外出探 索，并使个体之间充分交互信息，避免算法陷入局部最优，使算法不易早熟收 敛。
[0168]
精英狩猎策略原理：首先，建立一个仓库，在每次迭代中，根据取出15个 种群放入仓库中，并根据适应度值排序，确定精英个体(5个)，其余的为普通 个体(10个)，此时，若总的种群数量为30个，则仓库外的个体共有15个。按 照精英位置计算模型进行计算。
[0169]
具体地，通过贪婪策略获取最佳位置为本领域技术人员的常规技术手段， 本发明在此不再赘述。
[0170]
s26、当判断到未达到最大迭代次数时，计算更新后粒子种群中每个粒子的 适应度，并判断迭代次数是否不大于预设迭代次数。否则，输出最大迭代次数 时的最佳位置及其对应的适应度，以获得设计参数。
[0171]
本实施例提供的模型参数的计算方法，采用了改进的阿奎拉鹰模型，通过 在种群初始化阶段引入logistic混沌映射，提高种群分布多样性。在局部开发 阶段引入精英狩猎
策略，建立种群库以充分发挥精英个体与普通个体的潜在能 力，通过贪婪算法选择最佳位置进入下一次迭代，避免算法易于早熟收敛。最 后，通过惩罚函数将约束设计问题转换为无约束优化问题从而能够使用iao算 法求得最优解。能够获得更为准确的模型参数的设计参数。
[0172]
第二方面、
[0173]
本发明实施例提供了一种模型参数的计算装置，其包含：
[0174]
模型获取模块1，用于获取待求解的目标对象的计算模型。优选地，目标对 象为压力容器。
[0175]
模型求解模块2，用于通过改进的阿奎拉鹰模型对计算模型进行求解，以获 得目标对象的设计参数。其中，改进的阿奎拉鹰优化模型通过logistic混沌映射 计算初始种群，且在开发阶段采用精英狩猎策略进行局部开发，以避免过早收 敛。
[0176]
本实施例提供的模型参数的计算装置，采用了改进的阿奎拉鹰模型，通过 在种群初始化阶段引入logistic混沌映射，提高种群分布多样性。在局部开发 阶段引入精英狩猎策略，建立种群库以充分发挥精英个体与普通个体的潜在能 力，通过贪婪算法选择最佳位置进入下一次迭代，避免算法易于早熟收敛。最 后，通过惩罚函数将约束设计问题转换为无约束优化问题从而能够使用iao算 法求得最优解。能够获得更为准确的模型参数的设计参数。
[0177]
在一个可选地实施例中，模型求解模块2，具体包括：
[0178]
初始化单元，用于根据计算模型，随机生成初始粒子种群，并初始化迭代 次数和最大迭代次数。
[0179]
混沌映射单元，用于根据初始粒子种群，通过logistic混沌映射模型，计算 得到多样化粒子种群。
[0180]
适应度单元，用于计算多样化粒子种群中每个粒子的适应度，并判断迭代 次数是否不大于预设迭代次数。优选地，预设迭代次数为最大迭代次数的三分 之二。
[0181]
大范围单元，用于当迭代次数不大于预设迭代次数时，采用大范围搜索攻 击策略更新粒子种群的位置，并将迭代次数加一，和判断是否达到最大迭代次 数。
[0182]
局部单元，用于当迭代次数大于预设迭代次数时，采用局部攻击策略更新 粒子种群的位置，以及通过精英狩猎策略计算得到精英位置，然后通过贪婪策 略获取最佳位置，并将迭代次数加一，和判断是否达到最大迭代次数。
[0183]
输出单元，用于当判断到未达到最大迭代次数时，计算更新后粒子种群中 每个粒子的适应度，并判断迭代次数是否不大于预设迭代次数。否则，输出最 大迭代次数时的最佳位置及其对应的适应度，以获得设计参数。
[0184]
在一个可选地实施例中，初始化单元，具体用于：
[0185]
根据所述计算模型，获取搜索空间的上界和下界。
[0186]
根据随机函数，在搜索空间内生成所述初始粒子种群。具体地，初始种群 生成模型为x＝rand
×
(ub
‑
lb) lb，其中，x为粒子位置，rand为均匀分布 在0到1之间的随机数，ub为搜索空间的上边界，lb为搜索空间的下边界。
[0187]
在一个可选地实施例中，logistic混沌映射模型为：
[0188]
x(t 1)＝a
×
x(t)
×
(1
‑
x(t))
[0189]
其中，x(t 1)为更新后的粒子位置，a＝4，x(t)为当前的粒子位置。
[0190]
在一个可选地实施例中，粒子的适应度的适应度计算模型为：
[0191][0192]
其中，f(x)是惩罚函数值，f(x)是原函数值，m是惩罚因子,g
i
(x)为约束条 件，n为粒子数量。
[0193]
在一个可选地实施例中，采用大范围搜索攻击策略更新粒子种群的位置时， 大范围单元具体用于：
[0194]
获取第一随机数，并判断随机数是否小于第一预设值。
[0195]
当判断到所述第一随机数小于第一预设值时，采用高空飞行搜索并垂直俯 身攻击的策略更新粒子种群的位置。具体地，高空飞行搜索并垂直俯身攻击的 位置更新模型为：
[0196][0197][0198]
其中，x1(t 1)表示粒子的下一个更新位置，x
best
(t)代表目前已取得的最佳位 置，x
m
(t)是当前所有粒子的平均位置，t为当前迭代次数，t为最大迭代次数， n为种群数量，r1为0到1之间的随机数，x
i
(t)为当前第i个粒子位置。
[0199]
当判断到所述第一随机数不小于第一预设值时，采用环绕猎物飞行并短距 离滑翔攻击的策略更新粒子种群的位置。具体地，环绕猎物飞行并短距离滑翔 攻击的位置更新模型为：
[0200]
x2(t 1)＝x
best
(t)
×
lf(d) x
r
(t) (y
‑
x)
×
r2[0201][0202][0203][0204]
其中，x2(t 1)表示粒子的下一个更新位置，x
best
(t)代表目前已取得的最佳位 置，lf(d)代表莱维飞行函数，d代表种群维度，x
r
(t)代表任意粒子的位置，r2为 均匀分布在0到1之间的随机数，s＝0.01，u与v为0到1之间的随机数，β＝1.5。
[0205]
在一个可选地实施例中，采用局部攻击策略更新粒子种群的位置时，局部 单元具体用于：
[0206]
获取第二随机数，并判断随机数是否小于第二预设值。
[0207]
当判断到所述第二随机数小于所述第二预设值时，采用低空飞行并缓降攻 击的策略更新粒子位置。具体地，低空飞行并缓降攻击的位置更新模型为：
[0208]
x3(t 1)＝(x
best
(t)
‑
x
m
(t))
×
α
‑
r4 ((ub
‑
lb)
×
r5 lb)
×
δ
[0209][0210]
其中，x3(t 1)表示粒子的下一个更新位置，x
best
(t)代表目前已取得的最佳位 置，x
m
(t)是当前所有粒子的平均位置，α与δ为适应参数，其值为0.1，ub是上 边界，lb为下边界，r4和r5为0到1之间的随机数。x
i
(t)为当前第i个粒子位置， n粒子的数量。
[0211]
当判断到所述第二随机数不小于所述第二预设值时，采用地面近距离攻击 猎物的策略更新粒子位置。具体地，地面近距离攻击猎物的位置更新模型为：
[0212][0213]
其中，x4(t 1)表示粒子的下一个更新位置，qf(t)是质量函数，用于平衡其 搜索策略，x
best
(t)代表目前已取得的最佳位置，g1表示阿奎拉鹰追踪猎物的移动 参数，为[
‑
1,1]之间的随机数，x(t)为当前的粒子位置，r6、r7和r8为0到1之间 的随机数，g2在空中追逐猎物的飞行斜率，从2到0线性递减，lf(d)代表莱维 飞行函数，t是当前迭代次数，t为最大迭代次数。
[0214]
在一个可选地实施例中，通过精英狩猎策略计算得到精英位置的精英位置 计算模型为：
[0215]
x
e
(t 1)＝x
e
‑
best
cd
×
(x(t)
‑
x
e
‑
normal
(t)) σ
×
(x
rand1
‑
x
rand2
)
[0216][0217]
cd＝z0 0.1
×
tan(π
×
(r
‑
0.5))
[0218]
其中，x
e
(t 1)是基于精英狩猎策略产生的新位置，x
e
‑
best
为从5个精英个体 中随机选出的个体，cd为符合柯西分布的随机数，x(t)为当前的粒子位置，x
e
‑
normal
(t)为从仓库中随机选出的15个位置，x
rand1
与x
rand2
为从剩下的15个个体 中随机选出的，exponent＝2，σ
initial
＝1σ
final
＝0，z0＝2，t为当前迭代次数，t为 最大迭代次数，r为分布在[0,1]之间的随机数。
[0219]
第三方面、
[0220]
本发明实施例提供了一种模型参数的计算设备，其包括处理器、存储器， 以及存储在存储器内的计算机程序。计算机程序能够被处理器执行，以实现如 第一方面所说的模型参数的计算方法。
[0221]
第四方面、
[0222]
本发明实施例提供了一种计算机可读存储介质，计算机可读存储介质包括 存储
的计算机程序，其中，在计算机程序运行时控制计算机可读存储介质所在 设备执行如第一方面所说的模型参数的计算方法。
[0223]
在本发明实施例所提供的几个实施例中，应该理解到，所揭露的装置和方 法，也可以通过其它的方式实现。以上所描述的装置和方法实施例仅仅是示意 性的，例如，附图中的流程图和框图显示了根据本发明的多个实施例的装置、 方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程 图或框图中的每个方框可以代表一个模块、程序段或代码的一部分，所述模块、 程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。 也应当注意，在有些作为替换的实现方式中，方框中所标注的功能也可以以不 同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地 执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意 的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合， 可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专 用硬件与计算机指令的组合来实现。
[0224]
另外，在本发明各个实施例中的各功能模块可以集成在一起形成一个独立 的部分，也可以是各个模块单独存在，也可以两个或两个以上模块集成形成一 个独立的部分。
[0225]
所述功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时， 可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方 案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件 产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指 令用以使得一台计算机设备(可以是个人计算机，电子设备，或者网络设备等) 执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：u 盘、移动硬盘、只读存储器(rom，read
‑
only memory)、随机存取存储器(ram， random access memory)、磁碟或者光盘等各种可以存储程序代码的介质。需要 说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非 排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包 括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、 方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括 一个
……”
限定的要素，并不排除在包括所述要素的过程、方法、物品或者设 备中还存在另外的相同要素。
[0226]
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领 域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则 之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之 内。