范 蓉¹，王为群¹，姚娟²

¹南京理工大学理学院，江苏 南京

²南京理工大学自动化学院，江苏 南京

收稿日期：2018年6月6日；录用日期：2018年7月3日；发布日期：2018年7月10日

本文针对具有时变时滞项的多面体不确定离散时间系统，研究了鲁棒预见跟踪控制问题。首先对时变时滞项作两项近似，并将近似误差看作外部扰动，则将原系统转换成受干扰抑制的常时滞不确定离散系统。其次，基于误差扰动的方法，通过引入一个差分算子导出一个包含可预见信息和外部扰动的增广误差系统。然后利用LMI技术，给出原系统在控制器下，闭环系统渐近稳定的充分条件并通过解LMI得到具有预见跟踪性能的输出反馈控制器。最后，通过数值算例验证了本文结论的有效性。

关键词 :预见跟踪控制，增广误差系统，时变时滞，多面体不确定系统，线性矩阵不等式

Copyright © 2018 by authors and beplay安卓登录

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

预见控制是一种通过充分利用系统已知的未来信息来设计控制器使得系统稳定的控制方法。它能够设计出具有信息补偿功能的控制器，使得系统的静态误差有效减小，系统的跟踪水平和响应速度显著提高 [ 1 ] 。预见控制的概念首次于1966年提出 [ 2 ] 。利用可预见信息，来提高闭环系统的跟踪性能是预见控制理论的主要研究内容，被应用于各种控制系统中：对于离散系统 [ 3 ] [ 4 ] ，引入差分算子构造增广误差系统，并分别通过LQI (线性二次积分)和LMI (线性矩阵不等式)技术得到状态反馈预见控制器。类似的，对于连续时间系统 [ 5 ] ，引入微分算子来构造增广误差系统，基于最优控制理论，导出了最优预见控制器。进一步，预见控制理论被推广到广义线性系统中 [ 6 ] [ 7 ] 。对于随机系统 [ 8 ] [ 9 ] ，构造一个包含引入的积分器的增广误差系统，将跟踪问题转换成一个调节器问题，再利用随机控制理论中的动态规划的方法解决预见跟踪控制问题。除此之外， [ 10 ] [ 11 ] 将预见控制理论与滑模控制系统结合。 [ 12 ] [ 13 ] 利用可预见信息考虑了模糊控制。 [ 14 ] 提出了具有 H 2 性能指标的预见控制器的设计方案。 [ 15 ] 基于博弈论，给出了具有 H ∞ 性能指标的预见控制器的设计方法。预见跟踪控制在实际中的应用有：汽车驾驶，机器人路径控制，优化无人机航迹等 [ 16 ] [ 17 ] [ 18 ] 。

在实际工程领域中，不确定因素是难以避免的。参数依赖的不确定系统(即凸多面体系统)是不确定系统的一种。文献 [ 19 ] [ 20 ] 描述了凸多面不确定系统的预见控制问题并通过对系统状态之间取差分的方式来构造增广误差系统，但是这样的构造方式不能直接用于时变时滞不确定系统，本文在这些工作的基础上考虑时变时滞不确定系统的预见跟踪控制问题。因为在实际工程中，考虑到状态变量并不是都是可测的，所以本文考虑给系统设计一个包含可预见信息的输出反馈控制器，使得系统输出在设计的控制器下能够准确跟踪到目标轨迹。

本文研究了参数依赖的不确定离散时变时滞系统，提出闭环系统渐近稳定的充分条件，并基于LMI工具箱给出预见输出反馈控制器的设计方案。通过数值仿真验真本文结果的有效性

考虑如下的时变时滞不确定离散系统：

