粒子群优化算法分数阶PID控制器

1、目录 TOC o 1-5 h z HYPERLINK l bookmark8 o Current Document 摘要1 HYPERLINK l bookmark11 o Current Document 第1章绪论2 HYPERLINK l bookmark14 o Current Document 1.1研究背景和课题意义2 HYPERLINK l bookmark17 o Current Document 1.2基本的PID参数优化方法2 HYPERLINK l bookmark20 o Current Document 1.3常用的整定方法3 HYPERLINK l bookmark

2、23 o Current Document 1.4本文的主要工作5 HYPERLINK l bookmark26 o Current Document 第2章 粒子群算法的介绍6 HYPERLINK l bookmark29 o Current Document 2.1粒子群算法思想的起源6 HYPERLINK l bookmark33 o Current Document 2.2算法原理6 HYPERLINK l bookmark42 o Current Document 2.3算法流程7 HYPERLINK l bookmark51 o Current Document 2.4全局模型与局

3、部模型8 HYPERLINK l bookmark54 o Current Document 2.5算法特点9 HYPERLINK l bookmark57 o Current Document 2.6带惯性权重的粒子群算法9 HYPERLINK l bookmark60 o Current Document 2.7粒子群算法的研究现状10 HYPERLINK l bookmark63 o Current Document 第3章 用粒子群方法优化PID参数11 HYPERLINK l bookmark66 o Current Document PID控制原理11 HYPERLINK l bo

4、okmark72 o Current Document PID控制的特点12 HYPERLINK l bookmark75 o Current Document 3.3优化设计简介12 HYPERLINK l bookmark81 o Current Document 3.4目标函数选取13 HYPERLINK l bookmark87 o Current Document 第4章 系统仿真研究15 HYPERLINK l bookmark90 o Current Document MATLAB 中 simulink 仿真15 HYPERLINK l bookmark93 o Current

5、Document 粒子群算法参数整定18 HYPERLINK l bookmark96 o Current Document 4.3结果比较19结论20 HYPERLINK l bookmark103 o Current Document 致谢21 HYPERLINK l bookmark106 o Current Document 参考文献21 HYPERLINK l bookmark118 o Current Document 附录A （程序清单）22摘要本文主要研究基于粒子群算法控制系统PID参数优化设计方法以及对PID控 制的改进。PID参数的寻优方法有很多种，各种方法的都有各自的特点

6、，应按实 际的系统特点选择适当的方法。本文采用粒子群算法进行参数优化，主要做了如 下工作：其一，选择控制系统的目标函数，本控制系统选用时间乘以误差的绝对 值，通过对控制系统的逐步仿真，对结果进行分析。由于选取的这个目标函数的 解析式不能直接写出，故采用逐步仿真来实现；其二，本文先采用工程上的整定 方法（临界比例度法）粗略的确定其初始的三个参数疽K , Kd ,再利用粒子 群算法进行寻优，得到更好的PID参数；其三，采用SIMULINK的仿真工具对PID 参数优化系统进行仿真，得出系统的响应曲线。从中发现它的性能指标，都比原 来有了很大的改进。因此，采用粒子群算法的优越性是显而易见的。关键词目标

7、函数；PID参数；粒子群算法；优化设计；SIMULINK第1章绪论1.1研究背景和课题意义在现代工业控制领域，分数阶PID控制器由于其结构简单、鲁棒性好、可靠性高等优点 得到了广泛应用。分数阶PID的控制性能与控制器参数 %,K ,七 以及人,的优化整定直 接相关。在工业控制过程中，多数控制对象是高阶、时滞、非线性的，所以对分数阶PID控 制器的参数整定是较为困难的。优化问题是工业设计中经常遇到的问题，许多问题最后都可 以归结为优化问题。为了解决各种各样的优化问题，人们提出了许多优化算法，比较著名的 有爬山法、神经算法和遗传算法等。优化问题有两个主要问题。一是要求寻找全局最小点， 二是要求有较

8、高的收敛速度。爬山法精度较高，但是易于陷入局部极小。遗传算法、神经网 络算法等也还存在某些不足，前者要涉及到繁琐的编码解码过程和很大的计算量，后者的编 程和解码过程需要大量CPU时间，算法易早熟，收敛易陷入局部最优，往往不能同时满足控 制系统的速度和精度，且隐含层数目、神经元个数以及初始权值等参数选择都没有系统的方 法。1.2基本的PID参数优化方法目前PID参数整定优化方法有很多，比如单纯形法、最速下降法、误差积分准则STE最 优设定方法、遗传算法、蚁群算法等。单纯形法是一种求解多变量无约束最优化问题的直接 搜索法，是求解非线性函数的无约束极值的一种经验方法；最速下降法是一种以梯度法为基 础

