优秀毕业设计毕业论文实验锅炉炉温自校正PID控制系统设计

1、毕业设计（论文）题目：实验锅炉炉温自校正pid控制系统设计本文以递推最小二乘法为自适应规律，研究实验锅炉温度控制，用matlab语言编 程并进行系统仿真，在仿真结果的基础上进行分析研究.结果表明当采用自校止pid算 法时，系统有自适应能力，能根据被调节系统自动调节kp、ki、kd参数使系统达到稳定, 解决了长期以来大时滞实验锅炉系统p1d参数设定难的问题.为了使系统更精确，本文 采用了 “带遗忘因子的递推最小二乘法”估计算法.关键字：自校正pid;炉温控制；参数估计；最小二乘法abstractthis paper adopts recursive least square method to

2、research how to control the experimental boiler temperature and uses matlab to program and simulate the further research and analysis are made on the basis of system simulation results, the results show that when adopting the self-revised pid controller, the system can adapt the complex working cond

3、itions and the controller can select the pid parameters automatically morever, it has solved the difficult problem of setting large delay experimental boiler system pid parameters. in order to make the system more accurate, this paper adopts the estimation algorithm of recursive least-square method

4、with forgetting factor.keywords: self-revised pid controller; temperature control; recursive least square method; parameter estimation摘要abstract1绪论11选题背景及意义11.2论文的主要内容11.3本人主要工作12常规pid控制算法及其改进算法22常规pid控制原理22.1.1常规pid调节器算法22.1.2常规pid调节器的参数整定32.1.3常规pid调节器在实际应用中的局限 32.2数字pi d控制32.2.1位置式pid控制算法42.2.2增量

5、式pid控制算法52.2.3数字pid控制器的参数整定方法52.2.4采样周期的选择63自校正pid控制算法73自适应控制系统原理73.1.1 概述73.1.2模型参考自适应控制83.1.3自校正控制83.2自校正控制系统93.3.1递推最小二乘估计103.3.2带遗忘因子的递推最小二乘算法的递推算式123.3.3初值的确定133.4本论文所用自校正pid控制算法133.4.1具体框图和原理133.4.3带遗忘因子的递推最小二乘法174系统硬件的结构设计184系统硕件的结构184.2自校正pid实验锅炉控制系统原理184.3自校正控制算法设计204.4自校正pid算法设计流程图5 matlab

6、仿真及结果分析226结论26参考文献27致谢错误！未定义书签。附录281绪论1.1选题背景及意义锅炉炉温控制是典型的工业过程控制对象，例如：在冶金工业、化工生产、 电力工程、造纸行业、机械制造和食品加工等诸多领域中，人们都需要对各类加 热炉、热处理炉、反应炉和锅炉中的温度进行检测和控制。通常，电阻炉炉温控 制都采用偏差控制法。偏差控制的原理是先求岀实测炉温对所需炉温的偏差值， 然后对偏差值处理获得控制信号去调节电阻炉的加热功率，以实现对炉温的控 制。通常对偏差进行比例、积分和微分控制又称pid控制，是工业控制过程中 应 用最广泛的一种控制形式。但由于锅炉温度控制具有升温单向性，大惯性，纯滞 后

7、性等特点，很难用数学建立精确的模型和确定参数。当被控对象的结构和参数不能完全掌握，或得不到精确的数学模型时，控制 理论的其它技术难以采用时，系统控制器的结构和参数必须依靠经验和现场调试 来确定，多年来，人们一直在努力寻找一种方法，能够简单易行的在线实时自动 整定pid参数。自校正pid调节器就是其中的一种，自校正pid是自校正思想 和常规pid控制器相结合的产物，它吸收了两者的优点，设计参数少，能够在 线整定和校正pid控制器参数，具有较强的适应能力。1.2论文的主要内容本论文详细的阐述了常规pid控制优点和缺陷控制功能，以及如何把最小 二乘法的思维与常规pid相结合的，得出了基于最小二乘法的

8、自校正pid控制 算法。以及如何把自校正pid应用到实验锅炉稳定控制系统中，并对自校正pid 实验锅炉温度控制系统进行了 matlab仿真，并与常规pid控制的实验锅炉炉温 控制系统相比较，得岀了自校正pid能自适应实验锅炉系统且自动调节pid参 数以达到最优控制，解决了常规pid参数难以设定的问题。1.3本人主要工作(1) 对自适应控制技术和常规pid调节器的原理进行研究，熟悉其模型原理;(2) 确定整个实验锅炉自校正pid控制系统算法；(3) 设计实验锅炉对象模型组成框图和结构图； 对实验锅炉模型实施自校正pid控制，进行matlab仿真，得到多个设定 值的仿真波形，并与其他控制方法比较得

