Z. GONG AND M. ALDEEN aso been shown that æxed modes can be eiminated by vibrationa contro and by samping techniques èsee ë1ë,ë1ë,ë17ëè. However, ther

Size: px
Start display at page:

Download "Z. GONG AND M. ALDEEN aso been shown that æxed modes can be eiminated by vibrationa contro and by samping techniques èsee ë1ë,ë1ë,ë17ëè. However, ther"

Transcription

1 Journa of Mathematica Systems, Estimation, and Contro Vo. 7, No. 1, 1997, pp. 1í1 cæ 1997 Birkhíauser-Boston Stabiization of Decentraized Contro Systems æ Zhiming Gong Mohammad Adeen Abstract The probem of stabiization of inear time-invariant systems under genera decentraized feedback schemes is considered in this paper. Anove approach to the probem is advised, in which the interactions between the strongy connected subsystems of a system are treated as disturbances. A necessary and suæcient condition for the existence of a decentraized controer which stabiizes a given system is presented. Keywords: decentraized contro, æxed modes, arge-scae systems, stabiization, quotient æxed modes AMS Subject Cassiæcations: 9A1, 9A15, 9B5, 9D15 1 Introduction For neary two decades, much attention has been paid to the probem of stabiization of decentraized inear time-invariant systems. Wang and Davison ë1ë introduced the notion of ëæxed mode" and obtained a soution to the decentraized contro probem. They concuded that the absence of æxed modes is a necessary and suæcient condition for the existence of a decentraized inear time-invariant system which stabiizes a given system èaso see ë1ë,ëë,ëë,ë1ëè. This resut has been the most important contribution to the decentraized contro theory so far. Anderson and Moore ëë considered the case where a decentraized time-varying controer is empoyed. They concuded that the absence of unstabe æxed modes of a system is not a necessary condition for the existence of a stabiizing decentraized time-varying controer, and in fact æxed modes may be eiminated by time-varying controers. Wang ë0ë observed independenty the property of time-varying controer in the eimination of æxed modes. It has æ Received September 1, 199; received in æna form November 0, Summary appeared in Voume 7, Number 1,

2 Z. GONG AND M. ALDEEN aso been shown that æxed modes can be eiminated by vibrationa contro and by samping techniques èsee ë1ë,ë1ë,ë17ëè. However, there has been some confusion in the iterature as to what kind of æxed modes can actuay be eiminated, see, e.g., ë17ë,ë18ë. Wiems ëë cariæed the question to certain extent and argued that ëstructuray æxed modes" caused by the fact that the system is not strongy connected are aso æxed modes with respect to time-varying output feedback. It shoud be noted that ëstructuray æxed modes" meant by Wiems ëë are not the same structuray æxed modes as studied in ë15ë and ë1ë, where structuray æxed modes were deæned through the notions of structured matrices and structuray equivaent systems. Very recenty, Khargonekar and Ozguer ë11ë further cariæed the question by using ëifting" technique for periodic systems. They showed that a æxed modes, except those associated with the compementary subsystems èëë,ë1ëè having zero transfer function matrices, can be eiminated by a periodicay time-varying decentraized controer. Then, the foowing genera question may be asked. What kind of æxed modes can, or cannot, be eiminated by aneven arger cass of controers than periodicay time-varying decentraized controers, say genera timevarying, noninear, vibrationa, or any other kind of suitabe controers? To provide an answer to the question given above, we consider in this paper genera decentraized feedback schemes where no restrictions are imposed on the type and structure of the decentraized controers, as ong as the constraint of the decentraized information structure is satisæed. A necessary and suæcient condition for the existence of a decentraized controer which stabiizes a given system is presented in this paper. The condition is stated in a neat term of æxed modes of a quotient system, which were ærst identiæed in ë9ë and ë10ë and are termed Quotient Fixed Modes in this paper. It is shown that the quotient æxed modes are exacty those æxed modes associated with the compementary subsystems having zero transfer function matrices, which were identiæed by Khargonekar and Ozguer ë11ë. The resut of this paper provides a precise statement of the arguments in ëë and ë11ë and reveas the further fact that if a æxed modes cannot be eiminated by a decentraized periodicay time-varying controer, then there is no way this æxed mode can be eiminated by a decentraized controer, no matter what kind of controers are used. This provides a compete answer to the question of what kind of æxed modes can, or cannot, be eiminated. The ayout of the paper is as foows. Section introduces concept of time-varying vectors with a degree of exponentia stabiity and gives forma statement of the probem. The main resut of this paper is presented in Section and a proof of the main resut is given in Section. Section 5 is the concusion.

3 DECENTRALIZED CONTROL SYSTEMS Statement of the probem Consider a decentraized contro system described by _xètè =Axètè+ NX i=1 B i u i ètè è:1aè y i ètè =C i xètè; i=1;;æææ ;N è:1bè where N is the number of the contro stations, xètè R n is the state vector of the system, u i ètè R mi and y i ètè R ri are, respectivey, the input and the output vectors of the ith contro station, and A, B i and C i are rea constant matrices of appropriate dimensions. Assume that the information avaiabe to the ith contro station at time t is represented by I i ètè =fy i èè;u i èè:ë0;të; ë0;tèg: è:è It is aso assumed that the constraint of the decentraized information structure is such that the contro input u i ètè at the ith oca contro station can ony be cacuated from I i ètè, i.e., u i ètè =F i èi i ètè;tè; i=1;;æææ ;N è:è where F i èi i ètè;tè denotes a function of I i ètè and time t. Obviousy, the cass of the controers è.è incudes the decentraized inear time-invariant, time-varying, or even noninear controers. In fact, it is the argest cass of decentraized controers for the system è.1è. To state the probem studied in this paper, the foowing deænition is required. Deænition 1 A time-varying vector vètè is said to be stabe with a degree of exponentia stabiity èdesè æ, if there exists a bounded scaar function fèæè so that kvètèk fèt 0 èe,æèt,t0è ; 8t 0 0;tt 0 è:è where æ is a positive rea number and kvètèk denotes the Eucidean norm of vètè. In addition, if fèt 0 è=akvètèk; 8t 0 0 è:5è where a is a constant, then, the time-varying vector vètè is said to be uniformy stabe with a DES æ. Under Deænition 1, stabe time-varying vectors have the foowing usefu properties. Let G 1 ètè and G ètè denote two bounded function matrices with

4 Z. GONG AND M. ALDEEN appropriate dimensions. If time-varying vectors v 1 ètè and v ètè are stabe with DES æ 1 and æ, respectivey, then, èiè the vector G 1 ètèv 1 ètè is stabe with a DES æ 1 ; and èiiè the vector G 1 ètèv 1 ètè+g ètèv ètè is stabe with a DES æ=minfæ 1, æ g; and èiiiè the vector G1 ètèv 1 ètè G ètèv ètè is stabe with a DES æ=minfæ 1, æ g. In the seque, a system is said to be stabe with a DES æ if the state vector of this system, with zero externa input, is stabe with a DES æ for any possibe initia state. When the DES is not concerned, the system is said to be stabe for any positive rea numbers æ. The probem of stabiization of decentraized inear time-invariant systems can now be stated as foows: Find a decentraized controer è.è for the system è.1è so that the resuting cosed-oop system is stabe, or stabe with a prescribed DES. Main Resut In order to state the main resut of this paper, we need to introduce the concept of the digraphs and strongy connectedness of decentraized contro systems ëë,ë9ë. Consider the N-channe decentraized contro system è.1è, the digraph of this system is deæned as a set of N nodes and some directed arcs connecting these nodes. The nodes represent the contro stations of the system and directed arcs represent the connections between them. If C j èsi, Aè,1 B i = 0, then there exist a directed arc from node i to node j èi; j = 1; ; æææ ; Nè. If for each distinct pair of i and j, the digraph contains directed paths from node i to node j and from node j to node i, then the digraph is said to be strongy connected, and so is caed the system è.1è. If a digrah is not strongy connected, it can aways be decomposed uniquey into a number of strongy connected components, which is the argest strongy connected sub-digraphs, in the sense that whenever an extra node is added to such a strongy connected sub-digraph, it wi be no onger strongy connected. Corresponding to the decomposition of the digraph into a number of strongy connected components, the system equations è.1è can be written