9、的多维无约束最优化问题的数值计算法，它的基本思想是选取目标函数的负梯度方法（最 速下降方向）作为每步迭代的搜索方向，逐步逼近函数的极小值点；误差积分准则ISTE最优 设定方法是针对一类特定被控对象的，如果被控对象形式已知，可以考虑使用这种ISTE误 差积分准则作为目标函数进行参数优化；遗传算法借鉴了自然界优胜劣汰的进化思想，是模 拟达尔文生物进化论的自然选择和遗传学机理的生物进化过程的计算模型，通过模拟自然进 化过程搜索最优解的方法。其基本思想是：先初始化一个种群（种群是由经过基因编码的一 定数目的个体组成的，每个个体代表所求问题的一种解决方案），然后按照生物进化理论中 的适者生存和优胜劣汰的

10、原理，逐代演化产生出越来越好的个体。在每一代，根据个体的适 应度大小挑选出较好的个体，并借助于自然遗传学的遗传算子进行组合交叉和变异，产生出 代表新的解集的种群。经过数代的演化，将使得最终的种群更加适应环境，种群中的个体更 加优质，把最后种群中的最优个体经过解码后作为问题的近似最优解；蚁群算法是受到自然 界中真实蚁群集体行为的研究成果的启发而提出的基于种群的模拟进化算法。蚂蚁从蚁巢出 发寻找食物源，找到食物后在从食物源原路返回蚁巢的路上释放信息素，觅食的蚂蚁会跟随 这个信息素踪迹找到食物源。信息素按照一定的比例释放的。路径越短，释放的信息素越多， 浓度也越高；而信息素浓度越高，吸引的蚂蚁也越多

11、；吸引的蚂蚁越多，遗留下的信息素也 越多。最后所有的蚂蚁都集中到信息素浓度最高的一条路径上，这条路径就是从蚁巢到食物 源的最短路径。为解决最优化问题人们提出过许多新技术和新方法，但工业和科学领域大量 实际问题的困难程度正在日益增长，它们大多是根本无法在可接受的时间内找到解的问题。 这类优化问题的困难性不仅体现在具有极大的规模，更为重要的是，它们多数是非线性的、 动态的、多峰的、具有欺骗性的或者不具有任何导数信息。因此，发展通用性更强、效率更 高的优化算法总是需要的。1.3常用的整定方法这里列举在过程控制系统中常用的参数整定方法：经验法、衰减曲线法、临界比例度 法、反应曲线法。用衰减曲线法整定调

12、节器参数的方法是：在纯比例作用下，K为8，Kd 为0，目的是要得到4:1，衰减振荡过度过程曲线。根据所得曲线，若衰减大于4:1应调整 5朝小比例带方向；若小于4:1，应调整8朝大比例带方向。记下4:1的比例带5，并在记 录曲线上求得4:1衰减时的调节周期Tp，然后计算5，K , Kd各值。临界比例度法考虑的实质是通过现场试验找到等幅振荡的过渡过程，得到临界比例度和 等幅振荡周期。当操纵变量作阶跃变化时，被控变量随时间变化的曲线称为反应曲线。对有 自衡的非振荡过程，广义对象传递函数常可用GG)= Kexp (t s)件+1)近似。K，T和T 可用图解法等得出。调节器参数整定的反应曲线是依据广义对

13、象的K，T和T确定调节器参 数的方法。在这些指标中，不同的系统有不同的侧重：强调快速跟踪的系统要求调节时间尽可能短 些，强调稳定平稳的系统则要求超调量小，但基本上都要保证系统稳定收敛，即衰减比大于 1，超调量必须在允许值的范围内，另外余差尽可能小至最后为零。影响控制系统指标的因 素除了对象的时间常数、放大系数及滞后常数外，还有调节器的参数整定情况。调节器的参 数整定是一个复杂的问题，这是因为这些参数的整定要考虑控制对象的各种特性，以及一些 会影响系统运行过程的未知干扰；而且，调节器参数本身的调整也会对系统的特性产生重大 影响1-3。调节器的各参数对控制指标的具体影响主要体现在：比例带5 :比例

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

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

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

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

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

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

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

