计算机测控技术课件

退出授课教师：XXX计算机测控技术第1章绪论

本章主要介绍计算机测控制的一般概念，计算机测控系统的组成、分类、特点及其典型应用和发展。

退出1.1计算机测控系统的组成及特点1.1.1一般概念

模拟式自动控制系统也已达到相当完善的程度。但它的进一步发展受到了限制，在复杂控制规律的实现系统的最优化、可靠性等方面已不能满足更高的要求。现代控制理论的发展给自动控制系统增添了理论工具，而计算机技术的发展为新型控制规律的实现、构造高性能的控制系统提供了物质基础，两者的结合极大地推动了自动控制技术的发展。将模拟式自动控制系统中的控制器的功能用计算机来实现，就组成了一个典型的计算机控制系统，如图1-1所示。因此，计算机控制系统，简单地说，就是采用计算机来实现的工业自动控制(含管理)系统。1.1.2计算机控制系统的组成

(1)传感器和变送器：传感器的主要功能是将被检测的非电量参数转变成电学量，如热电偶把温度变成电压信号，压力传感器把压力变成电信号等等。变送器的作用是将传感器得到的电信号转变成适用于计算机接口使用的标准的电信号(如0-5V，0-10mA)。(2)输入通道：模拟量I：将经由传感器得到的工业对象的生产过程参数变换成二进制代码传送给计算机(A/D)数字量I通道：除完成编码数字输入输出外，还可将各种继电器、限位开关等的状态通过输入接口传送给计算机。(3)输出通道：将计算机输出的数字控制量变换为控制操作拉行机构的模拟信号，以实现对生产过程的控制(D/A)。或将计算机发出的开关动作逻辑信号经由输出接口传送给生产机械中的各个电子开关或电磁开关。(4)执行机构：为了控制生产过程，还需有执行机构。常用的执行机构有各种电动、液动、气动开关，电液伺服阀，交、直流电动机，步进电动机等等。二、软件从功能区分，软件可分为系统软件和应用软件。(1)系统软件是由计算机的制造厂商提供的·．用来管理计算机本身的资源、方便用户使用计算机的软件。常用的有操作系统、开发系统等，它们一般不需用户自行设计编程，只需掌握使用方法或根据实际需要加以适当改造即可。(2)应用软件是用户根据要解决的控制问题而编写的各种程序，比如各种数据采集、滤波程序、信号分析、控制量计算程序、生产过程监控程序等。1.2计算机在控制中的典型应用方式

一、操作指导控制系统(见图1-3)在操作指导控制系统中，计算机的输出不直接用来控制生产对象。计算机只是对生产过程的参数进行采集，然后根据一定的控制算法计算出供操作人员参考、选择的操作方案、最佳设定值等，操作人员根据计算机的输出信息去改变调节器的设定值，或者根据计算机输出的控制量执行相应的操作(如直接改变阀门开度)。三、监督计算机控制系统监督计算机控制(SupervisoryComputerContro1)系统简称SCC系统。在SCC系统中计算机根据工艺参数和过程参量检测值，按照所设计的控制算法进行计算，计算出最佳设定值直接传送给常规模拟调节器或者DDC计算机，，最后由模拟调节器或DDC计算机控制生产过程。SCC系统有两种类型，一种是SCC+模拟调节器，另一种是SCC+DDC控制系统。监督计算机控制系统构成示意图如图1-5所示。四、分级计算机控制系统第2章模拟量输入输出通道

2.1DAC及ADC的性能指标和选择要点2.2DAC及硬件电路2.3ADC及硬件电路2.4多路转换器退出

2.1.1性能指标 1．分辨率:转换位数n2．转换速率（A/D):采样定理稳定时间（D/A)3．输入电压范围（A/D),模拟量输出形式（D/A)4．供电电源5．工作环境2.1DAC及ADC的性能指标和选择要点DAC0832是带有两级数据输人缓冲锁存器的8位D／A转换器。其引脚如图所示。输入寄存器地址：ILE为1，/CS为0,/WR1为0DAC寄存器地址：/WR2为0，/XFER为0二、DAC0832应用接口电路单片机与DAC0832的接口，可根据需要按单级缓冲器方式、二级缓冲器方式联接。 1、单缓冲器连接方式DAC0832以单缓冲器方式与8051的接口电路图所示。２、双缓冲器连接方式采用双缓冲器连接方式时，DAC0832的数字量输入锁存和D／A转换输出分两步完成。首先，将数字量输入到各路D／A转换器的输入寄存器，然后，控制各路D／A转换器，使各路D／A转换器输入寄存器中的数据，同时进入DAC寄存器，并转换输出。所以，在这种工作方式下，DAC0832占用两个I／O地址，输入寄存器和DAC寄存器各占一个I／O地址。这种方式适用于几个模拟量同时输出的系统，每一路模拟量输出需一片DAC0832，也可用多片DAC0832构成多路模拟量同步输出的系统。

图为二路模拟量同步输出的8031系统。图中，1#DAC0832输入寄存器地址为ODFFFH，2#DAC0832输入寄存器的地址为OBFFFH，1#和2#DAC0832的第二级寄存器地址同为7FFFH。两片DAC0832的输出分别接图形显示器的X、Y偏转放大器的输入端。2.3.15G14433(双积分型)的应用一、5G14433A/D转换器的特性及结构