5 DECENTRALIZED CONTROL SYSTEMS 1 JQQk JJ QQs J J^ æ ææ æ æ Q QQQQs QQs 5 QQk Figure 1: Digraph of a Decentraized Contro System as _xètè =Axètè+ XN æ i=1 B æ i uæ i ètè è:1aè y æ i ètè =C æ i xètè; i=1;;æææ ;Næ è:1bè where N æ denotes the number of the strongy connected components of the system è.1è, u æ i ètè Rmæ i and yi æètè Rræ i denote the input and the output vectors consisting of a of the inputs and the outputs in the ith strongy connected component, respectivey, and Bi æ and Ci æ are the corresponding input and output matrices. Suppose that the N contro stations of the system è.1è are aggregated into the N æ contro stations as described by è.1è. Then the decentraized contro system è.1è with N æ contro stations is caed a quotient system of the system è.1è ë9ë. The foowing exampe iustrates the concept of quotient systems. Let i! j stands for the condition C j èsi, Aè,1 B i =0. Suppose that a decentraized contro system has 5 oca contro stations and ony the foowing conditions hod: 1!, 1!,!,! 1,!,! 5, 5!. Then, the digraph of the system is as shown in Figure 1. It is easy to see that this digraph has strongy connected components. Therefore, its quotient system has aggregated oca contro stations. For comparison, signa æow graphs of the origina decentraized contro system and its quotient system are iustrated in Figure and Figure, respectivey. Let æ æ denote the set of the æxed modes ë1ë of the quotient system è.1è, i.e. æ æ = ë K æ i Rmæ i æræ i è A + XN æ i=1 B æ i Kæ i Cæ i! è:è where èæè denotes the set of the eigenvaues of the argument matrix and the intersection takes over a of the matrices K æ i Rmæ i æræ i, i =1;;æææ ;N æ. Deænition The æxed modes of the quotient system è.1è are caed quotient æxed modes of the system è.1è. 5

6 Z. GONG AND M. ALDEEN 1 QQs QQk Q J æ æ JJ QQQQs QQs 5 J J^ æ ææ QQk???? Controer Controer Controer Controer 1 5? Controer Figure : A Decentraized Contro System 1 QQs Q JQQk JJ æ æ QQQQs QQs 5 J J^ æ ææ QQk????? Controer 1 Controer Figure : A Quotient System

7 DECENTRALIZED CONTROL SYSTEMS It is cear that the set of the quotient æxed modes of the system è.1è is a subset of the æxed modes of this system, i.e. æ æ æ, where æ denotes the set of the æxed modes of the system è.1è and is given by ë1ë è! ë NX æ= A+ B i K i C i : è:è K ir m i ær i The main resut of this paper is now stated in the foowing theorem. Theorem 1 Given the system è.1è, there exists a decentraized controer è.è for this system so that èiè the cosed-oop system is stabe with any prescribed DES if and ony if the system è.1è has no quotient æxed modes, i.e. æ æ = ;; and èiiè the cosed-oop system is stabe if and ony if the system è.1è has no unstabe quotient æxed modes, i.e. æ æ C,, where C, denotes the eft haf open region of the compex pane. There exists an obvious anaogy between Theorem 1 above and the ceebrated resut by Wang and Davison ë1, Theorem 1ë, which may be stated as foows: Given the system è.1è, there exists a decentraized inear time-invariant èdynamicè controer for this system so that the cosed-oop system is stabe if and ony if the system è.1è has no unstabe æxed modes, i.e. æ C,. Comparing with Wang and Davison's resut ë1, Theorem 1ë with Theorem 1 of this paper, it is cear that when the considered famiy of decentraized controers for the system è.1è is enarged from inear time-invariant ones to the most genera ones è.è, the set of æxed modes æ in Wang and Davison's resut ë1, Theorem 1ë shoud be repaced by the set of quotient æxed modes æ æ, an subset of æ. Using ëifting" technique, Khargonekar and Ozguer ë11, Theorem ë showed recenty that a inear time-invariant system is stabiizabe by a decentraized periodicay time-varying controer if and ony if there exists no unstabe ëincompeting zeros" in the compementary subsystems which have zero transfer function matrices. Their resut gives a precise statement of the argument caimed by Wiems ëë. As shown ate in this paper, the incompeting zero of the compementary subsystems having zero transfer matrices are exacty the quotient æxed modes deæned in this paper. Therefore, Theorem 1 given above extends the resuts of Wiems ëë and Kargonekar and Ozguer ë11ë, and reveas the further fact that if a decentraized inear time-invariant system cannot be stabiized by a decentraized periodicay time-varying controer, then, no matter what kind of structures of the controers are empoyed, there exists no decentraized controers which can stabiize the system. In the neat term of quotient 7 i=1

8 Z. GONG AND M. ALDEEN æxed modes, Theorem 1 of this paper gives a compete answer to the question on what kind of æxed modes can, and cannot, be eiminated under the constraint of decentraized contro. Proof of the Main Resut In this section, a proof of Theorem 1 is presented. We ærst introduce a emma. Lemma 1 Consider a decentraized contro system, which is subject to unmeasurabe disturbances and otherwise identica to the system è.1è, described by _xètè=axètè+ NX i=1 B i u i ètè+dètè è:1aè y i ètè =C i xètè+g i ètè; i=1;;æææ ;N è:1bè where dètè and g i ètè èi=1;;æææ ;Nè are disturbances and are stabe with a DES æ. If the system è.1è is strongy connected, then, there exists a decentraized controer è.è for the system è.1è so that èiè the state vector of the cosed-oop system is stabe with a prescribed DES æ æ if the system è.1è is centraized controabe and observabe; and èiiè the state vector of the cosed-oop system is stabe if the system è.1è has no unstabe centraized uncontroabe andèor unobservabe modes. Proof of Lemma 1: See the Appendix. Next, we introduce some preiminary deveopment of the proof of Theorem 1. From the deænitions of digraph and quotient systems, it foows that the transfer function matrix of the system è.1è is of bock trianguar structure when the contro stations of the system are ordered appropriatey. Assume that the contro stations of the system è.1è have been ordered appropriatey and the state vector has aso been chosen appropriatey. Then the matrices in the system equations è.1è and è.1è can assume the foowing trianguar structure ëë: A = ~A 1 æ æææ æ 0 A ~ æææ æ ; 0 0 æææ AN ~ æ 8

9 DECENTRALIZED CONTROL SYSTEMS B =ëb æ 1 Bæ æææbæ N æë= C = C æ 1 C æ.. C æ N æ 7 5 = ~B 1 æ æææ æ 0 B ~ æææ æ ; 0 0 æææ BN ~ æ ~C 1 æ æææ æ 0 C ~ æææ æ æææ CN ~ æ 7 5 è:è where C = ëc 0 1 C0 æææc 0 N ë0, B = ëb 1 B æææb N ë and ~ Ai, ~ Bi and ~ Ci èi = 1; ; æææ ;N æ è are nonzero matrices of appropriate dimensions. Based on the structure of the matrices è.è, the system equations è.1è and è.1è can equivaenty be written as _x i ètè = ~ Ai x i ètè+ ~ Bi u æ iètè+w i ètè è:aè y æ i ètè = ~ Ci x i ètè+v i ètè; i=1;;æææ ;N æ è:bè where x i ètè, w i ètè and v i ètè are vectors determined correspondingy by è.è and the system equations è.1è and è.1è. From equations è.è, it is cear that the case where the system è.1è is controed by a decentraized controer is equivaent to the case where each subsystem in è.è is controed by an individua decentraized controer. Aso, it can be seen from è.è that, for a certain ith subsystem, w i ètè and v i ètè wi not contain information about x i ètè, regardess of the structure and type of the decentraized controer. In other words, w i ètè and v i ètè can be considered as exogenic disturbances to their respective subsystems. Therefore, to anayze the decentraized stabiizabiity of the goba system, it suæces to anayze individuay the decentraized stabiizabiity of each subsystem in è.è, in the presence of the disturbances. Furthermore, the system è.1è is stabe or stabe with a DES æ if and ony if the state of each subsystem è.è is stabe or stabe with DES æ i æ, respectivey. It is easy to show that the subsystems described by the tripes è Ci ~, Ai ~, ~B i è, i =1;;æææ ;N æ, are strongy connected, with respect to the contro stations in their respective strongy connected components. These subsystems are therefore caed strongy connected subsystems of the system è.1è ë5ë. Let æ c i denote the set of centraized uncontroabe andèor unobservabe modes of the strongy connected subsystem è Ci ~ ; Ai ~ ; Bi ~ è. Then, the foowing hods ë1, Lemma 1ë. æ c i = ë K æ i Rm i ær i ~Ai + ~ Bi K æ i ~ C i ; i =1;;æææ ;N æ : è:è 9

