对具有纯滞后的一阶惯性环节的设计

1、 中央民族大学信息工程学院 计算机控制技术综合设计实验报告学生姓名： 学 号： 指导教师： 实验地点： 实验名称： 加热炉系统温度控制器设计目录一、设计题目及要求3二、设计方案与结构图31、计算机控制系统结构图32、硬件结构图4三、 电路硬件设计41、 电桥电路42、放大环节53、滤波电路64、A/D转换器65、D/A 转换电路7四、参数计算及仿真81、时数字调节器D（z）的实现8a、无控制作用下系统伯德图8b、最少拍下调节器函数9C、最少拍下系统伯德图11d、单位阶跃响应下系统输出12e、施加阶跃干扰信号13f、施加随机信号影响142、时数字调节器D（z）的实现16a、无控制作用下系统伯德图

2、16b、达林算法下调节器函数17C、达林算法下系统伯德图（未加增益）19d、达林算法下系统伯德图（加增益）20e、单位阶跃响应下系统输出21f、施加阶跃干扰信号22g、施加随机信号影响24五、心得与体会25一、设计题目及要求1、 针对一个具有纯滞后的一阶惯性环节的温度控制系统和给定的系统性能指标： 工程要求相角裕度为3060，幅值裕度6dB 要求测量范围-50200，测量精度0.5，分辨率0.22、 书面设计一个计算机控制系统的硬件布线连接图，并转化为系统结构图；3、 选择一种控制算法并借助软件工程知识编写程序流程图；4、 用MATLAB和SIMULINK进行仿真分析和验证; 对象确定：K=1

3、0*log(C*C-sqrt(C),rand(state,C),T=rang(1), 考虑=0或T/2两种情况。 C为学号的后3位数，如C=325，K=115.7，T=0.9824，=0或0.49125、 进行可靠性和抗干扰性的分析。二、设计方案与结构图1、计算机控制系统结构图被控过程保持器数字控制器w(t) + e(t) e(kT) u(kT) u(t) y(t)- T其控制过程可描述如下：1) 只有在采样开关闭合（即采样）的kT时刻，才对系统误差e(t)的瞬时值进行检测，也就是将整量化了的数字量e(kT)输入给计算机（数字控制器）。这一过程称为实时采集。2) 计算机对所采集的数据e(kT)

4、进行处理，即依给定的控制规律（数字控制器）确定该kT采样时刻的数字控制量u(kT)。这一过程称为实时决策。3) 将kT采样时刻决策给出的数字控制量u(kT)转换为kT时刻生效的模拟控制量u(t)控制被控对象。这一过程称为实时控制。2、硬件结构图被控对象温度控制DA转换单片机AD转换温度输出信号滤波信号放大温度检测设计方案：温度检测环节通过温度传感器实时测量温度并通过差分放大电路放大其信号，之后通过滤波环节除去信号中参杂的高频干扰，利用偏差信号进行控制。之后对信号进行采样，通过AD转换器将模拟信号转换成数字信号送入到单片机中去处理，单片机通过程序输出控制信号，控制信号通过DA转换器转换成模拟信号

5、之后对执行元件进行控制，使被控参数符合系统的要求。而该系统采用闭环控制，具有自动调节的能力，将温度控制在一定的范围内，由于单片机强大的算术运算和逻辑运算功能，使得高精度高性能的调节方案可以通过软件实现，并且具有良好扩展性(例如功能扩展)，故本系统理论上可以实现良好的控制效果并尽可能地智能化。三、 电路硬件设计1、 电桥电路电桥电路是用来把传感器的电阻、电容、电感变化转换为电压或电流。分直流电桥和交流电桥，交流电桥主要用于测量电容式传感器和电感式传感器的电容和电感的变化，直流电桥主要用于电阻式传感器。我们这里是用热敏电阻来测量温度的变化，所以先用直流电桥。而电桥可以分为单臂电桥，半桥差动和全桥差

6、动。全桥差动和板桥差动虽然输出信号对电源的影响减小，但增加传感器的个数，为了节约成本，可选择单臂电桥。电桥电路如下：图中R5,R6,R7,R8和R2构成电桥电路，其中R2起的作用是，在传感器变化为零时，调节R2使电桥平衡，输出为零。假设R5=R6=R7=R8=R，电桥输出。2、放大环节由于电桥输出的信号十分微弱，所以放大电路采用义用放大电路对信号进行放大，电路图如下：其中义用放大电路的放大倍数为：如上图所示的放大位数为79.5倍。由于放大信号很微弱，而且电桥输出电阻较高，所以义用放大电路的运算放大器要求输入电阻很高，这里可采用场效应管组成的输入电路T1084。3、滤波电路为减少或消除外界干扰，