5G14433是上海元件五厂生产的三位半ADC，它是一种双积分型ADC，具有精度高(精度相当于11位二进制ADC)、抗干扰性能好等优点。其缺点是转换速度慢，约1-10次／s。在不要求高速转换的场合，例如温度测控系统中，被广泛采用。5G14433ADC与Motorola公司的产品MCl4433可以互换。5G14433ADC的被转换电压量程为199.9mV或1.999V。转换结果以BCD码的形式分4次输出。(1)VAG：模拟地。(2)Vref：外接基准电压(2V或200mV)输入端。(3)Vx：被测电压输入端。(4)Rl、Rl／C、C：外接积分阻容元件端。外接元件典型值：①当量程为2V时，Cl=0.1luF，R1=470k；②当量程为200mV时，Cl=0.1uF，Rl=27k(5)C01、C02：外接失调补偿电容C0端，C0的典型值为0.1uF。(6)DU：更新转换结果输出的输入端。当DU与EOC连接时，每次转换结果都被更新。(7)CLKI、CLKO：时钟振荡器外接电阻Rc端。Rc的典型值为470k，时钟频率随着Rc的增加而下降。当CLKO为66kHz时，5G14433工作在最佳状况，CLKO最高不超过300kHz。(8)VEE：模拟部分的负电源端，接-5V。(9)Vss：数字地。VDD：正电源端。

(10)EOC：转换周期结束标志输出。每当转换周期结束，EOC端输出一个宽度为时钟周期一半的正脉冲。(11)/OR：过量称标志输出。当1Vx1>Vref时，/OR端输出低电平。(12)DSl—DS4：多路选通脉冲输出端。DSl对应千位，DS4对应个位。如图所示。

(13)Q0—Q3：BCD码数据输出线。其中Q0为最低位，Q3为最高位。当DS2、DS3、DS4选通期间，输出三位完整的BCD码(百位、十位、个位)；但在DSl选通期间，输出端Q0-Q3除了表示千位为0或1外，还表示了转换结果的正负极性和欠量程还是过量程，其含义：

DS1Q3Q2Q1Q0输出结果状态——————————————

11xx0千位数为010xx0千位数为11x1x0输出结果为正1x0x0输出结果为负———————————————————————10xx1输入信号过量程11xx1输入信号欠量程二、5G14433应用接口电路由于5G14433的结果输出是动态的，Q0—Q3和DSl—DS4都不是总线式的，因此必须通过并行接口和MCS—51相连。图为5G14433和8031P1口相连的接口逻辑。例：将上述电路中A／D转换结果存入8031内部RAM的20H、21H单元，按如下格式存放：设8031内部RAM可寻址位10H设为量程错标志位。读取A／D转换的结果。程序清单如下：PINT1:MOVA，P1；读P1口JNBACC．4，PINTl；查DSl

JBACC.0，PERR;判超量程JBACC.2，PLl;判极性SETB07H；为负，符号位置1AJMPPL2PLl：CLR07H；为正，符号位清0PL2：JBACC.3，PL3；判千位SETB04H；千位为1AJMPPL4PL3：CLR04H；千位为0PL4：MOVA，P1JNBACC.5，PL4；查DS2MOVR0，#20HXCHDA，@R0；保存百位PL5：MOVA，P1JNBACC.6，PL5；查DS3SWAPAINCR0MOV@R0，A；保存十位PL6：MOVA，P1JNBACC.7，PL6；查DS4XCHDA，@R0；保存个位RETPERR：SETBl0H；置1量程错标志RET2.3.3ADC0809(逐次逼近型)的应用一、ADC0809的性能指标及结构原理ADC0809是一种典型的A／D转换器，是8位8通道的A／D转换器，转换时间100us，输入电压范围：0-5V。其引脚如下图所示。（1）Vcc:+5V工作电压，GND:数字地（2）REF(+），REF(-):参考电压正负端（3）CLK:时钟信号输入端(内部500KHZ)(4)EOC:转换结束信号输出端。开始转换时为低电平，转换结束时为高电平。(5)数据输出线：D0---D7(6)地址线：START:A/D转换启动信号输入端ALE:地址锁存允许信号输入端OE:输出允许控制端A、B、C通道选择地址输入线二、ADC0809的应用接口电路ADC0809与805l之间的接口电路如下图所示。ADC0809时钟信号由单片机的ALE信号２分频获得（假设fosc=6MHZ)。ADC0809通道地址由P0口的低3位直接与ADC0809的A、B、C相连。ADC0809通道地址:IN0:7FF8HIN1:7FF9H…IN7:7FFFH下面采用查询的方法，对N(N≤8)路模拟信号进行A／D转换，转换后的N个数据顺序存放到起始地址为data_addr数据存区。ADST：MOVR1，#data_addr；数据区首地址指针MOVDPTR，#7FFF8H；指向第１个通道MOVR2，#08H；通道个数LOOP：MOVX@DPTR，A；启动A／D转换LOOP1:NOPJNBP3.3,LOOP1;查询转换结果EOC

