部分时变离散系统中的稳定性判据
1.Lyapunov稳定性理论
下面先给出Lyapunov稳定性的一些基本理论(网上资源较多这里不再过多赘述):
2.一类时变离散系统的稳定性
- 定理
对于离散时变系统 x ( k + 1 ) = A ( k ) x ( k ) x(k+1)=A(k)x(k) x(k+1)=A(k)x(k),若 lim k → ∞ A ( k ) = A \lim_{k \rightarrow \infty}A(k)=A limk→∞A(k)=A存在且 A A A的特征值 λ i \lambda_i λi满足 ∣ λ i ∣ < 1 |\lambda_i|<1 ∣λi∣<1,则系统 x ( k + 1 ) = A ( k ) x ( k ) x(k+1)=A(k)x(k) x(k+1)=A(k)x(k)渐进稳定到 0 \mathbf{0} 0。
-
证明
先引入以下引理:
若常阵 A A A特征值满足 ∣ λ i ∣ < 1 |\lambda_i|<1 ∣λi∣<1,则对任意的正定矩阵 Q Q Q,都存在唯一的正定阵 P P P,满足 A T P A − P = − Q A^TPA-P=-Q ATPA−P=−Q。
下面开始证明:
若 lim k → ∞ A ( k ) = A \lim_{k \rightarrow \infty} A(k)=A limk→∞A(k)=A,则令 A ( k ) = A + B ( k ) A(k)=A+B(k) A(k)=A+B(k),则有: lim k → ∞ B ( k ) = 0 \lim_{k\rightarrow \infty} B(k)=0 limk→∞B(k)=0。有引理知 ∃ 正定阵 P \exists 正定阵P ∃正定阵P满足 A T P A − P = − I A^TPA-P=-I ATPA−P=−I,其中 I I I为单位阵。
下面构造 V ( k ) = x ( k ) T P x ( k ) V(k)=x(k)^TPx(k) V(k)=x(k)TPx(k),则 V ( k ) V(k) V(k)正定。同时有:
V ( k + 1 ) = x ( k + 1 ) T P x ( k + 1 ) = x ( k ) T A ( k ) T P A ( k ) x ( k ) V(k+1)=x(k+1)^TPx(k+1)\\=x(k)^TA(k)^TPA(k)x(k) V(k+1)=x(k+1)TPx(k+1)=x(k)TA(k)TPA(k)x(k)
令 Δ V ( k ) = V ( k + 1 ) − V ( k ) \Delta V(k)=V(k+1)-V(k) ΔV(k)=V(k+1)−V(k),得到:
Δ V ( k ) = V ( k + 1 ) − V ( k ) = x ( k ) T ( A ( k ) T P A ( k ) − P ) x ( k ) = x ( k ) T ( ( A + B ( k ) ) T P ( A + B ( k ) ) − P ) x ( k ) = x ( k ) T ( A T P A − P + C ( k ) ) x ( k ) = − x ( k ) T x ( k ) + x ( k ) T C ( k ) x ( k ) \Delta V(k)=V(k+1)-V(k)\\=x(k)^T(A(k)^TPA(k)-P)x(k)\\=x(k)^T((A+B(k))^TP(A+B(k))-P)x(k)\\=x(k)^T(A^TPA-P+C(k))x(k)\\=-x(k)^Tx(k)+x(k)^TC(k)x(k) ΔV(k)=V(k+1)−V(k)=x(k)T(A(k)TPA(k)−P)x(k)=x(k)T((A+B(k))TP(A+B(k))−P)x(k)=x(k)T(ATPA−P+C(k))x(k)=−x(k)Tx(k)+x(k)TC(k)x(k)
其中: C ( k ) = A T P B ( k ) + B ( k ) T P A + B ( k ) T P B ( k ) C(k)=A^TPB(k)+B(k)^TPA+B(k)^TPB(k) C(k)=ATPB(k)+B(k)TPA+B(k)TPB(k)。有:
lim k → ∞ C ( k ) = 0 \lim_{k\rightarrow \infty}C(k)=0 k→∞limC(k)=0
此时有: lim k → ∞ x ( k ) T C ( k ) x ( k ) x ( k ) T x ( k ) = 0 \lim_{k \rightarrow \infty}\frac{x(k)^TC(k)x(k)}{x(k)^Tx(k)}=0 limk→∞x(k)Tx(k)x(k)TC(k)x(k)=0,由极限的定义有:
∃ k 0 , 当 k > k 0 时有 ∣ x ( k ) T C ( k ) x ( k ) ∣ < 1 2 x ( k ) T x ( k ) \exists k_0,当k>k_0时有|x(k)^TC(k)x(k)|<\frac{1}{2}x(k)^Tx(k) ∃k0,当k>k0时有∣x(k)TC(k)x(k)∣<21x(k)Tx(k)
因此 k > k 0 k>k_0 k>k0时有 Δ V ( k ) = − x ( k ) T x ( k ) + x ( k ) T C ( k ) x ( k ) < − 1 2 x ( k ) T x ( k ) \Delta V(k)=-x(k)^Tx(k)+x(k)^TC(k)x(k)<-\frac{1}{2}x(k)^Tx(k) ΔV(k)=−x(k)Tx(k)+x(k)TC(k)x(k)<−21x(k)Tx(k)。 故 Δ V ( k ) \Delta V(k) ΔV(k)在 k > k 0 k > k_0 k>k0时负定。
证毕。