9、出结论。2常规pid控制算法及其改进算法pid控制是模拟控制系统最常用的控制规律之一。模拟pid控制系统原理框图如图2所示，系统由模拟pid控制器和被控对象组成。r(t)+y(t)图2常规pid控制原理图2.1常规pid控制原理2.1.1常规pid调节器算法pid控制器是一种线性控制器，它根据给定值厂与实际输出值y(f)构成控 制偏差：(2-1)e(r) = r(r)-y(r)(2-2)将偏差的比例p、积分i和微分d通过线性组合构成控制量，对受控对象进 行控制。其控制规律为：u(t) = kp + e(t)d 杆1 f j°其传递函数为:=1 + 丄 +(2-3)e(s) 5 tts

10、 )式屮，kp比例系数，tt积分时间常数，td微分时间常数。pid控制器各校正环节的作用如下：(1) 比例环节：成比例地反映控制系统的偏差信号£()，偏差一旦产生，控 制器立即产生控制作用，以减小偏差。(2) 积分环节：主要用于消除误差，提高系统的精度。积分作用的强弱取决 于积分时间常数7； , 7；越大，积分作用越弱，反z则越强。(3) 微分环节：反映偏差信号的变化趋势(变化速率)，并能在偏差信号变得 太大之前，在系统中引入一个有效的早期修正信号，从而加快系统的动作速度， 减少调节时间。2.1.2常规pid调节器的参数整定pid控制器的参数整定是控制系统设计的核心内容。它是根据被控

11、过程的特 性确定pid控制器的比例系数、积分i寸间和微分i寸间的大小。pid控制器参数整定的方法很多，概括起来有两人类：一是理论计算整定法。 它主要是依据系统的数学模型，经过理论计算确定控制器参数。这种方法所得到 的计算数据未必可以直接用，还必须通过工程实际进行调整和修改。二是工程整 定方法，它主要依赖工程经验，直接在控制系统的试验中进行，并ii方法简单、 易于掌握，在工程实际中被广泛采用。pid控制器参数的工程整定方法，主要有 稳定边界法、动态特性参数法和衰减曲线法。三种方法各有其特点，其共同点都 是通过试验，然后按照工程经验公式对控制器参数进行整定。但无论采用哪-种 方法所得到的控制器参数

12、，都需要在实际运行中进行最后调整与完善。现在-般 采用的是稳定边界法。利用该方法进行pid控制器参数的整定步骤如下i(1) 首先预选择一个足够短的采样周期让系统工作；(2) 加入比例控制环节，直到系统对输入的阶跃响应出现临界振荡，记下这 吋的比例放大系数和临界振荡周期；(3) 在一定的控制度下通过公式讣算得到pid控制器的参数。2.1.3常规pid调节器在实际应用中的局限(1) 由于实现控制系统的器件的物理特性的限制，使得pid控制器获取得原 始信息偏离真实值，而其产生的控制作用偏离理论值。(2) 由于在系统的设计与整定过程中，要兼顾动态与稳定性能，只能采取折 衷方案，难以大幅度提高控制系统的

13、性能指标。(3) 对于存在强非线性、快速时变不确定性、强干扰等特性的对象，控制效 果较差。2.2数字pi d控制在连续时间控制系统中，pid控制器应用得非常广泛。其设计技术成熟，长 期以来形成了典型的结构，参数整定方便，结构更改灵活，能满足一般的控制要 求。计算机控制是一种采样控制，它只能根据采样时刻的偏差值计算控制量。因此连续pid控制算法不能直接使用，需要采用离散化方法。在计算机pid控制 中，使用的是数字pid控制器。数字pid控制比连续pid控制更为优越，因为 计算机程序的灵活性，很容易克服连续pid控制中存在的问题，经修正得到更 完善的数字pid算法。2.2.1位置式pid控制算法按

14、模拟pid控制算法，以一系列的采样吋刻点代表连续时间以矩形 法数值积分近似代替积分，以一阶向后差分近似代替微分，即：t = kt(2-4)(2-5)kk；=o&)幺伙)£伙一1) dt t"；=0可得离散pid表达式:“(/) = kp e(k)- e(j) + -e(k)-e伙一 1)k pe(k) + ( $>(力 + kde(k)- e(k 1)式中，k严kp*,k°=心半,t为采样周期，k = 1,2,粼)和幺伙-1)分别为第(kl)和第k时刻所得的偏差信号。根据z变换迭值定理和滞后定理可得：r1z工 </)= e(z),ze 伙1)