MOVXA，@DPTR；读取转换结果MOV@R1，A；结果转存到数据区INCDPTR；指向下一通道INCR1；修改数据区指针DJNZR2，LOOP；若8路未转换完则继续转换初始化程序如下：INT1：SETBIT1；外中断1初始化SETBEASETBEXlMOVDPTR,#7FF8H；启动0809对通道IN0转换MOVA，#00HMOVX@DPTR，A……中断服务程序清单如下：PRINTl：MOVDPTR，#7FF8H；读取A／D转换结果，送缓冲单元30HMOVXA，@DPTRMOV30H,AMOVA，#00H；启动0809对通道IN0转换MOVX@DPTR,ARETI2.4多路转换器一、多路转换器的作用和要求多路转换器又称多路开关。在实际数据处理系统或实际控制系统中，被测量或被控制量往往可能是几路或几十路。对这些回路的参量进行采样和A/D(或D/A)转换时，为了共用A/D(或D/A)转换器，以节省硬件，可以利用多路开关，轮流切换各被测量与A/D转换电路的通路，达到分时转换的目的。二、几种常用的多路开关集成电路芯片

(一)AD7501和AD7503AD7501和AD7503是具有8路输入通道、一路公共输出端的多路开关CMOS集成芯片，由三个地址线(A0、A1、A2)及EN的状态来选择8个通道中的一路。（二）A/D7052A/D7052是双四选一的多路模拟开关，由二个地址线(A0、A1)及EN的状态来选择8个通道中的二路。

（三）CD4051BCD4051B和AD7501相类似，有三位地址信号控制8路模拟信号的开关。与AD7501不同之处是，AD7501是单向的，只能用于多到一。因此，常用于多路模拟量的输入通道。而CD4051B是双向的，既可用于“多到一”，也可用于“一到多”的切换。因此它既常用于多路模拟量输入通道共享A/D转换器，也常用于多路模拟量输出通道，作为D/A转换器到多路模拟量输出的切换开关。第3章开关量输入输出通道

所谓开关量，就是输入信号为具有TTL电平的状态信号，如继电器的吸合与断开、光电门的导通与截止、限位开关、按钮、转换开关、接触器等电器的触点通断，其信号电平只有高、低两种电平。微处理器只要通过对输入信息分析是“1”还是“0”，即可知开关是合上还是断开。如果控制某个执行器的工作状态，只需送出“0”或“1”，即可由操作机构执行。

退出

但是由于工业现场存在着电、磁、振动、温度等各种干扰及各类执行器所要求的开关电压量级及功率不同，所以在接口电路中除根据需要选用不同的元器件外，还需要采用各种缓冲、隔离与驱动措施。如许多时候需要利用计算机控制一些高电压、大电流的负载。由于驱动能力所限，不能使用微机输出线直接驱动这些负载，此时就需要通过功率接口，将输入的低电压、小电流的信号转变成高电压、大电流的信号来进行驱动。另外，有些时候，也需要使用元器件将计算机机与强电进行隔离以防止现场强电磁干扰或工频电压通过输出通道反串到测控系统中，提高抗干扰能力。常用的功率接口有隔离组件光电耦合器、继电器与晶闸管。3.1光电耦合器光电耦合器是一种半导体光电器件，因为它体积小、寿命长、抗干扰能力强等优点，所以应用领域十分广泛。

1.原理：当输入端高电平时，发光二极管发光，光线照射到光敏三极管或光敏二极管上，产生相应的电流，从输出端输出，完成了电—光—电的转换。由于输入电路与输出电路是绝缘的，因此可起到对输出通道的隔离作用。2.应用：当单片机的P1.0口输出为高电平时，发光二极管处于截止状态，不发光，所以光电耦合器无电流输出；当P1.0口输出为低电平时，发光二极管导通并发光，光信号使光敏三极管产生电流输出。通过光电耦合器，可以间接地控制电路的开与断，达到了隔离的目的。3.2继电器当系统需要对负荷比较大的设备进行启动或停止操作时，光电耦合器的负载能力就不能满足要求了。此时便可以选择使用继电器来控制，这是因为继电器的负载能力远大于光电耦合器。继电器是一种电子控制组件，通常用于自动控制系统当中。在工作时是利用较小的电流去控制较大的电流。一般在驱动大型设备时，往往利用继电气作为测控系统输出至输出驱动级之间的第一级执行机构。通过该级继电器输出，可完成从低压电流到高压电流的过渡。

1.工作原理：继电器主要由铁芯、线圈、衔铁、触电、弹簧等组成。当在线圈两端加上一定的电压使线圈中形成一定方向的电流，由于电磁感应现象会产生相应的磁场，该磁场吸引衔铁，使衔铁克服弹簧的拉力与静触电吸合，使电路导通。当线圈两端电压消失后，电流产生的磁场也随之消失，由于没有了吸引力，因此衔铁在弹簧力的作用下将与触电分开，使电路断开。

2.应用：图为带有光电耦合器的继电器接口电路。使用光电耦合器可以进一步提高系统的抗干扰能力。在图中，当P1.0为低电平时，继电器K吸合；而当P1.0输出为高电平时，继电器K释放。带有光电耦合器的继电器接口电路3固态继电器固态继电器(SSR)是近年发展起来的一种新型电子继电器，其输入控制电流小，用TTL、HTL、CMOS等集成电路或加简单的辅助电路就可直接驱动，因此，它适宜于在微机测控系统中作为输出通道的控制元件；它其输出利用晶体管或晶闸管驱动，无触点。与普通的电磁式继电器和磁力开关相比，它具有无机械噪声、无抖动和回跳、开关速度快、体积小、重量轻、寿命长、工作可靠等特点，并且耐冲力、抗潮湿、抗腐蚀。因此，在微机测控等领域中，SSR已逐步取代了传统的电磁式继电器和磁力开关作为开关量输出的控制元件。

