基于MATLAB的PID控制器设计毕业论文正式

1、目 录 摘 要.错误！未定义书签。错误！未定义书签。 abstract.错误！未定义书签。错误！未定义书签。 第一章 绪 论.1 1.1 课题来源及 pid 控制简介.1 1.1.1 课题的来源和意义 .1 1.1.2 pid 控制简介.1 1.2 国内外研究现状及 matlab 简介.3 第二章 控制系统及 pid 调节.5 2.1 控制系统构成.5 2.2 pid 控制.5 2.2.1 比例、积分、微分 .5 2.2.2 、控制 .7 第三章 系统辨识.9 3.1 系统辨识.9 3.2 系统特性图.10 3.3 系统辨识方法.11 第四章 pid 最佳调整法与系统仿真.13 4.1 pid

2、 参数整定法概述.13 4.2 针对无转移函数的 pid 调整法.14 4.2.1 relay feedback 调整法.14 4.2.2 relay feedback 在计算机做仿真.14 4.2.3 在线调整法 .16 4.2.4 在线调整法在计算机做仿真 .17 4.3 针对有转移函数的 pid 调整方法.18 4.3.1 系统辨识法 .18 4.3.2 波德图法及根轨迹法 .20 第五章 油冷却机系统的 pid 控制器设计.21 5.1 油冷却机系统.21 5.1.1 油冷却机 .21 5.1.2 感测与转换器 .22 5.1.3 控制组件 .23 5.2 油冷却机系统之系统辨识.24

3、 5.3 油冷却机系统的 pid 参数整定.26 结 论.32 致 谢.33 参考文献.34 青岛农业大学机电工程学院本科毕业设计 1 第一章 绪 论 1.1 课题来源及 pid 控制简介 1.1.1 课题的来源和意义 任何闭环的控制系统都有它固有的特性，可以有很多种数学形式来描述它，如微分 方程、传递函数、状态空间方程等。但这样的系统如果不做任何的系统改造很难达到最 佳的控制效果，比如快速性稳定性准确性等。为了达到最佳的控制效果，我们在闭环系 统的中间加入 pid 控制器并通过调整 pid 参数来改造系统的结构特性，使其达到理想的 控制效果。 1.1.2 pid 控制简介 当今的自动控制技术

4、都是基于反馈的概念。反馈理论的要素包括三个部分：测量、 比较和执行。测量关心的变量，与期望值相比较，用这个误差纠正调节控制系统的响应。 这个理论和应用自动控制的关键是，做出正确的测量和比较后，如何才能更好地纠 正系统，pid （比例 - 积分 - 微分）控制器作为最早实用化的控制器已有 50 多年历史， 现在仍然是应用最广泛的工业控制器。 pid 控制器简单易懂，使用中不需精确的系统模 型等先决条件，因而成为应用最为广泛的控制器。 pid 控制器由比例单元（ p ） 、积分单元（ i ）和微分单元（ d ）组成。其输入 e (t) 与输出 u (t) 的关系为公式（1-1） 公式（1-1） 因

5、此它的传递函数为公式（1-2） 公式（1-2） 比例调节作用：是按比例反应系统的偏差,系统一旦出现了偏差,比例调节立即产生调 节作用用以减少偏差。比例作用大,可以加快调节,减少误差,但是过大的比例,使系统的稳 定性下降,甚至造成系统的不稳定。 积分调节作用：是使系统消除稳态误差,提高无差度。因为有误差,积分调节就进行,直 至无差,积分调节停止,积分调节输出一个常值。积分作用的强弱取决与积分时间常数 ti，ti 越小,积分作用就越强。反之 ti 大则积分作用弱,加入积分调节可使系统稳定性下降,动 青岛农业大学机电工程学院本科毕业设计 2 态响应变慢。积分作用常与另两种调节规律结合,组成 pi 调

6、节器或 pid 调节器。 微分调节作用：微分作用反映系统偏差信号的变化率,具有预见性,能预见偏差变化的 趋势,因此能产生超前的控制作用,在偏差还没有形成之前,已被微分调节作用消除。因此, 可以改善系统的动态性能。在微分时间选择合适情况下,可以减少超调,减少调节时间。微 分作用对噪声干扰有放大作用,因此过强的加微分调节,对系统抗干扰不利。此外,微分反应 的是变化率,而当输入没有变化时,微分作用输出为零。微分作用不能单独使用,需要与另外 两种调节规律相结合,组成 pd 或 pid 控制器。 pid 控制器由于用途广泛、使用灵活，已有系列化产品，使用中只需设定三个参数（ kp ， ki 和 kd ）

7、即可。在很多情况下，并不一定需要全部三个单元，可以取其中的 一到两个单元，但比例控制单元是必不可少的。 首先，pid 应用范围广。虽然很多控制过程是非线性或时变的，但通过对其简化可以 变成基本线性和动态特性不随时间变化的系统，这样 pid 就可控制了。 其次，pid 参数较易整定。也就是，pid 参数 kp，ki 和 kd 可以根据过程的动态特 性及时整定。如果过程的动态特性变化，例如可能由负载的变化引起系统动态特性变化， pid 参数就可以重新整定。 第三，pid 控制器在实践中也不断的得到改进，下面两个改进的例子，在工厂，总是 能看到许多回路都处于手动状态，原因是很难让过程在“自动”模式下

