同步(一致)行为是生物、生态、工程和社会科学等领域中最普遍的群聚现象之一. 在过去十几年里, 耦合系统中仅由局部交互引起的同步问题引起了大量研究者的关注^[1-5]. 在自然和工程系统中, 合作、竞争关系普遍存在, 且很多实际系统同时存在合作与竞争关系, 例如社会网络^[6]、存在合作与竞争的种群^[7]、竞争性细胞神经元^[8]和个性化推荐^[9]. 为了描述系统中的合作与竞争关系, 研究者们引入了符号图, 其中正数边表示合作关系, 负数边表示竞争关系.

目前, 越来越多的研究人员开始利用符号图来研究网络中的各种群聚现象^[10-16]. 在文献[10]中, Altafini研究了定义在符号图上的一个积分器网络, 并得到了关于双向一致的一些定理. 这里的双向一致表示所有的智能体都收敛到一个模量相等、符号不同的值. 其中, 作者假设符号图是结构平衡的, 即所有节点可以被分为两个阵营, 每个阵营内部是合作关系, 两个阵营之间是竞争关系. 这个假设对双向一致性结论的得出至关重要. 文献[10]的结论推广到了更一般的线性多智能体系统^[11-13], 其中每个智能体都由一个线性时不变系统表示. 例如对于有向图上的积分器网络, 文献[13] 在符号图含有生成树的情况下得到了达到双向一致的一些充分条件. 很多研究者陆续对各种特定网络展开了双向同步问题研究, 例如双向聚集^[14]、区间双向一致^[15]、含有时滞的双向一致^[16]等. 基于压缩性分析, 文献[17] 研究了耦合非线性网络的双向同步问题. 对于耦合离散系统构成的网络, 其双向同步问题也受到了很多研究者的关注^[18-19]. 对于更多的关于双向同步的研究, 可以参见综述文献[20-21].

在实际系统中, 随着时间的推移, 网络的拓扑结构可能会发生变化. 而且, 网络所形成的符号图可能不满足结构平衡特性. 例如, 在社会网络中, 个体之间的关系可能会由合作(友谊)到竞争(敌意)变化, 反之亦然; 在多党制的国家, 很多成员经常会从一个党派转向另一个党派. 当符号图不满足结构平衡性时, 网络不能达到双向同步. 在文献[22]中, 作者利用矩阵的最终为正性质, 分别研究了连续和离散时间舆论动力学模型的动力学行为. 当符号图随着时间变化的时候, 网络构成一个切换系统. 文献[23-24]考虑了所有符号图在结构上都是平衡的, 且敌对阵营的成员随着时间的推移是不变的情况. 具体地, 在文献[23]中, 作者得到了使非线性系统达到模同步的充分条件; 在文献[24]中, 作者设计了一种牵引控制, 使闭环系统实现双向同步. 如果这些符号图中的节点随着时间变化, 那么双向同步将不可能达到.

本文将研究含有对抗性关系和时变拓扑的耦合离散系统的有界双向同步(Bounded bipartite synchronization, BBS)问题. 考虑以下情形: 1)在某些时刻, 所有个体不能被分为两个敌对阵营; 2)虽然所有个体可以被划分为两个阵营. 但所形成敌对阵营中的成员会随时间改变. 当情形1)和2)出现时, 将这种耦合离散系统看成是一个特定网络的扰动, 在这个特定网络中, 所有的个体都可以被分成两个敌对阵营, 且二者中的成员随着时间的推移会保持不变. 在该特定网络的所有符号图都是连通的条件下, 本文得到了使系统达到有界双向同步的一些充分条件. 最后, 利用一个数值例子来说明所得结论的有效性.

本文符号说明如下: $ \vert x \vert $表示实数$ x $的绝对值, $ {\bf Z}^+ $表示正整数域, $ \vert\vert {\boldsymbol y}\vert\vert $表示向量$ {\boldsymbol y} $的范数, $ I_N $表示$ N $维单位矩阵, $ {\boldsymbol 1}_N $表示元素都为$ 1 $的$ N $维列向量, 运算符$ \otimes $表示Kronecker积. 对于矩阵$ A $, 符号$ \lambda _{\min}(A) $, $ \lambda_{\max}(A) $分别表示矩阵$ A $的最小特征值和最大特征值. ${\rm diag}\{{\cdot}\}$表示一个对角矩阵, $ {\rm sgn}(\cdot ) $代表符号函数. 如果对于每个固定的$ s $, 函数$ \beta \left( {r,\;s} \right) $是严格递增的且$ \beta \left( {\rm{0},\;s} \right)\equiv \rm{0} $, 对于每个固定的$ r, $ 函数$ \beta \left( {r,\;s} \right) $是严格递减的且$\lim\nolimits_{s\to \infty } \beta \left( {r,\;s} \right) = 0, $那么函数$ \beta \left( {r,\;s} \right) $称为 KL类函数.

