基于小生境遗传算法的高速动车组稳定性优化与灵敏度分析

自国内第一条高速铁路秦沈客运专线建成以来，经过引进-消化-吸收再创新，高速铁路技术得到了迅速发展。目前我国拥有世界最大规模与最高运营速度铁路网，并预计至2025年,铁路网规模达到17.5万km以上，其中高速铁路3.8万km左右[1], 基本完成了“四纵四横”等客运专线规划[2]。

伴随着高速铁路的发展与运用，车辆的稳定性要求越来越高，所暴露的问题也愈发严重，部分车辆出现了低锥度晃车和高锥度转向架失稳或抖车问题，即一次蛇行失稳与二次蛇行失稳问题[3]。如何对稳定性进行全面优化，达到同时兼顾车体稳定性和构架稳定性的目标，成为高速动车组发展的重要问题之一。

国内外针对铁道车辆的稳定性问题进行了大量的研究。Ahmadian[4] 在考虑轮轨接触和抗蛇行减振器非线性特性的基础上，研究了系统参数对转向架稳定性的影响，结果表明在一定范围内增大抗蛇行减振器阻尼反而会降低临界速度。Huang等[5]利用根轨迹法分析了车辆系统参数对车体一次蛇行稳定性的影响，研究结果表明轮轨接触参数和车辆悬挂参数对车体稳定性起决定性作用。孙建锋等[6]借助轮对运动参数的能量表达式，通过计算不同参数条件下的输入能量，对比了踏面等效锥度、轮对质量、一系悬挂刚度等参数对轮对稳定性的影响规律。此外，基于蛇行运动的稳定性理论，通过推导带抗蛇行减振器的转向架线性临界速度解析表达式，研究了不同等效锥度下抗蛇行减振器串联刚度和结构阻尼对临界速度的影响规律[7]。

在小鼠中，ACE基因完全失活是可存活的，而在人类则导致严重低血压、肾灌注不足和RTD，这是一种以出生时无尿或严重肾功能不全、缺乏肾近端小管分化和常见的围产期死亡为特征的严重肾脏疾病[5]。目前人类ACE基因突变的临床病例致病性报道比较有限。有报道称ACE基因的纯合突变或复合杂合突变会导致RTD的发生，临床表现主要为严重的羊水过少、出生时的肾功能不全及出生后短时间内的死亡[6]。在某些情况下，肾功能不全者需要出生之后的长期腹膜透析，例如发生Q1069R错义突变，导致ACE活性完全缺乏[7,8]。

上述研究仅是针对车体稳定性或构架稳定性开展的单参数单目标的优化分析，而参数对稳定性的影响具有交互作用，且部分参数针对车体稳定性和构架稳定性的优化方向不同，因此有必要开展多参数多目标优化。遗传算法对于求解非线性、多参数、多目标的复杂系统优化问题具有很好的效果，且在铁道车辆领域也有一定的应用基础。汤东胜等[8]将遗传算法应用于机车车辆横向动力学，以临界速度为衡量指标，平稳性指标为约束，设计了适用于动力学多参数的优化方案。沈文林[9]采用正交参数灵敏度分析方法与改进小生境遗传算法，对动车组悬挂参数进行了合理匹配，优化得到良好的动力学性能；Jiang等[10]采用遗传算法对铰接式单轨车辆系统的曲线通过动力学性能进行了悬挂参数的优化分析，目标函数中主要考虑了乘坐舒适性、轮重减载率、倾覆系数以及磨耗指数。

我们已赤裸地倒在她的床上，我的唇正要印上她那如花一样绽放的唇——忽然间，我感觉到一种不对的地方。我停住，仔细地观察、倾听、嗅闻……她仰卧在我的身边，那么完美，膝陇，热情，期待。她没有任何不对的地方。我猛然意识到，不对劲的是自己。

目前，针对车体稳定性与构架稳定性的多目标多参数优化研究还很少。本文将采用自适应小生境遗传算法，对高速动车组综合稳定性进行优化，并对多参数综合影响下的参数灵敏度进行分析，研究各参数对稳定性的影响程度。

1 高速动车组动力学模型与稳定性特征