8、平稳工作。由于这些 不足，采用 pid 的工业控制系统总是受产品质量、安全、产量和能源浪费等问题的困扰。 pid 参数自整定就是为了处理 pid 参数整定这个问题而产生的。现在，自动整定或自身整 定的 pid 控制器已是商业单回路控制器和分散控制系统的一个标准。 在一些情况下针对特定的系统设计的 pid 控制器控制得很好，但它们仍存在一些问 题需要解决：如果自整定要以模型为基础，为了 pid 参数的重新整定在线寻找和保持好 过程模型是较难的。闭环工作时，要求在过程中插入一个测试信号。这个方法会引起扰 动，所以基于模型的 pid 参数自整定在工业应用不是太好。 如果自整定是基于控制律的，经常难以

9、把由负载干扰引起的影响和过程动态特性变 化引起的影响区分开来，因此受到干扰的影响控制器会产生超调，产生一个不必要的自 适应转换。另外，由于基于控制律的系统没有成熟的稳定性分析方法，参数整定可靠与 否存在很多问题。 因此，许多自身整定参数的 pid 控制器经常工作在自动整定模式而不是连续的自身 整定模式。自动整定通常是指根据开环状态确定的简单过程模型自动计算 pid 参数。 但仍不可否认 pid 也有其固有的缺点： pid 在控制非线性、时变、耦合及参数和结构不确定的复杂过程时，工作地不是太 青岛农业大学机电工程学院本科毕业设计 3 好。最重要的是，如果 pid 控制器不能控制复杂过程，无论怎么

10、调参数都没用。 虽然有这些缺点，pid 控制器是最简单的有时却是最好的控制器。 1.2 国内外研究现状及 matlab 简介 pid 控制中最重要的是对其参数的控制，所以当今国内外 pid 控制技术的研究主要是 围绕如何对其参数整定进行的。 自 ziegler 和 nichols 提出 pid 参数整定方法起，有许多技术已经被用于 pid 控制器 的手动和自动整定.根据发展阶段的划分，可分为常规 pid 参数整定方法及智能 pid 参数 整定方法；按照被控对象个数来划分，可分为单变量 pid 参数整定方法及多变量 pid 参 数整定方法，前者包括现有大多数整定方法，后者是最近研究的热点及难点；

11、按控制量 的组合形式来划分，可分为线性 pid 参数整定方法及非线性 pid 参数整定方法，前者用 于经典 pid 调节器，后者用于由非线性跟踪-微分器和非线性组合方式生成的非线性 pid 控制器。 astrom 在 1988 年美国控制会议（acc）上作的面向智能控制 2的大会报告概述 了结合于新一代工业控制器中的两种控制思想自整定和自适应，为智能 pid 控制的 发展奠定了基础。他认为自整定控制器和自适应控制器能视为一个有经验的仪表工程师 的整定经验的自动化，在文3中继续阐述了这种思想,认为自整定调节器包含从实验中 提取过程动态特性的方法及控制设计方法，并可能决定何时使用 pi 或 pid

12、 控制，即自整 定调节器应具有推理能力。自适应 pid 的应用途径的不断扩大使得对其整定方法的应用 研究变得日益重要。目前，在众多的整定方法中，主要有两种方法在实际工业过程中应 用较好，一种是由福克斯波罗（foxboro）公司推出的基于模式识别的参数整定方法（基 于规则） ，另一种是基于继电反馈的参数整定方法（基于模型）.前者主要应用于 foxboro 的单回路 exact 控制器及其分散控制系统 i/a series 的 pide 功能块，其原理基于 bristol 在模式识别方面的早期工作11。后者的应用实例较多，这类控制器现在包括自整 定、增益计划设定及反馈和前馈增益的连续自适应等功能.

13、这些技术极大地简化了 pid 控 制器的使用，显着改进了它的性能，它们被统称为自适应智能控制技术。 4 自适应技术中最主要的是自整定。按工作机理划分，自整定方法能被分为两类：基 于模型的自整定方法和基于规则的自整定方法。 4 在基于模型的自整定方法中，可以通过暂态响应实验、参数估计及频率响应实验来 获得过程模型。 在基于规则的自整定方法中，不用获得过程实验模型，整定基于类似有经验的操作 者手动整定的规则。 青岛农业大学机电工程学院本科毕业设计 4 为了满足不同系统的要求，针对多变量和非线形的系统还分别采用了多变量 pid 参 数整定方法和非线性 pid 参数整定方法。 pid 控制算法是迄今为

14、止最通用的控制策略.有许多不同的方法以确定合适的控制器 参数.这些方法区分于复杂性、灵活性及使用的过程知识量。一个好的整定方法应该基于 合理地考虑以下特性的折衷：负载干扰衰减，测量噪声效果，过程变化的鲁棒性，设定 值变化的响应，所需模型，计算要求等.我们需要简单、直观、易用的方法，它们需要较 少的信息，并能够给出合适的性能。我们也需要那些尽管需要更多的信息及计算量，但 能给出较好性能的较复杂的方法。 从目前 pid 参数整定方法的研究和应用现状来看，以下几个方面将是今后一段时间 内研究和实践的重点。 4 对于单输入单输出被控对象，需要研究针对不稳定对象或被控过程存在较大干扰 情况下的 pid