1.   问题描述

考虑包含$ N $个离散系统的网络

其中, $ i = 1,2,\cdots ,N. $$ x_i \in {\bf R}^n $是第$ i $个节点的状态, A, B 是常数矩阵, $ u_i(k) $是控制输入. 假设网络的拓扑在$ p $个无向符号图$G\left( {E^k}\right)$ (符号图定义见附录A), $k = 1,2,\cdots ,p$之间切换, 其中切换信号是$\sigma (k):{\bf Z}^+\to $$ P: = \{1,\;2,\;\cdots \;,\rm{}p\}$, 它是一个分段右连续的函数. 控制输入$ u_i (k) $设计为

其中, $ K $是一个需要设计的增益矩阵, $ e_{ij} $是图$ G({E^k} ) $的边值. 令$x = [{x_1^{\rm T} \;\cdots\;x_N^{\rm T}} ]^{\rm T}$, $\{k_i:i = 0,1,\cdots\}$ 是$ \sigma (k) $的切换时刻. 存在正常数$ T>1 $, 使得$ k_{i+1} -k_i \ge T , $$ \forall i\ge 0. $

注1. 网络在切换信号下构成一个切换系统. 本文中要求存在正常数$ T>1 $, 使得$ k_{i+1} -k_i \ge T $, $ \forall i\ge 0 $. 这里的$ T>1 $可以看成是驻留时间. 如果没有驻留时间, 那么在有限时间内可能会有无限次切换, 对于系统的收敛性会有很大影响.

通常来说, 如果符号图结构平衡, 那么其所有节点可以划分为两个敌对阵营, 其中每个阵营中的个体之间的关系是合作的, 属于不同阵营的个体之间的关系是对立的. 对于符号图$ G({E^k} ),\,k = 1,2,\cdots,p $, 可能存在以下情况: 1)虽然每一个符号图都满足结构平衡, 即每个符号图都可以划分为两个敌对阵营, 但是每一个符号图的两个敌对阵营中的个体是不一样的, 例如在多党派执政的国家, 一些个体随着时间变化从一个阵营转移到另一个阵营; 2)可能存在某些不满足结构平衡的符号图. 在这些情况下, 网络很难达到双向同步. 为了研究这两种情况下的网络的同步问题, 将这些符号图看成是某些特定结构平衡符号图的扰动. 具体地, 假设符号图$ G( {E^k} ) $的邻接矩阵可以分为两个邻接矩阵, 即$ E^k = \bar{E}^k+w^k $, 其中, $ \bar {E}^k $是关于符号图$ G({\bar{E}^k}) $的一个邻接矩阵. 把控制输入(2)中的符号图改为$ G( {\bar {E}^k} ) $可以得到一个新的输入

因此, 由符号图$ G\left( {E^k}\right) $形成的耦合系统(1)和(2)可以看成是由符号图$ G( {\bar {E}^k}) $形成的耦合系统(1)和(3)的扰动. 而且, 假设符号图$ G( {\bar{E}^k} ) $, $ k = 1,2, \cdots , $$ p $的节点$ \{1,2,\cdots ,N\} $可以划分为两个敌对阵营$ V_1 $, $ V_2 $, 且存在一个符号矩阵$ \Psi\; (\Psi = {\rm diag}\{\sigma _1 ,\cdots ,\sigma _N \},$$ \,\sigma _i \in \{\pm 1\}) $, 使得矩阵$ \Psi \bar {E}^k\Psi $, $ k = 1,2,\cdots ,p $都是非负矩阵.

接下来, 本文将研究当控制输入为式(2)时, 网络(1)将在何种条件下达到有界双向同步. 双向同步和有界双向同步的定义分别如下.

