粒子群优化算法分数阶PID控制器精品文档完整版Word下载.docx

1、其二，本文先采用工程上的整定方法（临界比例度法）粗略的确定其初始的三个参数，再利用粒子群算法进行寻优，得到更好的PID参数；其三，采用SIMULINK的仿真工具对PID参数优化系统进行仿真，得出系统的响应曲线。从中发现它的性能指标，都比原来有了很大的改进。因此，采用粒子群算法的优越性是显而易见的。关键词 目标函数；PID参数；粒子群算法；优化设计；SIMULINK第 1 章 绪论1.1 研究背景和课题意义在现代工业控制领域，分数阶PID控制器由于其结构简单、鲁棒性好、可靠性高等优点得到了广泛应用。分数阶PID的控制性能与控制器参数以及的优化整定直接相关。在工业控制过程中，多数控制对象是高阶、时

2、滞、非线性的，所以对分数阶PID控制器的参数整定是较为困难的。优化问题是工业设计中经常遇到的问题，许多问题最后都可以归结为优化问题。为了解决各种各样的优化问题，人们提出了许多优化算法，比较著名的有爬山法、神经算法和遗传算法等。优化问题有两个主要问题。一是要求寻找全局最小点，二是要求有较高的收敛速度。爬山法精度较高，但是易于陷入局部极小。遗传算法、神经网络算法等也还存在某些不足，前者要涉及到繁琐的编码解码过程和很大的计算量，后者的编程和解码过程需要大量CPU时间，算法易早熟，收敛易陷入局部最优，往往不能同时满足控制系统的速度和精度，且隐含层数目、神经元个数以及初始权值等参数选择都没有系统的方法。

3、1.2 基本的PID参数优化方法目前PID参数整定优化方法有很多，比如单纯形法、最速下降法、误差积分准则ISTE最优设定方法、遗传算法、蚁群算法等。单纯形法是一种求解多变量无约束最优化问题的直接搜索法，是求解非线性函数的无约束极值的一种经验方法；最速下降法是一种以梯度法为基础的多维无约束最优化问题的数值计算法，它的基本思想是选取目标函数的负梯度方法（最速下降方向）作为每步迭代的搜索方向，逐步逼近函数的极小值点；误差积分准则ISTE最优设定方法是针对一类特定被控对象的，如果被控对象形式已知，可以考虑使用这种ISTE误差积分准则作为目标函数进行参数优化；遗传算法借鉴了自然界优胜劣汰的进化思想，是模

4、拟达尔文生物进化论的自然选择和遗传学机理的生物进化过程的计算模型，通过模拟自然进化过程搜索最优解的方法。其基本思想是：先初始化一个种群（种群是由经过基因编码的一定数目的个体组成的，每个个体代表所求问题的一种解决方案），然后按照生物进化理论中的适者生存和优胜劣汰的原理，逐代演化产生出越来越好的个体。在每一代，根据个体的适应度大小挑选出较好的个体，并借助于自然遗传学的遗传算子进行组合交叉和变异，产生出代表新的解集的种群。经过数代的演化，将使得最终的种群更加适应环境，种群中的个体更加优质，把最后种群中的最优个体经过解码后作为问题的近似最优解；蚁群算法是受到自然界中真实蚁群集体行为的研究成果的启发而提

5、出的基于种群的模拟进化算法。蚂蚁从蚁巢出发寻找食物源，找到食物后在从食物源原路返回蚁巢的路上释放信息素，觅食的蚂蚁会跟随这个信息素踪迹找到食物源。信息素按照一定的比例释放的。路径越短，释放的信息素越多，浓度也越高；而信息素浓度越高，吸引的蚂蚁也越多；吸引的蚂蚁越多，遗留下的信息素也越多。最后所有的蚂蚁都集中到信息素浓度最高的一条路径上，这条路径就是从蚁巢到食物源的最短路径。为解决最优化问题人们提出过许多新技术和新方法，但工业和科学领域大量实际问题的困难程度正在日益增长，它们大多是根本无法在可接受的时间内找到解的问题。这类优化问题的困难性不仅体现在具有极大的规模，更为重要的是，它们多数是非线性的

