多变量时滞系统解耦控制器解析设计方法

技术领域

本发明涉及的是一种工业过程控制技术领域的方法，具体是一种多变量时滞系统解耦控制器解析设计方法。

背景技术

化工生产中普遍存在具有传输时滞的多输入多输出过程，而且随着各种各样的先进生产工艺的快速发展，越来越多的生产过程被构造为高维多变量控制系统，从而实现高效率地生产高质量的产品。然而由于多变量过程的各输出通道之间存在交联耦合作用，使得大多数已发展的单变量控制方法很难用于多变量过程。目前针对多变量系统的研究方法主要包括线性二次最优(LQG)设计方法，定量反馈方法，奈奎斯特阵列方法，序列回差方法，并矢展开法，预测控制方法，多变量内模设计方法等。其中针对多变量系统的内模设计方法要求对控制对象进行内外分解。但是迄今为止，内外分解方法都是基于状态空间方法并使用数值运算，分解方法只能应用于有理对象，对含有时滞的多变量系统，不稳定系统都无法使用。其它的方法则不能很好的解决多变量系统鲁棒性设计问题。最重要的是上述的设计方法只针对不包含有时滞的多变量系统，对包含时滞的多变量对象则无能为力。目前，具有代表性的针对时滞多变量系统的解耦控制方法主要包括静态解耦器、动态解耦器、顺序解耦以及解耦控制器矩阵的设计方法，但是这些方法得到的控制器大都不便于在线调节和设定，并且所用到的相关专业理论知识较多，不便于被工程技术人员掌握和推广使用。

经对现有技术的文献检索发现，Dong等人在《Proceedings of the AmericanControl Conference》(美国控制会议)(1997年第5卷第3380-3384页)上发表了“Design of Robust Multivariable PID Controllers via IMC”(基于内模控制结构设计鲁棒多变量PID控制器)，文中提出一种解耦控制器设计方法，该方法是把控制对象模型和理论上基于对象模型分解得到的最优解耦控制器的一般形式都使用近似方法进行近似，然后通过对应元素相等来得到解耦控制器的具体表达式。其不足在于没有给出最优的控制对象分解步骤，这使该方法的实际应用受到了很大的限制。

发明内容

本发明的目的在于克服现有技术中的不足，提供一种针对包含时滞的多变量时滞系统解耦控制器解析设计方法。使其采用输入输出设计方法，设计方法容易理解，使用方便，不需要状态信息。本发明设计的控制器可以实现系统输出解耦，并且只需要调节有限的几个参数，就可以方便有效的定量调节系统性能和鲁棒性。

本发明是通过以下技术方案实现的，本发明具体步骤如下：

1)首先对化工多变量时滞过程的传递函数矩阵辨识模型

其中 $>>>g>ij>>>(>s>)>>=>>g>oij>>>(>s>)>>>e>>->>\theta >ij>>s>>>,>>>它是指从被控过程的第i个输入到第j个输出的传递函数，g$_oij(s)是其稳定正则的有理传递函数部分，θ_ij是其相应的过程传输时滞，i，j＝1，…p，进行分解，分解形式如下：

G(s)＝G_D(s)G_O(s)，                                            (1)

其中， $>>>G>D>>>(>s>)>>=>diag>\{>>e>>->>\theta >i>>s>>>,>.>.>.>,>>e>>->>\theta >p>>s>>>\}>,>->->->>(>2>)>>>>$

$>>>G>O>>>(>s>)>>=sup>>G>D>>->1>sup>>>(>s>)>>G>>(>s>)>>,>->->->>(>3>)>>>>$

其中θ_i取为G^-1(s)的第i列中最大预估值。

2)设计用于保证控制器可实现的两个调节因子：

