A short introduction to random trees

Size: px
Start display at page:

Download "A short introduction to random trees"

Transcription

1 Mongolian Mahemaical Journal (6) 5 A shor inroducion o random rees Chrisina Goldschmid Noes from wo lecures given a he CIMPA-DAAD Research School on Sochasic Processes and heir Applicaions, Naional Universiy of Mongolia, Ulaanbaaar, Mongolia, h July h Augus 5 Absrac Tree srucures are very common in naure, and we ofen wan o model rees which are somehow randomly generaed. In mahemaics, we hink of rees as conneced graphs (or neworks) wih no cycles. In hese lecures, I will discuss some simple models of random rees and wha we can say abou large insances of hem. The firs lecure will focus on Galon Wason branching processes. We hink of an idealised biological populaion (say of cells) in which each individual lives for a uni of ime, gives birh o a random number of children (wih a given disribuion) and hen dies. Wha can we say abou he way he populaion size evolves? Or abou he family ree of ha populaion? In he second lecure, we will ake a more combinaorial perspecive. Wha does a ree wih n verices, chosen uniformly a random from all he possibiliies, look like for large n? I urns ou ha we can use Galon Wason branching processes o help us answer his quesion. Along he way, we will encouner several beauiful bis of combinaorics and probabiliy, including Cayley s formula, random walks and (finally) Brownian moion. Conens. Branching processes Inroducion Generaing funcions Exincion Maringales Galon Wason rees Relaionship o random walks Uniform random rees Inroducion Cayley s formula Leaves The uniform random ree as a Galon Wason ree Taking limis The Brownian coninuum random ree References Mahemaics Subjec Classificaion 6C5, 6J. The auhor is graeful o Lady Margare Hall, Oxford, for funding which enabled her o ravel o Mongolia o deliver hese lecures.

2 54 CHRISTINA GOLDSCHMIDT Figure. The family ree of a branching process. We sar a he boom of he ree, wih a single individual in generaion. Then here are individuals in each of generaions and, 5 individuals in generaion, a single individual in generaion 4 and no individuals in subsequen generaions. (Noice ha in order o produce his picure, we mus have p() >.).. Inroducion. Branching processes Suppose we have a populaion (say of baceria). Each individual in he populaion lives a uni ime and, jus before dying, gives birh o a random number of children in he nex generaion. This number of children has probabiliy mass funcion p(k), k, called he offspring disribuion i.e. p(k) is he probabiliy ha an individual has k children. Differen individuals reproduce independenly in he same manner. Figure shows a possible family ree of such a populaion. Le X n be he size of he populaion in generaion n. The process (X n ) n is ofen called a Galon Wason branching process afer i was popularised by Galon and Wason [], alhough had already been sudied (much earlier) by Bienaymé [5]. As a simple model of he evoluion of a populaion, i has proved o be enormously powerful. I is now considered o be one of he fundamenal processes of probabiliy heory... Generaing funcions Suppose ha X =. Le C (n) i be he number of children of he ih individual in generaion n, so ha we may wrie X n+ = C (n) + C (n) + + C (n) X n. (We inerpre his sum as if X n =.) Noe ha C (n), C (n),... are independen and idenically disribued. Consider he probabiliy generaing funcions G(s) = k= p(k)sk (for s R such ha he righ-hand side is absoluely convergen) and G n (s) = E [ ] s Xn for n.