6、、动态的、多峰的、具有欺骗性的或者不具有任何导数信息。因此，发展通用性更强、效率更高的优化算法总是需要的。1.3 常用的整定方法 这里列举在过程控制系统中常用的参数整定方法：经验法、衰减曲线法、临界比例度法、反应曲线法。用衰减曲线法整定调节器参数的方法是：在纯比例作用下，为为0，目的是要得到，衰减振荡过度过程曲线。根据所得曲线，若衰减大于应调整朝小比例带方向；若小于，应调整朝大比例带方向。记下的比例带，并在记录曲线上求得衰减时的调节周期，然后计算各值。临界比例度法考虑的实质是通过现场试验找到等幅振荡的过渡过程，得到临界比例度和等幅振荡周期。当操纵变量作阶跃变化时，被控变量随时间变化的曲线称为反

7、应曲线。对有自衡的非振荡过程，广义对象传递函数常可用近似。K，和T可用图解法等得出。调节器参数整定的反应曲线是依据广义对象的K，和T确定调节器参数的方法。在这些指标中，不同的系统有不同的侧重：强调快速跟踪的系统要求调节时间尽可能短些，强调稳定平稳的系统则要求超调量小，但基本上都要保证系统稳定收敛，即衰减比大于1，超调量必须在允许值的范围内，另外余差尽可能小至最后为零。影响控制系统指标的因素除了对象的时间常数、放大系数及滞后常数外，还有调节器的参数整定情况。调节器的参数整定是一个复杂的问题，这是因为这些参数的整定要考虑控制对象的各种特性，以及一些会影响系统运行过程的未知干扰；而且，调节器参数本身

8、的调整也会对系统的特性产生重大影响1-3。调节器的各参数对控制指标的具体影响主要体现在：比例带：越小，上升时间减小，衰减比S减小，稳定度下降。在工程上，比例带 常用比例度P来描述。微分作用：微分作用的大小由微分时间来决定。越大，越能克服系统的容量滞后和测量滞后，对缩短调节时间有一定作用。积分作用：积分作用通过积分时间来体现。越小，消除余差越快，稳定度下降，振荡频率变高。要实现PID参数的自整定，首先要对被控制的对象有一个了解，然后选择相应的参数计算方法完成控制器参数的设计。据此，可将PID参数自整定分成两大类：辨识法和规则法。基于辨识法的PID参数自整定，被控对象的特性通过对被控对象数学模型的

9、分析来得到，在对象数学模型的基础上用基于模型的一类整定法计算PID参数。基于规则的PID参数自整定，则是运用系统临界点信息或系统响应曲线上的一些特征值来表征对象特性，控制器参数由基于规则的整定法得到4。尽管当今出现了许多高级控制方法，但是实际控制系统仍然是以比例积分微分（PID） 控制为主，即使已经有了一些行之有效的整定规则，但是手动整定PID控制器参数仍是一件复杂和费时的工作。因此出现了许多自整定算法5。无论那种整定方法，都不是万能的，它们各有长处和不足，都有一定的适应范围。为了提高传统PID整定技术的适应能力，好多新的方法，如遗传算法，模糊逻辑控制等在最近几年里获得了很快的发展，并广泛地应

10、用于PID控制器参数整定中6。每种控制方法都有各自的优点以及适用范围，在实际的操作中不同的方法来实现同一控制模型，其精确度也会有差别。在工程实践中，总希望所选的方案是一切可能的方案中最优的方案，这就是最优控制的问题。解决最优控制的数学方法称为最优化方法，近几十年来，它已经是一门迅速发展的学科。在自动控制方面，将优化技术用于系统设计，能使设计出来的控制系统在满足一定的约束条件下，达到某种性能指标的函数为最小（或最大），这就是控制系统的最优化问题。1.4 本文的主要工作本文采用粒子群算法对分数阶PID参数进行寻优。先选择控制系统的目标函数，本控制系统选用时间乘以误差的绝对值，通过对控制系统的逐步仿