10 Z. GONG AND M. ALDEEN Using è.è and è.è, it is not diæcut to show that æ æ = Në æ i=1 æ c i : è:5è Equation è.5è states that the set of quotient æxed modes of the system è.1è is just the union of the centraized uncontroabe andèor unobservabe modes of a of the strongy connected subsystems of this system. Now, we may appy the we-known Kaman decomposition ë8ë to the strongy connected subsystems è ~ Ci ; ~ Ai ; ~ Bi è and present straightforward proof of Theorem 1. Proof of the necessity part of Theorem 1: If the system è.1è has a quotient æxed mode,, then by è.5è, there is at east one strongy connected subsystem of the system è.1è, say, the subsystem è ~ Cp ; ~ Ap ; ~ Bp è, which has a centraized uncontroabe andèor unobservabe mode. Consider now the pth subsystem è.è corresponding to the strongy connected subsystem è ~ Cp ; ~ Ap ; ~ Bp è. When a decentraized controer èeven a centraized controerè for this subsystem is appied, the state response of the resuting cosed-oop system wi contain terms proportiona to e t when the initia state corresponding to the mode is nonzero. According to Kaman ë8ë, this statement is true regardess of the structure and type of the appied controer, i.e., regardess of whether the controer is inear or noninear, dynamic or static, time-varying or time-invariant, or any other kind of controer. This means that the state of the pth subsystem è.è can not be stabiized with an arbitrariy prescribed DES. By the same reason, if such aexists and is unstabe è = C, è, then, that subsystem can not be stabiized by any controer. This competes the proof of the necessity part of both statements èiè and èiiè in Theorem 1. Proof of the suæciency part of Theorem 1: For any prescribed positive rea number æ, et æ i, i =1;;æææ ;N æ, be numbers such that æ æ 1 æ æææ æ N æ. In the foowing, it is shown that if æ æ = ;, then there exists individua decentraized controers for each subsystem è.è so that the states of the resuting cosed-oop systems are stabe with DES æ 1 ;æ ;æææ ;æ N æ, respectivey. If æ æ = ;, then, by è.5è, each strongy connected subsystem of the system è.1è is centraized controabe and observabe. Consider ærst the N æ -th subsystem in è.è where w N æètè =v N æètè=0.according to Lemma 1, there exists a decentraized controer for this subsystem so that the state x N æètè of the resuting cosed-oop system is stabe with a DES æ N æ. Then, by noticing the matrix structure in è.è and using the properties of stabe time-invariant vectors, it foows that the disturbances w N æ,1ètè and the v Næ,1 ètè to èn æ, 1è-th subsystem in è.è are stabe with a DES æ N æ. Using Lemma 1 again, it foows that there exists a decentraized controer 10

11 DECENTRALIZED CONTROL SYSTEMS for the èn æ, 1è-th subsystem in è.è so that the state of the resuting cosed-oop system is stabe with a DES æ N æ,1. Repeating the above arguments for the rest of the subsystems in è.è, i.e. the subsystems in è.è with i = N æ,;n æ,;æææ ;;1;it is concuded that there exists N æ decentraized controers for each subsystem in è.è, so that the states of the resuting cosed-oop systems are stabe with DES æ N æ;æ N æ,1;æ N æ,;æææ ;æ ;æ 1 ; respectivey. This impies that there exists a decentraized controer è.è for system è.1è so that the resuting cosedoop system is stabe with the prescribed DES æ. This competes the proof of the suæciency part of statement èiè in Theorem 1. If the system è.1è has no unstabe quotient æxed modes, then, by è.5è, each strongy connected subsystem of the system è.1è has no unstabe centraized uncontroabe andèor unobservabe modes. Then, by using Kaman's decomposition methods ë8ë and arguments simiar to the above, it is not diæcut to show that there exists a decentraized controer for the system è.1è so that the resuting cosed-oop system is stabe. This competes the proof of the suæciency part of statement èiiè in Theorem 1, and consequenty the proof of Theorem 1. Khargonekar and Ozguar ë11ë showed that a inear time-invariant system can be stabiized by a decentraized periodicay time-varying controer if and ony if there exist no unstabe ëincompeting zero" in the compementary subsystems having zero transfer function matrices. It is shown in the foowing that these ëincompeting zero" are exacty the unstabe quotient æxed modes deæned in this paper. Consider the decentraized contro system è.1è with the matrices A, B and C in the form given by è.è. Then, according to the deænition of strongy connected digraphs, a the compementary subsystems with zero transfer function matrices have the form è ^Cj ;A; ^Bj è, j f1; ; æææ ;N æ g, with ^B j = æ B æ 1 Bæ æææ ; Bæ jæ and ^Cj = C æ j+1 C æ j+.; C æ N æ 7 5 è:è where Bi æ èi =1;;æææ ;jè and Ci æ èi = j +1;j +;æææ ;N æ è are matrices given in è.è. Consider now the system matrix I, A S j = ^Bj, ^Cj 0 : è:7è By the deænition given in ë11ë, a number is an incompeting zero of the compementary subsystems è ^Cj ;A; ^Bj è iæ the rank of the matrix S j is ess 11

12 Z. GONG AND M. ALDEEN than n, the dimension of the matrix A. Due to the trianguar structure of the matrices A, B and C as shown in è.è, we have where RankS j = Rank M 1 æ æææ æ 0 M æææ æ æææ M N æ æ M i = I, A ~ æ i Bi ~ ; i =1;;æææ ;j I, Ai ~ M i =, Ci ~ ; i = j +1;j+;æææ ;N æ : 7 5 è:8è Using a resut by Gong and Adeen ë5, Lemma A1ë, it can be shown easiy that the rank of the matrix S j is ess than n if and ony if there exists a matrix M i èi f1; ; æææ ;N æ gè which doesn't have fu rank. Then, by using the we-known rank tests for controabiity and observabiity of the system è Ci ~ ; Ai ~ ; Bi ~ è and the resut of è.5è, it is concuded that the set of the incompeting zeros of the compementary subsystems with zero transfer function matrices are exacty the set of the quotient æxed modes of the system è.1è. 5 Concusion The probem of stabiization of inear time-invariant systems under genera decentraized feedback schemes is investigated in this paper. By using the notion of quotient æxed modes, the paper presented a necessary and suæcient condition for the existence of a decentraized controer to stabiize a given inear time-invariant system ètheorem 1è. This resut gives a compete answer to the question on what kind of æxed modes can actuay be eiminated under the constraint of decentraized contro. The argument by Wiems ëë is now extended to the foowing: ëstructure æxed modes" caused by the fact that the system is not strongy connected are aso æxed modes with respect to decentraized time-varying output feedback controers; and in fact they are æxed modes with respect to a kind of decentraized controers, incuding noninear ones. Appendix Proof of Lemma 1. In this appendix, a proof of Lemma 1 is provided. Assume that the system è.1è is centraized controabe and observabe, and strongy connected. According to Anderson and Moore ëë, there exists a decentraized 1

13 DECENTRALIZED CONTROL SYSTEMS time-varying feedback controer, described by u i ètè = ^Ki ètèy i ètè+^u i ètè; i=1;;æææ ;N èa:1è where ^K are bounded matrices, for system è.1è so that the resuting timevarying cosed-oop system, described by _xètè = ^Axètè+ N X i=1 B i^u i ètè èa:aè y i ètè =C i xètè; i=1;;æææ ;N èa:bè where ^A = A + P N i=1 B i ^K i ètèc i, is uniformy competey controabe and observabe by an arbitrariy predetermined contro station of the system, say by contro station 1. Then, according to Ikeda et a. ë7ë, there exists a controer at contro station 1 for the system èa.è so that the cosed-oop system is uniformy stabe with any prescribed DES. Such a controer may be described by ^u 1 ètè =Fètèzètè èa:aè ^u i ètè =0; i=;;æææ ;N èa:bè where F ètè is a bounded feedback gain matrix and zètè is the estimate of the state vector xètè obtained by the observer _zètè =è^aètè,hètèc1 èzètè+b 1^u 1 ètè+hètèy 1 ètè where Hètè is a bounded observer gain matrix. Then the cosed-oop system is given by _x c ètè =A c ètèx c ètè èa:è where x c ètè = xètè zètè ; A c ètè= ^Aètè B1 Fètè HètèC 1 ^Aètè,HètèC1 +B 1 Fètè Now, appy the decentraized controer as described by èa.1è and èa.è to the system è.1è. This eads to the foowing cosed-oop system. where w c ètè = _x c ètè =A c ètèx c ètè+w c ètè dètè+ P N i=1 B i ^K i ètèg i ètè Hètèg 1 ètè : : èa:5è As dètè and g i ètè èi =1;;æææ ;Nè are a stabe with a DES æ, then, by noticing the property of stabe time-varying vectors, it foows that w c ètè is stabe with a DES æ. 1

