Analysis of the average depth in a suffix tree under a Markov model

Size: px
Start display at page:

Download "Analysis of the average depth in a suffix tree under a Markov model"

Transcription

1 Analysis of the average depth in a sffix tree nder a Markov model Jlien Fayolle, Mark Daniel Ward To cite this version: Jlien Fayolle, Mark Daniel Ward. Analysis of the average depth in a sffix tree nder a Markov model. Conrado Martínez International Conference on Analysis of Algorithms, 2005, Barcelona, Spain. Discrete Mathematics and Theoretical Compter Science, DMTCS Proceedings vol. AD, International Conference on Analysis of Algorithms, pp.95-04, 2005, DMTCS Proceedings. <hal > HAL Id: hal Sbmitted on 2 Ag 205 HAL is a mlti-disciplinary open access archive for the deposit and dissemination of scientific research docments, hether they are pblished or not. The docments may come from teaching and research instittions in France or abroad, or from pblic or private research centers. L archive overte plridisciplinaire HAL, est destinée a dépôt et à la diffsion de docments scientifiqes de nivea recherche, pbliés o non, émanant des établissements d enseignement et de recherche français o étrangers, des laboratoires pblics o privés.

2 2005 International Conference on Analysis of Algorithms DMTCS proc. AD, 2005, Analysis of the average depth in a sffix tree nder a Markov model Jlien Fayolle and Mark Daniel Ward 2 Projet Algorithmes, INRIA, F-7853 Rocqencort, France. jlien.fayolle@inria.fr 2 Department of Mathematics, Prde University, West Lafayette, IN, USA. mard@math.prde.ed In this report, e prove that nder a Markovian model of order one, the average depth of sffix trees of index n is asymptotically similar to the average depth of tries (a.k.a. digital trees bilt on n independent strings. This leads to an asymptotic behavior of (log n/h + C for the average of the depth of the sffix tree, here h is the entropy of the Markov model and C is constant. Or proof compares the generating fnctions for the average depth in tries and in sffix trees; the difference beteen these generating fnctions is shon to be asymptotically small. We conclde by sing the asymptotic behavior of the average depth in a trie nder the Markov model fond by Jacqet and Szpankoski ([4]. Keyords: Sffix trees, depth, average analysis, asymptotics, analytic methods Introdction The sffix tree is a data strctre that lies at the core of pattern matching, sed for example in the lossless data compression algorithm of Lempel and Ziv [7]. Sffix trees are also tilized in bio-informatics to track significant patterns. A thorogh srvey of the se of sffix trees in compter science can be fond in Apostolico ([]. The average depth in a sffix tree of index n + is the average time (nmber of letters read necessary to insert a ne sffix into a tree of index n; an alternate interpretation is the average time to discriminate beteen the crrent sffix and any previos one. A device called a sorce emits letters randomly, and independently of their emission date. In this report, e focs on Markov sorces of order one: the letter emitted at a given time is a random variable obeying a Markov dependency of order one i.e., the distribtion of each letter depends only on the letter actally emitted immediately beforehand. Increasing the order introdces no ne technical challenges comptations only become more intricate and or proof for a Markovian dependency of order can easily be extended to any sorce ith greater Markovian dependency. We assme that the sorce is stationary throghot this report. We denote the stationary probability vector as π = (π i, the transition probability matrix as P = (p i,j, and the probability the sorce starts emitting a text ith the ord as P(. We also introdce the stationary probability matrix Π hose ros are the stationary probability vector π. We make the sal assmptions that the nderlying Markov chain is irredcible and aperiodic. This paper focses on the asymptotic behavior of the average depth in sffix trees. There are to challenging parts in this stdy: first the sffixes on hich the sffix tree is bilt are mtally dependent. For example, if e jst fond the pattern = 0000 in the text, then it sffices to have the letter 0 next, in order to find once more in the text. The trie (a.k.a. digital tree is mch easier to analyze than the sffix tree, becase the strings in a trie are independent from each other. The second challenge lies in the probabilistic dependency beteen symbols (Markov model. The probability that a pattern occrs in the text depends not only on the pattern itself bt also on hat as previosly seen in the text. In Jacqet and Szpankoski ([4], inclsion-exclsion as sed to obtain the asymptotics of the average depth in an (independent trie ith nderlying Markovian sorce. Therefore, it sffices for s to compare the asymptotic average depths in sffix trees against that of independent tries (here the probability model is Markovian in both cases. We prove that, for a Markovian sorce of order one, the average depth of a sffix tree of index n has a similar asymptotic behavior as a trie bilt on n strings. In section 2, e present the main reslts of this paper. We give the precise asymptotics for the average depth in a sffix tree ith an nderlying Markovian sorce. Afterards, e give a sketch of the proof. In section 3 e prove that the atocorrelation polynomial S (z is approximately ith high probability c 2005 Discrete Mathematics and Theoretical Compter Science (DMTCS, Nancy, France

3 96 Jlien Fayolle and Mark Daniel Ward (for of large length. To prove that sffix trees and independent tries have similar average depths, e first derive bivariate generating fnctions for D n and D t n in section 4. Then in sbsections 5. and 5.2 e analyze the difference beteen the generating fnctions for D n and D t n by tilizing complex asymptotics. Ultimately, e conclde in sbsection 5.3 that the depth in a sffix tree has the same average p to the constant term as the depth in a trie. Or method is motivated by the analytic approach of Jacqet and Szpankoski ([3] for the sffix tree nder a memoryless (a.k.a. Bernolli sorce model. Or reslts are proved in the more general context of Markovian sorces. 2 Main Reslts Consider a tree (say it s binary and it shall be so throghot this paper ith n leaves nmbered from to n. The depth D n (i of the leaf nmber i is defined as the length of the path from the root to this leaf. Pick one of the n leaves at random (niformly; the typical depth D n of the tree is the depth of the picked leaf. The random variable D n informs s abot the typical profile of the tree, rather than its height (i.e., the maximm depth. A trie is defined recrsively on a finite set S of infinite ords on A = {0, } as if S = 0, trie(s = if S =,, trie(s\0, trie(s\ else, here S\α represents the set of ords starting ith the letter α hose first letter has been removed. In order to bild a sffix tree, one needs as inpt an infinite string T on A called text and an integer n. We rite T = T T 2 T 3... and then T (i = T i T i+ T i+2... denotes the ith sffix of T (for example, T = T (, namely the entire string itself, is the first sffix. The sffix tree of index n bilt on T is defined as the trie bilt on the set of the n first sffixes of T (namely, T (,..., T (n. In the context of sffix trees, D n (i is interpreted as the largest vale of k sch that there exists j i ith j n for hich T i+k i = T j+k j. Or discssion in this paper primarily concerns the comparison of the average depths in sffix trees verss independent tries. Throghot this discssion, D n alays denotes the typical depth in a sffix tree. In a trie bilt over n independent strings, e let D t n denote the typical depth. In [4], the asymptotic behavior of the average depth for a trie bilt on n independent strings nder a Markovian sorce of order as established sing inclsion-exclsion. In this paper, e prove the analogos reslt for sffix trees. Namely, e establish the asymptotic average depth for a sffix tree ith an nderlying Markovian sorce. Or proof techniqe consists of shoing that D n and D t n have similar asymptotic distribtions. To do this, e estimate the difference of their probability generating fnctions, and sho that the difference is asymptotically small. In order to prove that D n and D t n have asymptotically similar averages, e first discss the atocorrelation of a string. Since a string can overlap ith itself, the bivariate generating fnction D(z, for D n, here marks the depth and z the size, incldes the atocorrelation polynomial S (z. Fortnately, ith high probability, a random string has very little overlap ith itself. Therefore the atocorrelation polynomial of is close to ith high probability. After discssing atocorrelation, e present the bivariate generating fnctions for both D n and D t n, hich e denote by D(z, and D t (z,, respectively. We have D(z, = (z P( D (z 2 and D t (z, = A A zp( ( z + zp( 2, ( here D(z is a generating fnction. Or goal is to prove that the to generating fnctions are asymptotically very similar; it follos from there that Dn t and D n have the same average asymptotically p to the constant term. In order to determine the asymptotics of the difference D(z, D t (z,, e tilize complex analysis. From (, e see that D t (z, has exactly one pole of order 2 for each. Using Roché s theorem, e prove that, for z ρ (ith ρ > and for sfficiently large, the generating fnction D (z has exactly one dominant root of order 2 for each. Therefore D(z, has exactly one dominant pole of order 2 for each. We next se Cachy s theorem and singlarity analysis to qantify the contribtion

4 Analysis of the average depth in a sffix tree nder a Markov model 97 from each of these poles to the difference Q n ( := ( [z n ](D(z, D t (z,. Or analysis of the difference Q n ( ltimately relies on the Mellin transform. We conclde from there that the averages of the depths D n (in sffix trees and Dn t (in independent tries are asymptotically similar. Therefore, D n has mean h log n pls some flctations. Specifically, Theorem For a sffix tree bilt on the first n sffixes of a text prodced by a sorce nder the Markovian model of order one, the average typical depth is asymptotically E(D n = ( log n + γ + h 2 h 2h H + P (log n + O ( n ɛ, (2 for some positive ɛ, here γ is the Eler constant, h is the entropy of the Markov sorce, h 2 is the second order entropy, H is the stationary entropy and P is a fnction flctating arond zero ith a very small amplitde. In particlar h := π i p i,j log p i,j and H := π i log π i. (3 i,j A 2 i A 2 3 Atocorrelation The depth of the ith leaf in a sffix tree can be constred in the langage of pattern matching as the length of the longest factor starting at a position i in the text that can be seen at least once more in the text. Hence D n (i k means that there is at least another pattern in the text T starting ith T i+k i While looking for the longest pattern in T that matches the text starting in position i, there is a possibility that to occrrences of this longest pattern overlap; this possible overlap of the pattern ith itself cases complications in the enmeration of occrrences of in a text. The phenomenon exhibited hen overlaps ith itself is called atocorrelation. To accont for this overlapping, e se the probabilistic version of the atocorrelation polynomial. For a pattern of size k, the polynomial is defined as: k S (z := [ k i = i+ k ]P(k i+ k k i z i. (4 i=0 For simplicity, e introdce c i = [[ k i = i+ k ] here [.] is Iverson s bracket notation for the indicator fnction. We say there is an overlap of size k i if c i =. This means that the sffix and the prefix of size k i of coincide. Graphically, an overlap of a pattern (hite rectangles looks like this: here the to black boxes are the matching sffix and prefix. For example = has overlaps of sizes 3, 6 and 9; note that (here = 9 is alays a valid overlap since the pattern alays matches itself. We define a period of a pattern as any integer i for hich c i = (so a pattern often has several periods, the minimal period m( of is the smallest positive period, if one exists (if has no positive period for, e define m( = k. We no formlate precisely the intition that, for most patterns of size k, the atocorrelation polynomial is very close to. It stems from the fact that the sm of the probabilities P( over all loosely correlated patterns (patterns ith a large minimal period of a given size is very close to. Lemma There exist δ <, ρ > ith ρδ <, and θ > 0, sch that for any integer k A k [ S (ρ (ρδ k θ ]P( θδ k. (5 Proof: Note that S (z has a term of degree i ith i k if and only if c i =. If c i = the prefix of size i of ill repeat itself flly in as many times as there is i in k, hence knoing c i = and the first i letters of allos s to describe flly the ord. Therefore, given a fixed... i, there,

5 98 Jlien Fayolle and Mark Daniel Ward is exactly one ord i+... k sch that the polynomial S (z has minimal degree i k. We let p i,j = p i, j, π i = π i, and p = max ( p i,j, π i. Then, for fixed j and k, i,j 2 j i= A k [m( = i]p( = j i= j i=,..., i A i π p,2 p i,i,..., i A i π p,2 p 2,3 p i,j p k i, i+,..., k A k i [m( = i] p i,i+ p k,k bt e can factor p k i otside the inner sm, since it does not depend on,..., i. Next, e observe that,..., i A i π p,2 p 2,3 p i,i = ΠP i =, so A k [S (z has minimal degree j ]P( and this holds hen j = k/2. So j i= p k i pk j p, A k [all terms of S (z have degree > k/2 ]P( p k/2 p = θδk. (7 Remark that, if all terms of S (z have degree > k/2, then (6 k S (ρ i= k/2 (ρp i ρ k p k/2 + p = (ρδ k θ. (8 We select δ = p, θ = ( p and some ρ > ith δρ < to complete the proof of the lemma. The next lemma proves that, for sfficiently large and for some radis ρ >, the atocorrelation polynomial does not vanish on the disk of radis ρ. Lemma 2 There exist K, ρ > ith pρ <, and α > 0 sch that for any pattern of size larger than K and z in a disk of radis ρ, e have S (z α. Proof: Like for the previos lemma, e split the proof into to cases, according to the index i of the minimal period of the pattern of size k. Since the atocorrelation polynomial alays has c 0 =, e rite k S (z = + c j P(k j+ k k j z j. (9 j=i We introdce ρ > sch that pρ <. Therefore, if i > k/2, then k S (z c j P(k j+ k k j z j (pρ i pρ, (0 j=i in z ρ. Bt since i > k/2 and pρ <, e get S (z (pρ k/2 pρ. ( We observe that, for some K sfficiently large, any pattern of size larger than K satifies S (z α and this loer bond α is positive. We recall that if c i =, the prefix of size i of ill repeat itself flly in as many times as there is i in k i.e., q := k/i times, the remainder ill be the prefix v of of size r := k k/i i, hence = k/i v. We also introdce the ord v sch that vv = (of length t := i r = i k + k/i i.

6 Analysis of the average depth in a sffix tree nder a Markov model 99 If i k/2, e make the atocorrelation polynomial explicit: S (z = + P(v v k z i + P(v v k z 2i + + P(v q 2 v k z i(q S v (z. (2 Overall, e can rite P(v j v k = A j+ here A = p k,v p v,v 2... p v t,v... p v r,v r is the prodct of i transition probabilities, bt since k = v r e obtain S (z = + Az i + (Az i (Az i q S v (z = (Azi q Az i + (Az i q S v (z. (3 Then e provide a loer-bond for S (z : S (z (Az i q Az i (Az i q S v (z (pρ i(q + (pρ i (pρ i(q S v (z (pρ i(q + (pρ i (pρ i(q pρ. We have pρ < so that (pρ k tends to zero ith k. Since i(q is close to k (at orse k/3 if = v for some K 2 sfficiently large and patterns of size larger than K 2, only the term ( + (pρ i remains. Finally, e set K = max{k, K 2 }. 4 On the Generating Fnctions The probabilistic model for the random variable D n is the prodct of a Markov model of order one for the sorce generating the strings (one string in the case of the sffix tree, n for the trie, and a niform model on {,..., n} for choosing the leaf. Hence, if X is a random variable niformly distribted over {,..., n}, and T a random text generated by the sorce, the typical depth is D n := n D n (i(t [X = i]. i= For the rest of the paper, the t exponent on a qantity ill indicate its trie version. Or aim is to compare asymptotically the probability generating fnctions of the depth for a sffix tree (namely, D n ( := k P(D n = k k and a trie (namely, D t n( := k P(Dt n = k k. We provide in this section an explicit expression for these generating fnctions and their respective bivariate extensions D(z, := n nd n(z n and D t (z, := n ndt n(z n. We first derive D t n(. Each string from the set of strings S is associated to a niqe leaf in the trie. By definition of the trie, the letters read on the path going from the root of the trie to a leaf form the smallest prefix distingishing one string from the n others. We choose niformly a leaf among the n leaves of the trie. Let A k denote the prefix of length k of the string associated to this randomly selected leaf. We say that D t n < k if and only if the other n texts do not have as a prefix. It follos immediately that P(D t n(i < k = A k P(( P( n, conseqently D t n( = A P(( P( n and D t (z, = A for < and z <. The sffix tree generating fnction is knon from [5] and [3]: for < and z < D(z, = zp( ( z + zp( 2, (4 (z P( D (z 2, (5 A here D (z = ( zs (z + z P(( + ( zf (z and for z < P Π, F (z := Π π n 0(P n+ z n = [(P Π(I (P Πz π m, ] m,. (6

7 00 Jlien Fayolle and Mark Daniel Ward 5 Asymptotics 5. Isolating the dominant pole We prove first that for a pattern of size large enogh there is a single dominant root to D (z. Then e sho that there is a disk of radis greater than containing each single dominant root of the D (z s for any of size big enogh bt no other root of the D (z s. Lemma 3 There exists a radis ρ > and an integer K sch that for any of size larger than K, D (z has only one root in the disk of radis ρ. Proof: Let be a given pattern. We apply Roché s Theorem to sho the niqeness of the smallest modls root of D (z. The main condition e need to flfill is that, on a given circle z = ρ, f(z := ( zs (z > P(z k ( + ( zf (z =: g(z. (7 The fnction f is analytic since it is a polynomial, F is analytic for z < P Π (here P Π >, so g is too. For patterns of a size large enogh, P(z k ill be small enogh to obtain the desired condition. The main isse is the bonding from above of F (z on the circle of radis ρ. We note d = min a A π a ; this vale is positive (otherise a letter old never occr and since (P Π n+ = P n+ Π (remember that P Π = ΠP = Π and ΠΠ = Π, e have: F (z d n 0(P n+ Πz n [P n+ ] k, [Π] k, z n (8 d br n+ ρ n br d d n 0 k, n 0 rρ, (9 here b and r are constants (independent of the pattern ith 0 < r < and ρ sch that rρ <. Let K be an integer and ρ some radis satisfying Lemma 2, e ask ρ to be smaller than ρ so that pρ < and S (z α on z = ρ and for K. There exists K large enogh sch that for any k > K e verify the condition (pρ k ( + ( + ρ br d < α(ρ. (20 rρ On a disk of sch radis ρ and for k > max{k, K }, the assmptions of Roché s Theorem are satisfied, conseqently (f + g(z = D (z has exactly as many zeros in the centered disk of radis ρ as f(z, namely one zero, since S (z does not vanish. Frthermore the assmptions of Roché s Theorem are satisfied on z = ρ for any pattern of size larger than K := max{k, K }. So for any ith K, D (z has exactly one root ithin the disk z = ρ. 5.2 Compting resides The se of Roché s Theorem has established the existence of a single zero of smallest modls for D (z, e denote it A. We kno by Pringsheim s Theorem that A is real positive. Let also B and C be the vales of the first and second derivatives of D (z at z = A. We make se of a bootstrapping techniqe to find an expansion for A, B, and C along the poers of P( and e obtain A = + P( S ( + O(P2 (, [ B = S ( + P( k F ( 2 S ( S ( C = 2S ( + P( ] + O(P 2 (, and [ 3 S ( S ( + k(k 2F ( 2kF ( ] + O(P 2 (. We no compare D n ( and Dn( t to conclde that they are asymptotically close. We therefore introdce to ne generating fnctions ( z Q n ( := D 2 (z (D n( Dn( t and Q(z, := nq n (z n = P( n 0 A (2 z ( z( P( 2.

8 Analysis of the average depth in a sffix tree nder a Markov model 0 We apply Cachy s Theorem to Q(z, ith z rnning along the circle centered at the origin hose radis ρ as determined in 5.. There are only three singlarities ithin this contor: at z = 0, at A, and at ( P(. In order to jstify the presence of the third singlarity ithin the circle, e note that the condition (20 implies P(ρ < P(ρ k < (pρ k ( + ( + ρ br d < α(ρ < (ρ, (22 rρ since α is taken smaller than one. Ths ( P( is smaller than the radis ρ. For any of a size larger than K, e have I (ρ, := ( z P( 2iπ z =ρ z n+ D 2 (z z ( z( P( 2 dz = P(f( ( = Res(f(z; 0 + Res(f(z; A + Res f(z; P( ( ( ( = nq n ( + P( Res + Res z z n+ D 2 (z ; A z z n+ ( z( P( 2 ;. P( (23 Since z/( z( P( 2 is analytic at z = A it does not contribte to the reside of f(z at A and can be discarded in the comptation of the reside. The same is tre for the /D 2 (z part in the reside of f(z at z = /( P(. We compte the reside at A sing the expansion e fond for B and C throgh bootstrapping. We set k = to simplify the notation, and by a Taylor expansion near A e obtain ( z k (n+ Res D 2 (z ; A = A n ( (n + B 2 A C B 3. (24 In order to compte the reside at z = ( P(, e se the general formla [z n ]( az 2 = (n + a n ths [z ] z n ( z( P( 2 = [zn ] ( z( P( 2 = n( P(n. (25 Before e proceed to get the asymptotic behavior of Q n (, the folloing technical lemma is very sefl. Lemma 4 For any fnction f defined over the patterns on A, and any y e have A k P(f( y + f max P({ A k : f( > y}, here f max is the maximm of f over all patterns of size k. Proof: For the patterns of size k ith f( > y, f max is an pper bond of f(; the probability of the other patterns is smaller than. We also note that S (ρ ( pρ and D (z = O(ρ k for z ρ. Ths (details are omitted, e obtain the folloing bond for the sm of I (ρ, over all patterns of size k: I (ρ, = O((δρ k ρ n A k There are only finitely many patterns ith < K ; these terms provide a contribtion of at most O(B n to Q n ( for some B >. We have jst proved the folloing Lemma 5 For some β >, and for all β, there exists B > for hich e have Q n ( = ( ( (n + P( A n n B 2 C A A B 3 n( P( n + O(B n. (26

9 02 Jlien Fayolle and Mark Daniel Ward 5.3 Asymptotic behavior of Q n ( Lemma 6 For all β sch that < β < δ, there exists a positive ɛ sch that Q n ( = O(n ɛ niformly for all β. Proof: For n large enogh, the dominant term in eqation 26 is Q n ( = ( A n 2 P( B 2 ( P( n + O(n. (27 A We introdce the fnction [ ] [ A x 2 f (x = B 2 ( P( x A 2 B 2 ] exp( x, P( and perform a Mellin transform (see Flajolet, Gordon, and Dmas [2] for a thorogh overvie of this techniqe on P(f (x. It old have seemed natral to define f (x by its first term bt for simplicity s sake e force a O(x behavior for f (x near zero by sbtracting some vale. By an application of Lemma 4, ith y = ( δ k for some δ <, e prove that the sm is absoltely convergent for β. Hence the Mellin transform f (s, of the sm is defined and e have f (s = Γ(s ( (log A s A 2 B 2 ( log( P( s and f (s, = P( P(f (s. f (s, is analytical ithin an open strip < R(s < c for some positive c. In order to do so e split the sm over all patterns into to parts. For the patterns ith I(sP( small, e se an expansion of f(s and then apply Lemma 4. We make sre these patterns do not create any singlarity. For any given s there are only finitely many patterns ith I(sP( large so their contribtion to the sm f (s,, althogh each is individally large, do not create any singlarity. Therefore, e are able to take the inverse Mellin transform of f (s,, and since there is no pole in the strip e obtain the stated reslt. This last reslt is sfficient to prove or initial claim. By definition Dn( t = D n ( =, ths e have Q n ( = ( (D n( Dn( t Dn ( D n ( = Dt n( D t n(. (28 We have D n( = E(D n and (D t n ( = E(D t n, therefore hen tends to e obtain E(D t n E(D n = O(n ɛ. (29 This means that asymptotically the difference beteen the to averages is no larger than O(n ɛ for some positive ɛ. 6 Conclsion We have shon that the average depths of tries and sffix trees behave asymptotically likeise for a Markov model of order one. This reslt can be extended to any order of the Markov model. An extended analysis shold yield analogos reslts for both the variance and the limiting distribtion of typical depth. A normal distribtion is then expected for sffix trees. In the ftre, e also hope to extend or reslts to the more general probabilistic sorce model introdced by Vallée in [6]. Acknoledgements The athors thank Philippe Flajolet and Wojciech Szpankoski for serving as hosts dring the project. Both athors are thankfl to Wojciech Szpankoski for introdcing the problem, and also to Philippe Flajolet and Mireille Régnier for sefl discssions.

10 Analysis of the average depth in a sffix tree nder a Markov model 03 References [] Alberto Apostolico. The myriad virtes of sffix trees. In A. Apostolico and Z. Galil, editors, Combinatorial Algorithms on Words, volme 2 of NATO Advance Science Institte Series. Series F: Compter and Systems Sciences, pages Springer Verlag, 985. [2] Philippe Flajolet, Xavier Gordon, and Philippe Dmas. Mellin transforms and asymptotics: Harmonic sms. Theoretical Compter Science, 44( 2:3 58, Jne 995. [3] P. Jacqet and W. Szpankoski. Analytic approach to pattern matching. In M. Lothaire, editor, Applied Combinatorics on Words, chapter 7. Cambridge, [4] Philippe Jacqet and Wojciech Szpankoski. Analysis of digital tries ith Markovian dependency. IEEE Transactions on Information Theory, 37(5: , 99. [5] Mireille Régnier and Wojciech Szpankoski. On pattern freqency occrrences in a Markovian seqence. Algorithmica, 22(4:63 649, 998. [6] Brigitte Vallée. Dynamical sorces in information theory: Fndamental intervals and ord prefixes. Algorithmica, 29(/2: , 200. [7] Jacob Ziv and Abraham Lempel. A niversal algorithm for seqential data compression. IEEE Transactions on Information Theory, IT-23: , 977.

11 04 Jlien Fayolle and Mark Daniel Ward

Sources of Non Stationarity in the Semivariogram

Sources of Non Stationarity in the Semivariogram Sorces of Non Stationarity in the Semivariogram Migel A. Cba and Oy Leangthong Traditional ncertainty characterization techniqes sch as Simple Kriging or Seqential Gassian Simlation rely on stationary

More information

Constructive Root Bound for k-ary Rational Input Numbers

Constructive Root Bound for k-ary Rational Input Numbers Constrctive Root Bond for k-ary Rational Inpt Nmbers Sylvain Pion, Chee Yap To cite this version: Sylvain Pion, Chee Yap. Constrctive Root Bond for k-ary Rational Inpt Nmbers. 19th Annal ACM Symposim on

More information

On relative errors of floating-point operations: optimal bounds and applications

On relative errors of floating-point operations: optimal bounds and applications On relative errors of floating-point operations: optimal bonds and applications Clade-Pierre Jeannerod, Siegfried M. Rmp To cite this version: Clade-Pierre Jeannerod, Siegfried M. Rmp. On relative errors

More information

Analysis of the multiplicity matching parameter in suffix trees

Analysis of the multiplicity matching parameter in suffix trees 2005 International Conference on Analysis of Algorithms DMTCS proc. AD, 2005, 307 322 Analysis of the multiplicity matching parameter in suffix trees Mark Daniel Ward and Wojciech Szpankoski 2 Department

More information

Effects of Soil Spatial Variability on Bearing Capacity of Shallow Foundations

Effects of Soil Spatial Variability on Bearing Capacity of Shallow Foundations Geotechnical Safety and Risk V T. Schweckendiek et al. (Eds.) 2015 The athors and IOS Press. This article is pblished online with Open Access by IOS Press and distribted nder the terms of the Creative

More information

Upper Bounds on the Spanning Ratio of Constrained Theta-Graphs

Upper Bounds on the Spanning Ratio of Constrained Theta-Graphs Upper Bonds on the Spanning Ratio of Constrained Theta-Graphs Prosenjit Bose and André van Renssen School of Compter Science, Carleton University, Ottaa, Canada. jit@scs.carleton.ca, andre@cg.scs.carleton.ca

More information

Universal Scheme for Optimal Search and Stop

Universal Scheme for Optimal Search and Stop Universal Scheme for Optimal Search and Stop Sirin Nitinawarat Qalcomm Technologies, Inc. 5775 Morehose Drive San Diego, CA 92121, USA Email: sirin.nitinawarat@gmail.com Vengopal V. Veeravalli Coordinated

More information

Study on the impulsive pressure of tank oscillating by force towards multiple degrees of freedom

Study on the impulsive pressure of tank oscillating by force towards multiple degrees of freedom EPJ Web of Conferences 80, 0034 (08) EFM 07 Stdy on the implsive pressre of tank oscillating by force towards mltiple degrees of freedom Shigeyki Hibi,* The ational Defense Academy, Department of Mechanical

More information

Lecture Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2

Lecture Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2 BIJU PATNAIK UNIVERSITY OF TECHNOLOGY, ODISHA Lectre Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2 Prepared by, Dr. Sbhend Kmar Rath, BPUT, Odisha. Tring Machine- Miscellany UNIT 2 TURING MACHINE

More information

Spanning Trees with Many Leaves in Graphs without Diamonds and Blossoms

Spanning Trees with Many Leaves in Graphs without Diamonds and Blossoms Spanning Trees ith Many Leaes in Graphs ithot Diamonds and Blossoms Pal Bonsma Florian Zickfeld Technische Uniersität Berlin, Fachbereich Mathematik Str. des 7. Jni 36, 0623 Berlin, Germany {bonsma,zickfeld}@math.t-berlin.de

More information

CHARACTERIZATIONS OF EXPONENTIAL DISTRIBUTION VIA CONDITIONAL EXPECTATIONS OF RECORD VALUES. George P. Yanev

CHARACTERIZATIONS OF EXPONENTIAL DISTRIBUTION VIA CONDITIONAL EXPECTATIONS OF RECORD VALUES. George P. Yanev Pliska Std. Math. Blgar. 2 (211), 233 242 STUDIA MATHEMATICA BULGARICA CHARACTERIZATIONS OF EXPONENTIAL DISTRIBUTION VIA CONDITIONAL EXPECTATIONS OF RECORD VALUES George P. Yanev We prove that the exponential

More information

On relative errors of floating-point operations: optimal bounds and applications

On relative errors of floating-point operations: optimal bounds and applications On relative errors of floating-point operations: optimal bonds and applications Clade-Pierre Jeannerod, Siegfried M. Rmp To cite this version: Clade-Pierre Jeannerod, Siegfried M. Rmp. On relative errors

More information

The Lehmer matrix and its recursive analogue

The Lehmer matrix and its recursive analogue The Lehmer matrix and its recrsive analoge Emrah Kilic, Pantelimon Stănică TOBB Economics and Technology University, Mathematics Department 0660 Sogtoz, Ankara, Trkey; ekilic@etedtr Naval Postgradate School,

More information

Analysis of the Multiplicity Matching Parameter in Suffix Trees

Analysis of the Multiplicity Matching Parameter in Suffix Trees Analysis of the Multiplicity Matching Parameter in Suffix Trees Mark Daniel Ward 1 and Wojciech Szpankoski 2 1 Department of Mathematics, Purdue University, West Lafayette, IN, USA. mard@math.purdue.edu

More information

A Characterization of Interventional Distributions in Semi-Markovian Causal Models

A Characterization of Interventional Distributions in Semi-Markovian Causal Models A Characterization of Interventional Distribtions in Semi-Markovian Casal Models Jin Tian and Changsng Kang Department of Compter Science Iowa State University Ames, IA 50011 {jtian, cskang}@cs.iastate.ed

More information

Discontinuous Fluctuation Distribution for Time-Dependent Problems

Discontinuous Fluctuation Distribution for Time-Dependent Problems Discontinos Flctation Distribtion for Time-Dependent Problems Matthew Hbbard School of Compting, University of Leeds, Leeds, LS2 9JT, UK meh@comp.leeds.ac.k Introdction For some years now, the flctation

More information

Active Flux Schemes for Advection Diffusion

Active Flux Schemes for Advection Diffusion AIAA Aviation - Jne, Dallas, TX nd AIAA Comptational Flid Dynamics Conference AIAA - Active Fl Schemes for Advection Diffsion Hiroaki Nishikawa National Institte of Aerospace, Hampton, VA 3, USA Downloaded

More information

Introdction Finite elds play an increasingly important role in modern digital commnication systems. Typical areas of applications are cryptographic sc

Introdction Finite elds play an increasingly important role in modern digital commnication systems. Typical areas of applications are cryptographic sc A New Architectre for a Parallel Finite Field Mltiplier with Low Complexity Based on Composite Fields Christof Paar y IEEE Transactions on Compters, Jly 996, vol 45, no 7, pp 856-86 Abstract In this paper

More information

Frequency Estimation, Multiple Stationary Nonsinusoidal Resonances With Trend 1

Frequency Estimation, Multiple Stationary Nonsinusoidal Resonances With Trend 1 Freqency Estimation, Mltiple Stationary Nonsinsoidal Resonances With Trend 1 G. Larry Bretthorst Department of Chemistry, Washington University, St. Lois MO 6313 Abstract. In this paper, we address the

More information

An Investigation into Estimating Type B Degrees of Freedom

An Investigation into Estimating Type B Degrees of Freedom An Investigation into Estimating Type B Degrees of H. Castrp President, Integrated Sciences Grop Jne, 00 Backgrond The degrees of freedom associated with an ncertainty estimate qantifies the amont of information

More information

Non-Lecture I: Linear Programming. Th extremes of glory and of shame, Like east and west, become the same.

Non-Lecture I: Linear Programming. Th extremes of glory and of shame, Like east and west, become the same. The greatest flood has the soonest ebb; the sorest tempest the most sdden calm; the hottest love the coldest end; and from the deepest desire oftentimes enses the deadliest hate. Th extremes of glory and

More information

J.A. BURNS AND B.B. KING redced order controllers sensors/actators. The kernels of these integral representations are called fnctional gains. In [4],

J.A. BURNS AND B.B. KING redced order controllers sensors/actators. The kernels of these integral representations are called fnctional gains. In [4], Jornal of Mathematical Systems, Estimation, Control Vol. 8, No. 2, 1998, pp. 1{12 c 1998 Birkhaser-Boston A Note on the Mathematical Modelling of Damped Second Order Systems John A. Brns y Belinda B. King

More information

Setting The K Value And Polarization Mode Of The Delta Undulator

Setting The K Value And Polarization Mode Of The Delta Undulator LCLS-TN-4- Setting The Vale And Polarization Mode Of The Delta Undlator Zachary Wolf, Heinz-Dieter Nhn SLAC September 4, 04 Abstract This note provides the details for setting the longitdinal positions

More information

ON THE SHAPES OF BILATERAL GAMMA DENSITIES

ON THE SHAPES OF BILATERAL GAMMA DENSITIES ON THE SHAPES OF BILATERAL GAMMA DENSITIES UWE KÜCHLER, STEFAN TAPPE Abstract. We investigate the for parameter family of bilateral Gamma distribtions. The goal of this paper is to provide a thorogh treatment

More information

Imprecise Continuous-Time Markov Chains

Imprecise Continuous-Time Markov Chains Imprecise Continos-Time Markov Chains Thomas Krak *, Jasper De Bock, and Arno Siebes t.e.krak@.nl, a.p.j.m.siebes@.nl Utrecht University, Department of Information and Compting Sciences, Princetonplein

More information

Nonlinear parametric optimization using cylindrical algebraic decomposition

Nonlinear parametric optimization using cylindrical algebraic decomposition Proceedings of the 44th IEEE Conference on Decision and Control, and the Eropean Control Conference 2005 Seville, Spain, December 12-15, 2005 TC08.5 Nonlinear parametric optimization sing cylindrical algebraic

More information

Consistent Numerical Model for Wind Buffeting Analysis of Long-Span Bridges

Consistent Numerical Model for Wind Buffeting Analysis of Long-Span Bridges Consistent Nmerical Model for Wind Bffeting Analysis of Long-pan Bridges Dorian JANJIC Technical Director TDV GesmbH Graz, Astria Heinz PIRCHER Director TDV GesmbH Graz, Astria mmary The bffeting analysis

More information

4.2 First-Order Logic

4.2 First-Order Logic 64 First-Order Logic and Type Theory The problem can be seen in the two qestionable rles In the existential introdction, the term a has not yet been introdced into the derivation and its se can therefore

More information

Section 7.4: Integration of Rational Functions by Partial Fractions

Section 7.4: Integration of Rational Functions by Partial Fractions Section 7.4: Integration of Rational Fnctions by Partial Fractions This is abot as complicated as it gets. The Method of Partial Fractions Ecept for a few very special cases, crrently we have no way to

More information

A Theory of Markovian Time Inconsistent Stochastic Control in Discrete Time

A Theory of Markovian Time Inconsistent Stochastic Control in Discrete Time A Theory of Markovian Time Inconsistent Stochastic Control in Discrete Time Tomas Björk Department of Finance, Stockholm School of Economics tomas.bjork@hhs.se Agatha Mrgoci Department of Economics Aarhs

More information

Second-Order Wave Equation

Second-Order Wave Equation Second-Order Wave Eqation A. Salih Department of Aerospace Engineering Indian Institte of Space Science and Technology, Thirvananthapram 3 December 016 1 Introdction The classical wave eqation is a second-order

More information

Prandl established a universal velocity profile for flow parallel to the bed given by

Prandl established a universal velocity profile for flow parallel to the bed given by EM 0--00 (Part VI) (g) The nderlayers shold be at least three thicknesses of the W 50 stone, bt never less than 0.3 m (Ahrens 98b). The thickness can be calclated sing Eqation VI-5-9 with a coefficient

More information

Modelling by Differential Equations from Properties of Phenomenon to its Investigation

Modelling by Differential Equations from Properties of Phenomenon to its Investigation Modelling by Differential Eqations from Properties of Phenomenon to its Investigation V. Kleiza and O. Prvinis Kanas University of Technology, Lithania Abstract The Panevezys camps of Kanas University

More information

FOUNTAIN codes [3], [4] provide an efficient solution

FOUNTAIN codes [3], [4] provide an efficient solution Inactivation Decoding of LT and Raptor Codes: Analysis and Code Design Francisco Lázaro, Stdent Member, IEEE, Gianligi Liva, Senior Member, IEEE, Gerhard Bach, Fellow, IEEE arxiv:176.5814v1 [cs.it 19 Jn

More information

Concepts Introduced. Digital Electronics. Logic Blocks. Truth Tables

Concepts Introduced. Digital Electronics. Logic Blocks. Truth Tables Concepts Introdced Digital Electronics trth tables, logic eqations, and gates combinational logic seqential logic Digital electronics operate at either high or low voltage. Compters se a binary representation

More information

Prediction of Transmission Distortion for Wireless Video Communication: Analysis

Prediction of Transmission Distortion for Wireless Video Communication: Analysis Prediction of Transmission Distortion for Wireless Video Commnication: Analysis Zhifeng Chen and Dapeng W Department of Electrical and Compter Engineering, University of Florida, Gainesville, Florida 326

More information

CHANNEL SELECTION WITH RAYLEIGH FADING: A MULTI-ARMED BANDIT FRAMEWORK. Wassim Jouini and Christophe Moy

CHANNEL SELECTION WITH RAYLEIGH FADING: A MULTI-ARMED BANDIT FRAMEWORK. Wassim Jouini and Christophe Moy CHANNEL SELECTION WITH RAYLEIGH FADING: A MULTI-ARMED BANDIT FRAMEWORK Wassim Joini and Christophe Moy SUPELEC, IETR, SCEE, Avene de la Bolaie, CS 47601, 5576 Cesson Sévigné, France. INSERM U96 - IFR140-

More information

Predicting Popularity of Twitter Accounts through the Discovery of Link-Propagating Early Adopters

Predicting Popularity of Twitter Accounts through the Discovery of Link-Propagating Early Adopters Predicting Poplarity of Titter Acconts throgh the Discoery of Link-Propagating Early Adopters Daichi Imamori Gradate School of Informatics, Kyoto Uniersity Sakyo, Kyoto 606-850 Japan imamori@dl.soc.i.kyoto-.ac.jp

More information

Counting Markov Types

Counting Markov Types Counting Markov Types Philippe Jacquet, Charles Knessl, Wojciech Szpankowski To cite this version: Philippe Jacquet, Charles Knessl, Wojciech Szpankowski. Counting Markov Types. Drmota, Michael and Gittenberger,

More information

Experimental Study of an Impinging Round Jet

Experimental Study of an Impinging Round Jet Marie Crie ay Final Report : Experimental dy of an Impinging Rond Jet BOURDETTE Vincent Ph.D stdent at the Rovira i Virgili University (URV), Mechanical Engineering Department. Work carried ot dring a

More information

Ted Pedersen. Southern Methodist University. large sample assumptions implicit in traditional goodness

Ted Pedersen. Southern Methodist University. large sample assumptions implicit in traditional goodness Appears in the Proceedings of the Soth-Central SAS Users Grop Conference (SCSUG-96), Astin, TX, Oct 27-29, 1996 Fishing for Exactness Ted Pedersen Department of Compter Science & Engineering Sothern Methodist

More information

The Cryptanalysis of a New Public-Key Cryptosystem based on Modular Knapsacks

The Cryptanalysis of a New Public-Key Cryptosystem based on Modular Knapsacks The Cryptanalysis of a New Pblic-Key Cryptosystem based on Modlar Knapsacks Yeow Meng Chee Antoine Jox National Compter Systems DMI-GRECC Center for Information Technology 45 re d Ulm 73 Science Park Drive,

More information

arxiv: v1 [math.co] 25 Sep 2016

arxiv: v1 [math.co] 25 Sep 2016 arxi:1609.077891 [math.co] 25 Sep 2016 Total domination polynomial of graphs from primary sbgraphs Saeid Alikhani and Nasrin Jafari September 27, 2016 Department of Mathematics, Yazd Uniersity, 89195-741,

More information

Approach to a Proof of the Riemann Hypothesis by the Second Mean-Value Theorem of Calculus

Approach to a Proof of the Riemann Hypothesis by the Second Mean-Value Theorem of Calculus Advances in Pre Mathematics, 6, 6, 97- http://www.scirp.org/jornal/apm ISSN Online: 6-384 ISSN Print: 6-368 Approach to a Proof of the Riemann Hypothesis by the Second Mean-Vale Theorem of Calcls Alfred

More information

Characterizations of probability distributions via bivariate regression of record values

Characterizations of probability distributions via bivariate regression of record values Metrika (2008) 68:51 64 DOI 10.1007/s00184-007-0142-7 Characterizations of probability distribtions via bivariate regression of record vales George P. Yanev M. Ahsanllah M. I. Beg Received: 4 October 2006

More information

Lateral Load Capacity of Piles

Lateral Load Capacity of Piles Lateral Load Capacity of Piles M. T. DAVSSON, Department of Civil Engineering, University of llinois, Urbana Pile fondations sally find resistance to lateral loads from (a) passive soil resistance on the

More information

Multiple Choice Tries and Distributed Hash Tables

Multiple Choice Tries and Distributed Hash Tables Multiple Choice Tries and Distributed Hash Tables Luc Devroye and Gabor Lugosi and Gahyun Park and W. Szpankowski January 3, 2007 McGill University, Montreal, Canada U. Pompeu Fabra, Barcelona, Spain U.

More information

ON PREDATOR-PREY POPULATION DYNAMICS UNDER STOCHASTIC SWITCHED LIVING CONDITIONS

ON PREDATOR-PREY POPULATION DYNAMICS UNDER STOCHASTIC SWITCHED LIVING CONDITIONS ON PREDATOR-PREY POPULATION DYNAMICS UNDER STOCHASTIC SWITCHED LIVING CONDITIONS Aleksandrs Gehsbargs, Vitalijs Krjacko Riga Technical University agehsbarg@gmail.com Abstract. The dynamical system theory

More information

MEAN VALUE ESTIMATES OF z Ω(n) WHEN z 2.

MEAN VALUE ESTIMATES OF z Ω(n) WHEN z 2. MEAN VALUE ESTIMATES OF z Ωn WHEN z 2 KIM, SUNGJIN 1 Introdction Let n i m pei i be the prime factorization of n We denote Ωn by i m e i Then, for any fixed complex nmber z, we obtain a completely mltiplicative

More information

Bayes and Naïve Bayes Classifiers CS434

Bayes and Naïve Bayes Classifiers CS434 Bayes and Naïve Bayes Classifiers CS434 In this lectre 1. Review some basic probability concepts 2. Introdce a sefl probabilistic rle - Bayes rle 3. Introdce the learning algorithm based on Bayes rle (ths

More information

UNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL

UNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL 8th International DAAAM Baltic Conference "INDUSTRIAL ENGINEERING - 19-1 April 01, Tallinn, Estonia UNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL Põdra, P. & Laaneots, R. Abstract: Strength analysis is a

More information

Formules relatives aux probabilités qui dépendent de très grands nombers

Formules relatives aux probabilités qui dépendent de très grands nombers Formles relatives ax probabilités qi dépendent de très grands nombers M. Poisson Comptes rends II (836) pp. 603-63 In the most important applications of the theory of probabilities, the chances of events

More information

A generalized Alon-Boppana bound and weak Ramanujan graphs

A generalized Alon-Boppana bound and weak Ramanujan graphs A generalized Alon-Boppana bond and weak Ramanjan graphs Fan Chng Abstract A basic eigenvale bond de to Alon and Boppana holds only for reglar graphs. In this paper we give a generalized Alon-Boppana bond

More information

Collective Inference on Markov Models for Modeling Bird Migration

Collective Inference on Markov Models for Modeling Bird Migration Collective Inference on Markov Models for Modeling Bird Migration Daniel Sheldon Cornell University dsheldon@cs.cornell.ed M. A. Saleh Elmohamed Cornell University saleh@cam.cornell.ed Dexter Kozen Cornell

More information

1. INTRODUCTION. A solution for the dark matter mystery based on Euclidean relativity. Frédéric LASSIAILLE 2009 Page 1 14/05/2010. Frédéric LASSIAILLE

1. INTRODUCTION. A solution for the dark matter mystery based on Euclidean relativity. Frédéric LASSIAILLE 2009 Page 1 14/05/2010. Frédéric LASSIAILLE Frédéric LASSIAILLE 2009 Page 1 14/05/2010 Frédéric LASSIAILLE email: lmimi2003@hotmail.com http://lmi.chez-alice.fr/anglais A soltion for the dark matter mystery based on Eclidean relativity The stdy

More information

A Single Species in One Spatial Dimension

A Single Species in One Spatial Dimension Lectre 6 A Single Species in One Spatial Dimension Reading: Material similar to that in this section of the corse appears in Sections 1. and 13.5 of James D. Mrray (), Mathematical Biology I: An introction,

More information

Optimal search: a practical interpretation of information-driven sensor management

Optimal search: a practical interpretation of information-driven sensor management Optimal search: a practical interpretation of information-driven sensor management Fotios Katsilieris, Yvo Boers and Hans Driessen Thales Nederland B.V. Hengelo, the Netherlands Email: {Fotios.Katsilieris,

More information

EOQ Problem Well-Posedness: an Alternative Approach Establishing Sufficient Conditions

EOQ Problem Well-Posedness: an Alternative Approach Establishing Sufficient Conditions pplied Mathematical Sciences, Vol. 4, 2, no. 25, 23-29 EOQ Problem Well-Posedness: an lternative pproach Establishing Sfficient Conditions Giovanni Mingari Scarpello via Negroli, 6, 226 Milano Italy Daniele

More information

Evaluation of the Fiberglass-Reinforced Plastics Interfacial Behavior by using Ultrasonic Wave Propagation Method

Evaluation of the Fiberglass-Reinforced Plastics Interfacial Behavior by using Ultrasonic Wave Propagation Method 17th World Conference on Nondestrctive Testing, 5-8 Oct 008, Shanghai, China Evalation of the Fiberglass-Reinforced Plastics Interfacial Behavior by sing Ultrasonic Wave Propagation Method Jnjie CHANG

More information

Discussion of The Forward Search: Theory and Data Analysis by Anthony C. Atkinson, Marco Riani, and Andrea Ceroli

Discussion of The Forward Search: Theory and Data Analysis by Anthony C. Atkinson, Marco Riani, and Andrea Ceroli 1 Introdction Discssion of The Forward Search: Theory and Data Analysis by Anthony C. Atkinson, Marco Riani, and Andrea Ceroli Søren Johansen Department of Economics, University of Copenhagen and CREATES,

More information

Development of Second Order Plus Time Delay (SOPTD) Model from Orthonormal Basis Filter (OBF) Model

Development of Second Order Plus Time Delay (SOPTD) Model from Orthonormal Basis Filter (OBF) Model Development of Second Order Pls Time Delay (SOPTD) Model from Orthonormal Basis Filter (OBF) Model Lemma D. Tfa*, M. Ramasamy*, Sachin C. Patwardhan **, M. Shhaimi* *Chemical Engineering Department, Universiti

More information

A generalized Alon-Boppana bound and weak Ramanujan graphs

A generalized Alon-Boppana bound and weak Ramanujan graphs A generalized Alon-Boppana bond and weak Ramanjan graphs Fan Chng Department of Mathematics University of California, San Diego La Jolla, CA, U.S.A. fan@csd.ed Sbmitted: Feb 0, 206; Accepted: Jne 22, 206;

More information

Graphs and Their. Applications (6) K.M. Koh* F.M. Dong and E.G. Tay. 17 The Number of Spanning Trees

Graphs and Their. Applications (6) K.M. Koh* F.M. Dong and E.G. Tay. 17 The Number of Spanning Trees Graphs and Their Applications (6) by K.M. Koh* Department of Mathematics National University of Singapore, Singapore 1 ~ 7543 F.M. Dong and E.G. Tay Mathematics and Mathematics EdOOation National Institte

More information

A Regulator for Continuous Sedimentation in Ideal Clarifier-Thickener Units

A Regulator for Continuous Sedimentation in Ideal Clarifier-Thickener Units A Reglator for Continos Sedimentation in Ideal Clarifier-Thickener Units STEFAN DIEHL Centre for Mathematical Sciences, Lnd University, P.O. Box, SE- Lnd, Sweden e-mail: diehl@maths.lth.se) Abstract. The

More information

Decoder Error Probability of MRD Codes

Decoder Error Probability of MRD Codes Decoder Error Probability of MRD Codes Maximilien Gadolea Department of Electrical and Compter Engineering Lehigh University Bethlehem, PA 18015 USA E-mail: magc@lehighed Zhiyan Yan Department of Electrical

More information

Research Article Permanence of a Discrete Predator-Prey Systems with Beddington-DeAngelis Functional Response and Feedback Controls

Research Article Permanence of a Discrete Predator-Prey Systems with Beddington-DeAngelis Functional Response and Feedback Controls Hindawi Pblishing Corporation Discrete Dynamics in Natre and Society Volme 2008 Article ID 149267 8 pages doi:101155/2008/149267 Research Article Permanence of a Discrete Predator-Prey Systems with Beddington-DeAngelis

More information

Applying Laminar and Turbulent Flow and measuring Velocity Profile Using MATLAB

Applying Laminar and Turbulent Flow and measuring Velocity Profile Using MATLAB IOS Jornal of Mathematics (IOS-JM) e-issn: 78-578, p-issn: 319-765X. Volme 13, Isse 6 Ver. II (Nov. - Dec. 17), PP 5-59 www.iosrjornals.org Applying Laminar and Trblent Flow and measring Velocity Profile

More information

Lecture Notes: Finite Element Analysis, J.E. Akin, Rice University

Lecture Notes: Finite Element Analysis, J.E. Akin, Rice University 9. TRUSS ANALYSIS... 1 9.1 PLANAR TRUSS... 1 9. SPACE TRUSS... 11 9.3 SUMMARY... 1 9.4 EXERCISES... 15 9. Trss analysis 9.1 Planar trss: The differential eqation for the eqilibrim of an elastic bar (above)

More information

Move Blocking Strategies in Receding Horizon Control

Move Blocking Strategies in Receding Horizon Control Move Blocking Strategies in Receding Horizon Control Raphael Cagienard, Pascal Grieder, Eric C. Kerrigan and Manfred Morari Abstract In order to deal with the comptational brden of optimal control, it

More information

APPLICATION OF MICROTREMOR MEASUREMENTS TO EARTHQUAKE ENGINEERING

APPLICATION OF MICROTREMOR MEASUREMENTS TO EARTHQUAKE ENGINEERING 170 APPLICATION OF MICROTREMOR MEASUREMENTS TO EARTHQUAKE ENGINEERING I. M. Parton* and P. W. Taylor* INTRODUCTION Two previos papers pblished in this Blletin (Refs 1, 2) described the methods developed

More information

Decoder Error Probability of MRD Codes

Decoder Error Probability of MRD Codes Decoder Error Probability of MRD Codes Maximilien Gadolea Department of Electrical and Compter Engineering Lehigh University Bethlehem, PA 18015 USA E-mail: magc@lehigh.ed Zhiyan Yan Department of Electrical

More information

Stability of Model Predictive Control using Markov Chain Monte Carlo Optimisation

Stability of Model Predictive Control using Markov Chain Monte Carlo Optimisation Stability of Model Predictive Control sing Markov Chain Monte Carlo Optimisation Elilini Siva, Pal Golart, Jan Maciejowski and Nikolas Kantas Abstract We apply stochastic Lyapnov theory to perform stability

More information

PHASE STEERING AND FOCUSING BEHAVIOR OF ULTRASOUND IN CEMENTITIOUS MATERIALS

PHASE STEERING AND FOCUSING BEHAVIOR OF ULTRASOUND IN CEMENTITIOUS MATERIALS PHAS STRING AND FOCUSING BHAVIOR OF ULTRASOUND IN CMNTITIOUS MATRIALS Shi-Chang Wooh and Lawrence Azar Department of Civil and nvironmental ngineering Massachsetts Institte of Technology Cambridge, MA

More information

STURM-LIOUVILLE PROBLEMS

STURM-LIOUVILLE PROBLEMS STURM-LIOUVILLE PROBLEMS ANTON ZETTL Mathematics Department, Northern Illinois University, DeKalb, Illinois 60115. Dedicated to the memory of John Barrett. ABSTRACT. Reglar and singlar Strm-Lioville problems

More information

Chem 4501 Introduction to Thermodynamics, 3 Credits Kinetics, and Statistical Mechanics. Fall Semester Homework Problem Set Number 10 Solutions

Chem 4501 Introduction to Thermodynamics, 3 Credits Kinetics, and Statistical Mechanics. Fall Semester Homework Problem Set Number 10 Solutions Chem 4501 Introdction to Thermodynamics, 3 Credits Kinetics, and Statistical Mechanics Fall Semester 2017 Homework Problem Set Nmber 10 Soltions 1. McQarrie and Simon, 10-4. Paraphrase: Apply Eler s theorem

More information

Relativity II. The laws of physics are identical in all inertial frames of reference. equivalently

Relativity II. The laws of physics are identical in all inertial frames of reference. equivalently Relatiity II I. Henri Poincare's Relatiity Principle In the late 1800's, Henri Poincare proposed that the principle of Galilean relatiity be expanded to inclde all physical phenomena and not jst mechanics.

More information

Sareban: Evaluation of Three Common Algorithms for Structure Active Control

Sareban: Evaluation of Three Common Algorithms for Structure Active Control Engineering, Technology & Applied Science Research Vol. 7, No. 3, 2017, 1638-1646 1638 Evalation of Three Common Algorithms for Strctre Active Control Mohammad Sareban Department of Civil Engineering Shahrood

More information

Worst-case analysis of the LPT algorithm for single processor scheduling with time restrictions

Worst-case analysis of the LPT algorithm for single processor scheduling with time restrictions OR Spectrm 06 38:53 540 DOI 0.007/s009-06-043-5 REGULAR ARTICLE Worst-case analysis of the LPT algorithm for single processor schedling with time restrictions Oliver ran Fan Chng Ron Graham Received: Janary

More information

Department of Industrial Engineering Statistical Quality Control presented by Dr. Eng. Abed Schokry

Department of Industrial Engineering Statistical Quality Control presented by Dr. Eng. Abed Schokry Department of Indstrial Engineering Statistical Qality Control presented by Dr. Eng. Abed Schokry Department of Indstrial Engineering Statistical Qality Control C and U Chart presented by Dr. Eng. Abed

More information

Joint Transfer of Energy and Information in a Two-hop Relay Channel

Joint Transfer of Energy and Information in a Two-hop Relay Channel Joint Transfer of Energy and Information in a Two-hop Relay Channel Ali H. Abdollahi Bafghi, Mahtab Mirmohseni, and Mohammad Reza Aref Information Systems and Secrity Lab (ISSL Department of Electrical

More information

Network Coding for Multiple Unicasts: An Approach based on Linear Optimization

Network Coding for Multiple Unicasts: An Approach based on Linear Optimization Network Coding for Mltiple Unicasts: An Approach based on Linear Optimization Danail Traskov, Niranjan Ratnakar, Desmond S. Ln, Ralf Koetter, and Mriel Médard Abstract In this paper we consider the application

More information

Discussion Papers Department of Economics University of Copenhagen

Discussion Papers Department of Economics University of Copenhagen Discssion Papers Department of Economics University of Copenhagen No. 10-06 Discssion of The Forward Search: Theory and Data Analysis, by Anthony C. Atkinson, Marco Riani, and Andrea Ceroli Søren Johansen,

More information

The Dual of the Maximum Likelihood Method

The Dual of the Maximum Likelihood Method Department of Agricltral and Resorce Economics University of California, Davis The Dal of the Maximm Likelihood Method by Qirino Paris Working Paper No. 12-002 2012 Copyright @ 2012 by Qirino Paris All

More information

Convergence analysis of ant colony learning

Convergence analysis of ant colony learning Delft University of Technology Delft Center for Systems and Control Technical report 11-012 Convergence analysis of ant colony learning J van Ast R Babška and B De Schtter If yo want to cite this report

More information

CRITERIA FOR TOEPLITZ OPERATORS ON THE SPHERE. Jingbo Xia

CRITERIA FOR TOEPLITZ OPERATORS ON THE SPHERE. Jingbo Xia CRITERIA FOR TOEPLITZ OPERATORS ON THE SPHERE Jingbo Xia Abstract. Let H 2 (S) be the Hardy space on the nit sphere S in C n. We show that a set of inner fnctions Λ is sfficient for the prpose of determining

More information

A Survey of the Implementation of Numerical Schemes for Linear Advection Equation

A Survey of the Implementation of Numerical Schemes for Linear Advection Equation Advances in Pre Mathematics, 4, 4, 467-479 Pblished Online Agst 4 in SciRes. http://www.scirp.org/jornal/apm http://dx.doi.org/.436/apm.4.485 A Srvey of the Implementation of Nmerical Schemes for Linear

More information

Information Source Detection in the SIR Model: A Sample Path Based Approach

Information Source Detection in the SIR Model: A Sample Path Based Approach Information Sorce Detection in the SIR Model: A Sample Path Based Approach Kai Zh and Lei Ying School of Electrical, Compter and Energy Engineering Arizona State University Tempe, AZ, United States, 85287

More information

BLOOM S TAXONOMY. Following Bloom s Taxonomy to Assess Students

BLOOM S TAXONOMY. Following Bloom s Taxonomy to Assess Students BLOOM S TAXONOMY Topic Following Bloom s Taonomy to Assess Stdents Smmary A handot for stdents to eplain Bloom s taonomy that is sed for item writing and test constrction to test stdents to see if they

More information

On capital allocation by minimizing multivariate risk indicators

On capital allocation by minimizing multivariate risk indicators On capital allocation by minimizing mltivariate risk indicators Véroniqe Mame-Deschamps, Didier Rllière, Khalil Said To cite this version: Véroniqe Mame-Deschamps, Didier Rllière, Khalil Said. On capital

More information

Gradient Projection Anti-windup Scheme on Constrained Planar LTI Systems. Justin Teo and Jonathan P. How

Gradient Projection Anti-windup Scheme on Constrained Planar LTI Systems. Justin Teo and Jonathan P. How 1 Gradient Projection Anti-windp Scheme on Constrained Planar LTI Systems Jstin Teo and Jonathan P. How Technical Report ACL1 1 Aerospace Controls Laboratory Department of Aeronatics and Astronatics Massachsetts

More information

Pulses on a Struck String

Pulses on a Struck String 8.03 at ESG Spplemental Notes Plses on a Strck String These notes investigate specific eamples of transverse motion on a stretched string in cases where the string is at some time ndisplaced, bt with a

More information

L 1 -smoothing for the Ornstein-Uhlenbeck semigroup

L 1 -smoothing for the Ornstein-Uhlenbeck semigroup L -smoothing for the Ornstein-Uhlenbeck semigrop K. Ball, F. Barthe, W. Bednorz, K. Oleszkiewicz and P. Wolff September, 00 Abstract Given a probability density, we estimate the rate of decay of the measre

More information

Simplified Identification Scheme for Structures on a Flexible Base

Simplified Identification Scheme for Structures on a Flexible Base Simplified Identification Scheme for Strctres on a Flexible Base L.M. Star California State University, Long Beach G. Mylonais University of Patras, Greece J.P. Stewart University of California, Los Angeles

More information

Reaction-Diusion Systems with. 1-Homogeneous Non-linearity. Matthias Buger. Mathematisches Institut der Justus-Liebig-Universitat Gieen

Reaction-Diusion Systems with. 1-Homogeneous Non-linearity. Matthias Buger. Mathematisches Institut der Justus-Liebig-Universitat Gieen Reaction-Dision Systems ith 1-Homogeneos Non-linearity Matthias Bger Mathematisches Institt der Jsts-Liebig-Uniersitat Gieen Arndtstrae 2, D-35392 Gieen, Germany Abstract We describe the dynamics of a

More information

3.1 The Basic Two-Level Model - The Formulas

3.1 The Basic Two-Level Model - The Formulas CHAPTER 3 3 THE BASIC MULTILEVEL MODEL AND EXTENSIONS In the previos Chapter we introdced a nmber of models and we cleared ot the advantages of Mltilevel Models in the analysis of hierarchically nested

More information

Cuckoo hashing: Further analysis

Cuckoo hashing: Further analysis Information Processing Letters 86 (2003) 215 219 www.elsevier.com/locate/ipl Cckoo hashing: Frther analysis Lc Devroye,PatMorin School of Compter Science, McGill University, 3480 University Street, Montreal,

More information

(r2 + rf + k2r2f) that occurs when the space form of tne wave equation is expressed. Fl(U)--" [Ju(kr)u(ka)- Ju(ka)u(kr)]f(r)drr (1.

(r2 + rf + k2r2f) that occurs when the space form of tne wave equation is expressed. Fl(U)-- [Ju(kr)u(ka)- Ju(ka)u(kr)]f(r)drr (1. Internat. J. Math. & Math. Sci. VOL. ii NO. 4 (1988) 635-642 635 ON AN INTEGRAL TRANSFORM D. NAYLOR niversity of Weste Ontario London, Ontario, Canada (Received March 28, 1988) ABSTRACT. A formla of inversion

More information

arxiv: v1 [math.dg] 22 Aug 2011

arxiv: v1 [math.dg] 22 Aug 2011 LOCAL VOLUME COMPARISON FOR KÄHLER MANIFOLDS arxiv:1108.4231v1 [math.dg] 22 Ag 2011 GANG LIU Abstract. On Kähler manifolds with Ricci crvatre lower bond, assming the real analyticity of the metric, we

More information

On the scaling ranges of detrended fluctuation analysis for long-memory correlated short series of data

On the scaling ranges of detrended fluctuation analysis for long-memory correlated short series of data On the scaling ranges of detrended flctation analysis for long-memory correlated short series of data arxiv:126.17v1 [physics.data-an] 5 Jn 212 Darisz Grech (1) and Zygmnt Mazr (2) (1) Institte of Theoretical

More information

This is a repository copy of Enhanced Bounded Integral Control of Input-to-State Stable Nonlinear Systems.

This is a repository copy of Enhanced Bounded Integral Control of Input-to-State Stable Nonlinear Systems. This is a repository copy of Enhanced Bonded Integral Control of Inpt-to-State Stable Nonlinear Systems. White Rose Research Online URL for this paper: http://eprints.hiterose.ac.k/4854/ Version: Accepted

More information