The Korean Society of Marine Engineering
[ Original Paper ]
Journal of the Korean Society of Marine Engineering - Vol. 42, No. 7, pp.566-571
ISSN: 2234-7925 (Print) 2234-8352 (Online)
Print publication date 30 Sep 2018
Received 06 Nov 2017 Revised 25 May 2018 Accepted 18 Sep 2018
DOI: https://doi.org/10.5916/jkosme.2018.42.7.566

특이 섭동 시간지연 시스템의 안정성 분석

김범수1 ; 권재성
Stability analysis for singularly perturbed time delayed system
Beomsoo Kim1 ; Jaesung Kwon
1Associate Professor, Department of Mechanical System Engineering, Gyeongsang National University, Tel: 055-772-9101 kimbs@gnu.ac.kr

Correspondence to: Assistant Professor, Department of Mechanical System Engineering, Institute of Marine Industry, Gyeongsang National University, 38, Cheondaegukchi-gil, Tongyeong, Gyeongnam, 53064, Korea, E-mail: jkwon@gnu.ac.kr, Tel: 055-772-9102

Copyright ⓒ The Korean Society of Marine Engineering
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0), which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

초록

본 연구에서는 시간 지연을 갖는 특이 섭동 시스템의 안정성을 다루기 위해서 리아프노프-크라보스키 범함수에 기 반한 선형행렬 부등식 기법을 통해서 안정한 조건을 제시한다. 먼저 시간 지연이 있는 특이섭동 시스템 전체에 관해 안 정 조건을 구한다. 그리고 특이 섭동 시스템을 Chang 변환 기법을 이용해서 저차의 빠른 부시스템과 느린 부시스템으로 완전히 분리한 후에 각 부시스템에서의 안정 조건을 제시함으로써 적은 변수로 선형 행렬 부등식 문제를 효율적으로 풀 수 있음을 보인다. 마지막으로 예제 시스템을 통해 제안한 방법의 유효성과 장점을 보인다.

Abstract

In this paper, we propose the stability conditions for a singularly perturbed time-delayed system using Lyapunov–Krasovskii functionals and a linear matrix inequality (LMI). Further, we show that the singularly perturbed time-delayed system can be completely decoupled into a fast time-delayed subsystem and a slow time-delayed subsystem using the Chang transform and that the stability conditions for each subsystem can be efficiently derived using an LMI. This method reduces the variables in the LMI. Finally, the effectiveness and merits of the proposed approach are demonstrated through a numerical example.

Keywords:

Time delay, Singularly perturbed system, Stability, Lyapunov-Krasovskii functionals, Linear matrix inequality

키워드:

시간지연, 특이섭동시스템, 안정성, 리아프노프-크라보스키 범함수, 선형 행렬 부등식

1. 서 론

공학, 생물학 등 많은 분야의 시스템이 시간 지연을 갖는 미분방정식으로 표현된다[1]-[5]. 시간 지연성분이 지연 미분 방정식(DDE, Delay Differential Equation)으로 모델링 될 때 초월 특성(transcendental characteristics)에 의해 DDE는 무수히 많은 고유치를 가질 뿐만 아니라 시간 지연이 안정도에 미치는 영향을 고려해야 하며 이는 시스템 해석을 어렵게 한다.

특이 섭동 시스템은 기계적 요소와 전기적 요소가 결합된 시스템이나 매우 작은 시상수로 인한 와류(parasitic) 파라미터 등에 의해 차수가 증가된 시스템으로 느린 응답과 빠른 응답이 동시에 존재하는 다중 시간 척도(multi time scale) 특성을 갖는다. 특이 섭동 시스템은 조건수(condition number)가 큰 불량 조건(ill-condition)을 가져서 시스템 해석을 어렵게 한다. 따라서 이 문제를 회피하기 위해서 Chang 변환이나, 전체 시스템을 느린 부 시스템(slow subsystem)과 빠른(fast) 부 시스템으로 분리한 후 시스템 해석과 제어기를 설계하는 연구들이 진행되고 있다[1][2].

선박에서의 방향타 안정화 설계[6]와 전몰형 수중익선(fully submerged hydrofoil)의 회두각 유지(course keeping) 제어[7] 문제를 특이섭동의 시분리(time scale separation)를 이용해서 축차된 느린 시스템과 빠른 시스템으로 분리한 후 합성제어 기법으로 안정도를 해결하는 연구가 보고되었다.