15、= ce ;=0 i?则对离散pid的表达式(25)作z变换，可得"=kpeq + b+ 心比广1(初(2-6)1-1由(26)可得数字pid控制器的z传递函数为：g(c =络=kp + 吾才 + 心(1 厂)(2-7)e(s)1-z算法中，为了求和，必须将系统偏差的全部过去值心)。=1,2,灯都存储 起来。这种算法得出控制量的全量输出讥灯(的是控制量的绝对数值,ug的值 和执行机构的位置是一一对应的，因此称其为位置式pid控制算法。2.2 .2增量式pid控制算法位置式pid算法计算时需要对/伙)进行累加，计算机运算工作量很大。而 口，因为计算机输出的”对应的是执行机构的实际位置，

16、如果计算机出现故障, u(k)的大幅度变化会引起执行机构位置的大幅度变化，这就可能造成重人的生产 事故。并u有些执行机构需要的不是控制量的绝对值，而是控制量的增量(例如 去驱动步进电动机)时，需要的就是pid的增量算法。所谓增量式pid是指数 字控制器的输出只是控制量的增量皿伙)。根据位置式算法的公式(25河知，&一1u(k -1) = kp(k - v) + k工e(j) + k卜伙一 1) 一幺伙 一2 j=o则有：(2-8)(2-9)au(k) = u(k) - u(k -1)=kpw(k) -e(k -1)1 + k泮伙)+ kde(k) + e伙 + 2 )式(28)称为增量

17、式pid算法，对式(28)进行归并后，得:伙)=q()e(k) + qe(k -1) + q2e(k 一 2其屮,q° = k p + k + k d q=-kp-2kd% =心式(2-7)n经看不出是pid的表达式了，也看不出p、i、d作用的直接关系，只 表示了各次误差量对控制作用的影响。由式(29)看出，数字增量pid算法，只 要存储最近的三个误差采样值粼)，讹-1),锹-2)就足够了。由于增量式控制算 法不需要累加，控制增量伙)仅与最近3次釆样有关，所以误动作时影响小， 而且较容易通过加权处理获得比较好的控制效果。2.2.3数字pid控制器的参数整定方法数字pid控制器控制参数

18、的选择，可按连续一时间pid参数整定方法进行。 在选择数字pid参数之前，首先应该确定控制器的结构。对允许稳态误差的系 统，可以适当选择p或者pd控制器，使稳态误差在允许的范围内。对于必须消 除稳态误差系统，应选择包含积分控制的pi或者pid控制器。一般来说，pi、 pi和pd控制器应用较多。对于有滞后的对象，往往都加入微分控制。控制器结 构确定后，即可开始选择参数。参数的选择，是要根据受控对象的具体特性和对 控制系统的性能要求进行。工程上，一般要求整个闭环系统是稳定的，对给定量 的变化能迅速响应并平滑跟踪，超调量小；在不同干扰作用下，能保证被控量在 给定值；当环境参数发生变化时，整个系统能保

19、持稳定等等。这些要求，对于控 制系统自身性能来说，有些是卯盾的。我们必须满足主要方面的要求，兼顾其他 方而，适当地折中处理。pid控制器参数的整定，可以不依赖于受控对彖的数学 模型。工程上，pid控制器的参数常常是通过实验来确定，通过试凑，或者通过 实验经验公式来确定。2.2.4采样周期的选择在采样数据控制系统中，如果设采样周期为门则采样数率为1/t,采样角 频率co = 2tt/t采样周期丁是设计者要精心选择的重要参数，系统的性能与采样 周期的选择有密切关系。采样周期的选择受多方面因素的影响，主要考虑的因素 分析如下：(1) 香农(shannon)采样定理71t< (輕被采样信号的上限

20、角频率)%给岀了采样周期的上限。满足这一定理，采样信号方可模拟或者近似地恢复为原 模拟信号，而不丢失主要信息。在这个限制范围内，采样周期越小，采样一数据 控制系统的性能越接近于连续一时间控制系统。(2) 闭环系统对于给定信号的跟踪，要求采样周期要小。(3) 从抑制扰动的要求来说，采样周期应该选择得小一些。(4) 从执行元件的要求来看，有时要求输入控制信号要保持一定的宽度。(5) 从计算机精度考虑，采样周期不宜过短。(6) 从系统成本上考虑，希望采样周期越长越好。综合上述因素，选择采样周期，应在满足控制系统的性能耍求的条件下，尽 可能地选择低的采样速率。可以根据所设计的系统的具体情况，采用试凑的

21、方法, 在试凑过程中根据各种合理的建议来预选采样周期，多次试凑，选择性能较好地 一个作为最后的采样周期。3自校正pid控制算法3.1自适应控制系统原理3.1.1概述在许多工程中，被控对象和过程的数学模型事先难以确定，即使在某-条件 下已被确定了的数学模型，在条件改变后，其动态参数或其模型结构仍经常发生 变化。因而，常规调节器不能得到好的品质，这就需要设计一种特殊的控制系统, 它能够自动地补偿各方面的非预知变化，“自适应控制系统”也就应运而生。自 适应控制系统(adaptive control systems)是现代控制理论的一个重要分支。所谓 自适应控制系统是指能在系统和环境的信息不完备的情况