定义1. 如果存在依赖于非零初始条件的函数$ \zeta(k)\ne 0, $ 使得以下条件成立: $\lim\nolimits_{k\to \infty }( {x_i(k)-\zeta(k)} ) = $$ 0 ,$$ \forall i\in V_1, \lim\nolimits_{k\to \infty } \left({x_i (k)+\zeta(k)}\right) = 0, \forall i\in V_2, $ 那么控制输入为式(3)的网络(1)达到双向同步.

定义2. 如果满足以下两个条件, 那么控制输入为式(2)的网络(1)达到有界双向同步: 1)网络(1)在形式为式(3)的控制输入下达到双向同步; 2)存在一个正常数$ \xi $(依赖于非零初始条件), 一个KL类函数$ \beta (\cdot ,\cdot ) $(依赖于图$ G( {E^k}) $, $ k = 1,2,\cdots ,p )$, 使得${\vert \vert }\delta (k){\vert \vert }\le \beta $$ \left(\vert \vert \delta(0)\vert \vert , t\right)+\xi $成立, 其中$ \delta (k) = \;x(k)-$$\frac{1}{N}{\rm {\bf 1}}_N \otimes {\rm {\bf 1}}_N^{\rm T} \otimes I_n x(k). $

2.   主要结论

本节将研究以下两种情形: 1)在某些时刻, 所有个体不能划分为两个敌对阵营; 2)虽然所有个体可以划分为两个阵营, 但形成的敌对阵营中的成员会随时间改变. 如果符号图$ G( {\bar {E}^k} ) $, $ k = 1,2,\cdots ,\;p $都是连通的, 那么可以得到条件使得控制输入为式(2)的网络(1)达到有界双向同步. 为此, 给出以下假设:

假设1. 假设矩阵$ A $的所有特征值是模为1的半单特征值, 即所有约当块都是一维的.

进而, 针对存在对抗关系和时变拓扑的耦合离散系统, 可以得到定理1.

定理1. 考虑网络(1), 假定假设1成立且符号图$G( {\bar {E}^k})$, $ k = 1,2,\cdots ,p $连通. 如果存在$ \mu $使得不等式(4)成立(其中$ \Delta ^j = L^j-\bar {L}^j $),

那么控制输入为式(2)的网络(1)在$ K = \mu B^{\rm T}P^{\rm T}PA $时可以达到有界双向同步, 其中, $ \bar {A} = PAP^{-1} $, $ \bar {B} = $$ PB $, 可逆矩阵$ P $使得$ \bar {A} $是$ A $的约当标准型. 而且, 其最终界为$ \xi = \sqrt {\frac{\sigma _2 }{\sigma _1 }} \frac{\left\| {x(0)} \right\|( {\rm{1+}\sqrt {\rm{1+}\theta \alpha } } )}{\theta \alpha } ,$ 其中 $ 0< \theta < $$ 1 ,\; \alpha =\frac{\mu \lambda _2 \lambda _{\min } ( {\bar {A}^{\rm T}\bar {B}\bar {B}^{\rm T}\bar {A}} )}{2},\; \sigma _1 = \lambda _{\min },\; \sigma _2 = \lambda _{\max }\; (P^{\rm{T}}P), $ $ \lambda _2 = $$ \min _{k = 1,2,\cdots ,p} \lambda _2( {\bar {L}^k}). $

证明. 选择$ K = \mu B^{\rm T}P^{\rm T}PA $, 则控制输入为式(2)的网络(1)变为

其中, $ i = 1,2,\cdots ,N $. 式(5)可以写成如下所示的紧凑形式.

令$ \bar {x}(k) = \left( {I_N \otimes P} \right)x(k) $, 那么

令$ V_1 \left( {\bar {x}(k)} \right) = \frac{1}{2}\bar {x}^{\rm T}(k)\bar {x}(k) $, 那么$ V_1 $沿着式(7)的差分满足

基于条件(4), 可得

即得$ V_1 \left( {\bar {x}(k)} \right) $是非递增的, 且有$ \left\| {\bar {x}(k)} \right\|\le \left\| {\bar {x}(0)} \right\| $. 令$ \bar {y}(k) = \left( {\Psi \otimes P} \right)x(k) $, 在控制输入为式(3)时, 形成的闭环系统可表示为

由于图$ G( {\bar {E}^k} ) , $$k = 1,2,\cdots ,p$ 的节点$\{i = 1,2, \cdots , $$ N\}$可划分为两个敌对阵营$ V_1 $和$ V_2 $, 且图$ G( {\bar {E}^k}), $$ k = 1,$$ 2,\cdots ,p $是连通的, 基于定理1^[25], 可知网络(1)和(3)在任意切换信号下达到双向同步.

