• Keine Ergebnisse gefunden

The Koopman Operator

N/A
N/A
Protected

Academic year: 2021

Aktie "The Koopman Operator"

Copied!
17
0
0

Wird geladen.... (Jetzt Volltext ansehen)

Volltext

(1)

The Koopman Operator

Given(X,B, µ,T), we can take the space of complex-valued square-integrable observables L2(µ). This is a Hilbert space, equipped with inner product

hf,gi= Z

X f(x)·g(x) dµ.

TheKoopman operatoris dened as

UT :L2(µ)→L2(µ), UTf =f ◦T.

Exercise 8.1: Show that the Koopman operator is linear and positive.

(2)

By T -invariance ofµ, it is a unitary operator (it preserves the inner product). Indeed

hUTf,UTgi = Z

X f ◦T(x)·g ◦T(x) dµ

= Z

X(f ·g)◦T(x) dµ= Z

Xf ·g dµ=hf,gi, and therefore UTUT =UTUT =I .

Remark: This has several consequences, common to all unitary operators. The spectrumσ(UT) of UT is a closed subset of the unit circle, and it is in fact a group under multiplication.

(3)

The Transfer Operator

Denition: ThePerron-Frobenius or transfer operator of a transformation T :X →X is the dual of the Koopman operator:

Z

XPTf ·g dµ= Z

X f ·UTg dµ= Z

X f ·g ◦T dµ.

Note that, although UT is independent of the measure, PT is not.

Often it will be important to specify the measure explicitly, and this measure need not be invariant.

Exercise 8.2: Show that the Perron-Frobenius operator has the following properties:

1. PT is linear;

2. PT is positive: f ≥0 implies PTf ≥0.

3. R

PTf dµ=R f dµ. 4. PTk = (PT)k.

(4)

Lemma: Let T : [0,1]→[0,1]be a piecewise C1 interval map.

Then the Perron-Frobenius operator PT w.r.t. Lebesgue measure λ has the pointwise formula

PTf(x) = X

yT1(x)

f(y)

|T0(y)|.

Specically, a xed point of P is an invariant density for T . Proof Let 0=a0<a1 <· · ·<aN =1 be such that T is C1 monotone on each(ai1,ai). Let yi =T1(x)∩(ai1,ai).

(5)

The Transfer Operator

We obtain (PTf)(x) =

Fund. Th. Calc|{z}

d dx

Z x

0 PTf(s) ds= d dx

Z 1

0 (PTf)(s)·1[0,x](s)ds

=

def. P|{z}T

d dx

Z 1

0 f ·1[0,x]◦T(s) ds =

s=|{z}T(u)

d dx

Z

T1[0,x]

f(u) du

= X

T|

(ai1,ai)increasing

xT((ai1,ai))

d dx

Z yi(x) ai1

f(u) du

+ X

T|(ai1,ai)decreasing

xT((ai1,ai))

d dx

Z ai

yi(x)

f(u) du

+ X

T((ai1,ai))⊂[0,x]

d dx

Z ai

ai1

f(u) du

(6)

Continuing this:

(PTf)(x) = X

T|(ai1,ai)increasing

xT((ai1,ai))

f(yi)

T0(yi) + X

T|(ai1,ai)decreasing

xT((ai1,ai))

− f(yi) T0(yi) +0

= X

i

f(yi)

|T0(yi)|, as required.

(7)

The Transfer Operator

There is also a Perron-Frobenius operator with respect to a measureµLebesgue, instead of Lebesgue measure itself:

Lemma: If dµ=hdx, then the operator PTf = PT(f ·h)

h

acts as the Perron-Frobenius operator on(X,B,T, µ). Viewed dierently, if a function h≥0 is xed by PT (w.r.t.

Lebesgue) then dµ=h dx is an invariant measure.

(8)

Proof of the Lemma: Let A beµ-measurable and f ∈L1(µ). Then

Z

APTf dµ = Z

A

PT(f ·h) h h dx

= Z

APT(f ·h)dx

= Z

XPT(f ·h)1Adµ= Z

X f ·h·(1A◦T)dx

= Z

T1Af ·h dx = Z

T1Af dµ.

Because A is arbitrary, the lemma is proved.

(9)

Mixing

Denition: A probability measure preserving dynamical system (X,B, µ,T) ismixing (orstrong mixing) if

µ(Tn(A)∩B)→µ(A)µ(B) as n→ ∞ for every A,B ∈ B.

This says that the events A and B are asymptotically independent.

Exercise: Show that Lebesgue measureµis mixing for the doubling map. In fact, the n-thcorrelation coecient

Corn(A,B) :=µ(Tn(A)∩B)−µ(A)µ(B) =0 for every n≥1.

(10)

Proposition: A probability preserving dynamical system (X,B,T, µ) is mixing if and only if

Z

X f ◦Tn(x)·g(x) dµ→ Z

X f(x) dµ· Z

X g(x) dµas n→ ∞ for all f,g ∈L2(µ), or written in the notation of the Koopman operator UTf =f ◦T :

hUTnf,gi → hf,1ih1,gi as n→ ∞.

Proof: The if-direction follows by taking indicator functions f =1A and g =1B. For the only if-direction, general f,g ∈L2(µ) can be approximated by linear combinations of indicator functions.

(11)

Exponential Mixing

Assume without loss of generality thatR

g dµ=h1,gi=0. Then mixing means that

hUTnf,gi=hf,PTngi →0 as n→ ∞.

If the operator normkPTkrestricted to {g ∈L2(µ) :R

g dµ=0} is strictly less than 1, then

