算术平均值滤波.ppt

《算术平均值滤波.ppt》由会员分享，可在线阅读，更多相关《算术平均值滤波.ppt（81页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、1第七章 应用程序设计 数据处理程序 控制算法的程序设计本章要点本章要点本章要点本章要点2 微机控制系统中的软件设计，包括选择系统软件及设计应用软件。前者主要选择操作系统和算法语言，后者是根据控制对象的要求编制应用程序。若是有自己设计的单片机或是DSP系统，则不仅要编写应用程序，而且要根据系统要求编制专用的监控程序或系统程序。对单片机而言，则大部分是监控程序。3 所以，为控制，处理具体应用对象而编制的程序就叫应用程序。在微机控制系统中，为使程序紧凑，更重要的，为实时性强，常采用汇编语言编写。应用程序一般包括应用程序一般包括5部分：部分：1.数据处理程序：对计算机的输入输出数据进行处理或变换，如

2、数字滤波，线性化处理及工程量转换等。42.过程监视程序：为监视和控制生产过程，必须有相应的硬件与之配合，比如开关状态的检测及控制操作台，显示器等。应有巡回检测程序及控制操作台服务程序等，以及时反映现场情况和进行人机交互。3.控制算法程序：是为计算机能够按某种检测规律实现控制而编制的程序，如最优控制，PID控制，比值，前馈等。最常用的如PID控制。54.故障自诊断程序：即通过一定的软件去检查控制系统中的相关部件运转是否正常。已有各相关部件较成熟的自检程序，可供参考或借鉴。5.通用子程序：如制表打印格式程序，服务子程序库等。由于时间关系，重点介绍数据处理程序和控制算法程序。67.1 7.1 数据处

3、理程序数据处理程序数据处理程序数据处理程序7滤波：滤波：是将信号频谱中某些成分衰减或滤除。（滤除可用硬件如RC滤波器，也可用软件如数字滤波）高次谐波对电机产生有害影响。脉冲，铁损发热增加8数字滤波：实质是程序滤波。通过一定的计算方法，减少干扰在有用信号中的比重，使得输入计算机待处理的信号量可能是我们所要求的信号。程序限幅滤波 平均值滤波 中值滤波 一阶惯性滤波 复合数字滤波91.无需硬设备，各通道可公用，降低成本。2.各回路间不存在阻抗匹配等问题。故可靠性高，稳定性好。3.对频率很低（如0.01HZ以下）的信号也可滤波，克服了模拟滤波器的缺陷。4.改变滤波方法式滤波器参数只需修改程序，方便，灵

4、活。数字滤波的优点：数字滤波的优点：101.1.程序限幅滤波程序限幅滤波 经验说明，许多物理量的变化都需要一定的时间，相邻两次采样值之间的变化有一定的限度。当采样信号由于随机干扰，如大功率用电设备的启动或停止，造成电流的尖峰干扰或错误检测，以及变送器不稳定而引起的严重失真等现象时，可采用程序判断法进行滤波。程序判断滤波根据滤波方法的不同，可分为限幅滤波和限速滤波两种。11限幅滤波限速滤波 设顺序采样时刻t1，t2，t3所采集的参数分别为Y(1)，Y(2)，Y(3)，那么 显然，这类方法对滤除随机干扰特别有效，程序也很简单。122.平均值滤波平均值滤波 平均值滤波主要对压力，流量等周期脉动的采样

5、值进行平滑加工，但对脉冲性干扰的平滑不理想，它不适用于脉冲干扰比较严重的场合。算术平均值滤波加权平均值滤波法13（1）算术平均值滤波）算术平均值滤波公式为：可见，其实质实质是：把n个采样值相加，求其平均值，作为时刻k的值。显然n越大，结果越准确，但计算时间越长。算术平均值滤波是要寻找一个Y(k)，使该值与各采样值间误差的平方和为最小，14 算术平均值滤波主要用于对压力、流量等周期脉动参数的采样值进行平滑加工，但对脉冲性干扰的平滑作用尚不理想。因而它不适用于脉冲性干扰比较严重的场合。采样次数N的选取，取决于系统对于参数平滑度和灵敏度的要求。随着N值的增大，平滑度将提高，灵敏度将降低。所以要根据具