令$ z(k) = \left( {\Psi \otimes P} \right)x(k) $, 则控制输入为式(2)的网络(1)可表示为

令$ z_c (k) = \frac{1}{N}\sum\nolimits_{j = 1}^N {z_j (k)}, \bar {\delta }_i = z_i (k)-z_c (k) $, 则有$ \bar {\delta }(k) = $$ \left( {\Psi \otimes P} \right)\delta (k), $ 其中$\bar {\delta } =[ {\bar {\delta }_1^{\rm T}\, \cdots \,\bar {\delta }_N^{\rm T} } ]^{\rm{T}},$易得

令$ V_2( {\bar {\delta }(k)}) = \frac{1}{2}\bar {\delta }^{\rm T}(k)\bar {\delta }(k) $, 那么$ V_2 $沿着式(12)的差分满足

其中, 不等式第1部分可由条件(4)得到. 由于图$ G( {\bar {E}^k} ) $, $ k = 1,2,\cdots ,p $是连通的, 因而存在正交矩阵$ Q^{\sigma \left( k \right)}\in$$ {\bf R}^{N\times N} $, 使得

其中, $0 = \lambda _1^{\sigma (k)} < \lambda _2^{\sigma (k)}\,\le\,\cdots \,\le\, \lambda _N^{\sigma (k)}, Q^{\sigma (k)} = [q_1^{\sigma (k)} ,$$ q_2^{\sigma (k)} ,\cdots, q_N^{\sigma (k)}] $, $ q_1^{\sigma (k)} = \frac{{\rm {\bf 1}}_N }{\sqrt N } $是特征值$ \lambda _1^{\sigma (k)} \rm{ = 0} $对应的特征向量. 令$\bar {\delta }(k) = ( Q^{\sigma (k)}\otimes I_n )\hat {\delta }(k)$, 由$(( {Q^{\sigma (k)}} )^{\rm T}\otimes $$ I_n ) ( Q^{\sigma (k)}\otimes I_n )= I_{nN}$, 可得$ \bar {\delta }^{\rm T}\bar {\delta } = \hat {\delta }^{\rm T}\hat {\delta } $. 又由于$ \bar {\delta }_1 = $$ ( q_1^{\sigma (k)} \otimes I_n ) \hat {\delta }(k) = 0 $, 则可得

因此,

其中, $ 0<\theta <1 $. 所以下面的关系成立:

由$\delta (k) = ({\Psi \otimes P^{-1}})\bar {\delta }(k)z(k) = \left( {\Psi \otimes I_n } \right)\bar {x}(k) = $$\left(\Psi \;\otimes\; P \right) x(k) ,\;\left\| {\bar {x}(k)} \right\|\;\le \;\left\| {\bar {x}(0)} \right\|$, 可得 $\lambda _{\min }( {P^{\rm T}P})\times $$ \left\| {\delta (k)} \right\|^2\le$$ \| {\bar {\delta }(k)} \|^2 \le \lambda _{\max }( {P^{\rm T}P} )\left\|{\delta (k)}\right\|^2 $ 以及$\left\|{z(k)} \right\|^2\le $$ \lambda _{\max } ( {P^{\rm T}P} )\left\| {x(0)} \right\|^2$. 若不等式条件(15)成立, 则式(14)的左边部分成立.

当式(16)成立时, 式(15)成立.

因此, 对于$ \forall \left\| {\bar {\delta }(k)} \right\|\ge \frac{\sqrt {\sigma _2 } \left\| {x(0)} \right\|\left( {1+\sqrt {1+\theta \alpha } } \right)}{\theta \alpha } $,

为了应用引理1 (证明见附录B), 取$c_1 = c_2 = $$ {1}/{2},$$ c_3 = -\left( {1-\theta } \right)\alpha, $ $ c = \frac{\sqrt {\sigma _2 } \left\| {x(0)} \right\|\left( {1+\sqrt {1+\theta \alpha } } \right)}{\theta \alpha } $. 因此, 存在正常数$ \rho \ge 1, $$ 0<\gamma <1, $使得对于每个初始状态$ x(0) $, 网络(1)和(2)的解满足

