Control with Minimum Communication Cost per Symbol

Size: px
Start display at page:

Download "Control with Minimum Communication Cost per Symbol"

Transcription

1 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 constraints. We assume that the sensor that measures the state is connected to the actuator through a finite capacity communication channe over which an encoder at the sensor sends symbos from a finite aphabet to a decoder at the actuator. We consider a situation where one symbo from the aphabet consumes no communication resources, whereas each of the others consumes one unit of communication resources to transmit. This paper expores how the imposition of imits on an encoder s bit-rate and average resource consumption affect the encoder/decoder/controer s abiity to keep the process bounded. The main resut is a necessary and sufficient condition for a bounding encoder/decoder/controer which depends on the encoder s average bit rate, its average resource consumption, and the unstabe eigenvaues of the process. I. IRODUCTIO This paper addresses the probem of stabiizing a continuous-time inear time-invariant process under communication constraints. As in [1 7], we assume that the sensor that measures the state is connected to the actuator through a finite capacity communication channe. At each samping time, an encoder sends a symbo through the channe. The probem of determining whether or not it is possibe to bound the state of the process under this type of encoding scheme is not new; it was estabished in [2 4] that a necessary and sufficient condition for stabiity can be expressed as a simpe reationship between the unstabe eigenvaues of A and the average communication bit-rate. We expand upon this resut by considering the notion that encoders can effectivey save communication resources by not transmitting information, whie noting that the absence of an expicit transmission nevertheess conveys information. To capture this, we suppose that one symbo from the aphabet consumes no communication resources to transmit, whereas each of the others consumes one unit of communication resources. We then proceed to define the average cost per symbo of an encoder, which is essentiay the average fraction of non-free symbos emitted. This paper s main technica contribution is a necessary and sufficient condition for the existence of an encoder/decoder/controer that bounds the state of the process. This condition depends on the channe s average bit rate, the encoder s average cost per symbo, and the unstabe eigenvaues of A. Our resut extends [2] in the sense that as the constraint on the average cost per symbo is aowed to increase (becomes ooser), our necessary and sufficient condition becomes the condition from [2]. As with [2 4], our resut is constructive in the sense that we describe a famiy of encoder/decoder pairs that bound the process when our condition hods. A counterintuitive coroary to our main resut shows that if the process may be bounded with average bit-rate r, then there exists a bounding encoder/decoder/controer with average bit-rate r which uses no more than 50% non-free symbos in its symbo-stream. The remainder of this paper is organized as foows. ection III contains the main negative resut of the paper, namey that boundedness is not possibe when our condition does not hod. To prove this resut we actuay show that it is not possibe to bound the process with a arge cass of encoders which we ca M-of- encoders that incudes a the encoders with average cost per symbo not exceeding a given threshod. ection IV contains the positive resut of the paper, showing that when our condition does hod, there is an encoder/decoder pair that can bound the process; we provide the encoding scheme. II. PROBLEM TATEME Consider a stabiizabe inear time-invariant process 9x Ax Bu, x P R n, u P R m, (1) for which it is known that xp0q beongs to some bounded set X 0 R n. A sensor that measures the state xptq is connected to the actuator through a finite-data-rate, error-free, and deay-free communication channe. An encoder coocated with the sensor sampes the state once every T time units, and from this sequence of measurements txpkt q : k P 0 u causay constructs a sequence of symbos ts k P A : k P 0 u from a finite aphabet A. The encoder sends this symbo sequence through the channe at a rate of 1 symbo every T time units to a decoder/controer coocated with the actuator, which causay constructs the contro signa uptq, t 0 from the sequence of symbos ts k P A : k P 0 u that arrive at the decoder. We assume the channe faithfuy transmits each symbo without error. The positive time T between successive sampings is caed the samping period and has units of time units. The average bit-rate of an encoder/decoder pair is the amount of information that the encoder transmits in units of bits per time unit. For an encoder/decoder pair whose encoder transmits one symbo from an aphabet A every T time units, the pair s average bit-rate is given by r og 2 A. (2) T We consider encoder/decoder pairs whose aphabets each contain one specia symbo 0 P A that can be transmitted without consuming any communication resources, and

