2022年2022年计算机控制原理温度控制系统课程设计 .pdf

《2022年2022年计算机控制原理温度控制系统课程设计 .pdf》由会员分享，可在线阅读，更多相关《2022年2022年计算机控制原理温度控制系统课程设计 .pdf（23页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、计算机控制课程设计报告名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 23 页 - - - - - - - - - 目录1设计任务和控制要求1.1 设计任务1.2 控制要求2 设计方案2.1 温度控制系统组成2.2 设计步骤3 被控对象的参数设定3.1 数据记录32 T、 实验数据4算法原理5结果分析5.1 PID 调节结果5.2大林算法控制与处理6心得体会7 参考文献附录附录 1：原理图附录 2：程序名师资料总结 - - -精品资料欢迎下载 - - - - - - -

2、- - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 23 页 - - - - - - - - - 参考文献附录 1：电路原理图附录2：源程序温度控制系统设计1设计任务和控制要求1.1 设计任务计算机控制技术的课程设计是一个综合运用知识的过程，它不仅需要微型计算机控制理论、程序设计和接口技术等方面的基础知识，而且还需要具备一定的生产工艺知识。设计包括确定控制任务、系统总体方案设计、系统设计、控制软件的设计等，以便使学生掌握微型计算机控制系统设计的总体思路和方法，做到理论与实践的结合。名师资料总结 - - -精品资料欢迎下载 - - - - - -

3、 - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 23 页 - - - - - - - - - 1.2 控制要求1.2.1 被控对象：本设计的被控对象为电热丝。其被控对象温控数学模型为1T sess-KG）（1.2.2 温度范围：室温 100C ，控制精度正负 1C。1.2.3 输出形式：第一种 : 0-10mA 电流信号第二种：脉宽调制（ PWM）信号输出1.2.4 设计要求（1）数字显示温度，具有报警、参数设定、温升工艺曲线设置与修改等功能。（2）绘制温升曲线及其它处理；1.2.5 选用合适的控制算法（1）PID 控制PID控制器由