7、而且温度变化的频率很慢，所以我们使用幅频特性较好的二阶有源低通滤波电路，电路图如下：图中C1为0.2uF，C2为100pF,R13=R10=10K，低通滤波器的上限频率略等于314HZ，而且设定放大倍数为1。4、A/D转换器根据题目要求，温度测量范围为，分辨率为，则测量是量程为，再由，可得。可取，即采用12位的ADC和12位的DAC。当采用12位ADC进行模/数转换时,其分辨率为，远远满足设计要求。可采用ADC574A芯片。电路图如下：8051的P0口作为AD574A的地址线，P0口和P2.0、P2.1、P2.2、P2.3口作为数据线，用于接收获取AD574A的转换结果。P0口经地址锁存器74

8、373锁存，并经三-八译码器74138译码后的信号作为AD574A的片选信号输入。5、D/A 转换电路 由上分析，本设计需要采用12位的A/D转换器，与之相对应，可采用12位的D/A转换器DAC1208芯片进行模/数转换。8051的P0口作为DAC1208的地址线，P0口和P2.4、P2.5、P2.6、P2.7口作为数据线，用于传送经达林算法后的运算结果。P0口经地址锁存器74373锁存，并经三-八译码器74138译码后的信号作为DAC1208的片选信号输入。8051与DAC1208的连接图如下所示：四、参数计算及仿真由题目要求，其中C为学号后三位，即466，则有下式： 而；从而考虑 和两种情

9、况。1、时数字调节器D（z）的实现a、无控制作用下系统伯德图此时使用MATLAB画出其频率特性曲线图，程序如下： num=122.88; den=0.7479,1; w=logspace(-1,4); bode(num,den,w); grid on; title(系统122.88/(0.7479s+1)的Bode图)BODE图如下：无控制作用下系统伯德图程序：由上图可知，系统幅值裕度为无穷(Inf)，相角裕度为90.5deg。相角裕度和幅值裕度均不符合要求。需要对此进行改进。b、最少拍下调节器函数因为被控对象为一阶无纯滞后环节，因此可考虑使用数字 PID 算法或者最少拍设计方法等进行控制器设

10、计。现采用最少拍设计方法设计该对象的控制器 D(Z) ，采样时间Ts选择Ts=0.374s 。所为最少拍系统或最少调整时间系统，是指对下图所示的随动系统，在给定某种典型输入条件下，设计一个控制规律D（z），使闭环系统具有最快的响应速度，也就是能在最少个采样周期时间内准确跟踪输入给定值，即达到无静差的稳定状态。r(t) + e(k) u(k) u(t) y(t) -数字控制随动系统广义工业对象脉冲传递函数为：由上式知，所以解得所以则数学调节器为：因此可求得在单位阶跃输入下，系统期望输出序列：显然，当（即第一拍）后，。也就是经过一拍，系统响应输出序列已准确跟踪输入函数且在采样点上的稳态误差为零。系

11、统的开环传递函数为：C、最少拍下系统伯德图使用MATLAB进行Bode图仿真，其程序如下： Ts=0.374; num=1; den=1 -1; Zk=tf(num,den,Ts); margin(Zk)可得出伯德图如下：最少拍下系统伯德图实验程序：实验验证Bode图由上图可知，系统幅值裕度为无穷(Inf)，相角裕度为60deg。符合题目要求。d、单位阶跃响应下系统输出使用MATLAB进行SIMULINK仿真，仿真图如下：其中输入阶跃响应如下： 输入波形图系统输出如下：最后拍系统输出波形图由上图可知，系统在第一拍后。也就是经过一拍，系统响应输出序列已准确跟踪输入函数且在采样点上的稳态误差为零。

12、跟计算结果一样。则所设计系统成立。满足系统要求。e、施加阶跃干扰信号仿真图如下：施加干扰信号：施加干扰波形图图中输出波形如下：施加干扰后最少拍系统输出波形由上图可看出，系统的抗干扰性能较好，能在系统受到外界干扰后在较短时间内回到稳定值。f、施加随机信号影响为了研究系统跟随性能优劣，在输入处施加随机干扰信号，并观察其输出信号变化。仿真图如下： 施加随机信号仿真图随机信号波形如下： 随机信号输出波形系统输出波形图如下： 施加随机信号最少拍系统输出波形由上图可知，最少拍系统跟随性能好，能准确地跟上信号变化。由以上分析知当时，以最小拍来设计能够满足题目要求，系统能在较少时间内达到无静差，并且该系统具有

13、较好的抗干扰性和稳定性。2、时数字调节器D（z）的实现a、无控制作用下系统伯德图此时，使用MATLAB画出其频率特性曲线图，程序如下： k=122.88; den=0.749 1; G=tf(k,den,inputdelay,0.374); w=logspace(0,2,5); margin(G); title(系统122.88/(0.7479s+1)的Bode图)BODE图如下：实验结果：实验Bode图：由以上图知，系统幅值裕度为0.421dB，相角裕度为175deg，不符合题目要求。方案分析与选择：控制算法效果模拟化数字PID控制控制复杂且效果不一定好，当/Tm（最大时间常数）0.5时，采