固态继电器的内部逻辑图如图所示。固态继电器由光电耦合电路、触发电路、开关电路、过零控制电路和吸收电路5部分构成。这5部分被封装在一个6面体外壳内，成为一个整体，外面有4个引脚(图中的A、B、C、D)。如果是过零型SSR，就包括“过零控制电路”部分，而非过零型SSR则没有这部分电路。3.3晶闸管晶闸管又称可控硅，是一种大功率半导体器件，不但具有整流功能，还具有开关作用。利用晶闸管可以使用较小的功率控制较大的功率。同时，由于晶闸管的响应时间很短，远小于继电器，所以比较适用于高速系统。晶闸管本省体积小、重量轻、效率高，使用十分广泛。常用的晶闸管有单向晶闸管与双向晶闸管两类，还有光控晶闸管等。1.单向晶闸管：单向晶闸管由3个PN结组成。当晶闸管的A、K两端为反向电压时，晶闸管本身不会导通；当晶闸管的A、K两端接正向电压时，如果G端出现信号使G、K两端也为正向电压时，单向晶闸管将导通。需要注意的是，由于单向晶闸管导通后，并不以G端电平作为维持条件，所以G端可以作为触发信号。另外，单向晶闸管没有自我关断能力，在导通之后，只有当A、K两端的电压下降到不足以使其导通时单向晶闸管才会回到截止状态。应用：图示驱动脉冲直流负载的光电耦合—单向晶闸管接口电路。当P1.0=0时，光电耦合导通，三极管T1截止，晶闸管T2导通，直流负载加电；反之，P1.0=1时，由于工作电源是脉冲直流，所以T2关断，直流负载断电。LD为指示灯。2.双向晶闸管：如图1所示。主电极T1、T2无论加正向电压还是反向电压，其控制G的触发信号无论是正向还是反向，它都能被触发导通。一般用于交流控制。图2所示是驱动交流负载的光电耦合—双向晶闸管接口电路。T2为双向晶闸管。其工作原理同上。第4章信号的测量与分析算法一般在计算机应用系统的输入信号中，均含有种种噪声和干扰，它们来自被测信号源本身、传感器、外界干扰等。为了进行准确测量和控制，必须消除被测信号中的噪声和干扰。噪声有两大类：一类为周期性的；另一类为不规则的。前者的典型代表为50Hz的工频干扰。对于这类信号，采用积分时间等于20ms的整数倍的双积分A/D转换器，可有效的清除其影响。后者为随机信号，可以用数字滤波方法予以削弱或滤除。退出4.1数字滤波4.1.1一般概念数字滤波克服了模拟滤波器的不足，它与模拟滤波器相比有以下几个优点：1．数字滤波是用程序实现的，不需要增加设备，所以可靠性高、稳定性好。2．数字滤波可以对频率很低(如0.01Hz)的信号实行滤波，具有灵活、方便、功能强的特点。由于数字滤波具有以上优点，所以数字滤波在微机应用系统中得到了广泛的应用。（一）算术平均值法算术平均值法适用于对一般的具有随机干扰的信号滤波。它特别适于信号本身在一数值范围附近上下波动的情况，如流量、液平面等信号的测量。算术平均值法是要按输入的N个采样数据xi(i=1-N)，寻找y，使y各采样值间的偏差的平方和为最小：算术平均值法对信号的平滑滤波程度完全取决于N。当N较大时，平滑度高，但灵敏度低，即外界信号的变化对测量计算结果y的影响小，当N较小时，平滑度较低，但灵敏度高，应按具体情况选取N，如对流量测量，N可取8-16，对压力等测量，可取N=4。

（二）滑动平均值法

算术平均值法，每计算一次数据，需测量N次。对于测量速度较慢或要求数据计算速率较高的实时系统，该方法是无法使用的。例如A/D数据，数据采样速率为每秒10次，而要求每秒输入4次数据时，则N不能大于2。滑动平均值法只需进行一次测量，就能得到一个新的算术平均值。滑动平均值法采用队列作为测量数据存储器，队列的长度固定为N，每进行一次新的测量，把测量结果放人对尾，而扔掉原来队首的一个数据，这样在队列中始终有N个“最新”的数据。计算平均值时，只要把队列中的N个数据进行算术平均，就可得到新的算术平均值。这样每进行一次测量，就可计算得到一个新的算术平均值。（三）