4、比例控制单元（ P） 、积分控制单元（ I ）和微分控制单元（ D ）组成。其输入）（te与输出）（tu的关系为)/ )(*)(/ 1)()(dttdeTdtteTteKptuDI式中积分的上下限分别是 0 到 t 因此，它的传递函数为：) s*11 ()()()(DITsTKpsEsUsG名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 23 页 - - - - - - - - - 其中 Kp 为比例系数； TI为积分时间常数； TD 为微分时间常数（2）Dalin算法大

5、林算法是运用于自动控制领域中的一种算法，是一种根据被控对象先设计好理想闭环系统的控制器，再反过来综合调节的方法。设计的目标是使闭环传递函数)(s相当于一个纯滞后环节和一个惯性环节的串联，其中纯滞后环节的滞后时间与被控对象的纯滞后时间完全相同。这样就能保证使系统不产生超调，同时保证其稳定性。2 设计方案2.1 温度控制系统组成温度控制系统框图如图 1所示。该系统由以下部分组成：PID脉宽调制加温室温度变送器fUgUoU图1 温度控制系统框图工作原理在图1 所示系统中，加热室（加热丝）为该系统的被控对象，温度变送器用来检测加热室的温度，并转换为电压信号ou，检测到的电压信号经输出电压反馈单元作为系

6、统反馈值fu ，脉宽调制电路为加热室的加热执行单元。在参数给定值gu 的情况下，给定值gu 与反馈值fu 比较得到偏差，经过PID 调节器运算产生相名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 23 页 - - - - - - - - - 应的控制量， PID调节器的输出作为脉宽调制的输入信号，经脉宽调制电路产生占空比可调的脉冲信号，作为对加温室里电热丝的加热信号，使加温室（加热丝）里的温度稳定在给定值gu 。给定值gu 可由已知电源提供。温度变送器可采用热敏电阻等温度

7、传感器。2.2 设计步骤2.2.1被控对象（电热丝）模型通过实验方法确定被控对象数学模型（即电热丝的传递函数） ，该被控对象数学模型为1T sess-KG）（2.2.2 确定T、 T 为被控对象时间常数； 为滞后时间常数。2.2.3 确定被控对象数学模型在模拟控制系统中，参数整定方法较多，常用的实验整定方法有：临界比例度法、阶跃响应曲线法、试凑法等。数字控制时也可采用类似方法，如扩充的临界比例度法、扩充的阶跃响应曲线法与试凑法等等。本次设计采用扩充响应曲线法来确定被控对象模型。具体步骤如下：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - -

8、 - - 名师精心整理 - - - - - - - 第 6 页，共 23 页 - - - - - - - - - （1）首先不将数字控制器接入控制系统，让系统处于开环工作状态下，给定一阶跃信号，每10s 记录一次数据，直至系统输出稳定在给定值附近。（2）根据记录数据画出被调量在阶跃输入下的整个变化过程曲线。（3）在曲线最大斜率处做切线，求得滞后时间 ，被控对象时间常数1T以及它们的比值1T/ 。3 被控对象的参数设定3.1 数据记录如下 ：表1 第一次给定电压： 5V 每10s 记一数据1.64 3.54 4.21 4.69 4.98 5.19 5.34 5.43 5.49 1.67 3.58

9、 4.23 4.71 4.99 5.2 5.34 5.43 5.49 1.86 3.62 4.27 4.73 5 5.21 5.34 5.44 5.49 2.06 3.66 4.3 4.75 5.01 5.22 5.35 5.44 5.5 2.28 3.69 4.32 4.77 5.03 5.22 5.35 5.44 5.5 2.46 3.73 4.35 4.78 5.04 5.23 5.36 5.44 5.5 2.58 3.76 4.37 4.8 5.05 5.24 5.36 5.45 5.5 2.72 3.8 4.4 4.81 5.06 5.25 5.37 5.45 5.51 2.81 3

10、.84 4.42 4.83 5.07 5.26 5.37 5.46 5.51 2.92 3.87 4.45 4.84 5.08 5.26 5.38 5.46 5.51 3.01 3.91 4.48 4.85 5.09 5.27 5.38 5.46 5.51 3.07 3.94 4.51 4.87 5.1 5.28 5.39 5.46 5.52 3.13 3.97 4.53 4.89 5.11 5.29 5.39 5.47 5.52 3.19 4.01 4.56 4.9 5.12 5.29 5.4 5.47 5.52 3.25 4.04 4.58 4.92 5.13 5.3 5.4 5.48 3

11、.3 4.07 4.6 4.93 5.14 5.3 5.41 5.48 3.36 4.1 4.61 4.94 5.15 5.31 5.41 5.48 3.41 4.13 4.63 4.95 5.16 5.32 5.41 5.48 3.47 4.15 4.65 4.96 5.17 5.33 5.42 5.48 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页，共 23 页 - - - - - - - - - 3.51 4.18 4.67 4.97 5.18 5.33 5.42 5

12、.49 表2 第二次给定电压 4V 每10s 记一数据终值 4.75 2.11 3.36 3.85 4.17 4.37 4.49 4.57 4.66 4.72 2.13 3.39 3.87 4.18 4.38 4.49 4.58 4.66 4.72 2.26 3.42 3.89 4.19 4.38 4.49 4.58 4.66 4.72 2.39 3.45 3.91 4.21 4.4 4.51 4.58 4.66 4.73 2.53 3.48 3.92 4.22 4.4 4.51 4.59 4.66 4.73 2.64 3.5 3.94 4.24 4.41 4.51 4.59 4.69 4.7

13、3 2.73 3.53 3.96 4.24 4.41 4.52 4.59 4.69 4.73 2.81 3.56 3.97 4.25 4.43 4.52 4.62 4.69 4.73 2.87 3.59 3.99 4.25 4.43 4.52 4.62 4.69 4.73 2.93 3.62 4.01 4.27 4.44 4.54 4.62 4.69 4.73 2.98 3.64 4.02 4.27 4.44 4.54 4.62 4.69 4.73 3.03 3.67 4.04 4.3 4.45 4.54 4.62 4.71 4.73 3.07 3.69 4.05 4.3 4.45 4.55

14、4.62 4.71 4.73 3.11 3.71 4.07 4.31 4.46 4.55 4.64 4.71 3.15 3.74 4.08 4.31 4.46 4.55 4.64 4.71 3.2 3.76 4.1 4.34 4.47 4.56 4.64 4.71 3.24 3.78 4.11 4.34 4.47 4.56 4.64 4.71 3.27 3.8 4.13 4.35 4.47 4.56 4.64 4.72 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页，共 23

15、 页 - - - - - - - - - 3.31 3.81 4.14 4.35 4.47 4.57 4.64 4.72 3.33 3.83 4.16 4.37 4.47 4.57 4.66 4.72 表3 第三次给定电压 6V 每10s 记一数据终值 6.45 1.81 4.03 4.91 5.41 5.72 5.91 6.05 6.17 6.25 6.3 6.34 6.39 1.88 4.09 4.94 5.43 5.73 5.91 6.06 6.17 6.25 6.3 6.34 6.39 2.04 4.15 4.97 5.45 5.74 5.93 6.06 6.18 6.25 6.3 6

16、.35 6.41 2.23 4.2 5 5.46 5.75 5.93 6.06 6.18 6.25 6.3 6.35 6.41 2.42 4.26 5.03 5.48 5.76 5.94 6.08 6.18 6.25 6.3 6.35 6.41 2.59 4.32 5.06 5.5 5.77 5.94 6.08 6.19 6.25 6.32 6.35 6.41 2.73 4.37 5.09 5.51 5.78 5.96 6.08 6.19 6.26 6.32 6.35 6.41 2.88 4.42 5.11 5.53 5.79 5.96 6.1 6.19 6.26 6.32 6.35 6.41

17、 3.01 4.46 5.14 5.54 5.81 5.97 6.1 6.22 6.26 6.32 6.35 6.41 3.12 4.52 5.17 5.56 5.82 5.97 6.1 6.22 6.26 6.32 6.38 3.27 4.56 5.19 5.57 5.83 5.98 6.12 6.22 6.26 6.32 6.38 3.35 4.6 5.22 5.59 5.84 5.98 6.12 6.22 6.26 6.34 6.38 3.43 4.64 5.24 5.61 5.84 6 6.12 6.22 6.28 6.34 6.38 3.51 4.67 5.27 5.62 5.85

18、6 6.14 6.22 6.28 6.34 6.38 3.6 4.71 5.29 5.63 5.86 6 6.14 6.24 6.28 6.34 6.38 3.68 4.74 5.31 5.65 5.86 6.03 6.14 6.24 6.28 6.34 6.38 3.75 4.77 5.33 5.67 5.88 6.03 6.16 6.24 6.28 6.34 6.39 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 9 页，共 23 页 - - - - - - - - - 在曲线

19、最大斜率处做切线， 求得滞后时间 ，被控对象时间常数1T以及它们的比值1T/ 。数据如下：32 T、 实验数据次数1TK1T/第一次205 9.5 7.96 21.58 第二次220 10 7.87 21 第三次240 10.5 7.7 21.5 平均值221.6 10 7.8 21.5 1PID 参数整定3.82 4.8 5.35 5.68 5.88 6.03 6.16 6.24 6.28 6.34 6.39 3.9 4.84 5.37 5.69 5.89 6.05 6.16 6.24 6.3 6.34 6.39 3.98 4.88 5.39 5.7 5.89 6.05 6.17 6.24

20、6.3 6.34 6.39 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 10 页，共 23 页 - - - - - - - - - 根据实验得到的表 4， 初步确定控制器的 Kp， Ki， Kd 及采样周期 T 。控制度控制律采样周期Kp ITDT1.05 PI 0.1 0.841T3.4 PID 0.05 0.151T2.0 0.45 1.2 PI 0.2 0.781T3.6 PID 0.16 1.01T1.9 0.55 1.5 PI 0.5 0.681T3.9 PID 0.3

21、4 0.851T1.620.82 表5 PID 参数整定表在运行中，对上述参数作适当调整，以获得满意的性能。扩充响应曲线法通过测取阶跃响应曲线的t ，T 参数获得一个初步的PID 控制参数，然后在此基础上通过部分参数的调节（试凑）获得满意的控制性能。4算法原理计算机控制系统的方框图如图所示。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 11 页，共 23 页 - - - - - - - - - 图计算机控制系统框图1)常规数字 PID 控制算法常规的 PID 控制律为)()(dtt

22、deTdtteKtudP采用一阶差分法离散化后， 可以得到常规数字PID 控制位置式算法)1()()()()(1kekeTTieTTkeKkudkiiP简记为kikekeDieIkPeku1)1()()()()(这里 P、I 、D参数分别为TTKDTTKIKPPPPdi采用增量式形式有：2)-e(k1)-2e(k-De(k)Ie(k)1)-e(k-Pe(k)1)-u(ku(k)2)数字PID调节a.根据表 3，初步确定控制器的 Kp，Ki，Kd 及采样周期 T ，编写程序；将数字控制器接入控制系统，让系统处于闭环工作状态下，给定一阶跃信号，每10s 记录一次数据，直至系统输出稳定在给定值附近。