22、下改变自身特性来保 持良好工作品质的控制系统，又称适应控制系统。由于自适应控制的对象是那 些存在不定性的系统，所以这种控制应首先能在控制系统的运行过程中，通过不 断地量测系统的输入、状态、输出或性能参数，逐渐了解和掌握对象。然后根 据所得的过程信息，按一定的设计方法，做出控制决策去更新控制器的结构、参 数或控制作用，以便在某种意义下使控制效果达到最优或次最优，或达到某个预 期目标。按此设计思想建立的控制系统便是自适应控制系统。由此可见，一个 自适应控制系统必然具有下列三个基本特征：(1) 过程信息的在线积累在线积累过程信息的目的，是为了降低受控系统原有的不定性。为此可用系 统辨识的方法在线辨识

23、受控系统的结构和参数，直接积累过程信息；也可通过量 测能反映过程状态的某些辅助变量，间接积累过程信息。在系统辨识中，结构辨 识比参数估计困难得多。(2) 可调控制器可调控制器是指它的结构、参数或信号可以根据性能指标要求进行自动调 整。这种可调性要求是由受控系统的不定性决定，否则就无法对过程实现有效控 制。(3) 性能指标的控制性能指标的控制可分为开环控制方式和闭环控制方式两种。若与过程动态相 关联的辅助变量可测，而口此辅助变量与可调控制器参数之间的关系又可根据物 理学的知识和经验导出，这吋就可通过此辅助变量直接调整可调控制器，以期达 到预定的性能指标。这就是性能指标的开环控制，它的特点是没有根

24、据系统实际 达到的性能指标再作进一步的调整。与开环控制方式不同，在性能指标的闭环控 制中，还要获取实际性能与预定性能之间的偏差信息，将其反馈后修改可调控制 器，真到实际性能达到或接近预定性能为止。3.1.2模型参考自适应控制模型参考自适应控制系统(model reference adaptive system),简称mras, 由以下几个部分组成：参考模型、被控对象、反馈控制器和调整控制参数的自适 应机构等组成。在这个系统中，把参考输入r同吋加到参考模型和可调系统的输入端，用减 法器将参考模型和可调系统的输岀直接相减，得到广义误差信号e。自适应机构 按一定准则利用广义误差信号来修改可调系统的调

25、节器参数，或产生一个辅助输 入信号，使广义误差的某个指标达到极小。当可调系统特性与参考模型特性逐渐逼近, 广义误差趋于极小或下降为零，调节过程结束。当对象特性在运行中发生了变化，控制器参数的自适应调整过程与上述过程 完全一样。设计这类自适应控制系统的核心问题是如何综合自适应调整律，即自适应机 构所应遵循的算法。关于自适应律的设计目前存在两种不同的算法，-种称为局 部参数最优化的方法，即利用梯度或者其他参数优化的递推算法，求得一组控制 器的参数，使得某个预定的性能指标达到最小。这种方法的缺点是不能保证参数 调整过程中，系统总是稳定的。自适应律的另一种设计方法是基于稳定性理论的 方法，其基本思想是

26、保证控制器参数自适应调节过程是稳定的，然后再尽量使这 个过程收敛快一些。这类系统的难点在于系统稳定性分析，李雅普诺夫稳定性和 波波夫超稳定性理论都是设计自适应律的有效工具。3.1 .3自校正控制自校正控制是70年代发展起来的一种随机自适应控制，产生背景是：工业过 程控制由于强随机干扰、模型未知、参数时变、大时滞等因素，导致常规的控制 方法效果差。这类自适应控制系统的一个主要特点就是具有一个被控对象数学模型的在线 辨识环节，具体地说就是加入了一个对象参数的递推估计器。它根据系统的运行, 首先对被控对象进行在线辨识，然后再根据辨识出的模型参数和事先指定的性能 指标，在线地综合控制作用。当对象参数己

27、知时，对调节器的参数进行在线求解。 由于调节器的控制律是多样的，参数估计的方法也是多样的，因此自校正调节器 非常灵活，采用各种不同的控制方法和估计方法来搭配，能满足不同的性能要求, 其难点在于收敛性。本文主要采用自校正控制的pid调节器，所采用的控制律是pid调节，参 数估计方法是最小二乘法。3.2自校正控制系统图3.1 口校正控制系统结构图自校正控制系统结构图如图3所示：系统由三个部分组成：参数估计器、 控制器参数计算、调节器。其中，参数估计器又称辨识器，用来根据测量得到的 系统的输入输出信息，不断地在线辨识系统的模型结构和参数;控制器参数计算 是指按照一定的控制算法，由辨识得到的参数，取得