시간 지연을 갖는 특이섭동 시스템의 안정도 문제를 해결하기 위해 Lyapunov-Krasovskii 범함수를 이용하는 기법 등이 연구되어 왔다[8]. S. Pan et al. [9]은 다중 시간 지연을 갖는 동차 특이 섭동시스템을 느린 부 시스템과 빠른 부 시스템으로 분리한 후 각 시스템에서의 안정도 문제를 다루었으며, Z. Shao et al.은[10] 시간 지연을 갖는 동차 특이 섭동시스템을 Lyapunov 방법을 이용해서 안정도 문제를 다루었다. B. S. Kim et al. [11]은 합성제어 기법으로 시스템을 분리한 후 놈(norm)을 이용하여 안정화 영역을 구한 후 제어입력을 구하고 Lambert W 함수를 이용해서 시스템을 해석하는 방법을 제시하였다.

본 연구에서는 시간 지연을 갖는 특이 섭동 시스템을 Lyapunov-Krasovskii 정리와 선형 행렬 부등식을 이용하여 안정할 조건을 구한다. 선형 행렬 부등식 문제에 있어서 주요 관점 중 하나는 가능한 적은 변수로 문제를 해결하고자 하는 것이다. 따라서 본 연구에서는 특이 섭동 시스템을 Chang 변환[12]을 이용하여 저차의 부 시스템으로 완전 분리한 후, 각 부 시스템에서 LMI(Linear Matrix Inequality)를 이용해서 안정할 조건을 구한 후 이를 합성하여 전체 안정할 제어 입력을 구한다. 제안된 방법은 LMI를 통해서 전체 시스템의 안정 조건을 구할 때 보다 적은 변수를 이용해서 구할 수 있다.


2. 시간지연이 있는 특이 섭동 시스템

2.1 순수 시간지연을 갖는 특이 섭동 시스템

다음과 같은 시간 지연을 갖는 특이 섭동 시스템을 고려하자.

Eeẋt=Axt-h(1) 

여기서 Ee=In100ϵIn2,  A=A11A12A21A22,  xt=x1x2이며 x1tRn1, x2tRn2, n=n1+n2는 각각 느린 부 시스템과 빠른 부 시스템을 나타내고, 0<1 이고, h>0 는 시간지연 상수이다. 그리고 A11Rn1×n1,A12Rn1×n2,A21Rn2×n1,A22Rn2×n2 이다.

보조 정리 1.[13] 알맞은 차원을 갖는 행렬 Zi=ZiT i=1,2,3,4,5다음 선형 행렬 부등식 (2)를 만족하면

Z1>0Z1+ϵZ3ϵZ5TϵZ5ϵZ2>0Z1+ϵZ3ϵZ5TϵZ5ϵZ2+ϵ2Z4>0(2) 

선형 행렬 부등식 (3)이 성립된다.

EeZe=ZeEe>0(3) 

여기서 Ze=Z1+ϵZ3ϵZ5TϵZ5Z2+ϵZ4 이다.

보조 정리 1과 LKF(Lyapunov-Krasovskii functionals)을 이용해서 식 (1) 의 안정 조건을 구하면 다음과 같다.

정리 1. 대칭 양의 정부호 행렬(Positive Definite Matrix) Q가 존재해서 식 (2)선형 행렬 부등식 (4)를 만족하면 식 (1) 은 점근적으로 안정하다.

*100-Q+hATEe-1A<0*1=ATZe+ZeA+hZeAEe-1ATZe+Q(4) 

증명: 다음과 같이 LKF 을 선정하자.

Vxt=V1xt+V2xt+V3xt(5) 

여기서

V1xt=xTtEeZext(6) 
V2xt=t-htxTsQxsds(7) 
V3xt=-h0t+0tẋTsEeẋsdsdθ  (8) 

이다. Q가 양한정(positive definite) 행렬이고 보조정리 1 에 의해서 EeZe>0 이므로 Vxt>0 이 성립한다. 식(6)을 미분하면 다음과 같다.

V˙1xt=x˙TtEeTZext+xTtEeZex˙t=xTt-hATZext+xTtZeAxt-h(9) 

아래와 같은 Newton-Leibniz 모델 변환 xt-h=xt-t-htẋsds 을 이용하면 식 (9)는 다음과 같이 다시 쓸 수 있다.

V̇1xt=xt-t-htẋsdsTATZext+xTtZeAxt-t-htẋsds=xTtATZext+xTtZeAxt-2xTtZeAt-htẋsds(10) 

두 번째 항인 식 (7)을 미분하면 다음과 같다.