15、参数整定方法，使其在初始化、抗干扰和鲁棒性能方面进一步增强，使用 最少量的过程信息及较简单的操作就能较好地完成整定。 对于多入多出被控对象，需要研究针对具有显着耦合的多变量过程的多变量 pid 参数整定方法，进一步完善分散继电反馈方法，尽可能减少所需先验信息量，使其易于 在线整定。 4 智能 pid 控制技术有待进一步研究，将自适应、自整定和增益计划设定有机结合， 使其具有自动诊断功能；结合专家经验知识、直觉推理逻辑等专家系统思想和方法对原 有 pid 控制器设计思想及整定方法进行改进；将预测控制、模糊控制和 pid 控制相结合， 进一步提高控制系统性能，都是智能 pid 控制发展的极有前途的

16、方向。 4 matrix laboratory(缩写为 mat lab)软件包，是一种功能强、效率高、便于进行科学和 工程计算的交互式软件包。其中包括:一般数值分析、矩阵运算、数字信号处理、建模和 系统控制和优化等应用程序，并将应用程序和图形集于便于使用的集成环境中。在此环 境下所解问题的 mat lab 语言表述形式和其数学表达形式相同，不需要按传统的方法编程 并能够进行高效率和富有创造性的计算，同时提供了与其它高级语言的接口，是科学研 究和工程应用必备的工具。目前，在控制界、图像信号处理、生物医学工程等领域得到 广泛的应用。本论文设计中 pid 参数的整定用到的是 mat lab 中的 s

17、imulink,它是一个 强大的软件包 ，在液压系统仿真中只需要做数学模型的推导工作。用 simulink 对设 计好的系统进行仿真，可以预知效果，检验设计的正确性，为设计人员提供参考。其仿 真结果是否可用，取决于数学模型正确与否，因此要注意模型的合理及输入系统的参数 值要准确。 8 青岛农业大学机电工程学院本科毕业设计 5 第二章 控制系统及 pid 调节 2.1 控制系统构成 对控制对象的工作状态能进行自动控制的系统称为自动控制系统，一般由控制器与 控制对象组成，控制方式可分为连续控制与反馈控制，即一般所称，开回路与闭回路控 制。 连续控制系统的输出量对系统的控制作用没有任何影响，也就是说

18、，控制端与控制 对象为单向作用，这样的系统亦称开回路系统。 反馈控制是指将所要求的设定值与系统的输出值做比较，求其偏差量，利用这偏差 量将系统输出值使其与设定值调为一致。 反馈控制系统方块图一般如图 2-1 所示： 比较组件控制器被控对象 感测与转换 图 2-1 反馈控制系统方块图 2.2 pid 控制 将感测与转换器输出的讯号与设定值做比较，用输出信号源(2-10v 或 4-20ma)去控制 最终控制组件。在工程实际中，应用最为广泛的调节器控制规律为比例积分微分控制， 简称 pid 控制，又称 pid 调节。pid 控制器问世至今已有近 60 年的历史了，它以其结构 简单、稳定性好、工作可靠