2 additiona symbos that each require one unit of communication resources per transmission. One can think of the free symbo 0 as the absence of an expicit transmission. The communication resources at stake may be energy, time, or any other resource that may be consumed in the course of the communication process. In order to capture the average rate that an encoder consumes communication resources, we define the average cost per symbo of an encoder as foows: We say an encoder has average cost per symbo not exceeding γ max if for every symbo sequence ts k u the encoder may generate, we have 1 M 1 Ņ km I sk 0 γ max O 1 M 1 for a integers, M satisfying M 0, where I sk 0 1 if the kth symbo is not the free symbo, and 0 if it is. The summation in (3) captures the tota resources spent transmitting symbos s M, s M 1,..., s. Motivating this definition of average cost per symbo is the observation that the efthand side has the intuitive interpretation as the average cost per transmitted symbo between symbos s M and s. As M Ñ 8, the rightmost big-oh term vanishes, eaving γ max as an upper bound on the average ong-term cost per symbo of the symbo sequence. ote that the average cost per symbo of any encoder never exceeds 1 and does not depend on the samping period T. Whereas an encoder/decoder pair s average bit-rate r ony depends on its symbo aphabet A and samping period T, its average cost per symbo depends on every possibe symbo sequence it may generate, and therefore depends on the encoder/decoder pair, the controer, the process (1), and the initia condition xp0q. The specific question considered in this paper is: under what conditions on the bit-rate and average cost per symbo does there exist a controer and encoder/decoder pair that keep the state of process (1) bounded? III. ECEARY CODITIO FOR BOUDEDE WITH LIMITED-COMMUICATIO ECODER It is we known from [2 4] that it is possibe to construct a controer and encoder/decoder pair that bounds process (1) with average bit-rate r ony if r n 2 λ i ras, (4) where n denotes the base-e ogarithm, and the summation is over a eigenvaues of A with nonnegative rea part. The resut that foows shows that a arger bit-rate may be needed when one poses constraints on the encoder s average cost per symbo γ max. pecificay, when γ max {p 1q the (necessary) stabiity condition reduces to (4), but when γ max {p 1q a bit-rate arger than (4) is necessary to compensate for the diution of information content in the transmitted symbos due to the constraint on the average cost per symbo. Theorem 1: Consider an encoder/decoder pair with samping period T, an aphabet with nonfree symbos and one (3) free symbo, and an average cost per symbo not exceeding γ max. If this pair keeps the state of process (1) bounded for every initia condition x 0 P X 0, then we must have r fpγ max, q n 2 λ i ras, (5) where the bit-rate r is reated to and T via Equation (2), and the function f : r0, 1s r0, 8q Ñ r0, 8q is defined as # Hpγq γ og2 og fpγ, q 2 p 1q 0 γ γ 1, (6) and Hppq p og 2 ppq p1 pq og 2 p1 pq is the base-2 entropy of a Bernoui random variabe with parameter p. Remark 1: The function f is potted in Figure 1 for severa vaues of. It is worth making two observations regarding f. The first observation is that the function fpγ, q is monotone nonincreasing in for any fixed γ P r0, 1s, which impies that smaer aphabets are preferabe to arge ones when trying to satisfy (5) with a given fixed bit-rate. The second observation is that the average cost per time unit, which is γ{t, can be made arbitrariy sma in two different ways. 1) One coud pick T very arge, then everaging the fact that rf pγ, q is monotone increasing in, pick arge enough to satisfy (5). This approach has two downsides: First, with a arge samping period, the state, athough remaining bounded, can grow quite arge between transmissions. econd, arge means that the encoder/decoder pair must store and process a arge symbo ibrary, adding compexity to the pair s impementation. 2) Aternativey, one can make γ{t arbitrary sma by observing that the sequences γ k e k, T k e k? k, k P 0 have the property that γ k Ñ 0, T k Ñ 0, and γ k {T k Ñ 0, but Hpγ k q{t k Ñ 8, and hence r k fpγ k, q n 2 Ñ 8 (where r k og 2 p 1q{T k ). This means that one can find k P 0 to make the average cost per time unit γ k {T k arbitrariy sma, and aso satisfy the necessary condition (5). The drawback of this approach is that to achieve a sma samping period T in practice requires an encoder/decoder pair with a very precise cock. Remark 2: The addition of the free symbo effectivey increases the data rate without increasing the rate of resource consumption, as seen by the foowing two observations: Without the free symbos, the size of the aphabet woud be and the bit-rate woud be og 2 pq{t. It coud happen that this bit-rate is too sma to bound the pant, yet after the introduction of the free symbo, the condition (5) is satisfied. ince γ max is the fraction of non-free symbos, then the quantity r γ max is the number of bits per time unit spent transmitting non-free symbos. But since fpγ, q γ, again we see that the free symbos hep satisfy (5).

3 Fig. 1. A pot of f pγ, q for 1, 4, 20. A. Theorem 1 Proof etup We ead up to the proof of Theorem 1 by first estabishing three emmas centered around a restricted arge cass of encoders caed M-of- encoders. We first define M-of- encoders, which essentiay partition their symbo sequences into -ength codewords, each with M or fewer non-free symbos. Lemma 1 demonstrates that every encoder with a imited average cost per symbo is an M-of- encoder for appropriate and M. ext, in Lemma 2 we estabish a reationship between the number of codewords avaiabe to an M-of- encoder and the function f as defined in (6). Then, in Lemma 3 we everage previous work to estabish a necessary condition for an M-of- encoder to bound the state of the process. Finay, the proof of Theorem 1 everages these three resuts. To introduce the cass of M-of- encoders, we define an -symbo codeword to be a sequence ts 1, s 2,..., s u of consecutive symbos starting at an index k 1, with P 0. An M-of- encoder is an encoder for which every -symbo codeword has M or fewer non-free symbos, i.e., k 1 I sk 0 P 0. (7) The tota number of distinct -symbo codewords avaiabe to an M-of- encoder is thus given by Lp, M, q tm u i0 i i, (8) where the ith term in the summation counts the number of -symbo codewords with exacty i non-free symbos. ote that in keeping with the probem setup, the M-of- encoders considered here each draw their symbos from the symbo ibrary A t0, 1,..., u and transmit symbos with samping period T. An intuitive property of M-of- encoders is that an M-of- encoder has average cost per symbo not exceeding M{. To see this, suppose M and are fixed and et 1, 2 P 0 be arbitrary with The width of the interva r 1, 2 s may not be an integer mutipe of, but nevertheess it is between two integer mutipes of : we have pk 2q 2 1 ˆ 2 ˆ 1 k for some non-negative integer k P 0. From the first inequaity we obtain k , and from the second inequaity we have 2 k 1 I sk 0 km, because in each of the k -ength intervas containing r 1, 2 s there are at most M non-free symbos. Combining these yieds 2 k 1 I sk 0 M p 2 1 1q 0, (9) where 0 is some constant. Dividing both sides by p 2 1 1q and rearranging, (3) foows. The fact that an M-of- encoder refrains from sending expensive codewords effectivey reduces its abiity to transmit information. Indeed, since an M-of- encoder takes time units to transmit one of Lp, M, q codewords, the encoder transmits merey og 2 Lp,M,q bits per time unit. For M, this is stricty ess than the encoder s average bit-rate r. The first emma, proved in the appendix, shows that the set of M-of- encoders is compete in the sense that every encoder with average cost per symbo not exceeding a finite threshod γ max is actuay an M-of- encoder for sufficienty arge and M γ max. Lemma 1: For every encoder/decoder pair with average bit-rate r and average cost per symbo not exceeding γ max P r0, 1s, and every constant ɛ 0, there exist positive integers M and with M γ max p1 ɛq such that the encoder is an M-of- encoder. The next emma estabishes a reationship between the number of codewords Lp, M, q avaiabe to an M-of- encoder and the function f defined in (6). Lemma 2: For any, P 0 and γ P r0, 1s, the function L defined in (8) and the function f defined in (6) satisfy og 2 Lp, γ, q og 2 p 1qfpγ, q, (10) with equaity hoding ony when γ 0 or γ 1. Moreover, we have asymptotic equaity in the sense that og im 2 Lp, γ, q og Ñ8 2 p 1qfpγ, q. (11) Proof of Lemma 2. Let and be arbitrary positive integers. First we prove (10) for γ P p0, s. Appying the 1

4 Binomia Theorem to the identity 1 pγ p1 γqq, we obtain Ņ 1 γ i p1 γq i i i0 ince each term in the summation is positive, keeping ony the first tγu terms yieds the inequaity 1 tγu i0 i γ i p1 γq i (12) ext, a cacuation presented as Lemma 5 in the appendix reveas that γ i p1 γq i Hpγq i 2 γ (13) for a, P 0, γ P p0, {p 1qs, and i P r0, γs. Using this in (12) and simpifying yieds og 2 Lp, γ, q Hpγq γ og 2 (14) for a, P 0 and γ P p0, {p 1qs. ince og 2 p 1qfpγ, q Hpγq γ og 2 when γ P, (14) 0, 1 proves (10) for γ P p0, 1s. ext, suppose γ P p 1, 1q and observe from (8) that Lp, M, q is a sum of positive terms whose index reaches tmu, hence Lp, γ, q is stricty ess than Lp,, q for any γ 1. We concude that og 2 Lp, γ, q og 2 Lp,, q (15) og 2 p 1q (16) ince og 2 p 1qfpγ, q og 2 p 1q for γ P p 1, 1q, this concudes the proof of (10) for γ P p0, 1q. The proof of (10) for γ 0, foows merey from inspection of (10), and the γ 1 case foows from the equaity in (15). The proof of the asymptotic resut (11) appears in [8]. This concudes the proof of Lemma 2. The foowing emma provides a necessary condition for an M-of- encoder to bound process (1). Lemma 3: Consider an M-of- encoder/decoder pair using symbos t0,..., u with samping period T. If the pair keeps the state of (1) bounded for every initia condition, then we must have n Lp, M, q λ i ras. (17) Proof of Lemma 3. Consider an encoder/decoder/controer tripe that bounds the state of the process (1) and whose encoder is an M-of- encoder using symbos t0,..., u and samping period T. ince the encoder sends one of LpM,, q codewords every time units, the average bitrate of the encoder is og 2 LpM,, q{. Theorem 1 of [2] proves that if the state of the process (1) is bounded by an encoder/decoder/controer tripe under the communication constraints described in our probem setup, then the pair s average bit-rate must not be ess than 1 n 2 i:rλ λ iras 0 iras. Hence, we have og 2 LpM,, q 1 λ i ras, n 2 from which (17) foows. This proves the emma. ow we are ready to prove Theorem 1. Proof of Theorem 1. By Lemma 1, for any ɛ 0 there exist M, P 0 with M γ max p1 ɛq for which the encoder/decoder is an M-of- encoder. ince the state of the process is kept bounded, by Lemma 3 we have og λ i ras 2 Lp, M, q n 2. (18) ince L is monotonicay nondecreasing in its second argument and M γ max p1 ɛq, we have og 2 Lp, M, q og 2 Lp, γ max p1 ɛq, q. (19) Lemma 2 impies that og 2 Lp, γ max p1 ɛq, q rfpγ max p1 ɛq, q. (20) Combining these and etting ɛ Ñ 0, we obtain (5). This competes the proof of Theorem 1. IV. UFFICIE CODITIO FOR BOUDEDE WITH LIMITED-COMMUICATIO ECODER The previous section estabished a necessary condition (5) on the average bit-rate and average cost per symbo of an encoder/decoder pair in order to bound the state of (1). In this section, we show that that condition is aso sufficient for a bounding encoder/decoder to exist. The proof is constructive in that we provide the encoder/decoder. Theorem 2: Assume that A is diagonaizabe. For every P 0, T 0, and γ max P r0, 1s satisfying rfpγ max, q n 2 λ i ras, (21) where r is defined in (2) and the function f is defined in (6), there exists a controer and an encoder/decoder pair using nonfree symbos, samping period T, and average cost per symbo not exceeding γ max, that keeps the state of the process (1) bounded for every initia condition x 0 P X 0. The proof of Theorem 2 reies on the foowing emma, which provides a sufficient condition for the existence of an M-of- encoder to bound the state of process (1). Lemma 4: Assume that A is diagonaizabe. For every T 0 and, M, P 0 with M 0 satisfying n Lp, M, q λ i ras (22) there exists an M-of- encoder on aphabet t0,..., u with samping period T that keeps the state of the process (1) bounded for every initia condition.

5 Proof of Lemma 4. This proof buids on Theorem 2 from [2], which provides sufficient conditions on an encoder s average bit-rate to ensure the existence of a stabiizing controer and encoder/decoder pair. The resut states that if an encoder s average bit-rate r satisfies r 1 n 2 λ i ras, (23) where A is the continuous-time process matrix, then there exists a controer and encoder/decoder pair that bound the state of the process. An outine of the proof of this resut from [2] foows. The encoder works by pacing a bounding rectange around the voume where the state is known to ie, and partitions it into two pieces aong its argest axis. Then the encoder transmits a 0 or 1, depending on which subrectange the state ies in. ince the decoder can cacuate the bounding rectanges with its knowedge of X 0 and process dynamics in (1), it estimates the state to be the centroid of whichever sub-rectange the received symbo corresponds to. The decoder transmits this state estimate to the controer, which can be any stabiizing inear state-feedback controer. We now provide a summary of the proof of Lemma 4, and refer the reader to [8] for the detaied proof. uppose T,, M, satisfy the constraints in the statement of the emma as we as (22). Theorem 2 from [2] guarantees the existence of a stabiizing encoder/decoder/controer tripe with average bit-rate og 2 Lp, M, q{. A that remains is to adapt this encoder/decoder pair into our framework, i.e., an M-of- encoder that uses 1 symbos with samping period T. We do this by buiding an encoder which runs a copy of the encoder from [2] internay. The outer encoder feeds sampes of the state to the inner encoder, from which the inner encoder generates a string of 1 s and 0 s. Due to the bit-rate being og 2 Lp, M, q{, after T time units the inner encoder has generated a string of ength og 2 Lp, M, q. There are Lp, M, q possibe such strings, and so the outer encoder can map this string uniquey to one of the Lp, M, q possibe -ength codewords with M or fewer non-free symbos from the aphabet t0,..., u. The outer encoder transmits this codeword across the channe to the decoder, which performs the inverse mapping to recover the string of og 2 Lp, M, q 1 s and 0 s, which it deivers to its inner decoder. From this bit-string and knowedge of the process dynamics, the inner decoder computes a state estimate, which it deivers to the stabiizing inear state-feedback controer. This process repeats each T time units. ince the inner encoder/decboder pair bound the process, the outer M-of- encoder described here does as we. This proves the emma. ow we are ready to prove Theorem 2. Proof of Theorem 2. Assume that, T, and γ max satisfy (21), so that ɛ rfpγ max, q n 2 λ i ras 0. (24) n Lp,γmax,q Equation (11) estabishes that gets arbitrariy cose to rfpγ max, q as we increase, so we pick sufficienty arge to satisfy rfpγ max, q n 2 n Lp, γ max, q ɛ{2. (25) By (8), we have Lp, tγ max u, q Lp, γ max, q for every, γ max, and. etting M tγ max u, then by (24) and (25) we have found and M satisfying n Lp,M,q i:rλ λ iras 0 iras. Hence by Lemma 4, there exists an M-of- encoder/decoder which bounds the state of (1). ince a M-of- encoder/decoders have average cost per symbo not exceeding M{, and this encoder/decoder satisfies M{ γ max, we concude that this encoder has an average cost per symbo not exceeding γ max. This concudes the proof, since we have found the desired encoder. An unexpected consequence of Theorems 1 and 2 is that when it is possibe to keep the state of a process bounded with a given bit-rate r og 2 p 1q{T, one can aways find M-of- encoders that bound it for (essentiay) the same bit-rate and average cost per symbo not exceeding {p 1q, i.e., approximatey a fraction 1{p 1q of the symbos wi not consume communication resources. In the most advantageous case, the encoder/decoder use the aphabet t0, 1u and the encoder s symbo stream consumes no more than 50% of the communication resources. The price paid for using an encoder/decoder with average cost per symbo near {p 1q is that it may require prohibitivey ong codewords (arge ) as compared to an encoder with higher average cost per symbo. This is due to the fact that n Lp, γ max, q{ q is monotonicay nondecreasing in γ max, so that a smaer γ max generay requires a arger to satisfy (25). Coroary 1: If the process (1) can be bounded with an encoder/decoder pair with symbo aphabet t0, 1,..., u and samping period T, then for any ɛ 0 with T ɛ, there exists an M-of- encoder with aphabet t0, 1,..., u, samping period T ɛ, and average cost per symbo not exceeding {p 1q that bounds its state. Proof of Coroary 1. Let ɛ be arbitrary in p0, T q. ince the controer and encoder/decoder pair bound the system and the average cost per symbo never exceeds 1, by Theorem 1 we have rfp1, q n 2 λ i ras, (26) where r og 2 p 1q{T. ince fp{p 1q, q fp1, q for a P 0, we have og 2 p 1q f T ɛ 1, r f 1, rfp1, q n 2 λ i ras. (27) By Theorem 2, there exists an encoder/decoder pair with symbo aphabet t0,..., u, samping period T ɛ, and average cost per symbo not exceeding {p 1q which

6 bounds the state of (1). By Lemma 1, this encoder is an M-of- encoder for appropriatey chosen M and. V. COCLUIO AD FUTURE WORK In this paper we considered the probem of bounding the state of a continuous-time inear process under communication constraints. We considered constraints on both the average bit-rate and the average cost per symbo of an encoding scheme. Our main contribution was a necessary and sufficient condition on the process and constraints for which a bounding encoder/decoder/controer exists. In the absence of a imit on the average cost per symbo, the conditions recovered previous work. A surprising coroary to our main resut was the observation that one may impose a constraint on the average cost per symbo without necessariy needing to oosen the bit-rate constraint. pecificay, we proved that if a process may bounded with a particuar bit-rate, then there exists a (possiby very compex) encoder/decoder that can bound it using no more than 50% non-free symbos on average, yet obeying the same bit-rate. This was surprising because one woud expect the prohibition of some codewords woud require that the encoder aways compensate by transmitting at a higher bit-rate. We observed in Remark 1 that smaer aphabets incur a smaer penaty on the conditions for boundedness in Theorems 1 and 2. This suggests that encoding schemes with sma aphabets may be abe to bound the state of the process with bit-rates and average costs not far above the minimum theoretica bounds as estabished in Theorems 1 and 2. Eventbased contro strategies comprise one such cass of encoders; they use a sma number of non-free symbos to notify the decoder/controer about certain state-dependent events. We have preiminary resuts showing that event-based encoders can indeed be used to produce encoder/decoder pairs that amost achieve the minimum achievabe average cost per symbo that appears in Theorems 1 and 2. Finay, our probem setup considered merey wherether there exists a bounding encoder/decoder/controer tripe. It seems natura to extend this setup to finding stabiizing tripes. APPEDIX Proof of Lemma 1. By the definition of average cost per symbo not exceeding γ max in (3) there exists an integer 0 P 0 such that for any symbo sequence ts k u that the encoder generates, we have Ņ i1 I si0 0 γ P 0. (28) Pick P 0 arge enough to satisfy 0 2 ɛγ max and pick M t 0 2 γ max u. Combining these with (28), we obtain Ņ i1 I si0 0 γ max M 0 2 γ max γ max p1 ɛq, which estabishes that M γ max p1 ɛq. This competes the proof. Lemma 5: The foowing inequaity hods for a, P 0, q P p0, {p 1qs, and i P r0, qs: q i p1 qq i Hpqq i 2 q (29) where Hpqq q og 2 q p1 qq og 2 p1 qq is the the base- 2 entropy of a Bernoui random variabe with parameter q. Proof of Lemma 5. Let,, q, and i take arbitrary vaues from the sets described in the emma s statement. ince og 2 is a monotone increasing function, og 2 pq{p1 qqq for q 0 is maximized at the right endpoint vaue, q {p 1q, where it equas og 2. This eads to og 2 q og 2 p1 qq og 2 (30) for a P 0 and q P p0, {p 1qs. ext, since i P r0, qs by assumption, then i q 0. Mutipying (30) by i q, standard agebraic manipuations yied i og 2 q p iq og 2 p1 qq from which (29) foows. Hpqq pi qq og 2, REFERECE [1] R. Brockett and D. Liberzon, Quantized feedback stabiization of inear systems, Automatic Contro, IEEE Transactions on, vo. 45, no. 7, pp , Ju [2] J. P. Hespanha, A. Ortega, and L. Vasudevan, Towards the contro of inear systems with minimum bit-rate, in Proc. of the Int. ymp. on the Mathematica Theory of etworks and yst., Aug [3] G. air and R. Evans, Communication-imited stabiization of inear systems, in Decision and Contro, Proceedings of the 39th IEEE Conference on, vo. 1, 2000, pp vo.1. [4]. Tatikonda and. Mitter, Contro under communication constraints, Automatic Contro, IEEE Transactions on, vo. 49, no. 7, pp , juy [5] G.. air and R. J. Evans, Exponentia stabiisabiity of finite-dimensiona inear systems with imited data rates, Automatica, vo. 39, no. 4, pp , [6] A. ahai and. Mitter, The necessity and sufficiency of anytime capacity for stabiization of a inear system over a noisy communication ink part i: caar systems, Information Theory, IEEE Transactions on, vo. 52, no. 8, pp , Aug [7] A. Matveev and A. avkin, Mutirate stabiization of inear mutipe sensor systems via imited capacity communication channes, IAM Journa on Contro and Optimization, vo. 44, no. 2, pp , [8] J. Pearson, J. P. Hespanha, and D. Liberzon, Contro with minimum communication cost per symbo, University of Caifornia, anta Barbara, Tech. Rep., Mar [Onine]. Avaiabe: edu/ hespanha/techrep.htm

Control with Minimum Communication Cost per Symbol

Control with Minimum Communication Cost per Symbol Control with Minimum Communication Cost per ymbol Technical Report Justin Pearson, João P. Hespanha, Daniel Liberzon May 11, 2014 Abstract We address the problem of stabilizing a continuous-time linear

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

Iterative Decoding Performance Bounds for LDPC Codes on Noisy Channels

Iterative Decoding Performance Bounds for LDPC Codes on Noisy Channels Iterative Decoding Performance Bounds for LDPC Codes on Noisy Channes arxiv:cs/060700v1 [cs.it] 6 Ju 006 Chun-Hao Hsu and Achieas Anastasopouos Eectrica Engineering and Computer Science Department University

More information

FRST Multivariate Statistics. Multivariate Discriminant Analysis (MDA)

FRST Multivariate Statistics. Multivariate Discriminant Analysis (MDA) 1 FRST 531 -- Mutivariate Statistics Mutivariate Discriminant Anaysis (MDA) Purpose: 1. To predict which group (Y) an observation beongs to based on the characteristics of p predictor (X) variabes, using

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

Separation of Variables and a Spherical Shell with Surface Charge

Separation of Variables and a Spherical Shell with Surface Charge Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation

More information

Some Measures for Asymmetry of Distributions

Some Measures for Asymmetry of Distributions Some Measures for Asymmetry of Distributions Georgi N. Boshnakov First version: 31 January 2006 Research Report No. 5, 2006, Probabiity and Statistics Group Schoo of Mathematics, The University of Manchester

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

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

Pricing Multiple Products with the Multinomial Logit and Nested Logit Models: Concavity and Implications

Pricing Multiple Products with the Multinomial Logit and Nested Logit Models: Concavity and Implications Pricing Mutipe Products with the Mutinomia Logit and Nested Logit Modes: Concavity and Impications Hongmin Li Woonghee Tim Huh WP Carey Schoo of Business Arizona State University Tempe Arizona 85287 USA

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

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

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

Analysis of Emerson s Multiple Model Interpolation Estimation Algorithms: The MIMO Case

Analysis of Emerson s Multiple Model Interpolation Estimation Algorithms: The MIMO Case Technica Report PC-04-00 Anaysis of Emerson s Mutipe Mode Interpoation Estimation Agorithms: The MIMO Case João P. Hespanha Dae E. Seborg University of Caifornia, Santa Barbara February 0, 004 Anaysis

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

CS229 Lecture notes. Andrew Ng

CS229 Lecture notes. Andrew Ng CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view

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

4 Separation of Variables

4 Separation of Variables 4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE

More information

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is

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

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

Turbo Codes. Coding and Communication Laboratory. Dept. of Electrical Engineering, National Chung Hsing University

Turbo Codes. Coding and Communication Laboratory. Dept. of Electrical Engineering, National Chung Hsing University Turbo Codes Coding and Communication Laboratory Dept. of Eectrica Engineering, Nationa Chung Hsing University Turbo codes 1 Chapter 12: Turbo Codes 1. Introduction 2. Turbo code encoder 3. Design of intereaver

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

A proposed nonparametric mixture density estimation using B-spline functions

A proposed nonparametric mixture density estimation using B-spline functions A proposed nonparametric mixture density estimation using B-spine functions Atizez Hadrich a,b, Mourad Zribi a, Afif Masmoudi b a Laboratoire d Informatique Signa et Image de a Côte d Opae (LISIC-EA 4491),

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

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

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

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische

More information

Control with minimal cost-per-symbol encoding and quasi-optimality of event-based encoders Technical Report

Control with minimal cost-per-symbol encoding and quasi-optimality of event-based encoders Technical Report Control with minimal cost-per-symbol encoding and quasi-optimality of event-based encoders Technical Report Justin Pearson, João P. Hespanha, Daniel Liberzon : May 11, 2016 Abstract We consider the problem

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

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

STA 216 Project: Spline Approach to Discrete Survival Analysis

STA 216 Project: Spline Approach to Discrete Survival Analysis : Spine Approach to Discrete Surviva Anaysis November 4, 005 1 Introduction Athough continuous surviva anaysis differs much from the discrete surviva anaysis, there is certain ink between the two modeing

More information

On colorings of the Boolean lattice avoiding a rainbow copy of a poset arxiv: v1 [math.co] 21 Dec 2018

On colorings of the Boolean lattice avoiding a rainbow copy of a poset arxiv: v1 [math.co] 21 Dec 2018 On coorings of the Booean attice avoiding a rainbow copy of a poset arxiv:1812.09058v1 [math.co] 21 Dec 2018 Baázs Patkós Afréd Rényi Institute of Mathematics, Hungarian Academy of Scinces H-1053, Budapest,

More information

Efficiently Generating Random Bits from Finite State Markov Chains

Efficiently Generating Random Bits from Finite State Markov Chains 1 Efficienty Generating Random Bits from Finite State Markov Chains Hongchao Zhou and Jehoshua Bruck, Feow, IEEE Abstract The probem of random number generation from an uncorreated random source (of unknown

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

4 1-D Boundary Value Problems Heat Equation

4 1-D Boundary Value Problems Heat Equation 4 -D Boundary Vaue Probems Heat Equation The main purpose of this chapter is to study boundary vaue probems for the heat equation on a finite rod a x b. u t (x, t = ku xx (x, t, a < x < b, t > u(x, = ϕ(x

More information

Asynchronous Control for Coupled Markov Decision Systems

Asynchronous Control for Coupled Markov Decision Systems INFORMATION THEORY WORKSHOP (ITW) 22 Asynchronous Contro for Couped Marov Decision Systems Michae J. Neey University of Southern Caifornia Abstract This paper considers optima contro for a coection of

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

Statistical Learning Theory: A Primer

Statistical Learning Theory: A Primer Internationa Journa of Computer Vision 38(), 9 3, 2000 c 2000 uwer Academic Pubishers. Manufactured in The Netherands. Statistica Learning Theory: A Primer THEODOROS EVGENIOU, MASSIMILIANO PONTIL AND TOMASO

More information

Manipulation in Financial Markets and the Implications for Debt Financing

Manipulation in Financial Markets and the Implications for Debt Financing Manipuation in Financia Markets and the Impications for Debt Financing Leonid Spesivtsev This paper studies the situation when the firm is in financia distress and faces bankruptcy or debt restructuring.

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

Control with minimal cost-per-symbol encoding and quasi-optimality of event-based encoders

Control with minimal cost-per-symbol encoding and quasi-optimality of event-based encoders Control with minimal cost-per-symbol encoding and quasi-optimality of event-based encoders Justin Pearson, João P. Hespanha, Daniel Liberzon Abstract We consider the problem of stabilizing a continuous-time

More information

Asymptotic Properties of a Generalized Cross Entropy Optimization Algorithm

Asymptotic Properties of a Generalized Cross Entropy Optimization Algorithm 1 Asymptotic Properties of a Generaized Cross Entropy Optimization Agorithm Zijun Wu, Michae Koonko, Institute for Appied Stochastics and Operations Research, Caustha Technica University Abstract The discrete

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

14 Separation of Variables Method

14 Separation of Variables Method 14 Separation of Variabes Method Consider, for exampe, the Dirichet probem u t = Du xx < x u(x, ) = f(x) < x < u(, t) = = u(, t) t > Let u(x, t) = T (t)φ(x); now substitute into the equation: dt

More information

T.C. Banwell, S. Galli. {bct, Telcordia Technologies, Inc., 445 South Street, Morristown, NJ 07960, USA

T.C. Banwell, S. Galli. {bct, Telcordia Technologies, Inc., 445 South Street, Morristown, NJ 07960, USA ON THE SYMMETRY OF THE POWER INE CHANNE T.C. Banwe, S. Gai {bct, sgai}@research.tecordia.com Tecordia Technoogies, Inc., 445 South Street, Morristown, NJ 07960, USA Abstract The indoor power ine network

More information

FRIEZE GROUPS IN R 2

FRIEZE GROUPS IN R 2 FRIEZE GROUPS IN R 2 MAXWELL STOLARSKI Abstract. Focusing on the Eucidean pane under the Pythagorean Metric, our goa is to cassify the frieze groups, discrete subgroups of the set of isometries of the

More information

arxiv:quant-ph/ v3 6 Jan 1995

arxiv:quant-ph/ v3 6 Jan 1995 arxiv:quant-ph/9501001v3 6 Jan 1995 Critique of proposed imit to space time measurement, based on Wigner s cocks and mirrors L. Diósi and B. Lukács KFKI Research Institute for Partice and Nucear Physics

More information

8 Digifl'.11 Cth:uits and devices

8 Digifl'.11 Cth:uits and devices 8 Digif'. Cth:uits and devices 8. Introduction In anaog eectronics, votage is a continuous variabe. This is usefu because most physica quantities we encounter are continuous: sound eves, ight intensity,

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

Mat 1501 lecture notes, penultimate installment

Mat 1501 lecture notes, penultimate installment Mat 1501 ecture notes, penutimate instament 1. bounded variation: functions of a singe variabe optiona) I beieve that we wi not actuay use the materia in this section the point is mainy to motivate the

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

IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS, VOL. 15, NO. 2, FEBRUARY

IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS, VOL. 15, NO. 2, FEBRUARY IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS, VOL. 5, NO. 2, FEBRUARY 206 857 Optima Energy and Data Routing in Networks With Energy Cooperation Berk Gurakan, Student Member, IEEE, OmurOze,Member, IEEE,

More information

The basic equation for the production of turbulent kinetic energy in clouds is. dz + g w

The basic equation for the production of turbulent kinetic energy in clouds is. dz + g w Turbuence in couds The basic equation for the production of turbuent kinetic energy in couds is de TKE dt = u 0 w 0 du v 0 w 0 dv + g w q 0 q 0 e The first two terms on the RHS are associated with shear

More information

1 Heat Equation Dirichlet Boundary Conditions

1 Heat Equation Dirichlet Boundary Conditions Chapter 3 Heat Exampes in Rectanges Heat Equation Dirichet Boundary Conditions u t (x, t) = ku xx (x, t), < x (.) u(, t) =, u(, t) = u(x, ) = f(x). Separate Variabes Look for simpe soutions in the

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

Reichenbachian Common Cause Systems

Reichenbachian Common Cause Systems Reichenbachian Common Cause Systems G. Hofer-Szabó Department of Phiosophy Technica University of Budapest e-mai: gszabo@hps.ete.hu Mikós Rédei Department of History and Phiosophy of Science Eötvös University,

More information

Explicit overall risk minimization transductive bound

Explicit overall risk minimization transductive bound 1 Expicit overa risk minimization transductive bound Sergio Decherchi, Paoo Gastado, Sandro Ridea, Rodofo Zunino Dept. of Biophysica and Eectronic Engineering (DIBE), Genoa University Via Opera Pia 11a,

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

Primal and dual active-set methods for convex quadratic programming

Primal and dual active-set methods for convex quadratic programming Math. Program., Ser. A 216) 159:469 58 DOI 1.17/s117-15-966-2 FULL LENGTH PAPER Prima and dua active-set methods for convex quadratic programming Anders Forsgren 1 Phiip E. Gi 2 Eizabeth Wong 2 Received:

More information

Codes between MBR and MSR Points with Exact Repair Property

Codes between MBR and MSR Points with Exact Repair Property Codes between MBR and MSR Points with Exact Repair Property 1 Toni Ernva arxiv:1312.5106v1 [cs.it] 18 Dec 2013 Abstract In this paper distributed storage systems with exact repair are studied. A construction

More information

THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS

THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS ECCM6-6 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Sevie, Spain, -6 June 04 THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS M. Wysocki a,b*, M. Szpieg a, P. Heström a and F. Ohsson c a Swerea SICOMP

More information

A unified framework for Regularization Networks and Support Vector Machines. Theodoros Evgeniou, Massimiliano Pontil, Tomaso Poggio

A unified framework for Regularization Networks and Support Vector Machines. Theodoros Evgeniou, Massimiliano Pontil, Tomaso Poggio MASSACHUSETTS INSTITUTE OF TECHNOLOGY ARTIFICIAL INTELLIGENCE LABORATORY and CENTER FOR BIOLOGICAL AND COMPUTATIONAL LEARNING DEPARTMENT OF BRAIN AND COGNITIVE SCIENCES A.I. Memo No. 1654 March23, 1999

More information

Unconditional security of differential phase shift quantum key distribution

Unconditional security of differential phase shift quantum key distribution Unconditiona security of differentia phase shift quantum key distribution Kai Wen, Yoshihisa Yamamoto Ginzton Lab and Dept of Eectrica Engineering Stanford University Basic idea of DPS-QKD Protoco. Aice

More information

Do Schools Matter for High Math Achievement? Evidence from the American Mathematics Competitions Glenn Ellison and Ashley Swanson Online Appendix

Do Schools Matter for High Math Achievement? Evidence from the American Mathematics Competitions Glenn Ellison and Ashley Swanson Online Appendix VOL. NO. DO SCHOOLS MATTER FOR HIGH MATH ACHIEVEMENT? 43 Do Schoos Matter for High Math Achievement? Evidence from the American Mathematics Competitions Genn Eison and Ashey Swanson Onine Appendix Appendix

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

arxiv: v2 [cond-mat.stat-mech] 14 Nov 2008

arxiv: v2 [cond-mat.stat-mech] 14 Nov 2008 Random Booean Networks Barbara Drosse Institute of Condensed Matter Physics, Darmstadt University of Technoogy, Hochschustraße 6, 64289 Darmstadt, Germany (Dated: June 27) arxiv:76.335v2 [cond-mat.stat-mech]

More information

Reliability: Theory & Applications No.3, September 2006

Reliability: Theory & Applications No.3, September 2006 REDUNDANCY AND RENEWAL OF SERVERS IN OPENED QUEUING NETWORKS G. Sh. Tsitsiashvii M.A. Osipova Vadivosto, Russia 1 An opened queuing networ with a redundancy and a renewa of servers is considered. To cacuate

More information

Efficient Generation of Random Bits from Finite State Markov Chains

Efficient Generation of Random Bits from Finite State Markov Chains Efficient Generation of Random Bits from Finite State Markov Chains Hongchao Zhou and Jehoshua Bruck, Feow, IEEE Abstract The probem of random number generation from an uncorreated random source (of unknown

More information

Rate-Distortion Theory of Finite Point Processes

Rate-Distortion Theory of Finite Point Processes Rate-Distortion Theory of Finite Point Processes Günther Koiander, Dominic Schuhmacher, and Franz Hawatsch, Feow, IEEE Abstract We study the compression of data in the case where the usefu information

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

Pattern Frequency Sequences and Internal Zeros

Pattern Frequency Sequences and Internal Zeros Advances in Appied Mathematics 28, 395 420 (2002 doi:10.1006/aama.2001.0789, avaiabe onine at http://www.ideaibrary.com on Pattern Frequency Sequences and Interna Zeros Mikós Bóna Department of Mathematics,

More information

FFTs in Graphics and Vision. Spherical Convolution and Axial Symmetry Detection

FFTs in Graphics and Vision. Spherical Convolution and Axial Symmetry Detection FFTs in Graphics and Vision Spherica Convoution and Axia Symmetry Detection Outine Math Review Symmetry Genera Convoution Spherica Convoution Axia Symmetry Detection Math Review Symmetry: Given a unitary

More information

Throughput Optimal Scheduling for Wireless Downlinks with Reconfiguration Delay

Throughput Optimal Scheduling for Wireless Downlinks with Reconfiguration Delay Throughput Optima Scheduing for Wireess Downinks with Reconfiguration Deay Vineeth Baa Sukumaran vineethbs@gmai.com Department of Avionics Indian Institute of Space Science and Technoogy. Abstract We consider

More information

Nonlinear Output-Feedback Model Predictive Control with Moving Horizon Estimation

Nonlinear Output-Feedback Model Predictive Control with Moving Horizon Estimation Noninear Output-Feedback Mode Predictive Contro with Moving Horizon Estimation David A. Copp and João P. Hespanha Abstract We introduce an output-feedback approach to mode predictive contro that combines

More information

Coded Caching for Files with Distinct File Sizes

Coded Caching for Files with Distinct File Sizes Coded Caching for Fies with Distinct Fie Sizes Jinbei Zhang iaojun Lin Chih-Chun Wang inbing Wang Department of Eectronic Engineering Shanghai Jiao ong University China Schoo of Eectrica and Computer Engineering

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

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

Akaike Information Criterion for ANOVA Model with a Simple Order Restriction

Akaike Information Criterion for ANOVA Model with a Simple Order Restriction Akaike Information Criterion for ANOVA Mode with a Simpe Order Restriction Yu Inatsu * Department of Mathematics, Graduate Schoo of Science, Hiroshima University ABSTRACT In this paper, we consider Akaike

More information

Transcendence of stammering continued fractions. Yann BUGEAUD

Transcendence of stammering continued fractions. Yann BUGEAUD Transcendence of stammering continued fractions Yann BUGEAUD To the memory of Af van der Poorten Abstract. Let θ = [0; a 1, a 2,...] be an agebraic number of degree at east three. Recenty, we have estabished

More information

(1 ) = 1 for some 2 (0; 1); (1 + ) = 0 for some > 0:

(1 ) = 1 for some 2 (0; 1); (1 + ) = 0 for some > 0: Answers, na. Economics 4 Fa, 2009. Christiano.. The typica househod can engage in two types of activities producing current output and studying at home. Athough time spent on studying at home sacrices

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

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

VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES

VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES SARAH DAY, JEAN-PHILIPPE LESSARD, AND KONSTANTIN MISCHAIKOW Abstract. One of the most efficient methods for determining the equiibria of a continuous parameterized

More information

Small generators of function fields

Small generators of function fields Journa de Théorie des Nombres de Bordeaux 00 (XXXX), 000 000 Sma generators of function fieds par Martin Widmer Résumé. Soit K/k une extension finie d un corps goba, donc K contient un éément primitif

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

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

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, we as various

More information

Schedulability Analysis of Deferrable Scheduling Algorithms for Maintaining Real-Time Data Freshness

Schedulability Analysis of Deferrable Scheduling Algorithms for Maintaining Real-Time Data Freshness 1 Scheduabiity Anaysis of Deferrabe Scheduing Agorithms for Maintaining Rea-Time Data Freshness Song Han, Deji Chen, Ming Xiong, Kam-yiu Lam, Aoysius K. Mok, Krithi Ramamritham UT Austin, Emerson Process

More information

The EM Algorithm applied to determining new limit points of Mahler measures

The EM Algorithm applied to determining new limit points of Mahler measures Contro and Cybernetics vo. 39 (2010) No. 4 The EM Agorithm appied to determining new imit points of Maher measures by Souad E Otmani, Georges Rhin and Jean-Marc Sac-Épée Université Pau Veraine-Metz, LMAM,

More information

A Statistical Framework for Real-time Event Detection in Power Systems

A Statistical Framework for Real-time Event Detection in Power Systems 1 A Statistica Framework for Rea-time Event Detection in Power Systems Noan Uhrich, Tim Christman, Phiip Swisher, and Xichen Jiang Abstract A quickest change detection (QCD) agorithm is appied to the probem

More information

A NOTE ON QUASI-STATIONARY DISTRIBUTIONS OF BIRTH-DEATH PROCESSES AND THE SIS LOGISTIC EPIDEMIC

A NOTE ON QUASI-STATIONARY DISTRIBUTIONS OF BIRTH-DEATH PROCESSES AND THE SIS LOGISTIC EPIDEMIC (January 8, 2003) A NOTE ON QUASI-STATIONARY DISTRIBUTIONS OF BIRTH-DEATH PROCESSES AND THE SIS LOGISTIC EPIDEMIC DAMIAN CLANCY, University of Liverpoo PHILIP K. POLLETT, University of Queensand Abstract

More information

Chapter 7 PRODUCTION FUNCTIONS. Copyright 2005 by South-Western, a division of Thomson Learning. All rights reserved.

Chapter 7 PRODUCTION FUNCTIONS. Copyright 2005 by South-Western, a division of Thomson Learning. All rights reserved. Chapter 7 PRODUCTION FUNCTIONS Copyright 2005 by South-Western, a division of Thomson Learning. A rights reserved. 1 Production Function The firm s production function for a particuar good (q) shows the

More information

Arbitrary Throughput Versus Complexity Tradeoffs in Wireless Networks using Graph Partitioning

Arbitrary Throughput Versus Complexity Tradeoffs in Wireless Networks using Graph Partitioning University of Pennsyvania SchoaryCommons Departmenta Papers (ESE) Department of Eectrica & Systems Engineering November 2006 Arbitrary Throughput Versus Compexity Tradeoffs in Wireess Networks using Graph

More information

STABLE GRAPHS BENJAMIN OYE

STABLE GRAPHS BENJAMIN OYE STABLE GRAPHS BENJAMIN OYE Abstract. In Reguarity Lemmas for Stabe Graphs [1] Maiaris and Sheah appy toos from mode theory to obtain stronger forms of Ramsey's theorem and Szemeredi's reguarity emma for

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

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

Chemical Kinetics Part 2

Chemical Kinetics Part 2 Integrated Rate Laws Chemica Kinetics Part 2 The rate aw we have discussed thus far is the differentia rate aw. Let us consider the very simpe reaction: a A à products The differentia rate reates the rate

More information

Homework 5 Solutions

Homework 5 Solutions Stat 310B/Math 230B Theory of Probabiity Homework 5 Soutions Andrea Montanari Due on 2/19/2014 Exercise [5.3.20] 1. We caim that n 2 [ E[h F n ] = 2 n i=1 A i,n h(u)du ] I Ai,n (t). (1) Indeed, integrabiity

More information

A Fundamental Storage-Communication Tradeoff in Distributed Computing with Straggling Nodes

A Fundamental Storage-Communication Tradeoff in Distributed Computing with Straggling Nodes A Fundamenta Storage-Communication Tradeoff in Distributed Computing with Stragging odes ifa Yan, Michèe Wigger LTCI, Téécom ParisTech 75013 Paris, France Emai: {qifa.yan, michee.wigger} @teecom-paristech.fr

More information