6、体情况选取n，以便取得满意的滤波效果。15 算术平均滤波对每次采样给出了相同的加权系数1/n。实际应用场合需要加大新采样值在平均值中的比重，为此可采用：加权平均值滤波法l 流量信号，取n=10l 压力信号，取n=5l 温度等缓慢变化信号，取 n=216（2）加权平均值滤波法 运算关系式为：Ci为加权系数应满足：1)Ci=1 2)CnCn-1 C2 C1 因为加权系数常为小数，因此采用浮点运算更为合适。17 上述两种滤波程序，为确定一个采样值，需采集n个数，需要时间长，使得检测时间减慢，所以只能用在信号变化比较缓慢的场合。为了提高检测速度，可采用:递推平均滤波或叫数据滑动平均滤波基本思想：每采样

7、一个新值，即舍弃一个最早的采样值，然后求平均值，作为本次的采样。平均值滤波对周期脉动的干扰有较明确的抑制，但对偶尔出现的脉动干扰信号，滤波效果不理想，应采用中值滤波183.3.中值滤波中值滤波 对被测参数采样n次（n3）,再按大小顺序排列，再取中间值作为本次采样值。此法对防止脉冲干扰比较有效。对缓慢变化的过程变量采用此法有良好的效果，快速变化时要求采样点至少3，所以不适合用去除脉动性干扰。19程序：n=3为例204.一阶惯性滤波一阶惯性滤波 前面几种滤波方法基本上属于静态滤波，主要适用于变化比较慢的参数，但对于快速随机变化的参数，要用动态滤波方法，如：一阶惯性滤波。数学表达式：21 一阶惯性滤

8、波的基本意图是把本次采样值与上次滤波器平均值进行加权平均，因此在转入的过程中，任何快速的干扰均被滤除，仅保留缓慢变化的信号，为此称之为“低通滤波”。T为采样周期，一般事先确定好，Tf=RC为数字滤波器时间常数，通过实际运算可取最优值。方法：不断改变Tf，噪声最低或消除，即为该系统的Tf值。缺点：造成相位滞后，大小与滤波系数有关。225.复合数字滤波 为了进一步提高滤波效果，有时可以把两种或两种以上不同滤波功能的数字滤波器组合起来，构成复合数字滤波器，或称多级数字滤波器。例如：中值均值可由下式表示：若则上式也称为防脉冲干扰的平均值滤波。此外，也可采用双重滤波的方法，即把采样值经过低通滤波后，再经

9、过一次高通滤波，这样，结果更接近理想值，这实际上相当于多级RC滤波。23各种滤波性能的比较各种滤波性能的比较1.滤波效果滤波效果 一般来说，对于变化比较慢的参数，如温度，可选用程序判断滤波及一阶滞后滤波方法。对那些变化比较快的脉冲参数，如压力，流量等，则可选择算术平均和加权平均滤波法，特别是加权平均滤波法。至于要求比较高的系统，需要用复合滤波法。24 在算术平均滤波和加权平均滤波中，其滤波效果与所选择的采样次数N有关。N越大，则滤波效果越好，但花费的时间也越长。高通及低通滤波程序是比较特殊的滤波程序，使用时一定要根据其特点选用。2.滤波时间滤波时间 在考虑滤波效果的前提下，应尽量采用执行时间比

10、较短的程序，若计算机计算时间允许，可采用效果更好的复合滤波程序。25注意：数字滤波在热工和化工过程的DDC系统中并非一定需要，需根据具体情况，经过分析，实验加以选用。不适当的应用数字滤波（例如，可能将待控制的偏差值滤掉），反而会降低控制效果，甚至失控，因此必须给与注意。26线性化处理程序线性化处理程序 在工程实践中，把物理量转化为电信号的传感器，大多带有一定的非线性特征，不便于计算机处理。有的很难找出明确的数学表达式，需根据测量值采取一些特殊的方法进行处理，还有一些参数，相互之间虽然有明确的数学表达式，但计算起来相当麻烦，会占用较多的时间。因此，找出某种既方便又能满足实际功能要求的数据处理办法