14 Z. GONG AND M. ALDEEN For a given number æ æ, et æ be a number such that æ = æ and æ = minfæ; æg. Design the controer èa.è so that the system èa.è is uniformy stabe with a DES æ. Let æèæ; æè denote the transition matrix associated with the system èa.è. Since the system èa.è is uniformy stabe with a DES æ, it is easy to show that kæèt; t 0 èkae,æèt,t0è ; 8t 0 ;tt 0 èa:è where a is a constant. The state response of the system èa.5è is given by x c ètè =æèt; t 0 èx c èt 0 è+ Z t As w c ètè is stabe with a DES æ, i.e. t 0 æèt; èw c èèd; 8t 0 ;tt 0 : èa:7è kw c ètèk f 1 èt 0 èe,æèt,t0è ; 8t 0 ;tt 0 èa:8è where f 1 èæè is a bounded function, it foows that kx c ètèk akx c èt 0 èke,æèt,t0è +f 1 èt 0 è = akx c èt 0 èke,æèt,t0è + af 1èt 0 è æ, æ f èt 0 èe,æèt,t0è Z t where f èæè is a bounded function given by f èt 0 è=akx c èt 0 èk+ af 1èt 0 è jæ, æj : t 0 e,æèt,è e,æè,t0è d ne,æèt,t0è, e,æèt,t0èo èa:9è Equation èa.9è impies that the system èa.5è, and therefore the system è.1è, is stabe with the given DES æ. This competes the proof of part èiè of Lemma 1. For the case where the system è.1è is not centraized controabe and observabe, assume that the system has no unstabe centraized uncontroabe andèor unobservabe modes and strongy connected. Then, by using Kaman's canonica system decomposition form ë8ë and the simiar arguments to the above, it is not diæcut to show that there exists a decentraized controer which stabiizes the system è.1è. This competes the proof of Lemma 1. References ë1ë B.D.O. Anderson and D.J. Cements. Agebraic characterization of æxed modes in decentraized contro, Automatia 17 è1981è,

15 DECENTRALIZED CONTROL SYSTEMS ëë ëë ëë ë5ë ëë ë7ë ë8ë ë9ë B.D.O. Anderson and J.B. Moore. Time-varying feedback aws for decentraized contro, IEEE Trans. Automatic Contro è1981è, J.P. Corfmat and A.S. Morse. Decentraized contro of inear mutivariabe systems, Automatica 1 è197è, E.J. Davison and T.N. Chang. Decentraized stabiization and poe assignment for genera proper systems, IEEE Trans. Automat. Contr. 5 è1990è, 5-. Z. Gong and M. Adeen. On the characterization of æxed modes in decentraized contro, IEEE Trans. Automatic Contro 7 è199è, A.N. Gundes and C.A. Desoer. Agebraic Theory of Linear Feedback Systems with Fu and Decentraized Compensators. Berin: Springer- Verag, M. Ikeda, H. Maeda and S. Kodama. Estimation and feedback in inear time-varying systems: a deterministic theory, SIAM J. Contro 1 è1975è, 0-. R.E. Kaman. Mathematica description of inear dynamica systems, SIAM J. Contro 1 è19è, H. Kobayashi and T. Yoshikawa. Graph-theoretic approach to controabiity and ocaizabiity of decentraized contro systems, IEEE Trans. Automatic Contro 7 è198è, ë10ë H. Kobayashi, H. Hanafusa and T. Yoshikawa. Controabiity under decentraized information structure, IEEE Trans. Automatic Contro è1978è, ë11ë P.P. Khargonekar and A.B. Ozguer. Decentraized contro and periodic feedback, IEEE Trans. Automat. Contr. 9 è199è, ë1ë S. Lee, M. Meerkov and T. Runofsson. Vibrationa Feedback Contro: Zero pacement capabiities, IEEE Trans. Automat. Contr. è1987è, ë1ë A.B. Ozguer. Decentraized contro: A stabe proper fractiona approach, IEEE Trans. Automat. Contr. 5 è1990è, ë1ë U. Ozguner and E.J. Davison. Samping and decentraized æxed modes, Proc American Contro Conference è1985è,

16 Z. GONG AND M. ALDEEN ë15ë K. Reinschke. Graph-theoretic characterization of æxed modes in centraized and decentraized contro, Int. J. Contro 9 è198è, ë1ë M.E. Sezer and D.D. Sijak. Structuray æxed modes, Systems and Contro Letters 1 è1981è, 0-. ë17ë L. Trave, A. Tarras and A. Titi. An appication of vibrationa contro to cance unstabe decentraized æxed modes, IEEE Trans. Automat. Contr. 0 è1985è, 8-8. ë18ë L. Trave, A. Titi and A. Tarras. Large Scae Systems: Decentraization, Structure Constraints and Fixed Modes. Berin: Springer-Verag, ë19ë K.A. Unyeiogu, A.B. Ozguer and P.P. Khargonekar. Decentraized simutaneous stabisation and reiabe contro using periodic feedback, Systems and Contro Letters 18 è199è, -1. ë0ë S.H. Wang. An exampe in decentraized contro systems, IEEE Trans. Automat. Contr. è1989è, ë1ë S.H. Wang and E.J. Davison. On the stabiization of decentraized contro systems, IEEE Trans. Automatic Contro 18 è197è, ëë J.L. Wiems. Time-varying feedback for the stabiization of æxed modes in decentraization contro systems, Automatica 5 è1989è, Schoo of Eectrica and Eectronic Engineering, Nanyang Technoogica University, Nanyang Avenue, Singapore 9798 Department of Eectrica and Eectronic Engineering, The University of Mebourne, Parkvie, Victoria, Austraia 05 Communicated by Cyde F. Martin 1

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

