Simulated Annealing with Time dependent. Matthias Lowe. Universitat Bielefeld. Fakultat fur Mathematik. Postfach Bielefeld.
|
|
- Neal Greer
- 5 years ago
- Views:
Transcription
1 Simulated Annealing with Time dependent Energy Function via Sobolev inequalities Matthias Lowe Universitat Bielefeld Fakultat fur Mathematik Postfach Bielefeld Abstract We analyze the Simulated Annealing Algorithm with an energy function U t that depends on time. Assuming some regularity conditions on U t (especially that U t does not change too quickly in time), and choosing a logarithmic cooling schedule for the algorithm, we derive bounds on the Radon-Nikodym density of the distribution of the annealing algorithm at time t with respect to the invariant measure at time t. Moreover we estimate the entrance time of the algorithm into typical subsets V of the state space in terms of t (V ). c Keywords: Simulated Annealing, Sobolev inequalities, Spectral gap, Markov processes Introduction Let X be a nite set. The well known Simulated Annealing (SA) algorithm is an inhomogeneous Markov process Y t on X with the aim to minimize a given function U : X! R. The idea behind SA is to think of U as an energy function and to choose the Markov process in such a way that the transition kernel at time t has at its invariant measure t, the Gibbs distribution with
2 energy U at temperature = t. So, if t diverges to innity in the limit for large times the invariant measure becomes concentrated on the set of minima of U. On the other hand one has to "freeze" the process slowly enough to guarantee convergence of its distribution at time t towards equilibrium. Since the introduction of SA by Kirkpatrick et al. in [K83] it has found a wide range of applications (see [GG84],[V93]), and the problems mentioned above, i.e. the optimal choice of t and the asymptotic behavior of SA have been studied by several authors. The problem of an asymptotically optimal "cooling schedule" (i.e. the dependence of t from t) has been settled by Hajek [Ha88], who showed the optimality of a logarithmic freezing schedule t = log(a + Bt), where A and B are positive constants and m is, roughly m speaking, the largest hill to be climbed to reach the global minimum. The rate of convergence has been investigated by Gidas in [Gi85] and Mitra et al. in [M85] on the basis of Dobrushin coecients for Markov chains. Cantoni has obtained bounds of the type P (Y t E) O(t? ) in [C90] by using Freidlin-Vencel techniques (see [FV84]) on perturbed dynamical systems. In [HS88] Holley and Stroock were able to give bounds for P x (Y t = y)= t (y) for reversible SA processes based on Sobolev inequalities. Additionally they obtained sharp estimates for the spectral gap of the operators associated to the SA process. Based on rened spectral gap estimates Gotze in [G9] proved eective bounds for the waiting time of escaping from an unlikely set (i.e. a set of high energy). These techniques were generalized by Deuschel and Mazza in [DM9] to Markov processes on nite sets similar to SA. In this paper we are going to analyze that variant of the SA algorithm, where U depends on t, i.e. U = U t. These processes have rst been studied by Frigerio and Grillo in [FG93]. Not only they are a natural generalization of SA, but also they can be regarded as a discrete version of quantum diusions (see [FG93], Section 4). Moreover these processes become interesting when thinking of SA as a sort of genetic algorithm (the reader is referred to the book of Holland [H75] for an introduction and a survey on these algorithms). Then the SA process can be regarded as an individual undergoing a mutation, the mutation being governed by the "tness-function" U. In such a model the time dependence of U comes up naturally, if one e.g. thinks of the tness of an individual being given by its environment, that usually changes in time.
3 In [FG93] Frigerio and Grillo have proved that for SA with time dependent energy function the distribution of the process and its invariant measure become indistinguishable in the weak topology for large times t, if the following catalogue of conditions is fullled: the SA process is driven by a bistochastic matrix, U t (x) is uniformly bounded in x and t by a nonnegative number M and has a derivative with respect to t bounded by C t for some C; > 0, the generator L t associated to the process at time t admits a bound on its spectral gap of the form e?tm for some ; m > 0, and we choose a logarithmic cooling schedule t = log( + t) with h > m and m=h <. These h conditions, that are especially satised, when U t admits an uniform limit U look rather natural from a point of view of Sobolev inequalities and spectral gap estimates. On one hand one has chosen a cooling schedule that also is close to optimal, when U is constant in t. On the other hand in some sense m and M are the only values which enter in these techniques for time independent energy function. So if the variation of U t is moderate one should expect spectral gap estimates to work. The aim of this paper is to show that under the conditions formulated by Frigerio and Grillo in [FG93] this is indeed the case. By Sobolev inequalities we obtain bounds for the density of the distribution of the process with respect to the invariant distribution at time t and are able to analyze the speed of convergence towards equilibrium. It turns out the conditions formulated above cannot essentially be weakened. Our main tools for this approach are those of Holley and Stroock in [HS88] and Gotze in [G9]. We organize the rest of the paper in the following way. Section contains our basic denitions and a formal description of the model. In Section 3 we describe the behavior of the Radon-Nikodym derivative of the distribution of the process at time t with respect to the Gibbs-distribution at that time with the help of Sobolev inequalities and check, under which conditions these are fullled. Moreover we use spectral gap estimates to bound the probability of large entrance times into typical sets. The Model Let us now give a formal description of the algorithm we have mind. To this end let q(x; ) be a probability distribution on X. Assume that the 3
4 random walk associated with q(x; y) is irreducible (i.e. for all pairs x; y X there is a number n := n(x; y) such that q n (x; y) > 0 where we have put q n (x; y) = P zx q(x; z)q n? (z; y)). Moreover for simplicity we assume that q(x; y) is time-reversible with respect to a probability measure on X, i.e. (x; y) := (x)q(x; y) = (y)q(y; x) = (y; x) () for all x; y X. Note that is q-stationary and therefore irreducibility of q implies that charges every point in X. This assumption of reversibility often makes life easier in convergence theorems for stochastic processes. On the other hand a comparison of the spectral gap estimates for time homogeneous, reversible Markov chains by Diaconis and Stroock [DS9] with those for homogeneous, irreversible chains by Fill [F9] shows that there is some hope to drop the reversibility condition in future work. Let (U t ) t be a sequence (indexed by R + ) of functions U t : X! R. Moreover let us assume that t 7! t is a monotonely increasing function in t, which is strictly positive ( t will play the role of the inverse temperature at time t). Corresponding to q, U t, and t we may dene a transition probability q t (x; y) by q t (x; y) = ( exp(?t (U t (y)? U t (x)) + )q(x; y) if x 6= y? P z6=y q t (x; z) if x = y () where we have put (x) + := max(x; 0). It is easy to check that q t is reversible with respect to the Gibbs measure t (with reference measure ) with energy U t at inverse temperature t dened by t (y) := e?tut(y) Z t (y) (3) where Z t := P y e?tut(y) (y) is the partition function. Note also that q t again is irreducible and reversible and that t is its stationary distribution, or, equivalently, that the operator Q t given by [Q t ](x) := X yx (y)q t (x; y) (4) is a self-adjoint contraction on L ( t ). 4
5 If we now dene the operator L t := Q t? I on L ( t ), i.e. [L t ](x) := X yx((y)? (x))q t (x; y); x X (5) we are able to describe the continuous time SA algorithm. Namely let its transition probabilities P s;t (x; y), i.e the probabilities to arrive in y at time t when starting in x at time s, for s t be given by the [P s;t](x) = [P s;t [L t ]](x) for s < t (6) and [P s;s ](x) := (x) for : X! R. Here [P s;t ](x) := X yx (y)p s;t (x; y); x X: The intuitive meaning of.6 is that starting in x at time s the process waits for exponential time with mean one and then chooses a neighbor y according to q(x; y). If then U t (y) U t (x) (where t is the instantaneous time after having waited) the algorithm moves to y. Otherwise it moves to y with probability exp(? t (U t (y)? U t (x))) and stays in x with probability? exp(? t (U t (y)? U t (x))). Obviously this is the continuous time analogue of the ordinary discrete time SA process where we have chosen U to depend on the time. The advantage of using continuous time instead of discrete time becomes especially apparent when studying the time evolution of f t, the Radon-Nikodym derivative of the distribution of Y t with respect to t, which then can be described by means of dierential equations instead of dierence equations that are harder to solve. Let us remark that together with Q t also L t is a self-adjoint operator on L ( t ) with largest eigenvalue 0. Due to the irreducibility of q t the eigenspace of the eigenvalue 0 is one-dimensional and the constants are the only eigenfunctions. The key observation on which most of the following techniques are based is the relation of second smallest eigenvalue t of?l t and the Dirichlet form associated with?l t. This relation is given by t = inff E t(; ) V ar t () : : X! R non constantg: (7) 5
6 Here E t (; ) :=? < ; [L t ] > t = X x;yx ((x)? (y)) t (x; y) (8) is the discrete Dirichlet form associated to L t, < ; > t denotes the scalar product in L ( t ), V ar t () is the variance of with respect to t and t is given by t (x; y) := t (x)q t (x; y) = t (y)q t (y; x) = t (y; x): (9) Note that E t is related to the Dirichlet form E 0 (; ) := P x;yx ((x)? (y)) (x; y) via the relation E t (; ) = X x;yx ((x)? (y)) (x; y) exp(? t max(u t (x); U t (y))) Z t : (0) (.7) is derived by the well known variational characterization of eigenvalues (see e.g. [HJ85]). For the rest of the paper let us assume that the (U t ) t fulll the following conditions : First we may without loss of generality assume that min x U t (x) = 0. Moreover let us choose such sequences of energy functions that there exist universal constants 0 and ; C; m; M > 0 such that 0 U t (x) M for all x X, for all t 0 ; U t(x) C ( + t) for all x X; () t e?tm for all t 0: (3) We will conclude this section with the observation that assumption (.3) is automatically fullled if we choose m similar to [Ha88] as the maximal hill to be climbed at any time to reach to global minimum. More formally for xed t and any two points x; y X let us call a sequence (x i ) ik a path from x to y, if x = x, x k = y and q t (x i? ; x i ) > 0 for all i = ; : : : ; k?. Dene m t (x) := minfe 0 : there exists y X with U(y) < U(x) and a path (x i ) ik from x to y with U(x i ) U(x) + E for all i kg: 6
7 Finally put m t := max xx m t(x): (4) With these denitions we are able to give the precise order (up to a factor) of the spectral gap t : Theorem. There exist constants ;? > 0 such that e?tmt t?e?tmt (5) Proof: Note that the statement of Theorem. is the analogue to the spectral gap estimate Theorem. in [HS88] for the ordinary SA. As a matter of fact, since we analyze the spectral gap of the instantaneous time operator L t, which does not "feel" the change of U t in t, the proof stays the same for our case. Since we are not going to make use of the upper bound in the rest of the paper, we will skip its proof here. The proof of the lower bound follows the lines of the proof of Lemma.7 in [HS88]. By (.7) it suces to show the existence of a constant, such that for all t, and for all f : X! R E t (f; f) V ar t (f) =? P xx f(x)[l t f](x) t (x) V ar t (f) This is done by bounding the variance as follows : For any two x; y X take a path p x;y t with vertices p x;y t that where the % x;y t length of p x;y t max p x;y t (i) U t (p x;y t (i)) = min % x;y t max % x;y t (j) e?tmt : (6) U t (% x;y t (j)) are paths form x to y with vertices % x;y t and put N t := max x;yx n t(x; y): (i) from x to y such (j). Let n t (x; y) be the We introduce the indicator ( if there exists 0 i t nt (x; y) : p x;y t (i) = z and p x;y t (i + ) = w z;w(x; y) = 0 otherwise. In the following we put t z;w (x;y) X (x;y) V ar t (f) = x;yx = 0 if t z;w (x; y) = 0. Then for any f : X! R (f(x)? f(y)) t (x) t (y) 7
8 = X x;yx X x;yx 0 nt(x;y) i= n t(x;y) n t (x; y) f(p x;y t (i))? f(p x;y t (i? )) X i= (f(p x;y t A t (x) t (y) (i))? f(p x;y t (i? ))) t (x) t (y) by Jensen's inequality. Therefore X V ar t (f) N t t z;w(x; y)(f(z)? f(w)) e?t max(ut(z);ut(w)) (z; w) x;y;z;wx t (x) t (y) e?t max(ut(z);ut(w)) (z; w) X N t [ max t t (x) t (y)z t z;w(x; y) e?t max(ut(z);ut(w)) (z; w) ] z;wx x;yx [ X z;wx (f(z)? f(w)) e?t max(ut(z);ut(w)) (z; w) Z t ]: (7) Since N t is bounded by the size of the vertex set jxj (because of irreducibility) and the second term in brackets in (.7) is twice the Dirichlet form E t (f; f) (due to (.0)) it suces to calculate that t t (x) t (y)z t z;w(x; y) e?t max(ut(z);ut(w)) (z; w) = t z;w(x; y) (z; w) e tmt [ t z;w(x; y) (z; w) (x)(y) e?t(max(ut(z);ut(w))?ut(x)?ut(y)) Z t (x)(y) min vx (v) ]; since by denition m t (max(u t (z); U t (w))?u t (x)?u t (y)) whenever t z;w (x; y) = for all x; y; z; w X. Since [ t z;w (x;y) (x)(y) ] is bounded independently (z;w) minvx (v) of t, this proves the lower bound. Remark. Note that under condition (.) especially m := max t m t exists and therefore we obtain for all t. t e?tm 8
9 3 Convergence of the Algorithm via Sobolev inequalities In this section we are going to analyze the behavior of the density of the distribution of the process with respect to the invariant measure on the basis of a Sobolev inequality, which we rst assume and then show to hold true. The results of this section are similar to those for ordinary SA proved by Holley and Stroock in [HS88]. Like they we start with a result on a dierential inequality, on which the rest of the section is based. Lemma 3. Let a; b; c be non-negative numbers, let a > 0, let " > 0, and 0 <. Let u be a dierentiable function u : R +! R + 0 fullling d dt u(t) =: u0 (t) ( + t ) [?au +" + bu + cu ] for t > 0 (where we adopt the convention that primes denote derivatives with respect to t). Then u(t) where := max(b; a + c). 4 4" a? exp( " (? ( +? t)? )) Remark 3. Lemma 3. generalizes Lemma. in [HS88] where =. Our proof is similar. 3 5 " Proof: Put v(t) := (u(t)). Then v 0 (t) = u0 (t) (u(t)) ( + t ) [?a( (v(t) )+" ) + b v(t) + c] = ( + t ) [?a( + ( v(t) )+" + b v(t) + a + c] ( + t ) [?a v(t) ( + 4 " )+" + ( + v(t) )]: 9
10 Set w(t) := ( + v(t) ) exp(? (+t)? ). Then? w 0 (t) = exp(? ( + t)?? exp(? ( + t)?? =?a + t)? exp("( 4 "? Dene z(t) := (w(t))?" and obtain and therefore Hence )[ v0 (t)? )[? a v(t) ( + 4 " z 0 (t) =?"w0 (t) (w(t)) +" v(t) ( + ( + t) )] )+" ]( + t ) )(w(t)) +" ( + t ) : z 0 (t) a" + t)? exp("( )( 4 "? + t ) : z(t) a" Z t ( + s)? exp(" )( 4 " 0? + s ) ds?? = a" Z (+t) 4 " = a" 4? exp("s)ds + t)? [exp("( "? Recalling the denitions of z; w and v this yields exp(" which implies ( + t)?? u(t) )(u(t))?" a" 4 4 4" a )? exp("? )]: + t)? [exp("( "?? exp( " (? ( +? t)? )) )? exp("? )]; To make use of Lemma 3. let us assume that the following hypothesis on the SA are fullled: There is a 0 and " :
11 (A) there is a constant B < and such that where R t := sup x U 0 (x). 0 tm + t R t B ( + t) (A) there is a constant A and a < p < such that for all t 0 jjf? < f > t jj p;t A( + t) E t (f; f) for all f L p ( t ), where < > t denotes expectation with respect to t. The following Lemma shows, how to use Lemma 3. and assumption (A) as well as the Sobolev inequality (A) to bound the density of the distribution of the process at time t with respect to the invariant measure at time t. To this end set f t (y) := P x(y t = y) t (y) where x is the starting point of the SA process. Lemma 3.3 Assume that hypotheses (A) and (A) are fullled and set " := p? =p as well as K := 4 " (4AB + )? exp(? "(B + =A)) (? e))! :? Then after a waiting time of t e??. Proof: Therefore Note that jjf t jj ;t + K " d dt Z t =?Z t ( 0 t < U t > t + t < U 0 t > t): d dt t(y) =? t (y)( 0 t ^U t (y) + t ^U 0 t (y))
12 where ^Ut (y) := U t (y)? < U t > t. Using the Fokker-Planck equations (.6) we arrive at d dt f t(y) = X (f t (z)? f t (y)) t (z; y) + f t (y)( 0 t t (y) ^U 0 t (y) + t ^U t (y)): zx Putting u(t) := jjf t? jj ;t = jjf tjj ;t?,this implies u 0 (t) =?E t (f t ; f t ) + X y X f t (y)( 0 t ^U t (y) + t ^U 0 t (y)) t (y) =?E t (f t ; f t ) + X y X(f t (y)? ) ( 0 t ^U t (y) + t ^U 0 t (y)) t (y) + X y X(f t (y)? )( 0 t ^U t (y) + t ^U 0 t (y)) t (y)?e t (f t ; f t ) + X y X(f t (y)? ) ( 0 t M + t R t ) t (y) + X y X(f t (y)? )( 0 t M + t R t ) t (y) ( + t ) [? A u+" + Bu + Bu ]: Thus the estimate of the lemma follows from Lemma 3.. It remains to show that under conditions (.)-(.3) we have modeled the SA in such a way that assumptions (A) and (A) are fullled. This is the content of the following lemma. Lemma 3.4 Under the assumptions (.)-(.3) for any p > there exists a cooling schedule t such that the associated SA algorithm satises (A) and (A). Proof: We follow the ideas of the proof of Theorem. in [HS88]. So we choose < and according to assumption (.) and we put t = log( + t); m + " where " = p? M. Now for any xed t 0 on one hand we have p 0 t M + tr t M (m + ")( + t) const: ( + t) + log( + t) m + " C ( + t)
13 with a constant depending on for any 0 < <. On the other hand for every f : X! R X X < f > t ) yx(f(y)? p t (y) = (f(y)? < f > t ) j(f(y)? < f > t )j p? t (y) yx X (f(y)? < f > t ) t (y)jj(f? < f > t )jj p? yx emt E t (f; f)jjf? < f > t jj p? Since Z t we obtain Hence jjf? < f > t jj p e?tm min xx jjf? < f > t jj p min t(x) xx jjf? < f > t jj p p;t: jjf? < f > t jj p p;t M (p?) emt t jjf? < f > t jj p? p;t E t (f; f)e p : (min xx (x)) p? p Thus for the nite constant A = jjf? < f > t jj p;t (minx (x)) p?=p M (p?) Aet(m+ p E t (f; f): So the choice of t and " implies the statement of the lemma. Summarizing we have proved the following theorem Theorem 3.5 Let Y t be the SA process with time dependent energy function U t. Under conditions (.)-(.3) the cooling schedule t = log( + t) m+" with 0 < < implies the following bound on the Radon-Nikodym density f t (y) := Px(Yt=y) t(y) jjf t (y)jj + K " for t e??, and any starting point x X. Here the constants " and K are chosen according to lemmas above. 3
14 Remarks 3.6. Following the same ideas as in [HS88], Section, it is possible to derive uniform bounds for f t as well. We omit a detailed discussion here.. The conditions of Theorem 3.5 are not only sucient to bound f t in L ( t ) (and therefore to ensure convergence of the algorithm), but also are nearly necessary conditions to make SA sensible. Especially if m t and therefore any bound on U t diverges to innity with t! even the cooling schedule t = log( + t), which we expect to be optimal, is in fact a "heating" schedule, i.e. t! 0 and therefore the m equilibrium measures at time t converge to the uniform distribution on X. Obviously, in this situation SA does not make sense any longer. Also note that M being nite does not imply that Mt 0 = 0 (where M t = sup x U t (x)), but only that, roughly speaking, Mt 0 t? for some >, i.e M t changes noticeably slower than U t does. 3. Note that if in condition (.) > U t admits a limit (for t! ) U. As one would expect, SA with time dependent energy function U t then asymptotically behaves like SA with energy function U. Let us nally show, how to use Theorem 3.5 to derive a bound on the entrance time of the algorithm into typical sets. Note that the following theorem especially becomes interesting, when the minima of U t become concentrated on a subset V X. Theorem 3.7 Let V X and denote by V c the rst entrance time of the algorithm into V c when starting in x V. Then under the conditions and with the notations of Theorem 3.5 for t e??. P ( V c > t) ( + K " ) t (V ) + t (V ) Proof: We have P ( V c > t) = P (Y s V 8s t) inf st P (Y s V ); 4
15 because of fy r V 8r tg fy s V g if s t. Moreover (P (Y t V )? t (V )) = X t? ) V (y) t (y) A yx(f jjf t? jj ;t t (V ): Hence the assertion follows form Theorem 3.5. Surprisingly the SA with time dependent energy function (and < ) seems to resist calculations as in [G9] for the ordinary SA, where more applicable constants than our ( + K ") were obtained. On the other hand a rough estimate of K yields the following corollary. Corollary 3.8 With the notations of Theorem 3.7 we have 0 P ( V c > 56C + ( A e q (min x (x)) )=4 t (V ) + t (V ) for t e?=?. Proof: Estimating the constants in K with p = 4, i.e " = and = = yields the bound. Remark 3.9 Throughout the paper we have used a bound for the spectral gap of the form e?tm. As Theorem. shows these bounds are optimal up to a constant. On the other hand, since usually m, or m t resp. are as hard to determine as the argument where U t attains its minimum, these bounds as well as the associated cooling schedules are not very applicable. For practicle use Poincar-type or Cheeger-type bounds as have been obtained in [DS9] and [G9] seem to be preferable. References [C90] O. Cantoni, Asymptotics of simulated annealing algorithms, Ph.D. Thesis, Universit Paris-Sud, Orsey (990) 5
16 [DM9] J.D. Deuschel, C. Mazza, L Convergence of Time Non- Homogeneous Markov Processes I: Spectral Estimates, Preprint (99) [DS9] P. Diaconis, D. Stroock, Geometric Bounds for Eigenvalues of Markov chains, Ann. Appl. Prob., 36-6 (99) [F9] J.A. Fill, Eigenvalue bounds on Convergence to stationarity for nonreversible Markov chains, with an application to exclusion processes, Ann. Appl. Prob., 6-87 (99) [FV84] M.I. Freidlin, A.D. Vencel, Random Pertubations of Dynamical Systems, Springer, Grundlehren d. math. Wissenschaften 60 (984) [FG93] A. Frigerio, G. Grillo, Simulated annealing with time-dependent energy function, Math. Zeit. 3, 97-6 (993) [GG84] S. Geman, D. Geman, Stochastic Relaxation, Gibbs Distribution, and the Bayesian Restoration of images, IEEE Trans. on Pattern Anal. and Machine Intelligence, PAMI-6, 7-74 (984) [Gi85] [G9] B. Gidas Nonstationary chains and the convergence of the annealing algorithm, J. Stat. Phys. 39, 73-3 (985) F. Gotze, Rate of convergence for simulated annealing processes, Preprint 9-008, SFB 343, Bielefeld (99) [Ha88] B. Hajek, Cooling Schedules for optimal Annealing, Math. of Operations Research 3, 3-39 (988) [H75] J. Holland, Adaption in Natural and Articial Systems, Ann Arbor, University of Michigan Press (975) [HS88] R. Holley, D. Stroock, Simulated Annealing via Sobolev inequalities, Comm. Math. Phys. 5, (988) [HJ85] R. Horn, C. Johnson, Matrix analysis, Camb. Univ. Press (985) [K83] S. Kirkpatrick, C. Gelatt, M.P. Vecchi, Optimization by Simulated Annealing, Science 0, (983) 6
17 [M85] [V93] D. Mitra, F. Romeo, A. Sangiovanni-Vincentelli, Convergence and nite time behavior of simulated annealing, Proc. 4th Conf. Dec. a. Control, (985) R.V.V. Vidal (Ed.), Applied Simulated Annealing, Lecture Notes in Economics and Math. Syst. 369, Springer, Berlin (993) 7
OPTIMIZATION BY SIMULATED ANNEALING: A NECESSARY AND SUFFICIENT CONDITION FOR CONVERGENCE. Bruce Hajek* University of Illinois at Champaign-Urbana
OPTIMIZATION BY SIMULATED ANNEALING: A NECESSARY AND SUFFICIENT CONDITION FOR CONVERGENCE Bruce Hajek* University of Illinois at Champaign-Urbana A Monte Carlo optimization technique called "simulated
More informationOn the convergence rates of genetic algorithms
Theoretical Computer Science 229 (1999) 23 39 www.elsevier.com/locate/tcs On the convergence rates of genetic algorithms Jun He a;, Lishan Kang b a Department of Computer Science, Northern Jiaotong University,
More information1 Introduction This work follows a paper by P. Shields [1] concerned with a problem of a relation between the entropy rate of a nite-valued stationary
Prexes and the Entropy Rate for Long-Range Sources Ioannis Kontoyiannis Information Systems Laboratory, Electrical Engineering, Stanford University. Yurii M. Suhov Statistical Laboratory, Pure Math. &
More informationThe simple slice sampler is a specialised type of MCMC auxiliary variable method (Swendsen and Wang, 1987; Edwards and Sokal, 1988; Besag and Green, 1
Recent progress on computable bounds and the simple slice sampler by Gareth O. Roberts* and Jerey S. Rosenthal** (May, 1999.) This paper discusses general quantitative bounds on the convergence rates of
More informationON THE REGULARITY OF SAMPLE PATHS OF SUB-ELLIPTIC DIFFUSIONS ON MANIFOLDS
Bendikov, A. and Saloff-Coste, L. Osaka J. Math. 4 (5), 677 7 ON THE REGULARITY OF SAMPLE PATHS OF SUB-ELLIPTIC DIFFUSIONS ON MANIFOLDS ALEXANDER BENDIKOV and LAURENT SALOFF-COSTE (Received March 4, 4)
More informationPointwise convergence rate for nonlinear conservation. Eitan Tadmor and Tao Tang
Pointwise convergence rate for nonlinear conservation laws Eitan Tadmor and Tao Tang Abstract. We introduce a new method to obtain pointwise error estimates for vanishing viscosity and nite dierence approximations
More informationAnalysis on Graphs. Alexander Grigoryan Lecture Notes. University of Bielefeld, WS 2011/12
Analysis on Graphs Alexander Grigoryan Lecture Notes University of Bielefeld, WS 0/ Contents The Laplace operator on graphs 5. The notion of a graph............................. 5. Cayley graphs..................................
More informationCongurations of periodic orbits for equations with delayed positive feedback
Congurations of periodic orbits for equations with delayed positive feedback Dedicated to Professor Tibor Krisztin on the occasion of his 60th birthday Gabriella Vas 1 MTA-SZTE Analysis and Stochastics
More informationMax-Planck-Institut fur Mathematik in den Naturwissenschaften Leipzig Uniformly distributed measures in Euclidean spaces by Bernd Kirchheim and David Preiss Preprint-Nr.: 37 1998 Uniformly Distributed
More informationMonte Carlo methods for sampling-based Stochastic Optimization
Monte Carlo methods for sampling-based Stochastic Optimization Gersende FORT LTCI CNRS & Telecom ParisTech Paris, France Joint works with B. Jourdain, T. Lelièvre, G. Stoltz from ENPC and E. Kuhn from
More informationN.G.Bean, D.A.Green and P.G.Taylor. University of Adelaide. Adelaide. Abstract. process of an MMPP/M/1 queue is not a MAP unless the queue is a
WHEN IS A MAP POISSON N.G.Bean, D.A.Green and P.G.Taylor Department of Applied Mathematics University of Adelaide Adelaide 55 Abstract In a recent paper, Olivier and Walrand (994) claimed that the departure
More informationIntroduction to Machine Learning CMU-10701
Introduction to Machine Learning CMU-10701 Markov Chain Monte Carlo Methods Barnabás Póczos & Aarti Singh Contents Markov Chain Monte Carlo Methods Goal & Motivation Sampling Rejection Importance Markov
More informationLogarithmic Harnack inequalities
Logarithmic Harnack inequalities F. R. K. Chung University of Pennsylvania Philadelphia, Pennsylvania 19104 S.-T. Yau Harvard University Cambridge, assachusetts 02138 1 Introduction We consider the relationship
More informationQUASI-UNIFORMLY POSITIVE OPERATORS IN KREIN SPACE. Denitizable operators in Krein spaces have spectral properties similar to those
QUASI-UNIFORMLY POSITIVE OPERATORS IN KREIN SPACE BRANKO CURGUS and BRANKO NAJMAN Denitizable operators in Krein spaces have spectral properties similar to those of selfadjoint operators in Hilbert spaces.
More information1 Introduction We will consider traveling waves for reaction-diusion equations (R-D) u t = nx i;j=1 (a ij (x)u xi ) xj + f(u) uj t=0 = u 0 (x) (1.1) w
Reaction-Diusion Fronts in Periodically Layered Media George Papanicolaou and Xue Xin Courant Institute of Mathematical Sciences 251 Mercer Street, New York, N.Y. 10012 Abstract We compute the eective
More information215 Problem 1. (a) Define the total variation distance µ ν tv for probability distributions µ, ν on a finite set S. Show that
15 Problem 1. (a) Define the total variation distance µ ν tv for probability distributions µ, ν on a finite set S. Show that µ ν tv = (1/) x S µ(x) ν(x) = x S(µ(x) ν(x)) + where a + = max(a, 0). Show that
More information2 Section 2 However, in order to apply the above idea, we will need to allow non standard intervals ('; ) in the proof. More precisely, ' and may gene
Introduction 1 A dierential intermediate value theorem by Joris van der Hoeven D pt. de Math matiques (B t. 425) Universit Paris-Sud 91405 Orsay Cedex France June 2000 Abstract Let T be the eld of grid-based
More informationr( = f 2 L 2 (1.2) iku)! 0 as r!1: (1.3) It was shown in book [7] that if f is assumed to be the restriction of a function in C
Inverse Obstacle Problem: Local Uniqueness for Rougher Obstacles and the Identication of A Ball Changmei Liu Department of Mathematics University of North Carolina Chapel Hill, NC 27599 December, 1995
More informationON STATISTICAL INFERENCE UNDER ASYMMETRIC LOSS. Abstract. We introduce a wide class of asymmetric loss functions and show how to obtain
ON STATISTICAL INFERENCE UNDER ASYMMETRIC LOSS FUNCTIONS Michael Baron Received: Abstract We introduce a wide class of asymmetric loss functions and show how to obtain asymmetric-type optimal decision
More information6 Markov Chain Monte Carlo (MCMC)
6 Markov Chain Monte Carlo (MCMC) The underlying idea in MCMC is to replace the iid samples of basic MC methods, with dependent samples from an ergodic Markov chain, whose limiting (stationary) distribution
More informationConvex Optimization CMU-10725
Convex Optimization CMU-10725 Simulated Annealing Barnabás Póczos & Ryan Tibshirani Andrey Markov Markov Chains 2 Markov Chains Markov chain: Homogen Markov chain: 3 Markov Chains Assume that the state
More information2 The second case, in which Problem (P 1 ) reduces to the \one-phase" problem (P 2 ) 8 >< >: u t = u xx + uu x t > 0, x < (t) ; u((t); t) = q t > 0 ;
1 ON A FREE BOUNDARY PROBLEM ARISING IN DETONATION THEORY: CONVERGENCE TO TRAVELLING WAVES 1. INTRODUCTION. by M.Bertsch Dipartimento di Matematica Universita di Torino Via Principe Amedeo 8 10123 Torino,
More informationA converse Gaussian Poincare-type inequality for convex functions
Statistics & Proaility Letters 44 999 28 290 www.elsevier.nl/locate/stapro A converse Gaussian Poincare-type inequality for convex functions S.G. Bokov a;, C. Houdre ; ;2 a Department of Mathematics, Syktyvkar
More informationA Proof of the EOQ Formula Using Quasi-Variational. Inequalities. March 19, Abstract
A Proof of the EOQ Formula Using Quasi-Variational Inequalities Dir Beyer y and Suresh P. Sethi z March, 8 Abstract In this paper, we use quasi-variational inequalities to provide a rigorous proof of the
More informationgrowth rates of perturbed time-varying linear systems, [14]. For this setup it is also necessary to study discrete-time systems with a transition map
Remarks on universal nonsingular controls for discrete-time systems Eduardo D. Sontag a and Fabian R. Wirth b;1 a Department of Mathematics, Rutgers University, New Brunswick, NJ 08903, b sontag@hilbert.rutgers.edu
More informationComments on integral variants of ISS 1
Systems & Control Letters 34 (1998) 93 1 Comments on integral variants of ISS 1 Eduardo D. Sontag Department of Mathematics, Rutgers University, Piscataway, NJ 8854-819, USA Received 2 June 1997; received
More informationThe Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria
ESI The Erwin Schrodinger International Boltzmanngasse 9 Institute for Mathematical Physics A-19 Wien, Austria The Negative Discrete Spectrum of a Class of Two{Dimentional Schrodinger Operators with Magnetic
More informationOn It^o's formula for multidimensional Brownian motion
On It^o's formula for multidimensional Brownian motion by Hans Follmer and Philip Protter Institut fur Mathemati Humboldt-Universitat D-99 Berlin Mathematics and Statistics Department Purdue University
More informationMarkov Chain Monte Carlo. Simulated Annealing.
Aula 10. Simulated Annealing. 0 Markov Chain Monte Carlo. Simulated Annealing. Anatoli Iambartsev IME-USP Aula 10. Simulated Annealing. 1 [RC] Stochastic search. General iterative formula for optimizing
More informationPCMI LECTURE NOTES ON PROPERTY (T ), EXPANDER GRAPHS AND APPROXIMATE GROUPS (PRELIMINARY VERSION)
PCMI LECTURE NOTES ON PROPERTY (T ), EXPANDER GRAPHS AND APPROXIMATE GROUPS (PRELIMINARY VERSION) EMMANUEL BREUILLARD 1. Lecture 1, Spectral gaps for infinite groups and non-amenability The final aim of
More information2 optimal prices the link is either underloaded or critically loaded; it is never overloaded. For the social welfare maximization problem we show that
1 Pricing in a Large Single Link Loss System Costas A. Courcoubetis a and Martin I. Reiman b a ICS-FORTH and University of Crete Heraklion, Crete, Greece courcou@csi.forth.gr b Bell Labs, Lucent Technologies
More informationA NOTE ON THE POINCARÉ AND CHEEGER INEQUALITIES FOR SIMPLE RANDOM WALK ON A CONNECTED GRAPH. John Pike
A NOTE ON THE POINCARÉ AND CHEEGER INEQUALITIES FOR SIMPLE RANDOM WALK ON A CONNECTED GRAPH John Pike Department of Mathematics University of Southern California jpike@usc.edu Abstract. In 1991, Persi
More informationLecture two. January 17, 2019
Lecture two January 17, 2019 We will learn how to solve rst-order linear equations in this lecture. Example 1. 1) Find all solutions satisfy the equation u x (x, y) = 0. 2) Find the solution if we know
More informationLOGARITHMIC MULTIFRACTAL SPECTRUM OF STABLE. Department of Mathematics, National Taiwan University. Taipei, TAIWAN. and. S.
LOGARITHMIC MULTIFRACTAL SPECTRUM OF STABLE OCCUPATION MEASURE Narn{Rueih SHIEH Department of Mathematics, National Taiwan University Taipei, TAIWAN and S. James TAYLOR 2 School of Mathematics, University
More informationMARKOV CHAINS: STATIONARY DISTRIBUTIONS AND FUNCTIONS ON STATE SPACES. Contents
MARKOV CHAINS: STATIONARY DISTRIBUTIONS AND FUNCTIONS ON STATE SPACES JAMES READY Abstract. In this paper, we rst introduce the concepts of Markov Chains and their stationary distributions. We then discuss
More informationA Representation of Excessive Functions as Expected Suprema
A Representation of Excessive Functions as Expected Suprema Hans Föllmer & Thomas Knispel Humboldt-Universität zu Berlin Institut für Mathematik Unter den Linden 6 10099 Berlin, Germany E-mail: foellmer@math.hu-berlin.de,
More informationStrassen s law of the iterated logarithm for diusion processes for small time 1
Stochastic Processes and their Applications 74 (998) 9 Strassen s law of the iterated logarithm for diusion processes for small time Lucia Caramellino Dipartimento di Matematica, Universita degli Studi
More informationq-alg/ Mar 96
Integrality of Two Variable Kostka Functions Friedrich Knop* Department of Mathematics, Rutgers University, New Brunswick NJ 08903, USA knop@math.rutgers.edu 1. Introduction q-alg/9603027 29 Mar 96 Macdonald
More informationQuantum logics with given centres and variable state spaces Mirko Navara 1, Pavel Ptak 2 Abstract We ask which logics with a given centre allow for en
Quantum logics with given centres and variable state spaces Mirko Navara 1, Pavel Ptak 2 Abstract We ask which logics with a given centre allow for enlargements with an arbitrary state space. We show in
More information2. The averaging process
2. The averaging process David Aldous July 10, 2012 The two models in lecture 1 has only rather trivial interaction. In this lecture I discuss a simple FMIE process with less trivial interaction. It is
More informationGarrett: `Bernstein's analytic continuation of complex powers' 2 Let f be a polynomial in x 1 ; : : : ; x n with real coecients. For complex s, let f
1 Bernstein's analytic continuation of complex powers c1995, Paul Garrett, garrettmath.umn.edu version January 27, 1998 Analytic continuation of distributions Statement of the theorems on analytic continuation
More informationWITH REPULSIVE SINGULAR FORCES MEIRONG ZHANG. (Communicated by Hal L. Smith) of repulsive type in the sense that G(u)! +1 as u! 0.
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 127, Number 2, February 1999, Pages 41{47 S 2-9939(99)512-5 PERIODIC SOLUTIONS OF DAMPED DIFFERENTIAL SYSTEMS WITH REPULSIVE SINGULAR FORCES MEIRONG
More informationON THE ZERO-ONE LAW AND THE LAW OF LARGE NUMBERS FOR RANDOM WALK IN MIXING RAN- DOM ENVIRONMENT
Elect. Comm. in Probab. 10 (2005), 36 44 ELECTRONIC COMMUNICATIONS in PROBABILITY ON THE ZERO-ONE LAW AND THE LAW OF LARGE NUMBERS FOR RANDOM WALK IN MIXING RAN- DOM ENVIRONMENT FIRAS RASSOUL AGHA Department
More informationThe Pacic Institute for the Mathematical Sciences http://www.pims.math.ca pims@pims.math.ca Surprise Maximization D. Borwein Department of Mathematics University of Western Ontario London, Ontario, Canada
More informationLecture 5: Random Walks and Markov Chain
Spectral Graph Theory and Applications WS 20/202 Lecture 5: Random Walks and Markov Chain Lecturer: Thomas Sauerwald & He Sun Introduction to Markov Chains Definition 5.. A sequence of random variables
More informationA Chaotic Phenomenon in the Power Swing Equation Umesh G. Vaidya R. N. Banavar y N. M. Singh March 22, 2000 Abstract Existence of chaotic dynamics in
A Chaotic Phenomenon in the Power Swing Equation Umesh G. Vaidya R. N. Banavar y N. M. Singh March, Abstract Existence of chaotic dynamics in the classical swing equations of a power system of three interconnected
More informationThe Erwin Schrodinger International Pasteurgasse 6/7. Institute for Mathematical Physics A-1090 Wien, Austria
ESI The Erwin Schrodinger International Pasteurgasse 6/7 Institute for Mathematical Physics A-1090 Wien, Austria On the Point Spectrum of Dirence Schrodinger Operators Vladimir Buslaev Alexander Fedotov
More informationTHE LUBRICATION APPROXIMATION FOR THIN VISCOUS FILMS: REGULARITY AND LONG TIME BEHAVIOR OF WEAK SOLUTIONS A.L. BERTOZZI AND M. PUGH.
THE LUBRICATION APPROXIMATION FOR THIN VISCOUS FILMS: REGULARITY AND LONG TIME BEHAVIOR OF WEAK SOLUTIONS A.L. BERTOI AND M. PUGH April 1994 Abstract. We consider the fourth order degenerate diusion equation
More informationA slow transient diusion in a drifted stable potential
A slow transient diusion in a drifted stable potential Arvind Singh Université Paris VI Abstract We consider a diusion process X in a random potential V of the form V x = S x δx, where δ is a positive
More informationCOMPARISON THEOREMS FOR THE SPECTRAL GAP OF DIFFUSIONS PROCESSES AND SCHRÖDINGER OPERATORS ON AN INTERVAL. Ross G. Pinsky
COMPARISON THEOREMS FOR THE SPECTRAL GAP OF DIFFUSIONS PROCESSES AND SCHRÖDINGER OPERATORS ON AN INTERVAL Ross G. Pinsky Department of Mathematics Technion-Israel Institute of Technology Haifa, 32000 Israel
More information76 Griesemer, Lewis, Siedentop the Dirac operator in the gap with the number of negative eigenvalues of a corresponding Schrodinger operator (see [5])
Doc. Math. J. DMV 75 A Minimax Principle for Eigenvalues in Spectral Gaps: Dirac Operators with Coulomb Potentials Marcel Griesemer, Roger T. Lewis, Heinz Siedentop Received: March 9, 999 Communicated
More informationSIMU L TED ATED ANNEA L NG ING
SIMULATED ANNEALING Fundamental Concept Motivation by an analogy to the statistical mechanics of annealing in solids. => to coerce a solid (i.e., in a poor, unordered state) into a low energy thermodynamic
More informationSEMIGROUP APPROACH FOR PARTIAL DIFFERENTIAL EQUATIONS OF EVOLUTION
SEMIGROUP APPROACH FOR PARTIAL DIFFERENTIAL EQUATIONS OF EVOLUTION Istanbul Kemerburgaz University Istanbul Analysis Seminars 24 October 2014 Sabanc University Karaköy Communication Center 1 2 3 4 5 u(x,
More informationMarkov processes Course note 2. Martingale problems, recurrence properties of discrete time chains.
Institute for Applied Mathematics WS17/18 Massimiliano Gubinelli Markov processes Course note 2. Martingale problems, recurrence properties of discrete time chains. [version 1, 2017.11.1] We introduce
More informationA note on adiabatic theorem for Markov chains
Yevgeniy Kovchegov Abstract We state and prove a version of an adiabatic theorem for Markov chains using well known facts about mixing times. We extend the result to the case of continuous time Markov
More informationMaximum principle for the fractional diusion equations and its applications
Maximum principle for the fractional diusion equations and its applications Yuri Luchko Department of Mathematics, Physics, and Chemistry Beuth Technical University of Applied Sciences Berlin Berlin, Germany
More informationPARAMETER IDENTIFICATION IN THE FREQUENCY DOMAIN. H.T. Banks and Yun Wang. Center for Research in Scientic Computation
PARAMETER IDENTIFICATION IN THE FREQUENCY DOMAIN H.T. Banks and Yun Wang Center for Research in Scientic Computation North Carolina State University Raleigh, NC 7695-805 Revised: March 1993 Abstract In
More informationA Type of Shannon-McMillan Approximation Theorems for Second-Order Nonhomogeneous Markov Chains Indexed by a Double Rooted Tree
Caspian Journal of Applied Mathematics, Ecology and Economics V. 2, No, 204, July ISSN 560-4055 A Type of Shannon-McMillan Approximation Theorems for Second-Order Nonhomogeneous Markov Chains Indexed by
More informationThe Best Circulant Preconditioners for Hermitian Toeplitz Systems II: The Multiple-Zero Case Raymond H. Chan Michael K. Ng y Andy M. Yip z Abstract In
The Best Circulant Preconditioners for Hermitian Toeplitz Systems II: The Multiple-ero Case Raymond H. Chan Michael K. Ng y Andy M. Yip z Abstract In [0, 4], circulant-type preconditioners have been proposed
More informationComputing the Initial Temperature of Simulated Annealing
Computational Optimization and Applications, 29, 369 385, 2004 c 2004 Kluwer Academic Publishers. Manufactured in he Netherlands. Computing the Initial emperature of Simulated Annealing WALID BEN-AMEUR
More informationGerhard Rein. Mathematisches Institut der Universitat Munchen, Alan D. Rendall. Max-Planck-Institut fur Astrophysik, and
A regularity theorem for solutions of the spherically symmetric Vlasov-Einstein system Gerhard Rein Mathematisches Institut der Universitat Munchen, Theresienstr. 39, 80333 Munchen, Germany, Alan D. Rendall
More informationsemi-group and P x the law of starting at x 2 E. We consider a branching mechanism : R +! R + of the form (u) = au + bu 2 + n(dr)[e ru 1 + ru]; (;1) w
M.S.R.I. October 13-17 1997. Superprocesses and nonlinear partial dierential equations We would like to describe some links between superprocesses and some semi-linear partial dierential equations. Let
More informationBoundedly complete weak-cauchy basic sequences in Banach spaces with the PCP
Journal of Functional Analysis 253 (2007) 772 781 www.elsevier.com/locate/jfa Note Boundedly complete weak-cauchy basic sequences in Banach spaces with the PCP Haskell Rosenthal Department of Mathematics,
More informationReversible Markov chains
Reversible Markov chains Variational representations and ordering Chris Sherlock Abstract This pedagogical document explains three variational representations that are useful when comparing the efficiencies
More informationUniversity of Toronto Department of Statistics
Norm Comparisons for Data Augmentation by James P. Hobert Department of Statistics University of Florida and Jeffrey S. Rosenthal Department of Statistics University of Toronto Technical Report No. 0704
More informationON THE ASYMPTOTIC STABILITY IN TERMS OF TWO MEASURES FOR FUNCTIONAL DIFFERENTIAL EQUATIONS. G. Makay
ON THE ASYMPTOTIC STABILITY IN TERMS OF TWO MEASURES FOR FUNCTIONAL DIFFERENTIAL EQUATIONS G. Makay Student in Mathematics, University of Szeged, Szeged, H-6726, Hungary Key words and phrases: Lyapunov
More informationStability, Queue Length and Delay of Deterministic and Stochastic Queueing Networks Cheng-Shang Chang IBM Research Division T.J. Watson Research Cente
Stability, Queue Length and Delay of Deterministic and Stochastic Queueing Networks Cheng-Shang Chang IBM Research Division T.J. Watson Research Center P.O. Box 704 Yorktown Heights, NY 10598 cschang@watson.ibm.com
More informationThe average dimension of the hull of cyclic codes
Discrete Applied Mathematics 128 (2003) 275 292 www.elsevier.com/locate/dam The average dimension of the hull of cyclic codes Gintaras Skersys Matematikos ir Informatikos Fakultetas, Vilniaus Universitetas,
More informationExtremal Cases of the Ahlswede-Cai Inequality. A. J. Radclie and Zs. Szaniszlo. University of Nebraska-Lincoln. Department of Mathematics
Extremal Cases of the Ahlswede-Cai Inequality A J Radclie and Zs Szaniszlo University of Nebraska{Lincoln Department of Mathematics 810 Oldfather Hall University of Nebraska-Lincoln Lincoln, NE 68588 1
More informationKirsten Lackner Solberg. Dept. of Math. and Computer Science. Odense University, Denmark
Inference Systems for Binding Time Analysis Kirsten Lackner Solberg Dept. of Math. and Computer Science Odense University, Denmark e-mail: kls@imada.ou.dk June 21, 1993 Contents 1 Introduction 4 2 Review
More informationCompeting sources of variance reduction in parallel replica Monte Carlo, and optimization in the low temperature limit
Competing sources of variance reduction in parallel replica Monte Carlo, and optimization in the low temperature limit Paul Dupuis Division of Applied Mathematics Brown University IPAM (J. Doll, M. Snarski,
More information1. Introduction. Consider a single cell in a mobile phone system. A \call setup" is a request for achannel by an idle customer presently in the cell t
Heavy Trac Limit for a Mobile Phone System Loss Model Philip J. Fleming and Alexander Stolyar Motorola, Inc. Arlington Heights, IL Burton Simon Department of Mathematics University of Colorado at Denver
More informationDynamical Systems. August 13, 2013
Dynamical Systems Joshua Wilde, revised by Isabel Tecu, Takeshi Suzuki and María José Boccardi August 13, 2013 Dynamical Systems are systems, described by one or more equations, that evolve over time.
More informationOn averaging principle for diusion processes with null-recurrent fast component
Stochastic Processes and their Applications 93 (21) 229 24 www.elsevier.com/locate/spa On averaging principle for diusion processes with null-recurrent fast component R. Khasminskii a; ;1, N. Krylov b;2
More informationSBS Chapter 2: Limits & continuity
SBS Chapter 2: Limits & continuity (SBS 2.1) Limit of a function Consider a free falling body with no air resistance. Falls approximately s(t) = 16t 2 feet in t seconds. We already know how to nd the average
More informationLatent voter model on random regular graphs
Latent voter model on random regular graphs Shirshendu Chatterjee Cornell University (visiting Duke U.) Work in progress with Rick Durrett April 25, 2011 Outline Definition of voter model and duality with
More informationMARY ANN HORN be decoupled into three wave equations. Thus, we would hope that results, analogous to those available for the wave equation, would hold
Journal of Mathematical Systems, Estimation, and Control Vol. 8, No. 2, 1998, pp. 1{11 c 1998 Birkhauser-Boston Sharp Trace Regularity for the Solutions of the Equations of Dynamic Elasticity Mary Ann
More informationMGMT 69000: Topics in High-dimensional Data Analysis Falll 2016
MGMT 69000: Topics in High-dimensional Data Analysis Falll 2016 Lecture 14: Information Theoretic Methods Lecturer: Jiaming Xu Scribe: Hilda Ibriga, Adarsh Barik, December 02, 2016 Outline f-divergence
More informationLecture 12: Detailed balance and Eigenfunction methods
Miranda Holmes-Cerfon Applied Stochastic Analysis, Spring 2015 Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility),
More informationE.G. KALNINS AND WILLARD MILLER, JR. The notation used for -series and -integrals in this paper follows that of Gasper and Rahman [3].. A generalizati
A NOTE ON TENSOR PRODUCTS OF -ALGEBRA REPRESENTATIONS AND ORTHOGONAL POLYNOMIALS E.G. KALNINSy AND WILLARD MILLER, Jr.z Abstract. We work out examples of tensor products of distinct generalized s`) algebras
More information2 BAISHENG YAN When L =,it is easily seen that the set K = coincides with the set of conformal matrices, that is, K = = fr j 0 R 2 SO(n)g: Weakly L-qu
RELAXATION AND ATTAINMENT RESULTS FOR AN INTEGRAL FUNCTIONAL WITH UNBOUNDED ENERGY-WELL BAISHENG YAN Abstract. Consider functional I(u) = R jjdujn ; L det Duj dx whose energy-well consists of matrices
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationLinear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) 1.1 The Formal Denition of a Vector Space
Linear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) Contents 1 Vector Spaces 1 1.1 The Formal Denition of a Vector Space.................................. 1 1.2 Subspaces...................................................
More informationSimulated Annealing for Constrained Global Optimization
Monte Carlo Methods for Computation and Optimization Final Presentation Simulated Annealing for Constrained Global Optimization H. Edwin Romeijn & Robert L.Smith (1994) Presented by Ariel Schwartz Objective
More informationON THE DIAMETER OF THE ATTRACTOR OF AN IFS Serge Dubuc Raouf Hamzaoui Abstract We investigate methods for the evaluation of the diameter of the attrac
ON THE DIAMETER OF THE ATTRACTOR OF AN IFS Serge Dubuc Raouf Hamzaoui Abstract We investigate methods for the evaluation of the diameter of the attractor of an IFS. We propose an upper bound for the diameter
More informationRealization of set functions as cut functions of graphs and hypergraphs
Discrete Mathematics 226 (2001) 199 210 www.elsevier.com/locate/disc Realization of set functions as cut functions of graphs and hypergraphs Satoru Fujishige a;, Sachin B. Patkar b a Division of Systems
More informationG METHOD IN ACTION: FROM EXACT SAMPLING TO APPROXIMATE ONE
G METHOD IN ACTION: FROM EXACT SAMPLING TO APPROXIMATE ONE UDREA PÄUN Communicated by Marius Iosifescu The main contribution of this work is the unication, by G method using Markov chains, therefore, a
More informationLEBESGUE INTEGRATION. Introduction
LEBESGUE INTEGATION EYE SJAMAA Supplementary notes Math 414, Spring 25 Introduction The following heuristic argument is at the basis of the denition of the Lebesgue integral. This argument will be imprecise,
More informationIsoperimetric inequalities for cartesian products of graphs
Isoperimetric inequalities for cartesian products of graphs F. R. K. Chung University of Pennsylvania Philadelphia 19104 Prasad Tetali School of Mathematics Georgia Inst. of Technology Atlanta GA 3033-0160
More informationDeterministic Dynamic Programming
Deterministic Dynamic Programming 1 Value Function Consider the following optimal control problem in Mayer s form: V (t 0, x 0 ) = inf u U J(t 1, x(t 1 )) (1) subject to ẋ(t) = f(t, x(t), u(t)), x(t 0
More informationFACTOR MAPS BETWEEN TILING DYNAMICAL SYSTEMS KARL PETERSEN. systems which cannot be achieved by working within a nite window. By. 1.
FACTOR MAPS BETWEEN TILING DYNAMICAL SYSTEMS KARL PETERSEN Abstract. We show that there is no Curtis-Hedlund-Lyndon Theorem for factor maps between tiling dynamical systems: there are codes between such
More informationInstitute for Advanced Computer Studies. Department of Computer Science. On Markov Chains with Sluggish Transients. G. W. Stewart y.
University of Maryland Institute for Advanced Computer Studies Department of Computer Science College Park TR{94{77 TR{3306 On Markov Chains with Sluggish Transients G. W. Stewart y June, 994 ABSTRACT
More information3 Stability and Lyapunov Functions
CDS140a Nonlinear Systems: Local Theory 02/01/2011 3 Stability and Lyapunov Functions 3.1 Lyapunov Stability Denition: An equilibrium point x 0 of (1) is stable if for all ɛ > 0, there exists a δ > 0 such
More informationSTOCHASTIC DIFFERENTIAL EQUATIONS WITH EXTRA PROPERTIES H. JEROME KEISLER. Department of Mathematics. University of Wisconsin.
STOCHASTIC DIFFERENTIAL EQUATIONS WITH EXTRA PROPERTIES H. JEROME KEISLER Department of Mathematics University of Wisconsin Madison WI 5376 keisler@math.wisc.edu 1. Introduction The Loeb measure construction
More informationA global solution curve for a class of free boundary value problems arising in plasma physics
A global solution curve for a class of free boundary value problems arising in plasma physics Philip Korman epartment of Mathematical Sciences University of Cincinnati Cincinnati Ohio 4522-0025 Abstract
More informationDiscrete Halanay-type inequalities and applications
Nonlinear Analysis 55 (2003) 669 678 www.elsevier.com/locate/na Discrete Halanay-type inequalities and applications Eduardo Liz a;, Anatoli Ivanov b, Juan Bosco Ferreiro c a Departamento de Matematica
More informationSupplement A: Construction and properties of a reected diusion
Supplement A: Construction and properties of a reected diusion Jakub Chorowski 1 Construction Assumption 1. For given constants < d < D let the pair (σ, b) Θ, where Θ : Θ(d, D) {(σ, b) C 1 (, 1]) C 1 (,
More informationYadernaya Fizika, 56 (1993) Soviet Journal of Nuclear Physics, vol. 56, No 1, (1993)
Start of body part published in Russian in: Yadernaya Fizika, 56 (1993) 45-5 translated in English in: Soviet Journal of Nuclear Physics, vol. 56, No 1, (1993) A METHOD FOR CALCULATING THE HEAT KERNEL
More informationLecture Notes on Partial Dierential Equations (PDE)/ MaSc 221+MaSc 225
Lecture Notes on Partial Dierential Equations (PDE)/ MaSc 221+MaSc 225 Dr. Asmaa Al Themairi Assistant Professor a a Department of Mathematical sciences, University of Princess Nourah bint Abdulrahman,
More informationINRIA Rocquencourt, Le Chesnay Cedex (France) y Dept. of Mathematics, North Carolina State University, Raleigh NC USA
Nonlinear Observer Design using Implicit System Descriptions D. von Wissel, R. Nikoukhah, S. L. Campbell y and F. Delebecque INRIA Rocquencourt, 78 Le Chesnay Cedex (France) y Dept. of Mathematics, North
More information