防脉冲干扰平均值法在工业控制等应用场合中，经常会遇到尖脉冲干扰的现象。干扰通常只影响个别采样点的数据，此数据与其他采样点的数据相差比较大。如果采用一般的平均值法，则干扰将“平均”到计算结果上去，故平均值法不易消除由于脉冲干扰而引起的采样值的偏差。为此，可采取先对N个数据进行比较，去掉其中最大值和最小值，然后计算余下的N-2个数据的算术平均值。即可以滤去脉冲干扰又可滤去小的随机干扰。在实际应用中，N可取任何值，但为了加快测量计算速度，一般N不能太大，常取为4，为四取二再取平均值法。它具有计算方便速度快，存储量小等特点，故得到了广泛的应用。4.2信号分析处理方法一、时间域分析：1.幅值域分析:信号在幅值域上进行一些指标的统计计算:最大值、最小值、平均值、有效值、均方值、标准差、偏态系数、峰态系数。2.相关域分析:(1)自相关函数(2)互相关函数(3)波形分析(4)参数模型分析二、频域分析:(1)自谱分析(2)互谱分析(3)最大熵分析(4)倒谱分析(5)包络谱分析(6)三维谱阵分析三、时频分析：小波变换一、时间域分析1.幅值域分析:是信号在幅值域上进行一些指标的计算。(1)最大值：信号中最大的幅值（峰值）。(2)最大值时间：信号中幅值最大处的时间位置。(3)最小值：信号中最小的幅值谷值。(4)最小值时间：信号中幅值最小处的时间位置。(5)平均值:(6)有效值:(7)均方值:(8)标准差:标准差表示随机振动围绕均值“振荡”的强度，因而它是描述随机振动的重要物理参数。(9)偏态系数：＞0，正偏态；＜0，负偏态；＝0，中心对称。(10)峰态系数：幅域分析中常用的分析方法是振幅概率密度分布,它描述的是数据在幅值域的分布情况。例如轴承摩擦产生随机振动。轴承平稳和遍历随机振动的动态数据序列概率密度函数呈现出正态分布的钟形曲线。四种轴承摩擦状态动态数据的数字特征列于表中。在正常工况下，无摩擦产生，测试值的概率密度函数是标准正态的，，它表示动态数据集中在均值附近，并且变化范围对称于均值分布，如4#轴承；如果工况异常，分布状态便将发生变化。越大，峰态越差，摩擦越严重，如1#轴承。轴承号平均值标准差1#-0.0610.360.00272.502#0.143.23-0.0821.913#0.021.430.00921.364#0.000.56-0.0070.5172.相关域分析(1)自相关函数自相关函数是描述信号x(t)一个时刻的取值与另一个时刻的取值之间的依赖关系计算公式为其中R()为相关函数，x(t)为要分析的信号序列，为时间延迟。工程中通常使用相关系数来描述相关性，更具有对比性和方便性。相关系数函数定义：

式中，为均值，为方差。根据自相关函数图的形状可以判断原信号的性质。比如周期信号的自相关函数仍为周期相同的周期函数。因此，自相关函数可用于检测混于随机噪声中的确定性信号。因为周期信号或任意确定性数据在所有时间上都有其自相关函数，而随机信号则不是。

(2)

互相关函数互相关函数是对两个信号x(t)和y(t)进行分析的，描述x(t)一个时刻的取值与y(t)另一个时刻的取值之间的依赖关系，可以表示为：同样互相关系数的定义如下:互相关函数利用同频检测技术把振动信号中感兴趣的特定频率成分的振幅和相位检测出来，而将非同频信号滤掉。这在对研究平稳随机过程的动态特性和机械系统传递特性上起着十分重要的作用。

二、频域分析

1.功率谱分析

(1)自谱分析就是对一个信号进行频谱分析，FFT为快速傅立叶变换，定义为：正变换：逆变换：其中：，为时域信号，F为频域的频谱密度函数。自谱分析能够将实测的复杂工程信号分解成简单的谐波分量来研究，描述了信号的频率结构，因此，对设备的动态信号做功率谱分析相当于“透视”，从而了解设备各个部分的工作状况，功率谱分析在解决工程实际问题中获得了广泛应用。例如，图为某一大型航空发动机振动功率谱，在其功率谱图上，每个频率分量谱线都对应一定的零部件，利用这种对应关系，可以进行设备状态及故障的识别。(2)互谱分析对两个信号进行互功率谱计算。两路时域信号序列、经过FFT计算得到的幅值谱为X(k)、Y(k)，其乘积X*(k)Y(k)称为互功率谱。其数学定义为：=X*(k)Y(k)式中k为频域序列号，X*(k)表示X(k)的共轭函数。互功率谱函数一般和互相关函数具有同样的作用，但它提供的结果是频率的函数而不是时间的函数，并且具有相位信息，这就大大拓宽了其使用范围。(3)最大熵分析是一种现代谱分析方法，用这种方法可以观测区间以外的数据，以便填补出一个长得多的时区。因此最大熵法具有适合分析短时信号的特点，可用于分析采样点数较少的信号。信号的长度可以小于128点，如64点、32点、20点甚至16点，并且数据点数可以不是2的幂次方，其熵谱理论上可以有任意精细的频率分辨率。最大熵谱分析是利用已有的自相关函数，在保证最大熵的前提下构造出自相关函数的相邻值。重复这个步骤，把自相关函数向两边外推直至正负无穷处，最后作频域变换求得连续的功率谱估计。(4)倒谱分析是近代信号处理科学的一项新技术，它可以处理复杂频谱图上的周期结构。倒谱分析也称为二次频谱分析，是频谱的再次谱分析。它对具有同族或异族谐频以及多成分边频等复杂信号，找出功率谱上不易发现的问题，非常有效，因此常用于振动噪声源识别、机器故障诊断、语音分析、回声剔除等方面的研究。(5)包络谱分析对信号进行采集，确定信号中的共振频率，以该共振频率为中心进行合适带宽的带通滤波，对滤波后的信号进行包络线计算(通常使用Hilbert变换)，得到共振信号的包络线，最后对包络线进行频谱分析得到共振解调谱。由于包络线的频率一般较低，在进行频谱分析后可以再进行细化分析。