23、b.根据记录数据，画出系统在阶跃输入作用下的整个变化过名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 12 页，共 23 页 - - - - - - - - - 程曲线。2. 大林算法根据被控对象的S传递函数式（ 1） ，大林算法选定0( )1sesT snT0T 按控制要求选择。作为闭环控制的综合目标，与( ) s相对应的00/1/11(1)( )( )1T TTsnT TeezzZssez而包含零阶保持器被控对象的S传递函数为0111( )( )1TsTsseeKeG sG ss

24、sTs离散化后得到11/1/1(1)( )( )1T TnT TKezG zZ G sez于是可以得到大林算法控制器的Z传递函数01001/1/11(1)(1)( )(1)1(1)T TT TT TT TT TneezD zKeezez由此得到大林算法0000111/( )(1)(1) (1)(1)(1)( )(1)(1)(1)T TT TT TT TT TT TT Tu keu keu kneeee ke kKeKe5结果分析数据记录及结果分析5.1pid 调节结果名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理

25、- - - - - - - 第 13 页，共 23 页 - - - - - - - - - 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 14 页，共 23 页 - - - - - - - - - 根据记录数据及温升曲线得出PID参数对性能的影响：增大比例系数 Kp一般将加快系统的响应，在有静差的情况下有利于减小静差。 但过大会使系统有较大的超调， 并产生振荡，使系统稳定性变坏。增大积分时间 Ti 有利于减小超调，减小振荡，使系统更加稳定，但系统静差的消除将随之减慢。增大微分时间

