Teori, Fourieranalys
By Robin Andersson 2013-2014
Sammanfattning
Innehåller bevisen från teorilistan för kursen i fourieranalys på Chalmers för F2/TM/Kfkb 2014. För eventuella fel kontakta SNF.
Innehåll
1 Konvergenssatsen 1
1.1 Sats . . . 1 1.2 Bevis. . . 1
2 Termvis derivering av fourierserier 1
2.1 Sats . . . 1 2.2 Bevis. . . 2
3 Termvis integrering av fourierserier 2
3.1 Sats . . . 2 3.2 Bevis. . . 2
4 Sats 7.3 om faltning 2
4.1 Sats . . . 2 4.2 Bevis. . . 3
5 Fouriers inversionsformel då f och ˆf ∈ L1 3
5.1 Sats . . . 3 5.2 Bevis. . . 3
6 Plancherels formel f, g, ˆf , ˆg ∈ L1 4
6.1 Sats . . . 4 6.2 Bevis. . . 4 7 Definition av lågpassfilter och formulering av samplingssatsen i termer av låg-
passfilter 4
7.1 Definition . . . 4
8 Sats 3.8 om bästa approximationen 4
8.1 Sats . . . 4 8.2 Bevis. . . 4
10 Definition för ett reguljärt Sturm-Liouville-problem 6 10.1 Definition . . . 6 11 Sats 3.9 a) och b) om Sturm-Liouville-problem: Formulering och bevis 6 11.1 Sats . . . 6 11.2 Bevis. . . 6 12 Den genererande funktionen för Besselsfunktioner (formel 5.20) 6 12.1 Sats . . . 6 12.2 Bevis. . . 7 13 Bevis för Hermitepolynomens ortogonalitetsegenskap 7 13.1 Sats . . . 7 13.2 Bevis. . . 7 14 Härledning av den genererande funktionen för hermitepolynomen 7 14.1 Sats . . . 7 14.2 Bevis. . . 8
R. Andersson Teori, fourieranalysis F2/TM/Kfkb
1 Konvergenssatsen
1.1 Sats
Antag f 2π−periodisk och styckvis glatt i [−π,π]. Då konvergerar SNf (θ) i varje θ och gränsvärdet är f(θ) om f kontinuerlig i θ, annars 12(f (θ−) + f (θ+)).
1.2 Bevis
Inför en ny funktion
g(θ) = f (θ) − f (θ0)
eiθ− eiθ0 , θ 6= θ0.
wlog kan vi konstatera att g också är 2π−periodisk. Påstår att g är styckvis kontinuerlig, utanför θ0 är g styckvis glatt.
Nära θ0gäller
g(θ) = f (θ) − f (θ0) θ − θ0
θ − θ0 eiθ− eiθ0,
där vi får enligt medelvärdessatsen att det vänstra bråket har höger- respektive vänstergränsvärde då θ → θ0, samtidigt gäller
θ − θ0
eiθ− eiθ0 = 1 eiθ− eiθ0
θ − θ0
−→
θ→θ0
1 ieiθ0 .
Alltså är g Riemann-integrabel och korollariet till Bessel ger cn(g) → 0 då n → ±∞.
∴ f (θ) = eiθg(θ) − eiθ0g(θ) + f (θ0) , cn(eiθg(θ)) = 1 2π
Z π
−π
eiθg(θ)e−inθdθ = cn−1(g) .
Vi får alltså
cn(f ) = cn−1(g) − eiθ0cn(g) +
f (θ0) , n = 0 , 0 , n 6= 0 . Partialsumman är därmed enligt
SNf (θ0) =
N
X
n=−N
cn−1(g)einθ0−
N
X
n=−N
cn(g)ei(n+1)θ0+ f (θ0) =
= f (θ0) + c−N −1(g)e−iN θ0− cN(g)ei(N +1)θ0 −→
N →∞f (θ0) .
2 Termvis derivering av fourierserier
2.1 Sats
Antag att f är 2π−periodisk, kontinuerlig och styckvis glatt. Låt an, bn och cn
vara fourierkoefficienterna till f samt a0n b0n och c0n fourierkoefficienterna till f0. Då gäller
a0n= nbn, b0n = −nan, c0n= incn.
1
2.2 Bevis
Vi går efter definitionen och får c0n= 1
2π Z π
−π
f0(θ)e−inθdθ = 1
2πf (θ)e−inθπ
−π− 1 2π
Z π
−π
f (θ)(−ine−inθ) dθ =
eftersom f är 2π−periodisk samt einπ= e−inπ, så följer
= 1
2π · 0 + in 1 2π
Z π
−π
f (θ)e−inθdθ = incn.
Analogt följer proceduren för a0n och b0n.
3 Termvis integrering av fourierserier
3.1 Sats
Antag att f är 2π−periodisk och styckvis kontinuerlig med fourierkoefficienter an, bn och cn, låt F (θ) = R0θf (φ) dφ. Om c0=12a0= 0,
⇒ ∀θ , F (θ) = C0+X
n6=0
cn
ineinθ= 1 2A0+
∞
X
1
an
n sin θn − bn
n cos θn
.
3.2 Bevis
Det gäller att F är styckvis glatt, kontinuerlig och 2π−periodisk, ty, f saknar konstant term. Sats 2.2 gäller för F och ger om F har fourierserien
C0+X
n6=0
Cneinθ ,
att
cn= inCn, n 6= 0 : Cn= cn
in .
4 Sats 7.3 om faltning
4.1 Sats
Antag g ∈ L1(R), R g(x) dx = 1, g begränsad och x3g(x) → 0då x → ±∞. Låt f ∈ L1, om f kontinuerlig i x gäller f ∗ gε(x) → f (x) då ε → 0.
Vidare om f har vänster- respektive högergränsvärde i x och g är jämn gäller
f ∗ gε(x) −→
ε→0
f (x+) + f (x−)
2 .
R. Andersson Teori, fourieranalysis F2/TM/Kfkb
4.2 Bevis
f ∗ gε(x) − f (x) = Z
f (x − εz)g(z) dz − Z
f (x)g(z) dz = Z
(f (x − εz) − f (x))g(z) dz = Z
|z|<ε− 12
(f (x − εz) − f (x))g(z) dz + Z
|z|>ε−1/2
f (x − εz)g(z) dz − f (x) Z
|z|>ε−1/2
g(z) dz =
= I + II + III .
|I| ≤ sup
|y|<√ ε
|f (x − y) − f (x)|
Z
|g(z)| dz −→
ε→00 ,
|II| ≤ Z
|f (x − εz)| dz · sup
|z|>1/√ ε
|g(z)| = 1 ε
Z
|f (x − y)| dy · sup
|z|>1/√ ε
|g(z)| ≤ Z
|f | dy1
εcε32 −→
ε→00 , där c är en konstant. Vidare för den tredje integralen gäller
|III| ≤ |f (x)|
Z
|z|>ε−1/2
|g(z)| dz −→
ε→00
5 Fouriers inversionsformel då f och ˆ f ∈ L
15.1 Sats
Antag att f är integrerbar och styckvis kontinuerlig på R, definierad i sina dis- kontinuitetspunkter på ett sådant sätt att den uppfyller
f (x) =1
2[f (x−) + f (x+)] ∀x ∈R . Då gäller
f (x) = lim
ε→0
1 2π
Z
eiωxe−ε2ω2/2f (ω) dω .ˆ Vidare om ˆf ∈ Ł1, så är f kontinuerlig och
f (x) = 1 2π
Z
eiωxf (ω) dω,ˆ x ∈R .
5.2 Bevis
Eftersom ˆf ∈ Ł1kan vi applicera dominerade konvergenssatsen, alltså 1
2π Z
eiωxf (ω) dω = limˆ
ε→0
1 2π
Z
eiωxe−ε2ω2/2f (ω) dω =ˆ 1
2[f (x−) + f (x+)] .
Integralen i vänsterledet är 1/2π gånger fouriertransformen av ˆf i x−, som påpekat tidigare är fouriertransformen av integrerbara funktioner kontinuerliga. Alltså f är kontinuerlig och
1
2[f (x−) + f (x+)] = f (x) .
3
6 Plancherels formel f, g, ˆ f , ˆ g ∈ L
16.1 Sats
Antag att f, g ∈ Ł1 sådana att ˆf , ˆg ∈ Ł2. Då är också f, g, ˆf , ˆg ∈ Ł2 och vidare gäller 2πhf, gi = h ˆf , ˆgi .
6.2 Bevis
Ty,
2πhf, gi = 2π Z
f (x)g(x) dx = Z Z
f (x)eiωxˆg(ω) dω dx =
= Z Z
f (x)e−iωxg(ω) dω dx =ˆ
Z f (ω)ˆˆ g(ω) dω = h ˆf , ˆgi .
7 Definition av lågpassfilter och formulering av samplings- satsen i termer av lågpassfilter
7.1 Definition
Antag att f(t) är kontinuerlig med Fouriertransform ˆf (ω) = 0 för |ω| ≥ α. Om signalen samplas med frekvensen T1 ≥ απ så kan f(t) återvinnas ur den samplade signalen genom en lågpassfiltrering med avhuggningsfrekvesen α (LPα−filtrering) och multiplikation med T .
8 Sats 3.8 om bästa approximationen
8.1 Sats
Om f ∈ Ł2, är PN1 cn(f )φn den punkt i HN som ligger närmst f, dvs.
min ||f −
N
X
1
dnφn|| = ||f −
N
X
1
cn(f )φn||,
där minimum tas över alla skalärer d1, ..., dN, samt HN är det delrum av Ł2 som spänns upp av vektorerna {φn}N1 .
8.2 Bevis
Om vi väljer dn= cn(f )∀n blir normen i vänsterledet lika med högerledet, så det återstår bara att visa att den aldrig blir mindre än högerledet. Observera att vektorn f −PN
1 cn(f )φn är ortogonal mot varje φk, 1 ≤ k ≤ N , ty
hf −
N
X
1
cn(f )φn, φki = hf, φki − ck(f )hφk, φki = 0 ,
R. Andersson Teori, fourieranalysis F2/TM/Kfkb
det sista enl. def. av ck(f ). Den är därmed också ortogonal mot varje vektor i HN. Vi skriver nu
f −
N
X
1
dnφn= f −
N
X
1
cn(f )φn
! +
N
X
1
(cn(f ) − dn)φn,
och observerar att av de två termerna i högerledet ligger den andra i HN, och den är därför ortogonal mot den första termen. Då ger pythagoras sats
||f −
N
X
1
dnφn||2= ||f −
N
X
1
cn(f )φn||2+
N
X
1
|cn(f ) − dn|2||φn||2.
Eftersom den sista termen här aldrig är negativ, ser vi att
||f −
N
X
1
dnφn||2≥ ||f −
N
X
1
cn(f )φn||2.
Detta säger att normen i vänsterledet i satsens formulering alltid är minst lika stor som högerledet.
9 Satsen om fullständighet för ortogonalsystem
9.1 Sats
Låt {φn}∞1 vara ett ortogormalsystem i Ł2. Följande är ekvivalent:
a) {φn}∞1 är fullständigt.
b)
∞
X
1
|cn(f )|2||φn||2= ||f ||2∀ f ∈ Ł2.
c) om f ∈ Ł2 är ortogonal mot alla φn så gäller f = 0 .
9.2 Bevis
Vi vill visa att a) är ekvivalent med b), att a) implicerar c) och att c) implicerar a).
a) ⇔ b): Följer direkt av sats.
a) implicerar c): Vi utgår från att a) gäller och antar att f ∈ Ł2 är ortogonal mot alla φn. Då gäller att cn(f ) = 0, alltså fourierserienP∞
1 cn(f )φn = 0. Men på grund utav a) är denna serie f , alltså fås f = 0 och c) följer.
c) implicerar a): Tag en godtycklig funktion f ∈ Ł2. Vi vet redan att P∞
1 cn(f )φn kon- vergerar, enligt Bessels olikhet och av ett lemma, så vi kan bilda skillnaden f −P∞
1 cn(f )φn, skalärmultiplicerat med φk fås
hf, φki − ck(f )||φk||2= 0 , enligt definitionen av ck .
där vi utnyttjat skalärproduktens kontinuitet. Skillnaden är alltså ortogonal mot alla vektorer i ortogonalsystemet och därmed 0, ty, vi utgår ifrån c). Därmed är alla tre utsagor ekvivalenta.
5
10 Definition för ett reguljärt Sturm-Liouville-problem
10.1 Definition
Ett reguljärt Sturm-Liouville-problem i intervallet [a, b] ges av
• L f = (rf0)0+ pf, där r, r0 och p är kontinuerliga funktioner i [a, b] samt r > 0 i intervallet.
• Självadjungerade randvillkor.
• En viktfunktion w, kontinuerlig och positiv i [a,b].
11 Sats 3.9 a) och b) om Sturm-Liouville-problem: Formu- lering och bevis
11.1 Sats
Givet ett reguljärt Sturm-Liouville-problem på ett intervall [a, b], då gäller (a) All egenvärden är reella.
(b) Egenfunktioner tillhörande distinkta egenvärden är ortogonala med avse- ende på vikten w; d.v.s. om f och g är egenfunktioner med egenvärde λ och µ : λ 6= µ, så gäller
hf, giw= Z b
a
f (x)g(x)w(x) dx = 0 .
11.2 Bevis
A): Om λ är ett egenvärde med egenfunktion f så gäller
λ||f ||2w= hλwf, f i = −hL(f ), f i = −hf, L(f )i = hf, λwf i = ¯λhf, wf i = ¯λ||f ||2w,
eftersom f uppfyller själv-adjunkta randvillkor. Men eftersom ||f ||2w > 0 inses att ¯λ = λ, alltså λ ∈R .
B): Antag L(f ) + λwf = 0 och L(g) + µwg = 0, där f, g 6= 0. Med samma argument som i a) fås
λhf, giw= hλwf, gi = −hL(f ), gi = −hf, L(g)i = hf, µwgi = µhf, giw. Alltså om λ 6= µ måste vi ha hf, giw= 0.
12 Den genererande funktionen för Besselsfunktioner (for- mel 5.20)
12.1 Sats
∀x och ∀z 6= 0 ,
∞
X
−∞
Jn(x)zn= exp x 2
z −1
z
.
R. Andersson Teori, fourieranalysis F2/TM/Kfkb
12.2 Bevis
Vi vet taylorutvecklingen för exoch erhåller exphxz
2 i
=
∞
X
0
zj j!
x 2
j
, exp −x 2z
=
∞
X
0
(−1)k zkk!
x 2
k .
Vi vet att serierna är absolutkonvergenta, så vi kan multiplicera summorna och den resulterande dubbelsumman blir
exp x 2
z −1
z
=
∞
X
j, k=0
(−1)kzj−k j!k!
x 2
j+k .
Vi summerar serien genom att först addera alla termer som innehåller en given potens av zn och sen summerar vi över n. Vi sätter j − k = n eller j = k + n och erhåller mha
1
(k + n)! = 1
Γ(k + n + 1) = 0 , då k + n < 0 , alltså
exp x 2
z −1
z
=
∞
X
n=−∞
"∞ X
k=0
(−1)k k!(k + n)!
x 2
2k+n# zn=
∞
X
−∞
Jn(x)zn .
13 Bevis för Hermitepolynomens ortogonalitetsegenskap
13.1 Sats
Hermitepolynomen {Hn}∞0 är ortogonala på R med avseende på vikten w(x) = e−x2, och
||Hn||2w= 2nn!√ π .
13.2 Bevis
Om f är ett godtyckligt polynom, så gäller hf, Hniw=
Z ∞
−∞
f (x)Hn(x)e−x2dx = (−1)n Z ∞
−∞
f (x) dn
dxne−x2dx = Z ∞
−∞
f(n)(x)e−x2dx .
I sista ekvationen har vi integrerat partiell n gånger; termerna på randen blir 0 ty P (x)e−x2 → 0 då x → ±∞ ∀ polynom P . Vidare, om f är ett polynom av grad N < n, och speciellt om f = Hm där m < n så gäller f(n) ≡ 0, och alltså hf, Hniw = 0. Detta visar ortogonaliteten av Hermitepolynomen. Men, om f = Hn fås f (x) = (2x)n+ ... och vi får alltså f(n)≡ 2nn!, så
||Hn||2w= 2nn!
Z ∞
−∞
e−x2dx = 2nn!√ π .
14 Härledning av den genererande funktionen för hermite- polynomen
14.1 Sats
∀ x ∈R, ∀ z ∈ C,
∞
X
0
Hn(x)zn
n! = e2xz−z2 .
7
14.2 Bevis
Sätt u = x − z där x är fix, vi har då d/ du = − d/ dz och får alltså dn
dzne−(x−z)2 z=0
= (−1)n dn dune−u2
u=x
= e−u2Hn(u) u=x
= e−x2Hn(x) .
Vidare följer från Taylors formel
e−(x−z)2 =
∞
X
0
e−x2Hn(x)zn n! .