(6)三维谱阵分析信号的傅里叶频谱只能反映信号在整个时间过程中的平均频率情况，使用三维谱阵分析则可以反映长数据信号的频率特性随着时间变化而变化的情况。三维谱阵分析是对一个长信号在其不同时间位置上取一定长度的数据点(一般为1024点)分别进行傅里叶谱分析，然后将各次的谱分析结果在三维空间中依次排列起来以表现在不同时间位置上的信号频谱特性。相邻两条谱线所取时间的间隔通常相等，并使各次分析的时间位置在整个时间轴上均匀分布。分析谱线条数越多，则时间间隔越小，越能精确反映信号频率随时间变化的特性，但计算量也越大。

齿轮传动系统的振动加速度信号

概率密度函数分析结果

FFT功率谱分析

三维谱阵

三、时频分析：小波变换1小波变换在振动信号分析中的工程解释与应用

第5章数字控制器的模拟化设计5.1.引言:在微型计算机控制系统中，微型计算机作数字控制器。生产过程中，在大多数情况下被控对象具有连续的特性，而计算机的特性是离散的，由此组成了一个既有连续部分又有离散部分的混合系统。

对于这个系统，可以从两个不同的角度来看。

一方面，由于被控对象、A／D、微型计算机、D／A构成的组合体的输入和输出都是模拟量，所以，该系统可以看成是一个连续变化的模拟系统，因而可以用拉氏变换来进行分析。

Gp(s)被控对象的传递函数，D(s)校正装置的传递函数。

模拟化设计方法

另一方面，由于微型计算机与D／A、被控对象、A／D构成的组合体的输入和输出都是数字量，所以，这一系统又具有离散系统的特性，因而也可以用Z变换来进行分析。

H

(s)和Gp(s)构成了被控对象的广义对象传递函数G

(s)

D(s)微机校正装置的脉冲传递函数，

H

(s)零阶保持器的传递函数。

离散系统的设计方法。数字控制器的两种设计途径或两类设计方法。

一种是在一定条下，把计算机控制系统近似地看成模拟系统，用连续系统的理论来进行动态分析和设计，再将设计结果转变成数字计算机的控制算法，这种方法称为模拟化的设计方法，又称间接设计法。

另一种是，把计算机控制系统经过适当变换，变成纯粹的离散系统，用z变换等工具进行分析设计，直接设计出控制算法，该方法为离散化设计方法，又叫直接设计法。

模拟化设计方法的基本思路是．当系统的采样频率足够高时，采样系统的特性接近于连续变化的模拟系统，因而可以忽略采样开关和保持器，将整个系统看成是连续变化的模拟系统，从而用s域的方法设计校正装置D(z)，再使用s域到Z域的离散化方法求得离散传递函数D(z)。设计的实质是将一个模拟调节器离散化，用数字控制器取代模拟调节器。设计的基本步骤是，根据系统已有的连续模型，按连续系统理论设计模拟调节器，然后，按照一定的对应关系将模拟调节器离散化，得到等价的数字控制器，从而确定计算机的控制算法。

由于用经典方法设计连续系统已为工程技术人员所熟悉，且积累有一定经验，因此模拟化设计方法在实际中被广泛采用。

5.2离散化方法

模拟调节器的离散化方法有多种。由于数字控制器是在线进行控制的，对实时性要求较高，所以，必须采用简单、可靠和足够精确的方法。下面介绍几种常用的离散化方法。

一、差分变换法

模拟调节器若用微分方程的形式来表示，其导数可用差分近似。其变换的基本做法是，把原始的连续校正装置传递函数D(z)转换成微分方程，再用差分方程近似该微分方程。常用的差分近似方法有两种：后向差分和前向差分。为便于编程，离散化只采用后向差分法。

一阶后向差分一阶导数采用近似式:(5-1)(2)二阶后向差分二阶导数采用近似式:(5-2)例5-1求惯性环节的差分方程。

解：由

有化成微分方程:以采样周期T离散上述微分方程得：

即用一阶后向差分近似代替微分得:

代入上式得整理得例5-2求环节的差分方程。解：由，有即化成微分方程代入式(3-1)和式(3-2)得最后得到二、零阶保持器法：

又称阶跃响应不变法，其基本思想是：离散近似后的数字控制器的阶跃响应序列，必须与模拟调节器的阶跃响应的采样值相等。

或者(5-3)其中H(s)即式称为零阶保持器，T为采样周期。

零阶保持器法a)连续系统b)带采样和零阶保持c)等效离散系统例5-3用零阶保持器法求惯性环节的差分方程。

所以整理得

三、双线性变换法

双线性变换法又称突斯汀变换法，它是将s域函数与Z域函数进行转换的一种近似方法。

由Z变换定义，有将和展开成泰勒级数为：对上两式若只取前两项则有：即s近似为：所以，当已知连续传递函数D(s)时，则可计算D(z)：(5-5)，例5-4已知某连续控制器的传递函数，试用双线性变换法求出相应的数字控制器的脉冲传递函数D(z)，其中T=1s。