由于$ \sigma _1 \left\| {\delta (0)} \right\|^2\le \left\| {\bar {\delta }(0)} \right\|^2\le \sigma _2 \left\| {\delta (0)} \right\|^2 $, 可得

从而得到控制输入为式(2)的网络(1)达到有界双向同步. □

注2. 由定理1的证明过程可以看出, 最终界为$\sqrt \frac{\sigma _2 }{\sigma _1 }\frac{\left\| {x(0)} \right\|\left( {1+\sqrt {1+\theta \alpha } } \right)}{\theta \alpha }.$ 因此, 为了使最终界比较小, 可以选择使$ \left\| {x(0)} \right\| $很小或者$ \alpha $很大的初始条件.

注3. 在定理1中, 假设矩阵$ A $的所有特征值是模为1的半单特征值, 即所有约当块都是一维的. 在这种假设条件下, 矩阵$ A $是正交矩阵, 即$ A^{\rm T}A = I. $这时矩阵$ A $是中立稳定的.

3.   数值例子

本节将给出一个数值例子来验证所得结论的有效性.

例1. 对于网络(1), 令$ N = 4 $, 其中矩阵$ A, B $为

因为矩阵$ A $是正交的, 所以假设1成立. 定义切换信号$ \sigma(k) $如式(21), 其中$ s\in {\bf Z}^+. $

假设有两个无向图$ G({E^i}) $, $ i = 1,2, $如图1所示, 图$ G( {E^2}) $的节点不能划分为两个敌对阵营$ V_1 $和$ V_2 $. 假设$ G( {\bar{E}^1}) $和$ G({\bar{E}^2}) $分别对应于图2(a)和图2(b). 可知图$ G( {\bar {E}^i} ) $, $ i = 1,2 $的节点能划分为两个敌对阵营$ V_1 = \{1,2\} $, $ V_2 = \{3,4\} $.

图 1  无向图$G({E^i})$, $i = 1,2$

Fig. 1  The undirected signed graph $G({E^i})$, $i = 1,2$

下载: 全尺寸图片 幻灯片

图 2  无向图$G( {\bar {E}^i} )$, $i = 1,2$

Fig. 2  The undirected signed graph $G( {\bar {E}^i} )$, $i = 1,2$

下载: 全尺寸图片 幻灯片

对于图 $ G({\bar{E}^i}), $ $ i \;=\; 1,\;2, $可选择符号矩阵 $\Psi = $$ {\rm diag}\{1,1,-1,-1\} $使得$ \Psi\bar{E}^k\Psi $, $ k = 1,2 $是非负矩阵, 根据其拉普拉斯矩阵$ \bar{L}^1 $和 $ \bar{L}^2 $

可得$ \lambda _2 = \min _{k = 1,2} \lambda _2( {\bar {L}^k} ) = 1.5858 $. 又由于$ \sigma _1 =$$ \sigma _2 = $ $1 , \lambda _{\min } ({\bar {A}^{\rm T}\bar {B}\bar{B}^{\rm T}\bar{A}}) = 3.5570 $. 如果选择$ \mu = \rm{0.}1, \theta = 0.9,$ $ \alpha = 0.2820, \xi = 0.6272, $ 那么图$ G({\bar {E}^i}) $, $ i = 1,2 $描述的网络(1)和(3)在切换信号$ \sigma (k) $下的状态演变如图3所示, 可知该网络达到双向同步. 对于网络(1), 在图$ G({E^i}) $, $ i = 1,2 $和切换信号$ \sigma (k) $下的时间演变图如图4所示, 根据定理1, 控制输入为式(2)的网络(1)能达到有界双向同步, 且终值为$ \xi = 0.6272 $. 在图5中, 明确地描述了范数误差和最终界.

图 3  四智能体网络在拓扑为图2、切换信号为$ \sigma (k) $时的时间演变过程

Fig. 3  Time evolution of 4-agent network with topologies in Fig.2 and switching signal $ \sigma(k) $

下载: 全尺寸图片 幻灯片

图 4  四智能体网络在拓扑为图1、切换信号为$ \sigma (k) $时的时间演变过程

Fig. 4  Time evolution of 4-agent network with topologies in Fig.1 and switching signal $ \sigma(k) $

下载: 全尺寸图片 幻灯片

图 5  四智能体网络在切换信号$ \sigma(k) $下的范数误差和终值

Fig. 5  Norm error of the 4-agent network with switching signal $ \sigma(k) $