以国内某250 km/h级别的高速动车组为例，建立整车动力学模型，见图1。车辆由1个车体、2个构架、4个轮对与8个转臂轴箱组成。车体、构架、轮对均考虑伸缩、横移、浮沉、侧滚、点头和摇头6个自由度，转臂轴箱考虑点头自由度。模型总自由度为50个，垂向运动和横向运动耦合在一起。模型中建立了一系簧、一系垂向减振器、转臂节点、二系空簧、牵引拉杆、抗侧滚扭杆、二系横向减振器、抗蛇行减振器等力元。轮轨接触采用Hertz接触与Kalker简化理论进行计算，摩擦系数0.4。车辆模型自由度见表1，部分主要参数见表2。

图1 高速动车组整车动力学模型

表1 整车模型自由度

自由度纵向横向垂向侧滚点头摇头车体∨∨∨∨∨∨构架∨∨∨∨∨∨轮对∨∨∨∨∨∨轴箱∨

表2 高速动车组主要参数

参数取值参数取值轮对内侧距/m1.353轴距/m2.5定距/m17.8车轮滚动圆半径/m0.46轮对质量/t1.56构架质量/t2.5车体质量/t39转臂节点纵向刚度/(MN·m-1)15一系垂向刚度/(MN·m-1)2一系垂向阻尼系数/[kN·(m·s-1)-1]30二系水平刚度/(MN·m-1)0.2 二系垂向刚度/(MN·m-1)0.3二系横向阻尼/[kN·(m·s-1)-1]50抗蛇行减振器等效刚度/(MN·m-1)8抗蛇行卸荷速度/(m·s-1)0.006抗蛇行卸荷力/kN8

本文针对高速动车组的稳定性问题进行全面优化。高速动车组稳定性问题包括在实际运营过程中出现的一次蛇行失稳和二次蛇行失稳现象。一次蛇行失稳是指转向架蛇行与车体模态耦合现象，二次蛇行失稳是指转向架的蛇行失稳现象。为对该现象进行再现，选取3种轮轨匹配关系进行研究：标准踏面匹配过度打磨钢轨，对应锥度0.022低锥度，模拟一次蛇行现象；标准踏面匹配标准轨面，对应锥度0.038，模拟正常运行；实测磨耗到限踏面匹配标准轨面，对应锥度0.38高锥度，模拟二次蛇行失稳现象。3种轮轨关系的锥度图见图2。等效维度分别采用简化法、谐波法及UIC 519—2004[11]进行计算，结果见表3。

在大数据时代，对于国家审计部门而言，应当主动分析形势，加强互联网、物联网、云计算等现代信息技术的应用探索，借助现代信息技术，与审计工作深度融合，优化审计资源配置，提高动态监控水平和审计监管的连续性，加强审计整改跟踪，才能切实发挥审计应有的功能，具体建议如下。

食品工业作为改革开放过程中的受益者，在国民经济中占有重要位置。我国食品工业经历了哪几个阶段？各阶段得到哪些政策的扶植与推动？科技创新为食品工业发展起到哪些作用？消费者的需求又发生哪些变化？一个个问题的答案，正是食品工业飞速发展的佐证。

图2 不同轮轨关系的等效锥度

表3 不同方法计算所得等效锥度值

锥度分类计算方法简化法谐波法UIC 519—2004[11]低锥度左0.025,右0.0250.0220.023正常锥度左0.04,右0.040.0380.039高锥度左0.365,右0.3650.3820.388

在所采用的3种轮轨匹配条件下，对高速动车组进行稳定性分析，其分岔图见图3。改变一定参数模拟一次蛇行。从图3中可以看出，低锥度条件下，在速度为150 km/h之后会出现明显的一次蛇行现象；正常锥度条件下临界速度超过500 km/h；高锥度条件下临界速度为450 km/h。

图3 不同锥度下分岔图

2 小生境遗传算法原理与目标函数选取

为了对高速动车组的稳定性问题进行多目标多参数优化，选取自适应的小生境遗传算法进行实现。

这一天，有一半时间他一直躺着不动，尽力和昏迷斗争，等着那个要把他吃掉、而他也希望能吃掉的东西。有时候，疲倦的浪潮涌上来，淹没了他，他会做起很长的梦；然而在整个过程中，不论醒着或是做梦，他都在等着那种喘息和那条粗糙的舌头来舐他。

2.1 小生境遗传算法原理

遗传算法是模拟自然选择和遗传机制，采用进化形式达到最终目标的随机搜索技术。遗传算法的运算流程包括编码、初始种群、适应度、选择、交叉、变异六部分。整个计算过程为模拟自然选择、生物交配与生物变异过程，对随机生成的N个个体所组成的种群按照适应度进行选择，以交叉概率Pc随机交叉运算，以变异概率Pm随机选取个体中的基因位进行变异，从而获得适应度较高且具有一定多样性的种群。