V̇2xt=ddtt-htxTsQxsds=dtdtxTtQxt-dt-hdtxTt-hQxt-h=xTtQxt-xTt-hQxt-h

마지막 항인 식 (8)을 미분하면 다음과 같다.

V̇3xt=hẋTtEetẋt-t-htẋTsEeẋsds=hxTt-hATEe-1Axt-h-t-htẋTsEeẋsdθ

식 (10)의 적분 항은 아래 부등식 형태로 유도된다.

-2xTtZeAt-htẋsdshxTtZeAEe-1ATZeTxt+t-htẋTsEeẋsds

식 (5)로 주어진 LKF의 미분은 식 (6) - 식 (8) 항을 미분한 것을 더한 결과와 같으므로 식 (6)의 적분항을 식 (10)으로 대치하면 LKF 의 미분 결과는

V̇xtxTtATZext+xTtZeAxt+hxTtZeAEe-1ATZext+xTtQxt+xTt-hhATEe-1A-Qxt-h

로 된다. 여기서 yt=xTt xTt-hT 라 놓으면 V̇xt 의 우변은 다음과 같이 행렬로 표현할 수 있다.

yTt1,100hATEe-1A-Qyt1,1=ATZe+ZeA+hZeAEe-1ATZe+Q

여기서 V̇xt의 우변이 0 보다 작으면 시간 지연을 갖는 특이 섭동 시스템 식 (1)은 점근적으로 안정하다.

정리 1의 선형행렬부등식 문제에서 Zi, (i = 1,..,5) 와 Q가 대칭행렬이므로 3.5n2개의 변수를 구하면 된다.

2.2 시간지연을 갖는 특이 섭동 시스템의 분리

다음과 같은 시간지연 특이 섭동 시스템을 고려하자.

ẋt=Aext-h+But(11) 

여기서 B=B11ϵB2, Ae=A11A121ϵA211ϵA22, utRm 이다.

작은 파라미터 ϵ이 포함된 시스템 행렬 Aϵ 는 불량조건(ill-condition)으로서 시스템 해석을 어렵게 하고 있다 따라서 이 불량조건 문제를 회피하기 위해서 식 (12)로 구성된 Chang 변환 행렬을 이용하여 특이 섭동 시스템 식 (11)ϵ와 결합된 저차의 상태 방정식과 결합되지 않은 저차화된 상태 방정식으로 완전 분리한다.

T=I-ϵHL-ϵHLI.(12) 

식 (12)에서 행렬 L,H식 (11)의 시스템 행렬 Aϵiji,j = 1,2로 이루어진 비대칭 Riccati 방정식 (13) 및 (14) 의해이다.

Aϵ21-Aϵ22L+ϵLAϵ11-ϵLAϵ12L=0(13) 
Aϵ12-ϵAϵ11-Aϵ12LH-HAϵ22-ϵLAϵ12=0(14) 

식 (13)식 (14) 에서 LH 를 구하기 위해서는 오차가 주어진 값 보다 작아질 때 까지 뉴턴 반복법으로 구한다[14]. Chang 변환의 역행렬은 완전 분리된 시스템을 원래의 시스템으로 변환하기 위해 사용되는 것으로 다음과 같이 해석적으로 구할 수 있다.

T-1=IϵH-LI-ϵLH(15) 

Chang 변환을 통한 새로운 상태를 다음과 같이 정의하자.

η1tη2t=Tx1tx2t.(16) 

식 (16)을 이용하여 시스템 식 (11)를 Chang 변환하면 다음과 같다[2].

η̇1tϵη̇2t=Aϵ11-Aϵ12L00ϵLAϵ12+Aϵ22η1t-hη2t-h+B1-ϵHLB1-HB2ϵLB1+B2ut(17) 

따라서 시스템 식 (11)은 다음과 같이 저차의 느린 시스템 식 (18)과 저차의 빠른 시스템 식 (19)로 완전 분리된다.

η̇1t=Aϵ11-Aϵ12Lη1t-h+B1-ϵHLB1-HB2ut.(18) 
ϵη̇2t=ϵLAϵ12+Aϵ22η2t-h+ϵLB1+B2ut.(19) 

2.3 느린 부시스템에서의 안정도

식 (18)로 표현된 느린 부시스템에 상태 지연이 있는 제어 입력 식 (20)을 인가하면 느린 상태로만 구성된 시스템 식 (21)을 얻는다.

