Low order controller design for systems with time delays

Size: px
Start display at page:

Download "Low order controller design for systems with time delays"

Transcription

1 20 50th IEEE Conference on Deciion and Contro and European Contro Conference (CDC-ECC) Orando, FL, USA, December 2-5, 20 Low order controer deign for ytem with time deay A N Gündeş and H Özbay Abtract Finite-dimeniona controer ynthei method are deveoped for ome cae of inear, time-invariant, ingeinput inge-output, or muti-input muti-output ytem, which are ubject to time deay The propoed ynthei procedure give ow-order tabiizing controer that ao achieve integraaction o that contant reference input are tracked aymptoticay with zero teady-tate error I INTRODUCTION A wide range of dynamica phenomena cannot be modeed ufficienty accuratey a finite-dimeniona inear timeinvariant (LTI) ytem due to time deay The effect of thee deay often cannot be ignored and have to be incuded in the mode 4, 7 Thi paper preent finite-dimeniona tabiizing controer ynthei method for ome cae of LTI, inge-input inge-output (SISO) or muti-input mutioutput (MIMO) ytem that are ubject to time deay The propoed controer are impe and have ow-order, and they ao provide integra-action o that tep-input reference are tracked aymptoticay with zero teady-tate error The pant cae conidered in Section III-A and Section III-B have no retriction on their poe Thee pant may be tabe or untabe The (tranmiion) zero in the open eft-haf compex pane (OLHP) are unretricted and there may be any number of zero at infinity a we The dua cae in Section III-C conider pant with no retriction on the number or ocation of the (tranmiion) zero, but the poe are either in the OLHP or at the origin = 0 Section III-A examine SISO deayed pant of retarded type (eg, 4,, 2) and Theorem deveop a impe controer ynthei procedure, which i generaized and extended to MIMO ytem in Theorem 2, Section III-B Stabiity of deay ytem of retarded type and of neutra type wa tudied extenivey and many deay-independent and deay-dependent tabiity reut are avaiabe 7, 0 The tuning and interna mode contro technique ued in proce contro ytem generay appy to deay ytem 2 Infinite dimeniona integra action controer have been deigned in to maximize the aowabe controer gain uing the robut contro technique for infinite dimeniona ytem 5 For MIMO tabe pant ubject to inputoutput deay, proportiona-derivative (PD) and proportionaintegra-derivative (PID) controer were deigned in 8 for pant that have no more than two untabe poe coe Department of Eectrica and Computer Engineering, Univerity of Caifornia, Davi, CA 9566 angunde@ucdaviedu H Özbay i with the Department of Eectrica and Eectronic Engineering, Bikent Univerity, Ankara, Turkey hitay@bikentedutr //$ IEEE 5633 to the origin Arbitrary deay term in addition to inputoutput deay were conidered in 9 with decentraized controer tructure Retricting the deigned controer to be PD and PID impoe thee retriction on the number of untabe pant poe 3 The reut in thi paper appy to much wider cae of SISO and MIMO ytem by aowing the order of the controer to exceed that of PD or PID The advantage of integra-action and impe ow-order impementation are ti part of the ynthei technique Notation: Let C, R, R + denote compex, rea, and poitive rea number The extended coed right-haf compex pane i U = { C Re() 0} { }; R p denote rea proper rationa function (of ); S R p i the tabe ubet with no poe in U; M(S) i the et of matrice with entrie in S ; I r i the r r identity matrix The pace H i the et of a bounded anaytic function in C + For h H, the norm i defined a h = e up C + h(), where e up denote the eentia upremum A matrix-vaued function H i in M(H ) if a it entrie are in H ; in thi cae H = e up C + σ(h()), where σ denote the maximum inguar vaue Since a norm of interet here are H norm, we drop the norm ubcript, ie, From the induced L 2 gain point of view, a ytem whoe tranfer-matrix i H i tabe iff H M(H ) A quare tranfer-matrix H M(H ) i unimoduar iff H M(H ) We drop () in tranfer-matrice uch a G(); ue δ(n) to denote the degree of the poynomia n(); ue a m or a a 2 a m to denote the (m m) ona matrix, whoe ona entrie are a,,a m For G R m m p we ue coprime factorization over S ; ie, G = Y X denote a eft-coprime-factorization (LCF), where X,Y S m m, dety ( ) 0 For the deayed pant cae, we ue coprime factorization over H ; ie, Ĝ = Ŷ X denote a eft-coprime-factorization (LCF), where X,Ŷ H m m II PROBLEM DESCRIPTION Conider the feedback ytem Sy(Ĝ,C) in Fig ; C R m m p i the tranfer-function of the controer and Ĝ i the tranfer-function of the pant with time deay It i aumed that the feedback ytem i we-poed and the deay-free part of the pant (ie, the pant without the time deay term) and the controer have no untabe hiddenmode With u, v, w, y a the input and output vector, the coed-oop tranfer-matrix Ĥ from (u,v) to (w,y) i C(I + Ĥ = ĜC) C(I + ĜC) Ĝ () ĜC(I + ĜC) (I + ĜC) Ĝ Let the (input-error) tranfer-function from u to e be denoted

2 v u e C w Ĝ y Fig The feedback ytem Sy( b G, C) by H eu and et the (input-output) tranfer-function from u to y be denoted by H yu ; then H eu = (I +ĜC) = I ĜC(I +ĜC) = I H yu (2) Definition : a) The feedback ytem Sy(Ĝ,C) hown in Fig, i tabe if the coed-oop map Ĥ i in M(H ) b) The controer C tabiize Ĝ if C i proper and Sy(Ĝ,C) i tabe c) The ytem Sy(Ĝ,C) i tabe and ha integraaction if the coed-oop tranfer-function from (u,v) to (w,y) i tabe, and the (input-error) tranfer-function H eu ha bocking-zero at = 0 d) The controer C i aid to be an integra-action controer if C tabiize Ĝ and D(0) = 0 for any RCF C = ND Let Ĝ = Ŷ X, where Ŷ, X M(H ) Let C = ND be an RCF, where D,N S m m, detd( ) 0 Then C tabiize Ĝ if and ony if M M(H ), where M := Ŷ D + X N (3) Suppoe that the ytem Sy(Ĝ,C) i tabe and that tep input reference are appied at u(t) The teady-tate error e(t) due to tep input at u(t) goe to zero a t if and ony if H eu (0) = 0 Therefore, by Definition -(c), the tabe ytem Sy(Ĝ,C) achieve aymptotic tracking of contant reference input with zero teady-tate error if and ony if it ha integra-action By (3), write H eu = (I + ĜC) = DM Ŷ Then by Definition -(d), Sy(Ĝ,C) ha integraaction if C = ND i an integra-action controer ince D(0) = 0 impie H eu (0) = (DM Ŷ )(0) = 0 The ytem Sy(Ĝ,C) woud ao have integra-action if every entry of the MIMO pant ha poe at = 0 ince Ŷ (0) = 0 impie H eu (0) = 0 even if the controer D(0) 0 Therefore, it i not a neceary condition to have integra-action controer for the ytem to have integra-action when Ŷ (0) = 0 However, for robut deign, integra-action i achieved with poe dupicating the dynamic tructure of the exogenou igna that the reguator ha to proce; thee integra-action controer obey the we-known interna mode principe 6 We aume throughout that Ĝ ha no tranmiion-zero at = 0 Thi condition i a neceary condition for exitence of integra-action controer: Let the (m m) matrix Ĝ() have (norma) rankg() = m If Ĝ admit an integra-action controer, then it ha no tranmiion-zero at = 0 III CONTROLLER SYNTHESIS We propoe finite-dimeniona tabiizing controer ynthei for certain cae of pant that have time deay The dicuion in Section III-A appie to a ca of SISO deay ytem In Section III-B, the reut are extended to a ca of MIMO pant with poe anywhere in the compex pane, 5634 but zero retricted to be in the OLHP Section III-C incude MIMO deay ytem whoe zero are unretricted, but the poe are either at the origin or in the tabe region A SISO pant of retarded type We conider SISO deay pant decribed a x() ν Ĝ() = y() + q(), q() = e hi q i (), (4) i= where x(),y(),q i () are poynomia with rea coefficient, δ(x) δ(y) > δ(q i ), the integer h i > 0, i =,,ν We aume that the finite zero of Ĝ are in the OLHP, ie, the poynomia x() i tricty Hurwitz Let r := δ(y) δ(x) 0 Let ξ() be any monic r-th order tricty Hurwitz poynomia; for exampe, ξ = ( + a) r for any a R + Define Y n := y() x()ξ(), Y d := q() x()ξ(), Ŷ = Y n + Y d, X := ξ() (5) Then X,Y n S, Y d H, and Ĝ = (Ŷ ) X = (Y n + Y d ) X Theorem preent a finite-dimeniona controer ynthei for coed-oop tabiity Thi deign give integraaction controer of order r when the reative degree of x()/y() i r, or of order when the reative degree of x()/y() i zero Theorem : (SISO tabiizing controer ynthei): Let Ĝ() be a in (4) For any monic r-th order tricty Hurwitz poynomia ξ(), et Ĝ() = Ŷ X = (Y n + Y d ) X be a in (5) a) If r = 0, then chooe any g R + Let o R + be uch that o > gĝ = + + gŷ (6) Then the controer C o in (7) tabiize Ĝ: ( + g) C o = o (7) b) If r, then chooe any monic, tricty Hurwitz poynomia ξ() of order r Define Θ a Θ() := Y n ( ) ξ()ĝ() Let R + be uch that = (Y n () + Y d ())Y n ( ) (8) > r Θ() (9) Then the controer C in (0) tabiize Ĝ: r ξ() C = ( + ) r r Y n( ) (0) Remark: In Theorem, the SISO controer C o in (7) for r = 0, C in (0) for r are biproper, and each ha a poe at = 0 providing integra-action The remaining (r ) poe of r-th order controer C in (0) are a in the OLHP Proof of Theorem : a) If r = 0, then X = in Ĝ = Ŷ X Let N =, D = Co ; then C o = ND i a