SE-514 (OPTIMAL CONTROL) OPTIMAL CONTROL FOR SINGLE AND DOUBLE INVERTED PENDULUM. DONE BY: Fatai Olalekan ( Ayman Abdallah (973610)

SE-514 (OPTIMAL CONTROL) OPTIMAL CONTROL FOR SINGLE AND DOUBLE INVERTED PENDULUM. DONE BY: Fatai Olalekan ( Ayman Abdallah (973610) SE-54 (OPTIAL CONTROL OPTIAL CONTROL FOR SINGLE AND DOUBLE INVERTED PENDULU DONE BY: Fatai Oaekan (363 Ayman Abdaah (9736 PREPARED FOR: Dr. Sami E-Ferik Tabe of contents Abstract... 3 Introduction... 3

More information

A.V. SAVKIN AND I.R. PETERSEN uncertain systems in which the uncertainty satisæes a certain integral quadratic constraint; e.g., see ë5, 6, 7ë. The ad

A.V. SAVKIN AND I.R. PETERSEN uncertain systems in which the uncertainty satisæes a certain integral quadratic constraint; e.g., see ë5, 6, 7ë. The ad Journal of Mathematical Systems, Estimation, and Control Vol. 6, No. 3, 1996, pp. 1í14 cæ 1996 Birkhíauser-Boston Robust H 1 Control of Uncertain Systems with Structured Uncertainty æ Andrey V. Savkin

More information

Available online at ScienceDirect. IFAC PapersOnLine 50-1 (2017)

Available online at   ScienceDirect. IFAC PapersOnLine 50-1 (2017) Avaiabe onine at www.sciencedirect.com ScienceDirect IFAC PapersOnLine 50-1 (2017 3412 3417 Stabiization of discrete-time switched inear systems: Lyapunov-Metzer inequaities versus S-procedure characterizations

More information

Symbolic models for nonlinear control systems using approximate bisimulation

Symbolic models for nonlinear control systems using approximate bisimulation Symboic modes for noninear contro systems using approximate bisimuation Giordano Poa, Antoine Girard and Pauo Tabuada Abstract Contro systems are usuay modeed by differentia equations describing how physica

More information

THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES

THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES by Michae Neumann Department of Mathematics, University of Connecticut, Storrs, CT 06269 3009 and Ronad J. Stern Department of Mathematics, Concordia

More information

Research Article Numerical Range of Two Operators in Semi-Inner Product Spaces

Research Article Numerical Range of Two Operators in Semi-Inner Product Spaces Abstract and Appied Anaysis Voume 01, Artice ID 846396, 13 pages doi:10.1155/01/846396 Research Artice Numerica Range of Two Operators in Semi-Inner Product Spaces N. K. Sahu, 1 C. Nahak, 1 and S. Nanda

More information

STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION

STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION Journa of Sound and Vibration (996) 98(5), 643 65 STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM G. ERDOS AND T. SINGH Department of Mechanica and Aerospace Engineering, SUNY at Buffao,

More information

NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS

NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS TONY ALLEN, EMILY GEBHARDT, AND ADAM KLUBALL 3 ADVISOR: DR. TIFFANY KOLBA 4 Abstract. The phenomenon of noise-induced stabiization occurs

More information

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s

More information

arxiv: v1 [cs.sy] 31 May 2014

arxiv: v1 [cs.sy] 31 May 2014 Cooperative Contro of Linear Muti-Agent Systems via Distributed Output Reguation and Transient Synchronization Georg S. Seyboth, Wei Ren, Fran Agöwer arxiv:46.85v cs.sy 3 May 24 Abstract A wide range of

More information

Convergence Property of the Iri-Imai Algorithm for Some Smooth Convex Programming Problems

Convergence Property of the Iri-Imai Algorithm for Some Smooth Convex Programming Problems Convergence Property of the Iri-Imai Agorithm for Some Smooth Convex Programming Probems S. Zhang Communicated by Z.Q. Luo Assistant Professor, Department of Econometrics, University of Groningen, Groningen,

More information

XSAT of linear CNF formulas

XSAT of linear CNF formulas XSAT of inear CN formuas Bernd R. Schuh Dr. Bernd Schuh, D-50968 Kön, Germany; bernd.schuh@netcoogne.de eywords: compexity, XSAT, exact inear formua, -reguarity, -uniformity, NPcompeteness Abstract. Open

More information

Mode in Output Participation Factors for Linear Systems

Mode in Output Participation Factors for Linear Systems 2010 American ontro onference Marriott Waterfront, Batimore, MD, USA June 30-Juy 02, 2010 WeB05.5 Mode in Output Participation Factors for Linear Systems Li Sheng, yad H. Abed, Munther A. Hassouneh, Huizhong

More information

DIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM

DIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM DIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM MIKAEL NILSSON, MATTIAS DAHL AND INGVAR CLAESSON Bekinge Institute of Technoogy Department of Teecommunications and Signa Processing

More information

Observability of linear systems under adversarial attacks*

Observability of linear systems under adversarial attacks* Observabiity of inear systems under adversaria attacks* Michee S. Chong 1, Masashi Wakaiki 1 and João P. Hespanha 1 Abstract We address the probem of state estimation for muti-output continuous-time inear

More information

Conditions for Saddle-Point Equilibria in Output-Feedback MPC with MHE

Conditions for Saddle-Point Equilibria in Output-Feedback MPC with MHE Conditions for Sadde-Point Equiibria in Output-Feedback MPC with MHE David A. Copp and João P. Hespanha Abstract A new method for soving output-feedback mode predictive contro (MPC) and moving horizon

More information

Positive real and bounded real balancing for model reduction of descriptor systems

Positive real and bounded real balancing for model reduction of descriptor systems TECHNISCHE UNIVERSITÄT BERLIN Positive rea and bounded rea baancing for mode reduction of descriptor systems Timo Reis Tatjana Styke Preprint 28/25 Preprint-Reihe des Instituts für Mathematik Technische

More information

<C 2 2. λ 2 l. λ 1 l 1 < C 1

<C 2 2. λ 2 l. λ 1 l 1 < C 1 Teecommunication Network Contro and Management (EE E694) Prof. A. A. Lazar Notes for the ecture of 7/Feb/95 by Huayan Wang (this document was ast LaT E X-ed on May 9,995) Queueing Primer for Muticass Optima

More information

Section 6: Magnetostatics

Section 6: Magnetostatics agnetic fieds in matter Section 6: agnetostatics In the previous sections we assumed that the current density J is a known function of coordinates. In the presence of matter this is not aways true. The

More information

Systems & Control Letters. State and output feedback boundary control for a coupled PDE ODE system

Systems & Control Letters. State and output feedback boundary control for a coupled PDE ODE system Systems & Contro Letters 6 () 54 545 Contents ists avaiabe at ScienceDirect Systems & Contro Letters journa homepage: wwweseviercom/ocate/syscone State output feedback boundary contro for a couped PDE

More information

BASIC NOTIONS AND RESULTS IN TOPOLOGY. 1. Metric spaces. Sets with finite diameter are called bounded sets. For x X and r > 0 the set

BASIC NOTIONS AND RESULTS IN TOPOLOGY. 1. Metric spaces. Sets with finite diameter are called bounded sets. For x X and r > 0 the set BASIC NOTIONS AND RESULTS IN TOPOLOGY 1. Metric spaces A metric on a set X is a map d : X X R + with the properties: d(x, y) 0 and d(x, y) = 0 x = y, d(x, y) = d(y, x), d(x, y) d(x, z) + d(z, y), for a

More information

Cryptanalysis of PKP: A New Approach

Cryptanalysis of PKP: A New Approach Cryptanaysis of PKP: A New Approach Éiane Jaumes and Antoine Joux DCSSI 18, rue du Dr. Zamenhoff F-92131 Issy-es-Mx Cedex France eiane.jaumes@wanadoo.fr Antoine.Joux@ens.fr Abstract. Quite recenty, in

More information

A Brief Introduction to Markov Chains and Hidden Markov Models

A Brief Introduction to Markov Chains and Hidden Markov Models A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,

More information

An implicit Jacobi-like method for computing generalized hyperbolic SVD

An implicit Jacobi-like method for computing generalized hyperbolic SVD Linear Agebra and its Appications 358 (2003) 293 307 wwweseviercom/ocate/aa An impicit Jacobi-ike method for computing generaized hyperboic SVD Adam W Bojanczyk Schoo of Eectrica and Computer Engineering

More information

The Group Structure on a Smooth Tropical Cubic

The Group Structure on a Smooth Tropical Cubic The Group Structure on a Smooth Tropica Cubic Ethan Lake Apri 20, 2015 Abstract Just as in in cassica agebraic geometry, it is possibe to define a group aw on a smooth tropica cubic curve. In this note,

More information

Control with Minimum Communication Cost per Symbol

Control with Minimum Communication Cost per Symbol Contro with Minimum Communication Cost per ymbo Justin Pearson, João P. Hespanha, Danie Liberzon Abstract We address the probem of stabiizing a continuous-time inear time-invariant process under communication

More information

Disturbance decoupling by measurement feedback

Disturbance decoupling by measurement feedback Preprints of the 19th Word Congress The Internationa Federation of Automatic Contro Disturbance decouping by measurement feedback Arvo Kadmäe, Üe Kotta Institute of Cybernetics at TUT, Akadeemia tee 21,

More information

Tracking Control of Multiple Mobile Robots

Tracking Control of Multiple Mobile Robots Proceedings of the 2001 IEEE Internationa Conference on Robotics & Automation Seou, Korea May 21-26, 2001 Tracking Contro of Mutipe Mobie Robots A Case Study of Inter-Robot Coision-Free Probem Jurachart

More information

Limited magnitude error detecting codes over Z q

Limited magnitude error detecting codes over Z q Limited magnitude error detecting codes over Z q Noha Earief choo of Eectrica Engineering and Computer cience Oregon tate University Corvais, OR 97331, UA Emai: earief@eecsorstedu Bea Bose choo of Eectrica

More information

Problem set 6 The Perron Frobenius theorem.

Problem set 6 The Perron Frobenius theorem. Probem set 6 The Perron Frobenius theorem. Math 22a4 Oct 2 204, Due Oct.28 In a future probem set I want to discuss some criteria which aow us to concude that that the ground state of a sef-adjoint operator

More information

Coupling of LWR and phase transition models at boundary

Coupling of LWR and phase transition models at boundary Couping of LW and phase transition modes at boundary Mauro Garaveo Dipartimento di Matematica e Appicazioni, Università di Miano Bicocca, via. Cozzi 53, 20125 Miano Itay. Benedetto Piccoi Department of

More information

(f) is called a nearly holomorphic modular form of weight k + 2r as in [5].

(f) is called a nearly holomorphic modular form of weight k + 2r as in [5]. PRODUCTS OF NEARLY HOLOMORPHIC EIGENFORMS JEFFREY BEYERL, KEVIN JAMES, CATHERINE TRENTACOSTE, AND HUI XUE Abstract. We prove that the product of two neary hoomorphic Hece eigenforms is again a Hece eigenform

More information

ACTIVE electrical networks, which require energizing sources for their operation, are most widely

ACTIVE electrical networks, which require energizing sources for their operation, are most widely 1 Stabiization heory For Active Muti-port Networks Mayuresh Bakshi Member, IEEE,, Virendra R Sue and Maryam Shoejai Baghini, Senior Member, IEEE, arxiv:1606.03194v1 [cs.sy] 10 Jun 2016 Abstract his paper

More information

Global Optimality Principles for Polynomial Optimization Problems over Box or Bivalent Constraints by Separable Polynomial Approximations

Global Optimality Principles for Polynomial Optimization Problems over Box or Bivalent Constraints by Separable Polynomial Approximations Goba Optimaity Principes for Poynomia Optimization Probems over Box or Bivaent Constraints by Separabe Poynomia Approximations V. Jeyakumar, G. Li and S. Srisatkunarajah Revised Version II: December 23,

More information

Uniprocessor Feasibility of Sporadic Tasks with Constrained Deadlines is Strongly conp-complete

Uniprocessor Feasibility of Sporadic Tasks with Constrained Deadlines is Strongly conp-complete Uniprocessor Feasibiity of Sporadic Tasks with Constrained Deadines is Strongy conp-compete Pontus Ekberg and Wang Yi Uppsaa University, Sweden Emai: {pontus.ekberg yi}@it.uu.se Abstract Deciding the feasibiity

More information

The Symmetric and Antipersymmetric Solutions of the Matrix Equation A 1 X 1 B 1 + A 2 X 2 B A l X l B l = C and Its Optimal Approximation

The Symmetric and Antipersymmetric Solutions of the Matrix Equation A 1 X 1 B 1 + A 2 X 2 B A l X l B l = C and Its Optimal Approximation The Symmetric Antipersymmetric Soutions of the Matrix Equation A 1 X 1 B 1 + A 2 X 2 B 2 + + A X B C Its Optima Approximation Ying Zhang Member IAENG Abstract A matrix A (a ij) R n n is said to be symmetric

More information

Restricted weak type on maximal linear and multilinear integral maps.

Restricted weak type on maximal linear and multilinear integral maps. Restricted weak type on maxima inear and mutiinear integra maps. Oscar Basco Abstract It is shown that mutiinear operators of the form T (f 1,..., f k )(x) = R K(x, y n 1,..., y k )f 1 (y 1 )...f k (y

More information

Maximizing Sum Rate and Minimizing MSE on Multiuser Downlink: Optimality, Fast Algorithms and Equivalence via Max-min SIR

Maximizing Sum Rate and Minimizing MSE on Multiuser Downlink: Optimality, Fast Algorithms and Equivalence via Max-min SIR 1 Maximizing Sum Rate and Minimizing MSE on Mutiuser Downink: Optimaity, Fast Agorithms and Equivaence via Max-min SIR Chee Wei Tan 1,2, Mung Chiang 2 and R. Srikant 3 1 Caifornia Institute of Technoogy,

More information

Generalized Bell polynomials and the combinatorics of Poisson central moments

Generalized Bell polynomials and the combinatorics of Poisson central moments Generaized Be poynomias and the combinatorics of Poisson centra moments Nicoas Privaut Division of Mathematica Sciences Schoo of Physica and Mathematica Sciences Nanyang Technoogica University SPMS-MAS-05-43,

More information

MARKOV CHAINS AND MARKOV DECISION THEORY. Contents

MARKOV CHAINS AND MARKOV DECISION THEORY. Contents MARKOV CHAINS AND MARKOV DECISION THEORY ARINDRIMA DATTA Abstract. In this paper, we begin with a forma introduction to probabiity and expain the concept of random variabes and stochastic processes. After

More information

Steepest Descent Adaptation of Min-Max Fuzzy If-Then Rules 1

Steepest Descent Adaptation of Min-Max Fuzzy If-Then Rules 1 Steepest Descent Adaptation of Min-Max Fuzzy If-Then Rues 1 R.J. Marks II, S. Oh, P. Arabshahi Λ, T.P. Caude, J.J. Choi, B.G. Song Λ Λ Dept. of Eectrica Engineering Boeing Computer Services University

More information

UNIFORM CONVERGENCE OF MULTIPLIER CONVERGENT SERIES

UNIFORM CONVERGENCE OF MULTIPLIER CONVERGENT SERIES royecciones Vo. 26, N o 1, pp. 27-35, May 2007. Universidad Catóica de Norte Antofagasta - Chie UNIFORM CONVERGENCE OF MULTILIER CONVERGENT SERIES CHARLES SWARTZ NEW MEXICO STATE UNIVERSITY Received :

More information

A Solution to the 4-bit Parity Problem with a Single Quaternary Neuron

A Solution to the 4-bit Parity Problem with a Single Quaternary Neuron Neura Information Processing - Letters and Reviews Vo. 5, No. 2, November 2004 LETTER A Soution to the 4-bit Parity Probem with a Singe Quaternary Neuron Tohru Nitta Nationa Institute of Advanced Industria

More information

Distributed average consensus: Beyond the realm of linearity

Distributed average consensus: Beyond the realm of linearity Distributed average consensus: Beyond the ream of inearity Usman A. Khan, Soummya Kar, and José M. F. Moura Department of Eectrica and Computer Engineering Carnegie Meon University 5 Forbes Ave, Pittsburgh,

More information

(This is a sample cover image for this issue. The actual cover is not yet available at this time.)

(This is a sample cover image for this issue. The actual cover is not yet available at this time.) (This is a sampe cover image for this issue The actua cover is not yet avaiabe at this time) This artice appeared in a journa pubished by Esevier The attached copy is furnished to the author for interna

More information

Maximum likelihood decoding of trellis codes in fading channels with no receiver CSI is a polynomial-complexity problem

Maximum likelihood decoding of trellis codes in fading channels with no receiver CSI is a polynomial-complexity problem 1 Maximum ikeihood decoding of treis codes in fading channes with no receiver CSI is a poynomia-compexity probem Chun-Hao Hsu and Achieas Anastasopouos Eectrica Engineering and Computer Science Department

More information

Statistical Inference, Econometric Analysis and Matrix Algebra

Statistical Inference, Econometric Analysis and Matrix Algebra Statistica Inference, Econometric Anaysis and Matrix Agebra Bernhard Schipp Water Krämer Editors Statistica Inference, Econometric Anaysis and Matrix Agebra Festschrift in Honour of Götz Trenker Physica-Verag

More information

ORTHOGONAL MULTI-WAVELETS FROM MATRIX FACTORIZATION

ORTHOGONAL MULTI-WAVELETS FROM MATRIX FACTORIZATION J. Korean Math. Soc. 46 2009, No. 2, pp. 281 294 ORHOGONAL MLI-WAVELES FROM MARIX FACORIZAION Hongying Xiao Abstract. Accuracy of the scaing function is very crucia in waveet theory, or correspondingy,

More information

Investigation on spectrum of the adjacency matrix and Laplacian matrix of graph G l

Investigation on spectrum of the adjacency matrix and Laplacian matrix of graph G l Investigation on spectrum of the adjacency matrix and Lapacian matrix of graph G SHUHUA YIN Computer Science and Information Technoogy Coege Zhejiang Wani University Ningbo 3500 PEOPLE S REPUBLIC OF CHINA

More information

Guaranteed Cost H Controller Synthesis for Switched Systems Defined on Semialgebraic

Guaranteed Cost H Controller Synthesis for Switched Systems Defined on Semialgebraic Aaborg Universitet Guaranteed Cost H Controer Synthesis for Switched Systems Defined on Semiagebraic Sets Ahmadi, Mohamadreza; Mojaai, Hamed; Wisniewski, Rafa Pubished in: Noninear Anaysis: Hybrid Systems

More information

SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS

SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS ISEE 1 SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS By Yingying Fan and Jinchi Lv University of Southern Caifornia This Suppementary Materia

More information

Source and Relay Matrices Optimization for Multiuser Multi-Hop MIMO Relay Systems

Source and Relay Matrices Optimization for Multiuser Multi-Hop MIMO Relay Systems Source and Reay Matrices Optimization for Mutiuser Muti-Hop MIMO Reay Systems Yue Rong Department of Eectrica and Computer Engineering, Curtin University, Bentey, WA 6102, Austraia Abstract In this paper,

More information

An Algorithm for Pruning Redundant Modules in Min-Max Modular Network

An Algorithm for Pruning Redundant Modules in Min-Max Modular Network An Agorithm for Pruning Redundant Modues in Min-Max Moduar Network Hui-Cheng Lian and Bao-Liang Lu Department of Computer Science and Engineering, Shanghai Jiao Tong University 1954 Hua Shan Rd., Shanghai

More information

Integrating Factor Methods as Exponential Integrators

Integrating Factor Methods as Exponential Integrators Integrating Factor Methods as Exponentia Integrators Borisav V. Minchev Department of Mathematica Science, NTNU, 7491 Trondheim, Norway Borko.Minchev@ii.uib.no Abstract. Recenty a ot of effort has been

More information

Stochastic Automata Networks (SAN) - Modelling. and Evaluation. Paulo Fernandes 1. Brigitte Plateau 2. May 29, 1997

Stochastic Automata Networks (SAN) - Modelling. and Evaluation. Paulo Fernandes 1. Brigitte Plateau 2. May 29, 1997 Stochastic utomata etworks (S) - Modeing and Evauation Pauo Fernandes rigitte Pateau 2 May 29, 997 Institut ationa Poytechnique de Grenobe { IPG Ecoe ationae Superieure d'informatique et de Mathematiques

More information

Supplementary Appendix (not for publication) for: The Value of Network Information

Supplementary Appendix (not for publication) for: The Value of Network Information Suppementary Appendix not for pubication for: The Vaue of Network Information Itay P. Fainmesser and Andrea Gaeotti September 6, 03 This appendix incudes the proof of Proposition from the paper "The Vaue

More information

6 Wave Equation on an Interval: Separation of Variables

6 Wave Equation on an Interval: Separation of Variables 6 Wave Equation on an Interva: Separation of Variabes 6.1 Dirichet Boundary Conditions Ref: Strauss, Chapter 4 We now use the separation of variabes technique to study the wave equation on a finite interva.

More information

Lecture Note 3: Stationary Iterative Methods

Lecture Note 3: Stationary Iterative Methods MATH 5330: Computationa Methods of Linear Agebra Lecture Note 3: Stationary Iterative Methods Xianyi Zeng Department of Mathematica Sciences, UTEP Stationary Iterative Methods The Gaussian eimination (or

More information

On the Performance of Wireless Energy Harvesting Networks in a Boolean-Poisson Model

On the Performance of Wireless Energy Harvesting Networks in a Boolean-Poisson Model On the Performance of Wireess Energy Harvesting Networks in a Booean-Poisson Mode Han-Bae Kong, Ian Fint, Dusit Niyato, and Nicoas Privaut Schoo of Computer Engineering, Nanyang Technoogica University,

More information

An approximate method for solving the inverse scattering problem with fixed-energy data

An approximate method for solving the inverse scattering problem with fixed-energy data J. Inv. I-Posed Probems, Vo. 7, No. 6, pp. 561 571 (1999) c VSP 1999 An approximate method for soving the inverse scattering probem with fixed-energy data A. G. Ramm and W. Scheid Received May 12, 1999

More information

1D Heat Propagation Problems

1D Heat Propagation Problems Chapter 1 1D Heat Propagation Probems If the ambient space of the heat conduction has ony one dimension, the Fourier equation reduces to the foowing for an homogeneous body cρ T t = T λ 2 + Q, 1.1) x2

More information

Appendix of the Paper The Role of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Model

Appendix of the Paper The Role of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Model Appendix of the Paper The Roe of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Mode Caio Ameida cameida@fgv.br José Vicente jose.vaentim@bcb.gov.br June 008 1 Introduction In this

More information

DISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE

DISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE DISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE Yury Iyushin and Anton Mokeev Saint-Petersburg Mining University, Vasiievsky Isand, 1 st ine, Saint-Petersburg,

More information

C. Fourier Sine Series Overview

C. Fourier Sine Series Overview 12 PHILIP D. LOEWEN C. Fourier Sine Series Overview Let some constant > be given. The symboic form of the FSS Eigenvaue probem combines an ordinary differentia equation (ODE) on the interva (, ) with a

More information

12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes

12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes Maima and Minima 1. Introduction In this Section we anayse curves in the oca neighbourhood of a stationary point and, from this anaysis, deduce necessary conditions satisfied by oca maima and oca minima.

More information

PRIME TWISTS OF ELLIPTIC CURVES

PRIME TWISTS OF ELLIPTIC CURVES PRIME TWISTS OF ELLIPTIC CURVES DANIEL KRIZ AND CHAO LI Abstract. For certain eiptic curves E/Q with E(Q)[2] = Z/2Z, we prove a criterion for prime twists of E to have anaytic rank 0 or 1, based on a mod

More information

VI.G Exact free energy of the Square Lattice Ising model

VI.G Exact free energy of the Square Lattice Ising model VI.G Exact free energy of the Square Lattice Ising mode As indicated in eq.(vi.35), the Ising partition function is reated to a sum S, over coections of paths on the attice. The aowed graphs for a square

More information

u(x) s.t. px w x 0 Denote the solution to this problem by ˆx(p, x). In order to obtain ˆx we may simply solve the standard problem max x 0

u(x) s.t. px w x 0 Denote the solution to this problem by ˆx(p, x). In order to obtain ˆx we may simply solve the standard problem max x 0 Bocconi University PhD in Economics - Microeconomics I Prof M Messner Probem Set 4 - Soution Probem : If an individua has an endowment instead of a monetary income his weath depends on price eves In particuar,

More information

Stability Analysis and Application of Kalman Filtering with Irregularly Sampled Measurements

Stability Analysis and Application of Kalman Filtering with Irregularly Sampled Measurements 212 American Contro Conference Fairmont Queen Eizabeth, Montréa, Canada June 27-June 29, 212 Stabiity Anaysis and Appication of Kaman Fitering with Irreguary Samped Measurements Huazhen Fang and Raymond

More information

An explicit Jordan Decomposition of Companion matrices

An explicit Jordan Decomposition of Companion matrices An expicit Jordan Decomposition of Companion matrices Fermín S V Bazán Departamento de Matemática CFM UFSC 88040-900 Forianópois SC E-mai: fermin@mtmufscbr S Gratton CERFACS 42 Av Gaspard Coriois 31057

More information

Theory of Generalized k-difference Operator and Its Application in Number Theory

Theory of Generalized k-difference Operator and Its Application in Number Theory Internationa Journa of Mathematica Anaysis Vo. 9, 2015, no. 19, 955-964 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/10.12988/ijma.2015.5389 Theory of Generaized -Difference Operator and Its Appication

More information

Paper presented at the Workshop on Space Charge Physics in High Intensity Hadron Rings, sponsored by Brookhaven National Laboratory, May 4-7,1998

Paper presented at the Workshop on Space Charge Physics in High Intensity Hadron Rings, sponsored by Brookhaven National Laboratory, May 4-7,1998 Paper presented at the Workshop on Space Charge Physics in High ntensity Hadron Rings, sponsored by Brookhaven Nationa Laboratory, May 4-7,998 Noninear Sef Consistent High Resoution Beam Hao Agorithm in

More information

Physics 566: Quantum Optics Quantization of the Electromagnetic Field

Physics 566: Quantum Optics Quantization of the Electromagnetic Field Physics 566: Quantum Optics Quantization of the Eectromagnetic Fied Maxwe's Equations and Gauge invariance In ecture we earned how to quantize a one dimensiona scaar fied corresponding to vibrations on

More information

arxiv: v1 [cs.sy] 2 Apr 2019

arxiv: v1 [cs.sy] 2 Apr 2019 On the Existence of a Fixed Spectrum for a Multi-channel Linear System: A Matroid Theory Approach F Liu 1 and A S Morse 1 arxiv:190401499v1 [cssy] 2 Apr 2019 Abstract Conditions for the existence of a

More information

arxiv: v1 [math.co] 12 May 2013

arxiv: v1 [math.co] 12 May 2013 EMBEDDING CYCLES IN FINITE PLANES FELIX LAZEBNIK, KEITH E. MELLINGER, AND SCAR VEGA arxiv:1305.2646v1 [math.c] 12 May 2013 Abstract. We define and study embeddings of cyces in finite affine and projective

More information

General Decay of Solutions in a Viscoelastic Equation with Nonlinear Localized Damping

General Decay of Solutions in a Viscoelastic Equation with Nonlinear Localized Damping Journa of Mathematica Research with Appications Jan.,, Vo. 3, No., pp. 53 6 DOI:.377/j.issn:95-65...7 Http://jmre.dut.edu.cn Genera Decay of Soutions in a Viscoeastic Equation with Noninear Locaized Damping

More information

Math 124B January 17, 2012

Math 124B January 17, 2012 Math 124B January 17, 212 Viktor Grigoryan 3 Fu Fourier series We saw in previous ectures how the Dirichet and Neumann boundary conditions ead to respectivey sine and cosine Fourier series of the initia

More information

FOURIER SERIES ON ANY INTERVAL

FOURIER SERIES ON ANY INTERVAL FOURIER SERIES ON ANY INTERVAL Overview We have spent considerabe time earning how to compute Fourier series for functions that have a period of 2p on the interva (-p,p). We have aso seen how Fourier series

More information

Week 6 Lectures, Math 6451, Tanveer

Week 6 Lectures, Math 6451, Tanveer Fourier Series Week 6 Lectures, Math 645, Tanveer In the context of separation of variabe to find soutions of PDEs, we encountered or and in other cases f(x = f(x = a 0 + f(x = a 0 + b n sin nπx { a n

More information

Indirect Optimal Control of Dynamical Systems

Indirect Optimal Control of Dynamical Systems Computationa Mathematics and Mathematica Physics, Vo. 44, No. 3, 24, pp. 48 439. Transated from Zhurna Vychisite noi Matematiki i Matematicheskoi Fiziki, Vo. 44, No. 3, 24, pp. 444 466. Origina Russian

More information

Equilibrium of Heterogeneous Congestion Control Protocols

Equilibrium of Heterogeneous Congestion Control Protocols Equiibrium of Heterogeneous Congestion Contro Protocos Ao Tang Jiantao Wang Steven H. Low EAS Division, Caifornia Institute of Technoogy Mung Chiang EE Department, Princeton University Abstract When heterogeneous

More information

K a,k minors in graphs of bounded tree-width *

K a,k minors in graphs of bounded tree-width * K a,k minors in graphs of bounded tree-width * Thomas Böhme Institut für Mathematik Technische Universität Imenau Imenau, Germany E-mai: tboehme@theoinf.tu-imenau.de and John Maharry Department of Mathematics

More information

Numerical solution of one dimensional contaminant transport equation with variable coefficient (temporal) by using Haar wavelet

Numerical solution of one dimensional contaminant transport equation with variable coefficient (temporal) by using Haar wavelet Goba Journa of Pure and Appied Mathematics. ISSN 973-1768 Voume 1, Number (16), pp. 183-19 Research India Pubications http://www.ripubication.com Numerica soution of one dimensiona contaminant transport

More information

Position Control of Rolling Skateboard

Position Control of Rolling Skateboard osition Contro of Roing kateboard Baazs Varszegi enes Takacs Gabor tepan epartment of Appied Mechanics, Budapest University of Technoogy and Economics, Budapest, Hungary (e-mai: varszegi@mm.bme.hu) MTA-BME

More information

Smoothness equivalence properties of univariate subdivision schemes and their projection analogues

Smoothness equivalence properties of univariate subdivision schemes and their projection analogues Numerische Mathematik manuscript No. (wi be inserted by the editor) Smoothness equivaence properties of univariate subdivision schemes and their projection anaogues Phiipp Grohs TU Graz Institute of Geometry

More information

A unified framework for design and analysis of networked and quantized control systems

A unified framework for design and analysis of networked and quantized control systems 1 A unified framework for design and anaysis of networked and quantized contro systems Dragan ešić and Danie Liberzon Abstract We generaize and unify a range of recent resuts in quantized contro systems

More information

Discrete Techniques. Chapter Introduction

Discrete Techniques. Chapter Introduction Chapter 3 Discrete Techniques 3. Introduction In the previous two chapters we introduced Fourier transforms of continuous functions of the periodic and non-periodic (finite energy) type, as we as various

More information

$, (2.1) n="# #. (2.2)

$, (2.1) n=# #. (2.2) Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier

More information

King Fahd University of Petroleum & Minerals

King Fahd University of Petroleum & Minerals King Fahd University of Petroeum & Mineras DEPARTMENT OF MATHEMATICAL SCIENCES Technica Report Series TR 369 December 6 Genera decay of soutions of a viscoeastic equation Saim A. Messaoudi DHAHRAN 3161

More information

In this context, system identification is an important component

In this context, system identification is an important component 211 th IEEE Conference on Decision and Contro and European Contro Conference (CDC-ECC) Orando, FL, USA, December 12-1, 211 An improved strategy for NeuroMuscuar Bockade contro with parameter uncertainty

More information

Multiscale Domain Decomposition Preconditioners for 2 Anisotropic High-Contrast Problems UNCORRECTED PROOF

Multiscale Domain Decomposition Preconditioners for 2 Anisotropic High-Contrast Problems UNCORRECTED PROOF AQ Mutiscae Domain Decomposition Preconditioners for 2 Anisotropic High-Contrast Probems Yachin Efendiev, Juan Gavis, Raytcho Lazarov, Svetozar Margenov 2, 4 and Jun Ren 5 Department of Mathematics, TAMU,

More information

#A48 INTEGERS 12 (2012) ON A COMBINATORIAL CONJECTURE OF TU AND DENG

#A48 INTEGERS 12 (2012) ON A COMBINATORIAL CONJECTURE OF TU AND DENG #A48 INTEGERS 12 (2012) ON A COMBINATORIAL CONJECTURE OF TU AND DENG Guixin Deng Schoo of Mathematica Sciences, Guangxi Teachers Education University, Nanning, P.R.China dengguixin@ive.com Pingzhi Yuan

More information

Distributed Optimization With Local Domains: Applications in MPC and Network Flows

Distributed Optimization With Local Domains: Applications in MPC and Network Flows Distributed Optimization With Loca Domains: Appications in MPC and Network Fows João F. C. Mota, João M. F. Xavier, Pedro M. Q. Aguiar, and Markus Püsche Abstract In this paper we consider a network with

More information

A simple reliability block diagram method for safety integrity verification

A simple reliability block diagram method for safety integrity verification Reiabiity Engineering and System Safety 92 (2007) 1267 1273 www.esevier.com/ocate/ress A simpe reiabiity bock diagram method for safety integrity verification Haitao Guo, Xianhui Yang epartment of Automation,

More information

Adjustment of automatic control systems of production facilities at coal processing plants using multivariant physico- mathematical models

Adjustment of automatic control systems of production facilities at coal processing plants using multivariant physico- mathematical models IO Conference Series: Earth and Environmenta Science AER OEN ACCESS Adjustment of automatic contro systems of production faciities at coa processing pants using mutivariant physico- mathematica modes To

More information

Multiple Loop Self-Triggered Model Predictive Control for Network Scheduling and Control

Multiple Loop Self-Triggered Model Predictive Control for Network Scheduling and Control Mutipe Loop Sef-Triggered Mode Predictive Contro for Networ Scheduing and Contro Eri Henrisson, Danie E. Quevedo, Edwin G.W. Peters, Henri Sandberg, Kar Henri Johansson arxiv:5.38v [cs.sy] Feb 5 Abstract

More information

Efficient Visual-Inertial Navigation using a Rolling-Shutter Camera with Inaccurate Timestamps

Efficient Visual-Inertial Navigation using a Rolling-Shutter Camera with Inaccurate Timestamps Efficient Visua-Inertia Navigation using a Roing-Shutter Camera with Inaccurate Timestamps Chao X. Guo, Dimitrios G. Kottas, Ryan C. DuToit Ahmed Ahmed, Ruipeng Li and Stergios I. Roumeiotis Mutipe Autonomous

More information

Introduction to LMTO method

Introduction to LMTO method 1 Introduction to MTO method 24 February 2011; V172 P.Ravindran, FME-course on Ab initio Modeing of soar ce Materias 24 February 2011 Introduction to MTO method Ab initio Eectronic Structure Cacuations

More information

Homework #04 Answers and Hints (MATH4052 Partial Differential Equations)

Homework #04 Answers and Hints (MATH4052 Partial Differential Equations) Homework #4 Answers and Hints (MATH452 Partia Differentia Equations) Probem 1 (Page 89, Q2) Consider a meta rod ( < x < ), insuated aong its sides but not at its ends, which is initiay at temperature =

More information