11、就是本节要解决的问题。271.计算法 如果被测参数和转换的电信号之间有明确的数学表达式，且该表达式又便于计算，就应发挥微机计算机较强的优势，将其表达式转换为线性关系，再行计算。2.查表法 所谓查表法，就是把事先计算或测得的数据按一定的顺序编制成表格，查表程序的任务就是根据被测参数的值或者中间结果，查出最终所需要的结果。28 微机中表格表格的排列分为有序表和无序表两种。有序表指在表内存放的数有一定的规律（或按大小顺序；或按某种规律），无序表指数据的存放是随意的（或某种函数是无序的）。根据表格的编排方法，查表的方法有3种：顺序查表法计算查表法对分查表法29a.顺序查表法顺序查表法 顺序查表法用于无

12、序表格的查找。因为无序表中各项数据和排列无一定的规律可寻。当需要找某一个数（常被称为关键字），只能从表的第一项开始，逐项比较，直至找到关键字。在程序设计时，已知条件是表格的长度和起始地址。在查找的过程中要不断地修改地址和进行计数。直至找到关键字。若表中无此关键字必须给出提示。30b.计算查表法计算查表法 微机控制系统中使用的线性表，是若干个数据元素的集合，各数据元素在表格中的排列方法及所占用的存储单元个数是一样的。因为要搜索的内容和表格的排列方法有一定的关系，故各元素都严格地按某种顺序排列。在此前提条件下才可以使用计算查表法。31计算查表法的思路是：根据给定的要查找的元素，通过一定的计算，求出

13、该元素所对应的数值的地址。从而可从相应单元中取出所需要的量。计算查表法执行速度是很快的，使用这种查表法的关键是寻找一种计算表地址的公式，其缺点是对表格的要求非常挑剔。32 对分查表法的具体做法是：对于一个字节长度为n的线性表（设该表从小到大排列），先取n/2处的值与待搜索的关键字x进行比较，如果x大于n/2项的内容，下一次取n/2n的中间值，即与 项的内容进行比较。否则取0 n/2项的中间值，即与 项的内容进行比较。依次类推，可以逐次逼近待搜索的关键值所在位置，直至找到。若搜索完毕而未找到需给出提示。c.对分查表法对分查表法33 如果在应用中一些数据或函数能够满足从大到小或从小到大的顺序排列，

14、在这种情况下可使用对分查表法。显然对分查表其速度是很快的，它最多的查找次数为lg2n-1。例如，若查长为2068字节，对分查表最多只需要10次，而若采用顺序查表法，其平均查找次数1024次。343.折线法折线法 在微机控制系统中，有些参数输入和输出的关系只能通过一条曲线来表示，没有一个明确的表达式或一个有序表格。此种场合，只能够采用局部线性化对其进行处理。折线法实际上是线性插值的算法。这种计算程序的复杂程度和执行时间与折线的分段方程有关。35 若采用等距分段，可简化计算过程，但当函数的曲率和斜率变化较大时，引入误差较大。这可以通过细分，使之减小，但占用较多的内存，而且计算量也大。若采用不等距分

15、段，可根据函数曲率变化使分段直线 尽量接近该段曲线，使误差减小，但程序编制却比较麻烦。因此，如何分段视具体情况而定。36标度变换程序设计标度变换程序设计 生产过程中各个参数都有不同的量纲和数值。在计算机控制系统中所有这些物理参数经过变送器转换成A/D转换器能接收的电信号，又经A/D转换成一系列的数码。为了操作人员管理及生产过程的需要，必须把这些数码转换成各种工程量单位，使之便于显示，记录和打印。这种转换被称为标度变换或工程量转换。工程量转换有各种不同类型，它取决于被测参数和测量传感器的类型，要根据实际情况进行设计。371.线性参数的标度变换 线性标度变换时最常用的标度变换方式，其前提条件是被测

