logo资料库

新型PID控制及其应用(共六讲).pdf

第1页 / 共31页
第2页 / 共31页
第3页 / 共31页
第4页 / 共31页
第5页 / 共31页
第6页 / 共31页
第7页 / 共31页
第8页 / 共31页
资料共31页,剩余部分请下载后查看
第一讲_PID控制原理和自整定策略
第二讲_自适应PID控制
第三讲_智能PID控制
第四讲_模糊PID控制
第五讲_多变量PID控制器
第六讲新型PID控制器及其发展
·06· 工业仪表与自动化装置            1997 年第 4 期 新型 P ID 控制及其应用 第一讲 P ID 控制原理和自整定策略 陶永华 华东冶金学院 马鞍山: 243002   编者按 随着控制仪表的发展, 新的控制技术不断出现, 先进的控制算法不断产生和 发展。 本刊以前曾发表不少这类的专题文章。 为了使读者能系统了解这个领域的内容和 产品发展概况, 我们特请陶永华教授以讲座形式作全面介绍。讲座共分六讲: 第一讲, P ID 控制原理和自整定策略; 第二讲, 自适应 P ID 控制; 第三讲, 智能 P ID 控制; 第四讲, 模糊 P ID 控制; 第五讲, 多变量 P ID 控制; 第六讲, 新型控制器产品发展概述。   P ID 控制是最早发展起来的控制策略之 一, 由于其算法简单、鲁棒性好、可靠性高等优 点, 被广泛应用于工业过程控制。当用计算机实 现后, 数字 P ID 控制器更显示出参数调整灵 活、算法变化多样、简单方便的优点。 随着生产 的发展, 对控制的要求也越来越高, 随之发展出 许多以计算机为基础的新型控制算法, 如自适 应 P ID 控制、模糊 P ID 控制、智能 P ID 控制等 等。本讲座共分 6 讲, 将着重介绍这些新型 P ID 控制原理、方法及其应用。 我们期待着把 P ID 控制提到一个新的水平。 1 P ID 控制原理 1 模拟 P ID 控制器 1 模拟 P ID 控制系统原理框图如图 1- 1 所 示, 系统由模拟 P ID 控制器和受控对象组成。 图 1- 1 模拟 P ID 控制系统原理框图   P ID 控制器根据给定值 r (t) 与实际输出值 c (t) 构成的控制偏差:     e (t) = r ( t) - c ( t) (1- 1) 将偏差的比例 (P )、积分 ( I) 和微分 (D ) 通
1997 年第 4 期            工业仪表与自动化装置 ·16· 过线性组合构成控制量, 对受控对象进行控制。 其控制规律为: u ( t) = K p [ e ( t) + 1 T I t 0 e ( t) d t + T D ] d e ( t) d t (1- 2) 或写成传递函数形式: U (S ) E (S ) = K P (1+ G (s) = + T DS ) (1- 3) 式中, K P 为比例系数, T I 为积分时间常数, T D 为微分时间常数。 1 T IS 简单说来, P ID 控制器各校正环节的作用 是这样的: ●比例环节: 即时成比例地反应控制系统 的偏差信号 e ( t) , 偏差一旦产生, 控制器立即产 生控制作用以减小误差。 ●积分环节: 主要用于消除静差, 提高系统 的无差度, 积分作用的强弱取决于积分时间常 数 T I, T I 越大, 积分作用越弱, 反之则越强。 ●微分环节: 能反应偏差信号的变化趋势 (变化速率) , 并能在偏差信号值变得太大之前, 在系统中引入一个有效的早期修正信号, 从而 加快系统的动作速度, 减小调节时间。 1 2 数字 P ID 控制器 当计算机实现 P ID 控制时, 首先必须将上 述 P ID 控制规律的连续形式变成离散形式, 然 后才能编程实现。P ID 控制器控制算法的离散 形式为: u (k ) = K P {e (k ) + T T I e (k - 1) } 或 k ∑ j = 0 e ( j ) + T D T [ e (k ) - (1- 4) u (k ) = K P ·e (k ) + K I∑ e ( j ) + K D [ e (k ) k j= 0 - e (k - 1) ] 式中 T ——采样周期 (1- 5) k ——采样序号, k = 0, 1, 2, … u (k ) ——第 k 次采样时刻的计算机输出 值 由 Z 变换的性质: Z [ e (k - 1) = Z - 1E (Z ) Z [∑ k e ( j ) = E (Z ) 式 (1- 5) 的 Z 变换式为: j= 0 (1- Z - 1) U (Z ) = K P E (Z ) + K IE (Z ) K D [ E (Z ) - Z - 1E (Z ) ] (1- 6) 由式 (1- 6) 便可得到数字 P ID 控制器的 Z 传 递函数: (1- Z - 1) + G (Z ) = U (Z ) E (Z ) = K P + K I (1- Z - 1) + K D (1- Z - 1) (1- 7) 或者 G (Z ) = K P (1- Z - 1) + K I+ K D (1- Z - 1) 2 1- Z - 1 (1- 8) 数字 P ID 控制器如图 1- 2 所示。 图 1- 2 数字 P ID 控制器框图   由于计算机输出的 u (k ) 直接去控制执行 机构 (如阀门) , u (k ) 的值和执行机构的位置 (如 阀门开度) 是一一对应的, 所以我们通常称式 (1 - 4) 或式 (1- 5) 为位置式控制算法。 当执行机构需要的是控制量的增量 (例如 去驱动步进电机) 时, 可由式 (1- 5) 导出提供增 量的 P ID 控制算式。 根据递推原理可得: u (k - 1) = K P ·e (k - 1) + K I∑ e ( j ) k- 1 j= 0 e (k ) ——第 k 次采样时刻输入的偏差值 + K D [ e (k - 1) - e (k - 2) ] (1- 9) K I= K D = K P T T I K P T D T ——称为积分系数 ——称为微分系数 用式 (1- 5) 减式 (1- 9) 可得: u (k ) = K P [ e (k ) - e (k - 1) + K Ie (k ) + K D [ e (k ) - 2e (k - 1) + e (k - 2) ]
ΣÙ 工业仪表与自动化装置            1997 年第 4 期 e (k ) - (1—10) G (S) = K P (1+ + T DS) 有 Zieg ler- N icho ls 整定公式: 1 T IS K · 2T P K P = 1 T I= 2 T D = 0 实际应用时, 通常根据阶跃响应曲线 (图 2 (2—2) 5 - 1) , 人工测量出 K、T P、 参数, 然后按式 (2- 2) 计算 K P、T I、T D。用计算机进行辅助设计时, 一是可以用模式识别的方法识别出这些特证参 数; 一是可用曲线拟合的方法将阶跃响应数据 拟合成近似的一阶惯性加纯延迟环节的模型。 图 2- 1 阶跃响应曲线 1 2 2  IST E 最优设定方法 庄敏霞与 A therton 针对各种指标函数得 出了最优 P ID 参数整定的算法, 考虑下面给出 的最优指标通式。 ∞ ) = J n ( (2- 3) 这里 e ( t) 为进入 P ID 控制器的误差信号。 , t) ]2d t [ tne ( 0 根据设定点信号的最优自整定算法, 对式 (2- 3) 中给出的最优指标, 着重考虑 3 种情况, 即 n = 0, 简记作 ISE ( in teg ral squared erro r) 准则; n = 1, 简记作 IST E 准则; n = 2, 简记为 IST 2E 准则。 若已知系统的数学模型如式 (2- 1) 给出, 则对典型 P ID 结构可以建立经验公式: K P = a 1 K ( T P ) b1 T I= T P a 2+ b2 ( T D = a 3T P ( T P ) T P ) b3 (2- 4) Σ ·26· = K P e (k ) + K Ie (k ) + K D [ e (k - 1) ] 式中  e (k ) = e (k ) - e (k - 1) 式 (1- 10) 称为增量式 P ID 控制算法。 将 上式整理, 合并后得: u (k ) = A ·e (k ) - B ·e (k - 1) + C ·e (k - 2) (1- 11) ) 式中 A = K P (1+ T T I B = K P · (1+ C = K P T D T + T D T 2T D T ) 它们都是与采样周期, 比例系数、积分时间常 数、微分时间常数有关的系数。 增量式 P ID 和位置式 P ID 实质是一样的, 但增量式比位置式有许多优越之处: ● u (k ) 只与 k、k - 1、k - 2 时刻的偏差有 关, 节省内存和运算时间。 ●每次只作 u (k ) 计算, 而与位置式中积 分项∑e ( j ) 相比, 计算误差影响小。 ●若执行机构有积分能力 (如步进电机) , u (k ) , 即执行机构的变 则每次只需输出增量 化部分, 误动作造成的影响小。 ●手动—自动切换时冲击小, 便于实现无 扰动切换。 在具体应用数字 P ID 控制器时, 可采用一 些改进算法、如积分分离 P ID 控制算法、不完 全微分 P ID 控制算法, 带死区的 P ID 控制、变 速积分 P ID 算法等。 2 自整定 P ID 控制策略   P ID 控制器中 3 个参数的确定问题是实际 中的设计问题。 下面介绍几种较成熟的自整定 方法: 2 2 1 Z ieg ler- N icho ls 设定方法 受控对象大多可近似用一阶惯性加纯延迟 1 经验公式法 1 环节来表示, 传递函数为: s G P (s) = K e- 对于典型 P ID 控制器: (T PS + 1) (2- 1)
Ù 对不同的 T P 范围, 可以得出 (a, b) 参数 表如表 2- 1 所示。 由表中给出的 P ID 参数设 置可以通过M A TLAB 来简单地实现。 表 2- 1 设定点 P ID 控制器参数表 T P 范围 准则 a1 b1 a2 b2 a3 b3 0 1—1 1 1—2 ISE ISTE TS I2E ISE ISTE IST 2E 1 048 1 042 0 968 1 154 1 142 1 061 - 0 897 - 0 897 - 0 904 - 0 567 - 0 579 - 0 583 1 195 0 987 0 977 1 047 0 919 0 892 - 0 368 - 0 238 - 0 253 - 0 220 - 0 172 - 0 165 0 0 489 888 0 0 385 906 0 0 316 892 0 0 490 708 0 0 384 839 0 0 315 832 1 2 3 临界灵敏度法 当已知系统的临界比例增益 K C 和振荡周 期 T C 时, 也可用经验整定公式来确定 P ID 控 制器的参数。 例如: K P = 0 T I= 0 T D = 0 6K C 5T C 125T C (2- 5) 特征参数 T C 和 K C, 一般由系统整定试验确定。 或者用频率特性分析算法据受控过程 G P (S ) 直 接算得 K C 和 T C, 即由增益裕量 gm 确定 K C, 由 相位剪切频率 C 确定 T C: T C= 2 C K C= 10 (gm 20) 2 2 仿真试验法 (2- 6) 当受控过程的模型已知时, P ID 参数可通 过数字仿真试验来确定。常见的有两种方式: 半 自动的和全自动的。所谓半自动的方式, 就是由 人来设置并调整 P ID 参数, 由计算机来仿真系 统动态特性和计算系统的性能指标。 所谓全自 动的方式就是调整 P ID 参数 的任务也由计算 机来完成, 一般用最优化算法来整定 P ID 参 数。 对于单参数的 P 控制器常用黄金分割算 法。 对于多参数的 P I 和 P ID 控制器常用单纯 1997 年第 4 期            工业仪表与自动化装置 ·36· 形优化算法。优化性能指标常选 IA E 或 ISE, 一 个全自动的 P ID 参数优化 CAD 系统如图 2—2 所示。 图 2- 2 P ID 控制器参数优化 CAD 程序框图 m -Ha gglund 自整定 P ID 3  stro 控制结构   由前面的讨论可知, 若测出了系统的一阶 模型 (式 2—1) 或得出了系统的振荡频率 C 和 增益 K C, 则可以容易地设计出 P ID 控制器。以 往要想求出系统的这些特征参数, 需要使用离 线的方法来进行, 即首先通过试验测出系统的 特征参数, 然后再根据这些参数设计一个合适 的 P ID 控制器, 最后再将此控制器应用到原系 统的控制中。若系统的参数发生变化, 则应该再 重新开始这一过程。 m H a stro gg lund 提 出 了 一 种 继 电 型 P ID 自整定控制结构, 该方案的基本想法是在 控制系统中设置两种模态; 测试模态和调节模 态, 在测试模态下由一个继电非线性环节来测 试系统的振荡频率和增益, 而在调节模态下由 系统的特证参数首先得出 P ID 控制器, 然后由 此控制器对系统的动态性能进行调节。 如果系 统的参数发生变化时, 则需要重新进入测试模 态进行测试, 测试完成之后再回到调节模态进 行控制。 继电型 P ID 自整定控制结构如图 3—1 所 示。从图中可以看出, 两个模态之间的切换是靠
·46· 工业仪表与自动化装置            1997 年第 4 期 开关来实现的。 4 伪微分反馈控制策略   该控制策略如图 4—1 所示。   在算法中反馈环节为 K D (S ) = K dS + K P , 并引入了积分环节 K I S 来确保系统的闭环响 应中不存在静态误差, 系统外反馈回路的 K H 用来调节系统输出的幅值。 为了保证进入对象 环节的信号 m 2 ( t) 不至于过大, 在受控对象前 还加入了一个饱和限幅的非线性环节。 这样的 控制器有着比传统的 P ID 控制器更明显的优 越性。 因为它的前向通路是由单一的积分环节 驱动的。 而微分环节是对较为平滑的输出信号 作用的, 所以并不存在象传统 P ID 控制中微分 突变的现象, 有着较好的控制效果。 图 3- 1 继电型 P ID 自整定控制结构 5 基于启发式推理技术的参数自 整定   这种控制器参数自整定技术, 完全摒弃了 建模和基子模型的运算, 而是通过对闭环性能 的观测, 直接进行整定参数的操作, 现已成为一 种重要的在线自整定策略。 其商品化产品已从 独立地嵌入单回路控制器中, 发展到应用于计 算机控制系统中。P ID 参数整定规则是基于有 经验的操作者的直观判断和行之有效的整定公 式。 这种控制器在启动开环阶跃测试时预置 P ID 参数, 利用相关的模式识别技术在连续闭 环操作期间修改其参数值。 其中大多是在负荷 扰动或设定值变化后, 按照过程误差响应的振 荡周期 T、超调和衰减率, 整定其 P ID 参数, 如 图 5—1 所示。   还有一类不是在每一扰动或设定值变化后 进行整定, 而是执行迭代技术 (下转第 46 页) 图 4- 1 伪微分反馈控制策略
·64· 工业仪表与自动化装置            1997 年第 4 期 面, 其中包括原始参数画面、炉体总貌、棒状图、 历史趋势图、各段加热曲线及炉压曲线图。通过 这些画面, 操作人员可掌握系统的全局状况, 还 可以检索某个参数的历史记录和变化趋势。 操 作人员可以通过操作员键盘迅速调出所需的任 一画面。 3 2 系统带有完善的故障诊断与保护功能 当系统检测或调节器本身出现故障时, 系 统能自动脱离串级进行强迫手动, 指示灯闪烁, 可通过上位机或编程器获得出错信息以便及时 处理 本系统投入运行近一年时间, 运行可靠稳 定, 各项控制功能和指标完全满足工艺要求。 图 2- 5 二加热段温度设定曲线 (上接第 64 页)  在线搜索出最优整定参数。这 时, 控制系统连续监视被控过程的性能, 然后按 照专家整定规则, 修正当前的整定参数。 图 5- 1 直接确定整定参数的模式识别试验曲线 6 结束语   本讲介绍了 P ID 控制原理和自整定技术。 虽然, 大多数有实用价值的自整定技术是用于 单回路控制器, 只有少数用于分散型控制系统, 但其发展趋势却是以参数自整定算法软件形 式, 固化在计算机控制系统或 PL C 中, 使其作 为一种通用控制装置的内含部件。 对于对象为 时变系统或环境干扰变化较大时, 可采用自适 应、智能、模糊等方法实时对 P ID 参数进行整 定。 有关这些内容, 将在下面各讲分别进行介 绍。
·05· 工业仪表与自动化装置            1997 年第 5 期 新型 P ID 控制及其应用 第二讲 自适应 P ID 控制 陶永华 华东冶金学院 马鞍山: 243002 1 引言   所谓自适应控制 (A dap tive con tro l) 就是 在控制对象未知的情况下, 或者控制对象的参 数发生变化时, 调整控制器的控制方法或参数, 使控制系统达到预定的控制品质。 自适应控制的研究是从 50 年代开始的, 由 于计算机的迅速发展及随机控制理论、系统辨 识等学科的发展, 大大促进了自适应控制的研 究。 自适应控制与 P ID 控制器相结合, 形成了 所谓自适应 P ID 控制或自校正 P ID 控制技术 (人们统称为自适应 P ID 控制)。 自适应 P ID 控制具有自适应控制与普通 P ID 控制器两方面的优点: 首先, 它是自适应控制器, 就 是说, 它有自动辨识被控过程参数、自动整定控制 器参数、能够适应被控过程参数的变化等一系列优 点; 其次, 它又具有常规 P ID 控制器结构简单、鲁棒 性好、可靠性高、为现场工作人员和设计工程师们 所熟悉的优点。自适应P ID 控制所具有的这两大优 势, 使得它成为过程控制的一种较理想的自动化装 置, 成为人们竞相研究的对象和自适应控制发展的 一个方向。 自适应 P ID 控制器可分为 5 大类, 一类基 于被控过程参数辨识, 统称为参数自适应 P ID 控制器, 其参数的设计依赖于被控过程模型参 数的估计。 另一类基于被控过程的某些特征参 数, 诸如临界振荡增益 K C、临界振荡频率 等。这种类型的自适应 P ID 控制没有一个统一 名称, 我们姑且称为非参数自适应 P ID 控制 器。 其参数的设计直接依赖于过程的特征参数 C 和一些工程上常用的经验整定规则。 如果按控制器参数设计的原理来分, 自适 应 P ID 控制器又可分为五大类, 它们是: 极点 配置自适应 P ID 控制器、相消原理自适应 P ID 控制器、基于经验规则的自适应 P ID 控制器、 基于二次型性能指标的自适应 P ID 控制器和 智能或专家自适应 P ID 控制器。 本讲着重介绍极点配置自适应 P ID 控制 和一些具有应用前景的 P ID 自适应控制系统。 2 极点配置自适应 P ID 控制   极点配置自适应控制算法由W ellstead 等 人在 1979 年首先提出, 继而由 m 和 w it tenm ark, V ogel 和 Edgar, E llio tt 等人改进和 深化, 成为自适应控制中的一个重要组成部分, W ittenm ark 和 m 等人在此基础上提出 了极点配置自适应 P ID 控制算法。 stro stro 设 受 控 过 程 的 数 学 模 型 可 用 CA RM A M oving R eg ressive A u to ( Con tro lled A verage) 受控自回归滑动平均模型描述:  A (Z - 1) y (t) = Z - dB (Z - 1) u ( t) + C (Z - 1) ( t) (2—1)   其中 y ( t)、u ( t)、 (t) 分别为被控对象的输 出、输入和不可测的扰动噪声; Z - 1为滞后一步 算子, 即 Z - 1y ( t) = y ( t- 1) ; d 为滞后步数; A (Z - 1)、B (Z - 1)、C (Z - 1) 均为 Z - 1的多项式:   A (Z - 1) = 1+ a1Z - 1+ …+ anaZ - na B (Z - 1) = b 0+ b 1Z - 1+ …+ bnbZ - nb C (Z - 1) = 1+ c1Z - 1+ …+ cncZ - nc (b0≠0) (2—2)
1997 年第 5 期            工业仪表与自动化装置 ·15·    其中 na、nb、nc 分别为多项式 A (Z - 1)、B (Z - 1)、C (Z - 1) 的阶数; 一般 d 为已知, 但系数 a i (i= 1, 2, …, na) 和 bj ( j= 1, 2, …, n b) 为未知参 数, 需要在线辨识。 典型的闭环计算机控制系统结构如图 2— 1 所示。   其中 T (Z - 1) 称为期望特征多项式, 由设计 者根据实际工况和性能指标要求确定。 图 2—1 系统可直接写出自校正控制器输 出 u ( t) 的表达式: H (Z - 1) F (Z - 1) y r (t) - G (Z - 1) F (Z - 1) y (t) (2—5) u ( t) = 在数字系统中, 通常采用带数字滤波器的 P ID 控制器算法: u (Z - 1) = g 0+ g 1Z - 1+ g 2Z - 2 (1- Z - 1) (1+ f 1Z - 1) e (Z - 1) ) (2—6) (2—7) g 0= K P (1+ T T I 其中 g 1= - K P (1+ g 2= K P T D T + T D T 2T D T )   K P、T I 和 T D 分别为 P ID 控制器的比例系 数、积分时间和微分时间, T 为采样周期, f 1 为 滤波器参数。 为了将极点配置自校正控制器转换成增量 型 P ID 控 制 器, 根 据 式 ( 2—6) , 本 文 对 方 程 (2—5) 采用如下形式: F (Z - 1) = (1- Z - 1) (1+ f 1Z - 1)   G (Z - 1) = g 0+ g 1Z - 1+ g 2Z - 2 (2—8) (2—9) 图 2—1 典型的闭环计算机控制系统结构   系统闭环方程为: y (t) = Z - dB (Z - 1)H (Z - 1) y r (t) + F (Z - 1)C (Z - 1)V (t) A (Z - 1) F (Z - 1) + Z - dB (Z - 1)G (Z - 1) (2—3) 由式 (2—3) 可见, 按增广型自校正闭环极点配 置要求, 闭环特征多项式应满足: A (Z - 1) F (Z - 1) + Z - dB (Z - 1) G (Z - 1) = C (2—4) (Z - 1) T (Z - 1)   根据方程 (2—5) , 稳态时应有 y ( t) = y r ( t) , 系统才能达到伺服跟踪。 则应有: 图 2—2 极点配置自校正 P ID 控制器结构图
分享到:
收藏