28、响应控制算法的参数;调节 器则是形成控制量，送往被控过程或对象，达到控制目标。自校正控制系统有两个反馈回路，外环调节器以对象输入和输出的反馈为基 础，而内环调节器具有外环调整的可调参数。参数估计和控制设计必须在线实现, 因此参数估计必须采用递推算法，控制器设计必须采用计算尽量简单的设计方 法。该系统能自动校正自己的参数以得到希望的闭环系统特性。图3所示的自校正调节器的参数是经由参数估计和控制的设计计算而间接 进行更新的。但也可以将对象重新参数化，即利用调节器的参数来表示模型。这 时，就不需要进行设计计算这个环节。算法将大大简化，设计机构的方框将不复 存在，调节器参数将直接更新。本文研究的白校正

29、pid调节器就是直接更新参 数。自校正控制系统结构形式通常分为以下两种：(1) 显式结构。如图3所示，首先估计出控制对象的参数，然后进行控制 器设计，这样的机构称为显式结构。(2) 隐式结构。将控制器对象的参数估计和控制器设计这两个步骤结合在一 起，直接估计出控制器参数，从而大大简化白校正控制的算法，这样的结构形式 称为隐式结构。自校正调节器是不断地测量y伙)u(k),不断地辨识对象中的参数不断 修改控弘伙)过程。此算法由辨识和控制组成，其辨识和控制过程实际上是一系列 在线递推算法。由于需要对信息进行存储、处理和运算，而口要保证数据的实时 性，因此这种功能必须用计算机完成。3.3.1递推最小二

30、乘估计设一个线性系统，其输入为心心，兀3,£，输出y，输入输出均可测量，并 口它们的关系为：y = 0內+&2丘+仇兀其中，= 2,7?)未知，分别在时刻人心，®进行m次测量，可的一下 m个方程：y(l) = &內(1) + 02x2 (1) + + 3tjxn (1) y=臥+如2+仇捡ym = 0內(m) + 02x2 o) + + 0nxn (m)用矩阵表示上述m个方程，则有：y = x*e其中，& = &i020nlx|(l)x2(l)x/l)_%!(2)兀2兀“x =. 兀1(加)x2(m)xn(m)_y = y(l) y(2)y(

31、m)y以上方程表明，系统有n个参数，n个输入量。由于数据中有测量噪声或模 型误差影响，故上式应表示为：y = x*&+£其中，e = e(l) e(2)e(n)y称为误差向量。令j =ete = (y-x0)t(y-x0)为了估计未知参数，使丿最小，令® = (),计算可得：de人0 = (xtxyxty(3-1)a当xx为非奇异矩阵时，称式(31)为&的最小二乘法估计值。若已取得m组数据，可知其第m次辨识8结果的最小二乘估计何为：a进行第m+l次数据采集，则第m+1输入输出方程为：y(m + 1) = &內(m +1) + 2x2(m+ !) +

32、+ 0nxn (m + 1)(3-2)即:y(m + l) = x7'(加+ 1)0(3-3)式中：x 丁 (m + 1)=州(/? +1), x2 (m +1),xn (m +1)即：y(m + l) = xt(m + l)0a根据式(3-1),第m+1次辨识结果&的最小二乘估计可写为:a0(m) = (x：xmrxx(3-4)(3-5)(3-6)式中:西(1)兀2d)e.兀】(2)勺耳兀| (加)x2(m)x” (加)x,(77? + 1)(777 + 1)xmx爲x曲y(加)y(m +1)令矩阵p(m) = (x：xj'，则:pm +1) = (x爲 xqt =

33、p-1 (m) + x*x 爲r1根据矩阵求逆定理：(a + bcd)-1 = a_,- a 'ficc-1 + da lby da式(37)可写为:p(m +1) = p(2)- k (加 + 1)x 气加 + 1)p(加)式中，k(m 4-1) = p(m) - x (加 +1) / 1 + x 气加 + 1)p(加)x (加 +1) 则式(37)可写为:aaa0(jn + 1) = 0m + 1) + k(m + v)y(m 4-1) - x r (/n 4-1) 0(加)以上三式构成了递推最小二乘估计法。(3-7)(3-8)(3-9)(3-10)由式(310河知，新的估计值0(

34、m +1)可在上一步的估计值(加)的基础上加 以修正得到，而修正项是本次估计课差乘以加权系数，特别是式(39)中的分母 l + xr(m + l)p(m)x(m + l)是一标量。因此递推算法避免了矩阵求逆运算，可大 大缩短运算速度，适用于实时辨识和控制。3 .3 .2带遗忘因子的递推最小二乘算法的递推算式系统在线参数辨识中，人多数人以为数据越多，估计越准确，其实不然。当 数据增加到一定程度后参数辨识反而不准确，这就是所谓数据饱和问题。在递推 辨识中，最好利用当前数据，因为它真正反映系统当前的动态过程，数据越陈【口 则偏离当前动态特性越远，估计越不准确。遗忘因子算法的主要思路是对新老数据给予不