11、真，对结果进行分析。由于选取的这个目标函数的解析式不能直接写出，故采用逐步仿真来实现，然后采用工程上的整定方法（临界比例度法）粗略的确定其初始的五个参数以及积分和微分阶次，并以此进行寻优，得到较好的分数阶PID参数。再利用MATLAB编制粒子群算法寻优程序。通过粒子群算法优化系统性能最佳的分数阶PID参数后采用SIMULINK的仿真工具对分数阶PID参数优化系统进行仿真，得出系统的响应曲线。从中发现它的性能指标，都比原来的曲线有了很大的改进。第 2 章 粒子群算法的介绍2.1 粒子群算法思想的起源自然界中各种生物体均具有一定的群体行为，而人工生命的主要研究领域之一是探索自然界生物的群体行为，从

12、而在计算机上构建其群体模型。自然界中的鸟群和鱼群的群体行为一直是科学家的研究兴趣，生物学家Craig Reynolds在1987年提出了一个非常有影响的鸟群聚集模型7，在他的仿真中，每一个个体遵循：（1）避免与邻域个体相冲撞。（2）匹配邻域个体的速度。（3）飞向鸟群中心，且整个群体飞向目标。仿真中仅利用上面三条简单的规则，就可以非常接近的模拟出鸟群飞行的现象。1990年，生物学家Frank Heppner也提出了鸟类模型8，它的不同之处在于：鸟类被吸引飞到栖息地。在仿真中，一开始每一只鸟都没有特定的飞行目标，只是使用简单的规则确定自己的飞行方向和飞行速度（每一只鸟都试图留在鸟群中而又不相互碰撞

13、），当有一只鸟飞到栖息地时，它周围的鸟也会跟着飞向栖息地，这样，整个鸟群都会落在栖息地。1995年，美国社会心理学家James Kennedy和电气工程师Russell Eberhart共同提出了粒子群算法，其基本思想是受对鸟类群体行为进行建模与仿真的研究结果的启发。他们的模型和仿真算法主要对Frank Heppner的模型进行了修正，以使粒子飞向解空间并在最好解处降落。Kennedy在他的书中描述了粒子群算法思想的起源：自20世纪30年代以来，社会心理学的发展揭示：我们都是鱼群或鸟群聚集行为的遵循者。在人们的不断交互过程中，由于相互的影响和模仿，他们总会变得更相似，结果就形成了规范和文明。人

14、类的自然行为和鱼群及鸟群并不类似，而人类在高维认知空间中的思维轨迹却与之非常类似。思维背后的社会现象远比鱼群和鸟群聚集过程中的优美动作复杂的多：首先，思维发生在信念空间，其维数远远高于3；其次，当两种思想在认知空间会聚于同一点时，我们称其一致，而不是发生冲突。2.2 算法原理在一个D维的目标搜索空间中，有n个微粒组成一个粒子群，其中每个微粒是一个D维的向量，它的空间位置表示为xi =（xi1，xi2，xiD），i=1，2，n。微粒的空间位置是目标优化问题中的一个解，将它代入适应度函数可以计算出适应度值，根据适应度值的大小衡量微粒的优劣；第i个微粒的飞行速度也是一个D维的向量，记为vi=（vi1