19、、调整方便而成为工业控制主要和可靠的技术工具。当被控 对象的结构和参数不能完全掌握，或得不到精确的数学模型时，控制理论的其它设计技 术难以使用，系统的控制器的结构和参数必须依靠经验和现场调试来确定，这时应用 pid 控制技术最为方便。即当我们不完全了解一个系统和被控对象或不能通过有效的测量 青岛农业大学机电工程学院本科毕业设计 6 手段来获得系统的参数的时候，便最适合用 pid 控制技术。 2.2.1 比例、积分、微分 1. 比例 r1 r2 vi vo 图 2-2 比例电路 1 2 )( )( r r vi vo t t 公式（2-1） )( 1 2 )(tt vi r r vo 2. 积分

20、器 r1 vi vo 1/sc 图 2-3 积分电路 scrscrr sc vi vo t t111 1 111)( )( vi scr vo t 11 1 )( 公式（2-2）dtvi cr vo t 1 )( 1 vo r2 1/sc vi 青岛农业大学机电工程学院本科毕业设计 7 图 2-4 微分电路 3. 微分器 scr sc r vi vo t t 2 2 )( )( 1 )(2)(tt svicrvo (式 2-3) dt dvi crvo t2)( 实际中也有 pi 和 pd 控制器。pid 控制器就是根据系统的误差利用比例积分微分计 算出控制量，控制器输出和控制器输入（误差）之

21、间的关系在时域中如公式（2-4）和 （2-5）: u(t)=kp(e(t)+td+) 公式（2-4） dt tde )( dtte ti )( 1 u(s)=+e(s) 公式（2-5） p k s ki skd 公式中 u(s)和 e(s)分别为 u(t)和 e(t)的拉氏变换，其中、 p d d k k t i p k k ti p k 、分别为控制器的比例、积分、微分系数14 i k d k 2.2.2 、控制 比例（比例（p）控制）控制 比例控制是一种最简单的控制方式。其控制器的输出与输入误差讯号成比例关系。 当仅有比例控制时系统输出存在稳态误差（steady-state error）

22、。 积分（积分（i）控制）控制 在积分控制中，控制器的输出与输入误差讯号的积分成正比关系。 对一个自动控制系统，如果在进入稳态后存在稳态误差，则称这个控制系统是有稳 态误差的或简称有差系统（system with steady-state error） 。为了消除稳态误差，在控制器 中必须引入“积分项” 。积分项对误差取关于时间的积分，随着时间的增加，积分项会增 大。这样，即便误差很小，积分项也会随着时间的增加而加大，它推动控制器的输出增 大使稳态误差进一步减小，直到等于零。 因此，比例+积分(pi)控制器，可以使系统在进入稳态后无稳态误差。 微分（微分（d）控制）控制 青岛农业大学机电工程学

23、院本科毕业设计 8 在微分控制中，控制器的输出与输入误差讯号的微分（即误差的变化率）成正比关 系。 自动控制系统在克服误差的调节过程中可能会出现振荡甚至失稳。其原因是由于存 在有较大惯性的组件（环节）和（或）有滞后(delay)的组件，使力图克服误差的作用，其 变化总是落后于误差的变化。解决的办法是使克服误差的作用的变化要有些“超前” ，即 在误差接近零时，克服误差的作用就应该是零。这就是说，在控制器中仅引入“比例” 项往往是不够的，比例项的作用仅是放大误差的幅值，而目前需要增加的是“微分项” ， 它能预测误差变化的趋势，这样，具有比例+微分的控制器，就能够提前使克服误差的控 制作用等于零，甚

24、至为负值，从而避免了被控量的严重地冲过头。 所以对有较大惯性和（或）滞后的被控对象，比例+微分(pd)的控制器能改善系统在 调节过程中的动态特性。 12 青岛农业大学机电工程学院本科毕业设计 9 第三章系统辨识 3.1 系统辨识 (1) 所谓系统辨识即是在不知道系统转移函数时，根据系统特性辨识出来。 (2) 若被控对象的数学模式相当线性(linear)，且各项参数都可知道，则可用控制理论 来 设计 pid 控制器的系数大小。但实际的被控对象往往是非线性系统，且系统复杂， 难以精确地用数学式表达。所以工业上设计 pid 控制器时，常常使用实验方法而较少用 理论来设计。调整 pid 控制器的方法中

25、，最有名的是 ziegler-nichols 所提出的二个调整法 则。这个调整法测是基于带有延迟的一阶传递函数模型提出的，这种对象模型可以表示 为 1 )( ts ke sg ls 公式（3-1） 在实际的过程控制系统中，有大量的对象模型可以近似的由这样的一阶模型来表示， 如果不能物理的建立起系统的模型，我们还可以由实验提取相应的模型参数5。 (3) 将大小为 1 的阶跃信号加到被控对象如图 3-1 所示： 图 3-1 将阶跃信号加到被控对象 对大多数的被控对象，若输入为阶跃信号，则其输出 c(t)大多为 s 状曲线，如下图 3-2 所示。这个 s 状曲线称之为过程反应曲线（process r

26、eaction curve） 。 l 0.632k k t t” t c(t) 图 3-2 被控对象的阶跃响应图 青岛农业大学机电工程学院本科毕业设计 10 (4) 系统转移函数 空调方面： 图 3-3 空调系统示意图 图 3-4 空调系统方块图 由图 3-3 及图 3-4 可得知此系统的转移函数推导如下： qcqoq dt dt c r t q cst r t q 公式（3-2） 11 s r rcs r q t 3.2 系统特性图 (1) 系统为制热 使用最大信号去控制系统，直到稳定之后，也就是热到达无法再上升时，此时系统 特性就会出现，如下图 3-5 所示： 青岛农业大学机电工程学院本科

27、毕业设计 11 10 50 图 3-5 系统制热的特性图 (2) 系统为制冷 使用最大信号去控制系统，直到稳定之后，也就是冷到达无法再下降时，此时系统 特性就会出现，如下图 3-6 所示： 10 3010 30 图 3-6 系统制冷的特性图 3.3 系统辨识方法 （1）一阶系统带有延迟特性 l 0.632k k a t t” t c(t) 图 3-7 一阶系统带有延迟特性图 此切线 为最大斜率 青岛农业大学机电工程学院本科毕业设计 12 一阶系统加一个传递来近似被控对象，则其近似转移函数如公式 3-3 所示： 1ts k ls e 公式（3-3） 1 )( ts ke sg ls 其中 k、t

28、、l 可由上图 3-7 求得。 k：稳态时的大小。 t:时间常数。注：系统越大，时间常数越大。 l：延迟时间。 (2) k、t、l 的求法： k：如上图 3-3.1 所示，k 值相当于 c(t)在稳态时的大小。 t 与 l：求 t 及 l 必须在 s 形状曲线划一条切线(最大斜率)，画出切线之后，t 及 l 值 可以直接从图上得知。t 及 l 值与 c(t)及切线的关系如上图 3-7 所示。 青岛农业大学机电工程学院本科毕业设计 13 第四章pid 最佳调整法与系统仿真 4.1 pid 参数整定法概述 1.pid 参数整定方法 （1） relay feedback ：利用 relay 的 on

29、-off 控制方式，让系统产生一定的周期震 荡，再用 ziegler-nichols 调整法则去把 pid 值求出来。 （2） 在线调整：实际系统中在 pid 控制器输出电流信号装设电流表，调 p 值观察 电流表是否有一定的周期在动作，利用 ziegler-nichols 把 pid 求出来，pid 值求法与 relay feedback 一样。 （3） 波德图&跟轨迹：在 matlab 里的 simulink 绘出反馈方块图。转移函数在用 系统辨识方法辨识出来，之后输入指令算出 pid 值。13 2.pid 调整方式 pid 调整方式 有转移函数无转移函数 系统辨识法 波德图 根轨迹 rel

30、ay feedback在线调整 图 4-1 pid 调整方式 如上描述之 pid 调整方式分为有转函数和无转移函数，一般系统因为不知转移函数， 所以调 pid 值都会从 relay feedback 和在线调整去着手。波德图及根轨迹则相反，一定要 有转移函数才能去求 pid 值，那这技巧就在于要用系统辨识方法，辨识出转移函数出来， 再用 matlab 里的 simulink 画出反馈方块图，调出 pid 值。 15 所以整理出来，调 pid 值的方法有在线调整法、relay feedback、波德图法、根轨迹 法。前提是要由系统辨识出转移函数才可以使用波德图法和根轨迹法，如下图 4-2 所示。

31、 青岛农业大学机电工程学院本科毕业设计 14 4.2针对无转移函数的 pid 调整法 在一般实际系统中，往往因为过程系统转移函数要找出，之后再利用系统仿真找出 pid 值，但是也有不需要找出转移函数也可调出 pid 值的方法，以下一一介绍。 4.2.1relay feedback 调整法 图 4-3 relay feedback 调整法 如上图 4-3 所示，将 pid 控制器改成 relay，利用 relay 的 on-off 控制，将系统扰动， 可得到该系统于稳定状态时的震荡周期及临界增益（tu及u） ，在用下表 4-4 的 ziegler- nichols 第一个调整法则建议 pid 调

32、整值，即可算出该系统之p、ti、tv之值。 表 4-4 ziegler-nichols 第一个调整法则建议 pid 调整值9 4.2.2 relay feedback 在计算机做仿真 controller p k i t d t p0.5u pi0.45 u 0.83tu pid0.6u0.5tu0.125tu 图 4-2 由系统辨识法辨识出转移函数 青岛农业大学机电工程学院本科毕业设计 15 step 1： 以 matl ab 里的 simulink 绘出反馈方块，如下图 4-5 所示。 图 4-5 simulink 绘出的反馈方块图 stepstep 2 2：让 relay 做 on-of

33、f 动作，将系统扰动（on-off 动作，将以 做模拟） ， 如下图 4-6 所示。 图 4-6 step 3：即可得到系统的特性曲线，如下图 4-7 所示。 青岛农业大学机电工程学院本科毕业设计 16 图 4-7 系统震荡特性曲线 step 4：取得 tu 及 a，带入公式 3-1，计算出u。 以下为 relay feedback 临界震荡增益求法 公式（4-1） a d ku 4 ：振幅大小：电压值 4.2.3在线调整法 图 4 在线调整法示意图 在不知道系统转移函数的情况下，以在线调整法，直接于 pid 控制器做调整，亦即 pid 控制器里的 i 值与 d 值设为零，只调 p 值让系统产

34、生震荡，这时的 p 值为临界震荡 增益v，之后震荡周期也可算出来，只不过在线调整实务上与系统仿真差别在于在实务 上处理比较麻烦，要在 pid 控制器输出信号端在串接电流表，即可观察所调出的 p 值是 否会震荡，虽然比较上一个 relay feedback 法是可免除拆装 relay 的麻烦，但是就经验而 青岛农业大学机电工程学院本科毕业设计 17 言在实务上线上调整法效果会较 relay feedback 差，在线调整法也可在计算机做出仿真调 出 pid 值，可是前提之下如果在计算机使用在线调整法还需把系统转移函数辨识出来， 但是实务上与在计算机仿真相同之处是 pid 值求法还是需要用到调整法

35、则 ziegler-nichols 经验法则去调整，与 relay feedback 的经验法则一样，调出 pid 值。 4.2.4在线调整法在计算机做仿真 step 1：以 matlab 里的 simulink 绘出反馈方块，如下图 4-9 所示 图 4-9 反馈方块图 pid 方块图内为 图 4-10 pid 方块图 step 2：将 td 调为 0，ti 无限大，让系统为 p 控制，如下图 4-11 所示。 图 4-11 step 3：调整 kp使系统震荡，震荡时的 kp即为临界增益 ku，震荡周期即为 tv。 （使 在线调整时，不用看 a 求 ku） ，如下图 4-12 所示。 青岛农

36、业大学机电工程学院本科毕业设计 18 图 4-12 系统震荡特性图 step 4：再利用 ziegler-nichols 调整法则，即可求出该系统之p、ti，td之值。 4.3针对有转移函数的 pid 调整方法 4.3.1系统辨识法 系统反馈方块图在上述无转移函数 pid 调整法则有在线调整法与 relay feedback 调整 法之外，也可利用系统辨识出的转移函数在计算机仿真求出 pid 值，至于系统辨识转移 函数技巧在第三章已叙述过，接下来是要把辨识出来的转移函数用在反馈控制图，之后 应用系统辨识的经验公式 ziegler-nichols 第二个调整法求出 pid 值，13如下表 4-1

37、4 所 示。 图 4-13 由系统辨识法辨识出转移函数 青岛农业大学机电工程学院本科毕业设计 19 controller p k i t d t p a 1 pi ()* a 9 . 0 a 6 . 0 3.3l pid ()* a 2 . 1 a 9 . 0 2l 2 l 表 4-14 ziegler-nichols 第二个调整法则建议 pid 调整值9 *为本专题将经验公式修正后之值 上表 4.3.1 中，l 为延迟时间可参考图 4.3.1(b)。 上表 4.3.1 中，a 的解法可有以下 2 种： 解一：如下图 4-15 中可先观察系统特性曲线图，辨识出 a 值。 解二：利用三角比例法推

38、导求得 l 0.632k k a t t” t 图 4-15 利用三角比例法求出 a 值 ak a t l aak a lt l )( 公式（4-2） k a lt l lt l ka 用 ziegler-nichols 第一个调整法则求得之 pid 控制器加入系统后，一般闭环系统阶跃 响应最大超越的范围约在 10%60%之间。 所以 pid 控制器加入系统后往往先根据 ziegler-nichols 第二个调整法则调整 pid 值，然后再微调 pid 值至合乎规格为止。 青岛农业大学机电工程学院本科毕业设计 20 4.3.2波德图法及根轨迹法 利用系统辨识出来的转移函数，使用 matlab

39、软件去做系统仿真。由于本设计中 pid 参数的整定主要是基于系统辨识及 ziegler-nichols 调整法则，所以在此不用波德图法 及根轨迹法。 青岛农业大学机电工程学院本科毕业设计 21 第五章 油冷却机系统的 pid 控制器设计 5.1 油冷却机系统 本论文设计以油冷却机温度控制系统为被控对象进行 pid 控制器的参数整定及其设 计，下面介绍一下油冷却机系统以及各个组成部分： 5.1.1油冷却机 青岛农业大学机电工程学院本科毕业设计 22 图 5-1 油冷却机实物图 图 5-2油冷却机系统循环图 油冷却机系统循环主要可分为冷媒循环系统以及油循环系统。冷媒循环系统即为一 般常见之制冷循环

40、，而油循环则是将油打出后经过负载加热，再与冷媒循环的蒸发器作 热交换，再流回油槽做冷却用。16 5.1.2感测与转换器 青岛农业大学机电工程学院本科毕业设计 23 图 5-3pt100 实物图 电阻式温度检测器 (rtd,resistance temperature detector)一种物质材料作成的电阻,它 会随温度的上升而改变电阻值,如果它随温度的上升而电阻值也跟着上升就称为正电阻系 数,如果它随温度的上升而电阻值反而下降就称为负电阻系数6。pt100 温度传感器是一种 以白金(pt)作成的电阻式温度检测器,属于正电阻系数,其电阻和温度变化的关系式如下： r=ro(1+t) 其中 =0.

41、00392,ro 为 100(在 0的电阻值),t 为摄氏温度 vo=2.55ma 100(1+0.00392t)=0.255+t/1000 。6 图 5-4 感测与转换电路 由于一般电源供应较多零件之后,电源是带噪声的,因此我们使用齐纳二极管作为稳压 零件,由于 7.2v 齐纳二极管的作用,使得 1k 电阻和 5k 可变电阻之电压和为 6.5v,靠 5k 可 变电阻的调整可决定晶体管的射(集极)极电流,而我们须将集极电流调为 2.55ma,使得量测 电压 v 如箭头所示为 0.255+t/1000。其后的非反向放大器,输入电阻几乎无限大,同时又放 大 10 倍,使得运算放大器输出为 2.55

42、+t/100。 6v 齐纳二极管的作用如 7.2v 齐纳二极管的作用,我们利用它调出 2.55v,因此电压追 青岛农业大学机电工程学院本科毕业设计 24 随器的输出电压 v1 亦为 2.55v。 其后差动放大器之输出为 vo=10(v2-v1)=10(2.55+t/100-2.55)=t/10,如果现在室温为 25, 则输出电压为 2.5v。 5.1.3 控制组件 am+12vvingndx5revcomfwdtc1comta1tb1iin x1y2m-x3cmefm+y1ta2x2x6x4comtc2 图 5-5 变频器实物图 变频器功用： 控制器输出的信号连结到最终控制组件（变频器）的输入

43、端，使最终控制组件依控 制器所输出信号改变频率。 在前面介绍了系统辨识和基于 matlab 进行 pid 参数整定的几种方法，现在就来设 计油冷却机系统中的 pid 控制器，也就是进行 pid 参数的整定。 5.2 油冷却机系统之系统辨识 下图为整个油冷却机系统示意图： 电源端 负载 端 青岛农业大学机电工程学院本科毕业设计 25 图 5-6 油冷却机系统示意图 从第二章中可以知道变频器是个最终控制组件，其以最大讯号(12v)驱动油冷却机， 使油冷却机以全速运转，温度记录器感测油槽记录温度，过一段时间后，因负载与油温 达到平衡，所以温度会到一稳定值，从纪录器撷取出温度数据至计算机，利用 mat

44、lab 绘出系统特性图，从而辨识出转移函数。下面是油冷却机系统的系统特性求解图： 图 5-7 系统特性求解图 step 1：系统的起始温度及稳定温度差值即为。 k=24.1-15=9.1 step 2：将系统特性图画为近似曲线作为辅助，并画出最大切线斜率，此切线斜率交 于温度轴，该温度与起始温度的差即为。 step 3：最大切线斜率与系统特性曲线交点的时间为延迟时间。即为 l=190 step 4：先找出 0.632k 的温度值，画水平线与特性曲线相交，得知时间轴之值，此 值与延迟时间相减，即为时间常数。 t=2700-190=2510 所以辨识出转移函数为: 12510 1 . 9 )( 1

45、90 s e sg s 5.3 油冷却机系统的 pid 参数整定 我们已经运用系统辨识法通过系统特性图获得了部分参数和系统的传递函数 青岛农业大学机电工程学院本科毕业设计 26 ，我们只要再使用表 4-14 的 ziegler-nicholg 第二经验法则，便可以求出 12510 1 . 9 )( 190 s e sg s pid 值。 p=1.5 a 9 . 0 6 . 0 9 . 0 i=2l=2190=380 d=95 2 l 2 190 pid 控制系统调整结果如图 5-8： 图 5-8 pid 控制系统调整结果 现在再用 relay feedback 法来进行 pid 参数整定。 由

46、辨识出的传递函数为，以 matlab 算出延迟 190 秒的一阶方程 12510 1 . 9 )( 190 s e sg s 式， num,den=pade(190,1) num = -1.0000 0.0105 den = 1.0000 0.0105 在 simulink 中绘出方块图 5-9 青岛农业大学机电工程学院本科毕业设计 27 图 5-9 simulink 中绘出的方块图 让 relay 做 on-off 动作，将系统扰动（on-off 动作，将以 做仿真） ，系统特性 曲线如图 5-10 图 5-10 系统扰动后的特性曲线 由上图 a=0.5,tu=760.将 a 值代入所以 5

47、 . 1 5 . 0 14 ku 0 . 95 . 16 . 06 . 0kukp 3807605 . 05 . 0tuti 95760125 . 0 125 . 0 tutd 与用系统辨识法所得结果大体相同。在线调整法一样，在此就不再详解。 在 pid 控制器的设计中并不是只要整定完参数就完成任务，还要进行 pid 值的修正 订以及加干扰仿真实验，下面就以系统辨识法求得的 pid 值进行修正，其修正大多是根 据经验。pid 修正调整结果大体为： p=1.5 青岛农业大学机电工程学院本科毕业设计 28 i=400 d=120 图 5-11 pid 值修正前后系统响应曲线对比图 系统干扰：因为系

48、统中负载随时在变化，系统中负载突然变动，也能马上追上设定 值。 系统加干扰仿真实验： 图 5-12 干扰仿真实验响应曲线 由上图可知道，就算系统临时发生变化，也能马上恢复到设定值。 青岛农业大学机电工程学院本科毕业设计 29 在计算机做完仿真操作之后，接着要实际配线，将所求的 pid 值输入到 pid 控制器， 让机器实际运转16。 图 5-13 pid 反馈控制在油冷却机的使用 pid 控制器的操作面板15 图 5-14 pid 控制器的操作面板 青岛农业大学机电工程学院本科毕业设计 30 1 2 3 4 5 21 22 23 24 25 6 7 8 9 10 26 27 28 29 30

49、al1 al2 not used ct input b b a rtdtc ac85 264v sv2 input trnsmit nc no c main output 图 5-15 tm-109(30-i)接线图 按 mode 键约 2 秒，可在 pid 控制器做设定17： 青岛农业大学机电工程学院本科毕业设计 31 图 5-16 pid 控制器参数设定流程示意图 pid 控制器： 1.输入讯号端：电压或电流信号(420ma) 2.输出讯号端：电压或电流信号(420ma) 3.操作电源 青岛农业大学机电工程学院本科毕业设计 32 结 论 本论文是设计一个温度控制系统的 pid 控制器。pi

50、d 调节器从问世至今已历经了半个 多世纪，在这几十年中，人们为它的发展和推广做出了巨大的努力，使之成为工业过程 控制中主要的和可靠的技术工具。即使在微处理技术迅速发展的今天，过程控制中大部 分控制规律都未能离开 pid，这充分说明 pid 控制仍具有很强的生命力。pid 控制中一个 至关重要的问题，就是控制器三参数（比例系数、积分时间、微分时间）的整定。整定 的好坏不但会影响到控制质量，而且还会影响到控制器的鲁棒性1。所以本文重要的是来 介绍 pid 参数整定的方法。 在第一章绪论中说明了本课题的意义，matlab 软件的应用以及在这个方面的发展 趋势。第二章简单介绍了控制系统和 pid 调节

51、之后，在第三章中介绍了一种重要的求系 统传递函数的方法系统辨识法，这是借助 matlab 进行 pid 参数整定的前提。在第 四章中研究了基于 matlab 的 pid 参数整定的几种方法，并各举一例予以说明，主要有 relay feedback 法，在线调整法以及系统辨识法，波得图法及根轨迹法不做研究。以上几 章的知识是控制器设计的基础，在接下来的第五章里面便结合油冷却机系统应用以上几 章的知识进行 pid 控制器的参数整定及设计。 青岛农业大学机电工程学院本科毕业设计 33 致 谢 经过几个月的不断努力，毕业设计终于如期完成。从拿到设计题目到最后成设计并 定稿，其间经历了翻阅相关资料、熟悉

52、基础知识、学习 matlab 软件的使用，到开始写 论文以及最后的修改和装订成册这几个阶段。每个阶段工作的完成都使我在各个方面受 益匪浅。在这次毕业设计中，我的任务是完成基于 matlab 的温控系统 pid 控制器设计。 为了很好地完成设计任务，我经常去图书馆或者上网搜集各种资料文献，向指导老师和 各位同学请教，并且翻阅以前的课本、笔记，熟悉之前学过的相关知识。这些不仅仅巩 固了我以前所学的专业知识，而且使我接触了许多以前没接触过的新知识，大大地扩宽 了我的知识面。尤其是对于 pid 控制器的设计和应用，使我有了更加深入的了解，也使 知道了在现代的控制系统设计和建立中借助好的软件包的重要性及

53、未来的发展趋势。在 这次设计过程中，我明显感觉到自己在许多方面存在不足，譬如，对 word 的熟练使用， 对 matlab 软件的应用，对 pid 控制器的认识等等。我借此机会不断学习，努力提高多 方面的能力，弥补自己的不足。总的说来，通过这次毕业设计的完成，我在各方面都有 了很大的进步。特别是将大学所学的专业理论知识运用于实际设计中，与我们的生产、 生活联系在一起，让我对自己的专业有了更浓厚的兴趣，对自己的前途有了更充足的信 心和更美好的憧憬。 在整个设计过程中，我得到了陈俊东老师的悉心指导和帮助。在我遇到困难时，他 们总是及时地帮助我理清思路解决困惑，最终跨过了一个又一个障碍，顺利地完成了

54、毕 业论文的设计工作。在此，我向陈俊东老师表示衷心的感谢，向论文答辩组的各位老师 表示衷心的感谢。 由于本人的知识和能力有限，论文中难免存在许多不足的地方，敬请各位老师给予 批评指正. 青岛农业大学机电工程学院本科毕业设计 34 参考文献 1 夏红，赏星耀，宋建成. pid 参数自整定方法综述. 浙江科技学院学报, 2003,15(4)：236-240. 2 astrom k j. toward intelligent control. ieee control systems magazine,1989(april):6064. 3 astrom k j,hang c c,persson p

55、, ho w k.towards intelligent pid control.automatica,1992,28(1): 19 4 王 伟,张晶涛.pid 参数先进整定方法综述.自动化学报，2000,26(3):347-355. 5 薛定宇.反馈控制系统设计与分析-matlab 语言应用.清华大学出版社, 2000,250-251. 6 江明崇.感测与转换器.台北市：全华图书有限公司.1995（2）44-45 7 陈聪明.冷冻空调之自动控制. 台北市:全华图书有限公司.1996(3):89-91. 8应自炉. matlab 软件在控制系统仿真与分析中的应用.现代电子技术.1995(4):

56、12- 15. 9 aidan odwyer, handbook of pi and pid controller tuning rules, new jersey, world scientific, c2003:125-129. 10 李言俊,张科.系统辨识理论及应用.国防工业出版社,2002,16-19. 11 bristol e h.a simple adaptive system for industrial control. instrumentation technology, 1967(june), 156-159. 12 佚名.pid 调节概念及基本理论.2006,中国自动化

57、网. 13 原菊梅.参数整定的研究.北京工商大学学报(自然科学版) 2004,22(4),40-42. 14 刘金琨.先进 pid 控制及 matlab 仿真.电子工业出版社,2002 年 10 月,214-218. 15 沈金钟.pid 控制器 :理论.调整与实现. 台中市:沧海书局出版社 ,2001,68-74. 16 陶永华.新型 pid 控制及其应用.北京:电气自动化新技术丛书,2002 年 2 版,131-140. 17 aniruddha datta, ming-tzu ho, and shankar p. bhattacharyya, structure and synthesi

58、s of pid controllers, new york, 2000, 75-81. 09/20 11:46 102 机体齿飞面孔双卧多轴组合机床及 cad 设计 09/08 20:02 3kn 微型装载机设计 青岛农业大学机电工程学院本科毕业设计 35 09/20 15:09 45t 旋挖钻机变幅机构液压缸设计 08/30 15:32 5 吨卷扬机设计 10/30 17:12 c620 轴拨杆的工艺规程及钻 2-16 孔的钻床夹具设计 09/21 13:39 ca6140 车床拨叉零件的机械加工工艺规程及夹具设计 831003 08/30 15:37 cpu 风扇后盖的注塑模具设计 09

59、/20 16:19 gdc956160 工业对辊成型机设计 08/30 15:45 ls 型螺旋输送机的设计 10/07 23:43 ls 型螺旋输送机设计 09/20 16:23 p-90b 型耙斗式装载机设计 09/08 20:17 pe10 自行车无级变速器设计 10/07 09:23 话机机座下壳模具的设计与制造 09/08 20:20 t108 吨自卸车拐轴的断裂原因分析及优化设计 09/21 13:39 x-y 型数控铣床工作台的设计 09/08 20:25 yd5141syz 后压缩式垃圾车的上装箱体设计 10/07 09:20 zh1115w 柴油机气缸体三面粗镗组合机床总体及左主轴箱设计 09/21 15:34 zxt-06 型多臂机凸轮轴加工工艺及工装设计 10/30 16:04 三孔连杆零件的工艺规程及钻 35h6 孔的夹具设计 08/30 17:57 三层货运电梯曳引机及传动系统设计 10/29 14:08 上盖的工