在多目标遗传算法中，有向量评估遗传算法(VEGA)、非劣分层遗传算法(NSGA)、小生境Pareto遗传算法(NPGA),基于距离的Pareto遗传算法(DPGA)等等。在各种算法中，小生境遗传算法不仅种群多样性与全局优化效果较好，而且计算量较为经济[12]。

小生境遗传算法是基于生物倾向与相同个体或相似个体聚集的现象所提出的一种遗传算法。在优化过程中，将个体按照相似程度划分为若干类，在每一类中根据个体适应度挑选优秀个体再组成群[13]。具体做法为计算相似个体的海明距离，并对其中适应度较低的个体进行惩罚。

传统的小生境遗传算法需要预先设定好交叉概率Pc、变异概率Pm和小生境距离参数L，这些参数选取依赖经验并且不会随着进化过程而变化，降低运算效率。自适应的小生境遗传算法对这些问题进行了改进，使得小生境遗传算法的全局搜索能力和快速收敛能力得到最大的发挥。主要改进简要说明如下：

(1)动态的距离参数[14] 将个体间的距离判别参数设置为与进化代数有关的动态函数，保证种群进化初期个体间具有很大的差异性，也能保证在进化后期优良个体的相似性。

（3）拱部小导管注浆预加固。参考方案中的设计轮廓线，在其外侧布置小导管，从而为注浆提供通道，以便对顶部围岩区域进行预加固处理。小导管的间距有所讲究，彼此之间应达到30cm，同时每一根管的外倾角都应介于5～10°范围内，管长以4～6m为宜。运用风钻成孔的方法，将孔径控制在5Omm水平。当完成小导管安装作业后便可进行注浆处理，整个过程应在高压环境中进行，控制水灰比介于1：1～0.5：1范围内，讲求先稀后浓的原则，灌浆压力≥0.3MPa。当孔内不再吃浆时，应进行闭浆处理，整个过程时长≥5min，对应压力应＞O.5MPa。在严重渗水区域施工时，可采用双浆液注浆法，此时各管搭接长度≥1m。

(2)格雷码编码[15] 在格雷码下，海明距离可以很好的反应实际参数差异性。

(3)自适应交叉算子和变异算子[16] 在进化过程中根据个体适应度自动调整交叉概率和变异概率。

(4)精英保留策略[17] 将每一代中产生的最优秀个体进行保存直接进入下一代，特点不会被破坏。

2.2 目标函数选取与计算流程

采用加权求和的方式将多目标优化问题转化为单目标优化问题。对稳定性的优化分为一次蛇行性能优化与二次蛇行性能优化。此外由于涉及车体蛇行问题，还需要同时保证车辆的平稳性与舒适度。在评价函数设置中，采用临界速度，在速度200、250 km/h下0.3～3 Hz带通滤波后车体横向加速度最大值、横向平稳性、舒适度四类内容进行评估。因高速列车线路条件较好，曲线半径较大，曲线通过放在次要地位。临界速度采用轮对横移量大于0.5 mm失稳进行判断。

以各项指标的规定限值为基准，对各项动力学性能指标进行归一化处理。临界速度与其他指标取值方向相反，并且有一定约束条件：按照裕量设置速度最小为300 km/h，低于300 km/h个体予以惩罚；为防止性能过剩，临界速度最高设置为550 km/h，高于550 km/h临界速度均视为550 km/h。各项指标均取25%权重，最终得到的评价函数为

Obj=-Vcr/400+Spy/2.5+Com/2+Acc/0.2

(1)

式中：Obj为优化使用的评价函数值；Vcr为临界速度(采用速度400 km/h为归一化系数)；Spy为横向平稳性值；Com为舒适度值；Acc为0.3～3 Hz车体横向加速度最大值(根据原车计算结果取0.2 m/s2为归一化系数)。

通过上述方式，将高速动车组稳定性优化问题转化为了多参数与Obj最小值优化问题。在计算目标函数时，分别考虑异常低锥度、正常锥度和极限高锥度3种轮轨匹配关系下的动力学模型。以保证正常运营工况为主，异常低锥度与高锥度2种稳定性极限工况综合性能最优为目标；正常锥度考虑50%权重，其他锥度考虑为25%权重，作为对应参数综合目标值。对应自适应小生境遗传算法下稳定性优化计算流程见图4。