15、，vi2，viD）；第i个微粒所经历过的具有最好适应值的位置称为个体历史最好位置，记为pi=（pi1，pi2，piD）；整个微粒群所经历过的最好位置称为全局历史最好位置，记为pg=（pg1，pg2，pgD），粒子群的进化方程可描述为： （2.1） （2.2）其中：下标j表示微粒的第j维，下标i表示微粒i，t表示第t代，c1，c2为加速常量，通常在（0,2）间取值，r1 U（0,1），r2 U（0,1）为两个相互独立的随机函数。从上述微粒进化方程可以看出，c1调节微粒飞向自身最好位置方向的步长，c2调节微粒向全局最好位置飞行的步长。通过分析基本粒子群的一些特点，可以知道式（2.1）中其第一部分为

16、微粒先前的速度；其第二部分为“认知”部分，表示微粒本身的思考；其第三部分为“社会”部分，表示微粒间的社会信息共享。目前，虽然模型的社会部分和认知部分的相对重要性还没有从理论上给出结论，但有一些研究已经表明对一些问题，模型的社会部分显得对认知部分更重要。2.3 算法流程基本粒子群算法的流程如下：（1）初始化粒子群，随机初始化各粒子。（2）根据适应度函数计算各粒子的适应度值。（3）对每个粒子，将它的适应度值与它的历史最优的适应度值比较，如果更好，则将其作为历史最优。（4）对每个粒子，比较它的适应度值和群体所经历的最好位置的适应度值，如果更好，则将其作为群最优。（5）根据方程（2.1）和方程（2.2

17、）对粒子的速度和位置进行进化。（6）如果达到结束条件（足够好的解或最大迭代次数），则结束，否则转步骤（2）。粒子群优化算法的流程如图2.1所示图2.1 基本粒子群算法流程图2.4 全局模型与局部模型在2.2描述的算法中，粒子的行为是受自身最优pbest和全局最优gbest的影响，这种版本称为全局版本PSO算法，如图2.2所示。另一种为局部版本PSO算法，在该算法中，粒子的行为是不受全局最优gbest影响的，而是受自身最优pbest和拓扑结构中邻近粒子中的局部最优lbest影响的，如图2.3所示。对局部版本，式（2.1）改为： （2.3） 其中，pij为邻近粒子的局部最优。比较两种版本的算法，我

18、们可以发现：因为全局版本PSO算法中所有粒子信息是共享的，所以算法收敛到全局最优的速度比局部版本PSO算法快。但全局PSO算法易陷入局部最优；局部PSO算法允许粒子与邻近粒子比较，相互施加影响，虽然算法收敛速度慢，但不易陷入局部最优。 图2.2 gbest模型 图2.3 lbest模型2.5 算法特点粒子群算法具有以下主要优点： 易于描述 设置参数少 容易实现 收敛速度快粒子群算法很容易实现，计算代价低且占用计算机硬件资源少。粒子群算法已被证明能很好地解决许多全局优化问题。当然，PSO算法也和其它全局优化算法一样，有易陷入局部最优，收敛精度不高，后期收敛速度慢等缺点。2.6 带惯性权重的粒子群

19、算法探索是偏离原来的寻优轨迹去寻找一个更好的解，探索能力是一个算法的全局搜索能力。开发是利用一个好的解，继续原来的寻优轨迹去搜索更好的解，它是算法的局部搜索能力。如何确定局部搜索能力和全局搜索能力的比例，对一个问题的求解过程很重要。1998年，Yuhui Shi9提出了带有惯性权重的改进粒子群算法。其进化过程为： （2.4） （2.5） 在式（2.1）中，第一部分表示粒子先前的速度，用于保证算法的全局收敛性能；第二部分、第三部分则是使算法具有局部收敛能力。可以看出，式（2.4）中惯性权重w表示在多大程度上保留原来的速度。w较大，全局收敛能力强，局部收敛能力弱；w较小，局部收敛能力强，全局收敛能