Corn(f,g) =hf,PTngi ≤ kfk kgk kPTkn→0 exponentially fast.

Thusµis exponentially mixing: the events f and g are exponentially independent.

This motivates the spectral analysis of the transfer operator.

(12)

Denition: An operator P is called quasi-compact ifλ=1 is an eigenvalue of multiplicity one, and the rest of the spectrum (i.e., all λ∈Csuch that (P−λI)1 is not a wel-dened bounded

operator), lies in a disk{|λ| ≤σ} for someσ <1.

I Multiplicity one means that the space of eigenfunctions satisfying Pf =λf (forλ=1) is one-dimensional.

I Thisσ is called theessential spectral radius. The actual spectral radius is 1.

I Perron-Frobenius operators are seldom quasi-compact on L2(µ). The (non-trivial) task is to nd a Banach spaceB restricted to which P is quasi-compact.

(13)

Exponential Mixing

TheLasota-Yorke(orDoeblin-Fortet or Tulcea-Ionescu-Marinescu) inequality holds if there are two Banach spaces(Bs,k ks) and (Bw,k kw) (for strong and weak) and σ∈(0,1), L≥1 such that

I Bs ⊂Bw and the unit ball {f ∈Bs :kfks ≤1}is a compact subset of Bw.

I kPfks ≤σkfks+Lkfkw for all f ∈Bs.

I kPfkw ≤ kfkw for all f ∈Bw.

Under these conditions, P:Bs →Bs is quasi-compact, with essential spectral radius≤σ.

Theorem: If the Perron-Frobenius operator for T : [0,1]→[0,1]is quasi-compact with Bs ⊂L1(m)for Lebesgue measure m, then there is an absolutely continuous invariant measureµ which is exponentially mixing.

(14)

Exercise 8.3: Show that kPTfkL1 ≤ kfkL1.

Exercise 8.4: Suppose we have two sequences(an)nNand(bn)nN

satisfying: there areσ ∈(0,1) and L>0 such that:

an+1 ≤an and bn≤σbn+Lan.

Show that(bn)nN is bounded. Show that lim supnbn≤L/(1−σ). The Lasota-Yorke inequality got reinvented several times (Andrej Lasota and James Yorke being the most recent), and they applied this to expanding interval maps with specic choices of Banach spaces.

(15)

Exponential Mixing

Denition: Let g : [a,b]→R. The variation of g is dened to be Var[a,b]g =supXn

i=1

|g(xi)−g(xi1)|,

where the supremum runs over all nite partitions generated by points a=x0 <x1<· · ·<xn=b.

Note that Var is aseminorm, i.e.,Var(f) =06⇒f ≡0. In fact, Var f =Var(f +C) for every constant C. However,

BV :={f ∈L1:Var(f)<∞}with normkfkBV =Var(f) + Z

|f|dx is a Banach space. Its unit ball is compactly embedded in L1.

(16)

Lasota-Yorke Theorem: If T : [0,1]→[0,1]is a piecewise C2 map with infx∈[0,1]|T0(x)|>1, then T preserves an invariant measure µm=Lebesgue, and the density h= ddxµ is of bounded variation. Alsoµis exponentially mixing.

I Compared to the Folklore Theorem, there is an important dierence that the branches need not be onto.

I For instant, β-transformations:

Tβ : [0,1]→[0,1], Tβ(x) =βx mod 1

have an absolutely continuous measures for allβ >1, not just the integers β ≥2.

(17)

Exponential Mixing

The technical details are in the notes, but the key-steps are:

I Take an iterate Tnsuch that infx|(Tn)0(x)|>2, so σ :=2/infx|(Tn)0(x)|<1.

I For Bs =BV and Bw =L1(m) the Lasota-Yorke ineq. holds:

Var[0,1]PT,mg ≤σVar[0,1]g +LkgkL1(m).

I The limit µof Césaro means 1nPn1

k=0m◦Tk is invariant.

I The invariant densityh = ddxµ =limnPTn,m1∈BV.

I The eigenspace of PT,m for λ=1 is spanned by h.

I For all f ∈L1(m)and g ∈BV : Corn(f,g) =

Z

f ·PTn,mg h dx− Z

f dµ Z

g dµ because PTn,mh=h =

Z

f ·PT,m(g − Z

g dµ)h dx

≤ kfkL1kgkBVkPn|hk

≤ kfkL1kgkBVσn.

Referenzen

ÄHNLICHE DOKUMENTE

(Ahiran §ä'a-ra§ul al-muntazar, „Endlich ist der erwartete Mann gekommen&#34;, in: Ahmad H. az-Zaiyät, Wahy ar-Risäla, Bd. 1954 schreibt &#34;Abd an-Näsir in seiner „Philosophie

ences in the distribution of income and the local environmental public good, T CV Y , CV E ,ρ , and dierent correlations of income and endowment with the en- vironmental good at

The introduction of new financial assets increases the income share of financiers in GDP. It also pushes up the share of top wealth people in the country. Although the model is

Presence and being the centre of attraction are two different things and when one discusses about brands, there core purpose is to be the centre of attraction than the competition

8 In combination with the effects of a ⋃ -shaped e-i curve or corresponding effects in a more general equivalence scale specification, this trend has caused a procedure that have

Suppose the education level of the modern sector results in a wage income distribution that is fitted by a gamma pdf with the same shape parameter as that fitted to the wage

Mesa Callejas, Ramon Javier and Rhenals Monterrosa, Remberto and Londoño T., Carlos Mario.. Universidad

Con la salida de la convertibilidad a principios del 2002, y la declaración de default en los mercados internacionales de bonos, un grupo, cada vez más nutrido, de políticos y