16、参数值与A/D转换结果为线性关系。线性标度变换的公式为：38 上式为线性标度变换的通用公式，其中Am、A0、Nm、N0对某一固定的被测参数来说都是常数，不同的参数有着不同的值。为了使程序设计简单，一般把一次测量仪表的下限所A0所对应的A/D转换值置为0，即使N0=0。这样，式可写成：在很多测量系统中，仪表下限值A0=0,此时，其对应的N0=0，上式可进一步简化为：392.非线性参数的标度变换 上述标度变换公式仅适用于线性参数。对非线性参数的标度变换要具体问题具体分析，根据非线性关系求出相应的标度变换公式。流体的流量与被测流体流过节流装置前后产生的压力差的平方根成正比，于是得到测量流量时的标度变

17、换公式为：40 对于流量仪表，一般下限为0，即Q0=0，故上式可简化为：41 若在进行转换时，Q0所对应的数字量N0亦为0，则可进一步步简化为：对于那些不能用运算式表达的物理量（即非线性的传感设备），应根据其转换关系，用插值法或查表法进行标度变换。总之，标度变换与传感设备的性质有关，不总之，标度变换与传感设备的性质有关，不能简单地照搬公式。能简单地照搬公式。427.2 7.2 控制算法的程序设计控制算法的程序设计控制算法的程序设计控制算法的程序设计43 在微机控制系统分类中曾经讲过直接数字控制系统（Direct Digital Control）即DDC系统，这是工业生产中最常用的实时控制系统之

18、一，而DDC系统的基本问题是设计一个数字调节器：即按照实际系统确定控制规律，得出控制算法，编程实现这一算法。而按照偏差进行比例(P)，积分(I)，微分(D)调节的PID算法则是应用最广泛的一种控制算法，本节主要讨论PID算法的数字实现。44准连续准连续PID控制算法控制算法 实际上，对一个过程参数（连续变化量）而言，计算机是断续控制器。因为每一个采样周期T，计算机采样一次，算一次，输出一次，因为输出结果被保持器保持，所以不管此T内被控制量是否变化，控制输出量始终不变。所以，对控制作用来说，不是连续变化的。但是，当采样周期相对于过程时间常数是相当短时，这种断续变化的控制量就近似于连续变化的控制，

19、我们称之为 准连续控制准连续控制45 TTs时，就可用数字调节器代替模拟调节器。T和Ts可比拟时 经用采样控制理论设计数字调节器。具体实现：按照连续系统控制理论设计模拟调节器，然后对该调节器进行数字模拟。即：（同一数量级，10倍或以下）“准连续控制准连续控制”的的条件条件是：是：T Umax e 0，不加，维持原状 e 0，+e，增大 e 则 Kl=0 为e的门限值 则称式（则称式（6）为）为积分分离积分分离PID算式。算式。63其基本思想是：当偏差大于某个规定的门限值，删除积分作用，以便e不至于过大。只有当ej较小时，方引入积分作用，以消除静差。这样，控制量不易进入饱和区；即使进入了，也能较

20、快退出。所以系统的输出特性得到了改善。642.不完全微分的不完全微分的PID算法算法 在标准PID（位置，增量）算法中，即当有阶跃信号输入时，微分项输出急剧增加，容易引起调节过程的振荡，导致调节品质下降。为了克服这一点，又要使微分作用有效，可以采用不完全微分的PID算法。基本思想是：仿照模拟调节器的实际微分调节器，加入惯性环节，以克服完全微分的缺点。传递函数表达式：65不完全微分的PID位置式算式为：则(7)(8)其中：66 完全积分项对于阶跃信号将产生很大的微分输出信号，此信号急剧下降为0，因而引起系统的振荡。而在不完全微分系统中，其微分作用是逐渐下降的，微分输出信号按指数规律逐渐衰减为0，