21、第i个微 粒所经历过的具有最好适应值的位置称为个体历史最好位置，记为pi=（pi1，pi2,，PiD）； 整个微粒群所经历过的最好位置称为全局历史最好位置，记Xpg=（Pg1，Pg2,，PgD），粒子 群的进化方程可描述为：(2.1)(2.2)v G +1)= v Q+c, (t 舄 Qf Q)+c r (t 舄 Qf Q)ii1 1 ii2 2 gji七()+ % (t+1)其中：下标j表示微粒的箫维，下标i表示微粒i，t表示第t代，C，c2为加速常量，通常 在（0,2）间取值，U（0,1），r2U（0,1）为两个相互独立的随机函数。从上述微粒进化方程可 以看出，调节微粒飞向自身最好位置方向

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

23、较它的适应度值和群体所经历的最好位置的适应度值，如果更好， 则将其作为群最优。（5）根据方程（2.1）和方程（2.2）对粒子的速度和位置进行进化。（6）如果达到结束条件（足够好的解或最大迭代次数），则结束，否则转步骤（2）。粒子群优化算法的流程如图2.1所示图2.1基本粒子群算法流程图2.4全局模型与局部模型在2.2描述的算法中，粒子的行为是受自身最pbest和全局最优gbe st的影响，这种版本称为全局版本PSO算法，如图2.2所示。另一种为局部版本PSO算法，在该算法中，粒子的行 为是不受全局最优gbe st影响的，而是受自身最优pbest和拓扑结构中邻近粒子中的局部最优 lbest影响的

24、，如图2.3所示。对局部版本，式(2.1)改为：v (t +1)= v Q+c,(X Qr Q),c r(X Qr Q) (2.3)ijij 1 1 ij ij2 2 ij ij其中，Pjj为邻近粒子的局部最优。比较两种版本的算法，我们可以发现：因为全局版本PSO算法中所有粒子信息是共享的， 所以算法收敛到全局最优的速度比局部版本?$0算法快。但全局PS O算法易陷入局部最优;局部PSO算法允许粒子与邻近粒子比较，相互施加影响，虽然算法收敛速度慢，但不易陷入局部最优。2.5算法特点粒子群算法具有以下主要优点：易于描述设置参数少容易实现收敛速度快粒子群算法很容易实现，计算代价低且占用计算机硬件资

