实验锅炉炉温自校正pid控制系统设计

毕业设计(论文) 题目：实验锅炉炉温自校正PID控制系统设计 摘 要本文以递推最小二乘法为自适应规律, 研究实验锅炉温度控制, 用MATLAB语言编程并进行系统仿真, 在仿真结果的基础上进行分析研究. 结果表明当采用自校正PID算法时, 系统有自适应能力, 能根据被调节系统自动调节KP、KI、KD 参数使系统达到稳定, 解决了长期以来大时滞实验锅炉系统PID参数设定难的问题. 为了使系统更精确, 本文采用了“带遗忘因子的递推最小二乘法”估计算法.关键字： 自校正PID; 炉温控制; 参数估计; 最小二乘法 ABSTRACTThis paper adopts recursive least square method to 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 conditions 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 with forgetting factor.Keywords： Self-revised PID Controller; Temperature Control; Recursive Least Square Method; Parameter Estimation目 录摘 要IABSTRACTII1 绪论11.1 选题背景及意义11.2 论文的主要内容11.3 本人主要工作12 常规PID控制算法及其改进算法22.1 常规PID控制原理22.1.1 常规PID调节器算法22.1.2 常规PID调节器的参数整定32.1.3 常规PID调节器在实际应用中的局限32.2 数字PI D控制32.2.1 位置式PID控制算法42.2.2 增量式PID控制算法52.2.3 数字PID控制器的参数整定方法52.2.4 采样周期的选择63 自校正PID控制算法73.1 自适应控制系统原理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.1 系统硬件的结构184.2 自校正PID实验锅炉控制系统原理184.3 自校正控制算法设计204.4 自校正PID算法设计流程图215 MATLAB 仿真及结果分析226 结论26参考文献27致 谢28附 录291 绪论1.1 选题背景及意义锅炉炉温控制是典型的工业过程控制对象，例如：在冶金工业、化工生产、电力工程、造纸行业、机械制造和食品加工等诸多领域中，人们都需要对各类加热炉、热处理炉、反应炉和锅炉中的温度进行检测和控制。通常，电阻炉炉温控制都采用偏差控制法。偏差控制的原理是先求出实测炉温对所需炉温的偏差值，然后对偏差值处理获得控制信号去调节电阻炉的加热功率，以实现对炉温的控制。通常对偏差进行比例、积分和微分控制又称PID控制，是工业控制过程中 应用最广泛的一种控制形式。但由于锅炉温度控制具有升温单向性，大惯性，纯滞后性等特点，很难用数学建立精确的模型和确定参数。当被控对象的结构和参数不能完全掌握，或得不到精确的数学模型时，控制理论的其它技术难以采用时，系统控制器的结构和参数必须依靠经验和现场调试来确定，多年来，人们一直在努力寻找一种方法，能够简单易行的在线实时自动整定PID参数。自校正PID调节器就是其中的一种，自校正PID是自校正思想和常规PID控制器相结合的产物，它吸收了两者的优点，设计参数少，能够在线整定和校正PID控制器参数，具有较强的适应能力1。1.2 论文的主要内容本论文详细的阐述了常规PID控制优点和缺陷控制功能，以及如何把最小二乘法的思维与常规PID相结合的，得出了基于最小二乘法的自校正PID控制算法。以及如何把自校正PID应用到实验锅炉稳定控制系统中，并对自校正PID实验锅炉温度控制系统进行了matlab仿真，并与常规PID控制的实验锅炉炉温控制系统相比较，得出了自校正PID能自适应实验锅炉系统且自动调节PID参数以达到最优控制，解决了常规PID参数难以设定的问题。1.3 本人主要工作(1) 对自适应控制技术和常规PID调节器的原理进行研究，熟悉其模型原理；(2) 确定整个实验锅炉自校正PID控制系统算法；(3) 设计实验锅炉对象模型组成框图和结构图；(4) 对实验锅炉模型实施自校正PID控制，进行matlab仿真，得到多个设定值的仿真波形，并与其他控制方法比较得出结论。2 常规PID控制算法及其改进算法PID控制是模拟控制系统最常用的控制规律之一。模拟PID控制系统原理框图如图2.1所示，系统由模拟PID控制器和被控对象组成。比例R(t)+ +Y(t)微分积分被控对象+_U(t)图2.1常规PID控制原理图2.1 常规PID控制原理 常规PID调节器算法PID控制器是一种线性控制器，它根据给定值与实际输出值构成控制偏差： (2-1)将偏差的比例P、积分I和微分D通过线性组合构成控制量，对受控对象进行控制。其控制规律为： (2-2)其传递函数为： (2-3)式中，比例系数，积分时间常数，微分时间常数。PID控制器各校正环节的作用如下：(1) 比例环节：成比例地反映控制系统的偏差信号，偏差一旦产生，控制器立即产生控制作用，以减小偏差。(2) 积分环节：主要用于消除误差，提高系统的精度。积分作用的强弱取决于积分时间常数，越大，积分作用越弱，反之则越强。(3) 微分环节：反映偏差信号的变化趋势(变化速率)，并能在偏差信号变得太大之前，在系统中引入一个有效的早期修正信号，从而加快系统的动作速度，减少调节时间。 常规PID调节器的参数整定PID控制器的参数整定是控制系统设计的核心内容。它是根据被控过程的特性确定PID控制器的比例系数、积分时间和微分时间的大小。PID控制器参数整定的方法很多，概括起来有两大类：一是理论计算整定法。它主要是依据系统的数学模型，经过理论计算确定控制器参数。这种方法所得到的计算数据未必可以直接用，还必须通过工程实际进行调整和修改。二是工程整定方法，它主要依赖工程经验，直接在控制系统的试验中进行，并且方法简单、易于掌握，在工程实际中被广泛采用。PID控制器参数的工程整定方法，主要有稳定边界法、动态特性参数法和衰减曲线法。三种方法各有其特点，其共同点都是通过试验，然后按照工程经验公式对控制器参数进行整定。但无论采用哪一种方法所得到的控制器参数，都需要在实际运行中进行最后调整与完善。现在一般采用的是稳定边界法。利用该方法进行PID控制器参数的整定步骤如下2：(1) 首先预选择一个足够短的采样周期让系统工作；(2) 加入比例控制环节，直到系统对输入的阶跃响应出现临界振荡，记下这时的比例放大系数和临界振荡周期；(3) 在一定的控制度下通过公式计算得到PID控制器的参数。 常规PID调节器在实际应用中的局限(1) 由于实现控制系统的器件的物理特性的限制，使得PID控制器获取得原始信息偏离真实值，而其产生的控制作用偏离理论值。(2) 由于在系统的设计与整定过程中，要兼顾动态与稳定性能，只能采取折衷方案，难以大幅度提高控制系统的性能指标。(3) 对于存在强非线性、快速时变不确定性、强干扰等特性的对象，控制效果较差。2.2 数字PI D控制在连续时间控制系统中，PID控制器应用得非常广泛。其设计技术成熟，长期以来形成了典型的结构，参数整定方便，结构更改灵活，能满足一般的控制要求。计算机控制是一种采样控制，它只能根据采样时刻的偏差值计算控制量。因此连续PID控制算法不能直接使用，需要采用离散化方法。在计算机PID控制中，使用的是数字PID控制器。数字PID控制比连续PID控制更为优越，因为计算机程序的灵活性，很容易克服连续PID控制中存在的问题，经修正得到更完善的数字PID算法3。 位置式PID控制算法按模拟PID控制算法，以一系列的采样时刻点代表连续时间，以矩形法数值积分近似代替积分，以一阶向后差分近似代替微分，即： （2-4）可得离散PID表达式： (2-5)式中，为采样周期，和分别为第(k-1)和第k时刻所得的偏差信号。根据Z变换迭值定理和滞后定理可得：则对离散PID的表达式(2 -5)作z变换，可得 (2-6)由(2-6)可得数字PID控制器的Z传递函数为： (2-7)算法中，为了求和，必须将系统偏差的全部过去值都存储起来。这种算法得出控制量的全量输出(的是控制量的绝对数值，的值和执行机构的位置是一一对应的，因此称其为位置式PID控制算法4。2.2 .2 增量式PID控制算法位置式PID算法计算时需要对进行累加，计算机运算工作量很大。而且，因为计算机输出的对应的是执行机构的实际位置，如果计算机出现故障，的大幅度变化会引起执行机构位置的大幅度变化，这就可能造成重大的生产事故。并且有些执行机构需要的不是控制量的绝对值，而是控制量的增量(例如去驱动步进电动机)时，需要的就是PID的增量算法。所谓增量式PID是指数字控制器的输出只是控制量的增量。根据位置式算法的公式(2-5)可知，则有： (2-8)式(2-8)称为增量式PID算法，对式(2 -8)进行归并后，得： (2-9)其中，式(2-7)已经看不出是PID的表达式了，也看不出P、I、D作用的直接关系，只表示了各次误差量对控制作用的影响。由式(2-9)看出，数字增量PID算法，只要存储最近的三个误差采样值就足够了。由于增量式控制算法不需要累加，控制增量仅与最近3次采样有关，所以误动作时影响小，而且较容易通过加权处理获得比较好的控制效果5。2.2.3 数字PID控制器的参数整定方法数字PID控制器控制参数的选择，可按连续一时间PID参数整定方法进行。在选择数字PID参数之前，首先应该确定控制器的结构。对允许稳态误差的系统，可以适当选择P或者PD控制器，使稳态误差在允许的范围内。对于必须消除稳态误差系统，应选择包含积分控制的PI或者PID控制器。一般来说，PI、PI和PD控制器应用较多。对于有滞后的对象，往往都加入微分控制。控制器结构确定后，即可开始选择参数。参数的选择，是要根据受控对象的具体特性和对控制系统的性能要求进行。工程上，一般要求整个闭环系统是稳定的，对给定量的变化能迅速响应并平滑跟踪，超调量小；在不同干扰作用下，能保证被控量在给定值；当环境参数发生变化时，整个系统能保持稳定等等。这些要求，对于控制系统自身性能来说，有些是矛盾的。我们必须满足主要方面的要求，兼顾其他方面，适当地折中处理。PID控制器参数的整定，可以不依赖于受控对象的数学模型。工程上，PID控制器的参数常常是通过实验来确定，通过试凑，或者通过实验经验公式来确定6。2.2.4 采样周期的选择在采样数据控制系统中，如果设采样周期为r，则采样数率为1T，采样角频率采样周期丁是设计者要精心选择的重要参数，系统的性能与采样周期的选择有密切关系。采样周期的选择受多方面因素的影响，主要考虑的因素分析如下：(1)香农(Shannon)采样定理7 (被采样信号的上限角频率)给出了采样周期的上限。满足这一定理，采样信号方可模拟或者近似地恢复为原模拟信号，而不丢失主要信息。在这个限制范围内，采样周期越小，采样一数据控制系统的性能越接近于连续一时间控制系统。(2)闭环系统对于给定信号的跟踪，要求采样周期要小。(3)从抑制扰动的要求来说，采样周期应该选择得小一些。(4)从执行元件的要求来看，有时要求输入控制信号要保持一定的宽度。(5)从计算机精度考虑，采样周期不宜过短。(6)从系统成本上考虑，希望采样周期越长越好。综合上述因素，选择采样周期，应在满足控制系统的性能要求的条件下，尽可能地选择低的采样速率。可以根据所设计的系统的具体情况，采用试凑的方法，在试凑过程中根据各种合理的建议来预选采样周期，多次试凑，选择性能较好地一个作为最后的采样周期。3 自校正PID控制算法3.1 自适应控制系统原理3.1.1 概述在许多工程中，被控对象和过程的数学模型事先难以确定，即使在某一条件下已被确定了的数学模型，在条件改变后，其动态参数或其模型结构仍经常发生变化。因而，常规调节器不能得到好的品质，这就需要设计一种特殊的控制系统，它能够自动地补偿各方面的非预知变化，“自适应控制系统”也就应运而生。自适应控制系统(Adaptive Control Systems)是现代控制理论的一个重要分支。所谓自适应控制系统是指能在系统和环境的信息不完备的情况下改变自身特性来保持良好工作品质的控制系统，又称适应控制系统 。由于自适应控制的对象是那些存在不定性的系统，所以这种控制应首先能在控制系统的运行过程中，通过不断地量测系统的输入、状态、输出或性能参数，逐渐了解和掌握对象 。然后根据所得的过程信息，按一定的设计方法，做出控制决策去更新控制器的结构、参数或控制作用，以便在某种意义下使控制效果达到最优或次最优，或达到某个预期目标。按此设计思想建立的控制系统便是自适应控制系统 。由此可见，一个自适应控制系统必然具有下列三个基本特征8：(1) 过程信息的在线积累在线积累过程信息的目的，是为了降低受控系统原有的不定性。为此可用系统辨识的方法在线辨识受控系统的结构和参数，直接积累过程信息;也可通过量测能反映过程状态的某些辅助变量，间接积累过程信息。在系统辨识中，结构辨识比参数估计困难得多。(2) 可调控制器可调控制器是指它的结构、参数或信号可以根据性能指标要求进行自动调整。这种可调性要求是由受控系统的不定性决定，否则就无法对过程实现有效控制。(3) 性能指标的控制性能指标的控制可分为开环控制方式和闭环控制方式两种。若与过程动态相关联的辅助变量可测，而且此辅助变量与可调控制器参数之间的关系又可根据物理学的知识和经验导出，这时就可通过此辅助变量直接调整可调控制器，以期达到预定的性能指标。这就是性能指标的开环控制，它的特点是没有根据系统实际达到的性能指标再作进一步的调整。与开环控制方式不同，在性能指标的闭环控制中，还要获取实际性能与预定性能之间的偏差信息，将其反馈后修改可调控制器，真到实际性能达到或接近预定性能为止。3.1.2 模型参考自适应控制模型参考自适应控制系统(Model Reference Adaptive System)，简称MRAS，由以下几个部分组成：参考模型、被控对象、反馈控制器和调整控制参数的自适应机构等组成。在这个系统中，把参考输入r同时加到参考模型和可调系统的输入端，用减法器将参考模型和可调系统的输出直接相减，得到广义误差信号e。自适应机构按一定准则利用广义误差信号来修改可调系统的调节器参数，或产生一个辅助输入信号，使广义误差的某个指标达到极小。当可调系统特性与参考模型特性逐渐逼近，广义误差趋于极小或下降为零，调节过程结束。当对象特性在运行中发生了变化，控制器参数的自适应调整过程与上述过程完全一样。设计这类自适应控制系统的核心问题是如何综合自适应调整律，即自适应机构所应遵循的算法。关于自适应律的设计目前存在两种不同的算法，一种称为局部参数最优化的方法，即利用梯度或者其他参数优化的递推算法，求得一组控制器的参数，使得某个预定的性能指标达到最小。这种方法的缺点是不能保证参数调整过程中，系统总是稳定的。自适应律的另一种设计方法是基于稳定性理论的方法，其基本思想是保证控制器参数自适应调节过程是稳定的，然后再尽量使这个过程收敛快一些。这类系统的难点在于系统稳定性分析，李雅普诺夫稳定性和波波夫超稳定性理论都是设计自适应律的有效工具。3.1 .3 自校正控制自校正控制是70年代发展起来的一种随机自适应控制，产生背景是：工业过程控制由于强随机干扰、模型未知、参数时变、大时滞等因素，导致常规的控制方法效果差。这类自适应控制系统的一个主要特点就是具有一个被控对象数学模型的在线辨识环节，具体地说就是加入了一个对象参数的递推估计器。它根据系统的运行，首先对被控对象进行在线辨识，然后再根据辨识出的模型参数和事先指定的性能指标，在线地综合控制作用。当对象参数已知时，对调节器的参数进行在线求解。由于调节器的控制律是多样的，参数估计的方法也是多样的，因此自校正调节器非常灵活，采用各种不同的控制方法和估计方法来搭配，能满足不同的性能要求，其难点在于收敛性。本文主要采用自校正控制的PID调节器，所采用的控制律是PID调节，参数估计方法是最小二乘法。3.2 自校正控制系统U(t)调节参数估计参数设计机构估计器R(t)被控对象PID调节器Y(t)图3.1 自校正控制系统结构图自校正控制系统结构图如图3.1所示：系统由三个部分组成：参数估计器、控制器参数计算、调节器。其中，参数估计器又称辨识器，用来根据测量得到的系统的输入输出信息，不断地在线辨识系统的模型结构和参数;控制器参数计算是指按照一定的控制算法，由辨识得到的参数，取得响应控制算法的参数;调节器则是形成控制量，送往被控过程或对象，达到控制目标9。自校正控制系统有两个反馈回路，外环调节器以对象输入和输出的反馈为基础，而内环调节器具有外环调整的可调参数。参数估计和控制设计必须在线实现，因此参数估计必须采用递推算法，控制器设计必须采用计算尽量简单的设计方法。该系统能自动校正自己的参数以得到希望的闭环系统特性。图3.1所示的自校正调节器的参数是经由参数估计和控制的设计计算而间接进行更新的。但也可以将对象重新参数化，即利用调节器的参数来表示模型。这时，就不需要进行设计计算这个环节。算法将大大简化，设计机构的方框将不复存在，调节器参数将直接更新。本文研究的自校正PID调节器就是直接更新参数。自校正控制系统结构形式通常分为以下两种：(1) 显式结构。如图3.1所示，首先估计出控制对象的参数，然后进行控制器设计，这样的机构称为显式结构。(2) 隐式结构。将控制器对象的参数估计和控制器设计这两个步骤结合在一起，直接估计出控制器参数，从而大大简化自校正控制的算法，这样的结构形式称为隐式结构。自校正调节器是不断地测量 ，不断地辨识对象中的参数。不断修改控过程。此算法由辨识和控制组成，其辨识和控制过程实际上是一系列在线递推算法。由于需要对信息进行存储、处理和运算，而且要保证数据的实时性，因此这种功能必须用计算机完成。3.3.1 递推最小二乘估计设一个线性系统，其输入为，输出，输入输出均可测量，并且它们的关系为：其中，未知，分别在时刻进行m次测量，可的一下m个方程：用矩阵表示上述m个方程，则有：其中，以上方程表明，系统有n个参数，n个输入量。由于数据中有测量噪声或模型误差影响，故上式应表示为：其中，称为误差向量。令为了估计未知参数，使最小，令，计算可得： (3-1)当为非奇异矩阵时，称式(3-1)为的最小二乘法估计值。 若已取得m组数据，可知其第m次辨识结果的最小二乘估计10为：进行第m+1次数据采集，则第m+1输入输出方程为： (3-2)即： (3-3)式中： (3-4)即： (3-5)根据式(3-1)，第m+1次辨识结果的最小二乘估计可写为： (3-6)式中：令矩阵，则： (3-7)根据矩阵求逆定理：式(3-7)可写为： (3-8)式中， (3-9)则式(3-7)可写为： (3-10)以上三式构成了递推最小二乘估计法。由式(3-10)可知，新的估计值可在上一步的估计值的基础上加以修正得到，而修正项是本次估计误差乘以加权系数，特别是式(3-9)中的分母是一标量。因此递推算法避免了矩阵求逆运算，可大大缩短运算速度，适用于实时辨识和控制。3 .3 .2 带遗忘因子的递推最小二乘算法的递推算式系统在线参数辨识中，大多数人以为数据越多，估计越准确，其实不然。当数据增加到一定程度后参数辨识反而不准确，这就是所谓数据饱和问题。在递推辨识中，最好利用当前数据，因为它真正反映系统当前的动态过程，数据越陈旧则偏离当前动态特性越远，估计越不准确。遗忘因子算法的主要思路是对新老数据给予不同的对待，逐渐遗忘老数据的影响。具体做法是每取得一个新的时，将以前的所有数据 乘以一个小于1的加权因子则有：则式(3-7)可写为： (3-11)利用矩阵求逆公式，可得：其中：令则有： (3-12) (3-13)则被称为遗忘因子。此时：式(3-11)、(3-12)和(3-13)一起构成了带遗忘因子的递推最小二乘法推导公式。一般来说，的取值范围在0.95至0.995之间。对于遗忘因子的不同取值，可以得到不同的遗忘效果：值较小时的估计跟踪参数时变能力强，但是噪声干扰影响造成的估计波动大。3.3 .3 初值的确定初值的确定有两种方法11：(1) 先用一般最小二乘公式(3-6)求出初始估计值和，然后从第 (m+1)步开始递推计算。(2) 预先设定初始值为：=任意值，选一个很大的数，在103106之间，为单位矩阵。本文采用第二种方法设定初值。3.4 本论文所用自校正PID控制算法3 .4 .1 具体框图和原理本论文中所讨论的被控对象是已知的，不需要对其进行在线估计参数，使用PID调节器进行常规控制；所用的算法是带遗忘因子的递推最小二乘法。本文基本思想是：给出一个具体系统和一个闭环系统希望的传递函数，利用PID自校正调节器不断调节，使得输出y尽量达到期望值。其具体结构原理图如下图所示：PID调节器实验锅炉最小二乘参考模型Y+_Euyer_图3.2实验锅炉自校正PID结构原理图递推最小二乘法算法： 用差分方程描述的单输入、单输出过程为： (3-14)其中：表示过程输入和输出测量值，表示过程输出和输入测量值，为滞后时间，为后移位算子，假设是稳定的，且。 E(t)R(t)U(t)Y(t)图3.3单输入单输出闭环系统由公式(2-7)可知，数字PID调节器的单传递函数为：令 (3-15)则 (3-16)图3.3所示的单输入单输出闭环系统的传递函数是： (3-17)同时可以用指定闭环系统的极点和零点的方法来表示：化简得：又因故有： (3-18)令：则式(3-18)可化简为： (3-19)如果知道的值，就可以根据式(3-15)计算出PID调节器的三个参数、和的值。采用递推最小二乘法估计的数值。令 则式(3-19)变为： (3-20)定义误差矢量若满秩，为使最小，令则有： （3-21）有： (3-22)令则：对于 进行了N次观测，有即：其中，第N+1次的测量获得一个新的方程：则这N+1次的测量可表示为：其中：根据式(3-22)，有： (3-23)根据矩阵求逆公式对上式求逆，得： (3-24)根据式(3-23)有： (3-25)计算出后，根据，即知道了，将其直接带人式(3-16)得到t时刻的PID调节器： (3-26)根据计算出代入下式： 也可求得t时刻上的PID调节器12。3.4.3 带遗忘因子的递推最小二乘法参数的递推最小二乘估计与一次完成的最小二乘估计是数学等价的，它们都仅适用于估计时间定常过程的参数，而不适用于估计时变过程的参数。时变过程的特点是过程的参数可能随着时间变化而改变。因此，它的数学模型参数具有“时间性”，在利用动态过程的输入输出数据来辨识模型参数时，“老”的数据往往只能反映“老”的过程参数；而改变后的“新”参数，要靠用新的和比较新的实验数据来估计。因此时变过程参数估计的特点是，不同时段的实验数据的作用是有区别的。时变过程的参数估计有多种不同的算法，本文所用的是带遗忘因子的递推最小二乘法。考虑到遗忘因子，则 (3-28) (3-29)本论文采用的就是基于带遗忘因子的递推最小二乘法设计的自校正PID控制134 系统硬件的结构设计4.1 系统硬件的结构该系统以单片机为核心，包括模拟信号处理电路(热电偶放大电路、线性光电隔离电路、多路模拟开关、环境温度测量电路)、A/D 转换电路、RS232 接口电路、输出光电隔离电路、控制输出电路以及电源电路等几部分组成供电电源温度检测PT100线性光电隔离多路模拟开 关冷端温度测量电路A/D转换单片机AT89s52RS232计算机光电隔离控制输出变频器图4.1自校正PID实验锅炉温度控制系统控制图实验锅炉炉温的控制过程如下：测温元件将检测到的温度信号经过放大、隔断处理后，送到多路选择器，多路选择器在单片机 CPU 的控制下对这组信号进行选择，在某一时刻被选择的信号送到 A/D转换器，转换成数字量送给单片机，单片机对 A/D 结果处理(包括冷端补偿计算)，得出加热炉温度测量值，再将温度测量值与给定温度值进行 PID运算，产生控制脉宽输出，来控制加热炉的加热过程。4 .2 自校正PID实验锅炉控制系统原理被控对象为220V交流电阻加热实验锅炉。实验锅炉与加热器对象数学模型，水泵是用来使水活动这样可以均匀加热防止加热不均匀测量而产生误差。5# 锅炉调压器定值图4.21 锅炉温度定值控制TT101记录锅炉为大滞后系统，其一阶飞升曲线如下图：图4.22实验锅炉炉温测定飞升曲线设阶跃响应曲线为S形的单调曲线，如图4 .22所示，设阶跃输入幅值为，则增益可按下式求取： 时间常数及延迟时间可用作图法确定：在图4.22响应曲线的拐点做切线，切线与时间轴交与点，而与响应曲线稳态值交与点，则对应延迟时间，对应时间常数。响应传递函数为：测得实验锅炉为具有纯滞后的一阶惯性环节，即14：带零阶保持器的广义对象脉冲传递函数为 (4-1)取，则。将实测被控对象参数带入式(4-21)中得：从极点配置的观点出发，以典型的二阶系统闭环传递函数的离散特征多项式15 (4-2)为无阻尼自然震荡角频率，阻尼比。对于式(4-2)，当二阶系统最佳阻尼比时，在单位阶跃作用下的超调量，相角稳定裕量，它为二阶最佳动态响应模型。采样周期的关系 ： 当。当式(4 -22)的期望特征多项式为：建立期望的传递函数： (4-3)4.3 自校正控制算法设计实验锅炉的模型公式由得闭环系统希望的传递函数为：在递推估计调节参数时，设置遗忘因子。4.4 自校正PID算法设计流程图结束绘制出如 r y Y的图形，输出 和P的最终值YN图4.4 算法设计流程图5 Matlab 仿真及结果分析设输是阶跃信号输入，如图5.1图5.1系统输入信号图5.2系统期望输出图5.3系统经自校正PID调节后输出(初值为2，2，2)图5.4系统经常规PID调节输出的波形(初值为2，2，2)在matlab环境下可得上述闭环系统的期望输出和经过自校正PID调节后的系统实际输出。由图(5.3)可得自校正能够对期望值进行较好的跟踪。最终： 为了更好的进行对比，在matlab环境下进行常规PID调节后的输出为图(5.4)很明显，在在相同的初始值下，经常规PID调节的系统不稳定了把初始值设为则常规PID调节实验锅炉后的如图(5.5)。图5.5初始值设为常规PID调节的图形由图可知自校正有自动调节PID参数的能力，能自适应被控制的对像。为了验证自校正PID控制器的自动跟随能力，设系统输入该为方波信号为一个周期200s，幅值为1的矩形波，见下图：图5.6系统输入波形图5.7期望输出的波形图5.8系统经自校正PID输出的波形图5.9常规PID锅炉炉温控制系统输出波形由图(5.6)(5.7)(5.8)可得自校正PID控制的锅炉可以很好的跟随输入，而常规PID即使应用了自校正得到的值，虽然稳定了但从波形可以看出跟随系统输入的效果很差。为了检验是否自校正PID的初始值可以任意设定，改变初始值用图5.10系统自校正PID控制实验锅炉系统输出波形(初始值为) 由图(5.10)可以看出虽然在刚开始系统有很大的震荡，但在有限个震荡后系统输出为1跟随了系统的输入，最终的 说明自校正PID控制实现了实验锅炉温度的是最优控制。6 结论本论文认真地研究和分析了实验锅炉温度控制系统，利用自校正PID调节器，对递推最小二乘算法进行了推导。针对给实验锅炉炉温传递函数和期望传递函数的闭环系统，编制了自校正PID调节器参数调整过程的流程图，并利用MATLAB程序进行了编程和仿真实验。将其与普通自校正调节器进行了对比实验，验证了自校正PID调节器的优越性。并针对系统可能产生的一些情况变化，比如初值的改变、系统模型的改变、输入的改变，进行了仿真实验。通过实验证明，自校正PID调节器控制炉温温度和普通PID调节器控制炉温稳定相比，对于自校正系统具有很好的适应能力。参考文献1 王顺晃. 非线性PID算法及其在电阻加热炉集散控制系统中的应用J. 自动化学报, 1995. (6):675681 .2 闰秀英, 任庆昌, 盂庆龙. 一种自校正PID控制器设计与仿真研究J. 系统仿真报, 2006, (2): 753756.3 邓娟, 丛爽. 控制工程中几种自校正控制的应用及其性能对比研究J.电子学报, 2008, (3): 1418.4 张显库,贾欣乐.闭环增益成形控制新算法及其应用J.电子学报,1999,27(11): 1331355 Tjokro S, Shah S L. Adaptive PID ControlJ. Proceeding of the 1987 American Control Confernce, 1985(2): 15281534.6 刘伯春. 离散系统采样周期的选择J化工自动化及仪表, 1989, (1): 2833.7 Astrom K. J. Witten mark B. Adaptive Control, Addison-Wesley,1989李清泉等译, 自适应控制M. 北京: 科学出版社, 1992.8 刘东升. PID自校正调节器(J). 现代防御技术, 1991(4): 5659.9 夏天长(美)著, 熊光愣, 李芳芸译. 系统辨识最小二乘法M.北京: 清华大学出版社, 1983.2123.10 Hsia T.CSystem IdentificationLeast Square MethodsM.Le Xing Mass: Le XingTong Bokks, 1977.11 Shankar Sasty , Marc Bodeson. Adaptive ControlStabilityM.Convergence and Robustness, Prentice Hall.1989 .12 李培玉,谭大鹏, 张鹏, 倪笑裴, 林搏宇. 自校正PID在加弹机温控器中的应用研究. 浙江大学J. 电子学报, 2006, (2): 2548.13 周忠海, 惠超, 葬鹤超, 成园林, 张照文. 基于PID控制电阻加热炉炉温控制系统. 山东科学J , 2006, (10): 4547.14 唐永哲.自校正PID调节器的极点配置设计法J. 工业仪表与自动化装置, 1995, (1): 35 15 李开霞, 张金波, 曹爱华.一种实用的自校正PID控制器设计与仿真J.2008, (9): 167170.附 录自校正PID锅炉控制程序clear all;close all;e_0=0;e_1=0;u_0=0;u_1=0;u_01=0;y_0=0;y_1=0.5;y_2=0;y_3=0;Y_0=0;Y_1=0;Y_2=0;A_0=2,2, 2;A_1=A_0;P_0=106*eye(3,3);P_1=P_0;c=0.995;F_0=0;int k;for k = 1:1:100time(k)=k;%r(k)=sign(sin(0.01*pi*k);r(k)=sign(1)y(k)=-0.7987*y_0+0.5636*u_0;x(1,1)=1*y(k)+0.7987*y_0-0.5636*y_1;x(2,1)=1*y_0+0.7987*y_1-0.5636*y_2;x(3,1)=1*y_1+0.7987*y_2-0.5636*y_3;e(k)=r(k)-y(k);Y(k)=0 .7812*Y_0-0.2882*Y_1+0.506*r(k);E(k)=Y(k)-y(k);F(k)=abs(e(k)+e(k)2;F(k)=F(k)+F_0;F_0=F(k);if(abs(E(k)=0)u(k)=A_0(1)*e(k)+A_0(2)*e_0+A_0(3)*e_1+u_1;P_1=P_0/c-P_0*x*(c+x*P_0*x)(-1)*x*P_0/c;%v(k)=0.702*u_1-0 .0101*u_0-0.0601*u_01;v(k)=(0 .5636*u_0-0 .5636*u_01); %xiugaiA_1=A_0+P_0*x*(c+x*P_0*x)(-1)*(v(k)-(x)*A_0);A_0=A_1;P_0=P_1;endy_2=y_1;y_1=y_0;y_0=y(k);Y_1=Y_0;Y_0=Y(k);u_01=u_0;u_0=u_1;u_1=u(k);e_1=e_0;endfigure(1);plot(time, r, k);xlabel(time(s);ylabel(r);axis(-5 500,0 1.5);figure(2);plot(time, Y, k);xlabel(time(s);ylabel(Y);axis(-5 500，0 1.5);figure(3);plot(time，y，k);xlabel(time(s);ylabel(Y);axis(-5 500，0 1 .5);A_1，P_1，F(k)常规PID校正程序clear all;close all;y_1=0;y_0=0;u_1=0;u_2=0;e_1=0;e_2=0;x=0.42，-0.0707，-0.2419;int u100;int e1000;int r1000;for k=1:1:1000; time(k)=k; r(k)=sign(sin(0 .01*pi*k); %r(k)=sign(k); y(k)=-0.7987*y_0+0.5636*u_2; e(k)=r(k)-y(k); u(k)=(u_1+x(1)*e(k)+x(2)*e_1+x(3)*e_2); u_2=u_1;u_1=u(k); e_2=e_1;e_1=e(k); y_0=y(k);endplot(time, y, k)xlabel(time(s), ylabel(y)figure(2);plot(time, y, k);xlabel(time(s);ylabel(Y); axis(-5 100, 0 1.5);毕业设计（论文）管理规定及相关表格汇编目 录本科生毕业设计（论文）管理规定1毕业论文（设计）工作实施细则（草案）2优秀毕业论文评选办法 7毕业论文的写作与排版规范10毕业设计（论文）写作模板14毕业设计（论文）答辩提问记录表 21毕业设计 第 周 工作总结 22毕业设计（论文）过程跟踪表 23毕业设计（论文）课题申请表 24毕业设计（论文）任务书 25毕业设计（论文）开题报告 26毕业设计工作中期检查 27毕业设计工作中期检查28大学本科生毕业设计（论文）管理规定为提高本科生毕业设计（论文）质量，加强毕业生毕业设计管理工作，特制定如下规定：1. 本科毕业设计工作从第七学期（四年级第一学期）考试周前，完成教师选题、师生见面及指导教师向学生下达任务书（见附表 一）。毕业设计从第八学期开学正式开始，十六周内完成。2. 毕业设计开始时，教师必须填写毕业设计课题申请表、并由系汇总后交院毕业设计领导小组审核。（见附表）。3. 指导教师必须定期对学生辅导（每周至少两次），并将确定的周辅导时间上报教学办，教学办将定期和不定期进行检查。4. 指导教师自行组织学生的开题报告，并填写开题报告表（见附表二）5. 加强中期的监督与检查，第一次检查时间定在第五周，检查内容含学生资料阅读，方案论证情况。第二次检查定在第十周，形式为：由指导教师组织，院毕业设计领导小组成员旁听，学生口头汇报。主要检查学生的阶段性结果。两次检查学生均应填写中期检查报告表格（见附表三、四），指导教师或系对检查情况写出书面评语。6. 学生需上交一篇与本专业有关不少于5000字的外文文章翻译，并将原文和译文用A4打印加自设计封皮一起装订成册。7. 论文答辩前，指导教师必须给出评语及评分，然后由各系主任指定其他评阅人对论文进行评阅及评分（请参考毕业论文书面成绩评分表规定的评分标准评分）。8. 论文答辩由各系组织，并成立答辩委员会对所有学生论文进行答辩。答辩必须给出答辩委员会意见及评分并填写相应表格。答辩委员会成员必须按评分表中各项指标标准进行评分，然后由各成员给出的分数计算出得分。9. 毕业设计（论文）最终成绩计算：指导教师评分x 30% + 评阅人评分x 30% + 答辩委员会评分 x 40% = 总分10毕业设计（论文）提交的文档及装订要求a) 毕业论文一份（包含封皮、目录、中英文摘要、内容及参考文献）b) 不少于5000汉字的科技翻译资料一份（并将原文和译文用A4打印加自设计封皮一起装订成册）c) 毕业论文简介（A4纸打印12页）（包含题目、专业、年级、姓名、指导教师、毕业论文所做的工作、解决的问题、创新之处等）d) 毕业设计任务书（加在毕业论文目录之前和毕业论文一起装订，）e) 开题报告（加在毕业论文目录之前和毕业论文一起装订，）f) 毕业设计工作中期检查表、注（4）（5）（6）加在毕业论文目录之前和毕业论文一起装订。11毕业设计（论文）随正式论文一律附交电子文档（光盘或软盘）。被推荐参加学校优秀毕业设计（论文）3%评选的，还要另附（符合发表格式要求的）修改为5000字左右的电子文档。 2005年12月修订大学 毕业设计（论文）工作实施细则（试行）根据毕业设计（论文）工作暂行规定，为了进一步规范毕业设计（论文）工作的各个环节，结合我院实际情况，经院教学指导委员会讨论，院长办公联席会议通过，对我院毕业设计（论文）工作特制订以下实施细则。一、指导教师资格毕业设计（论文）的指导教师原则上应由各专业具有中级(含中级)以上专业技术职务的教师担任。也可以聘任具有相应职称、学术水平高的工程技术人员担任。指导教师名单由各系系主任核定，报分管教学的副院长批准后实施。二、课题的命题过程1命题的基本原则是：（1）目标原则：命题必须符合专业培养目标，具有一定的深度、广度与新颖性，能使学生得到比较全面的训练。（2）价值原则：命题应具有一定的学术研究或应用开发价值，尽可能做到理论与实际相结合，体现学科的发展性与应用性。（3）能力原则：命题的难度应考虑学生完成命题的条件与能力，因材施教。（4）兴趣原则：应考虑学生对命题的兴趣与探索研究的积极性。（5）多样化原则：命题应注意题目的不断更新与类型的多样化，同一指导教师的设计（论文）课题选题和近三届毕业设计课题的重复率不高于50%，题目数量应满足一人一题的基本要求。2课题的审核由各专业指导教师确定后的课题，填写毕业设计课题申请表。经系主任汇总后，报院毕业设计领导小组审核，由分管教学的副院长批准。3确定课题的时间要求一般在第七学期第15周前完成命题，对于经批准认可的命题，在第16周内向学生公布。指导教师在实际指导过程中，一般不应随意更改命题，确实需要更改命题时，应办理相关手续。三、选题与任务书1学生在教师指导下从公布的命题范围内选题。2学生自拟题目须经系主任审阅后指定指导教师，并报分管教学副院长批准。3选题原则上要求一人一题。4学生选题一般在第七学期的第17周前完成。5根据学生选题情况，填写大学毕业论文（设计）任务书，其中“设计所需的工作条件”应与课题的实际需要相符；“设计（论文）任务及目标”应指明毕业设计的详细任务（所需完成的工作、达到的设计要求等），不能简单了事。四、开题1指导学生书写开题报告，开题报告主要内容包括：课题研究的意义和背景、课题研究相关资料的研究情况（说明课题研究的历史、现状和发展趋势）、对课题研究的基本构思。2指导学生填写大学毕业论文（设计）开题报告，并要求学生在第八学期的第2周前提交不少于1000字的开题报告，基本内容应包括课题研究的意义和背景、课题研究相关资料的研究情况（说明课题研究的历史、现状和发展趋势）、对课题研究的基本构思、阶段性工作计划、参考文献等。参考文献的列写格式要符合要求。3第八学期的第3周前完成开题审核，通过开题审核的学生，一般不应更换选题；开题未获通过的学生，应及时根据开题指导小组的意见更换选题，并于开题后的2周内做好准备，向开题指导小组第二次提出开题申请；对再次开题未获通过的学生，取消其本次毕业设计（论文）资格，延至下届进行。指导教师要重视对开题未获通过的学生进行指导。4开题报告中，“指导教师意见”应对学生开题情况作出基本评价以及对学生毕业设计（论文）研究方案或毕业设计（论文）写作提纲提出改进意见，同时还应对学生是否可以进入研究和设计（