21、因而系统变化比较缓慢，不易引起振荡，其中延续时间的长短与KD（微分增益）有关。增量式：(9)6768 变速积分的PID算法其它其它PID算法算法 带死区的PID算法微分先行的PID算法 根据系统偏差大小改变积分速度，使其与偏差大小相对应，即偏差大时，积分累加速度减慢积分作用弱；反之积分累加速度加快，增强积分。B死区K死区增益693.纯滞后补偿算法：史密斯预估器纯滞后补偿算法：史密斯预估器 在工业控制中，许多控制对象往往具有纯滞后的性质。所谓纯滞后就是由于调节作用点或扰动作用点和被调参数相隔一定的距离，而导致被调参数的变化落后与调节作用或扰动的发生和变化。传递函数：比一般一阶惯性环节的前行通道上

22、，多了（乘）一个纯滞后环节）70解决办法：消除es使对象变为纯一阶惯性环节加一个史密斯预估计补偿纯滞后71可以看作两个环节，分别离散化：计算机实现：广义对象成为一阶惯性环节72另一种方法：采样采样-PID 这种补偿的方法对于模拟调节器几乎无法实现，用计算机则能很容易的实现。一阶惯性滤波的推到过程一样T是采样周期73PID调节器参数的选择调节器参数的选择 PID调节器设计的前提是：T Ts 准连续系统可仿照连续系统方法对PID参数整定稳定边界法；反应曲线法；衰减曲线法。741.稳定边界法稳定边界法 稳定边界法又称临界比例度法，即在生产工艺许可的情况下，用试验方法找出，当比例调节系统的被调量作等幅

23、振荡（即达到了稳定边界时的临界比例度Pm）时，按经验公式求出调节器的整定参数。调节器整定数据752.反应曲线法反应曲线法对象飞升特性 根据飞升特性曲线定出几个能代表该调节对象动态特性的参数，然后可直接按这个数据定出调节器最佳整定参数。76对于P调节器经验公式经验公式对于PI调节器对于PID调节器 根据代表对象动特性的参数K、T、，可以利用下面的经验公式，计算对应于衰减率 =0.75下调节器的最佳整定参数。77【例例】在某一蒸汽加热器的自动调节系统中，当调节器的输出从6mA改变到7mA时，温度记录仪的指针从85.0升到87.8，从原来的稳定状态达到新的稳定状态。仪表的刻度为50100，调节器输出

24、010mA，并通过飞升曲线测出 =1.2min，T=2.5min。如果采用PI和PID调节规律，试确定 0.75时的PID参数。783.衰减曲线法衰减曲线法 先把调节器改成比例作用，由手动投入自动，在达到稳定情况后，适当改变给定值，达到规定衰减比后，记下此时的比例度及周期，然后按表求得其它调节规律时的整定参数。调节器整定数据对象衰减比为1/10时的整定数据衰减振荡过渡过程79三种整定方法的比较三种整定方法的比较 稳定边界法：适用于一般的流量、压力、液面和温度调节系统，但实验时需要注意调节过程是否超限。反应曲线法：实验容易掌握，做实验所需时间比其它方法短些，但飞升特性曲线不太容易得到。衰减曲线法：安全，且容易掌握，能适用于各种类型的调节系统，但实验很费时间。80试凑法试凑法-瞎子爬山法，摸索可至山顶，杂凑法瞎子爬山法，摸索可至山顶，杂凑法 步骤是：先调P，只用P调节，有静差时，可考虑后加I，满意时，可用PI调节，动特性不好可考虑D，加D后反复试，满意为止。参数整定找最佳，从小到大顺序查。先是比例后积分，最后再把微分加。曲线振荡很频繁，比例度盘要放大。曲线漂浮绕大湾，比例度盘往小扳。曲线偏离回复慢，积分时间往下降。曲线波动周期长，积分时间再加长。曲线振荡频率快，先把微分降下来。动差大来波动慢。微分时间应加长。理想曲线两个波，前高后低 4比1。一看二调多分析，调节质量不会低。81