CN100428090C - 工业多输入输出过程的分布式pi和pid控制器定量整定方法 - Google Patents

工业多输入输出过程的分布式pi和pid控制器定量整定方法 Download PDF

Info

Publication number
CN100428090C
CN100428090C CNB2006101172782A CN200610117278A CN100428090C CN 100428090 C CN100428090 C CN 100428090C CN B2006101172782 A CNB2006101172782 A CN B2006101172782A CN 200610117278 A CN200610117278 A CN 200610117278A CN 100428090 C CN100428090 C CN 100428090C
Authority
CN
China
Prior art keywords
controller
sub
tjk
theta
matrix
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CNB2006101172782A
Other languages
English (en)
Other versions
CN1936737A (zh
Inventor
刘媛媛
张卫东
蔡云泽
高守伟
陈凤祥
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Jiaotong University
Original Assignee
Shanghai Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CNB2006101172782A priority Critical patent/CN100428090C/zh
Publication of CN1936737A publication Critical patent/CN1936737A/zh
Application granted granted Critical
Publication of CN100428090C publication Critical patent/CN100428090C/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Feedback Control In General (AREA)

Abstract

一种工业多输入输出过程的分布式PI和PID控制器的定量整定方法,用于工业过程控制技术领域。本发明首先根据辨识得到的被控对象的传递函数矩阵模型,利用有效相对耦合分析方法建立有效等价对角传递函数矩阵模型,然后根据建立的有效等价对角传递函数矩阵模型提出期望的闭环传递函数矩阵,进而获得分布式PI(或PID)控制器的定量整定公式。本发明充分考虑了系统耦合的动态特性,整定过程简单直观,数据运算量小,可方便应用于高维多输入输出过程。采取本发明方法可以实现在线单参数化调节各子控制器,因而操作简捷和方便,而且能够达到明显改进的控制效果。

Description

工业多输入输出过程的分布式PI和PID控制器定量整定方法
技术领域
本发明涉及的是一种用于工业过程控制技术领域的方法,特别是一种工业多输入输出过程的分布式PI和PID控制器的定量整定方法。
背景技术
在过程工业界,从20世纪40年代开始,采用PI和PID控制策略的简单反馈控制回路已成为过程控制的核心系统。目前,PI和PID控制器仍广泛应用,即便是在大量采用DCS控制的最现代化的装置中,这类回路仍占总回路数的80%-90%。然而,单回路PI和PID控制并不能适用于所有的过程,当所涉及的被控过程具有强耦合性、不确定性和大纯滞后等特征时,这类控制方法无法直接应用。为此,出现了针对具有复杂特性的多输入输出过程的PI和PID控制器矩阵设计方法。其中,最常用的方法是将多变量系统解耦为多个单变量系统,然后应用单变量系统的控制方法实现分布式(Decentralized)控制,也称为多环(Multiloop)控制结构。由于具有结构简单,成本低廉以及各闭环独立工作性强等优点,在工业生产实践中得到广泛采用。然而这种控制结构的稳定性容易受到实际被控过程不确定性的破坏,因此对该种控制结构中的多个标量PI(比例积分)或PID(比例积分微分)控制器的整定和调节要求比较高,为了保证这种控制结构的工作稳定性,通常将多个PI(或PID)控制器调节得比较松弛,但是由此造成控制系统的工作效率比较低,使得原材料和能源消耗较大,不利于经济生产和运行。
经对现有技术文献的检索发现,具有代表性的是国际著名教授I.L.Chien在文献《A Simple Multiloop Tuning Method for PID Controllers with NoProportional Kick》(无比例跳动的简单分布式PID控制器整定方法,发表在Industrial Engineering Chemical Research化工工程研究刊物,1999,38,1456-1468.)中,提出根据系统稳态时的耦合特性整定分布式闭环PID控制器的方法,这里简称之为Chien的方法。该方法简单,直观,方便实际工程应用,但是由于只考虑了系统耦合的稳态特性,可获得的改善的控制效果有限。需要指出,其他近期发表的文献给出的分布式闭环PI(或PID)控制器的整定方法虽然能够相对于实践中常用的简易整定方法取得显著的改进效果,但大多只适用于低维多变量过程。当这些方法被应用于高维多变量过程时,往往需采用数值化寻优设计方法来整定控制器,数据运算量比较大,不便于在线调节和设定,并且所用到的相关专业理论知识较多,不便于被工程技术人员掌握和推广使用。
发明内容
本发明的目的在于克服现有技术中的不足,提供一种针对工业多输入输出过程的分布式PI和PID控制器的定量整定方法。本发明控制器设计充分考虑了系统耦合的动态特性,整定过程简单直观,可以实现在线单参数化调节控制器,即对每个子控制器只需分别调节一个控制参数就能达到控制目的,因而操作简捷和方便,而且能够达到明显改进的控制效果。
本发明是通过以下技术方案实现的,本发明基于分布式闭环反馈控制结构,首先根据辨识得到的被控对象的传递函数矩阵模型,利用有效相对耦合分析方法建立有效等价对角传递函数矩阵模型,然后根据建立的有效等价对角传递函数矩阵模型提出期望的闭环传递函数矩阵,进而获得分布式PI(或PID)控制器的定量整定公式。实际运行控制系统时,在线单调地增减每个子控制器的单一调节参数,直至获得要求的控制系统标称性能及其鲁棒稳定性,并以最佳的方式在两者之间进行折中。
本发明在现有的电子控制设备和工控计算机中直接运行和实施,具体应用过程可以分为五个阶段:
第一、在组态界面上设置系统为“离线”调节状态,然后通过组态界面启动主机发出采样命令,当工控机的检测部分接到采样命令后,对被控制对象进行采样滤波,由模拟量输入通道将采样信号送入检测变送装置,再经A/D转换后得到数字信号后对对象进行辨识,对象辨识模块辨识出线性多变量过程的模型参数后,将模型参数送到主机的存储单元RAM中并由主机负责将数据显示在组态界面上。需要指出,这里的辨识模型是已经利用有效相对耦合阵列(ERGA)配对法则重新调整过的辨识模型。具体的有关有效相对耦合阵列(ERGA)配对法则的介绍和辨识方法可以参见相关技术文献,在这里不再详述并假设辨识,配对过程已经完成。工业多输入输出过程的传递函数矩阵辨识模型可用下式表示,
其中 g ij ( s ) = g oij ( s ) e - θ ij s , 它是指从被控过程的第i个输入到第j个输出的传递函数,goij(s)是其稳定正则的有理传递函数部分,θij是其相应的过程传输时滞。
第二、通过组态界面启动工控机的CPU调用事先编制好的“离线控制程序”解析的推导出分布式PI和PID控制器矩阵。详细的算法步骤如下:
①利用有效相对耦合分析方法建立有效对角传递函数矩阵,
G ~ ( s ) = diag { g ~ 0 ( s ) e - θ ~ jj s } = diag { g ojj e - θ ~ jj s } , ω → ∞ diag { 1 λ jj g oij e - θ ~ jj s } , ω → 0 - - - ( 2 )
其中
θ ~ jj = arg ( g ojj ( j ω ~ jj ) ) + π ω ~ jj - - - ( 3 )
ω ~ ij = λ ij ω ij ψ ij - - - ( 4 )
其中,λij是相对稳态增益矩阵 Λ = G ( 0 ) ⊗ G ( 0 ) - T 的元素,ψij是有效相对增益矩阵 Ψ = E ⊗ E - T 中的元素,这里
Figure C20061011727800098
E = G ( 0 ) ⊗ Ω - - - ( 6 )
其中
是相对关键频率矩阵,ωij是表达式arg[gij(jωij)]=-π的解。需要指出,在表示式(2)中,ω→∞表示高频段,ω→0表示低频段,但是具体作为区分高低频段的频率点并没有规定,这里采用该表达是方便PI/PID控制器参数的推导,具体步骤见第⑤步。从表达式(2)可以看出,建立有效对角传递函数矩阵包含了系统耦合引起的对角元素相位变化和稳态增益变化的所有信息,所以可以保证控制器的整定充分考虑系统耦合的动态特性,进而达到进一步改善系统控制效果的目的。
②根据建立的有效对角传递函数矩阵,利用鲁棒控制理论的H2最优性能目标,推导该分布式控制系统的期望闭环响应传递函数矩阵,
Figure C20061011727800101
其中
t jj ( s ) = 1 ( λ j s + 1 ) b Π k = 1 a ( z tjk - s z tjk + s ) e - θ ~ jj - - - ( 9 )
其中s=ztjk
Figure C20061011727800103
的右半平面零点,a是
Figure C20061011727800104
包含的右半平面零点的个数,b是
Figure C20061011727800105
的相对阶次,λj是用于得到实际期望的第j个控制闭环输出响应的调节参数,可以初始整定λj在(5-10)
Figure C20061011727800106
范围内。
③根据给出的期望闭环响应传递函数矩阵形式以及建立的有效对角传递函数矩阵,利用直接设计方法(direct design),推导分布式闭环控制器矩阵,
C ( s ) = diag { c jj ( s ) } = G ~ - 1 ( s ) T ( s ) I - T ( s )
= diag { 1 g ~ 0 ( s ) e - θ ~ jj s Π k = 1 a ( z tjk - s z tjk + s ) e - θ ~ jj ( λ j s + 1 ) b - Π k = 1 a ( z tjk - s z tjk + s ) e - θ ~ jj } - - - ( 10 )
④步骤③给出的分布式闭环控制器矩阵中子控制器的形式均以复杂的形式含有无理项,因此无法物理实现。这里采用数学Maclaurin展开级数对它们做有理逼近,从而方便实际构造分布式控制器矩阵。考虑到(10)中给出的分布式闭环控制器矩阵的子控制器的形式均具有积分特性,不能直接进行数学求导运算,所以先令
c jj ( s ) = c hj ( s ) = p hj ( s ) s , ω → ∞ c lj ( s ) = p lj ( s ) s , ω → 0 - - - ( 11 )
其中
c hj ( s ) = 1 g ojj e - θ ~ jj s Π k = 1 n ( z tjk - s z tjk + s ) e - θ ~ jj ( λ j s + 1 ) m - Π k = 1 n ( z tjk - s z tjk + s ) e - θ ~ jj - - - ( 12 )
c lj ( s ) = λ jj g ojj e - θ ~ jj s Π k = 1 n ( z tjk - s z tjk + s ) e - θ ~ jj ( λ j s + 1 ) m - Π k = 1 n ( z tjk - s z tjk + s ) e - θ ~ jj - - - ( 13 )
然后利用数学Maclaurin展开级数对其进行有理逼近,由此得到可以实际物理实现的控制器的有理逼近形式
c jj ( s ) = p hj ( s ) s = 1 s [ p hj ( 0 ) + p hj ′ ( 0 ) s + p hj ′ ′ ( 0 ) 2 ! s 2 + . . . ] , ω → ∞ ( a ) p lj ( s ) s = 1 s [ p lj ( 0 ) + p lj ′ ( 0 ) s + p lj ′ ′ ( 0 ) 2 ! s 2 + . . . ] , ω → 0 ( b ) - - - ( 14 )
其中
p hj ( n ) ( 0 ) = lim s → 0 d n ds [ s c hj ( s ) ] - - - ( 15 )
p lj ( n ) ( 0 ) = lim s → 0 d n ds [ sc lj ( s ) ] - - - ( 16 )
⑤工业生产中通常采用的分布式控制器矩阵中的子元素的具体形式为PI或PID,因此根据步骤④给出的可以实际物理实现的子控制器的有理逼近形式,取前两项就可以构造出子控制器的PI控制器的形式。同理,取其前三项就可以构成子控制器的PID控制器的形式。需要指出,因为PID控制算法中比例常数,微分时间常数在高频段起主导作用,而积分时间常数在低频段起主导作用,所以在通过给出的对分布式闭环控制器矩阵中子控制器的有理逼近获得PI或PID控制器参数时,比例常数和微分常数应采用表达式(14a),而积分常数应采用表达式(14b)。具体分布式控制器矩阵中的子PI或PID控制器的设计公式为,
c jj - PI ( s ) = K Cjj + 1 T Ijj s - - - ( 17 )
其中KCjj=p′hj(0)为比例增益,TIjj=plj(0)为积分时间常数。
C jj - PID ( s ) = ( K Cjj + T Ijj s + T Djj s ) 1 T Fjj s + 1 - - - ( 18 )
其中KCjj和TIjj如前所示, T Djj = p hj ′ ′ ( 0 ) 2 为微分时间常数,TFjj=(0.01-0.1)TDjj
第三、通过组态界面单调调节分布式控制器矩阵中的可调参数,并观察系统闭环响应曲线,由此离线确定最佳控制器参数。上述分布式控制器矩阵中的子PI(cjj-PI(s))和PID(Cjj-PID(s))控制器的设计公式(17-18)中的控制参数实质上均由单一调节参数λj整定,它对应着步骤②所示的期望控制闭环响应的形式和性能。调节参数λj的整定规则:调小λj可以加快对应的过程输出响应速度,提高控制系统的标称性能,但是相应所需的第j个控制器的输出能量要增大,并且它所对应的执行机构所需要提供的输出能量也要增大,会倾向于超出其容量范围,此外,在面临被控过程的未建模动态特性时,易于表现出过激行为,不利于控制系统的鲁棒稳定性;相反,增大λj会使对应的过程输出响应变缓,但是所要求的第j个控制器的输出能量减小,并且其所对应的执行机构所需要的输出能量也会减小,从而有利于提高控制系统的鲁棒稳定性。因此实际整定调节参数λj时,应在控制系统输出响应的标称性能与每个控制器及其执行机构的输出容量之间权衡。
第四、在组态界面上设置系统为“在线”调节状态,启动工控机的CPU读取最佳控制器参数,并执行“在线控制程序”得到当前时刻最优的控制量。“在线控制程序”可以按照提供的离散化公式编写,具体如下:
离散域PID控制算法公式为
Δ u jj ( n ) = K Cjj ( 1 + T s T Ijj + T Djj T s ) e jj ( n ) - K Cjj ( 1 + 2 T Djj T s ) e jj ( n - 1 ) + K Cjj T Djj T s e jj ( n - 2 ) - - - ( 19 )
ujj(n)=ujj(n-1)+Δujj(n)    (20)
式中KCjj-C(s)中第j个子控制器Cjj(s)的比例常数;
    TIjj-C(s)中第j个子控制器Cij(s)的积分时间常数,;
    TDjj-C(s)中第j个子控制器Cij(s)的微分时间常数;
    Ts-控制系统的采样周期;
    Δujj(n)-当前(n)时刻C(s)中第j个子控制器Cjj(s)的输出信号增量;
    ujj(n-1)-(n-1)时刻第j个子控制器Cjj(s)的控制器输出;
    ujj(n)-当前(n)时刻第j个子控制器Cjj(s)的控制器输出;
    ejj(n)-当前(n)时刻C(s)中第j个子控制器Cjj(s)输出与其给定值输入的偏差量;
    ejj(n-1)-前面第(n-1)时刻C(s)中第j个子控制器Cjj(s)的输出与其给定值输入的偏差量;
    ejj(n-2)-前面第(n-2)时刻C(s)中第j个子控制器Cjj(s)输出与其给定值输入的偏差量。
第五、对ujj(n)进行限幅,防止积分饱和,然后由D/A转换后输出至执行器,由执行器作用到被控对象,使被控对象运行在给定的范围内.此时组态界面上显示的是在线情况下的系统闭环响应曲线,然后用和离线调节方式下相同的方法进行在线微调。
全套调节过程可以在工控机组态界面上完成,此过程可以参看本说明书下文中提供的工业实例应用。与其它分布式控制器整定方法相比,本发明给出的分布式PI和PID控制器整定方法的最大特点是系统耦合的动态特性被充分考虑在控制器整定过程中,整定过程简单,直观,数据运算量小,并且可以实现系统性能的定量整定。在实际工业过程控制现场采用本发明给出的分布式PI和PID控制器的整定方法,只需在线单调地调节每个子控制器中的单一调节参数即可使控制系统达到指定的工作性能,而且便于实现控制系统的标称性能和鲁棒稳定性之间的最佳折中,因而操作简便,使得控制系统输出响应快速平稳,从而显著地克服了传统整定方法的主要缺点。本发明给出的分布式PI和PID控制器定量整定方法可广泛应用于能源、冶金、石化、轻工、医药、建材、纺织等行业中的多输入输出生产过程的控制和调节。
附图说明
图1为本发明在设计整定分布式PI和PID控制器矩阵时所基于的闭环控制结构图。其中C为控制器,G为被控对象,r和y分别为闭环系统的输入和输出,u为控制器输出,d为扰动信号。
图2为本发明实际运行过程的结构示意图。
图3为实际聚合反应堆对象在标称情况下对于两路单位阶跃给定值输入信号所得到的输出闭环响应。
其中,图3(a)示出了第一个过程输出的响应曲线,图3(b)示出了第二个过程输出的响应曲线。实线表示本发明,其中粗实线表示分布式PI,细实线表示分布式PID,点线表示Chien方法。
图4为实际聚合反应堆对象在有乘性不确定性情况下对于两路单位阶跃给定值输入信号所得到的输出闭环响应。
其中,图4(a)示出了第一个过程输出的响应曲线,图4(b)示出了第二个过程输出的响应曲线。实线表示系统在有乘性输入不确定性情况下,没有调节控制器参数的输出响应曲线,点线表示系统在有乘性输入不确定性情况下,调节控制器参数后的输出响应曲线。
具体实施方式
下面结合附图对本发明的实施例作详细说明:本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
对于一个广泛研究采用的工业聚合反应堆对象
G ( s ) = 22.89 e - 0.2 s 4.572 s + 1 - 11.64 e - 0.4 s 1.807 s + 1 4.689 e - 0.2 s 2.174 s + 1 5.80 e - 0.4 s 1.801 s + 1
本实施例给出的分布式PI和PID控制器的整定方法,第一步:按照附图1所示的结构框图组建一个分布式闭环控制系统,其中具体配对情况根据有效相对耦合配对法则来确定,本例中最佳配对为1-1/2-2;
第二步:在组态界面上点击“运行”键,启动工控机的CPU调用事先编制好的“离线控制程序”解析的推导出最优控制器。具体计算过程如下:
(1)应用公式(2-7),建立有效对角传递函数矩阵为
G ~ ( s ) = 22.89 e - 0.2014 s 4.572 s + 1 0 0 5.8 e - 0.4012 s 1.801 s + 1 , ω → ∞ 22.89 e - 0.2014 s 0.7087 ( 4.572 s + 1 ) 0 0 5.8 e - 0.4012 s 0.7087 ( 1.801 s + 1 ) , ω → 0
(2)应用公式(8-9)求解出期望闭环传函,
T ( s ) = e - 0.2014 s λ 1 s + 1 0 0 e - 0.4012 s λ 2 s + 1
(3)应用公式(10)得出理想期望的分布式闭环控制器的形式
C ( s ) = G ~ - 1 ( s ) e - 0.2014 s λ 1 s + 1 - e - 0.2014 s 0 0 e - 0.4012 s λ 2 s + 1 - e - 0.4012 s
(4)应用公式(11-16)得出理想期望的两个分布式闭环控制器有理逼近实现的级数形式;
(5)初始设置两个调节参数λ1=0.3,λ2=0.87,应用公式(17-18)得出分布式PI和PID控制器的形式,即
c 11 - PI ( s ) = 0.4019 + 0.2529 s ; c 22 - PI ( s ) = 0.0618 + 0.0961 s
c 11 - PID ( s ) = 0.4019 + 0.2529 s + 0.0114 s
c 22 - PID ( s ) = 0.0618 + 0.0961 s - 0.0105 s
需要说明,如上初始设置分布式PI和PID控制器的调节参数λ1和λ2的值是为了得到与Chien的方法下的系统给定值响应的上升时间相同,以便比较。
第三步:通过组态界面观察系统闭环响应。输入信号按以下方式加入:先在t=0秒时刻给第一路输入量加入单位阶跃输入信号:r1=1/s,而第二路输入信号为r2=0,然后在t=25秒时刻给第二路输入也都加入单位阶输入信号。所得到的系统闭环响应曲线为图3中所示。从图3中可以看出,本实施例给出的分布式PI和PID控制器矩阵调节的系统输出响应(实线)明显优于Chien的方法(点线)。然后假设实际被控过程G(s)存在乘性输入不确定性ΔI=diag{(s+0.3)/(s+1),(s+0.3)/(s+1)}。这里ΔI可以近似地物理解释为,被控过程的两个输入调节阀在高频段具有高达100%的不确定性,并且在低频段工作范围具有将近30%的不确定性。在这种严重的过程输入不确定性下进行如上所述仿真实验,本发明给出的控制器的整定方法所得到的过程输出响应的计算机仿真结果如附图4所示。由图4可以看到,本实施例给出的控制器的整定方法(实线)能够良好地保证系统的鲁棒稳定性。此外可以看到,单调地增大控制器中的调节参数λj就可以使过程输出的给定值响应的振荡减小,如图4中的点线所示。因此,采用本实施例给出的控制器的整定方法可以很方便地在线进行单调地调节系统输出响应,从而达到实际要求的工作指标。
第四步:在组态界面上设置系统为“在线”调节状态,启动工控机的CPU读取最佳控制器参数,并执行“在线控制程序”得到当前时刻最优的控制量。
第五步:对u(k)进行限幅,防止积分饱和,然后由D/A转换后输出至执行器,由执行器作用到被控对象,使被控对象运行在给定的范围内.此时组态界面上显示的是在线情况下的系统闭环响应曲线,观察曲线进行在线微调,如此周而复始实现控制。
以上阐述的是本发明给出的一个实施例所表现出的优良控制效果。本发明针对工业过程中的一般多输入输出过程模型给出分布式PI和PID控制器的定量整定方法,所以适用于各种不同的工业生产过程。本发明给出的分布式PI和PID控制器的定量整定方法可广泛应用于石化、冶金、医药、建材和纺织等行业的生产过程。

Claims (2)

1、一种工业多输入输出过程的分布式PI和PID控制器定量整定方法,其特征在于,包括如下具体步骤:
第一、采用已经利用有效相对耦合阵列配对法则重新调整过的辨识模型,工业多输入输出过程的传递函数矩阵辨识模型用下式表示,
Figure C2006101172780002C1
其中 g ij ( s ) = g oij ( s ) e - θ ij s , 它是指从被控过程的第i个输入到第j个输出的传递函数,goij(s)是其稳定正则的有理传递函数部分,θij是其相应的过程传输时滞;
第二、通过组态界面启动工控机的CPU调用事先编制好的“离线控制程序”解析的推导出分布式PI和PID控制器矩阵,步骤如下:
①利用有效相对耦合分析方法建立有效对角传递函数矩阵,
G ~ ( s ) = diag { g ~ o ( s ) e - θ ~ jj s } = diag { g ojj e - θ ~ jj s } , ω → ∞ diag { 1 λ jj g oij e - θ ~ jj s } , ω → ∞ - - - ( a )
其中
θ ~ jj = arg ( g ojj ( j ω ~ jj ) ) + π ω ~ jj
ω ~ ij = λ ij ω ij ψ ij
其中,λij是相对稳态增益矩阵 Λ = G ( 0 ) ⊗ G ( 0 ) - T 的元素,ψij是有效相对增益矩阵 ψ = E ⊗ E - T 中的元素,这里
Figure C2006101172780002C8
E = G ( 0 ) ⊗ Ω
其中
是相对关键频率矩阵,ωij是表达式arg[gij(jωij)]=-π的解,在表示式(a)中,ω→∞表示高频段,ω→0表示低频段,但是具体作为区分高低频段的频率点并没有规定,这里采用该表达是方便PI/PID控制器参数的推导,具体步骤见第⑤步;
②根据建立的有效对角传递函数矩阵,利用鲁棒控制理论的H2最优性能目标,推导分布式控制系统的期望闭环响应传递函数矩阵,
其中
t jj ( s ) = 1 ( λ j s + 1 ) b Π k = 1 a ( z tjk - s z tjk + s ) e - θ ~ jj
其中s=ztjk
Figure C2006101172780003C5
的右半平面零点,a是
Figure C2006101172780003C6
包含的右半平面零点的个数,b是
Figure C2006101172780003C7
的相对阶次,λj是用于得到实际期望的第j个控制闭环输出响应的调节参数,初始整定λj
Figure C2006101172780003C8
范围内;
③根据给出的期望闭环响应传递函数矩阵形式以及建立的有效对角传递函数矩阵,利用直接设计方法,推导分布式闭环控制器矩阵,
C ( s ) = diag { c jj ( s ) } = G ~ - 1 ( s ) T ( s ) I - T ( s )
= diag { 1 g ~ o ( s ) e - θ ~ jj s Π k = 1 a ( z tjk - s z tjk + s ) e - θ ~ jj ( λ j s + 1 ) b - Π k = 1 a ( z tjk - s z tjk + s ) e - θ ~ jj }
④采用数学Maclaurin展开级数对步骤③给出的分布式闭环控制器矩阵中子控制器的形式做有理逼近,由此得到可以实际物理实现的控制器的有理逼近形式
c jj ( s ) = p hj ( s ) s = 1 s [ p hj ( 0 ) + p hj ′ ( 0 ) s + p hj ′ ′ ( 0 ) 2 ! s 2 + . . . ] , ω → ∞ - - - ( b ) p lj ( s ) s = 1 s [ p lj ( 0 ) + p lj ′ ( 0 ) s + p lj ′ ′ ( 0 ) 2 ! s 2 + . . . ] , ω → 0 - - - ( c )
其中
p hj ( n ) ( 0 ) = lim s → 0 d n ds [ sc hj ( s ) ]
p lj ( n ) ( 0 ) = lim s → 0 d n ds [ sc lj ( s ) ]
这里
c hj ( s ) = 1 g ojj e - θ ~ jj s Π k = 1 n ( z tjk - s z tjk + s ) e - θ ~ jj ( λ j s + 1 ) m - Π k = 1 n ( z tjk - s z tjk + s ) e - θ ~ jj
c lj ( s ) = λ jj g ojj e - θ ~ jj s Π k = 1 n ( z tjk - s z tjk + s ) e - θ ~ jj ( λ j s + 1 ) m - Π k = 1 n ( z tjk - s z tjk + s ) e - θ ~ jj
⑤根据步骤④给出的能实际物理实现的子控制器的有理逼近形式,取前两项构造出子控制器的PI控制器的形式,同理,取其前三项构成子控制器的PID控制器的形式,具体分布式控制器矩阵中的子PI或PID控制器的设计公式为,
c jj - PI ( s ) = K Cjj + 1 T Ijj s - - - ( d )
其中KCjj=p′hj(0)为比例增益,TIjj=plj(0)为积分时间常数;
C jj - PID ( s ) = ( K Cjj + T Ijj s + T Djj s ) 1 T Fjj s + 1 - - - ( e )
其中KCjj和TIjj如前所示, T Djj = p hj ′ ′ ( 0 ) 2 为微分时间常数,TFjj=(0.01-0.1)TDjj
第三、通过组态界面单调调节分布式控制器矩阵中的可调参数,并观察系统闭环响应曲线,由此离线确定最佳控制器参数,上述分布式控制器矩阵中的子PI(cjj-PI(s))和PID(Cjj-PID(s))控制器的设计公式(d、e)中的控制参数实质上均由单一调节参数λj整定,它对应着步骤②所示的期望控制闭环响应的形式和性能,实际整定调节参数λj时,应在控制系统输出响应的标称性能与每个控制器及其执行机构的输出容量之间权衡;
第四、在组态界面上设置系统为“在线”调节状态,启动工控机的CPU读取最佳控制器参数,并执行“在线控制程序”得到当前时刻最优的控制量;
第五、对ujj(n)进行限幅,防止积分饱和,然后由D/A转换后输出至执行器,由执行器作用到被控对象,使被控对象运行在给定的范围内,此时组态界面上显示的是在线情况下的系统闭环响应曲线,然后用和离线调节方式下相同的方法进行在线微调。
2、如权利要求1所述的工业多输入输出过程的分布式PI和PID控制器定量整定方法,其特征是,所述的“在线控制程序”,按照以下提供的离散化公式编写:
离散域PID控制算法公式为
Δ u jj ( n ) = K Cjj ( 1 + T s T Ijj + T Djj T s ) e jj ( n ) - K Cjj ( 1 + 2 T Djj T s ) e jj ( n - 1 ) + K Cjj T Djj T s e jj ( n - 2 )
ujj(n)=ujj(n-1)+Δujj(n)
以上两个控制信号增量算式中:
KCjj-C(s)中第j个子控制器Cjj(s)的比例常数;
TIjj-C(s)中第j个子控制器Cjj(s)的积分时间常数;
TDjj-C(s)中第j个子控制器Cjj(s)的微分时间常数;
Ts-控制系统的采样周期;
Δujj(n)-当前(n)时刻C(s)中第j个子控制器Cjj(s)的输出信号增量;
ujj(n-1)-(n-1)时刻第j个子控制器Cjj(s)的控制器输出;
ujj(n)-当前(n)时刻第j个子控制器Cjj(s)的控制器输出;
ejj(n)-当前(n)时刻C(s)中第j个子控制器Cjj(s)输出与其给定值输入的偏差量;
ejj(n-1)-前面第(n-1)时刻C(s)中第j个子控制器Cjj(s)的输出与其给定值输入的偏差量;
ejj(n-2)-前面第(n-2)时刻C(s)中第j个子控制器Cjj(s)输出与其给定值输入的偏差量。
CNB2006101172782A 2006-10-19 2006-10-19 工业多输入输出过程的分布式pi和pid控制器定量整定方法 Expired - Fee Related CN100428090C (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2006101172782A CN100428090C (zh) 2006-10-19 2006-10-19 工业多输入输出过程的分布式pi和pid控制器定量整定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2006101172782A CN100428090C (zh) 2006-10-19 2006-10-19 工业多输入输出过程的分布式pi和pid控制器定量整定方法

Publications (2)

Publication Number Publication Date
CN1936737A CN1936737A (zh) 2007-03-28
CN100428090C true CN100428090C (zh) 2008-10-22

Family

ID=37954287

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2006101172782A Expired - Fee Related CN100428090C (zh) 2006-10-19 2006-10-19 工业多输入输出过程的分布式pi和pid控制器定量整定方法

Country Status (1)

Country Link
CN (1) CN100428090C (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102007032484A1 (de) * 2007-07-12 2009-01-22 Beckhoff Automation Gmbh Regelverfahren und Regelvorrichtung mit mehrkanaliger Rückführung
CN102426417B (zh) * 2011-12-13 2013-10-02 中冶南方(武汉)自动化有限公司 Pi参数混合整定法
CN103576711B (zh) * 2013-11-14 2015-09-09 上海交通大学 基于定量单参数pid控制的化工反应器温度控制方法
CN105988437A (zh) * 2015-02-10 2016-10-05 张家港市宏鑫源科技有限公司 一种采用分布式控制系统的多级万能轧机
CN104793496B (zh) * 2015-04-16 2017-06-06 北京化工大学 二输入二输出聚合反应器的解耦控制系统
CN104950725A (zh) * 2015-07-31 2015-09-30 北京化工大学 一种二元精馏塔分布式控制系统
CN105337298A (zh) * 2015-11-10 2016-02-17 南方电网科学研究院有限责任公司 交直流系统附加控制器通道间相互影响的评估方法和系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6445980B1 (en) * 1999-07-10 2002-09-03 Mykrolis Corporation System and method for a variable gain proportional-integral (PI) controller
CN1445629A (zh) * 2003-03-06 2003-10-01 上海交通大学 定量整定鲁棒性的智能比例积分微分控制方法
CN1588245A (zh) * 2004-07-22 2005-03-02 上海交通大学 化工多变量生产过程的解耦控制系统
CN1645276A (zh) * 2005-02-03 2005-07-27 上海交通大学 化工双输入输出过程的分布式pi和pid控制器的定量整定方法
JP2005258717A (ja) * 2004-03-10 2005-09-22 Toshiba Corp 制御装置のパラメータ設定方法及び装置並びにプログラム
CN1794118A (zh) * 2005-12-29 2006-06-28 上海交通大学 单输入单输出系统的极限pid控制方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6445980B1 (en) * 1999-07-10 2002-09-03 Mykrolis Corporation System and method for a variable gain proportional-integral (PI) controller
CN1445629A (zh) * 2003-03-06 2003-10-01 上海交通大学 定量整定鲁棒性的智能比例积分微分控制方法
JP2005258717A (ja) * 2004-03-10 2005-09-22 Toshiba Corp 制御装置のパラメータ設定方法及び装置並びにプログラム
CN1588245A (zh) * 2004-07-22 2005-03-02 上海交通大学 化工多变量生产过程的解耦控制系统
CN1645276A (zh) * 2005-02-03 2005-07-27 上海交通大学 化工双输入输出过程的分布式pi和pid控制器的定量整定方法
CN1794118A (zh) * 2005-12-29 2006-06-28 上海交通大学 单输入单输出系统的极限pid控制方法

Non-Patent Citations (14)

* Cited by examiner, † Cited by third party
Title
A Simple Multiloop Tuning Method for PID Controllerswith No Proportional Kick. I-Lung Chien,Hsiao-Ping Huang and Jen-Chien Yang.INDUSTRIAL & ENGINEERING CHEMISTRY RESEARCH,Vol.38 No.4. 1999
A Simple Multiloop Tuning Method for PID Controllerswith No Proportional Kick. I-Lung Chien,Hsiao-Ping Huang and Jen-Chien Yang.INDUSTRIAL & ENGINEERING CHEMISTRY RESEARCH,Vol.38 No.4. 1999 *
A Simple Multiloop Tuning Method for PID Controllerswith No Proportional Kick. I-Lung Chien,Hsiao-Ping Huang and Jen-Chien Yang.INDUSTRIAL &amp *
Decentralized ControlSystem Design for MultivariableProcessess -A Novel MethodBased on Effective RelativeGain Array. Qiang Xiong,WenJian Cai,MaoJun He,Ming He.INDUSTRIAL & ENGINEERING CHEMISTRY RESEARCH,Vol.45 No.8. 2006
Decentralized ControlSystem Design for MultivariableProcessess -A Novel MethodBased on Effective RelativeGain Array. Qiang Xiong,WenJian Cai,MaoJun He,Ming He.INDUSTRIAL & ENGINEERING CHEMISTRY RESEARCH,Vol.45 No.8. 2006 *
Decentralized ControlSystem Design for MultivariableProcessess -A Novel MethodBased on Effective RelativeGain Array. Qiang Xiong,WenJian Cai,MaoJun He,Ming He.INDUSTRIAL &amp *
ENGINEERING CHEMISTRY RESEARCH,Vol.38 No.4. 1999 *
ENGINEERING CHEMISTRY RESEARCH,Vol.45 No.8. 2006 *
典型不稳定时滞对象的H ∞最优灵敏度控制. 张卫东,许晓鸣.上海交通大学学报,第35卷第9期. 2001
典型不稳定时滞对象的H ∞最优灵敏度控制. 张卫东,许晓鸣.上海交通大学学报,第35卷第9期. 2001 *
多变量时滞过程的解耦控制设计. 刘涛,张卫东,顾诞英.自动化学报,第31卷第6期. 2005
多变量时滞过程的解耦控制设计. 刘涛,张卫东,顾诞英.自动化学报,第31卷第6期. 2005 *
给定值和干扰响应解耦的新型控制器设计. 张卫东,许晓鸣,何星.自动化学报,第27卷第1期. 2001
给定值和干扰响应解耦的新型控制器设计. 张卫东,许晓鸣,何星.自动化学报,第27卷第1期. 2001 *

Also Published As

Publication number Publication date
CN1936737A (zh) 2007-03-28

Similar Documents

Publication Publication Date Title
CN100428090C (zh) 工业多输入输出过程的分布式pi和pid控制器定量整定方法
CN1332274C (zh) 化工双输入输出过程的分布式pi和pid控制器的定量整定方法
CN100462877C (zh) 工业过程中非方系统的解耦控制方法
Ramasamy et al. Control of ball mill grinding circuit using model predictive control scheme
CN100458602C (zh) 多输入多输出系统的极限pid控制方法
CN100337169C (zh) 化工多变量生产过程的两自由度解耦控制系统
Zhao et al. Performance-based parameter tuning method of model-driven PID control systems
Zhang et al. Dominant pole placement for multi-loop control systems
CN101950156A (zh) 一种自适应串级pid控制方法
Silva et al. Controller design via Padé approximation can lead to instability
Isa et al. Comparative study of PID controlled modes on automatic water level measurement system
CN1765373A (zh) 中药生产浓缩过程中蒸发速度的控制方法
CN100422883C (zh) 线性多变量工业过程的全局最优控制器整定方法
Evald et al. A hybrid robust model reference adaptive controller and proportional integral controller without reference model for partially modeled systems
Helian et al. Adaptive robust control of a pump control hydraulic system
Poulsen et al. Constrained predictive control and its application to a coupled-tanks apparatus
Dawane et al. Original research Ar-ticle Performance improvement of DC servo motor using sliding mode controller
Tian et al. Double-controller scheme for separating load rejection from set-point tracking
Xiong et al. Adaptive control of a synchronizing servo-system
CN100449433C (zh) 分布式内模控制系统的参数定量整定方法
Nagaraj et al. Optimum tuning algorithms for PID controller—a soft computing approach
Vogel et al. Industrial experience with state-space model predictive control
Huang et al. Feedback linearization and equivalent-disturbance compensation control strategy for piezoelectric stage
CN108319141A (zh) 结合线性二次调节器与分数阶滑模控制的复合控制方法
CN115780006B (zh) 一种基于案例推理的高压辊磨机料重智能控制方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20081022

Termination date: 20111019