26、Td有利于加快系统响应，使超调量减小，稳定性增加，但系统对扰动的抑制能力减弱， 对扰动有较敏感的响应。5.2大林算法控制与处理温度时间关系如下图示：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 15 页，共 23 页 - - - - - - - - - 6 心得体会计算机控制技术是一门以电子技术、自动控制技术、 计算机应用技术为基础，以计算机控制技术为核心，综合可编程控制技术、单片机技术、 计算机网络技术， 从而实现生产技术的精密化、生产设备的信息化、生产过程的自动化及机电控制系统

27、的最佳化的专门学科。通过本次课程设计，让我受益良多。原来我学习计算机控制的时候，基本上不清楚学习的理论知识如何应用在实践中，所以觉得计控学习起来比较难，经过这次课程设计，让我清楚了计算机控制的一些理论在实践中的应用。课程设计一开始，老师让我们先参观实验室，这是很重要的一个环节，参观让我们对实验的仪器有了个感性的认识，让我们明白我们需要的实验仪器是怎么样的，是我们做到心中有数。接下来的一段日子就是课程设计了，时间三个星期，课程设计一开始，大家都非常积极主动的思考怎么完成任务，我们是三个人一组，一起讨论该怎么做，如何测数据，商量好了，再进行下一步。在这个过程中并不总是一帆风顺名师资料总结 - -

28、-精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 16 页，共 23 页 - - - - - - - - - 的。实验一开始就遇到了问题：我的实验箱最开始不能用，经过了几次调换之后发现坏箱子还真是不少啊，最后一个管子烂了，但是老师说不影响实验的结果，于是我们就用到那个不好的做。在了解实验的理论的过程中也遇到了很多困惑，首先，vb 的脚本语言我们都不是很熟悉，我们就在网上下载了关于机器类似的编程的语言的讲解，在结合老师发的实验指导书，我发现对这语言也算是有了一定的了解。接下来就是具体实验了，说实话，实验内容并

