• 대한전기학회
Mobile QR Code QR CODE : The Transactions of the Korean Institute of Electrical Engineers
  • COPE
  • kcse
  • 한국과학기술단체총연합회
  • 한국학술지인용색인
  • Scopus
  • crossref
  • orcid




Stability, time-delay, LKF, constrained quadratic function, equivalent transform, LMI

1. 서 론

시간 지연은 실제 시스템에서는 피할 수 없는 존재이기고 또한 시스템의 안정성에도 지대한 영향을 미치기에, 지난 수십 년 동안 이론 및 응용에서 매우 활발히 다루어진 분야이다[11,16,17].

다음으로 기술되는 시변 시간지연을 갖는 시간지연 선형 시스템을 생각하자.

(1)
$$ \left\{\begin{array}{l} \dot{x}(t)=A x(t)+A_{1} x(t-d(t)) \\ x(\theta)=\phi(\theta), \theta \in[-h, 0] \end{array}\right. $$

여기서 $x(t)\in R^{n}$은 시스템의 상태, $A ,\: A_{1}\in R^{n\times n}$은 상수 시스템 행렬들, $\phi(\theta)$는 상태의 초기 조건, 그리고 시변 시간 지연 $d(t)$는 다음을 만족한다.

(2)
$0\le d(t)\le h,\:\mu_{1}\le\dot d(t)\le\mu_{2}< 1 .$

시간 지연(2)를 갖는 시스템(1)의 안정성을 보장하는 좀 더 나은 조건을 구하기 위해 다음의 세 가지 요인이 중심 요소이다: (i) 적당한 Lyapunov-Krasovskii 함수(LKF)의 선정. (ii) 적당한 적분 부등식의 선정. (iii) LMI형태로의 변환. 특히 적분 부등식의 경우 많은 진보가 있어왔다. 초기의 Jensen 부등식[17]으로 시작하여 Wirtinger 부등식 기초한 적분부등식[2]의 출현으로 많은 진보가 있었고, 이어 Bessel-Legendre 부등식에 기초한 적분부등식[5]의 출현으로 더 많은 진보가 있었다. 그 외에 reciprocally convex 부등식[1]과 이를 4개의 추가 변수를 이용한 적분부등식[5], 그리고 Free-weighting matrices 적분부등식[4]등이 안정성 해석에 많은 기여를 하였다.

또한, 고차의 Bessel-Legendre 부등식에 기초한 적분부등식[10]을 사용하거나, LKF에 포함된 특수 항으로 인하여 LKF의 시간 미분 항이 음이 되는 LMI 조건에 시간지연 $d(t)$에 대한 affine 함수 형식, 즉 $d(t)X_{1}+ X_{0}\le 0$, 이 아닌 2차 다항식 형식, 즉 $d^{2}(t)X_{2}+ d(t)X_{1}+ X_{0}\le 0$, 의 조건이 요구되어 이에 대한 충분조건이 제시되었다[6]. 또한 이를 개선하기 위한 여러 시도가 있었고, 최근에는 이의 필요충분조건이 제시되어[12-14] 안정성에 좀 더 개선이 이루어졌다.

본 논문에서는 첫째, 시스템 궤적에 따른 시간-미분이 시간지연 $d(t)$에 대한 2차 다항식의 형태와 2차함수의 적분의 합 형태를 갖는 LKF를 제시한다. 둘째, 2차함수의 적분항에 Bessel-Legendre 부등식[10]과 reciprocally convex 부등식[1]을 연속 적용한다. 셋째, 새로이 제시된 구속된 2차함수를 LMI로 등가 변환하는 필요충분조건을 이용하여 LKF의 시간-미분의 상한을 LMI형태로 변환하여, 안정성을 보장하는 LMI 형태의 충분조건을 제시한다. 끝으로, 대표적인 예제 2개를 통하여, 본 논문의 결과는 기존의 안정성 결과보다 향상된 결과임을 보인다.

2. 예비 결과

다음의 보조정리들은 주요 결과 유도에 사용될 예비 결과이다.

보조정리 1 [2]: $t_{1}< t_{3},\: t_{2}\in(t_{1},\: t_{3})$라 하고, 양확정 행렬 $R\in R^{n\times n}$에 대하여 $\hat R =\mathrm{d iag}\{R ,\: 3 R ,\: 5R\}$라 하면, 다음이 항상 성립한다.

(3)
$-\int_{t_{1}}^{t_{3}}\dot x^{T}(s)R\dot x(s)ds\le -\dfrac{1}{t_{3}- t_{2}}q^{T}(t_{2},\: t_{3})\hat R q(t_{2},\: t_{3}) $ $-\dfrac{1}{t_{2}- t_{1}}q^{T}(t_{1},\: t_{2})\hat R q(t_{1},\: t_{2})$,

여기서

$$q(a ,\:b)=\begin{bmatrix}x(b)-x(a)\\ x(b)+x(a)-\dfrac{2}{b-a}\int_{a}^{b}x(s)ds\\ x(b)-x(a)+\dfrac{6}{b-a}\int_{a}^{b}x(s)ds-\dfrac{12}{(b-a)^{2}}\int_{a}^{b}(s-a)x(s)ds\end{bmatrix}.$$

보조정리 2 [7]: 양확정 행렬 $R_{1},\: R_{2}\in R^{q\times q}$, 대칭행렬 $\hat X_{1},\:\hat X_{2}\in R^{q\times q}$, 그리고 사각행렬 $\hat Y_{1},\:\hat Y_{2}\in R^{q\times q}$에 대하여 다음을 만족하면

$$ \left[\begin{array}{rr} R_{1}-\hat{X}_{1} & \hat{Y}_{1} \\ \cdot & R_{2} \end{array}\right] \geq 0,\left[\begin{array}{cc} R_{1} & \hat{Y}_{2} \\ \cdot & R_{2}-\bar{X}_{2} \end{array}\right] \geq 0 $$

스칼라 $\alpha\in(0,\: 1)$에 대하여 항상 다음이 성립한다.

$$-\dfrac{1}{\alpha}\xi_{1}^{T}R_{1}\xi_{1}-\dfrac{1}{1-\alpha}\xi_{2}^{T}R_{2}\xi_{2}\le -\xi_{1}^{T}[R_{1}+(1-\alpha)\hat X_{1}]\xi_{1}$$ $$-\xi_{1}^{T}[\alpha\hat Y_{1}+(1-\alpha)\hat Y_{2}]\xi_{2}-\xi_{2}^{T}[R_{2}+\alpha\hat X_{2}]\xi_{2}$$.

최근에는 이차함수를 LMI로 변환하는 필요충분조건이 제시되어[12-14] 안정성 결과에 크게 기여하였다. 다음의 보조정리 3은 기존의 결과와 다른 형태로 주어지는 필요충분조건이다.

보조정리 3: 스칼라 $y\in R ,\:$ 벡터 $\xi\in R^{q}$ 그리고 대칭행렬 $A_{2},\: A_{1},\: A_{0}\in R^{q\times q}$와 사각행렬 $M,\: Y_{1}\in R^{q\times q}$에 다음의 동치관계가 성립한다.

증명: 먼저 $\mathrm{g}(y)=(y-h_{1})(y - h_{2})\le 0,\:\forall y\in[h_{1},\: h_{2}]$임과 음확정성의 정의, 그리고 잘 알려진 S-procedure 성질[15]을 적용하면 다음을 얻는다.

이것으로 증명을 마친다.

Remark: 보조정리 3의 특징은 $y$항에 있는 행렬의 일반적인 형태인 대칭행렬($A_{1}= A_{1}^{T}$)과 비대칭행령의 대칭행렬 항($Y_{1}^{T}+ Y_{1}$)의 합 형태를 그대로 살려서 구속조건이 있는 2차함수의 음에 관한 필요충분조건을 S-procedure를 이용하여 얻은 것이다. 구하여진 조건이 필요충분조건이 되는 것은 S-procedure의 성질에 의한 것이다.

3. 주요 결과

먼저, 주요 결과를 유도하는데 편의를 위하여, 다음을 정의하자.

(4)
$$ \begin{aligned} &t_{d}=t-d(t), t_{h}=t-h, h_{d}(t)=h-d(t)\\ &\left\{\begin{array}{l} u_{1}(t)=\frac{1}{d(t)} \int_{t_{d}}^{t} x(s) d s, u_{2}(t)=\frac{1}{d^{2}(t)} \int_{t_{d}}^{t}\left(s-t_{d}\right) x(s) d s \\ v_{1}(t)=\frac{1}{h_{d}(t)} \int_{t_{d}}^{t} x(s) d s, v_{2}(t)=\frac{1}{h_{d}^{2}(t)} \int_{t_{h}}^{t_{d}}\left(s-t_{h}\right) x(s) d s \end{array}\right.\\ &\left\{\begin{array}{l} \eta_{0}(t)=\operatorname{col}\left\{x(t), x\left(t_{d}\right), x\left(t_{h}\right)\right\} \\ \eta_{1}(t)=\operatorname{col}\left\{\eta_{0}(t), u_{1}(t), u_{2}(t)\right\} \\ \eta_{2}(t)=\operatorname{col}\left\{\eta_{0}(t), v_{1}(t), v_{2}(t)\right\} \\ \left.\eta_{3}(t)=\operatorname{col}\left\{\eta_{0}(t), d(t) \mid u_{1}(t), u_{2}(t)\right], h_{d}(t)\left[v_{1}(t), v_{2}(t)\right]\right\}, \\ w_{1}(t, s)=\operatorname{col}\left\{\int_{0}^{t} x(\theta) d \theta, x(s), \dot{x}(s)\right\} \\ w_{2}(t, s)=\operatorname{col}\left\{\int_{0}^{t_{d}} x(\theta) d \theta, x(s), \dot{x}(s)\right\} \end{array}\right. \\ &\begin{array}{r} \xi_{t}=\operatorname{col}\left\{x(t), x\left(t_{d}\right), x\left(t_{h}\right), \dot{x}\left(t_{d}\right), \dot{x}\left(t_{h}\right)\right. \\ \left.u_{1}(t), v_{1}(t), u_{2}(t), v_{2}(t)\right\} \end{array} \end{aligned} $$

다음으로 새로운 형태의 LKF 후보함수인 다음 quadratic 함수를 생각하자.

(5)
$$ \begin{aligned} V\left(x_{t}\right)=& d(t) \eta_{1}^{T}(t) S_{1} \eta_{1}(t)+h_{d}(t) \eta_{2}^{T}(t) S_{2} \eta_{2}(t) \\ &+\eta_{3}^{T}(t)\left[d(t) P_{1}+h_{d}(t) P_{2}\right] \eta_{3} \\ &+\int_{t_{d}}^{t} w_{1}^{T}(t, s) Q_{1} w_{1}(t, s) d s+\int_{t_{h}}^{t_{d}} w_{2}^{T}(t, s) Q_{2} w_{2}(t, s) d s \\ &+\int_{t_{h}}^{t}(h-t+s) x^{T}(s) R \dot{x}(s) d s \end{aligned} $$

그러면 대칭행렬들의 조건 $S_{1},\: S_{2,\:}P_{1},\: P_{2,\:}Q_{1},\: Q_{2},\: R>0$하에서(5)의 $V(x_{t})$는 LKF 후보함수로써의 조건을 갖춤을 알 수 있다[16].

다음으로는 위의 LKF 후보함수(5)를 이용하여 시간지연(2)를 갖는 시간지연 선형 시스템(1)의 안정성을 보장하는 주요 결과를 제시한다.

정리 1: 7개의 확정 대칭행렬들 $S_{1},\: S_{2}\in bold R^{5n\times 5n},\:$ $P_{1},\:$ $P_{2}\in bold R^{7n\times 7n},\:$ $ Q_{1},\:$ $Q_{2}\in bold R^{3n\times 3n},\:$ $R\in bold R^{n\times n}$, 2개의 대칭행렬 $\hat X_{1},\:\hat X_{2}\in bold R^{3n\times 3n},\:$ 그리고 4개의 일반행렬 $\hat Y_{1},\:$ $\hat Y_{2}\in bold R^{3n\times 3n},\:$ $M_{1},\: $ $M_{2}\in bold R^{9n\times 9n}$이 존재하여, 다음의 LMI들을 만족하면

(6)
$$ \left[\begin{array}{cc} \bar{R}-\hat{X}_{1} & \hat{Y}_{1} \\ \cdot & \bar{R} \end{array}\right] \geq 0,\left[\begin{array}{cc} \bar{R} & \hat{Y}_{2} \\ \cdot & \bar{R}-\hat{X}_{2} \end{array}\right] \geq 0, M_{i}^{T}+M_{i}>0, i=1,2 $$

(7)
$$ \left[\begin{array}{cc} A_{0}\left(\mu_{1}\right) & \frac{1}{2} A_{1}+Y_{1}^{T}\left(\mu_{1}\right)+h M_{1} \\ \cdot & A_{2}\left(\mu_{1}\right)-\left(M_{1}^{T}+M_{1}\right) \end{array}\right]<0 $$

(8)
$$ \left[\begin{array}{cc} A_{0}\left(\mu_{2}\right) & \frac{1}{2} A_{1}+Y_{1}^{T}\left(\mu_{2}\right)+h M_{1} \\ \cdot & A_{2}\left(\mu_{2}\right)-\left(M_{2}^{T}+M_{2}\right) \end{array}\right]<0 $$

시간지연(2)를 갖는 시간지연 선형 시스템(1)은 점근적으로 안정하다. 여기서

(9)
$$ \begin{array}{l} A_{0}(\dot{d}(t))=\dot{d}(t)\left[C_{1}^{T} S_{1} C_{1}-C_{2}^{T} S_{2} C_{2}\right]+H e\left[C_{1}^{T} S_{1} C_{3}+C_{2}^{T} S_{2} C_{5}\right] \\ -\dot{d}(t) E_{5}^{T}\left(P_{1}-P_{2}\right) E_{5}+h H e\left[E_{5}^{T} P_{2} E_{7}\right]+E_{8}^{T} Q_{1} E_{8} \\ -(1-\dot{d}(t)) E_{9}^{T} Q_{1} E_{9}+H e\left[E_{11} Q_{1} E_{12}\right]+(1-\dot{d}(t)) E_{15}^{T} Q_{2} E_{15} \\ \quad-E_{10}^{T} Q_{2} E_{10}+H e\left[E_{18} Q_{2} E_{19}\right]+h^{2} A_{c}^{T} R A_{c} \\ \quad-U_{1}^{T}\left[\hat{R}+X_{1}\right] U_{1}-H e\left[U_{1}^{T} Y_{2} U_{2}\right]-U_{2}^{T} \hat{R} U_{2} \end{array} $$

(10)
$A_{1}=\dfrac{1}{h}U_{1}^{T}X_{1}U_{1}-\dfrac{1}{h}U_{2}^{T}X_{2}U_{2},\:$

(11)
$$ \begin{aligned} Y_{1}(\dot{d}(t))=& C_{1}^{T} S_{1} C_{4}-C_{2}^{T} S_{2} C_{4}-\dot{d}(t) E_{5}^{T}\left(P_{1}-P_{2}\right) E_{6} \\ &+h E_{0}^{T} P_{2} E_{7}+E_{5}^{T}\left(P_{1}-P_{2}\right) E_{7}-(1-\dot{d}(t)) E_{9}^{T} Q_{1} E_{10} \\ &+E_{11} Q_{1} E_{13}-E_{10}^{T} Q_{2} E_{17}+E_{16}^{T} Q_{2} E_{20} \\ &-\frac{1}{h} U_{1}^{T}\left(Y_{1}-Y_{2}\right) U_{2} \end{aligned} $$

(12)
$$ \begin{aligned} A_{2}(\dot{d}(t))=&-\dot{d}(t) E_{0}^{T}\left(P_{1}-P_{2}\right) E_{\sigma}+H e\left[E_{0}^{T}\left(P_{1}-P_{2}\right) E_{7}\right] \\ &-(1-\dot{d}(t)) E_{10}^{T} Q_{1} E_{10}+H e\left[E_{11}^{T} Q_{1} E_{14}\right] \\ &-E_{17}^{T} Q_{2} E_{17}+H e\left[E_{16}^{T} Q_{2} E_{21}\right] \end{aligned} $$

이며, 또한 사용된 벡터 또는 행렬들은 다음과 같다.

$$ \begin{aligned} &e_{i}=\left[\begin{array}{lll} 0_{n} \times(i-1) n & I_{n} \times n & 0_{n \times(9-i) n} \end{array}\right], i=1,2, \cdots, 9, \quad e_{0}=0_{n \times 9 n}\\ &\tilde{e}_{2}=(1-\dot{d}(t)) e_{2}\\ &\tilde{e}_{6}=(1-d(t)) e_{6}\\ &A_{c}=A e_{1}+A_{1} e_{2}\\ &\hat{R}=\operatorname{diag}\{R, 3 R, 5 R\}\\ &C_{1}=\operatorname{col}\left\{e_{1}, e_{2}, e_{3}, e_{6}, e_{8}\right\}\\ &C_{2}=\operatorname{col}\left\{e_{1}, e_{2}, e_{3}, e_{7}, e_{9}\right\}\\ &C_{3}=\operatorname{col}\left\{e_{0}, e_{0}, e_{0,} e_{1}-\tilde{e}_{2}-d(t) e_{6,} e_{1}-\tilde{e}_{6}-2 d(t) e_{8}\right\}\\ &C_{4}=\operatorname{col}\left\{A_{c,} \tilde{e}_{4}, e_{5}, e_{0}, e_{0}\right\}\\ &C_{5}=\operatorname{col}\left\{h A_{c}, h e_{4}, h e_{5}, \tilde{e}_{2}-e_{3}+\dot{d}(t) e_{7}, \tilde{e}_{2}-e_{7}+2 d(t) e_{9}\right\}\\ &E_{5}=\operatorname{col}\left\{e_{1}, e_{2}, e_{3}, e_{0}, e_{0}, h e_{7}, h e_{9}\right\}\\ &E_{6}=\operatorname{col}\left\{e_{0,} e_{0}, e_{0}, e_{6}, e_{8},-e_{7},-e_{9}\right\}\\ &E_{7}=\operatorname{col}\left\{A_{c}, \tilde{e}_{4}, e_{5}, e_{1}-\tilde{e}_{2}, e_{1}-\tilde{e}_{6}-\dot{d}(t) e_{8}, \tilde{e}_{2}-e_{3}, \tilde{e}_{2}-e_{7}+\dot{d}(t) e_{9}\right\}\\ &E_{8}=\operatorname{col}\left\{e_{0}, e_{1}, A_{c}\right\}\\ &E_{9}=\operatorname{col}\left\{e_{0}, e_{2}, e_{4}\right\} \\ &E_{10}=\operatorname{col}\left\{e_{6}, e_{0}, e_{0}\right\} \\ &E_{11}=\operatorname{col}\left\{e_{1}, e_{0}, e_{0}\right\} \\ &E_{12}=\operatorname{col}\left\{e_{0}, e_{0}, e_{1}-e_{2}\right\} \\ &E_{13}=\operatorname{col}\left\{e_{0}, e_{6}, e_{0}\right\} \\ &E_{14}=\operatorname{col}\left\{e_{8}, e_{0}, e_{0}\right\} \\ &E_{15}=\operatorname{col}\left\{e_{0}, e_{2}, e_{4}\right\} \\ &E_{16}=\operatorname{col}\left\{h e_{7}, e_{3}, e_{5}\right\} \\ &E_{17}=\operatorname{col}\left\{-e_{7}, e_{0}, e_{0}\right\} \\ &E_{18}=\operatorname{col}\left\{\tilde{e}_{2}, e_{0}, e_{0}\right\} \\ &E_{19}=\operatorname{col}\left\{h^{2} e_{9}, h e_{7}, e_{2}-e_{3}\right\} \\ &E_{20}=\operatorname{col}\left\{-2 h e_{9},-e_{7}, e_{0}\right\} \\ &E_{21}=\operatorname{col}\left\{e_{9}, e_{0}, e_{0}\right\} \\ &U_{1}=\operatorname{col}\left\{e_{1}-e_{2}, e_{1}+e_{2}-2 e_{6}, e_{1}-e_{2}+6 e_{6}-12 e_{8}\right\} \\ &U_{2}=\operatorname{col}\left\{e_{2}-e_{3}, e_{2}+e_{3}-2 e_{7}, e_{2}-e_{3}+6 e_{7}-12 e_{9}\right\} \end{aligned} $$

증명: 먼저 위의(5)에 정의된 증강된 LKF 후보함수의 시스템 궤적(1)에 따른 시간 미분을 구하면 다음을 얻는다.

(13)
$$ \begin{aligned} \dot{V}\left(x_{t}\right)=& d(t)\left[\eta_{1}^{T}(t) S_{1} \eta_{1}(t)-\eta_{2}^{T}(t) S_{2} \eta_{2}(t)\right]+2 \eta_{1}^{T}(t) S_{1} d(t) \dot{\eta}_{1}(t) \\ +& 2 \eta_{2}(t) S_{2} h_{d}(t) \dot{\eta}_{2}(t)-\dot{d}(t) \eta_{3}^{T}(t)\left[P_{1}-P_{2}\right] \eta_{3}(t) \\ +& 2 \eta_{3}^{T}(t)\left[h P_{2}+d(t)\left(P_{1}-P_{2}\right)\right] \dot{\eta}_{3}(t)+w_{1}^{T}(t, t) Q_{1} w_{1}(t, t) \\ &-(1-\dot{d}(t)) w_{1}^{T}\left(t, t_{d}\right) Q_{1} w_{1}\left(t, t_{d}\right) \\ &+2 \frac{\partial}{\partial t} w_{1}(t, s) Q_{1} \int_{t_{d}}^{t} w_{1}(t, s) d s \\ &+(1-\dot{d}(t)) w_{2}^{T}\left(t, t_{d}\right) Q_{2} w_{2}\left(t, t_{d}\right) \\ &-w_{2}^{T}\left(t, t_{h}\right) Q_{2} w_{2}\left(t, t_{h}\right)+2 \frac{\partial}{\partial t} w_{2}(t, s) Q_{1} \int_{t_{h}}^{t_{d}} w_{2}(t, s) d s \\ &+h^{2} x^{T}(t) \operatorname{Ri}(t)-h \int_{t_{h}}^{t} x^{T}(s) R \dot{x}(s) d s \\ =& \xi_{t}^{T}\left\{\dot{d}(t)\left[C_{1}^{T} S_{1} C_{1}-C_{2}^{T} S_{2} C_{2}\right]+2 C_{1}^{T} S_{1}\left[C_{3}+d(t) C_{4}\right]\right.\\ &+2 C_{2}^{T} S_{2}\left[C_{5}-d(t) C_{4}\right] \\ &-\dot{d}(t)\left[E_{5}+d(t) E_{0}\right]^{T}\left(P_{1}-P_{2}\right)\left[E_{5}+d(t) E_{0}\right] \\ &+2\left[E_{5}+d(t) E_{0}\right]\left[h P_{2}+d(t)\left(P_{1}-P_{2}\right)\right] E_{7} \\ &+E_{8}^{T} Q_{1} E_{8}-(1-\dot{d}(t))\left[E_{9}+d(t) E_{10}\right]^{T} Q_{1}\left[E_{9}+d(t) E_{10}\right] \\ &+2 E_{11}^{T} Q_{1}\left[E_{12}+d(t) E_{13}+d^{2}(t) E_{14}\right] \\ &+(1-\dot{d}(t)) E_{15}^{T} Q_{2} E_{15} \\ &-\left[E_{10}+d(t) E_{17}\right]^{T} Q_{2}\left[E_{10}+d(t) E_{17}\right] \\ &+2 E_{18}^{T} Q_{2}\left[E_{19}+d(t) E_{20}+d^{2}(t) E_{21}\right] \\ &\left.+h^{2} A_{c}^{T} R A_{c}\right\} \xi_{t}-h \int_{t_{h}}^{t} x^{T}(s) R \dot{x}(s) d s \end{aligned} $$

여기서 벡터 $\xi_{t}$는(4)에 정의된 것이고, 적분항의 상한을 구하기 위하여 보조정리 1과 보조정리 2를 연속적으로 적용하고, $\rho_{1}= U_{1}\xi_{t},\:\rho_{2}= U_{2}\xi_{t}$라 하면 다음을 얻는다.

(14)
$$ \begin{array}{l} -h \int_{t_{h}}^{t} x^{T}(s) R \dot{x}(s) d s \leq-\rho_{1}^{T}\left[\bar{R}+\left(1-\frac{d(t)}{h}\right) \hat{X}_{1}\right] \rho_{1} \\ -2 \rho_{1}^{T}\left[\hat{Y}_{2}+\frac{d(t)}{h}\left(\hat{Y}_{1}-\hat{Y}_{2}\right)\right] \rho_{2}-\rho_{2}^{T}\left[\bar{R}+\frac{d(t)}{h} \widehat{X}_{2}\right] \rho_{2} \end{array} $$

다음으로(14)(13)에 적용하면 다음을 얻는다.

$$ \begin{aligned} \dot{v}\left(x_{t}\right) \leq & \xi_{t}^{T}\left\{d^{2}(t) A_{2}(\dot{d}(t))\right.\\ &\left.+d(t)\left[A_{1}+Y_{1}^{T}(d(t))+Y_{1}(\dot{d}(t))\right]+A_{0}(\dot{d}(t))\right\} \xi_{t} \\ &=\xi_{t}^{T}\{\Psi[d(t), \dot{d}(t)]\}\xi_{t} \end{aligned} $$

여기서 $A_{2}(\dot d(t)),\: A_{1},\: Y_{1}(\dot d(t)),\: A_{0}(\dot d(t))$는(9)-(12)에 정의된 대칭행렬들이고 이들은 모두 스칼라 $\dot d(t)$에 대한 affine 함수이다. 끝으로 보조정리 3을 적용하면 다음을 얻으므로

$$ \begin{aligned} (6)-(8) &\Leftrightarrow\left[\begin{array}{cc} A_{0}\left(\mu_{i}\right) & \frac{1}{2} A_{1}+Y_{1}^{T}\left(\mu_{i}\right)+h M_{i} \\ \cdot & A_{2}\left(\mu_{i}\right)-\left(M_{i}^{T}+M_{i}\right) \end{array}\right]<0, i=1,2 \\ & \Leftrightarrow \Psi[d(t), \dot{d}(t)]<0, \forall d(t) \in[0, h], \mu_{1} \leq \dot{d}(t) \leq \mu_{2}<1 \\ & \Rightarrow \dot{v}\left(x_{t}\right)<0, \forall x_{t} \neq 0 \text { and under }(2). \end{aligned} $$

지연 조건(2)를 갖는 시간지연 시스템(1)은 점근적으로 안정하다. 이것으로 증명을 마친다.

4. 수치 예제

다음은 위에서 제시된 결과의 유용성을 잘 알려진 대표적인 2개의 수치예제를 통하여 보이도록 한다.

예제 1: 다음의 같은 행렬을 갖는 시간지연 시스템(1)을 고려하자[2-4,7-10,12-14].

$$ A=\begin{bmatrix}-2&&0\\0&&-0.9\end{bmatrix},\:A_{1}=\begin{bmatrix}-1&&0\\-1&&-1\end{bmatrix}$$

다음의 표 1은 $-\mu_{1}=\mu_{2}$의 여러 값에 따른 예제 1의 시스템의 안정도를 보장하는 최대 허용 $h$를 기존의 결과들과 비교 정리한 표이다.

표 1. 예제 1 시스템의 안정도를 보장하는 최대 시간지연

Table 1. Maximal allowable delay for Example 1

$-\mu_{1}=\mu_{2}$

0.1

0.5

0.8

[2]

4.7038

2.4208

2.1377

[3]

4.8117

3.1282

2.6982

[4]

4.7883

3.0557

2.6150

[8]

4.809

3.109

2.710

[10]

5.01

3.19

2.70

[12]

5.044

3.443

2.983

[14]

5.026

3.428

2.977

[13]

5.084

3.482

3.005

본 논문

5.094

3.499

3.049

표 1에서 보듯이 예제 1의 경우에 새로이 제시된 정리1의 결과는 기존의 결과보다 우수함을 알 수 있다.

다음 표 2표 1의 계산을 위해 필요한 변수의 개수를 최근의 결과들과 비교한 것이다.

표 2. 최근 결과에 사용된 변수 수

Table 2. Number of variables for recent results

[6]

[14]

[12]

[13]

본 논문

$27 n^{2}+ 4n$

\begin{align*} 252.5 n^{2}\\+ 0.5 n \end{align*}

$235n^{2}+ 34n$

\begin{align*} 249.5n^{2}\\+ 13.5n \end{align*}

$272 n^{2}+19n$

예제 2: 다음의 같은 행렬을 갖는 시간지연 시스템(1)을 고려하자[2,6,12-14].

$$A=\begin{bmatrix}0&&1\\-1&&-2\end{bmatrix},\:A_{1}=\begin{bmatrix}0&&0\\-1&&1\end{bmatrix}$$

다음의 표 3은 예제 1과 마찬가지로 여러 경우의 $-\mu_{1}=\mu_{2}$ 값에 따른 예제 2의 안정도를 보장하는 최대 허용 $h$를 기존의 결과와 비교 정리한 표이다.

표 3. 예제 2 시스템의 안정도를 보장하는 최대 시간지연

Table 3. Maximal allowable delay for Example 2

$-\mu_{1}=\mu_{2}$

0.1

0.2

0.5

0.8

[2]

6.590

3.672

1.411

1.275

[6]

6.727

3.920

1.923

1.367

12]

7.685

4.969

2.774

2.117

[14]

7.651

4.936

2.764

2.114

[13]

7.714

5.003

2.809

2.149

본 논문

7.723

5.024

2.838

2.165

표 3에서 보듯이 예제 2의 경우도 새로이 제시된 정리1의 결과는 기존의 결과보다 우수함을 알 수 있다.

5. 결 론

본 논문에서는 새로운 형태의 증강된 LKF 후보함수를 제안하고, 먼저 Bessel-Legendre 부등식, 개선된 reciprocally convex 부등식을 적용하여 이의 시간 미분의 상한을 구하였다. 다음으로 2차 행렬 다항식이 정하여진 구간에서 음이 되는 새로운 형태의 필요충분조건을 이용하여 안정성을 보장하는 결과를 LMI로 제시하였다. 끝으로 제시된 결과는 두 개의 대표적 예제를 통하여 기존의 결과보다 우수함을 보였다.

References

1 
P. Park, J. W. Ko, C. Jeong, 2011, Reciprocally convex approach to stability of systems with time-varying delays, Automatica, Vol. 47, pp. 235-238DOI
2 
A. Seuret, F. Gouaisbaut, 2013, Wirtinger-based integral inequality: Application to time-delayed systems, Automatica, Vol. 49, pp. 2860-2866DOI
3 
O. M. Kwon, M. J. Park, J. H. Park, S. M. Lee, E. J. Cha, 2014, Improved results on stability of linear systems with time-varying delays via Wirtinger-based integral inequality, J. of the Franklin Institute, Vol. 351, pp. 5382-5398DOI
4 
H. -B. Zeng, Y. He, M. Wu, J. She, 2015, Free matrix based integral inequality for stability analysis of systems with time-varying delay, IEEE Trans. on Automatic Control, Vol. 60, pp. 2768-2772DOI
5 
A. Seuret, F. Gouaisbaut, 2015, Hierarchy of LMI conditions for the stability analysis of time-delay systems, Systems & Control Letters, Vol. 81, pp. 1-7DOI
6 
J. -H. Kim, 2016, Further improvement of Jensen inequality and application to stability of time-delayed systems, Automatica, Vol. 64, pp. 121-125DOI
7 
A. Seuret, F. Gouaisbaut, 2016, Delay dependent reciprocally convex combination lemma, Rapport LAAS n16006Google Search
8 
C. -K. Zhang, Y. He, L. Jiang, M. Wu, 2017, Notes on stability of time-delay systems: bounding inequalities and augmented Lyapunov-Krasovskii functionals, IEEE, Trans. Automatic Control, Vol. 62, pp. 5331-5336DOI
9 
X. -M. Zhang, Q. -L. Han, A. Seuret, F. Gouaisbaut, 2017, An improved reciprocally convex inequality and an augmented Lyapunov-Krasovskii functional for stability of linear systems with time-varying delay, Automatica, Vol. 84, pp. 221-226DOI
10 
A. Seuret, F. Gouaisbaut, 2018, Stability of linear systems with time-varying delays using Bessel-Legendre inequalities, IEEE Trans. Automatic Control, Vol. 63, pp. 225-232DOI
11 
X. -M. Zhang, Q. -L. Han, A. Seuret, F. Gouaisbaut, Y. He, 2019, Overview of recent advances in stability of linear systems with time-varying delays, IET Control Theory and Applications, Vol. 12, pp. 1-16DOI
12 
F. de Oliveira, M. C., F. Souza, 2020, Further refinements in stability conditions for time-varying delay systems, Applied Math. Comput., Vol. 359DOI
13 
B. L. Zhang, L. Cheng, K. Pan, X -M, Zhang, 2020, Reducing conservatism of stability criteria for linear systems with time-varying delay using an improved tripple-integral ine- quality, Applied mathematics and computation, Vol. 380, pp. 125-254DOI
14 
J. M. Park, P. G. Park, 2020, Finite-interval quadratic polynomial inequalities and their application to time-delay systems, J. of Franklin Institute, Vol. 357, pp. 4316-4327DOI
15 
S. Boyd, L. E. Ghaoui, E. Feron, V. Balakrishhnan, 1994, Linear Matrix Inequalities in System and Control Theory, Studies in Applied mathematicsGoogle Search
16 
K. Gu, V.L. Kharitonov, J. Chen, 2003, Stability of time-delay systems, BirkhausserGoogle Search
17 
E. Fridman, 2014, Introduction to time-delay systems: Analysis and control, BirkhauserGoogle Search

저자소개

Jin-Hoon Kim
../../Resources/kiee/KIEE.2020.69.12.1904/au1.png

He received the B.S. from Seoul National University in 1985, and M.Sc and Ph.D. from KAIST in 1989 and 1983, respectively.

From 1993 to 1994, he was a full-time lecturer with the department of control and instrument engineering, Gyeongsang National University.

In 1998 and 2008, he was a visiting professor at UCI and UTA, respectively.

Since 1995, he is currently a professor with the school of electronics engineering, Chungbuk National University, and a researcher at the Research Institute for Computer and Information Communication.

Tel: 043-261-2387, Fax: 043-268-2386.

E-mail : jinhkim@cbnu.ac.kr