3 coprime factorization of the propoed controer in (7) By (3), C o tabiize Ĝ if and ony if M H, where M = X N + Ŷ D = + Ŷ () o ( + g) o (+g) Since o atifie (6), Ŷ () <, which i a ufficient condition for M H Therefore, C o tabiize Ĝ ince M H b) Define := ( + ) r r Let N = r, D = () ξ() Y n( ) Then C = ND i a coprime factorization of the propoed controer in (0) By (3), C tabiize Ĝ if and ony if M M(H ), where M = X N + Ŷ D = r ξ() + Ŷ ()() ξ() Y n( ) r = ( ( + ) r + () ( + ) r Ŷ ()Y n( ) ( + )r ) ξ() = ( + ( + ) r Ŷ ()Y n( ) ( + )r ) ξ() + )r = ( + Θ())( ( + ) r ξ() A ufficient condition for M to be in M(H ) i that (+) Θ() < We firt how that r ( + ) r = ( + )r r ( + ) r = r () ( ) For r, (+) r r =0 (+) = r r impie that the norm in () i greater than or equa to r/ We prove the norm in () i e than or equa to r/ by iteration: For r =, () hod ince (+) = / For r = 2, 2 (+)2 (+) = (+) = 2/ For r = 3, 2 + = / impie ( + )3 3 ( + ) 3 ( + )( + ) 2 2 ( + ) 2 = ( + ) + ( + )2 2 ( + ) 2 + ( + )2 2 ( + ) 2 + = 2 + = 3, hence, () hod Continuing imiary, uppoe that () hod for r and how that it hod for (r + ): ( + )r+ r+ ( + ) r+ = ( + )( + )r r ( + ) ( + ) r + ( + ) r + ( + ) r + = r + = r +, hence, () hod In (8), Y n ()Y n ( ) = impie Y n ()Y n ( ) I S Since δ(q) < δ(y) = δ(xξ), we have Y d ()Y n ( ) H Therefore, Θ() = Y n ()Y n ( ) I + Y d ()Y n ( ) H Since atifie (9), (+) Θ() r r Θ() < Therefore, C in (0) tabiize Ĝ ince M H 5635 B MIMO pant with unretricted poe We conider (m m) MIMO pant with deay, where the deay are a in the denominator matrix Ŷ M(H ) of Ĝ = Ŷ X Therefore, X i deay-free and we denote it by X M(S) We aume that Ĝ can be written a Ĝ = Ŷ X ; Ŷ = Y n + Y d, Y n () M(S), ν dety n ( ) 0, Y d = e hi Q i (), Q i ( ) = 0 (2) i= We aume that the tranmiion-zero of Ĝ are a in the OLHP and at infinity, ie, rankx( ) < m but rankx() = m for a C + With n k and d k a poynomia, write X nk () () = (3) d k () k, {,,m} Since the tranmiion-zero of Ĝ are a in the OLHP, X ha no poe in the coed right-haf compex pane C + (ie, d k are tricty Hurwitz) but may have poe at infinity, ie, X may be improper Define the integer r k and r a { δ(nk ) δ(d r k := k ), if δ(n k ) > δ(d k ) 0, if δ(n k ) δ(d k ) r := max k m r k, =,,m (4) Let ξ () be any monic r -th order tricty Hurwitz poynomia, =,,m; eg, ξ () = ( + a) r for a R + Define () := ξ () ξ 2 () ξ m () (5) If r = 0, then ξ = Athough X may be improper, X n i tabe ince k () d k ()ξ () S Define Ŷ ( ) := (X()Ĝ() ) ; ie, Y j ( ) = (Ĝ()X ()) By (2), Ŷ ( ) = Y n( ) For thi ca of MIMO pant, Theorem 2 preent a finite-dimeniona controer ynthei with integra-action Theorem 2: (MIMO tabiizing controer ynthei): Let Ĝ = Ŷ X = (Y n + Y d ) X be a in (2) Define Θ a Θ() := Ŷ ()Y n( ) I (6) For =,,m, define ρ a {, if r = 0 ρ := r, if r Let R + be uch that For =,,m, define a (7) > max ρ Θ (8) () := ( + ) ρ ρ (9) Then the controer C in (20) tabiize Ĝ: C = X ρ ρ2 () () 2 () ρm m () Y n ( ) (20)

4 m Remark: In Theorem 2, the ona term of ρ in the controer C in (20) a have poe at = 0 and hence, C ha integra-action The term correponding to r = 0 are in the form The term correponding to r are in the form ( r (+) r ), with one poe at = 0, and the r remaining (r ) poe a in the OLHP Proof of Theorem 2: For =,,m, define η a η := ( + ) ρ r ; (2) η = (+) if r = 0, and η = if r If r = 0, then et ξ = If r, then chooe any monic, tricty Hurwitz poynomia ξ () of order r With a in (5), et N = X η ρ η 2 ρ2 η m ρm m = X ρ ξ ()( + ) ρ, r D = Y n ( ) η () m = Y n ( ) () ξ ()( + ) ρ r m (22) Then C = ND i a right-factorization of the propoed controer in (20) By (3), C tabiize Ĝ if and ony if M M(H ), where M = X N + Ŷ D = XX ρ η m + Ŷ ()Y n( ) η m ρ m = ξ ()( + ) ρ r m + Ŷ ()Y n( ) () ξ ()( + ) ρ r ρ m = ( ( + ) ρ m + Ŷ ()Y n( ) () ( + ) r m ( + ) r ) ξ () m = (I + Ŷ ()Y n( ) () I ( + ) ρ ) ( + ) r m ξ () m () ( + ) r m = (I+Θ() ( + ) ρ ) ξ m The entrie of () (+) for r ρ = 0 have norm (+) = The entrie for r have norm ( + )r r ( + ) r = r = ρ a hown in the proof of Theorem Therefore, m () ρ ( + ) ρ = max m In (6), Y n ()Y n ( ) χ M(S) ince Y n ()Y n ( ) = I Since Q i ( ) = 0 by (2), 5636 Y d ()Y n ( ) M(H ); hence, Θ H Therefore C tabiize Ĝ ince M M(H ) if m () ρ Θ() ( + ) ρ Θ() max <, which hod ince atifie (8) C MIMO pant with unretricted tranmiion-zero We conider (m m) MIMO pant with deay, where the deay are a in the numerator matrix X M(H ) of Ĝ = Ŷ X, ie, the denominator matrix Ŷ i deay-free and we denote it by Y M(S) Therefore, we aume that a eft-coprime factorization of Ĝ can be written a Ĝ = Y X ; Xij = e hij X ij ; i,j =,,m; (23) the integer h ij 0; Y M(S) i deay-free; X M(H ) and X ij denote it ij-th entry Suppoe that each ij-th entry X ij of X may contain any arbitrary deay term and that the deay are known If the finite-dimeniona part Y of the deayed pant Ĝ i tabe, then (23) impie that the entrie of Ĝ may contain a different arbitrary known deay term Let Ĝ have fu (norma) rank m Let Ĝ have no tranmiion zero at = 0, equivaenty, rank X(0) = m We ao aume that Y may have poe anywhere in the OLHP, but the ony U-poe of are a at = 0, ie, the ony C + -poe of Y are at the origin The entrie of Y may have different mutipicitie of poe at = 0 and ome entrie may have ony poe in the tabe region C\U Write Y () = Y k () k,,,m (24) For =,,m, define the integer γ k 0 be the number of poe of Y k () at = 0, and define γ a γ := max k m γ k ; (25) ie, γ 0 i the arget number of poe at = 0 of the entrie in the -th coumn of Y () For =,,m, athough Y k () S, (Y k () γ (+β) ) S for any β R γ + For thi ca of (MIMO or SISO) pant, Theorem 3 preent a finite-dimeniona controer ynthei; Coroary incude integra-action in the tabiizing controer ynthei Theorem 3: (MIMO tabiizing controer ynthei): Let Ĝ = Y X be a in (23) Define Φ a Φ() := X()X(0) I (26) Chooe β R + uch that β < Φ() (27) max γ For =,,m, define ψ () := ( + β) γ γ (28) Then the controer C in (29) tabiize Ĝ: ψ () 0 0 γ ψ 0 2() C = X(0) 0 γ 2 Y () (29) 0 0 ψ m() γm

5 Coroary : (Integra-action controer ynthei): Under the aumption of Theorem 3, chooe β R + uch that β < For =,,m, define Φ() (30) + max γ ψ () := ( + β) +γ +γ (3) Then the integra-action controer C in (32) tabiize Ĝ: ψ () 0 0 +γ ψ C = X(0) 0 2() 0 +γ 2 Y () (32) Proof of Theorem 3: Let Ψ() := 0 0 ψ() (+β) γ Define N := X(0) Ψ(), D := Y γ (+β) γ ψ 2() (+β) γ 2 ψ m() +γm ψ m() (+β) γm (33) γ 2 (+β) γm γ 2 (+β) (34) γm Since the order of the poynomia ψ () i (γ ), the tricty-proper term S are tabe and hence, ψ () (+β) γ N = X(0) Ψ() M(S) Since (Y k () γ (+β) γ ) S, the matrix D M(S) Therefore, C = ND i a rightfactorization of the propoed controer in (29) By (3), C tabiize Ĝ if and ony if M M(H ), where M = X N + Y D = X()X(0) Ψ() + Y ()Y () γ (+β) γ = X()X(0) Ψ + m γ (+β) γ m = X()X(0) IΨ() = X()X(0) IΨ() = Φ()Ψ() A ufficient condition for M to be in M(H ) i that Φ()Ψ() < To find Ψ(), we firt how that ψ () ( + β) γ = ( + β)γ γ ( + β) γ = β γ (35) For ( γ = 0, (35) obviouy hod For γ, (+β) γ γ ) (+β) = βγ γ impie that the norm in (35) = i greater than or equa to β γ We prove the norm in (35) i e than or equa to β γ For γ =, (35) hod ince β +β = β For γ = 2, (+β)2 2 (+β) = 2β +β2 2 (+β) = 2 2β and hence, (35) hod For γ = 3, +β = impie (+β)3 3 (+β) 3 +β (+β)2 2 (+β) +β = (2β+β) = 2 3β and hence, (35) hod Continuing imiary, uppoe that (35) hod for γ and how that it hod for (γ + ): ( + β)γ+ γ+ ( + β) γ + + β ψ () ( + β) γ + β = (βγ + β) = β(γ + ) and hence, (35) hod Now (35) impie Ψ() = β max γ Since β atifie (27), Φ()Ψ() Φ() Ψ() = β max γ Φ() < Therefore, C tabiize Ĝ ince M M(H ) Proof of Coroary : Let ψ Ψ() := () ψ 2() ψ m() (+β) +r (+β) +r 2 (+β) +rm (36) With D a in (34), et Ñ = X(0) Ψ(), D = (+β) D Then C = Ñ D i a right-factorization of the propoed controer in (32); ince D(0) = 0, by Definition -(d), C i an integra-action controer We how that by (3), C tabiize Ĝ if and ony if M M(H ), where M = X Ñ + Y D = X Ñ + (+β) Y D = X()X(0) Ψ() + (+β) r m (+β) r = X()X(0) I Ψ() = 5637 X()X(0) I Ψ() = Φ() Ψ() A ufficient condition for M M(H ) i Φ() Ψ() <, where, by (35), Ψ() = β ( + max γ ) Since β atifie (30), Φ() Ψ() Φ() Ψ() = β ( + max γ ) Φ() < Therefore, the integra-action controer C tabiize Ĝ ince M M(H ) Exampe : Conider IV EXAMPLES ( + ) Ĝ() = ( ) + 2( )e h + 5e h2 x() = y() + q(), h 2 = π 2 h (37) The pant Ĝ i in the ca conidered in Section III-A Since the reative degree r =, the controer a in (0) i a firt order controer with integra action (ie a PI controer) Let ξ() = ( + b) for a free parameter b > 0; define Y n () = ( ) ( + )( + b), Y d() = 2( )e h + 5e h2 ( + )( + b) With Θ a in (8), Θ() = (2( )e h + 5e h2 ) (3 + b) b + 2, ( + )( + b) the minimum vaue Θ() of atifying (9) i hown in Fig 2 for variou h For h 0, 25, if we chooe b = 2, then = 8 atifie (9) The controer in (0) i C() = 8(+2) Thi feedback ytem i tabe if the tranformed characteritic equation + (+) Θ() = 0 ha no root in the coed right haf pane Since Θ < and ( + ) = /, the ma gain theorem impie tabiity Now change the pant Ĝ in (37) to Ĝ w = x() ( p) w y() + q(), w > 0, p R, and x(),y(),q() are the ame a in (37) The reative degree become r = w+, and the deayed part

6 min increaing direction of h h = 0, 05, 0, 20, b Fig 2 Exampe : b veru Θ for variou h Y d () of Ŷ () remain the ame The new Y n i Y nw = ( p)w y() ( + )( + b) w, and Y nw ( ) = Y n ( ) = We re-cacuate Θ() = Y nw ()Y n ( ) I + Y d ()Y n ( ) Y nw ()Y n ( ) I + Y d ()Y n ( ) =: Θ n () + Θ d () A condition on more conervative than (9) i > r Θ n () + r Θ d () For exampe, if p = 05 and w =, then r = 2 and Θ() = 45 for h 0, 25 Therefore, (9) hod if we chooe = 30 The controer in (0) i then given by C() = 900(+2) (+60) Exampe 2: For h 2 = π 2 h, conider Ĝ() = Ŷ () X() = + 2e h e 0 h The MIMO pant Ĝ i in the ca in Section III-B We have Ŷ ( ) = Y n ( ) =, r 0 = r 2 = We compute Θ() = 05 2e h (+e h ) + (2+3e h 2 ) + It can be hown that for a h 00, 25, (8) i atified for = 5 Hence C() a in (20) tabiize the given Ĝ: C() = (+05)(+2) ( ) +05 (+2)( 2 ++5) ( ) Exampe 3: Conider the foowing pant Ĝ, which i in the ca conidered in Section III-C: Ĝ() = K( + z)e h γ ( + p )( + p 2 ), where γ, p,p 2 R +, z R \ {0}, and h 0; Note that z may be poitive or negative, ie, Ĝ may have a finite zero in the right-haf compex pane Write Ĝ = Y X a Ĝ = Y X γ K( + z)e h = ( ( + b) ) ( γ ( + b) γ ( + p )( + p 2 ) ) 5638 for any b R + With X(0) = bγ p p 2 Kz, and Φ a in (26), et β R + atify (27), ie, β < γ Φ, where Φ = b γ p p 2 K( + z)e h Kz( + b) γ ( + p )( + p 2 ) Then C = X(0) (+β) γ γ Y = a p p2 (+β)γ γ γ K z (+b) a γ in (29) i a controer that tabiize Ĝ The controer C i tabe, and it order i γ, the ame a the number of poe of Ĝ at = 0, which i e than the pant order Let β R + atify (30), ie, β < +γ Φ Then an integraaction controer C a in (32) that tabiize Ĝ i C = X(0) (+β) γ + γ + Y = bγ p p 2(+β) γ+ γ+ γ + K z (+b) For γ exampe, if γ =, then C and C become C = X(0) β b p p2 β Y = K z (+b), where β < Φ, and C = X(0) (2β +β 2 ) Y = b p p2(2β +β2 ) 2 K z (+b), where β < 2 Φ V CONCLUSIONS We propoed finite-dimeniona controer deign for certain cae of SISO and MIMO ytem ubject to deay Thee deign achieve coed-oop tabiity and integraaction The controer order matche the reative degree of the finite-dimeniona part of the pant for the pant in Section III-A-III-B or the number of pant poe at the origin in Section III-C Performance pecification beyond aymptotic tracking of contant reference are not within the cope of thi tudy Future work wi focu on expanding the pant cae to thoe that aow finite right-haf pane zero whie not retricting the ocation of untabe poe REFERENCES C Bonnet, J R Partington, Anayi of fractiona deay ytem of retarded and neutra type, Automatica, vo 38, pp 3338, C Bonnet, JR Partington, Stabiization of ome fractiona deay ytem of neutra type, Automatica, vo 43, pp , RF Curtain, K Gover, Robut tabiization of infinite-dimeniona ytem by finite-dimeniona controer, Sytem Contro Letter, 7 (), pp 447, RF Curtain, H Zwart, An Introduction to Infinite-dimeniona Linear Sytem Theory, Text in Appied Mathematic, Vo 2, Springer, New York, C Foia, H Özbay, A Tannenbaum, Robut Contro of Infinite Dimeniona Sytem, LNCIS 209, Springer-Verag, London, B A Franci and W A Wonham, The interna mode principe for inear mutivariabe reguator, Appied Mathematic & Optimization, 2:2, pp 70-95, K Gu, V L Kharitonov, J Chen, Stabiity of Time-Deay Sytem, Birkhäuer, Boton, A N Gündeş, H Özbay, and A B Özgüer, PID controer ynthei for a ca of untabe MIMO pant with I/O Deay, Automatica, vo 43, no, pp 35-42, A N Gündeş, H Özbay, Reiabe decentraized contro of deayed MIMO pant, Int Jour Contro, vo 83, no 3, pp , S I Nicuecu, Deay Effect on Stabiity: A Robut Contro Approach, LNCIS, vo 269, Heideberg: Springer-Verag, 200 H Özbay, A N Gündeş, Integra action controer for ytem with time deay, in Topic in Time Deay Sytem Anayi, Agorithm and Contro, LNCIS 388, J J Loieau, W Michie, S-I Nicuecu, R Sipahi (Ed), pp , Springer-Verag, London, S Skogetad, Simpe anaytic rue for mode reduction and PID controer tuning, Jour Proce Contro, vo 3, pp , G J Siva, A Datta, S P Bhattacharyya, PID Controer for Time- Deay Sytem, Birkhäuer, Boton, 2005

Then C pid (s) S h -stabilizes G(s) if and only if Ĉpid(ŝ) S 0 - stabilizes Ĝ(ŝ). For any ρ R +, an RCF of Ĉ pid (ŝ) is given by

Then C pid (s) S h -stabilizes G(s) if and only if Ĉpid(ŝ) S 0 - stabilizes Ĝ(ŝ). For any ρ R +, an RCF of Ĉ pid (ŝ) is given by 9 American Control Conference Hyatt Regency Riverfront, St. Loui, MO, USA June -, 9 WeC5.5 PID Controller Synthei with Shifted Axi Pole Aignment for a Cla of MIMO Sytem A. N. Gündeş and T. S. Chang Abtract

More information

Well-Posedness of Feedback Loop:

Well-Posedness of Feedback Loop: ntena Stabiity We-oedne of Feedback Loop: onide the foowing feedback ytem - u u p d i d y Let be both pope tanfe function. Howeve u n d di 3 3 ote that the tanfe function fom the extena igna n d d to u

More information

VIII. Addition of Angular Momenta

VIII. Addition of Angular Momenta VIII Addition of Anguar Momenta a Couped and Uncouped Bae When deaing with two different ource of anguar momentum, Ĵ and Ĵ, there are two obviou bae that one might chooe to work in The firt i caed the

More information

THE PARAMETERIZATION OF ALL TWO-DEGREES-OF-FREEDOM SEMISTRONGLY STABILIZING CONTROLLERS. Tatsuya Hoshikawa, Kou Yamada and Yuko Tatsumi

THE PARAMETERIZATION OF ALL TWO-DEGREES-OF-FREEDOM SEMISTRONGLY STABILIZING CONTROLLERS. Tatsuya Hoshikawa, Kou Yamada and Yuko Tatsumi International Journal of Innovative Computing, Information Control ICIC International c 206 ISSN 349-498 Volume 2, Number 2, April 206 pp. 357 370 THE PARAMETERIZATION OF ALL TWO-DEGREES-OF-FREEDOM SEMISTRONGLY

More information

Research Article Simplicity and Commutative Bases of Derivations in Polynomial and Power Series Rings

Research Article Simplicity and Commutative Bases of Derivations in Polynomial and Power Series Rings ISRN Agebra Voume 2013 Artice ID 560648 4 page http://dx.doi.org/10.1155/2013/560648 Reearch Artice Simpicity and Commutative Bae of Derivation in Poynomia and Power Serie Ring Rene Batazar Univeridade

More information

EE Control Systems LECTURE 14

EE Control Systems LECTURE 14 Updated: Tueday, March 3, 999 EE 434 - Control Sytem LECTURE 4 Copyright FL Lewi 999 All right reerved ROOT LOCUS DESIGN TECHNIQUE Suppoe the cloed-loop tranfer function depend on a deign parameter k We

More information

CHAPTER 4 DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL

CHAPTER 4 DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL 98 CHAPTER DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL INTRODUCTION The deign of ytem uing tate pace model for the deign i called a modern control deign and it i

More information

Linear System Fundamentals

Linear System Fundamentals Linear Sytem Fundamental MEM 355 Performance Enhancement of Dynamical Sytem Harry G. Kwatny Department of Mechanical Engineering & Mechanic Drexel Univerity Content Sytem Repreentation Stability Concept

More information

Active Sway Control of a Single Pendulum Gantry Crane System using Output-Delayed Feedback Control Technique

Active Sway Control of a Single Pendulum Gantry Crane System using Output-Delayed Feedback Control Technique 21 11th Int. Conf. Contro, Automation, Robotic and Viion Singapore, 7-1th December 21 Active Sway Contro of a Singe Penduum Gantry Crane Sytem uing Output-Deayed Feedback Contro echnique Rajeeb Dey, Nihant

More information

Root Locus Contents. Root locus, sketching algorithm. Root locus, examples. Root locus, proofs. Root locus, control examples

Root Locus Contents. Root locus, sketching algorithm. Root locus, examples. Root locus, proofs. Root locus, control examples Root Locu Content Root locu, ketching algorithm Root locu, example Root locu, proof Root locu, control example Root locu, influence of zero and pole Root locu, lead lag controller deign 9 Spring ME45 -

More information

SIMPLE LOW-ORDER AND INTEGRAL-ACTION CONTROLLER SYNTHESIS FOR MIMO SYSTEMS WITH TIME DELAYS

SIMPLE LOW-ORDER AND INTEGRAL-ACTION CONTROLLER SYNTHESIS FOR MIMO SYSTEMS WITH TIME DELAYS Appl. Comput. Math., V.10, N.2, 2011, pp.242-249 SIMPLE LOW-ORDER AND INTEGRAL-ACTION CONTROLLER SYNTHESIS FOR MIMO SYSTEMS WITH TIME DELAYS A.N. GÜNDEŞ1, A.N. METE 2 Abtact. A imple finite-dimenional

More information

Independent Joint Control

Independent Joint Control ME135 ADANCED OBOICS Independent oint Contro ee-hwan yu Schoo of Mechanica Engineering Introduction to Contro Contro: deterining the tie hitory of joint input to do a coanded otion Contro ethod are depend

More information

Chapter 13. Root Locus Introduction

Chapter 13. Root Locus Introduction Chapter 13 Root Locu 13.1 Introduction In the previou chapter we had a glimpe of controller deign iue through ome imple example. Obviouly when we have higher order ytem, uch imple deign technique will

More information

Methods for calculation of the coupling coefficients in the Coupling Cavity Model of arbitrary chain of resonators

Methods for calculation of the coupling coefficients in the Coupling Cavity Model of arbitrary chain of resonators Method for cacuation of the couping coefficient in the Couping Cavity Mode of arbitrary chain of reonator M.I. Ayzaty V.V.Mytrocheno Nationa Science Center Kharov Intitute of Phyic and echnoogy (NSC KIP)

More information

March 18, 2014 Academic Year 2013/14

March 18, 2014 Academic Year 2013/14 POLITONG - SHANGHAI BASIC AUTOMATIC CONTROL Exam grade March 8, 4 Academic Year 3/4 NAME (Pinyin/Italian)... STUDENT ID Ue only thee page (including the back) for anwer. Do not ue additional heet. Ue of

More information

LOW ORDER MIMO CONTROLLER DESIGN FOR AN ENGINE DISTURBANCE REJECTION PROBLEM. P.Dickinson, A.T.Shenton

LOW ORDER MIMO CONTROLLER DESIGN FOR AN ENGINE DISTURBANCE REJECTION PROBLEM. P.Dickinson, A.T.Shenton LOW ORDER MIMO CONTROLLER DESIGN FOR AN ENGINE DISTURBANCE REJECTION PROBLEM P.Dickinon, A.T.Shenton Department of Engineering, The Univerity of Liverpool, Liverpool L69 3GH, UK Abtract: Thi paper compare

More information

= 25 ohms, and line #2 has R c2

= 25 ohms, and line #2 has R c2 Soution for Aignment #3.A tranmiion-ine circuit i driven by a tep-function generator with V = 0 vot and RS =0 ohm. ine # ha ength of = cm. Both ine have the ame peed of trave, u = 0 cm/n. ine # ha characteritic

More information

CHAPTER 3. FUZZY LOGIC DIRECT TORQUE CONTROL.

CHAPTER 3. FUZZY LOGIC DIRECT TORQUE CONTROL. CHAPTER 3. FUZZY LOGIC DIRECT TORQUE CONTROL. 3. - Introduction. In DTC induction otor drive there are torque and fux rippe becaue none of the VSI tate i abe to generate the exact votage vaue required

More information

Automatic Control Systems. Part III: Root Locus Technique

Automatic Control Systems. Part III: Root Locus Technique www.pdhcenter.com PDH Coure E40 www.pdhonline.org Automatic Control Sytem Part III: Root Locu Technique By Shih-Min Hu, Ph.D., P.E. Page of 30 www.pdhcenter.com PDH Coure E40 www.pdhonline.org VI. Root

More information

Gain and Phase Margins Based Delay Dependent Stability Analysis of Two- Area LFC System with Communication Delays

Gain and Phase Margins Based Delay Dependent Stability Analysis of Two- Area LFC System with Communication Delays Gain and Phae Margin Baed Delay Dependent Stability Analyi of Two- Area LFC Sytem with Communication Delay Şahin Sönmez and Saffet Ayaun Department of Electrical Engineering, Niğde Ömer Halidemir Univerity,

More information

Chapter 4 Interconnection of LTI Systems

Chapter 4 Interconnection of LTI Systems Chapter 4 Interconnection of LTI Sytem 4. INTRODUCTION Block diagram and ignal flow graph are commonly ued to decribe a large feedback control ytem. Each block in the ytem i repreented by a tranfer function,

More information

Bogoliubov Transformation in Classical Mechanics

Bogoliubov Transformation in Classical Mechanics Bogoliubov Tranformation in Claical Mechanic Canonical Tranformation Suppoe we have a et of complex canonical variable, {a j }, and would like to conider another et of variable, {b }, b b ({a j }). How

More information

Mathematical modeling of control systems. Laith Batarseh. Mathematical modeling of control systems

Mathematical modeling of control systems. Laith Batarseh. Mathematical modeling of control systems Chapter two Laith Batareh Mathematical modeling The dynamic of many ytem, whether they are mechanical, electrical, thermal, economic, biological, and o on, may be decribed in term of differential equation

More information

MEM 355 Performance Enhancement of Dynamical Systems Root Locus Analysis

MEM 355 Performance Enhancement of Dynamical Systems Root Locus Analysis MEM 355 Performance Enhancement of Dynamical Sytem Root Locu Analyi Harry G. Kwatny Department of Mechanical Engineering & Mechanic Drexel Univerity Outline The root locu method wa introduced by Evan in

More information

The Root Locus Method

The Root Locus Method The Root Locu Method MEM 355 Performance Enhancement of Dynamical Sytem Harry G. Kwatny Department of Mechanical Engineering & Mechanic Drexel Univerity Outline The root locu method wa introduced by Evan

More information

A Simplified Methodology for the Synthesis of Adaptive Flight Control Systems

A Simplified Methodology for the Synthesis of Adaptive Flight Control Systems A Simplified Methodology for the Synthei of Adaptive Flight Control Sytem J.ROUSHANIAN, F.NADJAFI Department of Mechanical Engineering KNT Univerity of Technology 3Mirdamad St. Tehran IRAN Abtract- A implified

More information

An Inequality for Nonnegative Matrices and the Inverse Eigenvalue Problem

An Inequality for Nonnegative Matrices and the Inverse Eigenvalue Problem An Inequality for Nonnegative Matrice and the Invere Eigenvalue Problem Robert Ream Program in Mathematical Science The Univerity of Texa at Dalla Box 83688, Richardon, Texa 7583-688 Abtract We preent

More information

Control Systems Engineering ( Chapter 7. Steady-State Errors ) Prof. Kwang-Chun Ho Tel: Fax:

Control Systems Engineering ( Chapter 7. Steady-State Errors ) Prof. Kwang-Chun Ho Tel: Fax: Control Sytem Engineering ( Chapter 7. Steady-State Error Prof. Kwang-Chun Ho kwangho@hanung.ac.kr Tel: 0-760-453 Fax:0-760-4435 Introduction In thi leon, you will learn the following : How to find the

More information

SOME RESULTS ON INFINITE POWER TOWERS

SOME RESULTS ON INFINITE POWER TOWERS NNTDM 16 2010) 3, 18-24 SOME RESULTS ON INFINITE POWER TOWERS Mladen Vailev - Miana 5, V. Hugo Str., Sofia 1124, Bulgaria E-mail:miana@abv.bg Abtract To my friend Kratyu Gumnerov In the paper the infinite

More information

Evolutionary Algorithms Based Fixed Order Robust Controller Design and Robustness Performance Analysis

Evolutionary Algorithms Based Fixed Order Robust Controller Design and Robustness Performance Analysis Proceeding of 01 4th International Conference on Machine Learning and Computing IPCSIT vol. 5 (01) (01) IACSIT Pre, Singapore Evolutionary Algorithm Baed Fixed Order Robut Controller Deign and Robutne

More information

Multivariable Control Systems

Multivariable Control Systems Lecture Multivariable Control Sytem Ali Karimpour Aociate Profeor Ferdowi Univerity of Mahhad Lecture Reference are appeared in the lat lide. Dr. Ali Karimpour May 6 Uncertainty in Multivariable Sytem

More information

TRIPLE SOLUTIONS FOR THE ONE-DIMENSIONAL

TRIPLE SOLUTIONS FOR THE ONE-DIMENSIONAL GLASNIK MATEMATIČKI Vol. 38583, 73 84 TRIPLE SOLUTIONS FOR THE ONE-DIMENSIONAL p-laplacian Haihen Lü, Donal O Regan and Ravi P. Agarwal Academy of Mathematic and Sytem Science, Beijing, China, National

More information

Problem Set 8 Solutions

Problem Set 8 Solutions Deign and Analyi of Algorithm April 29, 2015 Maachuett Intitute of Technology 6.046J/18.410J Prof. Erik Demaine, Srini Devada, and Nancy Lynch Problem Set 8 Solution Problem Set 8 Solution Thi problem

More information

A PLC BASED MIMO PID CONTROLLER FOR MULTIVARIABLE INDUSTRIAL PROCESSES

A PLC BASED MIMO PID CONTROLLER FOR MULTIVARIABLE INDUSTRIAL PROCESSES ABCM Sympoium Serie in Mechatronic - Vol. 3 - pp.87-96 Copyright c 8 by ABCM A PLC BASE MIMO PI CONOLLE FO MULIVAIABLE INUSIAL POCESSES Joé Maria Galvez, jmgalvez@ufmg.br epartment of Mechanical Engineering

More information

Hyperbolic Partial Differential Equations

Hyperbolic Partial Differential Equations Hyperbolic Partial Differential Equation Evolution equation aociated with irreverible phyical procee like diffuion heat conduction lead to parabolic partial differential equation. When the equation i a

More information

CHEAP CONTROL PERFORMANCE LIMITATIONS OF INPUT CONSTRAINED LINEAR SYSTEMS

CHEAP CONTROL PERFORMANCE LIMITATIONS OF INPUT CONSTRAINED LINEAR SYSTEMS Copyright 22 IFAC 5th Triennial World Congre, Barcelona, Spain CHEAP CONTROL PERFORMANCE LIMITATIONS OF INPUT CONSTRAINED LINEAR SYSTEMS Tritan Pérez Graham C. Goodwin Maria M. Serón Department of Electrical

More information

List coloring hypergraphs

List coloring hypergraphs Lit coloring hypergraph Penny Haxell Jacque Vertraete Department of Combinatoric and Optimization Univerity of Waterloo Waterloo, Ontario, Canada pehaxell@uwaterloo.ca Department of Mathematic Univerity

More information

Control of Delayed Integrating Processes Using Two Feedback Controllers R MS Approach

Control of Delayed Integrating Processes Using Two Feedback Controllers R MS Approach Proceeding of the 7th WSEAS International Conference on SYSTEM SCIENCE and SIMULATION in ENGINEERING (ICOSSSE '8) Control of Delayed Integrating Procee Uing Two Feedback Controller R MS Approach LIBOR

More information

Multicolor Sunflowers

Multicolor Sunflowers Multicolor Sunflower Dhruv Mubayi Lujia Wang October 19, 2017 Abtract A unflower i a collection of ditinct et uch that the interection of any two of them i the ame a the common interection C of all of

More information

Given the following circuit with unknown initial capacitor voltage v(0): X(s) Immediately, we know that the transfer function H(s) is

Given the following circuit with unknown initial capacitor voltage v(0): X(s) Immediately, we know that the transfer function H(s) is EE 4G Note: Chapter 6 Intructor: Cheung More about ZSR and ZIR. Finding unknown initial condition: Given the following circuit with unknown initial capacitor voltage v0: F v0/ / Input xt 0Ω Output yt -

More information

Duality Model of TCP/AQM

Duality Model of TCP/AQM Optimization and Contro of Network Duaity Mode of TCP/AQM Lijun Chen 02/08/2016 Agenda Utiity maimization and dua decompoition An introduction to TCP congetion contro A genera dynamica mode of TCP/AQM

More information

A Constraint Propagation Algorithm for Determining the Stability Margin. The paper addresses the stability margin assessment for linear systems

A Constraint Propagation Algorithm for Determining the Stability Margin. The paper addresses the stability margin assessment for linear systems A Contraint Propagation Algorithm for Determining the Stability Margin of Linear Parameter Circuit and Sytem Lubomir Kolev and Simona Filipova-Petrakieva Abtract The paper addree the tability margin aement

More information

OBSERVER DESIGN FOR DISCRETE-TIME LINEAR SWITCHING SYSTEMS 1

OBSERVER DESIGN FOR DISCRETE-TIME LINEAR SWITCHING SYSTEMS 1 OBSERVER DESIGN FOR DISCRETE-TIME LINEAR SWITCHING SYSTEMS 1 E. De Santi, M.D. Di Benedetto Department of Electrical Engineering and Computer Science, Univerity of L Aquila. Email: (deanti,dibenede@ing.univaq.it

More information

Minimal state space realization of MIMO systems in the max algebra

Minimal state space realization of MIMO systems in the max algebra KULeuven Department of Electrical Engineering (ESAT) SISTA Technical report 94-54 Minimal tate pace realization of MIMO ytem in the max algebra B De Schutter and B De Moor If you want to cite thi report

More information

Root Locus Diagram. Root loci: The portion of root locus when k assume positive values: that is 0

Root Locus Diagram. Root loci: The portion of root locus when k assume positive values: that is 0 Objective Root Locu Diagram Upon completion of thi chapter you will be able to: Plot the Root Locu for a given Tranfer Function by varying gain of the ytem, Analye the tability of the ytem from the root

More information

SOLUTIONS TO ALGEBRAIC GEOMETRY AND ARITHMETIC CURVES BY QING LIU. I will collect my solutions to some of the exercises in this book in this document.

SOLUTIONS TO ALGEBRAIC GEOMETRY AND ARITHMETIC CURVES BY QING LIU. I will collect my solutions to some of the exercises in this book in this document. SOLUTIONS TO ALGEBRAIC GEOMETRY AND ARITHMETIC CURVES BY QING LIU CİHAN BAHRAN I will collect my olution to ome of the exercie in thi book in thi document. Section 2.1 1. Let A = k[[t ]] be the ring of

More information

PI control system design for Electromagnetic Molding Machine based on Linear Programing

PI control system design for Electromagnetic Molding Machine based on Linear Programing PI control ytem deign for Electromagnetic Molding Machine baed on Linear Programing Takayuki Ihizaki, Kenji Kahima, Jun-ichi Imura*, Atuhi Katoh and Hirohi Morita** Abtract In thi paper, we deign a PI

More information

EE Control Systems LECTURE 6

EE Control Systems LECTURE 6 Copyright FL Lewi 999 All right reerved EE - Control Sytem LECTURE 6 Updated: Sunday, February, 999 BLOCK DIAGRAM AND MASON'S FORMULA A linear time-invariant (LTI) ytem can be repreented in many way, including:

More information

Copyright 1967, by the author(s). All rights reserved.

Copyright 1967, by the author(s). All rights reserved. Copyright 1967, by the author(). All right reerved. Permiion to make digital or hard copie of all or part of thi work for peronal or claroom ue i granted without fee provided that copie are not made or

More information

Convex Hulls of Curves Sam Burton

Convex Hulls of Curves Sam Burton Convex Hull of Curve Sam Burton 1 Introduction Thi paper will primarily be concerned with determining the face of convex hull of curve of the form C = {(t, t a, t b ) t [ 1, 1]}, a < b N in R 3. We hall

More information

EXTENDED STABILITY MARGINS ON CONTROLLER DESIGN FOR NONLINEAR INPUT DELAY SYSTEMS. Otto J. Roesch, Hubert Roth, Asif Iqbal

EXTENDED STABILITY MARGINS ON CONTROLLER DESIGN FOR NONLINEAR INPUT DELAY SYSTEMS. Otto J. Roesch, Hubert Roth, Asif Iqbal EXTENDED STABILITY MARGINS ON CONTROLLER DESIGN FOR NONLINEAR INPUT DELAY SYSTEMS Otto J. Roech, Hubert Roth, Aif Iqbal Intitute of Automatic Control Engineering Univerity Siegen, Germany {otto.roech,

More information

Laplace Examples, Inverse, Rational Form

Laplace Examples, Inverse, Rational Form Lecture 20 Outine: Lapace Eampe, Invere, Rationa Form Announcement: HW 5 poted, due Friday ore Lapace Tranform Eampe Invere Lapace Tranform Rationa Lapace Tranform ROC for Right/Left/Toided Signa agnitude/phae

More information

7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281

7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281 72 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 28 and i 2 Show how Euler formula (page 33) can then be ued to deduce the reult a ( a) 2 b 2 {e at co bt} {e at in bt} b ( a) 2 b 2 5 Under what condition

More information

Robust Decentralized Design of H -based Frequency Stabilizer of SMES

Robust Decentralized Design of H -based Frequency Stabilizer of SMES International Energy Journal: Vol. 6, No., Part, June 005-59 Robut Decentralized Deign of H -baed Frequency Stabilizer of SMES www.erd.ait.ac.th/reric C. Vorakulpipat *, M. Leelajindakrirerk *, and I.

More information

Non-linear Control and Observer Design for the Self-energizing Electro-Hydraulic Brake

Non-linear Control and Observer Design for the Self-energizing Electro-Hydraulic Brake 9. Non-inear Contro and Oberver Deign for the Sef-energizing Eectro-Hydrauic Brake Michae Kühnein, Juian Ewad, Hubertu Murrenhoff RWTH achen Univerity, Germany, Intitute for uid Power Drive and Contro

More information

NUMERICAL SOLUTION OF THE SYSTEM OF LINEAR FREDHOLM INTEGRAL EQUATIONS BASED ON DEGENERATING KERNELS

NUMERICAL SOLUTION OF THE SYSTEM OF LINEAR FREDHOLM INTEGRAL EQUATIONS BASED ON DEGENERATING KERNELS TWMS J Pure App Math V6, N1, 2015, pp111-119 NUMERICAL SOLUTION OF THE SYSTEM OF LINEAR FREDHOLM INTEGRAL EQUATIONS BASED ON DEGENERATING KERNELS S KARIMI 1, M JOZI 1 Abtract In thi paper, a new numerica

More information

Analysis and Design of a Third Order Phase-Lock Loop

Analysis and Design of a Third Order Phase-Lock Loop Analyi Deign of a Third Order Phae-Lock Loop DANIEL Y. ABRAMOVITCH Ford Aeropace Corporation 3939 Fabian Way, MS: X- Palo Alto, CA 94303 Abtract Typical implementation of a phae-lock loop (PLL) are econd

More information

CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS

CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.1 INTRODUCTION 8.2 REDUCED ORDER MODEL DESIGN FOR LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.3

More information

To describe a queuing system, an input process and an output process has to be specified.

To describe a queuing system, an input process and an output process has to be specified. 5. Queue (aiting Line) Queuing terminology Input Service Output To decribe a ueuing ytem, an input proce and an output proce ha to be pecified. For example ituation input proce output proce Bank Cutomer

More information

MODERN CONTROL SYSTEMS

MODERN CONTROL SYSTEMS MODERN CONTROL SYSTEMS Lecture 1 Root Locu Emam Fathy Department of Electrical and Control Engineering email: emfmz@aat.edu http://www.aat.edu/cv.php?dip_unit=346&er=68525 1 Introduction What i root locu?

More information

Solutions. Digital Control Systems ( ) 120 minutes examination time + 15 minutes reading time at the beginning of the exam

Solutions. Digital Control Systems ( ) 120 minutes examination time + 15 minutes reading time at the beginning of the exam BSc - Sample Examination Digital Control Sytem (5-588-) Prof. L. Guzzella Solution Exam Duration: Number of Quetion: Rating: Permitted aid: minute examination time + 5 minute reading time at the beginning

More information

1 Routh Array: 15 points

1 Routh Array: 15 points EE C28 / ME34 Problem Set 3 Solution Fall 2 Routh Array: 5 point Conider the ytem below, with D() k(+), w(t), G() +2, and H y() 2 ++2 2(+). Find the cloed loop tranfer function Y () R(), and range of k

More information

STOCHASTIC DIFFERENTIAL GAMES:THE LINEAR QUADRATIC ZERO SUM CASE

STOCHASTIC DIFFERENTIAL GAMES:THE LINEAR QUADRATIC ZERO SUM CASE Sankhyā : he Indian Journal of Statitic 1995, Volume 57, Serie A, Pt. 1, pp.161 165 SOCHASIC DIFFERENIAL GAMES:HE LINEAR QUADRAIC ZERO SUM CASE By R. ARDANUY Univeridad de Salamanca SUMMARY. hi paper conider

More information

ROBUST CONTROLLER DESIGN WITH HARD INPUT CONSTRAINTS. TIME DOMAIN APPROACH. Received August 2015; revised November 2015

ROBUST CONTROLLER DESIGN WITH HARD INPUT CONSTRAINTS. TIME DOMAIN APPROACH. Received August 2015; revised November 2015 International Journal of Innovative Computing, Information and Control ICIC International c 2016 ISSN 1349-4198 Volume 12, Number 1, February 2016 pp. 161 170 ROBUST CONTROLLER DESIGN WITH HARD INPUT CONSTRAINTS.

More information

Advanced D-Partitioning Analysis and its Comparison with the Kharitonov s Theorem Assessment

Advanced D-Partitioning Analysis and its Comparison with the Kharitonov s Theorem Assessment Journal of Multidiciplinary Engineering Science and Technology (JMEST) ISSN: 59- Vol. Iue, January - 5 Advanced D-Partitioning Analyi and it Comparion with the haritonov Theorem Aement amen M. Yanev Profeor,

More information

ONLINE APPENDIX: TESTABLE IMPLICATIONS OF TRANSLATION INVARIANCE AND HOMOTHETICITY: VARIATIONAL, MAXMIN, CARA AND CRRA PREFERENCES

ONLINE APPENDIX: TESTABLE IMPLICATIONS OF TRANSLATION INVARIANCE AND HOMOTHETICITY: VARIATIONAL, MAXMIN, CARA AND CRRA PREFERENCES ONLINE APPENDIX: TESTABLE IMPLICATIONS OF TRANSLATION INVARIANCE AND HOMOTHETICITY: VARIATIONAL, MAXMIN, CARA AND CRRA PREFERENCES CHRISTOPHER P. CHAMBERS, FEDERICO ECHENIQUE, AND KOTA SAITO In thi online

More information

MULTIPLE POSITIVE SOLUTIONS OF BOUNDARY VALUE PROBLEMS FOR P-LAPLACIAN IMPULSIVE DYNAMIC EQUATIONS ON TIME SCALES

MULTIPLE POSITIVE SOLUTIONS OF BOUNDARY VALUE PROBLEMS FOR P-LAPLACIAN IMPULSIVE DYNAMIC EQUATIONS ON TIME SCALES Fixed Point Theory, 5(24, No. 2, 475-486 http://www.math.ubbcluj.ro/ nodeacj/fptcj.html MULTIPLE POSITIVE SOLUTIONS OF BOUNDARY VALUE PROBLEMS FOR P-LAPLACIAN IMPULSIVE DYNAMIC EQUATIONS ON TIME SCALES

More information

Multi-dimensional Fuzzy Euler Approximation

Multi-dimensional Fuzzy Euler Approximation Mathematica Aeterna, Vol 7, 2017, no 2, 163-176 Multi-dimenional Fuzzy Euler Approximation Yangyang Hao College of Mathematic and Information Science Hebei Univerity, Baoding 071002, China hdhyywa@163com

More information

DYNAMIC MODELS FOR CONTROLLER DESIGN

DYNAMIC MODELS FOR CONTROLLER DESIGN DYNAMIC MODELS FOR CONTROLLER DESIGN M.T. Tham (996,999) Dept. of Chemical and Proce Engineering Newcatle upon Tyne, NE 7RU, UK.. INTRODUCTION The problem of deigning a good control ytem i baically that

More information

Manprit Kaur and Arun Kumar

Manprit Kaur and Arun Kumar CUBIC X-SPLINE INTERPOLATORY FUNCTIONS Manprit Kaur and Arun Kumar manpreet2410@gmail.com, arun04@rediffmail.com Department of Mathematic and Computer Science, R. D. Univerity, Jabalpur, INDIA. Abtract:

More information

Hybrid Projective Dislocated Synchronization of Liu Chaotic System Based on Parameters Identification

Hybrid Projective Dislocated Synchronization of Liu Chaotic System Based on Parameters Identification www.ccenet.org/ma Modern Applied Science Vol. 6, No. ; February Hybrid Projective Dilocated Synchronization of Liu Chaotic Sytem Baed on Parameter Identification Yanfei Chen College of Science, Guilin

More information

CONTROL OF INTEGRATING PROCESS WITH DEAD TIME USING AUTO-TUNING APPROACH

CONTROL OF INTEGRATING PROCESS WITH DEAD TIME USING AUTO-TUNING APPROACH Brazilian Journal of Chemical Engineering ISSN 004-6632 Printed in Brazil www.abeq.org.br/bjche Vol. 26, No. 0, pp. 89-98, January - March, 2009 CONROL OF INEGRAING PROCESS WIH DEAD IME USING AUO-UNING

More information

Reliability Analysis of Embedded System with Different Modes of Failure Emphasizing Reboot Delay

Reliability Analysis of Embedded System with Different Modes of Failure Emphasizing Reboot Delay International Journal of Applied Science and Engineering 3., 4: 449-47 Reliability Analyi of Embedded Sytem with Different Mode of Failure Emphaizing Reboot Delay Deepak Kumar* and S. B. Singh Department

More information

A Comparative Study on Control Techniques of Non-square Matrix Distillation Column

A Comparative Study on Control Techniques of Non-square Matrix Distillation Column IJCTA, 8(3), 215, pp 1129-1136 International Science Pre A Comparative Study on Control Technique of Non-quare Matrix Ditillation Column 1 S Bhat Vinayambika, 2 S Shanmuga Priya, and 3 I Thirunavukkarau*

More information

CONTROL SYSTEMS, ROBOTICS AND AUTOMATION Vol. VIII Decoupling Control - M. Fikar

CONTROL SYSTEMS, ROBOTICS AND AUTOMATION Vol. VIII Decoupling Control - M. Fikar DECOUPLING CONTROL M. Fikar Department of Proce Control, Faculty of Chemical and Food Technology, Slovak Univerity of Technology in Bratilava, Radlinkého 9, SK-812 37 Bratilava, Slovakia Keyword: Decoupling:

More information

Simple Observer Based Synchronization of Lorenz System with Parametric Uncertainty

Simple Observer Based Synchronization of Lorenz System with Parametric Uncertainty IOSR Journal of Electrical and Electronic Engineering (IOSR-JEEE) ISSN: 78-676Volume, Iue 6 (Nov. - Dec. 0), PP 4-0 Simple Oberver Baed Synchronization of Lorenz Sytem with Parametric Uncertainty Manih

More information

CDMA Signature Sequences with Low Peak-to-Average-Power Ratio via Alternating Projection

CDMA Signature Sequences with Low Peak-to-Average-Power Ratio via Alternating Projection CDMA Signature Sequence with Low Peak-to-Average-Power Ratio via Alternating Projection Joel A Tropp Int for Comp Engr and Sci (ICES) The Univerity of Texa at Autin 1 Univerity Station C0200 Autin, TX

More information

Lecture 10 Filtering: Applied Concepts

Lecture 10 Filtering: Applied Concepts Lecture Filtering: Applied Concept In the previou two lecture, you have learned about finite-impule-repone (FIR) and infinite-impule-repone (IIR) filter. In thee lecture, we introduced the concept of filtering

More information

General System of Nonconvex Variational Inequalities and Parallel Projection Method

General System of Nonconvex Variational Inequalities and Parallel Projection Method Mathematica Moravica Vol. 16-2 (2012), 79 87 General Sytem of Nonconvex Variational Inequalitie and Parallel Projection Method Balwant Singh Thakur and Suja Varghee Abtract. Uing the prox-regularity notion,

More information

New bounds for Morse clusters

New bounds for Morse clusters New bound for More cluter Tamá Vinkó Advanced Concept Team, European Space Agency, ESTEC Keplerlaan 1, 2201 AZ Noordwijk, The Netherland Tama.Vinko@ea.int and Arnold Neumaier Fakultät für Mathematik, Univerität

More information

Input-to-state stability for a class of Lurie systems

Input-to-state stability for a class of Lurie systems Automatica 38 (2002) 945 949 www.esevier.com/ocate/automatica Brief Paper Input-to-state stabiity for a cass of Lurie systems Murat Arcak a;, Andrew Tee b a Department of Eectrica, Computer and Systems

More information

Riser Dynamic Analysis Using WKB-Based Dynamic Stiffness Method

Riser Dynamic Analysis Using WKB-Based Dynamic Stiffness Method Rier ynamic Anayi Uing WKB-Baed ynamic Stiffne Method The MIT Facuty ha made thi artice openy avaiabe. Peae hare how thi acce benefit you. Your tory matter. Citation A Pubihed Pubiher Cheng, Yongming,

More information

IEOR 3106: Fall 2013, Professor Whitt Topics for Discussion: Tuesday, November 19 Alternating Renewal Processes and The Renewal Equation

IEOR 3106: Fall 2013, Professor Whitt Topics for Discussion: Tuesday, November 19 Alternating Renewal Processes and The Renewal Equation IEOR 316: Fall 213, Profeor Whitt Topic for Dicuion: Tueday, November 19 Alternating Renewal Procee and The Renewal Equation 1 Alternating Renewal Procee An alternating renewal proce alternate between

More information

THINKING IN PYRAMIDS

THINKING IN PYRAMIDS ECS 178 Course Notes THINKING IN PYRAMIDS Kenneth I. Joy Institute for Data Anaysis and Visuaization Department of Computer Science University of Caifornia, Davis Overview It is frequenty usefu to think

More information

arxiv:patt-sol/ v1 10 Jun 1997

arxiv:patt-sol/ v1 10 Jun 1997 Ampitude Equation for Eectrotatic Wave: mutipe pecie John David Crawford and Anandhan Jayaraman Department of Phyic and Atronomy Univerity of Pittburgh Pittburgh, Pennyvania 15260 (November 15, 2018) arxiv:patt-o/9706001v1

More information

c n b n 0. c k 0 x b n < 1 b k b n = 0. } of integers between 0 and b 1 such that x = b k. b k c k c k

c n b n 0. c k 0 x b n < 1 b k b n = 0. } of integers between 0 and b 1 such that x = b k. b k c k c k 1. Exitence Let x (0, 1). Define c k inductively. Suppoe c 1,..., c k 1 are already defined. We let c k be the leat integer uch that x k An eay proof by induction give that and for all k. Therefore c n

More information

Lecture 4. Chapter 11 Nise. Controller Design via Frequency Response. G. Hovland 2004

Lecture 4. Chapter 11 Nise. Controller Design via Frequency Response. G. Hovland 2004 METR4200 Advanced Control Lecture 4 Chapter Nie Controller Deign via Frequency Repone G. Hovland 2004 Deign Goal Tranient repone via imple gain adjutment Cacade compenator to improve teady-tate error Cacade

More information

PIPELINING AND PARALLEL PROCESSING. UNIT 4 Real time Signal Processing

PIPELINING AND PARALLEL PROCESSING. UNIT 4 Real time Signal Processing PIPELINING AND PARALLEL PROCESSING UNIT 4 Real time Signal Proceing Content Introduction Pipeling of FIR Digital Filter Parallel proceing Low power Deign FIR Digital Filter A FIR Filter i defined a follow:

More information

CONTROL SYSTEMS. Chapter 5 : Root Locus Diagram. GATE Objective & Numerical Type Solutions. The transfer function of a closed loop system is

CONTROL SYSTEMS. Chapter 5 : Root Locus Diagram. GATE Objective & Numerical Type Solutions. The transfer function of a closed loop system is CONTROL SYSTEMS Chapter 5 : Root Locu Diagram GATE Objective & Numerical Type Solution Quetion 1 [Work Book] [GATE EC 199 IISc-Bangalore : Mark] The tranfer function of a cloed loop ytem i T () where i

More information

696 Fu Jing-Li et al Vol. 12 form in generalized coordinate Q ffiq dt = 0 ( = 1; ;n): (3) For nonholonomic ytem, ffiq are not independent of

696 Fu Jing-Li et al Vol. 12 form in generalized coordinate Q  ffiq dt = 0 ( = 1; ;n): (3) For nonholonomic ytem, ffiq are not independent of Vol 12 No 7, July 2003 cfl 2003 Chin. Phy. Soc. 1009-1963/2003/12(07)/0695-05 Chinee Phyic and IOP Publihing Ltd Lie ymmetrie and conerved quantitie of controllable nonholonomic dynamical ytem Fu Jing-Li(ΛΠ±)

More information

7 Minimal realization and coprime fraction

7 Minimal realization and coprime fraction 7 Minimal realization and coprime fraction 7 Introduction If a tranfer function i realizable, what i the mallet poible dimenion? Realization with the mallet poible dimenion are called minimal-dimenional

More information

Social Studies 201 Notes for March 18, 2005

Social Studies 201 Notes for March 18, 2005 1 Social Studie 201 Note for March 18, 2005 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the

More information

Unbounded solutions of second order discrete BVPs on infinite intervals

Unbounded solutions of second order discrete BVPs on infinite intervals Available online at www.tjna.com J. Nonlinear Sci. Appl. 9 206), 357 369 Reearch Article Unbounded olution of econd order dicrete BVP on infinite interval Hairong Lian a,, Jingwu Li a, Ravi P Agarwal b

More information

Pacific Journal of Mathematics

Pacific Journal of Mathematics Pacific Journal of Mathematic OSCILLAION AND NONOSCILLAION OF FORCED SECOND ORDER DYNAMIC EQUAIONS MARIN BOHNER AND CHRISOPHER C. ISDELL Volume 230 No. March 2007 PACIFIC JOURNAL OF MAHEMAICS Vol. 230,

More information

Exercises for lectures 19 Polynomial methods

Exercises for lectures 19 Polynomial methods Exercie for lecture 19 Polynomial method Michael Šebek Automatic control 016 15-4-17 Diviion of polynomial with and without remainder Polynomial form a circle, but not a body. (Circle alo form integer,

More information

A Simple Approach to Synthesizing Naïve Quantized Control for Reference Tracking

A Simple Approach to Synthesizing Naïve Quantized Control for Reference Tracking A Simple Approach to Syntheizing Naïve Quantized Control for Reference Tracking SHIANG-HUA YU Department of Electrical Engineering National Sun Yat-Sen Univerity 70 Lien-Hai Road, Kaohiung 804 TAIAN Abtract:

More information

Control Systems Analysis and Design by the Root-Locus Method

Control Systems Analysis and Design by the Root-Locus Method 6 Control Sytem Analyi and Deign by the Root-Locu Method 6 1 INTRODUCTION The baic characteritic of the tranient repone of a cloed-loop ytem i cloely related to the location of the cloed-loop pole. If

More information

Assessment of Performance for Single Loop Control Systems

Assessment of Performance for Single Loop Control Systems Aement of Performance for Single Loop Control Sytem Hiao-Ping Huang and Jyh-Cheng Jeng Department of Chemical Engineering National Taiwan Univerity Taipei 1617, Taiwan Abtract Aement of performance in

More information

EFFECT ON PERSISTENCE OF INTRA-SPECIFIC COMPETITION IN COMPETITION MODELS

EFFECT ON PERSISTENCE OF INTRA-SPECIFIC COMPETITION IN COMPETITION MODELS Electronic Journal of Differential Equation, Vol. 2007(2007, No. 25, pp. 0. ISSN: 072-669. URL: http://ejde.math.txtate.edu or http://ejde.math.unt.edu ftp ejde.math.txtate.edu (login: ftp EFFECT ON PERSISTENCE

More information

Lecture 21. The Lovasz splitting-off lemma Topics in Combinatorial Optimization April 29th, 2004

Lecture 21. The Lovasz splitting-off lemma Topics in Combinatorial Optimization April 29th, 2004 18.997 Topic in Combinatorial Optimization April 29th, 2004 Lecture 21 Lecturer: Michel X. Goeman Scribe: Mohammad Mahdian 1 The Lovaz plitting-off lemma Lovaz plitting-off lemma tate the following. Theorem

More information