{ x ( k + 1 ) = A ( θ ) x ( k ) + A d ( θ ) x ( k − d ( k ) ) + B ( θ ) u ( k ) y ( k ) = C ( θ ) x ( k ) (1)

其中 x ( k ) ∈ R p 是系统的状态向量， u ( k ) ∈ R m 系统的输入控制， y ( k ) ∈ R q 是系统的输出。 A ( θ ) ， A d ( θ ) ， B ( θ ) ， C ( θ ) 是恰当维数的不确定系数矩阵，它们不依赖于时间变量k。 d ( k ) 表示时变时滞，且满足： d 1 ≤ d ( k ) ≤ d 2 ，其中 d 1 和 d 2 是 d ( k ) 的上下界。

假设1：系统(1)的不确定系数矩阵满足如下关系式：

[ A ( θ )             A d ( θ )             B ( θ )           C ( θ ) ] = ∑ i = 1 s θ i [ A i         A d i         B i           C i ] (2)

其中 A i ， A d i ， B i ， C i ， ( i = 1 , 2 , ⋯ , s ) 为适当维数的常矩阵。 θ = ( θ 1     θ 2       ⋯         θ s ) T ∈ R s 为不确定参数向量，满足：

θ ∈ Θ : { θ ∈ R s | θ i ≥ 0 , ( i = 1 , 2 , ⋯ , s ) , ∑ i = 1 s θ i = 1 } (3)

假设2：参考信号 r ( k ) 的预见长度为 M r 即在时刻k， M r 步未来值 r ( k + 1 ) , r ( k + 2 ) , ⋯ , r ( k + M r ) 以及参考信号的现在和过去的值均可获得。并假定超过 M r 步的未来值是不变的，即

j = M r + 1 ,         M r + 2 ,       M r + 3 , ⋯ (4)

注1：假设1表明系统 r ( k + j ) = r ( k + M r ) (1)的不确定系数矩阵和时间变量k无关，且参数向量 θ 在s维凸多面体中取值。假设2是预见控制问题中的标准假设。理论研究和实际例子验证可预见信号仅在某一段步长内对系统的性能产生影响，超过预见长度的参考信号对系统性能的影响是比较弱的，因此本文中，假设超过预见长度的参考信号是个定常数。则显然有：

Δ r ( k + M r ) = r ( k + M r + 1 ) − r ( k + M r ) = 0 (5)

定义误差信号 e ( k ) ∈ R q 为 e ( k ) = y ( k ) − r ( k ) 。本文针对系统(1)，研究预见跟踪问题，我们期望通过设计的控制器，使得系统的输出 y ( k ) 能够实现对参考信号 r ( k ) 的跟踪，即使得 lim k → ∞ e ( k ) = lim k → ∞ ( y ( k ) − r ( k ) ) = 0 能够成立。

本文需要如下引理：

引理1( [ 21 ] )：考虑矩阵 S = ( S 11 S 12 S 12 T S 22 ) ，其中 S 11 和 S 22 是对称矩阵并且可逆；则以下条件等价：

(1) S < 0 。

(2) S 11 < 0 ,     S 22 − S 21 S 11 − 1 S 21 T < 0 。

(3) S 22 < 0 ,       S 11 − S 12 S 22 − 1 S 12 T < 0 。

引理2( [ 22 ] )：对于适当维数的矩阵T，R，S，N及数 β ， T + S T R T + R S < 0 成立如果满足不等式 :

( T ( β R T + N S ) T β R T + N S − β N − β N T ) < 0

引理3( [ 23 ] )：对于矩阵L， H ∈ R s × s and H > 0 ，存在关系 − L T H − 1 L ≤ − L − L T + H 。

首先采用 [ 24 ] 中两项近似法，对系统(1)变形。然后采用预见控制理论中误差系统的方法来构造增广误差系统。

对导出的增广误差系统(15)，设计如下的输出反馈控制器：

Δ u ( k ) = K ( a ) Z ( k ) = ( ∑ i = 1 s a i K i ) Z ( k ) (16)

其中 K i 和 a i ( i = 1 , 2 , ⋯ , s ) 分别是不确定参数矩阵和可调节变量且 a i ≥ 0 ， ∑ i = 1 s a i = 1 。

基于(15)和(16)导出闭环系统：

{ x ¯ ( k + 1 ) = [ A ¯ ( θ ) + B ¯ ( θ ) K ( a ) C Z ( θ ) ] x ¯ ( k ) + A ¯ d ( θ ) x ¯ ( k − d 1 ) + A ¯ d ( θ ) x ¯ ( k − d 2 ) + η ¯ ( θ ) w ( k ) Z ( k ) = C Z ( θ ) x ¯ ( k ) (17)

基于设计的输出反馈控制器(16)，如果闭环系统是渐近稳定的，则有 lim k → ∞ x ¯ ( k ) = 0 ，注意到 e ( k ) 是闭环系统状态 x ¯ ( k ) 的一部分，即可得到 lim k → ∞ e ( k ) = 0 。因此我们设计一个包含参考信号预见行为的输出反馈控制器能够实现系统(1)的输出对参考信号的鲁棒跟踪目的。在(16)中，增益矩阵 能够通过解LMI方法获得。

为导出闭环系统渐近稳定的充分条件，考虑(17)的标称系统：

{ x ¯ ( k + 1 ) = [ A ¯ ( θ ) + B ¯ ( θ ) K ( a ) C Z ( θ ) ] x ¯ ( k )   + A ¯ d ( θ ) x ¯ ( k − d 1 ) + A ¯ d ( θ ) x ¯ ( k − d 2 ) Z ( k ) = C Z ( θ ) x ¯ ( k ) (18)

系统(18)渐近稳定的充分条件由如下的定理1给出。

定理1：给定正数 β ，可逆矩阵W和矩阵Q，系统(18)是鲁棒渐近稳定的如果存在一个正定矩阵 P ( θ ) ， S 1 ( θ ) ， S 2 ( θ ) ，矩阵 L ( a ) 和可逆矩阵 G ( a ) ，U使得

Π ( a , θ ) < 0 (19)

其中 Π ( a , θ ) 为下列矩阵

( − G ( a ) T − G ( a ) + P − 1 ( θ ) ∗ ∗ ∗ 0 − S 1 ( θ ) ∗ ∗ 0 0 − S 2 ( θ ) ∗ A ¯ ( θ ) G ( a ) + B ¯ ( θ ) L ( a ) Q A ¯ d ( θ ) A ¯ d ( θ ) − P − 1 ( θ ) G ( a ) 0 0 0 G ( a ) 0 0 0 C Z ( θ ) G ( a ) − U Q 0 0 0 ∗ ∗ ∗ ∗ ∗ ∗ ∗ ∗ ∗ ∗ ∗ ∗ − S 1 ( θ ) − 1 ∗ ∗ 0 − S 2 ( θ ) − 1 ∗ 0 0 − β U W − β W T U T )

证明：令 L ( a ) U − 1 = K ( a ) ， A ^ ( θ ) = A ¯ ( θ ) + B ¯ ( θ ) K ( a ) C Z ( θ ) ，由引理2，条件(19)能够保证：

( − G ( a ) T − G ( a ) + P − 1 ( θ ) ∗ ∗ ∗ ∗ ∗ 0 − S 1 ( θ ) ∗ ∗ ∗ ∗ 0 − S 2 ( θ ) ∗ ∗ ∗ A ^ ( θ ) G ( a ) A ¯ d ( θ ) A ¯ d ( θ ) − P ( θ ) − 1 ∗ ∗ G ( a ) 0 0 0 − S 1 ( θ ) − 1 ∗ G ( a ) 0 0 0 0 − S 2 ( θ ) − 1 ) < 0 (20)

又因为 − G ( a ) T − G ( a ) + P ( θ ) − 1 ≥ − G ( a ) T P ( θ ) G ( a ) ，所以

( − G ( a ) T P ( θ ) G ( a ) ∗ ∗ ∗ ∗ ∗ 0 − S 1 ( θ ) ∗ ∗ ∗ ∗ 0 0 − S 2 ( θ ) ∗ ∗ ∗ A ^ ( θ ) G ( a ) A ¯ d ( θ ) A ¯ d ( θ ) − P ( θ ) − 1 ∗ ∗ G ( a ) 0 0 0 − S 1 ( θ ) − 1 ∗ G ( a ) 0 0 0 0 − S 2 ( θ ) − 1 ) < 0 (21)

不等式(21)能够保证下式(22)成立(注：(22)式不等式左边作合同变换：前乘可逆对称阵 d i a g { G ( a ) , I , I , I , I , I } ，后乘它的转置阵，则可得到不等式(21)，故(21)与(22)等价)：

( P ( θ ) ∗ ∗ ∗ ∗ ∗ 0 − S 1 ( θ ) ∗ ∗ ∗ ∗ 0 0 − S 2 ( θ ) ∗ ∗ ∗ A ^ ( θ ) A ¯ d ( θ ) A ¯ d ( θ ) − P ( θ ) − 1 ∗ ∗ I 0 0 0 − S 1 ( θ ) − 1 ∗ I 0 0 0 0 − S 2 ( θ ) − 1 ) < 0 (22)

利用Schur补引理，由(23)式能够得到：

Ω = ( A ^ ( θ ) T P ( θ ) A ^ ( θ ) − P ( θ ) + S 1 ( θ ) + S 2 ( θ ) A ¯ d ( θ ) T P ( θ ) A ^ ( θ ) A ¯ d ( θ ) T P ( θ ) A ^ ( θ ) A ^ ( θ ) T P ( θ ) A ¯ d ( θ ) A ^ ( θ ) T P ( θ ) A ¯ d ( θ ) A ¯ d ( θ ) T P ( θ ) A ¯ d ( θ ) − S 1 ( θ ) A ¯ d ( θ ) T P ( θ ) A ¯ d ( θ ) A ¯ d ( θ ) T P ( θ ) A ¯ d ( θ ) A ¯ d ( θ ) T P ( θ ) A ¯ d ( θ ) − S 2 ( θ ) ) < 0 (23)

取Lyapunov函数：

V ( x ¯ ( k ) ) = x ¯ ( k ) T P ( θ ) x ¯ ( k ) + ∑ i = k − d 1 k − 1 x ¯ ( i ) T S 1 ( θ ) x ¯ ( i ) + ∑ i = k − d 2 k − 1 x ¯ ( i ) T S 2 ( θ ) x ¯ ( i ) (24)

其中 P ( θ ) = ∑ i = 1 s θ i P i ， S 1 ( θ ) = ∑ i = 1 s θ i S 1 i ， S 2 ( θ ) = ∑ i = 1 s θ i S 2 i 。

定义 Δ V ( x ¯ ( k ) ) = V ( x ¯ ( k + 1 ) ) − V ( x ¯ ( k ) ) 则有：

Δ V ( x ¯ ( k ) ) = ( x ¯ T ( k ) x ¯ T ( k − d 1 ) x ¯ T ( k − d 2 ) ) Ω ( x ¯ T ( k ) x ¯ T ( k − d 1 ) x ¯ T ( k − d 2 ) ) T (25)

由(23)，得到 Δ V ( x ¯ ( k ) ) < 0 ，则证得系统(17)是鲁棒渐近稳定的。

注意到定理1依赖不确定参数 θ ，以下定理2是在假设1的基础上，将定理1转换成一般LMI的形式，以便通过LMI工具箱得到相关矩阵的可行性解。

定理2：给定正数 β ，可逆矩阵W和矩阵Q，如果存在正定矩阵 X i ， Y 1 i ， Y 2 i ，矩阵 L j 和可逆矩阵 ，U ( i = 1 , 2 , 3 , ⋯ , s ; j = 1 , 2 , 3 , ⋯ , s ) 使得

Γ i j < 0             ( i , j ∈ { 1 , 2 , 3 , ⋯ , s } ) (26)

其中， Γ i j = ( − G j T − G j + X i ∗ ∗ ∗ 0 − Y 1 i ∗ ∗ 0 0 − Y 2 i ∗ A ¯ i G j + B ¯ i L j Q A ¯ d i Y 1 i A ¯ d i Y 2 i − X i G j 0 0 0 G j 0 0 0 C Z i G j − U Q 0 0 β W T L j T B ¯ i T ∗ ∗ ∗ ∗ ∗ ∗ ∗ ∗ ∗ ∗ ∗ ∗ − Y 1 i ∗ ∗ 0 − Y 2 i ∗ 0 0 − β U W − β W T U T )

则系统(17)是鲁棒渐近稳定，且增益矩阵 K i = L i U − 1 。控制器由下式给出：

Δ u ( k ) = K ( a ) Z ( k ) = ∑ i = 1 s a i L i U − 1 Z ( k ) (27)

证明：令 X i = P i − 1 ， Y 1 i = S 1 i − 1 ， Y 2 i = S 2 i − 1 。对(26)进行合同变换：左乘一个可逆对称阵 d i a g { I , S 1 i , S 2 i , I , I , I , I } 右乘它的转置阵，则有：

Π ( a , θ ) = ∑ i = 1 s ∑ j = 1 s θ i a j Π i j < 0 (28)

再根据定理1可以保证系统(17)渐近稳定，因此定理2得证。

基于定理1和定理2，我们给系统(1)设计具有预见跟踪性能的控制器由下面定理3给出。

定理3：如果假设1和假设2满足，且LMI(26)有一组可行性解，则设计的控制器为：

u ( k ) = K e ∑ i = 0 k − 1 e ( i ) + K y y ( k ) + ∑ j = 0 M r K r ( j ) r ( k + j ) (29)

其中，

K e = ∑ i = 1 s a i K e i ， K y = ∑ i = 1 s a i K y i ， K r ( j ) = ∑ i = 1 s a i K r i ( j ) ， K e i = ∑ j = 1 p K j i ， K i = L i U − 1 。

证明：基于定理1和定理2，可以得到期望的预见控制器，其中增益矩阵 K i 可通过解LMI(26)中的 L i     ( i = 1 , 2 , ⋯ , s ) 和U得到 ( K i = L i U − 1 ) ，即有：

Δ u ( k ) = K ( a ) Z ( k )     = ∑ i = 1 s a i   K i   Z ( k ) = ∑ i = 1 s a i ( K e i e ( k ) + K y i Δ y ( i ) + ∑ j = 0 M r K r i ( j ) Δ r ( k + j ) ) (30)

其中，

K i = [ K 1 i K 2 i ⋯ K p i K y i           K r i ( 0 ) K r i ( 1 ) ⋯ K r i ( M r ) ] = L i U − 1 (31)

注意到 K e i = ∑ j = 1 p K j i ， K e = ∑ i = 1 s a i K e i ， K y = ∑ i = 1 s a i K y i ， K r ( j ) = ∑ i = 1 s a i K r i ( j ) ，令增益矩阵 K = [ K e K y K r ( 0 ) K r ( 1 ) ⋯ K r ( M r ) ] ，则有 Δ u ( k ) = K e e ( k ) + K y Δ y ( k ) + ∑ j = 0 M r K r ( j ) Δ r ( k + j ) ，由 Δ u ( k ) = u ( k + 1 ) − u ( k ) ，知道 u ( k ) = ∑ i = − M r − 1 k − 1 Δ u ( i ) + u ( − M r − 1 ) 。进一步，在假定 x ( i ) = 0 ， u ( i ) = 0 ， r ( i ) = 0 ， ( i < 0 ) 下得到：

u ( k ) = K e ∑ i = 0 k − 1 e ( i ) + K y y ( k ) + ∑ j = 0 M r K r ( j ) r ( k + j ) (32)

定理3得证。

注2：矩阵 L i 和可逆矩阵U由不等式(26)决定。等式(31)决定着 K 1 i ， K 2 i ， ⋯ ， K p i ， K y i ， K r i ( 0 ) ， K r i ( 1 ) ， ⋯ ， K r i ( M r ) 以及 K i 之间的关系。在控制器(29下，系统(1)的闭环系统的输出 y ( k ) 能够准确的跟踪到参考信号 r ( k ) 。

注3：根据(29)，我们知道系统(1)的预见控制器包括三部分：第一部分是跟踪误差补偿，第二部分是输出反馈，最后一部分是可预见信息，表明控制器具备预见行为。

考虑系统(1)的鲁棒预见跟踪控制问题。令

A ( θ ) = ( 1 − 0.5 0.8 1 0 0 0 1 0 ) θ 1 + ( 1 − 0.8 0.5 1 0 0 0 1 0 ) θ 2 , A d ( θ ) = ( 0 0.01 0 0 0 0.1 0 0 0 ) θ 1 + ( 0.01 0 0 0 − 0.1 0 0 0 0 ) θ 2

B ( θ ) = ( 1 0.5 0 ) θ 1 + ( 1 1 0 ) θ 2 , C ( θ ) = [ 1 0 0 ] θ 1 + [ 1 0 0 ] θ 2 .

这里我们取 β = 0.6 ， a 1 = 0.3 ， a 2 = 0.7 ，并取矩阵 Q = 0.9 C z ( θ ) T C Z ( θ ) ， W = 0.6 C z ( θ ) 。不确定参数取 θ 1 = 0.4 ， θ 2 = 0.6 。本文仿真将从两种参考目标分别来考虑三种预见长度下的跟踪情况(① M r = 0 ; ② M r = 2 ; ③ M r = 6 )。

基于定理2，我们通过解LMI(24)得到矩阵变量 L 1 , L 2 以及U，则增益矩阵 K i = L i U − 1 ( i = 1 , 2 )能够计算得出。于是能够得到输出反馈增益矩阵 K ( a ) = a 1 K 1 + a 2 K 2 ，其中 K 1 , K 2 能从LMI(24)中得到。

若 M r = 0 ，则 K ( a ) = [ 0.0000 0.0000 − 0.3655 − 1.4296 1.0631 ] 从而

K = [ K e K y K r ( 0 ) ] = [ − 0.3655 − 1.4296 1.0631 ]

若 M r = 2 ，则 K ( a ) = [ − 0.0001 − 0.0001 − 0.3366 − 1.3943       0.2949 0.2717 0.2028 ]

K = [ K e K y K r ( 0 ) K r ( 1 ) K r ( 2 ) ] = [ − 0.3368 − 1.3943 0.2949 0.2717 0.2028 ]

若 M r = 6 ，则

K ( a ) = [ 0.0000 0.0000 − 0.3413 − 1.3946   0.3102 0.2971                   0.2414 0.1732 0.0825 0.0285 0.0100 ]

K = [ K e K y K r ( 0 ) K r ( 1 ) K r ( 2 )         K r ( 3 ) K r ( 4 ) K r ( 5 ) K r ( 6 ) ] = [ − 0.3413 − 1.3946 0.3102 0.2971 0.2414         0.1732 0.0825 0.0285 0.0100 ]

假定 x ( k ) 的初始状态为 x ( 0 ) = ( 0 0 0 ) T ，我们可以得到如下的仿真结果。

情形1：参考信号取为

r ( k ) = { 1 , k ≥ 30 0 , k < 30

设计输出反馈控制器为： u ( k ) = K e ∑ i = 0 k − 1 e ( i ) + K y y ( k ) + ∑ j = 0 M r K r ( j ) r ( k + j ) ，则仿真结果：

图1表示的是预见长度分别取为：① M r = 0 ；② M r = 2 ；③ M r = 6 时，闭环系统的输出跟踪参考信号(情形1)的情况；

图2表示的是预见长度分别取为：① M r = 0 ；② M r = 2 ；③ M r = 6 时闭环系统的控制输入；

图3是预见长度分别取为：① M r = 0 ；② M r = 2 ；③ M r = 6 时闭环系统的跟踪误差。

图1和图3表明在设计的预见控制器下，系统(1)的输出能够有效地跟踪到参考信号。同时动态特征(峰值 y max ，峰值时间 K p ，上升时间 K r 以及稳定时间 K s )能够反映出系统性能，仿真过程中可以发现闭环系统有一个期望的稳态响应特征：随着预见长度的增加，动态性能指数反而下降，即参考信号的预见长度越长，闭环系统的输出越能快速的跟踪到参考信号。

M r = 0 : y max = 1.2589 , K p = 34 , K r = 31 , K s = 58

M r = 2 : y max = 1.0292 , K p = 38 , K r = 29 , K s = 54

M r = 6 : y max = 1.0238 , K p = 38 , K r = 25 , K s = 53

针对同一时变时滞，我们重新选择参考信号：

情形2：参考信号取为

r ( k ) = { 0 ,       k ≤ 10 0.1 k − 1 , 10 < k < 30 5 ,       k ≥ 30

类似的，我们能够得到如下的仿真结果：

图4表示的是预见长度分别取为：① M r = 0 ；② M r = 2 ；③ M r = 6 时，闭环系统的输出跟踪参考信号(情形2)的情况；

图5表示的是预见长度分别取为：① M r = 0 ；② M r = 2 ；③ M r = 6 时闭环系统的控制输入；

图6是预见长度分别取为：① M r = 0 ；② M r = 2 ；③ M r = 6 时闭环系统的跟踪误差。

对于预见跟踪控制问题，利用渐近稳定性的定义得到闭环系统的状态渐近趋于平衡点0。对于系统(17)，条件(26)保证系统渐近稳定，则可以得到 lim k → ∞ x ¯ ( k ) = 0 。误差信号 e ( k ) 作为状态 x ¯ ( k ) 的一部分必然有 lim k → ∞ e ( k ) = 0 ，于是实现闭环系统输出跟踪到参考信号的目的。图3和图6是闭环系统的输出与参考信号的跟踪误差轨迹图，可以看出在达到一定时间，跟踪误差趋于0，实现了预见跟踪控制。

本文考虑时变时滞的情况，研究了预见跟踪控制问题。通过设计一个预见控制器使得闭环系统的输

图1. 参考信号(情形1)以及闭环系统的输出

图2. 闭环系统的控制输入

图3. 闭环系统的跟踪误差

图4. 参考信号(情形2)以及闭环系统的输出

图5. 闭环系统的控制输入

图6. 闭环系统的跟踪误差

出能够准确跟踪到参考信号。对于跟踪问题，通过引入误差信号来构造增广误差系统，使其闭环系统渐近稳定来达到跟踪目的，并利用LMI工具箱求出增益矩阵得到预见控制器。但是时变时滞的影响会破坏系统的稳定性，且不能直接构造增广误差系统。本文给出了在时变时滞的影响下闭环系统渐近稳定的充分条件，并提出了具有预见行为的输出反馈控制器的设计方法。最后数值仿真说明了本文结论的有效性。

在此对“国家自然科学基金”，基金号61603188的资助表示感谢。

范 蓉,王为群,姚 娟. 一类时变时滞不确定离散系统的预见跟踪控制Preview Tracking Control for a Class of Uncertain Discrete-Time System with Time-Varying Delay[J]. 动力系统与控制, 2018, 07(03): 201-213. https://doi.org/10.12677/DSC.2018.73022