35、同的对待，逐渐遗忘老数据的影响。具体做法是每取得一个新的)5 + 1)时，将以前的所有数据 乘以一个小于1的加权因子/(0</<1),则有:yxm儿+i则式(37河写为:p® +1)=(x爲 x=r2p- (m) + x讪 x爲一】(3-11)利用矩阵求逆公式，可得:+1) =- kqn + l)xr(/z? + l)p(m)-/其中:k (加 + 1)=p(加)x0 + 1)/2 + xr(m4-l)p(m)x(/7? + l)令则有:p = 72p(m + 1) = p(m) 一 k(m + l)xt(m + 1)p(m)p(3-12)+ 1)p + xt (m- )

36、p(m)x (m +1)(3-13)则。被称为遗忘因子(°s<1)。此时：式(3-11). (312)和(3一起构成了带遗忘因子的递推最小二乘法推 导公式。一般来说，°的取值范圉在0.95至0.995之间。对于遗忘因子p的不同a取值，可以得到不同的遗忘效果：°值较小时&(加)的估计跟踪参数时变能力强, 但是噪声干扰影响造成的估计波动大。3.3 .3初值的确定初值的确定有两种方法"3a(1) 先用一般最小二乘公式(36)求出初始估计值3(加)和p(m),然后从第 (m+1)步开始递推计算。(2) 预先设定初始值为：0二任意值，p(0) = a

37、1!/选一个很大的数，在103106之间，/为单位矩阵。本文采用第二种方法设 定初值。3.4本论文所用自校正pid控制算法3 .4.1具体框图和原理本论文中所讨论的被控对象是已知的，不需要对其进行在线估计参数，使用 pid调节器进行常规控制；所用的算法是带遗忘因子的递推最小二乘法。木文基木思想是：给出一个具体系统和一个闭环系统希望的传递函数，利用 pid自校正调节器不断调节，使得输出y尽量达到期望值。其具体结构原理图如下图所示：图3. 2实验锅炉自校正ptd结构原理图递推最小二乘法算法：用差分方程描述的单输入、单输出过程为：(3-14)其中：y上表示过程输入和输出测量值，）心）上表示过程输出和

38、输入测量 值，£为滞后时间，z"为后移位算子，a（zl） = l + alzl +a2z2 '- + anzn,bzx） = % + 方忆“ + 乞疋2 + + bnzn假设是稳定的，且z>1。图3.3单输入单输出闭环系统由公式（27）可知，数字pid调节器的单传递函数为：g = kp +- + k“（l-厂）e 1-z'1°二kp + kj + kj -區 + 2忑门厂+ k皿彳 t（l-z,）令 g = k» a. = kpt + 2kdt, a3 = kpt + k/ + kdt 则 g(z) =a3-a2zl +a,z2卩(

39、1一厂)图3.3所示的单输入单输岀闭环系统的传递函数是:叫=zkb(zl)a3+a2zl +axz2l4(z") + f b(z)43 + a?+ &忆边(3-15)(3-16)(3-17)%同时可以用指定闭环系统的极点和零点的方法来表示：w （-）二 0（?）二z"b（zt）如 + 仅2厂1 +可厂2 ' _ f（z“）_（厂|） + 疋*8（厂|）&3 + &2厂'+少厂2化简得：厂（巾+勺厂+%3戶（已-q（）b（zs（/）"（l 厂）0（旷）仏为呛）又因故有:(勺+如厂+%-2)円厂)-0(厂)曲)"(1-

40、厂)0(厂)叨)(3-18)令：亍=p(门)_0(广)曲)则式(3 可化简为：u(t) = a3 y(t) + a2 y(t -1) + cr, y(t - 2 :(3-19)如果知道可。2。3的值，就可以根据式（3计算出pid调节器的三个参数 kp、k,和心的值。采用递推最小二乘法估计少、勺 巾的数值。令9t =a3,a2,axyxt =刃)，刃一 1),刃一 2) 则式(3 变为：u(t) = xt 0定义误差矢量e = u(t)-xt0j =st£= lu(t)-xr0irlu(t)-xt3=(u(r)2 - 2u(t)xt0 + otxxt0若(x/f1满秩，为使丿最小，令竺

41、=0 de 则有：刃 _ 0(历()2-2历(0*0 + &7%”0) 5060a=-22w(r)xr + 2xxt 0=0有：a& =( xxtyxtu(t)令 p =(xxty 则：a e = pxtu(t)对于 w(r) = xt0(3-20)(3-21)(3-22)进行了 n次观测，有即:其中,瓦(0 = x & u2(t) = x0un(t) = x0 ,x：=_<(0_ 理(/) _4(r)_un第n+1次的测量获得一个新的方程：% (0 =易+则这n+1次的测量可表示为:其中:百弘)%(/) 屮n+i() 根据式(3-22),有：u兀"+1

42、(')“n+13+1 =(心 +=府+兀柏边+ef根据矩阵求逆公式对上式求逆，得：pn+ = pn - pnxn+1(01 + 4+1(0pnxn+1 (r)'1 xj+1 (t)pn(3-23)(3-24)根据式(3-23) w：a0n+1 = pn+mn+mn+1aa=on + pnxn+l (r)l + 兀+1 (t)pnxn+ (/)爲+i -珂+i (0 on (3-25)aaaaaa计算出&后,根据& = *2心,即知道了 &3o2s将其直接带人式(3-16)得到t时刻的pid调节器:u(f)q3+q2 z l +qi(3-26)根据 a=k

43、d,a2 = kpt + 2kdt,a3 = kpt + k,t + kdt 计算出 k。、k,心代入下式：ku(t)=kp+-+kd(l-z-')e(t)l-z也可求得t时刻上的pid调节器他。3.4.3带遗忘因子的递推最小二乘法参数的递推最小二乘估计与一次完成的最小二乘估计是数学等价的，它们都 仅适用于估计时间定常过程的参数，而不适用于估计时变过程的参数。时变过程 的特点是过程的参数可能随着吋间变化而改变。因此，它的数学模型参数具有“吋 间性”，在利用动态过程的输入输出数据来辨识模型参数时，“老”的数据往往 只能反映“老”的过程参数；而改变后的“新”参数，要靠用新的和比较新的实 验

44、数据来估计。因此时变过程参数估计的特点是，不同时段的实验数据的作用是有区别的。 时变过程的参数估计有多种不同的算法，本文所用的是带遗忘因子的递推最小二 乘法。考虑到遗忘因子，则(3-28)(3-29)&n+1 = &n + pgg (t)p + 对+ (r)恥n+1 广叵n+1(0 兀：+1 (0 °n 本论文采用的就是基于带遗忘因子的递推最小二乘法设计的自校正pid控制4系统硬件的结构设计4.1系统硬件的结构该系统以单片机为核心，包括模拟信号处理电路（热电偶放大电路、线性光 电隔离电路、多路模拟开关、环境温度测量电路）、a/d转换电路、rs232接口 电路、输岀光电

45、隔离电路、控制输岀电路以及电源电路等几部分组成图4.1自校正pid实验锅炉温度控制系统控制图实验锅炉炉温的控制过程如下：测温元件将检测到的温度信号经过放大、 隔断处理后，送到多路选择器，多路选择器在单片机cpu的控制下对这组信号 进行选择，在某一时刻被选择的信号送到a/d转换器，转换成数字量送给单片 机，单片机对a/d结果处理（包括冷端补偿计算），得出加热炉温度测量值，再 将温度测量值与给定温度值进行pid运算，产生控制脉宽输出，来控制加热炉 的加热过程。4.2自校正pid实验锅炉控制系统原理被控对象为220v交流电阻加热实验锅炉。实验锅炉与加热器对象数学模型, 水泵是用来使水活动这样可以均匀

46、加热防止加热不均匀测量而产生误差。图4.22实验锅炉炉温测定飞升曲线设阶跃响应曲线为s形的单调曲线，如图4 .22所示，设阶跃输入幅值为,则增益kp可按下式求取：© =气 x0)时间常数厂及延迟时间：可用作图法确定：在图4.22响应曲线的拐点”做切 线，切线与时间轴交与a点，而与响应曲线稳态值交与b点，贝i()4对应延迟时 间一 对应时间常数t。响应传递函数为：gp =心+ 1测得实验锅炉为具有纯滞后的一阶惯性环节，bpl4：2&"178s + 1带零阶保持器的广义对象脉冲传递函数为gp(c) = z(严,k&严=kpd。)严)=h(z'd *tps

47、 + l 一1 + 计1一为(4-1)取厂=，则n = l,d = 2。将实测被控对象参数带入式(4-21)中得:/?0 =0.5636, =-0.7987从极点配置的观点出发，以典型的二阶系统闭环传递函数的离散特征多项式丽八厂)=1- 2不叫cos cdt-£2z + 严才2(4-2)©为无阻尼自然震荡角频率，£阻尼比。对于式(4-2),当二阶系统最佳阻尼比£ = 0.707吋，在单位阶跃作用下的超调 量 =4.3%,相角稳定裕量/(.) = 65.5它为二阶最佳动态响应模型。采样周期卩和哄£的关系：t =丁 =10- 20)当 n丁 =10

48、。当 £=0.7077 = 40$时，得©=0.022/ 式(422)的期望特征多项式为：)=1 - 0.7812厂 + 0.2882厂 建立期望的传递函数：时)0.509p(厂)1-0.7812厂+0.2882 厂(4-3)4.3自校正控制算法设计实验锅炉的模型公式由w(宁*吿亠鴿右y(/) + 0.7987(/ 一 1)二 0.5636«(r- 2) 闭环系统希望的僦函数为：0.509”厂)u(z'1)1-0.7812 厂 + 0.2882z-2在递推估计调节参数时，设置遗忘因子q = 0.995。4.4自校正pid算法设计流程图图4.4算法设计流程图

49、5 matlab仿真及结果分析设输是阶跃信号输入，如图5.1图5系统输入信号图5.2系统期望输出图5.3系统经自校正pid调节后输出（初值为2, 2, 2）图5.4系统经常规pid调节输出的波形(初值为2, 2, 2)在matlab环境下可得上述闭环系统的期望输出和经过自校正pid调节后的 系统实际输出。由图(5.3河得自校正能够对期望值进行较好的跟踪。最终：_ 0.4248 _ae= -0.07070.2419为了更好的进行对比，在matlab环境下进行常规pid调节后的输出为图(54) 很明显，在在相同的初始值下，经常规pid调节的系统不稳定了把初始值设为_ 0.4248 _a0= -0.

50、0707则常规pid调节实验锅炉后的如图(5.5)。0.2419图5.5初始值设为&常规pid调节的图形由图可知自校正有自动调节pid参数的能力，能自适应被控制的对像。为 了验证自校正pid控制器的自动跟随能力，设系统输入该为方波信号为一个周 期200s,幅值为1的矩形波，见下图：图5.6系统输入波形图5.7期望输出的波形图5.8系统经自校正pid输岀的波形图5.9常规pid锅炉炉温控制系统输出波形由图(5.6)(5.7)(5.8)可得自校正pid控制的锅炉可以很好的跟随输入，而常规apid即使应用了自校正得到的&值，虽然稳定了但从波形可以看出跟随系统输入 的效果很差。_20_

51、为了检验是否自校正pid的初始值可以任意设定，改变初始值用155图5.10系统自校正pid控制实验锅炉系统输出波形(初始值为g )由图(5.10)可以看出虽然在刚开始系统有很大的震荡，但在有限个震荡后系统输 _ 0.4247 _aa a出为1跟随了系统的输入，最终的-0.0707 ee说明自校正pid控制实0.2419现了实验锅炉温度的是最优控制。6结论本论文认真地研究和分析了实验锅炉温度控制系统，利用自校正pid调节 器，对递推最小二乘算法进行了推导。针对给实验锅炉炉温传递函数和期望传递 函数的闭环系统，编制了自校正pid调节器参数调整过程的流程图，并利用 matlab程序进行了编程和仿真实

52、验。将其与普通自校正调节器进行了对比实 验，验证了自校正pid调节器的优越性。并针对系统可能产生的一些情况变化, 比如初值的改变、系统模型的改变、输入的改变，进行了仿真实验。通过实验证 明，自校正pid调节器控制炉温温度和普通pid调节器控制炉温稳定相比，对 于自校正系统具有很好的适应能力。参考文献fl王顺晃.非线性pid算法及其在电阻加热炉集散控制系统屮的应用卩自动化学报， 1995. (6):675681 .2 闰秀英，任庆昌，盂庆龙.一种自校正pid控制器设计与仿真研究j.系统仿真报， 2006, (2): 753756.邓娟，从爽.控制工程中几种自校正控制的应用及其性能对比研究jj.电

53、子学报，2008,(3) : 1418.4 张显库，贾欣乐.闭环增益成形控制新算法及其应用j.电子学报，1999, 27(11):133135.f5 tjokro s, shah s l. adaptive pid controlfj. proceeding of the 1987 american control confemce, 1985(2): 15281534. 刘伯春.离散系统采样周期的选择j.化工自动化及仪表，1989, (1): 2833.7 astrom k. j. witten mark b. adaptive control, addison-wesley, 1989 李清泉等译，自适 应控制m.北京：科学出版社，1992.8 刘东升.pid自校正调节器(j).现代防御技术，1991