图4 自适应小生境算法稳定性优化流程图

3 稳定性优化与灵敏度分析

根据自适应小生境算法，对所选高速动车组进行稳定性优化。以原始参数为中心，结合各悬挂元件常用参数范围，选取横向减振器等效刚度K_Csy (3～10 MN/m)、阻尼Csy (20～80 kN/(m·s-1))，抗蛇行减振器等效刚度K_Csx (3～12 MN/m)、卸荷速度dV (0.002～0.02 m/s)、卸荷力dF(3～12 kN)共5个对稳定性影响较大的参数进行优化。

初始种群随机个体取25，共迭代15代。精英个体数为3个，初始小生境参数L为8.4。适应度最低个体的交叉概率Pc、变异概率Pm分别为0.7、0.1。进行小生境遗传算法迭代计算。计算代数目标函数最小值变化见图5。

图5 目标函数最小值变化

最终，本文经迭代后获得了该高速列车稳定性优化的最优解与次优解，见表4。

X-431 PAD III 是元征公司研发的一款基于Android 系统，支持“智能诊断”和“在线编程”的高端汽车故障诊断设备。本文将介绍使用X-431 PAD III 进行2015年宝马528Li全景摄像机系统标定和学习。

他预感到，在情感上有着某种相似孤独感的伍亦苒同样也在期待，期待在故事的中段因为情节激荡而生发出的高潮。

表4 小生境遗传算法最终优化参数

参数组合K_Csy /(MN·m-1)Csy/(kN·s·m-1)K_Csx/(MN·m-1)dV /(m·s-1)dF /N组合15.426.77.30.003 39 900组合23.642.98.40.003 310 700组合36.137.2100.005 211 700

3种方案最终目标函数值相近，代表了不同的性能倾向。对3种参数组合下，高速动车组的主要性能进行对比分析。各等效锥度下的分岔图与平稳性见图6、图7。分岔图采用升速结果进行各组参数对比。

扭矩限制器预紧螺栓上紧时的摩擦系数，以及摩擦片的摩擦系数对扭矩限制器装配后的实际打滑扭矩有很大的影响。如果上述摩擦系数离散度比较大，会造成扭矩限制器装配后的实际打滑扭矩与理论值偏差很大。实际打滑扭矩比要求值偏小较多，会造成联轴器频繁打滑；实际打滑扭矩比要求值偏大，又起不到过载保护的作用。

图6 各参数组合分岔图对比

表5为提取的失稳速度与失稳后的轮对最大横移量。由图6与表5可以看出，3组参数组合在稳定性上表现各异，一次蛇行稳定性参数选取与二次蛇行出现了一定的矛盾。正常锥度下各参数组合稳定性相似；低锥度条件下，参数组合1整体轮对横移量最小，参数组合3最大，但各组参数都优于原车参数；高锥度条件下，参数组合1临界速度最低，而参数组合3临界速度最高。可见，参数组合1条件下偏向于一次蛇行优化，参数组合3偏向于二次蛇行优化，而参数组合2介于了两者之间。

表5 各参数组合下临界速度与失稳后最大横移量

参数低锥度正常锥度高锥度临界速度/(km·h-1)横移量/mm临界速度/(km·h-1)横移量/mm临界速度/(km·h-1)横移量/mm原始参数1401.6955004608.97参数组合11301.2455004408.2参数组合21601.3455004608.19参数组合31401.5455004807.97

由图7平稳性分析可知，各参数组合的横向平稳性得到了明显的优化。在各锥度条件下，各参数组合平稳性都优于原车参数，其中参数组合1最优，参数组合2与参数组合3横向平稳性相似。

其次，作者对翻译过程原创性地进行了回溯性研究。书中对诗歌翻译的体验性进行了描述，将诗歌翻译与诗歌创作相提并论，认为诗歌翻译不是静态地复制原诗，而是一种诗学挪用和创造，是译者诗学能力和诗学爱好的展示。书中既有形而上的考察，也有形而下的举例，涵括了诗歌翻译中对本体的认识、对原作的解读、译作中的语言表达等阶段，无论是在原创性上还是系统性上都值得称道。

图7 各参数组合平稳性对比