$>>>G>N>>>(>s>)>>=>diag>\{>>\Pi >>j>=>1>>>r>z>>>>>(>>>->s>+>>z>j>>>>s>+sup>>z>j>*sup>>>>)>>>k>ij>>>,>.>.>.>,>>\Pi >>j>=>1>>>r>z>>>>>(>>>->s>+>>z>j>>>>s>+sup>>z>j>*sup>>>>)>>>k>pj>>>\}>->->->>(>4>)>>>>$

$>>J>>(>s>)>>=>diag>\{>>1>>>(>>\lambda >i>>s>+>1>)>>>p>i>>>>,>.>.>.>,>>1>>>(>>\lambda >n>>s>+>1>)>>>p>n>>>>\}>->->->>(>5>)>>>>$

其中z_j是G^-1(s)第i列中不稳定的极点，k_ij是z_j的最大个数。λ_i为控制器可调参数，p_i取G(s)中第i列元素的最大相对阶次。

3)设计解耦的控制器C为如下形式：

$>>C>>(>s>)>>=>>>>G>>->1>>>>(>s>)>>>G>D>>>(>s>)>>>G>N>>>(>s>)>>J>>(>s>)>>>>I>->>G>D>>>(>s>)>>>G>N>>>(>s>)>>J>>(>s>)>>>>.>->->->>(>6>)>>>>$

需要说明，在上述设计过程中，如果对象传递函数矩阵分解后G_O(s)中仍然含有时滞项，这将导致设计的控制器是无穷高阶，这时使用控制器降阶技术对该控制器进行降阶。

本发明在现有工控计算机上可以直接运行实施，具体实施过程如下：

a)在组态界面上设计可滑动的滑块，编写程序使滑块的移动与解耦控制器中可调参数的变化保持一致。整定时先确定部分参数的初始值，主要包括根据实际工况对控制对象的模型进行辨识得到对象模型参数；根据系统动态性能和鲁棒稳定性的要求确定期望的性能指标参数，控制器可调参数初始值，通常可取1。然后由工控机将数据送到存储单元RAM中；并在组态界面上设置系统为“离线”调节状态。其中多变量对象在线辨识方法有很多种，如继电反馈法等。

b)在组态界面上点击“运行”键，由此启动工控机的CPU执行事先编制好的实现解耦控制器的程序。此程序算法正是在上述新提出的因式分解规则的基础上创新得出的。

c)单调调节组态界面上的滑块来调节解耦控制器的可调参数，并观察系统闭环响应曲线，由此确定最佳控制器参数。调节参数λ_i的整定规则：调小λ_i可以加快对应的过程输出响应速度，提高控制系统的标称性能，但是相应所需的第i个控制器的输出能量要增大，并且它所对应的执行机构所需要提供的输出能量也要增大，会倾向于超出其容量范围，此外，在面临被控过程的未建模动态特性时，易于表现出过激行为，不利于控制系统的鲁棒稳定性；相反，增大λ_i会使对应的过程输出响应变缓，但是所要求的第i个控制器的输出能量减小，并且其所对应的执行机构所需要的输出能量也会减小，从而有利于提高控制系统的鲁棒稳定性。因此实际整定调节参数λ_i时，应在控制系统输出响应的标称性能与每个控制器及其执行机构的输出容量之间权衡。

d)在组态界面上设置系统为“在线”调节状态，由此启动工控机的CPU读取最佳控制器参数，然后按照离散化公式计算离散域解耦控制算式，得到当前时刻离线状态下最优的控制量。

e)对u_j(k)进行限幅，防止积分饱和，然后由D/A转换后输出至执行器，由执行器作用到被控对象，使被控对象运行在给定的范围内。此时组态界面上显示的是在线情况下的系统闭环响应曲线，然后用和离线调节方式下相同的方法进行在线微调，如此周而复始实现控制。

所述的d)中，按照离散化公式计算离散域解耦控制算式，得到当前时刻离线状态下最优的控制量，具体如下：

首先对控制器中每个元素进行离散化，并化成形如式(7)的标准形式：