3 A SHORT INTRODUCTION TO RANDOM TREES 55 Proposiion.. For n, G n+ (s) = G n (G(s)) = G(G(... G(s)...)) = G(G n (s)). }{{} n+ imes Proof. Since X =, we have G (s) = s. Also, we ge X = C () which has probabiliy mass funcion p(k), k. So G (s) = E [ s X ] = G(s). Since we ge G n+ (s) = E [ s Xn+] [ = E X n X n+ = C (n) i, = k= i= s Xn i= C(n) i [ E s Xn ] i= C(n) i ] Xn = k P (X n = k) and, since C (n), C (n),... are independen of X n, his is equal o k= E [s ] k i= C(n) i P (X n = k). Since C (n), C (n),... are independen and idenically disribued, his in urn equals Hence, by inducion, for n, ] k E [s C(n) P (Xn = k) k= = (G(s)) k P (X n = k) k= = G n (G(s)). G n (s) = G(G(... G(s)...)) = G(G n (s)). }{{} n imes Corollary.. Suppose ha he mean number of children of a single individual is µ i.e. k= kp(k) = µ, for some µ <. Then E [X n ] = µ n. Proof. Noe ha µ = G ( ) and ha E [X n ] = G n( ). By he chain rule, Leing s, we ge G n(s) = d ds G(G n (s)) = G n (s)g (G n (s)). E [X n ] = E [X n ] µ = = µ n.

4 56 CHRISTINA GOLDSCHMIDT Exercise.. Show ha if σ := var (X ) < hen var (X n ) = { σ µ n (µ n ) µ if µ nσ if µ =. Corollary. ells us ha we ge exponenial growh on average if µ > and exponenial decrease if µ <. This raises an ineresing quesion: does he populaion die ou?.. Exincion In his secion, we invesigae he exincion probabiliy: q := P (populaion dies ou). If p() = hen every individual has a leas one child and so he populaion clearly grows forever. If p() >, on he oher hand, hen he populaion dies ou wih posiive probabiliy because q P (X = ) = p() >. Noice ha his holds even in he cases where E [X n ] grows as n. Proposiion.4. We have q = lim n P (X n = ). Proof. The even ha he populaion dies ou can be wrien as {X m = }. m= Noice ha he evens in his union are no disjoin: indeed, X m = implies ha X n = for all n > m. I follows ha n {X m = } = {X n = }. m= m= We have ha n m= {X m = } n+ m= {X m = } m= {X m = } and so by monoone convergence, ( ) ( n ) P {X m = } = lim P {X m = } = lim P (X n = ). n n m= Example.5. Suppose ha p(k) = (/) k+, k, so ha each individual has a geomeric number of offspring. We ll find he disribuion of X n. Firs calculae ( ) k+ G(s) = s k = s. k=

5 A SHORT INTRODUCTION TO RANDOM TREES 5 By plugging his ino iself a couple of imes, we ge G (s) = s s, G (s) = s 4 s. A naural guess is ha G n (s) = n (n )s (n+) ns which is, in fac, he case, as can be proved by inducion. If we wan he probabiliy mass funcion of X n, we need o expand his quaniy ou in powers of s. We have (n + ) ns = (n + ) Muliplying by n (n )s, we ge n k+ s k G n (s) = (n + ) k+ k= k= ( ns/(n + )) = k= n k s k (n + ) k+. n k (n )s k (n + ) k = n n + + We can read off he coefficiens now o see ha { n n+ if k = P (X n = k) = n k if k. (n+) k+ k= n k s k (n + ) k+. Noice ha P (X n = ) as n, which indicaes ha he populaion dies ou evenually in his case. Exercise.6. Suppose now ha p(k) = p( p) k, k, for some p (, ). (a) Find he probabiliy generaing funcion G(s) of his disribuion. Wha is is mean, µ? (b) Suppose ha p /. Prove, by inducion, ha he probabiliy generaing funcion of X n is G n (s) = p [( p)n p n ] ( p)s[( p) n p n ] [( p) n+ p n+ ] ( p)s[( p) n p n ], for n. { (c) Find G n () = P (X n = ) and hence show ha q = min p p }., Deduce ha he populaion always dies ou if µ whereas i has posiive probabiliy of surviving forever if µ >. Le us now reurn o he general case. We can find an equaion saisfied by q by condiioning on he number of children of he firs individual. q = P (populaion dies ou X = k) P (X = k) = k= P (populaion dies ou X = k) p(k). k= Remember ha each of he k individuals in he firs generaion behaves exacly like he paren. In paricular, we can hink of each of hem saring is own family, which is an independen copy of he original family. Moreover, he whole populaion dies ou if and only if all of hese sub-populaions die ou. If we had k families, his occurs wih

6 5 CHRISTINA GOLDSCHMIDT probabiliy q k. So q = q k p(k) = G(q). k= Exercise.. Check ha for he offspring disribuion in Exercise.6 he exincion probabiliy q you obained does, indeed, solve s = G(s). The equaion q = G(q) doesn quie enable us o deermine q: noice ha is always a soluion, bu i s no necessarily he only soluion in [, ]. I urns ou ha when here are muliple soluions, he one we wan is he smalles one in [, ]. Theorem.. Suppose ha p() > and p() + p() <. (a) The equaion s = G(s) has a mos wo soluions in [, ]. The exincion probabiliy q is he smalles non-negaive roo of he equaion G(s) = s. (b) Suppose ha he offspring disribuion has mean µ. Then if µ hen q = if µ > hen q <. Proof. have We observe ha G() = p() >, G() = and ha µ = G ( ). We also G (s) = j(j )p(j)s j for s <. j= Since p() + p() <, we have G (s) > for s [, ). Hence, G is sricly convex. Consider he case µ. Then G() > and G (s) G ( ) = µ for all s [, ]. I follows ha G(s) > s for all s [, ) and so he unique roo in [, ] is s =, which mus necessarily be q. See Figure. Consider now he case µ >. Then since G() >, G() =, G ( ) > and G is sricly convex, here exiss precisely one roo α of s = G(s) in [, ). See Figure. Now G() G(α) = α. So G(G()) G(α) = α and, ieraing, G n () α for all n. I follows ha q = lim n G n () α. Bu since q mus be a roo of G(s) = s, we deduce ha q = α. We refer o he case µ < as subcriical, he case µ = as criical and he case µ > as supercriical..4. Maringales We have jus seen ha, in he subcriical and criical cases, we have X n as n. In he supercriical case, on he oher hand, here is posiive probabiliy ha X n. We also have ha he populaion grows exponenially on average. Wha can we say almos surely? For n, define W n = µ n X n and le F n = σ(x m, m n).

7 A SHORT INTRODUCTION TO RANDOM TREES 5 G(s) G(s) s s G(s) s Figure. Solving s = G(s) for µ <, µ = and µ > respecively. Proposiion.. (W n, n ) is a non-negaive maringale (wih respec o he filraion (F n ) n ) and so exiss almos surely. W := lim n W n Proof. We have already shown ha E [X n ] = µ n so ha E [W n ] = for all n. Hence, W n is inegrable. Moreover, [ Xn ] E [W n+ F n ] = µ (n+) E [X n+ F n ] = µ (n+) E F n [ = µ W n E i= C (n) C (n) i ],

8 6 CHRISTINA GOLDSCHMIDT since X n is F n -measurable, [ ] C (n) i, i are independen of F n and are idenically disribued. Since E = µ, we ge C (n) E [W n+ F n ] = W n, as required. The convergence resul now follows from he almos sure maringale convergence heorem. Proposiion. ells us ha, in a rough sense, W n µ n W. Using Faou s lemma, we obain ha E [W ] lim inf n E [W n ] =. I is, however, sill possible ha W. (Indeed, Proposiion. holds even if µ, in which case we mus have W.) Proposiion.. or o. Suppose ha < µ <. Then P (W = ) is eiher equal o q Proof. We have X X n+ = X n (i), where X n (), X n (),... are i.i.d. copies of X n, independen of X. So, similarly, i= X W n+ = µ W n (i). Leing n and using Proposiion., we see ha W has he same disribuion as µ X i= W (i). (.) i= ] In paricular, P (W = ) = E [P (W = ) Z = G(P (W = )). So P (W = ) is a roo of s = G(s). Bu, by Theorem. he possible roos are q and and he resul follows. Exercise.. Consider again he Geomeric offspring disribuion of Exercise.6. Assume ha p < /, so ha µ = p p >. Suppose ha W has he disribuion specified by P (W = ) = p p ( p)w, P (W > w) = p p e p for w >. Show ha W solves he recursive disribuional equaion (.). [Hin: Firs show ha he number of erms in he sum (.) which are non-zero has a Geomeric disribuion, wih parameer p/( p). Then use momen generaing funcions.] A necessary and sufficien condiion for E [W ] = is given in he following heorem, which we won prove.

9 A SHORT INTRODUCTION TO RANDOM TREES 6 Theorem. (Kesen and Sigum []). Suppose ha < µ <, and wrie log + x = max{log x, }. Then he following are equivalen: E [W ] = P (W > non-exincion) = k= kp(k) log+ k <..5. Galon Wason rees We have so far jus considered he populaion sizes X, X, X,... in each of he successive generaions. Bu i s also naural o wan o hink abou a richer srucure: he family ree (also called he genealogical ree). In order o do his, i will urn ou o be useful o consider a canonical labelling for rees (called he Ulam Harris labelling). Le N = {,,...} and le U = { } n= Nn. In general, an elemen u U is a sequence of naural numbers represening a poin in an infiniary ree. We will hink of u as a sring and wrie u u... u n. The concaenaion of srings u = u u... u n and v = v v... v m is wrien uv = u u... u n v v... v m. The label of a verex u U indicaes is genealogy: u has paren p(u) = u u... u n. u has children u, u,.... We wrie u = n for he generaion of u. Definiion.. A rooed ordered ree is a subse of U such ha for all u such ha u, p(u) for all u, here exiss c(u) {,,,...} such ha for j N, uj if and only if j c(u). The ree is rooed a. The number of children of u in is given by c(u). We wrie #() for he size (number of verices) of and noe ha #() = + u c(u). We can hink of he family ree of a branching process as a random rooed ordered ree which is obained by sampling independen and idenically disribued random variables C(u), u U, each having he offspring disribuion. The corresponding Galon Wason ree T is hen obained by finding he larges subse of U such ha he condiions of Definiion. are fulfilled. Noe ha he ree is finie if and only if he branching process becomes exinc..6. Relaionship o random walks Trees are slighly awkward objecs o manipulae and so we ofen find i useful o encode hem in erms of funcions. This will be useful in he nex lecure. There are several differen ways o do his we will concenrae on one which works well wih our labelling. See Secions. and. of Le Gall [] for more deails. For a rooed ordered ree, le v = and, for i, le v i be he ih verex of when he verices are lised in lexicographic order of label.

10 6 CHRISTINA GOLDSCHMIDT Definiion.4. The Lukasiewicz pah associaed wih a ree of size #() = n is he funcion l : {,,..., n} {,,,...} defined by l() = and, for i n, l(i) = l(i ) + c(v i ) i.e. l(i) = i j= (c(v j) ). Noe ha by (.), we mus have l(n) =. Moreover, l(i) for i n. Exercise.5. unique ree. Show ha a given Lukasiewicz pah (l(i), i n) encodes a In oher words, here is a bijecion beween rooed ordered rees and Lukasiewicz pahs. In he case of Galon Wason rees, his correspondence gives us somehing paricularly nice. Wrie (L(i), i #(T)) for he Lukasiewicz pah associaed wih a Galon Wason ree T. Le (S(i), i ) be a random walk wih sep-sizes having disribuion p(k + ), k and le τ = inf{i : S(i) = }. Proposiion.6. he same disribuion. The processes (L(i), i #(T)) and (S(i), i τ) have Proof (skech). We have L(i) = i j= (C(v j) ), so ha L(i) for i #(T) and L(#(T)) =. Since he random variables C(u), u U are independen, ake values in {,,,...} and have probabiliy mass funcion p(k + ), k, he resul seems inuiively clear. The problem is ha he labels v, v,... are random and depend on (C(u), u U), so i is no, in fac, obvious ha C(v ), C(v ),... are independen and idenically disribued. The poin is ha, a sep n + of he lexicographic exploraion of he ree, he labels v,..., v n and he numbers C(v ),..., C(v n ) enable us o deermine (a) if we have reached he end of he ree and (b) if no, which is he label v n+. Bu, in case (b), he informaion we already possess can influence he value of C(v n+ ). See Proposiion.5 of Le Gall [] for he deails. Noe ha we can (and will) hink of he random walk as coninuing on beyond he firs ime i his. Remark.. Proposiion.6 enails ha we can hink of he exincion probabiliy of a branching process as he probabiliy ha he random walk ever his. We observe ha he sep sizes of he random walk have mean µ. Assuming addiionally ha p() > and p() + p() <, we see ha i is, indeed, he case ha he walk his wih probabiliy if µ (corresponding o exincion) and has posiive probabiliy of never hiing if µ >.

11 A SHORT INTRODUCTION TO RANDOM TREES 6.. Inroducion. Uniform random rees In his lecure, we re going o ake a more combinaorial perspecive on random rees, bu i will evenually connec back o wha we discussed in he firs lecure. Here we will concenrae on labelled rees. We le T n be he se of (unordered) labelled rees on n verices (labelled by [n] := {,,..., n}. For example, T consiss of he rees Unordered means ha hese rees are all he same: bu his one is differen: 4 5 Suppose we ake all of he rees on n labelled verices and pick one uniformly a random call i T n. Wha does T n look like? Phrased anoher way: wha does a ypical ree on n verices look like?.. Cayley s formula A good saring place is o ask how many elemens T n has. Theorem. (Cayley s formula). For n, T n = n n. This heorem was firs proved by Borchard [6] and laer exended by Cayley []. Proof. We give a beauiful proof due o Piman [4]. Call he quaniy we re afer a(n). I s useful for he purposes of his proof o hink of rooed rees wih direced edges, where he edges are all direced away from he roo. (We can simply forge he roo and he direcions of he arrow a he end.) Sar from he empy graph on our n

12 64 CHRISTINA GOLDSCHMIDT verices and coun he number of differen sequences of direced edges ha we can add o form a rooed ree. Firs way o coun. Sar wih one of he a(n) unrooed rees on n verices. There are n ways o pick a roo. Direc all edges away from he roo. There are (n )! differen orders in which we migh add he n direced edges. So, we have a(n) n (n )! = a(n)n! differen sequences of direced edges. Second way o coun. Sar by hinking of each of he n verices as a roo. A each sep, we have a fores of rooed rees, wih edges direced away from he roos. Add an edge from an arbirary verex v o he roo of a differen ree, and roo he resuling ree a he roo of v s old ree. If we have added k edges already, hen we have a fores conaining n k rees (we reduce he number of rees by a each sep). So here are n(n k ) choices for he edge we add: he saring verex can be any one of he n verices and here are n k roos of oher rees o which i can connec. So he oal number of choices is n k= n(n k ) = n n (n )! = n n n!. So a(n)n! = n n n! and cancelling gives a(n) = n n. There are many quesions we migh ask abou our uniform random ree T n. For example, How many leaves (i.e. verices wih only one neighbour) are here? How many verices have wo neighbours? Or hree? Or more? Wha is he diameer of he ree (i.e. he lengh of he longes pah beween wo poins in he ree)? Pick a verex uniformly a random o be he roo. How many verices are here a disance d from he roo? Because he ree is random, all of hese quaniies are random variables... Leaves We firs invesigae he proporion of leaves in a uniform random ree. Theorem.. as n. Le N n denoe he number of leaves in T n. Then N n n p e Proof. where Noe ha I i = N n = n i= { if he verex labelled i is a leaf oherwise. I i

13 A SHORT INTRODUCTION TO RANDOM TREES 65 We have E [I i ] = P (i is a leaf). In order o calculae his probabiliy, we need o deermine how many rees on [n] have i as a leaf. Each such ree can be hough of as a ree on [n] \ {i} plus an edge from one of he oher verices o i. Since here are (n ) n rees on n labelled verices, i follows ha here are (n ) n rees which have i as a leaf. So P (i is a leaf) = (n )n n n = ( n) n e, as n. If he indicaor random variables I, I,..., I n were independen, we would now be able o apply he weak law of large numbers. Bu his is no quie he case. Le i j. Then P (i is a leaf and j is a leaf) = (n ) (n ) n 4 ( n n = n) n e, so ha I i and I j are asympoically independen. In paricular, ( cov (I i, I j ) = n) n ( n) (n ), as n. We also have ( var (I i ) = Hence, using symmery, var ( Nn n ( ) ) n ( ) n ( e )e. n n ) = n n var (I i ) + n cov (I i, I j ) i= i<j = n var (I ) + n n cov (I, I ), as n. Finally, for any ɛ >, we can apply Chebyshev s inequaliy o obain ( ) N n P n e > ɛ var (N n/n) ɛ, as n, as desired..4. The uniform random ree as a Galon Wason ree We now make he link beween our sudy of he uniform random ree and he las lecure on branching processes. Proposiion.. The uniform random ree T n has he same disribuion as a ree generaed as follows: Take a Galon Wason ree wih Poisson() offspring disribuion Condiion i o have oal progeny precisely n Assign he roo he label and he oher verices random labels chosen from {,..., n}. Forge he roo and he original ordering.

14 66 CHRISTINA GOLDSCHMIDT Proof. Le T be a Poisson() Galon Wason ree, and le be a paricular ree wih #() = n. Then P (T = ) = e c(v)! = e n c(v)!. v v Now observe ha P (#(T) = n) is a funcion only of n. Hence, P (T = #(T) = n) = f(n) v c(v)! for some funcion f. Moreover, here are (n )! differen ways o label wih he roo labelled and he oher labels chosen from {,..., n}. Now consider a paricular unordered labelled ree T n, and hink of i as rooed a. Then we can say which are he children of a verex v (i.e. hose which are furher away from he roo han he verex iself). Then here are v c(v)! differen possible orderings of he labelled ree. In paricular, here are v c(v)! differen labelled ordered rees which yield he same labelled ree when we forge he ordering. So he probabiliy ha our condiioned Galon Wason ree gives a paricular ree T n is c(v)! v f(n) (n )! v c(v)! = f(n) (n )!. Since his depends only on n, and no on any oher feaure of he ree, i mus be he case ha he ree is uniformly disribued on T n. Remark.4. Oher classes of combinaorial rees may also be obained as Galon Wason rees condiioned on heir oal progeny. For example, aking p(k) = ( ) k+, k, gives a uniform plane ree aking p() = / and p() = / gives a uniform complee binary ree (as long as he oal progeny is odd)..5. Taking limis We ve already seen ha if we look a he proporion of he verices in T n which are leaves hen we obain a meaningful limi as n. Wha can we say more generally abou T n as n ges large? We will make crucial use of he link o Galon Wason rees oulined in he las secion. We won prove he resuls in his secion or he nex see Aldous [] and Le Gall [] for more deails. Recall ha if T is a Galon Wason ree wih offspring disribuion p(k), k, hen is Lukasiewicz pah (L(i), i #(T)) has he same disribuion as a random walk (S(k), k ) wih sep-size disribuion p(k + ), k, sopped a τ, he firs ime S his. We will concenrae on he criical case, where µ =. The fundamenal idea ha we will make use of in his secion comes from he following heorem. Theorem.5 (Donsker s heorem). Suppose ha (S(k), k ) is a random walk wih sep-sizes of mean and finie variance σ. Then, as n, ( ) σ n S( n ), d (B(), ),

15 A SHORT INTRODUCTION TO RANDOM TREES 6 Figure. A simulaion of a sandard Brownian excursion. where (B(), ) is a sandard Brownian moion. In oher words, in he criical case, here is a good way o rescale he random walk in such a way ha we obain a non-rivial limiing process. If we have a Poisson() offspring disribuion, we ge ha he sep-sizes precisely have mean and variance. Raher han a criical Galon Wason of random size, we would like o consider a criical Galon Wason ree condiioned o have size n. The corresponding Lukasiewicz pah is a random walk condiioned o hi for he firs ime a ime n. Le (S (n) (k), k n) have he same disribuion as (S(k), k n) condiioned on he even {τ = n}. Theorem.6 (Kaigh []). As n, ( ) σ n S(n) d ( n ), (e(), ), where (e(), ) is a sandard Brownian excursion. Informally, e is a Brownian moion condiioned o say posiive in he inerval (, ) and o come back o a ime (see Figure ). This heorem suggess ha here should be a limi for he ree iself, and ha limi should somehow be encoded by e. This is, indeed, he case, as was proved by Aldous in a series of papers in he early s [,, ]. The limi is known as he Brownian coninuum random ree..6. The Brownian coninuum random ree In order o see how he limi ree arises, we will go hrough he consrucion in he special case where he offspring disribuion is given by p(k) = ( ) k+, k i.e. when

16 6 CHRISTINA GOLDSCHMIDT we have a uniform plane ree. We will use an encoding for our ree which is somewha easier o visualise han he Lukasiewicz pah. The conour funcion is obained by simply racing he conour of he ree from lef o righ a speed, so ha we pass along each edge wice. Record he disance from he roo a each ime o ge (C(k), k n )

17 A SHORT INTRODUCTION TO RANDOM TREES

18 CHRISTINA GOLDSCHMIDT The conour funcion is a sor of expanded version of he ree. For convenience, we will define a slighly shifed version: le C() =, C(n) = and, for k n, C(k) = + C(k ). Exercise.. Show ha ( C(k), k n) has he same disribuion as a simple symmeric random walk condiioned o reurn o he origin for he firs ime a ime n. I s easy o linearly inerpolae o give a coninuous funcion C : [, n] [, ), and similarly for C. Since he variance of he offspring disribuion is, i follows from Theorem.6 ha n ( C(n), ) d (e(), ) as n and, hence, n (C(n), ) d (e(), ). I s also relaively sraighforward o see how o ge back a ree from a conour funcion. Inuiively speaking, we pu glue on he underside of he funcion and hen push he wo sides ogeher unil hey mee (see Figure 4). The Brownian coninuum random ree is he objec we obain by doing his gluing procedure o he funcion (e(), ) see Figure 5.

19 A SHORT INTRODUCTION TO RANDOM TREES Figure 4. Trees from excursions Figure 5. A simulaion of a Brownian coninuum random ree For a general criical Galon Wason ree wih finie offspring variance σ condiioned o have oal progeny n, i urns ou ha he convergence of he (suiably rescaled) Lukasiewicz pah o a sandard Brownian excursion acually implies he convergence of he conour funcion of he ree (alhough his is far from rivial o see). Hence, any such condiioned ree converges in disribuion as n o he Brownian coninuum random ree. See Aldous [] and Le Gall [] for more deails, in paricular he specificaion of he opology in which convergence occurs. Furher reading In hese lecures, I have jus scraped he surface of a large and fascinaing opic. For an accoun of he hisory of branching processes, see Kendall []. For much more on he heory of branching processes, see he classic accoun of Ahreya and Ney [4] or he recen lecure course of Shi [5]. To learn abou maringales, see he beauiful book of

20 A SHORT INTRODUCTION TO RANDOM TREES Williams [6]. The accoun I have presened here of he relaionship beween Galon Wason rees and random walks owes much o he excellen survey paper of Le Gall []. Tha paper also gives he rigorous deails of he convergence of a criical Galon Wason ree wih finie offspring variance o he Brownian coninuum random ree. For overview of he Brownian coninuum random ree, see Aldous aricle []. Finally, for an accoun of a variey of models of random rees, see he book of Drmoa []. Acknowledgemens The simulaion of a large Galon Wason ree (approximaing he Brownian coninuum random ree) in Figure 5 was done using Mahemaica code wrien by Igor Korchemski. I am very graeful o him for leing me use his program. I would like o express my graiude o he organisers of he research school in Ulaanbaaar a which hese lecures were given for heir ireless effors and wonderful hospialiy. Finally, I would like o hank he audience for heir warm welcome and infecious enhusiasm, which made he research school an unforgeable experience. References. D. Aldous, The coninuum random ree. I, Ann. Probab. (), no.,.., The coninuum random ree. II. An overview, Sochasic analysis (Durham, ), London Mah. Soc. Lecure Noe Ser., vol. 6, Cambridge Universiy Press, Cambridge,, pp..., The coninuum random ree. III, Ann. Probab. (), no., K. Ahreya and P. Ney, Branching processes, Springer-Verlag, New York,, Die Grundlehren der mahemaischen Wissenschafen, Band I. J. Bienaymé, De la loi de muliplicaion e de la durée des familles, Soc. Philomah. Paris Exrais, Sér. 5 (45),. 6. C. W. Borchard, Über eine der Inerpolaion ensprechende Darsellung der Eliminaions- Resulane, J. Reine Angew. Mah. 5 (6),.. A. Cayley, A heorem on rees, Q. J. Pure Appl. Mah. (), 6.. M. Drmoa, Random rees: An inerplay beween combinaorics and probabiliy, Springer-Verlag, Vienna,.. F. Galon and H.W. Wason, On he probabiliy of exincion of families, J. Anhropol. Ins. 4 (4), 44.. W. D. Kaigh, An invariance principle for random walk condiioned by a lae reurn o zero, Ann. Probabiliy 4 (6), no., 5.. D.G. Kendall, Branching processes since, J. London Mah. Soc. 4 (66), H. Kesen and B. Sigum, A limi heorem for mulidimensional Galon Wason processes, Ann. Mah. Sais. (66),.. J.-F. Le Gall, Random rees and applicaions, Probab. Surv. (5), J. Piman, Coalescen random foress, J. Combin. Theory Ser. A 5 (), no., Zhan Shi, Random walks and rees, X Symposium on Probabiliy and Sochasic Processes and he Firs Join Meeing France-Mexico of Probabiliy, ESAIM Proc., vol., EDP Sci., Les Ulis,, pp.. 6. David Williams, Probabiliy wih maringales, Cambridge Mahemaical Texbooks, Cambridge Universiy Press, Cambridge,. Chrisina Goldschmid Deparmen of Saisics and Lady Margare Hall Universiy of Oxford Unied Kingdom goldschm@sas.ox.ac.uk

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still. Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in

More information

The Arcsine Distribution

The Arcsine Distribution The Arcsine Disribuion Chris H. Rycrof Ocober 6, 006 A common heme of he class has been ha he saisics of single walker are ofen very differen from hose of an ensemble of walkers. On he firs homework, we

More information

Notes for Lecture 17-18

Notes for Lecture 17-18 U.C. Berkeley CS278: Compuaional Complexiy Handou N7-8 Professor Luca Trevisan April 3-8, 2008 Noes for Lecure 7-8 In hese wo lecures we prove he firs half of he PCP Theorem, he Amplificaion Lemma, up

More information

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t... Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger

More information

KINEMATICS IN ONE DIMENSION

KINEMATICS IN ONE DIMENSION KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec

More information

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes Half-Range Series 2.5 Inroducion In his Secion we address he following problem: Can we find a Fourier series expansion of a funcion defined over a finie inerval? Of course we recognise ha such a funcion

More information

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities: Mah 4 Eam Review Problems Problem. Calculae he 3rd Taylor polynomial for arcsin a =. Soluion. Le f() = arcsin. For his problem, we use he formula f() + f () + f ()! + f () 3! for he 3rd Taylor polynomial

More information

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals

More information

Unit Root Time Series. Univariate random walk

Unit Root Time Series. Univariate random walk Uni Roo ime Series Univariae random walk Consider he regression y y where ~ iid N 0, he leas squares esimae of is: ˆ yy y y yy Now wha if = If y y hen le y 0 =0 so ha y j j If ~ iid N 0, hen y ~ N 0, he

More information

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

More information

5. Stochastic processes (1)

5. Stochastic processes (1) Lec05.pp S-38.45 - Inroducion o Teleraffic Theory Spring 2005 Conens Basic conceps Poisson process 2 Sochasic processes () Consider some quaniy in a eleraffic (or any) sysem I ypically evolves in ime randomly

More information

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution Physics 7b: Saisical Mechanics Fokker-Planck Equaion The Langevin equaion approach o he evoluion of he velociy disribuion for he Brownian paricle migh leave you uncomforable. A more formal reamen of his

More information

We just finished the Erdős-Stone Theorem, and ex(n, F ) (1 1/(χ(F ) 1)) ( n

We just finished the Erdős-Stone Theorem, and ex(n, F ) (1 1/(χ(F ) 1)) ( n Lecure 3 - Kövari-Sós-Turán Theorem Jacques Versraëe jacques@ucsd.edu We jus finished he Erdős-Sone Theorem, and ex(n, F ) ( /(χ(f ) )) ( n 2). So we have asympoics when χ(f ) 3 bu no when χ(f ) = 2 i.e.

More information

Biol. 356 Lab 8. Mortality, Recruitment, and Migration Rates

Biol. 356 Lab 8. Mortality, Recruitment, and Migration Rates Biol. 356 Lab 8. Moraliy, Recruimen, and Migraion Raes (modified from Cox, 00, General Ecology Lab Manual, McGraw Hill) Las week we esimaed populaion size hrough several mehods. One assumpion of all hese

More information

ON THE NUMBER OF FAMILIES OF BRANCHING PROCESSES WITH IMMIGRATION WITH FAMILY SIZES WITHIN RANDOM INTERVAL

ON THE NUMBER OF FAMILIES OF BRANCHING PROCESSES WITH IMMIGRATION WITH FAMILY SIZES WITHIN RANDOM INTERVAL ON THE NUMBER OF FAMILIES OF BRANCHING PROCESSES ITH IMMIGRATION ITH FAMILY SIZES ITHIN RANDOM INTERVAL Husna Hasan School of Mahemaical Sciences Universii Sains Malaysia, 8 Minden, Pulau Pinang, Malaysia

More information

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes Some common engineering funcions 2.7 Inroducion This secion provides a caalogue of some common funcions ofen used in Science and Engineering. These include polynomials, raional funcions, he modulus funcion

More information

An Introduction to Malliavin calculus and its applications

An Introduction to Malliavin calculus and its applications An Inroducion o Malliavin calculus and is applicaions Lecure 5: Smoohness of he densiy and Hörmander s heorem David Nualar Deparmen of Mahemaics Kansas Universiy Universiy of Wyoming Summer School 214

More information

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales Advances in Dynamical Sysems and Applicaions. ISSN 0973-5321 Volume 1 Number 1 (2006, pp. 103 112 c Research India Publicaions hp://www.ripublicaion.com/adsa.hm The Asympoic Behavior of Nonoscillaory Soluions

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

Lecture Notes 2. The Hilbert Space Approach to Time Series

Lecture Notes 2. The Hilbert Space Approach to Time Series Time Series Seven N. Durlauf Universiy of Wisconsin. Basic ideas Lecure Noes. The Hilber Space Approach o Time Series The Hilber space framework provides a very powerful language for discussing he relaionship

More information

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance

More information

Morning Time: 1 hour 30 minutes Additional materials (enclosed):

Morning Time: 1 hour 30 minutes Additional materials (enclosed): ADVANCED GCE 78/0 MATHEMATICS (MEI) Differenial Equaions THURSDAY JANUARY 008 Morning Time: hour 30 minues Addiional maerials (enclosed): None Addiional maerials (required): Answer Bookle (8 pages) Graph

More information

EXERCISES FOR SECTION 1.5

EXERCISES FOR SECTION 1.5 1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler

More information

Some Basic Information about M-S-D Systems

Some Basic Information about M-S-D Systems Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,

More information

18 Biological models with discrete time

18 Biological models with discrete time 8 Biological models wih discree ime The mos imporan applicaions, however, may be pedagogical. The elegan body of mahemaical heory peraining o linear sysems (Fourier analysis, orhogonal funcions, and so

More information

Echocardiography Project and Finite Fourier Series

Echocardiography Project and Finite Fourier Series Echocardiography Projec and Finie Fourier Series 1 U M An echocardiagram is a plo of how a porion of he hear moves as he funcion of ime over he one or more hearbea cycles If he hearbea repeas iself every

More information

3.6 Derivatives as Rates of Change

3.6 Derivatives as Rates of Change 3.6 Derivaives as Raes of Change Problem 1 John is walking along a sraigh pah. His posiion a he ime >0 is given by s = f(). He sars a =0from his house (f(0) = 0) and he graph of f is given below. (a) Describe

More information

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check

More information

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal

More information

Instructor: Barry McQuarrie Page 1 of 5

Instructor: Barry McQuarrie Page 1 of 5 Procedure for Solving radical equaions 1. Algebraically isolae one radical by iself on one side of equal sign. 2. Raise each side of he equaion o an appropriae power o remove he radical. 3. Simplify. 4.

More information

Cash Flow Valuation Mode Lin Discrete Time

Cash Flow Valuation Mode Lin Discrete Time IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728,p-ISSN: 2319-765X, 6, Issue 6 (May. - Jun. 2013), PP 35-41 Cash Flow Valuaion Mode Lin Discree Time Olayiwola. M. A. and Oni, N. O. Deparmen of Mahemaics

More information

Some Ramsey results for the n-cube

Some Ramsey results for the n-cube Some Ramsey resuls for he n-cube Ron Graham Universiy of California, San Diego Jozsef Solymosi Universiy of Briish Columbia, Vancouver, Canada Absrac In his noe we esablish a Ramsey-ype resul for cerain

More information

Stochastic models and their distributions

Stochastic models and their distributions Sochasic models and heir disribuions Couning cusomers Suppose ha n cusomers arrive a a grocery a imes, say T 1,, T n, each of which akes any real number in he inerval (, ) equally likely The values T 1,,

More information

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t Exercise 7 C P = α + β R P + u C = αp + βr + v (a) (b) C R = α P R + β + w (c) Assumpions abou he disurbances u, v, w : Classical assumions on he disurbance of one of he equaions, eg. on (b): E(v v s P,

More information

non -negative cone Population dynamics motivates the study of linear models whose coefficient matrices are non-negative or positive.

non -negative cone Population dynamics motivates the study of linear models whose coefficient matrices are non-negative or positive. LECTURE 3 Linear/Nonnegaive Marix Models x ( = Px ( A= m m marix, x= m vecor Linear sysems of difference equaions arise in several difference conexs: Linear approximaions (linearizaion Perurbaion analysis

More information

Let us start with a two dimensional case. We consider a vector ( x,

Let us start with a two dimensional case. We consider a vector ( x, Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our

More information

5.1 - Logarithms and Their Properties

5.1 - Logarithms and Their Properties Chaper 5 Logarihmic Funcions 5.1 - Logarihms and Their Properies Suppose ha a populaion grows according o he formula P 10, where P is he colony size a ime, in hours. When will he populaion be 2500? We

More information

ACE 562 Fall Lecture 4: Simple Linear Regression Model: Specification and Estimation. by Professor Scott H. Irwin

ACE 562 Fall Lecture 4: Simple Linear Regression Model: Specification and Estimation. by Professor Scott H. Irwin ACE 56 Fall 005 Lecure 4: Simple Linear Regression Model: Specificaion and Esimaion by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Simple Regression: Economic and Saisical Model

More information

Question 1: Question 2: Topology Exercise Sheet 3

Question 1: Question 2: Topology Exercise Sheet 3 Topology Exercise Shee 3 Prof. Dr. Alessandro Siso Due o 14 March Quesions 1 and 6 are more concepual and should have prioriy. Quesions 4 and 5 admi a relaively shor soluion. Quesion 7 is harder, and you

More information

Final Spring 2007

Final Spring 2007 .615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o

More information

20. Applications of the Genetic-Drift Model

20. Applications of the Genetic-Drift Model 0. Applicaions of he Geneic-Drif Model 1) Deermining he probabiliy of forming any paricular combinaion of genoypes in he nex generaion: Example: If he parenal allele frequencies are p 0 = 0.35 and q 0

More information

Introduction to Probability and Statistics Slides 4 Chapter 4

Introduction to Probability and Statistics Slides 4 Chapter 4 Inroducion o Probabiliy and Saisics Slides 4 Chaper 4 Ammar M. Sarhan, asarhan@mahsa.dal.ca Deparmen of Mahemaics and Saisics, Dalhousie Universiy Fall Semeser 8 Dr. Ammar Sarhan Chaper 4 Coninuous Random

More information

Homework sheet Exercises done during the lecture of March 12, 2014

Homework sheet Exercises done during the lecture of March 12, 2014 EXERCISE SESSION 2A FOR THE COURSE GÉOMÉTRIE EUCLIDIENNE, NON EUCLIDIENNE ET PROJECTIVE MATTEO TOMMASINI Homework shee 3-4 - Exercises done during he lecure of March 2, 204 Exercise 2 Is i rue ha he parameerized

More information

ACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H.

ACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H. ACE 56 Fall 005 Lecure 5: he Simple Linear Regression Model: Sampling Properies of he Leas Squares Esimaors by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Inference in he Simple

More information

Fishing limits and the Logistic Equation. 1

Fishing limits and the Logistic Equation. 1 Fishing limis and he Logisic Equaion. 1 1. The Logisic Equaion. The logisic equaion is an equaion governing populaion growh for populaions in an environmen wih a limied amoun of resources (for insance,

More information

) were both constant and we brought them from under the integral.

) were both constant and we brought them from under the integral. YIELD-PER-RECRUIT (coninued The yield-per-recrui model applies o a cohor, bu we saw in he Age Disribuions lecure ha he properies of a cohor do no apply in general o a collecion of cohors, which is wha

More information

CHAPTER 12 DIRECT CURRENT CIRCUITS

CHAPTER 12 DIRECT CURRENT CIRCUITS CHAPTER 12 DIRECT CURRENT CIUITS DIRECT CURRENT CIUITS 257 12.1 RESISTORS IN SERIES AND IN PARALLEL When wo resisors are conneced ogeher as shown in Figure 12.1 we said ha hey are conneced in series. As

More information

THE GENERALIZED PASCAL MATRIX VIA THE GENERALIZED FIBONACCI MATRIX AND THE GENERALIZED PELL MATRIX

THE GENERALIZED PASCAL MATRIX VIA THE GENERALIZED FIBONACCI MATRIX AND THE GENERALIZED PELL MATRIX J Korean Mah Soc 45 008, No, pp 479 49 THE GENERALIZED PASCAL MATRIX VIA THE GENERALIZED FIBONACCI MATRIX AND THE GENERALIZED PELL MATRIX Gwang-yeon Lee and Seong-Hoon Cho Reprined from he Journal of he

More information

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]

More information

The Strong Law of Large Numbers

The Strong Law of Large Numbers Lecure 9 The Srong Law of Large Numbers Reading: Grimme-Sirzaker 7.2; David Williams Probabiliy wih Maringales 7.2 Furher reading: Grimme-Sirzaker 7.1, 7.3-7.5 Wih he Convergence Theorem (Theorem 54) and

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence MATH 433/533, Fourier Analysis Secion 6, Proof of Fourier s Theorem for Poinwise Convergence Firs, some commens abou inegraing periodic funcions. If g is a periodic funcion, g(x + ) g(x) for all real x,

More information

Predator - Prey Model Trajectories and the nonlinear conservation law

Predator - Prey Model Trajectories and the nonlinear conservation law Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories

More information

( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively:

( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively: XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of

More information

STA 114: Statistics. Notes 2. Statistical Models and the Likelihood Function

STA 114: Statistics. Notes 2. Statistical Models and the Likelihood Function STA 114: Saisics Noes 2. Saisical Models and he Likelihood Funcion Describing Daa & Saisical Models A physicis has a heory ha makes a precise predicion of wha s o be observed in daa. If he daa doesn mach

More information

Lab #2: Kinematics in 1-Dimension

Lab #2: Kinematics in 1-Dimension Reading Assignmen: Chaper 2, Secions 2-1 hrough 2-8 Lab #2: Kinemaics in 1-Dimension Inroducion: The sudy of moion is broken ino wo main areas of sudy kinemaics and dynamics. Kinemaics is he descripion

More information

Oscillation of an Euler Cauchy Dynamic Equation S. Huff, G. Olumolode, N. Pennington, and A. Peterson

Oscillation of an Euler Cauchy Dynamic Equation S. Huff, G. Olumolode, N. Pennington, and A. Peterson PROCEEDINGS OF THE FOURTH INTERNATIONAL CONFERENCE ON DYNAMICAL SYSTEMS AND DIFFERENTIAL EQUATIONS May 4 7, 00, Wilmingon, NC, USA pp 0 Oscillaion of an Euler Cauchy Dynamic Equaion S Huff, G Olumolode,

More information

An Introduction to Backward Stochastic Differential Equations (BSDEs) PIMS Summer School 2016 in Mathematical Finance.

An Introduction to Backward Stochastic Differential Equations (BSDEs) PIMS Summer School 2016 in Mathematical Finance. 1 An Inroducion o Backward Sochasic Differenial Equaions (BSDEs) PIMS Summer School 2016 in Mahemaical Finance June 25, 2016 Chrisoph Frei cfrei@ualbera.ca This inroducion is based on Touzi [14], Bouchard

More information

Welcome Back to Physics 215!

Welcome Back to Physics 215! Welcome Back o Physics 215! (General Physics I) Thurs. Jan 19 h, 2017 Lecure01-2 1 Las ime: Syllabus Unis and dimensional analysis Today: Displacemen, velociy, acceleraion graphs Nex ime: More acceleraion

More information

1 Review of Zero-Sum Games

1 Review of Zero-Sum Games COS 5: heoreical Machine Learning Lecurer: Rob Schapire Lecure #23 Scribe: Eugene Brevdo April 30, 2008 Review of Zero-Sum Games Las ime we inroduced a mahemaical model for wo player zero-sum games. Any

More information

SOME MORE APPLICATIONS OF THE HAHN-BANACH THEOREM

SOME MORE APPLICATIONS OF THE HAHN-BANACH THEOREM SOME MORE APPLICATIONS OF THE HAHN-BANACH THEOREM FRANCISCO JAVIER GARCÍA-PACHECO, DANIELE PUGLISI, AND GUSTI VAN ZYL Absrac We give a new proof of he fac ha equivalen norms on subspaces can be exended

More information

Chapter 7: Solving Trig Equations

Chapter 7: Solving Trig Equations Haberman MTH Secion I: The Trigonomeric Funcions Chaper 7: Solving Trig Equaions Le s sar by solving a couple of equaions ha involve he sine funcion EXAMPLE a: Solve he equaion sin( ) The inverse funcions

More information

Stability and Bifurcation in a Neural Network Model with Two Delays

Stability and Bifurcation in a Neural Network Model with Two Delays Inernaional Mahemaical Forum, Vol. 6, 11, no. 35, 175-1731 Sabiliy and Bifurcaion in a Neural Nework Model wih Two Delays GuangPing Hu and XiaoLing Li School of Mahemaics and Physics, Nanjing Universiy

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and

More information

4.1 - Logarithms and Their Properties

4.1 - Logarithms and Their Properties Chaper 4 Logarihmic Funcions 4.1 - Logarihms and Their Properies Wha is a Logarihm? We define he common logarihm funcion, simply he log funcion, wrien log 10 x log x, as follows: If x is a posiive number,

More information

Lecture 10: The Poincaré Inequality in Euclidean space

Lecture 10: The Poincaré Inequality in Euclidean space Deparmens of Mahemaics Monana Sae Universiy Fall 215 Prof. Kevin Wildrick n inroducion o non-smooh analysis and geomery Lecure 1: The Poincaré Inequaliy in Euclidean space 1. Wha is he Poincaré inequaliy?

More information

arxiv: v1 [math.ca] 15 Nov 2016

arxiv: v1 [math.ca] 15 Nov 2016 arxiv:6.599v [mah.ca] 5 Nov 26 Counerexamples on Jumarie s hree basic fracional calculus formulae for non-differeniable coninuous funcions Cheng-shi Liu Deparmen of Mahemaics Norheas Peroleum Universiy

More information

!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)

!!#$%&#'()!#&'(*%)+,&',-)./0)1-*23) "#"$%&#'()"#&'(*%)+,&',-)./)1-*) #$%&'()*+,&',-.%,/)*+,-&1*#$)()5*6$+$%*,7&*-'-&1*(,-&*6&,7.$%$+*&%'(*8$&',-,%'-&1*(,-&*6&,79*(&,%: ;..,*&1$&$.$%&'()*1$$.,'&',-9*(&,%)?%*,('&5

More information

Seminar 4: Hotelling 2

Seminar 4: Hotelling 2 Seminar 4: Hoelling 2 November 3, 211 1 Exercise Par 1 Iso-elasic demand A non renewable resource of a known sock S can be exraced a zero cos. Demand for he resource is of he form: D(p ) = p ε ε > A a

More information

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4.

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4. PHY1 Elecriciy Topic 7 (Lecures 1 & 11) Elecric Circuis n his opic, we will cover: 1) Elecromoive Force (EMF) ) Series and parallel resisor combinaions 3) Kirchhoff s rules for circuis 4) Time dependence

More information

13.3 Term structure models

13.3 Term structure models 13.3 Term srucure models 13.3.1 Expecaions hypohesis model - Simples "model" a) shor rae b) expecaions o ge oher prices Resul: y () = 1 h +1 δ = φ( δ)+ε +1 f () = E (y +1) (1) =δ + φ( δ) f (3) = E (y +)

More information

15. Vector Valued Functions

15. Vector Valued Functions 1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,

More information

The expectation value of the field operator.

The expectation value of the field operator. The expecaion value of he field operaor. Dan Solomon Universiy of Illinois Chicago, IL dsolom@uic.edu June, 04 Absrac. Much of he mahemaical developmen of quanum field heory has been in suppor of deermining

More information

The average rate of change between two points on a function is d t

The average rate of change between two points on a function is d t SM Dae: Secion: Objecive: The average rae of change beween wo poins on a funcion is d. For example, if he funcion ( ) represens he disance in miles ha a car has raveled afer hours, hen finding he slope

More information

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow 1 KEY Mah 4 Miderm I Fall 8 secions 1 and Insrucor: Sco Glasgow Please do NOT wrie on his eam. No credi will be given for such work. Raher wrie in a blue book, or on our own paper, preferabl engineering

More information

= ( ) ) or a system of differential equations with continuous parametrization (T = R

= ( ) ) or a system of differential equations with continuous parametrization (T = R XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of

More information

arxiv: v1 [math.pr] 19 Feb 2011

arxiv: v1 [math.pr] 19 Feb 2011 A NOTE ON FELLER SEMIGROUPS AND RESOLVENTS VADIM KOSTRYKIN, JÜRGEN POTTHOFF, AND ROBERT SCHRADER ABSTRACT. Various equivalen condiions for a semigroup or a resolven generaed by a Markov process o be of

More information

LECTURE 1: GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS

LECTURE 1: GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS LECTURE : GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS We will work wih a coninuous ime reversible Markov chain X on a finie conneced sae space, wih generaor Lf(x = y q x,yf(y. (Recall ha q

More information

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.

More information

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t M ah 5 2 7 Fall 2 0 0 9 L ecure 1 0 O c. 7, 2 0 0 9 Hamilon- J acobi Equaion: Explici Formulas In his lecure we ry o apply he mehod of characerisics o he Hamilon-Jacobi equaion: u + H D u, x = 0 in R n

More information

INDEPENDENT SETS IN GRAPHS WITH GIVEN MINIMUM DEGREE

INDEPENDENT SETS IN GRAPHS WITH GIVEN MINIMUM DEGREE INDEPENDENT SETS IN GRAPHS WITH GIVEN MINIMUM DEGREE JAMES ALEXANDER, JONATHAN CUTLER, AND TIM MINK Absrac The enumeraion of independen ses in graphs wih various resricions has been a opic of much ineres

More information

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous

More information

Math 334 Fall 2011 Homework 11 Solutions

Math 334 Fall 2011 Homework 11 Solutions Dec. 2, 2 Mah 334 Fall 2 Homework Soluions Basic Problem. Transform he following iniial value problem ino an iniial value problem for a sysem: u + p()u + q() u g(), u() u, u () v. () Soluion. Le v u. Then

More information

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon 3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of

More information

Math 333 Problem Set #2 Solution 14 February 2003

Math 333 Problem Set #2 Solution 14 February 2003 Mah 333 Problem Se #2 Soluion 14 February 2003 A1. Solve he iniial value problem dy dx = x2 + e 3x ; 2y 4 y(0) = 1. Soluion: This is separable; we wrie 2y 4 dy = x 2 + e x dx and inegrae o ge The iniial

More information

dy dx = xey (a) y(0) = 2 (b) y(1) = 2.5 SOLUTION: See next page

dy dx = xey (a) y(0) = 2 (b) y(1) = 2.5 SOLUTION: See next page Assignmen 1 MATH 2270 SOLUTION Please wrie ou complee soluions for each of he following 6 problems (one more will sill be added). You may, of course, consul wih your classmaes, he exbook or oher resources,

More information

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3. Mah Rahman Exam Review Soluions () Consider he IVP: ( 4)y 3y + 4y = ; y(3) = 0, y (3) =. (a) Please deermine he longes inerval for which he IVP is guaraneed o have a unique soluion. Soluion: The disconinuiies

More information

Robust estimation based on the first- and third-moment restrictions of the power transformation model

Robust estimation based on the first- and third-moment restrictions of the power transformation model h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,

More information

SMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15.

SMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15. SMT Calculus Tes Soluions February 5,. Le f() = and le g() =. Compue f ()g (). Answer: 5 Soluion: We noe ha f () = and g () = 6. Then f ()g () =. Plugging in = we ge f ()g () = 6 = 3 5 = 5.. There is a

More information

THE MATRIX-TREE THEOREM

THE MATRIX-TREE THEOREM THE MATRIX-TREE THEOREM 1 The Marix-Tree Theorem. The Marix-Tree Theorem is a formula for he number of spanning rees of a graph in erms of he deerminan of a cerain marix. We begin wih he necessary graph-heoreical

More information

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3 and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or

More information

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan Ground Rules PC11 Fundamenals of Physics I Lecures 3 and 4 Moion in One Dimension A/Prof Tay Seng Chuan 1 Swich off your handphone and pager Swich off your lapop compuer and keep i No alking while lecure

More information

Christos Papadimitriou & Luca Trevisan November 22, 2016

Christos Papadimitriou & Luca Trevisan November 22, 2016 U.C. Bereley CS170: Algorihms Handou LN-11-22 Chrisos Papadimiriou & Luca Trevisan November 22, 2016 Sreaming algorihms In his lecure and he nex one we sudy memory-efficien algorihms ha process a sream

More information

Undetermined coefficients for local fractional differential equations

Undetermined coefficients for local fractional differential equations Available online a www.isr-publicaions.com/jmcs J. Mah. Compuer Sci. 16 (2016), 140 146 Research Aricle Undeermined coefficiens for local fracional differenial equaions Roshdi Khalil a,, Mohammed Al Horani

More information

Displacement ( x) x x x

Displacement ( x) x x x Kinemaics Kinemaics is he branch of mechanics ha describes he moion of objecs wihou necessarily discussing wha causes he moion. 1-Dimensional Kinemaics (or 1- Dimensional moion) refers o moion in a sraigh

More information

Asymptotic Equipartition Property - Seminar 3, part 1

Asymptotic Equipartition Property - Seminar 3, part 1 Asympoic Equipariion Propery - Seminar 3, par 1 Ocober 22, 2013 Problem 1 (Calculaion of ypical se) To clarify he noion of a ypical se A (n) ε and he smalles se of high probabiliy B (n), we will calculae

More information

Math 116 Second Midterm March 21, 2016

Math 116 Second Midterm March 21, 2016 Mah 6 Second Miderm March, 06 UMID: EXAM SOLUTIONS Iniials: Insrucor: Secion:. Do no open his exam unil you are old o do so.. Do no wrie your name anywhere on his exam. 3. This exam has pages including

More information

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is UNIT IMPULSE RESPONSE, UNIT STEP RESPONSE, STABILITY. Uni impulse funcion (Dirac dela funcion, dela funcion) rigorously defined is no sricly a funcion, bu disribuion (or measure), precise reamen requires

More information