14、用常规PID控制难以得到良好的控制Smith纯滞后补偿控制将模型加入到反馈控制系统中，有延迟的一部分用于抵消被延迟了的被控量，无延迟部分反映到调节器，让调节器提前动作，从而可明显地减少超调量和加快调节过程。但是该方式对系统受到的负荷干扰无补偿作用，且控制效果严重依赖于对象的动态模型精度，特别是纯滞后时间达林算法针对多数工业过程具有纯滞后特性，解决/Tm0.5时常规PID控制的难题，能使被控对象输出没有超调量或超调量很小由于被控对象具有纯延时特性，导致控制系统的稳定性降低，过渡过程变坏，且此对象的纯延时时间与对象惯性时间常数T之比为0.5，采用常规PID算法难以获得良好的系统性能，因此使用达林算

15、法。当使用达林算法求取字调节器时，必须求得系统脉冲传递函数。而达林算法目标是设计一个数字调节器，使得整个闭环系统的传递函数为具有纯时延特性的一阶惯性环节，其中，纯时延时间等于工业对象的纯时延时间，因此，我们可以设，其中为系统的期望闭环时间常数。因达林算法也是一种零、极点配置的设计方法，它适用于广义对象含有时延环节、且要求系统没有超调量或超调量很小的应用场合，适当放宽采样周期T和的取值范围或取值为同一数量级时，有利于消除或削弱振铃现象。则可取=0.8。b、达林算法下调节器函数由以上分析知，带纯时延特性的一阶惯性环节如下：其中，0.374，0.7479，K=122.88所设计数字调节器传递函数如下

16、：，式中，0.8。系统脉冲传递函数式中，为采样周期，取L=2，可得。为期望闭环系统的时间常数。则广义工业对象脉冲传递函数为：则所设计数字调节器形式为：因此可求得在单位阶跃输入下，系统期望输出序列：系统的开环传递函数为：C、达林算法下系统伯德图（未加增益）由以上传递函数，通过MATLAB求其伯德图，程序如下： Ts=0.187; num=0 0 0 0.2084; den=1 -0.7916 0 0.2084; Zk=tf(num,den,Ts); margin(Zk)由上图可知，在未加入适当的比例增益，系统的幅值裕度为12.6dB，相角裕度为72.5deg，可见该系统并不符合设计要求，故需要加

17、入适当的比例增益，降低系统的相角裕度。d、达林算法下系统伯德图（加增益）对系统加入一个比例增益为之后，即传递函数变为：通过MATLAB进行仿真，其仿真程序如下： Ts=0.187; num=0 0 0 0.4168; den=1 -0.7916 0 0.2084; Zk=tf(num,den,Ts); margin(Zk)可得伯德图如下：实验bode图：由上图可知，系统此时的幅值裕度为6.57dB，相角裕度为54.1deg，均符合题目要求。e、单位阶跃响应下系统输出通过Simulink仿真，其仿真图如下：得出输出波形如下：由以上图分析可知，系统相对于无滞后系统时会有一定延时时间才能达到无静差，

18、并且系统伴有一定的超调。但系统最后能在较短时间内跟上输出，达到无静差。f、施加阶跃干扰信号通过Simulink仿真，其仿真图如下： 达林算法下系统仿真图（含增益）干扰信号波形：施加干扰信号波形图系统输出波形图：由以上分析可知，达林算法下系统抗干扰性较好，在施加干扰后，能在较短时间回到稳定状态。g、施加随机信号影响通过Simulink仿真，其仿真图如下：达林算法施加随机信号仿真图随机信号波形：系统输出波形图：由以上分析知，达林算法下所设计系统在随机信号输入下，能很好地跟随其变化。综上时，使用达林算法设计的系统在加入增益后，系统的相角裕度，幅值裕度均符合题目所要求，并且该系统具有较好的抗干扰性和稳定性，能在较短时间内使输出达到无静差。五、心得与体会通过本次课程设计，我掌握了计算机控制设计的基本方法。通过独立完成方案分析、整体方案设计、模块设计、控制对象分析、控制算法的选择和控制器的设计与仿真，更好地提高了我系统思考的能力，提高了学以致用、理论联系实际能力，提高了动手的能力。初次接触该设计时，由于一开始T时计算错误以及对无延时系统的分析不合理，导致一直不能达到题目要求。经过查找资料和与同学讨论，最后明白T值是需要使用Matlab计算所得，而并非是一个简单的随机数。而无延时时，一开始并未考虑过多，