20、力弱。当w=1时，式（2.4）与式（2.1）完全一样，表明带惯性权重的粒子群算法是基本粒子群算法的扩展。实验结果表明，w在0.8,1.2之间时，PSO算法有更快的收敛速度，而当w1.2时，算法则易陷入局部极值。2.7 粒子群算法的研究现状在算法的理论研究方面。目前PSO算法还没有成熟的理论分析，少部分研究者对算法的收敛性进行了分析，大部分研究者在算法的结构和性能改善方面进行研究，包括参数分析，拓扑结构，粒子多样性保持，算法融合和性能比较等。PSO由于有简单、易于实现、设置参数少、无需梯度信息等特点，其在连续非线性优化问题和组合优化问题中都表现出良好的效果。第 3 章 用粒子群方法优化PID参数

21、分数阶PID控制是最早发展起来的控制策略之一，是指将偏差的比例（P）、积分（I）和微分（D）以及微积分的系数通过线性组合构成控制量，对被控对象进行控制。随着计算机的普及，分数阶PID控制在生产过程中已成为一种最常用的控制方法，在机电、冶金、机械、化工等诸多行业中获得了广泛的应用。3.1 分数阶PID控制原理图3.1给出分数阶PID控制系统的原理框图，该控制系统由模拟分数阶PID控制器和被控对象组成。图3.1 分数阶PID控制器原理图分数阶PID控制器的闭环控制系统如图所示，控制器主要是由分数阶PID和被控对象模型组成的反馈，闭环控制系统。根据给定的输入值Input与实际的输出值Output直间

22、的误差e（t），分数阶PID控制器通过对误差非线性控制，形成控制量u（t），对被控对象模型plant进行控制，以达到期望的输出目的。控制器的输入输出关系在时域中的表达：= （3.1）其中，为比例系数，为积分时间常数，为微分时间常数， 0为积分次数，0为微分次数。分数阶PID控制器中的各个校正环节的作用如下：1、比例环节：成比例地反映控制系统的偏差信号e（t），偏差一旦产生，控制器立即产生控制作用，以减小偏差；2、积分环节：主要用于消除静差，提高系统的无差度。积分作用的强弱取决于积分时间常数凡，凡越大，积分作用越强，反之则越弱；3、微分环节：反映偏差信号的变化趋势（变化速率），并能在误差信号变得

23、太大之前，在系统中引入一个有效的早期修正信号，从而加快系统的动作速度，减少调节时间。调节积分阶次和微分阶次使系统调节优化效果更佳明显。3.2 分数阶PID控制的特点分数阶PID控制器原理简单、鲁棒性好、可靠性高，因此一直是工业过程控制中应用最广的策略，尤其适用于可建立精确数学模型的确定性系统。在实际工业生产过程往往具有非线性、时变不确定性等困难性，难以建立精确的数学模型，应用常规PID控制器不能达到理想的控制效果。此外，在实际生产的现场中，常规PID控制器往往会受到参数整定过程繁杂的困扰，出现整定不良、性能欠佳的情况，对运行工况的适应性也很差。所以应用分数阶PID控制器。3.3 优化设计简介所

24、谓优化设计就是一种对问题寻优的过程，人们所从事的任何工作都希望尽可能做好，以期得到一个理想的目标。在日常的设计过程中，常常需要根据产品设计的要求，合理地确定各种参数，以达到最佳的设计目标。实际上，在任何一项设计工作中都包含着寻优过程，但这种寻优在很大程度上带有经验性，多根据人们的直觉、经验及不断试验而实现的，由于受到经验、时间、环境等条件的限制，往往难以得到最佳的结果。优化设计是20世纪60年代发展起来的一门新的学科，它是最优化技术和计算机技术在设计领域应用的结果。优化设计为工程设计提供了一种重要的科学设计方法，在解决复杂设计问题时，它能从众多的设计方案中找到尽可能完善的设计方案。要实现问题的