5.3PID数字控制器的设计

按偏差的比例(Proportional)、积分(Integral)和微分(Derivative)进行控制的调节器。简称PID调节器，是连续系统中技术成熟，应用最为广泛的一种调节器。它的结构简单，参数易于调整。实际运行经验及理论分析证明，运用PID调节器在对相当多的工业对象进行控制能取得满意效果。

一、PID位置式算式的递推形式在模拟调节系统中，PID控制算法的模拟表达式为：(5-6)式中，u(t)为调节器的输出信号；e(t)为偏差信号，它等于给定量与输出量之差；Kp为比例系数；TI为积分时间常数；TD为微分时间常数。

如果采样周期T取得足够小，该算式可以很好地逼近模拟PID算式，因而使被控过程与连续控制过程十分接近。表示的控制算法提供了执行机构的位置u(k)，即其输出值与阀门开度的位置一一对应，所以，通常把式(5-9)称为PID的位置式控制算式或位置式PID控制算法。

连续的时间离散化，即：

积分用累加求和近似得：微分用一阶后向差分近似得：式中，T为采样周期；e(k)为系统第k次采样时刻的偏差值；e(k-1)为系统第(k-1)次采样时刻的偏差值；k为采样序号，k=0,1,2…将式(5-7)和式(5-8)代入式(5-6)，可得离散的PID表达式：(5-7)(5-8)(5-9)如果在式(5-9)中，令(称为积分系数)

(称为微分系数)

则

此即为离散化的位置式PID控制算法的编程表达式。当进行控制时，KP、KI、KD可先分别求出并放在指定的内存单之中，则可实现式(5-9)。(5-10)由式(5-9)可以看出，每次输出与过去的所有状态有关，计算复杂，浪费内存。下面，我们来推导计算较为简单的递推算式。为此，对式(5-9)作如下变动

考虑到第k-1次采样时有使式(5-9)两边对应减去式(5-11)，得

(5-11)整理得：

式中，

(5-12)式(5-12)就是PID位置式算式的递推形式，是编程时常用的形式之一。二、增量式PID控制算法

如果令Δu(k)=u(k)-u(k-1)，则Δu(k)=a0e(k)-a1e(k-1)+a2e(k-2)(5-13)式中a0、a1、a2同式(5-12)中一样。因为在计算机控制中，a0、a1、a2都可事先求出，所以，实际控制时只须获得e(k)、e(k-1)、e(k-2)三个有限的偏差值即可求出控制增量。由于其控制输出对应于执行机构的位置的增量，故式(5-13)常被称为PID控制的增量式算式。实际中，当执行机构需要用控制量的增量来进行控制如控制步进电动机时，可按式(5-13)来计算，在执行机构中常用步进电动机来完成位置的累计。PID增量式控制的原理图如图。PID增量式控制的原理图增量式PID控制算法程序框图增量式PID控制算法与位置式PID控制算法相比较，有如下一些优点：

1)位置式算法每次输出与整个过去状态有关，算式中要用到过去偏差的累加值，容易产生较大的积累误差。而增量式中只需计算增量，算式中不需要累加，控制增量的确定仅与最近几次偏差采样值有关，当存在计算误差或精度不足时，对控制量计算的影响较小，且较容易通过加权处理获得比较好的控制效果。

2)由于计算机只输出控制增量，所以误动作时影响小，且必要时可用逻辑判惭的方法去掉，对系统安全运行有利。

3)手动——自动切换时冲击比较小。

由于上述优点，增量式PID控制算法得到了广泛的应用。

三、速度式的控制算法它采用位置式的导数形式，也就是

即

由于在一般计算机控制系统中，采样周期T是一个常数，所以，速度式与增量式在算法上没有本质区别。在实际应用中，PID数字控制各种算式形式的选择视执行器的形式、被控对象的特性而定。若执行机构不带积分部件，其位置和计算机输出的数字量是一一对应的话(如电液伺服阀)，就要采用位置式算式。若执行机构带积分部件(如步进电动机或步进电动机带动阀门或带动多圈电位器)时，就可选用增量式算式。若执行机构要求速度设定就可选用速度式算式，通常，速度式算式用得较少。第5章数字控制器的直接设计5.1.引言:前一章讨论了数字控制器的模拟化设计方法即间接设计法，它是立足于连续系统的设计，并在计算机上采用数字模拟的方法来实现，其优点是可以充分运用工程设计者熟悉的各种连续系统的设计方法和经验，将它移植到数字计算机上予以实现，从而达到满意的控制效果。但是，模拟化设计方法选定的采样周期必须足够小，除必须满足采样定理外，还要求采样周期的变化对系统性能影响不大。当所选择的采样周期较大或对控制的质量要求较高以及用一台计算机实现多回路控制时，很难满足要求。因而，往往是从被控对象的特性出发，直接根据采样系统理论来设计数字控制器，这种方法称为直接设计法。

如图所示的离散控制系统。Gc(s)为被控对象，G(z)为广义对象的脉冲传递函数，D(z)代表被设计的数字控制器，Φ(z)为系统闭环脉冲传递函数，其表达式。(1)

系统设计的目标是要设计一个数字控制器的脉冲传递函数D(z)，利用它来控制被控对象，达到期望的性能指标。由式(1)可得:

(2)

由式（2）可知，当已知G(z)时，只要根据设计要求选择好Φ(z)

就可求得D(z)。因此，在已知对象特性的前提下，设计步骤为：

1)求得带零阶保持器的被控对象的广义脉冲传递函数G(z)。

2)根据系统的性能指标要求以及实现的约束条件构造闭环z传递函数Φ(z)。

3)依据式(2)确定数字控制器的传递函数D(z)。

4)由D(z)确定控制算法并编制程序5.2.最少拍无差系统的设计

最少拍无差系统，是指在典型的控制输入信号作用下能在最少几个采样周期内达到稳态无静差的系统。q最少拍数

123对最少拍控制系统设计的具体要求如下：

1.准确性要求对典型的参考输入信号，在到达稳态后，系统在采样点的输出值能准确跟踪输入信号，不存在静差。

2.快速性要求在各种使系统在有限拍内到达稳态的设计中，系统准确跟踪输入信号所需的采样周期数应为最少。

3．稳定性要求数字控制器D(z)必须在物理上可以实现且闭环系统必须是稳定的。

由图，系统的误差传递函数为根据准确性要求，系统无稳态误差，而对于时间t为幂函数的典型输入函数其z变换的一般形式为其中A(z)为不包含(1-z-1)因子的关于z-1的多项式1）物理上的可实现性要求:控制器当前的输出信号，只能与当前时刻的输入信号，以前的输入信号和输出信号有关，而与未来的输入信号无关。这就要求数字控制器的输入信号E(z)不能有z的正幂项。2）根据快速性要求：使系统的稳态误差尽快为零，故必然有p=q,F(z)=1单位阶跃：q=1单位速度：q=2单位加速度：q=3必须含有因子(1-z-1)q,即

其中，q为对应于典型输入函数R(z)中分母(1-z-1)因子的阶次。

所以，对于典型的输入来说，有1.单位阶跃输入，即由式(4-11)、(4-12),有即说明系统只需一拍(一个采样周期)，输出就能跟随输入。此时用长除法可得阶跃输入时的输出

2．单位速度输入说明系统只需两拍(两个采样周期)，在采样点上偏差即为零，输出就跟随输入。此时，输出为，即有即单位速度输入时的输出3．单位加速度输入

输入函数即有即说明系统的过渡过程共需三拍(三个采样周期)，此时，输出为

q最少拍数

1235.3最少拍快速有波纹系统设计的一般方法

综合最少拍系统设计中需满足的准确性要求、快速性要求、物理上的可实现性以及稳定性要求。设广义对象的脉冲传递函数为

式中，Gc(s)为对象S的传递函数，当Gc(s)中不含有延迟环节时m=1；当Gc(s)中有延滞环节时，一般m>1。G1(z)是G(z)中不包含单位圆外或圆上的零极点，以及不包含延滞环节z-m的部分。是广义对象在单位圆外或圆上的u个零点，是广义对象在单位圆外或圆上的v个极点。

注意：当G(z)中有z=1的极点时，稳定性和准确性取得一致，可降阶处理。考虑上述条件后，数字控制器中显然不再包含G(z)在单位圆上和圆外的零极点，在物理上具有可实现性。即

1)设定，把G(z)中单位圆上或圆外的极点作为自己的零点，即F1(z)是关于z-1的多项式，且不包含G(z)中不稳定极点ai。2)设定，把G(z)中所有单位圆上或圆外的零点作为自己的零点，即F2(z)是关于z-1的多项式，且不包含G(z)中在单位圆上或圆外的零点bi前面讲过其中，q为对应于典型输入函数R(z)中分母因子的阶次。综合考虑系统的准确性、快速性和稳定性要求，闭环脉冲传递函数中必须选择为

单位阶跃：q=1单位速度：q=2单位加速度：q=3式中，m为广义对象的瞬变滞后；bi为G(z)在z平面单位圆外或圆上的零点，u为G(z)在z平面单位圆外或圆上的零点数；v为G(z)在z平面单位圆外或圆上的极点数。由准确性条件知，包含有因子；由稳定性条件知，

必须包含G(z)在z平面单位圆外和圆上的极点，即有

的因子，其中，ai为G(z)在z平面单位圆外或圆上的非重极点；v为非重极点个数。而例计算机系统，被控对象的传递函数，知：输入为单位速度函数，试设计快速有波纹系统的解：所以带入有显然根据稳定性要求，G(z)中z=1的极点应包含在的零点中，另一方面，对于单位速度输入设计时，由准确性条件知必须包含这一因子。所以，准确性条件中已经满足或包含了稳定性条件要求，故可降一阶设计，可设解得：故闭环脉冲传递函数为所以：此即所求的数字控制器的脉冲传递函数。另外，由求得系统的输出为

由求得数字控制器的输出为

5.4最少拍无波纹系统的设计

波纹产生的原因及设计要求，按上一节所述设计控制器，系统输出在采样点之间存在着波纹，是由控制量输出序列的波动引起的。其根源在于控制变量的z变换有非零的极点。要使控制量u(t)在稳态过程中为零或常数值，必须使多项式即

为关于z-1的有限多项式，因此，此时闭环脉冲传递函数，必须包含G(z)的分子多项式P(z)