u1t=K1η1t-h(20) 
η̇1t=Aη1η1t-hAη1=B1-ϵHLB1-HB2K1+Aϵ11-Aϵ12L(21) 

정리 1과 유사하게 시스템 식 (21)에 대한 점진적 안정도 조건을 Newton-Leibniz 모델 변환을 이용하면 다음과 같이 구할 수 있다.

보조 정리 2. Q=QT>0, P=PT>0, S=ST>0 인 행렬이 존재해서 선형 행렬 부등식 (22)을 만족하면 시스템 식 (21)은 점진적으로 안정하다.

*200hAη1TSAη1-Q<0*2=Aη1TP+PAη1+hPAη1SAη1TP+Q(22) 

증명: 다음과 같이 LKF 를 선정하자.

V1η1=η1TtPη1t(23) 
V2η1=t-htη1TsQη1sds(24) 
V3η1=-h0t+0tṅ1TsSη̇1sdsdθ(25) 

식 (23)을 미분하면 다음과 같다.

V̇1η1=η1Tt-hATPη1t+η1TtPAη1t-h

이다. 여기에 Newton-Leibniz 변환 행렬

η1t-h=η1t-t-htη̇1sds

을 이용하면 V̇η1은 다음과 같이 다시 쓸 수 있다.

V̇1η1=η1TtAη1TPη1t+η1TtPAη1η1t-2η1TtPAη1t-htη̇1sds

여기서

t-htη1tη1a-hTPAη1Aη1TSPAη1Aη1TTη1tη1a-hda0

를 이용하면

-2η1TtPAη1t-htṅ1sdshη1TtPAη1SAη1Tη1t+t-htη̇1TsSη̇1sds

이다. 식 (24)식 (25)의 미분은 각각 다음과 같다.

V̇2η1=η1TtQη1t-η1Tt-hQη1t-hV̇3η1=hη1Tt-hATSAη1t-h-t-htẋTsSẋsds

따라서 LKF 의 미분을 더한 결과는 다음과 같다.

V̇1η1+V̇2η1+V̇3η1η1TtAη1TPη1t+η1TtPAη1η1t+η1TtQη1t-η1Tt-hQη1t-h+hη1Tt-hATSAη1t-h+hη1TtPAη1SAη1Tη1t

따라서 식 (22)가 성립하면 시간 지연 시스템 식 (21)은 점진적으로 안정하다.

2.4 빠른 부시스템에서의 안정도

식 (19)로 표현된 빠른 부시스템에 제어입력 식 (26)을 인가하면 빠른 상태로만 구성된 시스템 식 (27)을 구축할 수 있다.

u2t=K2η2t-h(26) 
ϵη̇2tϵLAϵ12+Aϵ22+ϵLB1+B2K2η2t-h(27) 

느린 부시스템에서와 같이 빠른 부시스템에서도 다음 보조정리 3 을 만족하면 시스템 식 (27)은 점진적으로 안정하다.

보조 정리 3. Q=QT>0, P=PT>0 인 행렬이 존재해서 선형 행렬 부등식 (28)을 만족하면 시스템 식 (27)은 점진적으로 안정하다.

*300hϵAη2TAη2-ϵQ<0*3=Aη2TP+PAη2+hϵPAη2Aη2TP+ϵQ(28) 

여기서 Aη2=ϵLAϵ12+Aϵ22+(ϵLB1+B2)K2 이다.

증명: 여기서는 Lyapunov-Krasovskii 범함수를 다음과 같이 선정한다.

V1η2=η2TtPη2t(29) 
V2η2=t-htη2TsQη2sds(30) 
V3η2=-h0t+θtη̇2TsSη̇2sdsdθ(31) 

식 (31)에서 S=I 로 놓은 후 보조정리 2의 증명 방법을 적용하면 부등식 (28)을 구할 수 있다.

2.5 합성 제어 시스템에서의 안정도

식 (20), 식 (26), 그리고 식 (16)에 의해서 합성 제어 u(t)는 다음과 같다.

ut=u1t+u2t=K1K2Tx1t-hx2t-h(32) 

각 부시스템에서 구한 제어 입력과 전체 시스템에 인가할 제어 입력과의 관계는 다음과 같이 정리된다.

보조 정리 4. 보조정리 2 와 3을 만족하면 특이섭동 시스템 식 (11)은 점근적으로 안정하다.