25、源少。粒子群算法已被证明能 很好地解决许多全局优化问题。当然，PSO算法也和其它全局优化算法一样，有易陷入局部 最优，收敛精度不高，后期收敛速度慢等缺点。2.6带惯性权重的粒子群算法探索是偏离原来的寻优轨迹去寻找一个更好的解，探索能力是一个算法的全局搜索能 力。开发是利用一个好的解，继续原来的寻优轨迹去搜索更好的解，它是算法的局部搜索能 力。如何确定局部搜索能力和全局搜索能力的比例，对一个问题的求解过程很重要。1998(2.4)(2.5)年，Yuhui Shi9提出了带有惯性权重的改进粒子群算法。其进化过程为: v (t +1)二心 Q+c, (t 禹(t)- x (t )+c, (t 禹(t

26、)- x (t) jj1 1 jj2 2 gji七(t+1)二七(t)+% (t+1)在式(2.1)中，第一部分表示粒子先前的速度，用于保证算法的全局收敛性能；第二部分、 第三部分则是使算法具有局部收敛能力。可以看出，式：2.4)中惯性权重可表示在多大程度上 保留原来的速度。w较大，全局收敛能力强，局部收敛能力弱；w较小，局部收敛能力强， 全局收敛能力弱。当w=1时，式(2.4)与式(2.1)完全一样，表明带惯性权重的粒子群算法是基本粒子群算法 的扩展。实验结果表明，w在0.8,1.2之间时，PSO算法有更快的收敛速度，而当w1.2时， 算法则易陷入局部极值。2.7粒子群算法的研究现状在算法的

27、理论研究方面。目前?$。算法还没有成熟的理论分析，少部分研究者对算法的 收敛性进行了分析，大部分研究者在算法的结构和性能改善方面进行研究，包括参数分析， 拓扑结构，粒子多样性保持，算法融合和性能比较等PS O由于有简单、易于实现、设置参 数少、无需梯度信息等特点，其在连续非线性优化问题和组合优化问题中都表现出良好的效 果。第3章用粒子群方法优化PID参数分数阶PID控制是最早发展起来的控制策略之一，是指将偏差的比例(P)、积分(I)和微分 (D)以及微积分的系数人,日通过线性组合构成控制量，对被控对象进行控制。随着计算机的 普及，分数阶PID控制在生产过程中已成为一种最常用的控制方法，在机电、

28、冶金、机械、 化工等诸多行业中获得了广泛的应用。3.1分数阶PID控制原理图3.1给出分数阶PID控制系统的原理框图，该控制系统由模拟分数阶PID控制器和被控 对象组成。图3.1分数阶PID控制器原理图分数阶PID控制器的闭环控制系统如图所示，控制器主要是由分数阶PID和被控对象模型组成 的反馈，闭环控制系统。根据给定的输入值Input与实际的输出值Output直间的误差e(t)， 分数阶PID控制器通过对误差非线性控制，形成控制量口(t)，对被控对象模型plant进行控 制，以达到期望的输出目的。控制器的输入输出关系在时域中的表达：u(t) = k e(t) + K D2(t) + K Du

29、e(t)(3.1)其中，Kp为比例系数，K为积分时间常数，Kd为微分时间常数，人0为积分次数， 0 为微分次数。分数阶PID控制器中的各个校正环节的作用如下：1、比例环节：成比例地反映控制系统的偏差信号e(t)，偏差一旦产生，控制器立即产生 控制作用，以减小偏差；2、积分环节：主要用于消除静差，提高系统的无差度。积分作用的强弱取决于积分时 间常数凡，凡越大，积分作用越强，反之则越弱；3、微分环节：反映偏差信号的变化趋势(变化速率)，并能在误差信号变得太大之前， 在系统中引入一个有效的早期修正信号，从而加快系统的动作速度，减少调节时间。调节积分阶次人和微分阶次使系统调节优化效果更佳明显。3.2分

30、数阶PID控制的特点分数阶PID控制器原理简单、鲁棒性好、可靠性高，因此一直是工业过程控制中应用最 广的策略，尤其适用于可建立精确数学模型的确定性系统。在实际工业生产过程往往具有非线性、时变不确定性等困难性，难以建立精确的数学模 型，应用常规PID控制器不能达到理想的控制效果。此外，在实际生产的现场中，常规PID 控制器往往会受到参数整定过程繁杂的困扰，出现整定不良、性能欠佳的情况，对运行工况 的适应性也很差。所以应用分数阶PID控制器。3.3优化设计简介所谓优化设计就是一种对问题寻优的过程，人们所从事的任何工作都希望尽 可能做好，以期得到一个理想的目标。在日常的设计过程中，常常需要根据产品

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

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

33、变量（一维）优化和多变量优化；若按约束条件的不同，可分为无约束优 化和有约束优化；若按目标函数数量的不同，又有单目标优化和多目标优化103.4目标函数选取在参数最优化的问题中要涉及性能指标函数，性能指标函数是被寻参数的函 数，称为目标函数。选择不同的目标函数的出发点是使它即能比较明确的反映系 统的品质，又便于计算。当然选择不同的目标函数，即使对于同一系统，寻优最 后得到的优化参数也是会有所不同的。目标函数的选择分为两大类：第一类是特征型目标函数，它是按照系统的输 出响应的特征提出的。第二类是误差型目标函数，它是采用期望响应和实际响应 之差的某个函数作为目标函数。这种目标函数实际上是对第一类目标

34、函数的几个 特征向量做数学分析，把它们包含在一个目标函数的表达式中。因此它反映整个 系统的性能。几种常用的误差型目标函数：（1）误差平方的积分型。这种目标函数的表达式为J = j e2(t )dt(3.4)0其中e(t)=r(t)-y(t)表示系统误差。一般要求e(t )越小越好，即要求控制系统的 输出响应y(t)尽可能的接近输入r(t)。由于在过度过程中e(t)时正时负，故取误差 的平方进行积分。这种目标函数在数学上是很容易实现的，常常可以得到比较简 单的解析式。但是在过度过程中，不同时期的误差是不完全相同的，如果全部用 误差的平方再积分显然是不怎么合理的，不能很好的反映系统的最终品质指标的

35、 要求。时间乘以平方误差型。这种目标函数的表达式为J = jt te dt(3.5)0由于在误差平方上乘以了 t，相当加上了时间权。这样过度过程的初始误差 考虑比较少，而着重权衡过度过程中后期出现的误差。这种目标函数的选取不止 一种方法可以更精确地反映系统的最终品质要求。误差绝对值积分型。这种目标函数的表达式为J牛(t )也 或者为 J = jtt|e(t)dt(3.6)00其寻优方法显然要比其他两种方法优点突出。一方面加了绝对值，它克服了在过度过程 中e (t)时正时负的缺点，另外加了时间t，这样过度过程中后期出现的误差也基本上能消除。 因此本文在选择目标函数的表达式取J =tt|e(tdt

36、。0第4章系统仿真研究4.1 MATLAB 中SIMULINK 仿真1.选用函数为整定的函数。利用matlab中的simulink进行仿真模拟通过取不2s 2 + 3s +1团的参数值获得不同的波形图，进而可以得到不同参数对系统的影响。然后通过粒子群优化算法得到最优的参数。图形如下：itapl图4.1P模型控制器输出波形如下:图4.2为P模型输出波形图由图4.2可知，增大七会增大系统的超调量，但降低了峰值时间和调节时间；而减少K则相对地降低了超调量，但增加了峰值时间和调节时间。图4.3PI模型控制器输出的波形如下：图4.4为PI模型的输出波形图由图4.4可知，增大K增加了系统了超调量；而减少K

37、.则相对地降低了系统超调量;无K则系统存在余差。i图4.5PD模型图输出的波形如下:图4.6为PD模型的输出波形图由图4.6可知，增大K会降低系统的超调量；减小K会相对地增大系统超调量4.2粒子群算法参数整定用粒子群算法优化后得到的参数为七=5.。17,K =.229,如，0.6,R =1.2.。如下图给出整定前后输出的波形的:与1=_图4.7整定前的波形图-dluq图4.8整定后的波形图4.3结果比较粗略的确定其初始的参数气,K ,Kd以及人,，并通过仿真得到响应曲线，曲线的 上升时间虽然比较快，但是过度时间比较长，超调量也过大，这对工程实践是不利的。再采 用粒子群算法整定参数，通过粒子群算

38、法对参数优化后的曲线，曲线的各个指标也有了明显 的提高，尤其是超调量有了明显的减少，上升时间也有了明显的缩短。这正是我们所期待的。 通过上图可以说明粒子群算法可以有效的优化分数阶PID模型。结论本设计采用多变量寻优的粒子群算法对控制系统的PID参数进行优化设计，通过采用一般 的整定方法粗略的确定其初始的参数七,K , % 以及人,，并采用粒子群算法用SIMULINK的仿真工具对PID参数进行优化，得出系统的响应曲线。从中发现它的性能指标， 都比原来的曲线有了很大的改进。曲线的各个指标也有了明显的提高，尤其是超调量有了明 显的减少，上升时间也有了明显的缩短。因此，采用粒子群算法的优越性是显而易见

39、的。除此之外，本设计还研究了?、I、D各参数对系统的影响以及K 对该系统的作用，叙 述如下：增大Kp会增大系统的超调量，但降低了峰值时间和调节时间；而减少Kp则相对地 降低了超调量，但增加了峰值时间和调节时间。增大K增加了系统了超调量；而减少K.则相对地降低了系统超调量；无K.则系统 存在余差。增大K会降低系统的超调量；减小K会相对地增大系统超调量。加了屋后，系统的超调量有了明显进步，稳定时间也相对降低，有效的改善了控 制品质。致谢本次毕业设计论文全程都是在黄宜庆老师一个人的指导下完成的。黄老师治学严谨，学 识渊博，为人和蔼，诲人不倦。在我写论文陷入山穷水尽时，给予及时的帮助和指导，并在 关键

40、处给予细致的讲解，使我时有柳暗花明又一村的感觉，同时及时指正了我论文中粗心马 虎的缺点，使我在做论文的同时又养成了严谨的治学态度，这对我以后的学业和工作都是有 很大的帮助的。参考文献邱黎辉，等.模糊PID控制在中央空调系统中的应用J.计算机测量与控制，2004, 12(1)： 15-26.Ho Ming-Jzu, Lin Chia-Yin. PID controller design for robust performanceJ. IEEE Transactions on Automatic Control, 2003, 48(8)： 3-8.Ho M.T.Robust and non-fragile PID controller design J. Robust