下载: 全尺寸图片 幻灯片

4.   结论

当存在对抗关系和切换拓扑时, 本文研究了耦合离散线性系统的同步问题. 针对实际中可能存在的两种情形, 研究了耦合离散系统的有界双向同步问题, 得到了使闭环系统在任意切换信号下达到有界双向同步的充分条件. 数值仿真验证了本文所得理论的正确性. 本文的结论对于系统矩阵有一定的要求, 后续工作将考虑更一般的情况.

附录 A.   符号图

符号图$ G\left( {V,\varepsilon } \right) $由一个有限节点集和一个边集组成, 节点集记为$ V = \{1,2,\cdots ,N\} $, 边集记为$\varepsilon = \{ \left( {i,j} \right):i\ne j, i, j\in $$ V \}\subseteq V\times V$. 令$ E = ( {e_{ij} } ) $是图$ G $的一个邻接矩阵, 利用$ G\left( E \right) $来表示邻接矩阵为E 的符号图, 图$ G\left( E \right) $的拉普拉斯矩阵定义为$ L= C_r-E , $其中$C_r = {\rm diag}\{ \sum\nolimits_{j = 1}^N{| {e_{1j} } |, \cdots , \sum\nolimits_{j = 1}^N {| {e_{Nj} }|} } \}$. 由 i到j 的边$ \left( {i,j} \right)\in \varepsilon $是有向边, 其中节点$ i ,j$分别称为父节点和子节点. 如果$ \left( {j,i} \right),\left( {i,j} \right)\in \varepsilon, $ 那么图$ G\left( E \right) $是无向图. 文中定义$\varepsilon ^+ \,=\, \{ \left( {i,j} \right)\vert e_{ij} \, > \,0\},\ \varepsilon ^- = \{ \left( {i,j} \right) \vert e_{ij} < 0\} , \varepsilon = $$ \varepsilon ^+\cup\varepsilon^-.$ 由不同节点$ \left( {i_1 ,i_2 } \right), \left( {i_2 ,i_3 } \right), \cdots, $$ \left( i_{l-1} , \right. $$ \left. i_l \right) $所组成的边的一个序列称为路径(路径长度为$ l-1 )$. 若符号图中的任意两个不同节点之间存在路径, 则该图称为是强连通的. 已知包含相同节点集的$ p $个符号图$ G\left( {E^k} \right) = $$ ( {V,\varepsilon _k ,E^k}) $, $ k = 1,2,\cdots, p $, 则在切换信号$ \sigma(k) $下, 可以定义一个时变符号图, 即$G( {E^{\sigma (k)}} ) = $$ ( {V,\varepsilon ^{\sigma (k)},E^{\sigma(k)}} ) $.

附录 B.   预备引理

考虑如下差分方程

其中, $ x\in {\bf R}^n $, $ f:{\bf R}^n\to {\bf R}^n $是连续的, $ f(0) = 0 $.

引理1. 令$ V:{\bf R}^n\to {\bf R}^n $是一个连续函数, 且满足

其中, $ \forall k\ge 0 $, $ \forall x\in {\bf R}^n $, $ c,c_1 ,c_2 ,c_3 $是正常数. 那么, 对每个初始状态$ x(0) $, 存在正常数$ \rho \ge 1 $, $ 0<\gamma <1 $, 有$ T\ge 0 $(取决于$ x(0) $和$ c) $, 使得系统(B1)的解满足

证明. 本引理证明类似于定理4.18^[26]的证明. 令$\Omega _c = $$ \{ x\in {\bf R}^n\vert V(x)\le $$ c \},$若初始$ x(0)\in \Omega $, 则系统 (B1) 的解依赖于$\Omega _c $, 这是因为$ V(x(k)) $在边界上是负的. 对于$ {\bf R}^n-\Omega _c $内部的某个解, 令$ T $是它进入$ \Omega _c $的起始时刻, 则对于所有的$ k\in \left[ {0,T} \right]\cap {\bf Z}^+ $, 有下式成立:

因此,

又由于$ V\left( {x(k)} \right)\ge 0, $易得${c_3 }/{c_2 } < 1.$ 所以$( {1-{c_3 }/{c_2 }} ) < 1 .$ 可以得到

令$\rho = \sqrt {{c_2 }/{c_1 }}$, $\gamma = \sqrt {1-{c_3 }/{c_2 }}$, 则可以得到