25、优化必须具备两个条件，一是存在一个优化目标；另一是具有多个方案可供选择。工程设计问题的最优化，可以表达为一组优选的设计参数，在满足一系列限制条件下，使设计指标达到最优。因而，优化设计的数学模型可由设计变量、目标函数和设计约束条件三部分组成。（1） 设计变量：在工程设计中，为区别不同的设计方案，通常是以被称为设计变量的不同参数来表示。（2） 目标函数：每一个设计问题，都有一个或多个设计中所追求的目标，它们可以用设计变量的函数来表示，被称为目标函。 （3） 设计约束：优化设计不仅要使所选择方案的设计指标达到最佳值，同时还必须满足一些附加的设计条件，这些附加设计条件都构成对设计变量取值的限制，在优化

26、设计中被称为设计约束。工程设计中的优化方法有多种类型，有不同的分类方法。若按设计变量数值的不同，可将优化设计分为单变量（一维）优化和多变量优化；若按约束条件的不同，可分为无约束优化和有约束优化；若按目标函数数量的不同，又有单目标优化和多目标优化10。3.4 目标函数选取在参数最优化的问题中要涉及性能指标函数，性能指标函数是被寻参数的函数，称为目标函数。选择不同的目标函数的出发点是使它即能比较明确的反映系统的品质，又便于计算。当然选择不同的目标函数，即使对于同一系统，寻优最后得到的优化参数也是会有所不同的。目标函数的选择分为两大类：第一类是特征型目标函数，它是按照系统的输出响应的特征提出的。第二

27、类是误差型目标函数，它是采用期望响应和实际响应之差的某个函数作为目标函数。这种目标函数实际上是对第一类目标函数的几个特征向量做数学分析，把它们包含在一个目标函数的表达式中。因此它反映整个系统的性能。几种常用的误差型目标函数：（1）误差平方的积分型。这种目标函数的表达式为 （3.4）其中e（t）=r（t）-y（t）表示系统误差。一般要求e（t）越小越好，即要求控制系统的输出响应y（t）尽可能的接近输入r（t）。由于在过度过程中e（t）时正时负，故取误差的平方进行积分。这种目标函数在数学上是很容易实现的，常常可以得到比较简单的解析式。但是在过度过程中，不同时期的误差是不完全相同的，如果全部用误差的

28、平方再积分显然是不怎么合理的，不能很好的反映系统的最终品质指标的要求。（2）时间乘以平方误差型。 （3.5）由于在误差平方上乘以了t，相当加上了时间权。这样过度过程的初始误差考虑比较少，而着重权衡过度过程中后期出现的误差。这种目标函数的选取不止一种方法可以更精确地反映系统的最终品质要求。（3）误差绝对值积分型。或者为 （3.6）其寻优方法显然要比其他两种方法优点突出。一方面加了绝对值，它克服了在过度过程中e（t）时正时负的缺点，另外加了时间t，这样过度过程中后期出现的误差也基本上能消除。因此本文在选择目标函数的表达式取。第 4 章 系统仿真研究4.1 MATLAB中SIMULINK仿真1.选用

29、函数为整定的函数。利用matlab中的simulink进行仿真模拟通过取不团的参数值获得不同的波形图，进而可以得到不同参数对系统的影响。然后通过粒子群优化算法得到最优的参数。图形如下：图4.1P模型控制器输出波形如下：图4.2为P模型输出波形图由图4.2可知，增大会增大系统的超调量，但降低了峰值时间和调节时间；而减少则相对地降低了超调量，但增加了峰值时间和调节时间。2.下图4.3给出了PI的模型图图4.3PI模型控制器输出的波形如下：图4.4为PI模型的输出波形图由图4.4可知，增大增加了系统了超调量；则相对地降低了系统超调量；无则系统存在余差。3下图4.5给出了PD控制模型图4.5PD模型图图4.6为PD模型的输出波形图由图4.6可知，增大会降低系统的超调量；减小会相对地增大系统超调量4.2 粒子群算法参数整定 用粒子群算法优化后得到的参数为=5.0017， =0.0229， =65.7811， 0.6，=1.2.。如下图给出整定前后输出的波形的：图4.7整定前的波形图图