증명: 각 부시스템 식 (18)식 (19) 가 안정하면 Chang 변환 행렬 T가 상수 행렬이면서 비특이행렬(nonsingular matrix) 이므로 x1Ttx2TtT=T-1η1Ttη2tTT 이다. η1Ttη2tT<α, α>0 이면 T-1η1Ttη2tTTT-1 η1Ttη2tTT 이므로 x1Ttx2TtT α 이 성립하여 전체 시스템은 안정하다.

시간 지연이 포함된 특이섭동 시스템의 안정조건을 구한 정리 1 의 경우에는 LMI를 풀 때 3.5n2개의 변수를 구해야 하지만 두 개의 부 시스템으로 분리된 접근법에서 보조정리 2 및 3을 이용할 때는 3n2개의 변수만을 구하면 된다.


3. 수치 예

다음과 같이 시간 지연을 갖는 특이 섭동 시스템을 고려하자.

Eeẋt=Axt-h+But(33) 

여기서

Ee=I20202ϵI2, A=5020900.524100.3450.4652000.262-5, B=1.20.010.10.0226 

이고, h=0.1, ϵ=0.09 이다.

식 (15)의 Chang 변환 행렬은 비대칭 Riccati 방정식 (13)(14)를 풀면

L=3.75451.49091.49060.3749,  H=0.00750.00030.05360.0028

로 구해지고 식 (21)로 표현된 느린 부시스템의 파라미터는 다음과 같이 구해진다.

Aη1a=Aϵ11-Aϵ12L=5.09.0-1.2953-0.5144Aη1b=B1-ϵHLB1-HB2=1.1961-0.0177

부등식 (22)을 슈어 보수를 이용하면

*000-QAη1T0Aη1-1hS-1<0*=Aη1TP+PAη1+hPAη1SAη1TP+Q(34) 

와 같이 되고, Z1=P-1로 정의한 후 부등식 (34) 양변에 diag{Z1, Z1} 를 곱하면 다음과 같다.

**000-Z1QZ1Z1Aη1a+Aη1bK1Z10Z1Aη1a+Aη1bK1TZ1-1hZ1S-1Z1<0**=Z1Aη1a+Aη1bK1T+Aη1a+Aη1bK1Z1+hPAη1a+Aη1bK1SAη1a+Aη1bK1T+Z1QZ1(35) 

여기서 F1=K1Z1, L1=Z1QZ1, J1=K1S, M1=K1SK1T, N1=Z1S-1Z1, V1=Z1Aη1bF1, W1=Z1Aη1aZ1와 같이 정의하여 변수 변경하면 부등식 (35)

000-L1W1T+V1T0W1+V1-1hN1<0=Z1Aη1aT+Aη1aZ1+hSAη1aT+hJ1TAη1bT+Aη1bhJ1Aη1aT+F1+hM1TAη1bT+L1+F1T+Aη1bT(36) 

와 같이 LMI 형태로 변환된다. 선형 행렬 부등식 (36)에서 F1, Z1을 구한 후 식 (20)으로 표현된 느린 부시스템의 제어입력 이득은 다음과 같이 구해진다.

K1=F1Z1-1=-8.2433-5.7770(37) 

마찬가지로 위에서 구해진 L,H를 적용하면 빠른 부시스템 식 (27)의 파라미터는 다음과 같이 구해진다.

Aη2a=ϵLAϵ12-Aϵ22=5.21302.911122.2339-55.5556Aη2b=ϵLB1+B2=0.50680.1839

합동 변환(Z2=P-1)과 F2=K2Z2, L2=ϵZ2QZ2, M2=K2K2T,N2=Z2Z2, V2=Z2Aη2bF2, W2=Z2Aη2aZ2 와 같이 변수 변경하면 부등식 (28)을 슈어 보수한 행렬은

 Δ000-L2W2T+V2T0W2+V2-ϵhN2<0Δ=Z2Aη2aT+F2TAη2aT+Aη2aZ2+Aη2bF2+hAη2aAη2aT+hAη2bK2Aη2aT+hAη2aK2TAη2bT+hAη2bM2TAη2bT+L2(38) 

와 같이 LMI 형태로 변환된다. 선형행렬 부등식 (38)에서 F2, Z2를 구한 후 식 (26)으로 표현된 느린 부시스템의 제어입력 이득은 다음과 같이 구해진다.

K2=F2Z2-1=-68.8366-64.1008(39) 

식 (32)로 표현된 전체 시스템의 제어입력은 식 (37)식 (39)에서 구해진 부시스템들의 제어입력을 이용하면 K=[-277.62 9.55 128.01 -511.99]와 같이 구해진다.