上述参数优化中，稳定性与平稳性均随着参数变化表现出了一定的趋势性。为了能够更加清晰的反映所选取的各参数对于稳定性与平稳性各项性能的影响，进行参数灵敏度计算分析。统计所有的计算工况，将各参数与各项指标对应关系列出，见图8～图12。

对所有数据进行线性拟合，提取拟合系数。低锥度工况关注200 km/h横向车体加速度，正常锥度与高锥度工况关注临界速度。各工况都关注横向平稳性。

采用各参数的百分比差分灵敏度进行各参数影响对比

(2)

式中：Sp(i)为第i个参数值下的百分比灵敏度；xi和xi+1分别为第i和第i+1个参数值；yi和yi+1分别为第i和第i+1个参数值对应的指标值。

提取到的各参数对各项指标的百分比灵敏度最大值见表6。

表6 各参数对各项指标的百分比灵敏度最大值 %

参数车体横向加速度/(m·s-2)低锥度临界速度/(km·h-1)横向平稳性正常锥度高锥度低锥度正常锥度高锥度K_Csy /(MN·m-1)6.47-5.38-17.030.150.194.45Csy /(kN.s·m-1)35.81-15.8629.5125.4924.0228.04K_Csx/(MN·m-1)52.705.7421.636.648.02-5.44dV/(m·s-1)67.60-17.0611.2317.0818.126.56dF/N-97.4738.0943.10-15.86-14.68-14.39

由表6可见：

（2）使用建筑和装修去完成一体化的设计，以及施工，理其在想状态下可以说是装修能够伴随着主体施工去同步完成的；

(1)对于所选用的车型，在所选取的参数优化范围内，K_Csy的增大会使高锥度轮轨关系临界速度有所降低；Csy的增大会恶化车体横向加速度(35.81%)并致使横向平稳性(28.04%)更加恶劣，但高锥度临界速度会有所增加(29.51%)；K_Csx增大，车体横向加速度会更加恶劣(52.70%)，但临界速度有所增加(21.63%)；dV增大会恶化车体横向加速度(67.60%)与正常锥度下的临界速度(-17.06%)；dF增大，会优化车体横向加速度(-97.47%)与临界速度(38.09%)。

(2)对于所选用的车型，在所选取的参数优化范围内，对于车体一次蛇行，最为有效的措施为增大dF，百分比达到-97.10%，其次为减小K_Csx、dV，分别为52.70%、67.60%；对于正常轮轨关系临界速度，最为有效的措施同样为增大dF，百分比达到38.09%；对于高锥度轮轨关系临界速度，增大dF同样最为有效，百分比达到43.10%，其次为增大Csy，为29.51%；对于横向平稳性，影响最大为Csy，灵敏度最高达到28.04%，建议在保证稳定性的前提下取3 000～4 000 N/(m/s)较小值。

图8 K_Csy(横向减振器等效刚度)对不同目标的影响规律与灵敏度统计拟合

图9 Csy(横向减振器阻尼)对不同目标的影响规律与灵敏度统计拟合

图10 K_Csx(抗蛇行减振器等效刚度)对不同目标的影响规律与灵敏度统计拟合

图11 dV(抗蛇行减振器卸荷速度)对不同目标的影响规律与灵敏度统计拟合

图12 dF(抗蛇行减振器卸荷力)对不同目标的影响规律与灵敏度统计拟合

4 结论

本文通过建立国内某250 km/h级别的高速动车组动力学模型，采用自适应的小生境遗传算法，对车辆的一次蛇行与二次蛇行稳定性问题进行了多目标多参数优化，并对各参数灵敏度进行了分析。通过对比分析，可以得到以下结论：

(1)考虑一次蛇行条件下车体横向加速度，二次蛇行条件下临界速度，兼顾列车运行平稳性与舒适度，对异常低锥度、正常锥度与高锥度3种轮轨关系进行考虑，将各目标值进行综合，给出了稳定性问题优化的目标函数。

(2)对横向减振器等效刚度、阻尼，抗蛇行减振器等效刚度、卸荷速度、卸荷力共5个对稳定性影响较大的参数进行优化。选取算法参数，进行优化迭代。最终获得目标函数值相似的3组优化参数组合。3组参数分别偏向于一次蛇行，偏向于二次蛇行优化与介于两者之间。各参数综合指标都优于原车参数。以速度250 km/h运行为例，优化后，低锥度下横移量由1.6 mm降低为1 mm左右，高锥度下横向平稳性由2.80降为2.48左右。