$>>>C>ij>>>(>z>)>>=>>>>b>1>>+>>b>2>>>z>>->1>>>+>·>·>·>>b>>m>->1>>>>z>>->>(>\beta >->1>)>>>>>>>a>1>>+>>a>2>>>z>>->1>>>+>·>·>·>>a>>n>->1>>>>z>>->>(>\alpha >->1>)>>>>>>=>>>>u>ij>>>(>z>)>>>>>e>ij>>>(>z>)>>>>->->->>(>7>)>>>>$

由表达式(7)得控制器输出控制量的分量表达式：

a₁u_ij(z)+a₂u_ij(z-1)+…a_n-1u_ij(z-α+1)

＝b₁e_ij(z)+b₂e_ij(z-1)+…b_m-1e_ij(z-β+1)

                                                        (8)

将上式写成时间递推形式为如下形式：

a₁u_ij(k)+a₂u_ij(k-1)+…a_n-1u_ij(k-α+1)

＝b₁e_ij(k)+b₂e_ij(k-1)+…b_m-1e_ij(k-β+1)

                                                        (9)

由此得到控制器输出控制量的表达式如下式(10)所示：

$>>>u>j>>>(>k>)>>=>>\Sigma >>i>=>1>>n>>>u>ij>>>(>k>)>>->->->>(>10>)>>>>$

上述表达式(8)中α，β分别表示控制器每一个元素分子分母的阶次。表达式(8-10)式中：

u_ij(k)-当前(k)时刻第C_ij(z)控制器的输出控制量

e_ij(k)-当前(k)时刻第C_ij(z)控制器的输入偏差量

u_ij(k-α+1)-当前(k-α+1)时刻第C_ij(z)控制器的输出控制量

e_ij(k-β+1)-当前(k-β+1)时刻第C_ij(z)控制器的输入偏差量

u_j(k)-当前(k)时刻控制器第j个输出控制量。

本发明全套调节过程可以在工控机组态界面上完成，与传统的多变量设计方法相比，本发明给出的针对多变量时滞系统工程模型的控制器设计方法的最大特点是具有很强的实用性。设计的控制器可以实现系统闭环响应解耦，控制器设计具有最优，解析的特点。用传统解耦控制方法进行控制时，很难分析系统的鲁棒性，使用本发明则可以实现系统性能和鲁棒性的定量调节。本发明除了可以用于普通的含有时滞的多变量系统，还可以用于含有时滞的不稳定非最小相位系统以及含有状态时滞的多变量过程。

已经证明，本发明给出的设计方法是最优的，本发明与Morari and Zafiriou，1989中内模控制方法的区别在于，内模控制是基于状态空间理论的针对稳定有理对象的非解耦设计方法。

附图说明

图1为本发明给出的解耦控制器的设计方案示意图。

图2为本发明实际运行时采用的结构示意图。

图3为本发明给出的解耦控制器设计方法所采用的闭环控制结构图。

图4为本发明给出的组态界面。

其中，图4(a)示出了初始状态下的组态界面，图4(b)示出了工作状态下的组态界面。

图5为精馏塔对象的输出闭环响应。

其中，图5(a)示出了在第一个阶跃输入作用下系统输出的响应曲线，图5(b)示出了第二个阶跃输入作用下系统输出的响应曲线，由本发明给出未经过近似的控制器作用下的系统输出响应曲线用点线表示，由本发明给出并经过近似的控制器作用下的系统输出响应曲线用实线表示。

图6为在有乘性不确定性作用下，精馏塔对象的输出闭环响应。

其中，图6(a)示出了在第一个阶跃输入作用下系统输出的响应曲线，图6(b)示出了第二个阶跃输入作用下系统输出的响应曲线，点线表示系统在有乘性不确定性情况下，没有调节控制器参数的输出响应曲线，实线表示系统在有乘性不确定性情况下，调节控制器参数后的输出响应曲线。