4. 결 론

본 논문에서는 시간지연을 갖는 특이 섭동 시스템의 안정화 조건을 전체 시스템을 고려한 경우와 Chang 변환을 통해 시스템을 저차의 느린 부시스템과 빠른 부시스템으로 완전 분리한 후 각 부시스템에서 안정화 조건을 제시했다. 시간지연을 고려하여 Lyapunov-Krasovskii 범함수에 이차 적분을 포함한 덜 보수적인 방법으로 안정화 조건을 보였다. 시스템을 저차의 시스템으로 분리하고 안정화 조건을 선형 행렬 부등식을 이용해서 구할 때 변수의 개수를 줄일 수 있음을 보였고 수치 예를 통해 그 타당성을 보였다.

Acknowledgments

이 연구는 2016년도 경상대학교 신임교원 연구기반 조성 연구비에 의해 수행되었음.

References

  • P. Kokotovic, H. Khali, and J. O'Reilly, Singular Perturbation Methods in Control: Analysis and Design, Society for Industrial and Applied Mathematics, (1987).
  • B. S. Kim, “Stabilizing controller design for time-delay singularly perturbed systems by H norm and lambert W function”, The Transactions of the Korean Institute of Electrical Engineers, 62(8), p1144-1150, (2014), (in Korean). [https://doi.org/10.5370/kiee.2013.62.8.1144]
  • H. S. Seong, and H. S. Choi, “Robust control of a heave compensation system for offshore cranes considering the time-delay”, Journal of the Korean Society of Marine Engineering, 41(1), p105-110, (2017), (in Korean). [https://doi.org/10.5916/jkosme.2017.41.1.105]
  • S. L. Lee, Y. H. Lee, J. K. Ahn, J. K. Son, K. T. Ryu, and M. O. So, “PID control of unstable processes with time delay”, Journal of the Korean Society of Marine Engineering, 33(5), p721-728, (2009), (in Korean).
  • J. Richard, “Time-delay systems: An overview of some recent advances and open problems”, Automatica, 39(10), p1667-1694, (2003). [https://doi.org/10.1016/s0005-1098(03)00167-5]
  • R. Ren, Z. Zou, and X. Wang, “A two-time scale control law based on singular perturbations used in rudder roll stabilization of ships”, Ocean Engineering, 88, p488-498, (2014). [https://doi.org/10.1016/j.oceaneng.2014.07.006]
  • S. Liu, C. Xu, and L. Zhang, “Robust course keeping control of a fully submerged hydrofoil vessel without velocity measurement: An iterative learning approach,”, Mathematical Problems in Engineering, (2017), [Online]. Available: https://www.hindawi.com/journals/mpe/2017/7979438. [https://doi.org/10.1155/2017/7979438]
  • P. G. Park, W. L. Lee, and S. Y. Lee, “Stability on time delay systems: A survey”, Journal of Institute of Control, Robotics and Systems, 20(3), p289-297, (2014), (in Korean).
  • S. Pan, C. Chen, and J. Hsieh, “Stability analysis for a class of singularly perturbed systems with multiple time delays”, Journal of Dynamic Systems Measurement and Control, 126(3), p462-466, (2004). [https://doi.org/10.1115/1.1793172]
  • Z. Shao, and J. Rowland, “Stability of time-delay singularly perturbed systems”, IEE Proceedings - Control Theory and Applications, 142(2), p111-113, (1995). [https://doi.org/10.1049/ip-cta:19951739]
  • B. S. Kim, and S. W. Ahn, “Stabilizing controller design for time-delay pure singularly perturbed systems via lambert W function”, Journal of the Korean Society for Power System Engineering, 18(1), p120-127, (2014), (in Korean). [https://doi.org/10.9726/kspse.2014.18.1.120]
  • K. W. Chang, “Singular perturbations of a general boundary value problem”, “Singular perturbations of a general boundary value problem,”, 3(3), p520-526, (1972). [https://doi.org/10.1137/0503050]
  • F. Sun, L. Zhou, Q. Zhang, and Y. Shen, “Stability bound analysis and synthesis for singularly perturbed systems with time-varying delay”, Mathematical Problems in Engineering, (2013), [Online]. Available: https://www.hindawi.com/journals/mpe/2013/517258. [https://doi.org/10.1155/2013/517258]
  • Z. Gajic, and M. T. Lim, Optimal Control of Singularly Perturbed Linear Systems and Applications : High-Accuracy Techniques, New York: Taylor & Francis Inc, (2001).