29、不复杂，但是问什么要用三个星期的时间去完成呢，我觉得最主要的原因是因为这个实验确实需要花时间，因为每一组数据都需要差不多一个小时的时间去记录，而且要测量很多组的数据。实验结束后，我觉得最大的感受就是，做这个实验，最最需要的是耐心。经过了这次课程设计， 我发现学到了很多的东西， 首先我觉得能把自己所学到的东西应用到实践中去，这是一个很大的收获。其次，仅仅是所学到的东西根本就不能满足工程或者说一个具体的项目的要求，要完成一个具体的项目，需要额外的学习很多新的东西，因此不能满足于自己的所学， 应该养成终身学习的态度。附录 1：原理图名师资料总结 - - -精品资料欢迎下载 - - - - - - -

30、 - - - - - - - - - - - 名师精心整理 - - - - - - - 第 17 页，共 23 页 - - - - - - - - - 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 18 页，共 23 页 - - - - - - - - - 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 19 页，共 23 页 - - - - - - - - - 附录 2

31、：程序、1. PID 整定option explicit public KP,KI,KD public T,Ti,Td public e0,e1,e2,u public x2,x1,x0 public sume 初始化通道参数Sub Init() Daq.SetChannelXNumber 1 X 通道使用 A/D 通道 1 Daq.SetChannelYNumber 2 Y 通道使用 A/D 通道 2 Daq.SetSignal1State true D/A 通道 1 开通Daq.SetSignal2State true D/A 通道 2 关闭Daq.SetExpLastTime 0 持续时

32、间， 0 表示持续Daq.SetExpSamplePeriod 1000 采样周期， 10ms T=0.5 Ti=40 Td=50 Kp=5 KI=Kp*T/Ti KD=Kp*Td/T e0=0 e1=0 e2=0 sume=0 u=0 End Sub 测试信号发生函数Sub Signal() Dim eWave(1000) 测试信号定义， 1000 为每个周期点数dim i dim A,Ts Ts=1000 A=4 Signal1.SetPeriod Ts 测试信号 1 的周期， 单位 ms， 周期 1000ms Signal1.SetAmp A 测试信号 1 的幅值，单位 V，幅值为 4

33、For i = 0 To ubound(ewave)-1 eWave(i) =A Next Signal1.setwave eWave End Sub 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 20 页，共 23 页 - - - - - - - - - 运行 文件Sub Run() e0=e1 e1=e2 e2=Signal1.getcursignal()-daq.getchannelXvalue(0) sume=sume+e2 u=KP*e2+KI*sume+KD*(e2-e

34、1) Daq.SetSignal1Value u Daq.SetSignal2Value e2 End Sub 2. 大林算法option explicit public e0,e1,e2,u0,u1,u2 public T,T1,T0,K,m1,m2,k1,k2 初始化通道参数Sub Init() Daq.SetChannelXNumber 1 X 通道使用 A/D 通道 1 Daq.SetChannelYNumber 2 Y 通道使用 A/D 通道 2 Daq.SetSignal1State true D/A 通道 1 开通Daq.SetSignal2State false D/A 通道

35、2 关闭Daq.SetExpLastTime 0 持续时间， 0 表示一值持续Daq.SetExpSamplePeriod 4000 采样周期， 4000ms T=4000/1000 T1=320 被控对象惯性环节时间常数K=0.8 放大倍数T0=150 理想闭环系统一节惯性环节时间常数e0=0 e1=0 e2=0 u0=0 u1=0 u2=0 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 21 页，共 23 页 - - - - - - - - - m1=exp(-T/T0) m2

36、=1-m1 k1=m2/(K*(1-exp(-T/T1) k2=-exp(-T/T1)*k1 End Sub 测试信号发生函数Sub Signal() Dim eWave(1000) 测试信号定义， 1000为每个周期点数dim i dim A,Ts Ts=1000 A=4 Signal1.SetPeriod Ts 测试信号 1 的周期，单位 ms， 周期1000ms Signal1.SetAmp A 测试信号 1 的幅值，单位 V，幅值为 4 For i = 0 To ubound(ewave)-1 eWave(i) =A Next Signal1.setwave eWave End Sub

37、 运行 文件Sub Run() u0=u1 u1=u2 e0=e1 e1=e2 e2=Signal1.getcursignal()-daq.getchannelXvalue(0) u2=m1*u1+m2*u0+k1*e2+k2*e1 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 22 页，共 23 页 - - - - - - - - - Daq.SetSignal1Value u1 延迟一个时间，故用u1 来输出End Sub 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 23 页，共 23 页 - - - - - - - - -