(3)分析各参数百分比灵敏度。横向减振器阻尼主要影响横向平稳性，抗蛇行减振器等效刚度主要影响一次蛇行加速度，抗蛇行减振器卸荷力与卸荷速度对一次蛇行与二次蛇行都影响较大。对于所选用的车型，在所选取的参数优化范围内，优化车辆稳定性最为有效的措施为增大卸荷力。

参考文献：

[1] 中华人民共和国国家发展和改革委员会. 中长期铁路网规划[EB/OL]. (2016-07-13)[2019-12-17].https://www.ndrc.gov.cn/xxgk/zcfb/ghwb/201607/t20160720_962188_ext.html.

[2] 搜狐网. 四纵四横高铁网最后一“横”收官 [EB/OL]. (2018-05-14)[2019-12-17].http://www.sohu.com/a/231565908_618586.

[3] 孟素英. 高速列车蛇行运动半主动控制研究[D].成都：西南交通大学,2018.

[4] AHMADIAN M, YANG S P. Effect of System Nonlinearities on Locomotive Bogie Hunting Stability [J]. Vehicle System Dynamics, 1998, 29(6): 365-384.

[5] HUANG C H, ZENG J, LIANG S L. Carbody Hunting Investigation of a High Speed Passenger Car[J]. Journal of Mechanical Science and Technology,2013,27(8):2283-2292.

[6] 孙建锋,池茂儒,吴兴文,等.基于能量法的轮对蛇行运动稳定性[J].交通运输工程学报,2018,18(2):82-89.

SUN Jianfeng, CHI Maoru, WU Xingwen, et al. Hunting Motion Stability of Wheelset Based on Energy Method [J]. Journal of Traffic and Transportation Engineering, 2018, 18(2): 82-89.

[7] 孙建锋,池茂儒,吴兴文,等.抗蛇行减振器参数对车辆稳定性的影响分析[J].振动.测试与诊断,2018,38(6):1155-1160，1291.

SUN Jianfeng, CHI Maoru, WU Xingwen, et al. Analysis of the Influence of Anti-snaking Shock Absorber Parameters on Vehicle Stability [J]. Journal of Vibration. Measurement & Diagnosis, 2018, 38(6): 1155-1160,1291.

[8] 汤东胜,范佩鑫.基于遗传算法的机车车辆横向动力学多参数优化设计[J].内燃机车,2000(12):9-13,1.

TANG Dongsheng, FAN Peixin. Multi-parameter Optimization with Genetic Algorithm for Lateral Dynamics of Railway Vehicles [J]. Diesel Locomotives, 2000(12): 9-13,1.

[9] 沈文林. 动车组参数灵敏度分析及参数优化[D].成都：西南交通大学,2013.

[10] JIANG Y Z, WU P B, ZENG J, et al. Multi-parameter and Multi-objective Optimisation of Articulated Monorail Vehicle System Dynamics Using Genetic Algorithm [J]. Vehicle System Dynamics, 2020, 58(1):74-91.

[11] International Union of Railways. Method for Determining the Equivalent Conicity: UIC 519—2004 [S]． 1st ed． Paris: International Union of Railways，2004．

[12] 徐磊. 基于遗传算法的多目标优化问题的研究与应用[D]. 长沙：中南大学, 2007.

[13] 沈文林,池茂儒.改进的小生境遗传算法在铁道车辆优化设计中的应用[J].机械,2014,41(1):5-8,31.

SHEN Wenlin, CHI Maoru. Improved Niche Genetic Algorithm in Optimization Design of Railway Vehicle [J]. Machinery, 2014, 41(1): 5-8, 31.

[14] LI B H, CHOI S S, VILATHGAMUWA D M. A New Control Strategy For Energy-saving Dynamic Voltage Restoration[C]//Power Engineering Society Summer Meeting. New York: IEEE, 2000.

[15] 雷英杰, 张善文. MATLAB遗传算法工具箱及应用[M]. 西安:西安电子科技大学出版社, 2014.

[16] 黄聪明,陈湘秀.小生境遗传算法的改进[J].北京理工大学学报,2004,24(8):675-678.

HUANG Congming, CHEN Xiangxiu. Improvements on Niche Genetic Algorithm [J]. Transactions of Beijing Institute of Technology, 2004, 24(8): 675-678.

[17] 王小平, 曹立明. 遗传算法——理论、应用与软件实现[M]. 西安:西安交通大学出版社, 2002.