08.08.2013 Views

G. Gripenberg (TKK) Mat-1.1531 Gk3-I 20 augusti 2009 1 / 89

G. Gripenberg (TKK) Mat-1.1531 Gk3-I 20 augusti 2009 1 / 89

G. Gripenberg (TKK) Mat-1.1531 Gk3-I 20 augusti 2009 1 / 89

SHOW MORE
SHOW LESS

Transform your PDFs into Flipbooks and boost your revenue!

Leverage SEO-optimized Flipbooks, powerful backlinks, and multimedia content to professionally showcase your products and significantly increase your reach.

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 1 / <strong>89</strong>


Komplexa tal<br />

”In mathematics you don’t understand things. You just get used to them.”<br />

Reell och imaginär del, konjugering, absolutbelopp<br />

z = x + iy = x + y i, x, y ∈ R, i 2 = −1<br />

C är mängden av komplexa tal<br />

John von Neumann<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 2 / <strong>89</strong>


Komplexa tal<br />

”In mathematics you don’t understand things. You just get used to them.”<br />

Reell och imaginär del, konjugering, absolutbelopp<br />

z = x + iy = x + y i, x, y ∈ R, i 2 = −1<br />

C är mängden av komplexa tal<br />

Reell del: Re (x + iy) = x<br />

Imaginär del: Im (x + iy) = y s˚a Im (z) är allts˚a ett reellt tal<br />

John von Neumann<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 2 / <strong>89</strong>


Komplexa tal<br />

”In mathematics you don’t understand things. You just get used to them.”<br />

Reell och imaginär del, konjugering, absolutbelopp<br />

z = x + iy = x + y i, x, y ∈ R, i 2 = −1<br />

C är mängden av komplexa tal<br />

Reell del: Re (x + iy) = x<br />

Imaginär del: Im (x + iy) = y s˚a Im (z) är allts˚a ett reellt tal<br />

Konjugering: x + iy = x − iy (= x + i(−y))<br />

John von Neumann<br />

Absolutbelopp (eller modul) |x + iy| = mod (x + iy) = x 2 + y 2<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 2 / <strong>89</strong>


Komplexa tal<br />

”In mathematics you don’t understand things. You just get used to them.”<br />

Reell och imaginär del, konjugering, absolutbelopp<br />

z = x + iy = x + y i, x, y ∈ R, i 2 = −1<br />

C är mängden av komplexa tal<br />

Reell del: Re (x + iy) = x<br />

Imaginär del: Im (x + iy) = y s˚a Im (z) är allts˚a ett reellt tal<br />

Konjugering: x + iy = x − iy (= x + i(−y))<br />

John von Neumann<br />

Absolutbelopp (eller modul) |x + iy| = mod (x + iy) = x 2 + y 2<br />

Räkneregler<br />

|z| 2 = zz, z1 + z2 = z1 + z2, z1 − z2 = z1 − z2<br />

<br />

z1<br />

z = z, z1z2 = z1 z2, = z2<br />

z1<br />

z2<br />

˛ ˛<br />

|z1 + z2| ≤ |z1| + |z2|, ˛ |z1| − |z2| ˛ ≤ |z1 − z2|<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 2 / <strong>89</strong>


Exempel<br />

Vid addition och subtraktion av komplexa tal adderar och subtraherar man<br />

de reella och imaginära delarna var för sig s˚a att tex.<br />

(8 + 2i) + (−3 − 4i) = (8 + (−3)) + (2 + (−4))i = 5 − 2i.<br />

Vid multiplikation gäller det bara att komma ih˚ag att i 2 = −1:<br />

(8 + 2i)(−3 − 4i) = 8 · (−3) + 8 · (−4)i + 2 · (−3)i + 2 · (−4)i 2<br />

= −24 − 32i − 6i − 8 · (−1) = −16 − 38i.<br />

Division av komplexa tal kan räknas s˚a att man förlänger med nämnarens<br />

konjugat s˚a att man i nämnaren f˚ar ett reellt tal, tex.:<br />

8 + 2i (8 + 2i)(−3 + 4i) −24 + 32i − 6i + 8i2<br />

= =<br />

−3 − 4i (−3 − 4i)(−3 + 4i) (−3) 2 − (4i) 2<br />

= −24 − 8 + 26i<br />

9 − 16i2 = −32 + 26i 26<br />

= −32 +<br />

9 + 16 25 25 i.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 3 / <strong>89</strong>


Kommentar<br />

Ett annat, formellt mera korrekt, sätt att definiera de komplexa talen är<br />

att inte alls (explicit) tala om den imaginära konstanten i utan tala om<br />

punkter (eller vektorer) (x, y) i planet R 2 och definiera räkneoperationer<br />

för dem som motsvarar räkneoperationerna för vanliga reella tal. Addition<br />

är inget problem eftersom det enda förnuftiga är att definiera<br />

(x1, y1) + (x2, y2) = (x1 + x2, y1 + y2),<br />

vilket är addition av vektorer. Ett viktigt villkor som multiplikationen skall<br />

uppfylla är att produkten av tv˚a ”punkter” endast f˚ar vara ”noll” (dvs.<br />

(0, 0)) om ˚atminstone den ena faktorn är ”noll”. Detta uppn˚as om man<br />

definierar<br />

(x1, y1) · (x2, y2) = (x1x2 − y1y2, x1y2 + x2y1)<br />

och man kan d˚a visa att ”alla räkneregler gäller”.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 4 / <strong>89</strong>


Argument eller fasvinkel<br />

Ifall Re (z) = x och Im (z) = y s˚a är argumentet θ = arg(z) av z<br />

⎧ <br />

y<br />

<br />

arctan (+2kπ), x > 0, x + iy<br />

⎪⎨<br />

x<br />

•.<br />

<br />

y<br />

<br />

θ .<br />

θ = arctan + π (+2kπ), x < 0,<br />

x<br />

⎪⎩<br />

y π<br />

(+2kπ), x = 0<br />

|y| 2<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 5 / <strong>89</strong>


Argument eller fasvinkel<br />

Ifall Re (z) = x och Im (z) = y s˚a är argumentet θ = arg(z) av z<br />

⎧ <br />

y<br />

<br />

arctan (+2kπ), x > 0, x + iy<br />

⎪⎨<br />

x<br />

•.<br />

<br />

y<br />

<br />

θ .<br />

θ = arctan + π (+2kπ), x < 0,<br />

x<br />

⎪⎩<br />

y π<br />

(+2kπ), x = 0<br />

|y| 2<br />

atan2<br />

I de flesta programmeringsspr˚ak finns en funktion atan2 som räknar ut det<br />

argument av x + iy som ligger i intervallet (−π, π] med kommandot<br />

atan2(y,x). Observera att i tex. Excel och OOCalc skall man skriva<br />

atan2(x;y) dvs. byta ordning p˚a argumenten.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 5 / <strong>89</strong>


Polär framställning<br />

z = r(cos(θ) + i sin(θ)) = re iθ , r ≥ 0<br />

⇔ |z| = r och arg(z) = θ<br />

⇔ Re (z) = r cos(θ) och Im (z) = r sin(θ)<br />

Kommentar<br />

Om x är ett reellt tal kan man skriva x = |x|sign (x) vilket motsvarar den<br />

polära framställningen z = |z|e iθ med den skillnaden att teckenfunktionen<br />

sign (x) bara f˚ar tv˚a värden (eftersom man inte behöver bry sig om<br />

sign (0)).<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 6 / <strong>89</strong>


Exempel<br />

arg(−3) =?<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 7 / <strong>89</strong>


Exempel<br />

arg(−3) =?<br />

Eftersom den reella delen är negativ är argumentent<br />

) + π (+2kπ) = π (+2kπ).<br />

arctan( 0<br />

−3<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 7 / <strong>89</strong>


Exempel<br />

arg(−3) =?<br />

Eftersom den reella delen är negativ är argumentent<br />

) + π (+2kπ) = π (+2kπ).<br />

arctan( 0<br />

−3<br />

arg(2 − 2i) =?<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 7 / <strong>89</strong>


Exempel<br />

arg(−3) =?<br />

Eftersom den reella delen är negativ är argumentent<br />

) + π (+2kπ) = π (+2kπ).<br />

arctan( 0<br />

−3<br />

arg(2 − 2i) =?<br />

Argumentet är arctan( −2<br />

2<br />

) (+2kπ) = − π<br />

4<br />

(+2kπ).<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 7 / <strong>89</strong>


Exempel<br />

arg(−3) =?<br />

Eftersom den reella delen är negativ är argumentent<br />

) + π (+2kπ) = π (+2kπ).<br />

arctan( 0<br />

−3<br />

arg(2 − 2i) =?<br />

Argumentet är arctan( −2<br />

2<br />

arg(−3e −i 0.1234 ) =?<br />

) (+2kπ) = − π<br />

4<br />

(+2kπ).<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 7 / <strong>89</strong>


Exempel<br />

arg(−3) =?<br />

Eftersom den reella delen är negativ är argumentent<br />

) + π (+2kπ) = π (+2kπ).<br />

arctan( 0<br />

−3<br />

arg(2 − 2i) =?<br />

Argumentet är arctan( −2<br />

2<br />

) (+2kπ) = − π<br />

4<br />

(+2kπ).<br />

arg(−3e −i 0.1234 ) =?<br />

Argument är arg(−3) + arg(e −i 0.1234 ) = π − 0.1234 (+2kπ).<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 7 / <strong>89</strong>


Räkneregler<br />

|z1z2| = |z1||z2|, arg(z1z2) = arg(z1) + arg(z2)<br />

|z n | = |z| n <br />

z1<br />

<br />

, <br />

|z1|<br />

z2<br />

=<br />

|z2|<br />

arg (z n <br />

z1<br />

) = n arg(z), arg = arg(z1) − arg(z2)<br />

z2<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 8 / <strong>89</strong>


Räkneregler<br />

|z1z2| = |z1||z2|, arg(z1z2) = arg(z1) + arg(z2)<br />

|z n | = |z| n <br />

z1<br />

<br />

, <br />

|z1|<br />

z2<br />

=<br />

|z2|<br />

arg (z n <br />

z1<br />

) = n arg(z), arg = arg(z1) − arg(z2)<br />

z1 = z2 ⇔ Re (z1) = Re (z2), Im (z1) = Im (z2)<br />

⇔ |z1| = |z2|, θ1 = θ2 + 2kπ<br />

där θ1 = arg(z1) och θ2 = arg(z2)<br />

z2<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 8 / <strong>89</strong>


Exponentfunktionen<br />

exp(x + iy) = e x+iy = e x cos(y) + i sin(y) <br />

e z1+z2 = e z1 e z2<br />

|e z | = e Re (z) , arg(e z ) = Im (z)<br />

e z = 0, z ∈ C, |e iθ | = 1 θ ∈ R<br />

d’Moivres formel<br />

cos(nt) + i sin(nt) = e int = e it n = cos(t) + i sin(t) n<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 9 / <strong>89</strong>


Logaritmfunktionen z = ln(w) ⇔ w = e z<br />

Om z = x + iy s˚a är |e z | = e x och arg(e z ) = y<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 10 / <strong>89</strong>


Logaritmfunktionen z = ln(w) ⇔ w = e z<br />

Om z = x + iy s˚a är |e z | = e x och arg(e z ) = y<br />

och om w = e z m˚aste |w| = |e z | = e x och<br />

arg(w) + 2kπ = arg(e z ) = y<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 10 / <strong>89</strong>


Logaritmfunktionen z = ln(w) ⇔ w = e z<br />

Om z = x + iy s˚a är |e z | = e x och arg(e z ) = y<br />

och om w = e z m˚aste |w| = |e z | = e x och<br />

arg(w) + 2kπ = arg(e z ) = y<br />

dvs. x = ln(|w|) s˚a att z = ln(w) = ln(|w|) + i(arg(w) + 2kπ).<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 10 / <strong>89</strong>


Logaritmfunktionen z = ln(w) ⇔ w = e z<br />

Om z = x + iy s˚a är |e z | = e x och arg(e z ) = y<br />

och om w = e z m˚aste |w| = |e z | = e x och<br />

arg(w) + 2kπ = arg(e z ) = y<br />

dvs. x = ln(|w|) s˚a att z = ln(w) = ln(|w|) + i(arg(w) + 2kπ).<br />

För att f˚a en ”ordentlig” logaritmfunktion med bara ett värde i varje<br />

punkt kan man tex. definiera Ln(w) = ln(|w|) + iArg(w) där Arg(w)<br />

är argumentet valt s˚a att −π < Arg(z) ≤ π s˚a att tex. ln(|w|)<br />

egentligen är Ln(|w|)<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 10 / <strong>89</strong>


Rötter: z = w 1<br />

n ⇔ w = z n<br />

Om z = |z|e iϕ , dvs. ϕ = arg(z) s˚a är |z n | = |z| n och arg(z n ) = nϕ<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 11 / <strong>89</strong>


Rötter: z = w 1<br />

n ⇔ w = z n<br />

Om z = |z|e iϕ , dvs. ϕ = arg(z) s˚a är |z n | = |z| n och arg(z n ) = nϕ<br />

och om w = z n s˚a är |w| = |z| n och arg(w) + 2kπ = nϕ s˚a att om<br />

arg(w) = θ s˚a är<br />

|z| = |w| 1<br />

n och ϕ = θ 2kπ<br />

n + n<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 11 / <strong>89</strong>


Rötter: z = w 1<br />

n ⇔ w = z n<br />

Om z = |z|e iϕ , dvs. ϕ = arg(z) s˚a är |z n | = |z| n och arg(z n ) = nϕ<br />

och om w = z n s˚a är |w| = |z| n och arg(w) + 2kπ = nϕ s˚a att om<br />

arg(w) = θ s˚a är<br />

|z| = |w| 1<br />

n och ϕ = θ<br />

n<br />

z = w 1<br />

n = n√ w = n |w|<br />

2kπ + n dvs.<br />

<br />

θ+2kπ<br />

θ+2kπ<br />

cos + i sin<br />

<br />

,<br />

n<br />

där k = 0, 1, . . . , n − 1 eftersom man f˚ar samma värden för k + n som<br />

för k.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 11 / <strong>89</strong><br />

n


Exempel<br />

L˚at w = 1 + i . Bestäm den lösning till ekvationen z 4 = w, vars argument<br />

ligger i intervallet [π, 3<br />

2 π].<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 12 / <strong>89</strong>


Exempel<br />

L˚at w = 1 + i . Bestäm den lösning till ekvationen z4 = w, vars argument<br />

ligger i intervallet [π, 3<br />

2π]. Lösning: Absolutbeloppet av talet w är |w| = √ 12 + 12 = √ 2 ≈ 1.4142,<br />

och w:s argument är arctan( 1 π<br />

1 ) = 4 .<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 12 / <strong>89</strong>


Exempel<br />

L˚at w = 1 + i . Bestäm den lösning till ekvationen z4 = w, vars argument<br />

ligger i intervallet [π, 3<br />

2π]. Lösning: Absolutbeloppet av talet w är |w| = √ 12 + 12 = √ 2 ≈ 1.4142,<br />

och w:s argument är arctan( 1 π<br />

1 ) = 4 . Ifall |z| = r och arg(z) = ϕ, s˚a är<br />

|z4 | = r 4 och arg(z 4 ) = 4ϕ. Om nu z4 = w s˚a är r 4 = |w| = √ 2 och<br />

+ 2kπ där k är ett heltal.<br />

4ϕ = π<br />

4<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 12 / <strong>89</strong>


Exempel<br />

L˚at w = 1 + i . Bestäm den lösning till ekvationen z4 = w, vars argument<br />

ligger i intervallet [π, 3<br />

2π]. Lösning: Absolutbeloppet av talet w är |w| = √ 12 + 12 = √ 2 ≈ 1.4142,<br />

och w:s argument är arctan( 1 π<br />

1 ) = 4 . Ifall |z| = r och arg(z) = ϕ, s˚a är<br />

|z4 | = r 4 och arg(z 4 ) = 4ϕ. Om nu z4 = w s˚a är r 4 = |w| = √ 2 och<br />

4ϕ = π<br />

4 + 2kπ där k är ett heltal. Av detta följer att r = 8√ 2 ≈ 1.0905 och<br />

ϕ = π π<br />

16 + 2 k = 0.19635 + 1.5708k. Nu f˚ar man olika lösningar d˚a<br />

k = 0, 1, . . . , 3 eftersom man d˚a tex. k = 4 f˚ar samma tal som d˚a k = 0<br />

osv.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 12 / <strong>89</strong>


Exempel<br />

L˚at w = 1 + i . Bestäm den lösning till ekvationen z4 = w, vars argument<br />

ligger i intervallet [π, 3<br />

2π]. Lösning: Absolutbeloppet av talet w är |w| = √ 12 + 12 = √ 2 ≈ 1.4142,<br />

och w:s argument är arctan( 1 π<br />

1 ) = 4 . Ifall |z| = r och arg(z) = ϕ, s˚a är<br />

|z4 | = r 4 och arg(z 4 ) = 4ϕ. Om nu z4 = w s˚a är r 4 = |w| = √ 2 och<br />

4ϕ = π<br />

4 + 2kπ där k är ett heltal. Av detta följer att r = 8√ 2 ≈ 1.0905 och<br />

ϕ = π π<br />

16 + 2 k = 0.19635 + 1.5708k. Nu f˚ar man olika lösningar d˚a<br />

k = 0, 1, . . . , 3 eftersom man d˚a tex. k = 4 f˚ar samma tal som d˚a k = 0<br />

osv. Eftersom argumenten för dehär lösningarna är π π π π<br />

16 , 16 + 2 , 16 + π och<br />

π 3π<br />

3<br />

16 + 2 s˚a ser vi att den lösning vars argument ligger i intervallet [π, 2π] f˚as d˚a k = 2 och är allts˚a<br />

<br />

<br />

z2 = 1.0905 cos(0.19635 + 1.5708 · 2) + i sin(0.19635 + 1.5708 · 2)<br />

= −1.0696 − i 0.21275.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 12 / <strong>89</strong>


Trigonometriska funktioner<br />

e iz − e −iz <br />

sin(z) = 1<br />

2i<br />

cos(z) = 1<br />

<br />

2 eiz + e−iz <br />

sin(x + iy) = sin(x) cosh(y) + i cos(x) sinh(y)<br />

cos(x + iy) = cos(x) cosh(y) − i sin(x) sinh(y)<br />

Möbius-avbildning<br />

z ↦→ az+b<br />

cz+d , ad − bc = 0<br />

Cirklar och linjer avbildas p˚a cirklar eller linjer<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 13 / <strong>89</strong>


Öppna och slutna mängder<br />

En mängd Ω ⊂ C är öppen ifall<br />

för varje z0 ∈ Ω finns ett tal δ > 0 s˚a att<br />

{ z ∈ C : |z − z0| < δ } ⊂ Ω<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 14 / <strong>89</strong>


Öppna och slutna mängder<br />

En mängd Ω ⊂ C är öppen ifall<br />

för varje z0 ∈ Ω finns ett tal δ > 0 s˚a att<br />

{ z ∈ C : |z − z0| < δ } ⊂ Ω<br />

En mängd A ⊂ C är sluten ifall<br />

C \ A = { z ∈ C : z /∈ A } är öppen<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 14 / <strong>89</strong>


Öppna och slutna mängder<br />

En mängd Ω ⊂ C är öppen ifall<br />

för varje z0 ∈ Ω finns ett tal δ > 0 s˚a att<br />

{ z ∈ C : |z − z0| < δ } ⊂ Ω<br />

En mängd A ⊂ C är sluten ifall<br />

C \ A = { z ∈ C : z /∈ A } är öppen<br />

En mängd Ω ⊂ C är öppen när<br />

den inte inneh˚aller n˚agon randpunkt, dvs. ∂Ω ∩ Ω = ∅<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 14 / <strong>89</strong>


Öppna och slutna mängder<br />

En mängd Ω ⊂ C är öppen ifall<br />

för varje z0 ∈ Ω finns ett tal δ > 0 s˚a att<br />

{ z ∈ C : |z − z0| < δ } ⊂ Ω<br />

En mängd A ⊂ C är sluten ifall<br />

C \ A = { z ∈ C : z /∈ A } är öppen<br />

En mängd Ω ⊂ C är öppen när<br />

den inte inneh˚aller n˚agon randpunkt, dvs. ∂Ω ∩ Ω = ∅<br />

där randen ∂Ω best˚ar av de punkter z ∈ C för vilka det<br />

för varje δ > 0 finns en punkt zi ∈ Ω och en punkt zu ∈ C \ Ω<br />

s˚a att |zi − z| < δ och |zu − z| < δ<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 14 / <strong>89</strong>


Öppna och slutna mängder<br />

En mängd Ω ⊂ C är öppen ifall<br />

för varje z0 ∈ Ω finns ett tal δ > 0 s˚a att<br />

{ z ∈ C : |z − z0| < δ } ⊂ Ω<br />

En mängd A ⊂ C är sluten ifall<br />

C \ A = { z ∈ C : z /∈ A } är öppen<br />

En mängd Ω ⊂ C är öppen när<br />

den inte inneh˚aller n˚agon randpunkt, dvs. ∂Ω ∩ Ω = ∅<br />

där randen ∂Ω best˚ar av de punkter z ∈ C för vilka det<br />

för varje δ > 0 finns en punkt zi ∈ Ω och en punkt zu ∈ C \ Ω<br />

s˚a att |zi − z| < δ och |zu − z| < δ<br />

En mängd Ω ⊂ C är sluten när<br />

den inneh˚aller alla randpunkter, dvs. ∂Ω ⊂ Ω<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 14 / <strong>89</strong>


Sammanhängande mängder<br />

En öppen mängd Ω ⊂ C är sammanhängande ifall det<br />

för alla z0, z1 ∈ Ω finns en kontinuerlig funktion<br />

γ : [0, 1] → Ω s˚a att γ(0) = z0 och γ(1) = z1<br />

(och allts˚a γ(t) ∈ Ω, t ∈ [0, 1], dvs. en kurva i Ω fr˚an z0 till z1).<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 15 / <strong>89</strong>


Sammanhängande mängder<br />

En öppen mängd Ω ⊂ C är sammanhängande ifall det<br />

för alla z0, z1 ∈ Ω finns en kontinuerlig funktion<br />

γ : [0, 1] → Ω s˚a att γ(0) = z0 och γ(1) = z1<br />

(och allts˚a γ(t) ∈ Ω, t ∈ [0, 1], dvs. en kurva i Ω fr˚an z0 till z1).<br />

Mera allmänt: A ⊂ C är sammanhängande om följande gäller (som<br />

allts˚a är ekvivalent med ovanst˚aende villkor d˚a A är öppen):<br />

Ifall A ⊂ Ω1 ∪ Ω2 där Ωj ⊂ C, j = 1, 2 är öppen och<br />

Ω1 ∩ Ω2 = ∅ s˚a är A ∩ Ω1 = ∅ eller A ∩ Ω2 = ∅<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 15 / <strong>89</strong>


Sammanhängande mängder<br />

En öppen mängd Ω ⊂ C är sammanhängande ifall det<br />

för alla z0, z1 ∈ Ω finns en kontinuerlig funktion<br />

γ : [0, 1] → Ω s˚a att γ(0) = z0 och γ(1) = z1<br />

(och allts˚a γ(t) ∈ Ω, t ∈ [0, 1], dvs. en kurva i Ω fr˚an z0 till z1).<br />

Mera allmänt: A ⊂ C är sammanhängande om följande gäller (som<br />

allts˚a är ekvivalent med ovanst˚aende villkor d˚a A är öppen):<br />

Ifall A ⊂ Ω1 ∪ Ω2 där Ωj ⊂ C, j = 1, 2 är öppen och<br />

Ω1 ∩ Ω2 = ∅ s˚a är A ∩ Ω1 = ∅ eller A ∩ Ω2 = ∅<br />

Kontinuerliga funktioner<br />

Om Ω ⊂ C s˚a är funktionen f : Ω → C kontinuerlig i Ω om<br />

lim<br />

z→z0<br />

z∈Ω<br />

f (z) = f (z0) z0 ∈ Ω.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 15 / <strong>89</strong>


Derivator<br />

En funktion f är deriverbar i punkten z0 ifall<br />

limz→z0<br />

f (z)−f (z0)<br />

z−z0<br />

= f ′ (z0)<br />

för n˚agot komplext tal f ′ (z0) dvs.,<br />

för varje ɛ > 0 finns det ett tal δ > 0 s˚a att om 0 < |z − z0| < δ s˚a är<br />

<br />

f definierad för z och z0 och <br />

f (z) − f (z0)<br />

<br />

− f<br />

z − z0<br />

′ <br />

<br />

(z0) <br />

< ɛ<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 16 / <strong>89</strong>


Derivator<br />

En funktion f är deriverbar i punkten z0 ifall<br />

limz→z0<br />

f (z)−f (z0)<br />

z−z0<br />

= f ′ (z0)<br />

för n˚agot komplext tal f ′ (z0) dvs.,<br />

för varje ɛ > 0 finns det ett tal δ > 0 s˚a att om 0 < |z − z0| < δ s˚a är<br />

<br />

f definierad för z och z0 och <br />

f (z) − f (z0)<br />

<br />

− f<br />

z − z0<br />

′ <br />

<br />

(z0) <br />

< ɛ<br />

Alla normala deriveringsregler gäller och tex.<br />

d<br />

dz zm = mz m−1 ( d˚a m inte är ett heltal)<br />

d<br />

dz ez = e z<br />

d<br />

1<br />

dz ln(z) = z<br />

d<br />

d<br />

dz sin(z) = cos(z)<br />

d<br />

dz gf (z) = g ′f (z) f ′ (z)<br />

<br />

d<br />

dz f (z)g(z) = f ′ (z)g(z) + f (z)g ′ (z)<br />

dz cos(z) = − sin(z)<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 16 / <strong>89</strong>


Analytiska funktioner<br />

En funktion f är analytisk i mängden A ⊂ C ifall<br />

det finns en öppen mängd Ω ⊂ C s˚a att A ⊂ Ω och<br />

f är deriverbar i varje punkt i Ω<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 17 / <strong>89</strong>


Analytiska funktioner<br />

En funktion f är analytisk i mängden A ⊂ C ifall<br />

det finns en öppen mängd Ω ⊂ C s˚a att A ⊂ Ω och<br />

f är deriverbar i varje punkt i Ω<br />

Cauchy-Riemann ekvationerna<br />

Ifall f (z) = u(x, y) + iv(x, y) d˚a z = x + iy s˚a gäller<br />

ux(x, y) = vy (x, y)<br />

uy (x, y) = −vx(x, y)<br />

i de punkter där f är deriverbar<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 17 / <strong>89</strong>


Stig<br />

En stig, eller bitvis differentierbar glatt kurva är en kontinuerlig funktion<br />

γ : [a, b] → C s˚a att det finns ändligt m˚anga punkter<br />

a = t0 < t1 < . . . < tn = b s˚a att γ är kontinuerligt deriverbar i varje<br />

intervall [tj−1, tj] och derivatan (höger- eller vänster- i ändpunkterna) är<br />

aldrig 0.<br />

γ∗ def<br />

= { γ(t) : t ∈ [a, b] }<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 18 / <strong>89</strong>


Stig<br />

En stig, eller bitvis differentierbar glatt kurva är en kontinuerlig funktion<br />

γ : [a, b] → C s˚a att det finns ändligt m˚anga punkter<br />

a = t0 < t1 < . . . < tn = b s˚a att γ är kontinuerligt deriverbar i varje<br />

intervall [tj−1, tj] och derivatan (höger- eller vänster- i ändpunkterna) är<br />

aldrig 0.<br />

γ∗ def<br />

= { γ(t) : t ∈ [a, b] }<br />

Stigintegraler<br />

Ifall γ är en stig är<br />

<br />

γ<br />

f (z) dz =<br />

b<br />

a<br />

f (γ(t))γ ′ (t) dt<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 18 / <strong>89</strong>


Stig<br />

En stig, eller bitvis differentierbar glatt kurva är en kontinuerlig funktion<br />

γ : [a, b] → C s˚a att det finns ändligt m˚anga punkter<br />

a = t0 < t1 < . . . < tn = b s˚a att γ är kontinuerligt deriverbar i varje<br />

intervall [tj−1, tj] och derivatan (höger- eller vänster- i ändpunkterna) är<br />

aldrig 0.<br />

γ∗ def<br />

= { γ(t) : t ∈ [a, b] }<br />

Stigintegraler<br />

Ifall γ är en stig är<br />

<br />

f (z) dz =<br />

b<br />

f (γ(t))γ ′ (t) dt<br />

γ<br />

a<br />

Ofta skriver man <br />

<br />

f (z) dz istället för f (z) dz där C är en ”riktad kurva” om<br />

C γ<br />

det är klart hur parameterframställningen γ för C (som allts˚a är γ∗ ”med<br />

riktning”) skall väljas.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 18 / <strong>89</strong>


Exempel p˚a integraler<br />

Om C är ”sträckan fr˚an z0 till z1” s˚a är<br />

<br />

<br />

f (z) dz = f (z) dz =<br />

C<br />

γ<br />

1<br />

0<br />

f (γ(t))γ ′ (t) dt,<br />

där γ(t) = (1 − t)z0 + tz1, t ∈ [0, 1]. Alternativt,<br />

<br />

C f (z) dz = b<br />

a f (γ(t))γ′ (t) dt där γ(t) = b−t<br />

b−a z0 + t−a<br />

b−a z1<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 19 / <strong>89</strong>


Exempel p˚a integraler<br />

Om C är ”sträckan fr˚an z0 till z1” s˚a är<br />

<br />

<br />

f (z) dz = f (z) dz =<br />

C<br />

γ<br />

1<br />

0<br />

f (γ(t))γ ′ (t) dt,<br />

där γ(t) = (1 − t)z0 + tz1, t ∈ [0, 1]. Alternativt,<br />

<br />

C f (z) dz = b<br />

a f (γ(t))γ′ (t) dt där γ(t) = b−t<br />

b−a z0 + t−a<br />

b−a z1<br />

Om C är ”cirkelb˚agen |z − z0| = r fr˚an z0 + reiα till z0 + reiβ motsols” (där α < β) s˚a är<br />

<br />

C<br />

<br />

f (z) dz =<br />

där γ(t) = z0 + re it , t ∈ [α, β].<br />

γ<br />

β<br />

f (z) dz = f (γ(t))γ ′ (t) dt,<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 19 / <strong>89</strong><br />

α


Lemma<br />

Tv˚a stigar γ1 : [a, b] → C och γ2 : [c, d] → C där γ1(b) = γ2(c) kan<br />

kombineras till en stig γ = γ1 ⊕ γ2, γ : [a, b + d − c] → C, s˚a att<br />

(<br />

γ1(t), t ∈ [a, b],<br />

γ(t) =<br />

γ2(c + t − b), t ∈ [b, b + d − c],<br />

och d˚a är <br />

γ<br />

<br />

f (z) dz =<br />

γ1<br />

<br />

f (z) dz +<br />

(P˚a motsvarande sätt kan en stig delas upp.)<br />

γ2<br />

f (z) dz.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 <strong>20</strong> / <strong>89</strong>


Lemma<br />

Tv˚a stigar γ1 : [a, b] → C och γ2 : [c, d] → C där γ1(b) = γ2(c) kan<br />

kombineras till en stig γ = γ1 ⊕ γ2, γ : [a, b + d − c] → C, s˚a att<br />

(<br />

γ1(t), t ∈ [a, b],<br />

γ(t) =<br />

γ2(c + t − b), t ∈ [b, b + d − c],<br />

och d˚a är <br />

γ<br />

<br />

f (z) dz =<br />

γ1<br />

<br />

f (z) dz +<br />

(P˚a motsvarande sätt kan en stig delas upp.)<br />

Lemma<br />

γ2<br />

f (z) dz.<br />

Om γ : [a, b] → C är en stig och −γ eller γ↢ är stigen<br />

γ↢ (t) = γ(a + b − t), t ∈ [a, b], dvs. stigen γ ”i omvänd riktning”, s˚a är<br />

<br />

<br />

f (z) dz = − f (z) dz.<br />

γ ↢<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 <strong>20</strong> / <strong>89</strong><br />

γ


Slutna stigar<br />

En stig γ : [a, b] → C är sluten ifall γ(a) = γ(b) och <br />

att γ är sluten.<br />

γ<br />

f (z) dz betyder<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 21 / <strong>89</strong>


Slutna stigar<br />

En stig γ : [a, b] → C är sluten ifall γ(a) = γ(b) och <br />

att γ är sluten.<br />

Lemma<br />

Ifall |f (z)| ≤ M d˚a z ∈ γ∗ , (a < b) och<br />

längden av γ är ≤ L s˚a gäller<br />

<br />

<br />

<br />

<br />

f (z) dz<br />

≤<br />

b<br />

|f (γ(t))||γ ′ (t)| dt ≤ ML.<br />

γ<br />

a<br />

γ<br />

f (z) dz betyder<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 21 / <strong>89</strong>


Slutna stigar<br />

En stig γ : [a, b] → C är sluten ifall γ(a) = γ(b) och <br />

att γ är sluten.<br />

Lemma<br />

Ifall |f (z)| ≤ M d˚a z ∈ γ∗ , (a < b) och<br />

längden av γ är ≤ L s˚a gäller<br />

<br />

<br />

<br />

<br />

f (z) dz<br />

≤<br />

b<br />

|f (γ(t))||γ ′ (t)| dt ≤ ML.<br />

Lemma<br />

γ<br />

a<br />

γ<br />

f (z) dz betyder<br />

Om γ : [a, b] → C är en stig och f är en kontinuerlig funktion s˚a att det<br />

finns en deriverbar funktion F s˚a att F ′ (z) = f (z) d˚a z ∈ γ ∗ s˚a är<br />

<br />

γ<br />

f (z) dz = F (γ(b)) − F (γ(a)).<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 21 / <strong>89</strong>


Vridningstal<br />

Ifall γ är en sluten stig s˚a är<br />

ν(γ, w) = 1<br />

<br />

1<br />

2πi γ z−w dz<br />

stigens vridningstal i förh˚allande till w /∈ γ∗ .<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 22 / <strong>89</strong>


Vridningstal<br />

Ifall γ är en sluten stig s˚a är<br />

ν(γ, w) = 1<br />

<br />

1<br />

2πi γ z−w dz<br />

stigens vridningstal i förh˚allande till w /∈ γ∗ .<br />

Lemma<br />

ν(γ, w) ett heltal som anger hur m˚anga varv γ g˚ar runt w i positiv<br />

riktning och är konstant i varje öppen sammanhängande delmängd av<br />

C \ γ ∗ .<br />

ν(γ ↢ , w) = −ν(γ, w), w /∈ γ ∗<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 22 / <strong>89</strong>


Cauchys integralteorem<br />

Om γ är en sluten stig och f är analytisk p˚a γ∗ och i alla punkter<br />

”innanför” γ∗ <br />

s˚a är<br />

f (z) dz = 0<br />

γ<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 23 / <strong>89</strong>


Cauchys integralteorem<br />

Om γ är en sluten stig och f är analytisk p˚a γ∗ och i alla punkter<br />

”innanför” γ∗ <br />

s˚a är<br />

f (z) dz = 0<br />

γ<br />

Antagandet mera exakt: γ är en sluten stig, f är analytisk i en öppen<br />

mängd Ω s˚a att γ ∗ ⊂ Ω och ν(γ, z) = 0 för alla z ∈ C \ Ω<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 23 / <strong>89</strong>


Exempel<br />

Visa hur man kan räkna ut integralen ∞<br />

0<br />

integralteorem.<br />

sin(t)<br />

t dt med hjälp av Cauchys<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 24 / <strong>89</strong>


Exempel<br />

Visa hur man kan räkna ut integralen ∞ sin(t)<br />

0 t dt med hjälp av Cauchys<br />

integralteorem.<br />

Lösning: Vi bildar en sluten stig γ som är sammansatt av stigarn γS,<br />

γ [−S,−r], γr och av γ [r,S] där 0 < r < S och γS är den stig som best˚ar av<br />

sträckorna fr˚an S till S + i √ S, fr˚an S + i √ S till −S + i √ S och fr˚an<br />

−S + i √ S till −S, γ [−S,−r] är sträckan fr˚an −S + 0i till −r + 0i, γr är<br />

cirkelb˚agen |z| = r i negativ riktning fr˚an −r to r och γ [r,S] är sträckan<br />

fr˚an r + 0i till S + 0i.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 24 / <strong>89</strong>


Exempel<br />

Visa hur man kan räkna ut integralen ∞ sin(t)<br />

0 t dt med hjälp av Cauchys<br />

integralteorem.<br />

Lösning: Vi bildar en sluten stig γ som är sammansatt av stigarn γS,<br />

γ [−S,−r], γr och av γ [r,S] där 0 < r < S och γS är den stig som best˚ar av<br />

sträckorna fr˚an S till S + i √ S, fr˚an S + i √ S till −S + i √ S och fr˚an<br />

−S + i √ S till −S, γ [−S,−r] är sträckan fr˚an −S + 0i till −r + 0i, γr är<br />

cirkelb˚agen |z| = r i negativ riktning fr˚an −r to r och γ [r,S] är sträckan<br />

fr˚an r + 0i till S + 0i. Av Cauchys integralteorem följer att <br />

dvs. <br />

γS<br />

e iz<br />

z<br />

<br />

e dz + γ [−S,−r]<br />

iz <br />

z dz + γr<br />

e iz<br />

z<br />

<br />

e dz + γ [r,S]<br />

iz<br />

z dz = 0<br />

γ<br />

e iz<br />

z<br />

dz = 0<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 24 / <strong>89</strong>


Exempel<br />

Visa hur man kan räkna ut integralen ∞ sin(t)<br />

0 t dt med hjälp av Cauchys<br />

integralteorem.<br />

Lösning: Vi bildar en sluten stig γ som är sammansatt av stigarn γS,<br />

γ [−S,−r], γr och av γ [r,S] där 0 < r < S och γS är den stig som best˚ar av<br />

sträckorna fr˚an S till S + i √ S, fr˚an S + i √ S till −S + i √ S och fr˚an<br />

−S + i √ S till −S, γ [−S,−r] är sträckan fr˚an −S + 0i till −r + 0i, γr är<br />

cirkelb˚agen |z| = r i negativ riktning fr˚an −r to r och γ [r,S] är sträckan<br />

fr˚an r + 0i till S + 0i. Av Cauchys integralteorem följer att <br />

e<br />

γ<br />

iz<br />

z dz = 0<br />

dvs. <br />

e<br />

γS<br />

iz <br />

e<br />

z dz + γ [−S,−r]<br />

iz <br />

e<br />

z dz + γr<br />

iz <br />

e<br />

z dz + γ [r,S]<br />

iz<br />

z dz = 0 s˚a det gäller<br />

att visa att ∞ sin(t)<br />

1<br />

0 t dt = lim r→0 2<br />

S→∞<br />

Im<br />

e<br />

γ [−S,−r]<br />

iz <br />

e<br />

z dz + γ [r,S]<br />

iz<br />

z dz<br />

<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 24 / <strong>89</strong>


Exempel<br />

Visa hur man kan räkna ut integralen ∞ sin(t)<br />

0 t dt med hjälp av Cauchys<br />

integralteorem.<br />

Lösning: Vi bildar en sluten stig γ som är sammansatt av stigarn γS,<br />

γ [−S,−r], γr och av γ [r,S] där 0 < r < S och γS är den stig som best˚ar av<br />

sträckorna fr˚an S till S + i √ S, fr˚an S + i √ S till −S + i √ S och fr˚an<br />

−S + i √ S till −S, γ [−S,−r] är sträckan fr˚an −S + 0i till −r + 0i, γr är<br />

cirkelb˚agen |z| = r i negativ riktning fr˚an −r to r och γ [r,S] är sträckan<br />

fr˚an r + 0i till S + 0i. Av Cauchys integralteorem följer att <br />

e<br />

γ<br />

iz<br />

z dz = 0<br />

dvs. <br />

e<br />

γS<br />

iz <br />

e<br />

z dz + γ [−S,−r]<br />

iz <br />

e<br />

z dz + γr<br />

iz <br />

e<br />

z dz + γ [r,S]<br />

iz<br />

z dz = 0 s˚a det gäller<br />

att visa att ∞ sin(t)<br />

1<br />

0 t dt = lim r→0 2<br />

S→∞<br />

Im<br />

e<br />

γ [−S,−r]<br />

iz <br />

e<br />

z dz + γ [r,S]<br />

iz<br />

z dz<br />

<br />

och<br />

<br />

<br />

sedan räkna ut limS→∞ dz och limr→0<br />

γS<br />

e iz<br />

z<br />

γr<br />

e iz<br />

z dz.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 24 / <strong>89</strong>


Cauchys integralteorem, ver. 1.1<br />

Om γ1 och γ2 är slutna stigar och f är analytisk p˚a γ∗ 1 och γ∗ 2<br />

punkter ”mellan” γ∗ 1 och γ∗ 2 s˚a är<br />

<br />

γ1<br />

<br />

f (z) dz =<br />

γ2<br />

f (z) dz.<br />

och i alla<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 25 / <strong>89</strong>


Cauchys integralteorem, ver. 1.1<br />

Om γ1 och γ2 är slutna stigar och f är analytisk p˚a γ∗ 1 och γ∗ 2<br />

punkter ”mellan” γ∗ 1 och γ∗ 2 s˚a är<br />

<br />

γ1<br />

<br />

f (z) dz =<br />

γ2<br />

f (z) dz.<br />

och i alla<br />

Antagandet mera exakt: γ1 och γ2 är slutna stigar, f är analytisk i en<br />

öppen mängd Ω s˚a att γ ∗ 1 ⊂ Ω, γ∗ 2 ⊂ Ω och ν(γ1, z) = ν(γ2, z) för alla<br />

z ∈ C \ Ω.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 25 / <strong>89</strong>


Cauchys integralformel, ver 1.0<br />

Om γ är en sluten stig som g˚ar ett varv runt punkten w och f är analytisk<br />

p˚a γ ∗ och i alla punkter ”innanför” γ ∗ s˚a är<br />

f (w) = 1<br />

2πi<br />

<br />

γ<br />

f (z)<br />

z − w dz.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 26 / <strong>89</strong>


Cauchys integralformel, ver 1.0<br />

Om γ är en sluten stig som g˚ar ett varv runt punkten w och f är analytisk<br />

p˚a γ ∗ och i alla punkter ”innanför” γ ∗ s˚a är<br />

f (w) = 1<br />

2πi<br />

Cauchys integralformel, ver 1.1<br />

<br />

γ<br />

f (z)<br />

z − w dz.<br />

Om γ är en sluten stig, f är analytisk i en öppen mängd Ω s˚a att γ∗ ⊂ Ω<br />

och ν(γ, z) = 0 för alla z ∈ C \ Ω s˚a är<br />

f (w)ν(γ, w) = 1<br />

<br />

f (z)<br />

2πi γ z − w dz,<br />

för alla w ∈ Ω \ γ ∗ .<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 26 / <strong>89</strong>


Exempel<br />

<br />

z + 2<br />

Beräkna integralen<br />

dz d˚a γ är randen (i positiv riktning)<br />

γ (z + 1)(z − 2)<br />

av rektangeln med hörn i punkterna −2 + 3i, −2 − 3i, 1 − 3i och 1 + 3i.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 27 / <strong>89</strong>


Exempel<br />

<br />

z + 2<br />

Beräkna integralen<br />

dz d˚a γ är randen (i positiv riktning)<br />

γ (z + 1)(z − 2)<br />

av rektangeln med hörn i punkterna −2 + 3i, −2 − 3i, 1 − 3i och 1 + 3i.<br />

Lösning: Vi konstaterar att funktionen som skall integreras är analytisk i<br />

alla punkter utom −1 och 2 och av dessa är det bara −1 som ligger<br />

innanför γ∗ z + 2<br />

. Därför väljer vi f (z) = och eftersom f är analytisk i<br />

z − 2<br />

alla punkter utom 2 och eftersom denna punkt ligger utanför rektangeln<br />

kan vi använda Cauchys integralteorem.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 27 / <strong>89</strong>


Exempel<br />

<br />

z + 2<br />

Beräkna integralen<br />

dz d˚a γ är randen (i positiv riktning)<br />

γ (z + 1)(z − 2)<br />

av rektangeln med hörn i punkterna −2 + 3i, −2 − 3i, 1 − 3i och 1 + 3i.<br />

Lösning: Vi konstaterar att funktionen som skall integreras är analytisk i<br />

alla punkter utom −1 och 2 och av dessa är det bara −1 som ligger<br />

innanför γ∗ z + 2<br />

. Därför väljer vi f (z) = och eftersom f är analytisk i<br />

z − 2<br />

alla punkter utom 2 och eftersom denna punkt ligger utanför rektangeln<br />

kan vi använda Cauchys integralteorem. Punkten −1 ligger inne i<br />

rektangeln och vi f˚ar därför<br />

<br />

<br />

z + 2<br />

f (z)<br />

−1 + 2<br />

dz =<br />

dz = 2πif (−1) = 2πi = −2π<br />

(z + 1)(z − 2) z − (−1) −1 − 2 3 i.<br />

γ<br />

γ<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 27 / <strong>89</strong>


Cauchys integralformel, ver 1.2<br />

Om f är en analytisk funktion i en öppen mängd Ω ⊂ C, γ1, γ2, . . . , γn är<br />

slutna stigar s˚a att γ ∗ j ⊂ Ω och n<br />

j=1 ν(γj, z) = 0 för varje z ∈ C \ Ω s˚a är<br />

f (w)<br />

n<br />

j=1<br />

n<br />

<br />

j=1<br />

γj<br />

ν(γj, w) = 1<br />

2πi<br />

f (z) dz = 0<br />

n<br />

<br />

j=1<br />

γj<br />

f (z)<br />

z − w dz, w ∈ Ω \ ∪n j=1γ ∗ j<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 28 / <strong>89</strong>


Teorem<br />

Om f är analytisk i den öppna mängden Ω s˚a är ocks˚a f ′ deriverbar i Ω,<br />

dvs. f oändligt m˚anga g˚anger deriverbar i Ω.<br />

Cauchys olikheter<br />

|f (n) (w)| ≤ n!M<br />

r n<br />

ifall f är analytisk i { z : |z − w| ≤ r } och |f (z)| ≤ M d˚a |z − w| = r<br />

eller ifall f är analytisk i { z : |z − w| < r } och |f (z)| ≤ M d˚a<br />

|z − w| < r.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 29 / <strong>89</strong>


Liouvilles teorem<br />

Om f är analytisk i C och f är begränsad (dvs. |f (z)| ≤ M < ∞ för alla<br />

z ∈ C) s˚a är f en konstant.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 30 / <strong>89</strong>


Liouvilles teorem<br />

Om f är analytisk i C och f är begränsad (dvs. |f (z)| ≤ M < ∞ för alla<br />

z ∈ C) s˚a är f en konstant.<br />

Moreras teorem<br />

<br />

Om f är kontinuerlig i den öppna mängden Ω ⊂ C och<br />

för varje triangel T ⊂ Ω s˚a är f analytisk i Ω.<br />

∂T<br />

f (z) dz = 0<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 30 / <strong>89</strong>


Konvergensradie<br />

Om z0 ∈ C och an ∈ C, n ≥ 0 s˚a finns det ett tal R, 0 ≤ R ≤ ∞, s˚a att<br />

serien ∞<br />

n=0 an(z − z0) n<br />

konvergerar (absolut) d˚a |z − z0| < R<br />

divergerar d˚a |z − z0| > R.<br />

Talet R kallas seriens ∞<br />

n=0 an(z − z0) n konvergensradie.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 31 / <strong>89</strong>


Exempel<br />

Bestäm konvergensradien för potensserien<br />

∞<br />

n=1<br />

n!<br />

n n z2n+2 .<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 32 / <strong>89</strong>


Exempel<br />

Bestäm konvergensradien för potensserien<br />

Lösning: Vi använder kvottestet och räknar<br />

lim<br />

n→∞<br />

<br />

<br />

<br />

<br />

<br />

<br />

(n+1)!<br />

(n+1) n+1 z2(n+1)+2 n!<br />

n n z 2n+2<br />

<br />

<br />

<br />

<br />

<br />

<br />

= lim<br />

n→∞<br />

∞<br />

n=1<br />

n + 1<br />

(n + 1) (n+1)n<br />

n n<br />

n!<br />

n n z2n+2 .<br />

|z| 2n+2−2−2n−2<br />

1<br />

= lim<br />

n→∞ (1 + 1<br />

n )n |z|2 = 1<br />

e |z2 |.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 32 / <strong>89</strong>


Exempel<br />

Bestäm konvergensradien för potensserien<br />

Lösning: Vi använder kvottestet och räknar<br />

lim<br />

n→∞<br />

<br />

<br />

<br />

<br />

<br />

<br />

(n+1)!<br />

(n+1) n+1 z2(n+1)+2 n!<br />

n n z 2n+2<br />

<br />

<br />

<br />

<br />

<br />

<br />

= lim<br />

n→∞<br />

∞<br />

n=1<br />

n + 1<br />

(n + 1) (n+1)n<br />

n n<br />

n!<br />

n n z2n+2 .<br />

|z| 2n+2−2−2n−2<br />

1<br />

= lim<br />

n→∞ (1 + 1<br />

n )n |z|2 = 1<br />

e |z2 |.<br />

Serien konvergerar d˚a gränsvärdet är < 1 dvs. d˚a |z 2 | < e eller |z| < √ e.<br />

P˚a motsvarande sätt divergerar serien d˚a gränsvärdet av kvoterna är > 1,<br />

dvs. d˚a |z| > √ e. Detta innebär att seriens konvergensradie är √ e.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 32 / <strong>89</strong>


Potensserie<br />

Ifall f är analytisk innanför cirkeln |z − z0| = r s˚a är<br />

f (z) =<br />

∞<br />

an(z − z0) n , |z − z0| < r,<br />

n=0<br />

där an = 1<br />

n! f (n) (z0) och konvergensradien är ≥ r.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 33 / <strong>89</strong>


Potensserie<br />

Ifall f är analytisk innanför cirkeln |z − z0| = r s˚a är<br />

f (z) =<br />

∞<br />

an(z − z0) n , |z − z0| < r,<br />

n=0<br />

där an = 1<br />

n! f (n) (z0) och konvergensradien är ≥ r.<br />

Omvänt, ifall f (z) = ∞<br />

n=0 an(z − z0) n s˚a är f analytisk i mängden<br />

{ z : |z − z0| < R } där R är konvergensradien och<br />

f ′ (z) =<br />

∞<br />

nan(z − z0) n−1 , |z − z0| < R<br />

n=1<br />

dvs. serien kan deriveras (och integreras) termvis.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 33 / <strong>89</strong>


Laurent-serie<br />

Om f är analytisk i mängden { z : r < |z − z0| < R } s˚a är<br />

f (z) =<br />

∞<br />

n=−∞<br />

an(z − z0) n ,<br />

där serien konvergerar (absolut) för alla z s˚a att r < |z − z0| < R och kan<br />

deriveras och integreras termvis.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 34 / <strong>89</strong>


Laurent-serie<br />

Om f är analytisk i mängden { z : r < |z − z0| < R } s˚a är<br />

f (z) =<br />

∞<br />

n=−∞<br />

an(z − z0) n ,<br />

där serien konvergerar (absolut) för alla z s˚a att r < |z − z0| < R och kan<br />

deriveras och integreras termvis.<br />

an = 1<br />

2πi<br />

<br />

Cr∗,z 0<br />

f (z)<br />

dz,<br />

(z − z0) n+1<br />

där Cr∗,z0 är cirkeln |z − z0| = r∗, r < r∗ < R,<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 34 / <strong>89</strong>


l’Hopitals regel<br />

Om f och g är analytiska i z0 och f (z0) = g(z0) = 0 och gränsvärdet<br />

f<br />

lim<br />

z→z0<br />

′ (z)<br />

g ′ existerar s˚a är<br />

(z)<br />

f (z) f<br />

lim = lim<br />

z→z0 g(z) z→z0<br />

′ (z)<br />

g ′ (z)<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 35 / <strong>89</strong>


l’Hopitals regel<br />

Om f och g är analytiska i z0 och f (z0) = g(z0) = 0 och gränsvärdet<br />

f<br />

lim<br />

z→z0<br />

′ (z)<br />

g ′ existerar s˚a är<br />

(z)<br />

f (z) f<br />

lim = lim<br />

z→z0 g(z) z→z0<br />

′ (z)<br />

g ′ (z)<br />

Exempel p˚a serieutvecklingar<br />

1<br />

1−z = 1 + z + z2 + z3 + . . . = ∞ n=0 zn , |z| < 1<br />

ez = 1 + z + z2 z3<br />

2! + 3! + . . . = ∞ n=0 zn<br />

n!<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 35 / <strong>89</strong>


l’Hopitals regel<br />

Om f och g är analytiska i z0 och f (z0) = g(z0) = 0 och gränsvärdet<br />

f<br />

lim<br />

z→z0<br />

′ (z)<br />

g ′ existerar s˚a är<br />

(z)<br />

f (z) f<br />

lim = lim<br />

z→z0 g(z) z→z0<br />

′ (z)<br />

g ′ (z)<br />

Exempel p˚a serieutvecklingar<br />

1<br />

1−z = 1 + z + z2 + z3 + . . . = ∞ n=0 zn , |z| < 1<br />

ez = 1 + z + z2 z3<br />

2! + 3! + . . . = ∞ n=0 zn<br />

n!<br />

sin(z) = z − z3 z5<br />

3! + 5! − . . . = ∞ n=0<br />

(−1) n z 2n+1<br />

(2n+1)!<br />

cos(z) = 1 − z2 z4<br />

2! + 4! − . . . = ∞ (−1)<br />

n=0<br />

nz 2n<br />

(2n)!<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 35 / <strong>89</strong>


z-transformationer<br />

Om x(n), n = 0, 1, 2, är en talföljd s˚a är z-transformationen av x<br />

Z(x)(z) =<br />

∞<br />

x(n)z −n<br />

n=0<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 36 / <strong>89</strong>


z-transformationer<br />

Om x(n), n = 0, 1, 2, är en talföljd s˚a är z-transformationen av x<br />

Konvergensomr˚adet<br />

Z(x)(z) =<br />

∞<br />

x(n)z −n<br />

Det finns ett tal R, 0 ≤ R ≤ ∞ s˚a att serien ∞<br />

n=0 x(n)z−n konvergerar<br />

d˚a |z| > R och divergerar d˚a |z| < R.<br />

n=0<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 36 / <strong>89</strong>


Överföringsfunktioner<br />

Antag att H är en funktion som avbildar en talföljd x = (x(n)) ∞ n=0<br />

talföljd H(x) = (H(x)(n)) ∞ n=0 s˚a att<br />

H är linjär, dvs H(λx + µy) = λH(x) + µH(y);<br />

p˚a en<br />

H är translationsinvariant, dvs. H(τ(x)) = τ(H(x)) där τ(x)(0) = 0<br />

och τ(x)(n) = x(n − 1) d˚a n ≥ 1;<br />

Det finns ett tal C < ∞ s˚a att |H(δ)(n)| ≤ C n+1 , n ≥ 0 där δ(0) = 1<br />

och δ(n) = 0, n ≥ 1.<br />

d˚a finns det en funktion H(z), den sk. överföringsfunktionen s˚a att<br />

Y (z) = H(z)X (z)<br />

där X = Z(x) och Y = Z(H(x)). Funktionen H(z) är z-transformationen<br />

av H(δ).<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 37 / <strong>89</strong>


Residy<br />

Om f (z) =<br />

∞<br />

n=−∞<br />

residyn av f i punkten z0.<br />

an(z − z0) n d˚a 0 < |z − z0| < R s˚a är a−1 = Res(f , z0)<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 38 / <strong>89</strong>


Residy<br />

Om f (z) =<br />

∞<br />

n=−∞<br />

residyn av f i punkten z0.<br />

an(z − z0) n d˚a 0 < |z − z0| < R s˚a är a−1 = Res(f , z0)<br />

Lemma<br />

Om<br />

<br />

lim (z − z0)f (z)<br />

z→z0<br />

= w existerar<br />

eller d˚a m > 1<br />

1 d<br />

lim<br />

z→z0 (m − 1)!<br />

m−1<br />

dzm−1 <br />

(z − z0) m f (z) = w existerar<br />

s˚a är Res(f , z0) = w.<br />

(m m˚aste vara minst ordningen av polen i z0.)<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 38 / <strong>89</strong>


Residyteoremet<br />

Ifall γ är en sluten stig och f är analytisk p˚a och innanför γ ∗ bortsett fr˚an<br />

punkterna zj för vilka ν(γ, zj) = 1, j = 1, . . . , k s˚a är<br />

<br />

γ<br />

f (z) dz = 2πi<br />

k<br />

Res(f , zj).<br />

j=1<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 39 / <strong>89</strong>


Varför?<br />

Om γj är en cirkel med mittpunkt zj och s˚a liten radie rj att f är analytisk<br />

i { z : 0 < |z − zj| ≤ rj } ⊂ Ω s˚a kan vi skriva<br />

Eftersom d<br />

dz<br />

f (z) =<br />

∞<br />

n=−∞<br />

aj,n(z − zj) n , 0 < |z − zj| ≤ rj.<br />

1<br />

n+1zn+1 = zn d˚a n = −1 blir <br />

γj aj,n(z − zj) n dz = 0 d˚a<br />

n = −1 och<br />

<br />

<br />

f (z) dz = aj,−1<br />

γj<br />

γj<br />

1<br />

dz = 2πiaj,−1 = 2πiRes(f , zj).<br />

z − zj<br />

Nu följer residyteoremet av att man tillämpar Cauchys integralteorem, ver<br />

1.2 p˚a stigarna γ och γ ↢ j , j = 1, . . . , k.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 40 / <strong>89</strong>


Beräkning av integraler<br />

Om f är analytisk p˚a reella axeln och i övre halvplanet (Im (z) > 0)<br />

|zf (z)| = 0, (dvs. om<br />

bortsett fr˚an punkterna zj, j = 1, . . . , k, och limIm z≥0<br />

|z|→∞<br />

f är rationell, ifall skillnaden mellan nämnarens och täljarens gradtal är minst 2)<br />

s˚a är<br />

∞<br />

−∞<br />

R<br />

f (x) dx = lim f (x) dx = 2πi<br />

R→∞ −R<br />

k<br />

Res(f , zj).<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 41 / <strong>89</strong><br />

j=1


Exempel<br />

∞<br />

Beräkna med hjälp av residyteoremet integralen<br />

1<br />

−∞<br />

1<br />

(x 2 + 1)(x 2 + 4) dx.<br />

Lösning: L˚at f (z) = (z2 +1)(z2 . Denna funktion är analytisk i alla<br />

+4)<br />

punkter utom nämnarens nollställen vilka är<br />

z = ±i och z = ±2i.<br />

Av dessa är det endast +i och +2i som ligger i övre halvplanet. Eftersom<br />

nämnarens gradtal minus täljarens gradtal är 4 s˚a ser vi ocks˚a att<br />

lim |z||f (z)| = 0.<br />

|z|→∞<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 42 / <strong>89</strong>


forts.<br />

Vi kan allts˚a använda residyteoremet och vi f˚ar<br />

∞ 1<br />

−∞ (x 2 + 1)(x 2 + 4) dx = 2πiRes(f , i) + Res(f , 2i) <br />

<br />

= 2πi lim (z − i)f (z) + 2πi lim (z − 2i)f (z)<br />

z→i<br />

z→2i<br />

z − i<br />

= 2πi lim<br />

z→i (z − i)(z + i)(z − 2i)(z + 2i)<br />

z − 2i<br />

+ 2πi lim<br />

z→2i (z − i)(z + i)(z − 2i)(z + 2i)<br />

1<br />

1<br />

= 2πi lim<br />

+ 2πi lim<br />

z→i (z + i)(z − 2i)(z + 2i) z→2i (z − i)(z + i)(z + 2i)<br />

2πi 2πi π π π<br />

= + = − =<br />

2i (−i) 3i i 3i 4i 3 6 6 .<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 43 / <strong>89</strong>


Teorem<br />

Ifall a > 0, g är analytisk p˚a reella axeln och i övre halvplanet (Im (z) > 0)<br />

|g(z)| = 0 (när g är<br />

bortsett fr˚an punkterna zj, j = 1, . . . , k, och limIm z≥0<br />

|z|→∞<br />

rationell, ifall skillnaden mellan nämnarens och täljarens gradtal är minst 1) s˚a är<br />

∞<br />

−∞<br />

e iax R<br />

g(x) dx = lim<br />

R→∞<br />

e<br />

−R<br />

iax g(x) dx = 2πi<br />

k<br />

Res(e iaz g(z), zj).<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 44 / <strong>89</strong><br />

j=1


Teorem<br />

Ifall a > 0, g är analytisk p˚a reella axeln och i övre halvplanet (Im (z) > 0)<br />

|g(z)| = 0 (när g är<br />

bortsett fr˚an punkterna zj, j = 1, . . . , k, och limIm z≥0<br />

|z|→∞<br />

rationell, ifall skillnaden mellan nämnarens och täljarens gradtal är minst 1) s˚a är<br />

∞<br />

−∞<br />

Obs!<br />

e iax R<br />

g(x) dx = lim<br />

R→∞<br />

e<br />

−R<br />

iax g(x) dx = 2πi<br />

k<br />

Res(e iaz g(z), zj).<br />

j=1<br />

R<br />

−R eiax <br />

g(x) dx och<br />

Om g(x) är reell s˚a är R<br />

−R sin(ax)g(x) dx = Im<br />

<br />

R<br />

R<br />

−R cos(ax)g(x) dx = Re −R eiax <br />

g(x) dx , och det är i b˚ada fallen<br />

enklare att räkna ut limR→∞<br />

R<br />

−R eiax g(x) dx och sedan ta imaginära eller<br />

reella delen än att skriva tex. sin(ax) = 1<br />

2i (eiax − e iax ). (Att försöka<br />

tillämpa residyteoremet p˚a funktionen sin(az)g(z) eller cos(az)g(z) lyckas<br />

vanligen inte alls.)<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 44 / <strong>89</strong>


Lemma<br />

2π<br />

0<br />

F (cos(t), sin(t)) dt<br />

=<br />

ifall f (z) = F<br />

2π<br />

0<br />

F <br />

1 it −it<br />

2 e + e , 1<br />

it −it<br />

2i e − e 1<br />

ieit ieit <br />

dt = f (z) dz<br />

γ<br />

<br />

1<br />

2 z + 1<br />

<br />

, 1<br />

<br />

2i z − 1<br />

<br />

1<br />

z<br />

z<br />

iz<br />

och γ är enhetscirkeln<br />

|z| = 1. Integralen kan räknas ut med residyteoremet (ifall förutsättningarna<br />

är uppfyllda).<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 45 / <strong>89</strong>


Nollställen, poler<br />

∞<br />

Om f (z) = an(z − z0) n d˚a 0 < |z − z0| < R där am = 0 och R > 0,<br />

n=m<br />

(och f (z0) = a0 om m ≥ 0) s˚a är<br />

z0 ett nollställe av ordningen m om m > 0 (dvs. med multipliciteten m),<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 46 / <strong>89</strong>


Nollställen, poler<br />

∞<br />

Om f (z) = an(z − z0) n d˚a 0 < |z − z0| < R där am = 0 och R > 0,<br />

n=m<br />

(och f (z0) = a0 om m ≥ 0) s˚a är<br />

z0 ett nollställe av ordningen m om m > 0 (dvs. med multipliciteten m),<br />

z0 är en pol av ordningen |m| om m < 0 (dvs. med multipliciteten |m|).<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 46 / <strong>89</strong>


Nollställen, poler<br />

∞<br />

Om f (z) = an(z − z0) n d˚a 0 < |z − z0| < R där am = 0 och R > 0,<br />

n=m<br />

(och f (z0) = a0 om m ≥ 0) s˚a är<br />

z0 ett nollställe av ordningen m om m > 0 (dvs. med multipliciteten m),<br />

z0 är en pol av ordningen |m| om m < 0 (dvs. med multipliciteten |m|).<br />

Om f (z) =<br />

∞<br />

an(z − z0) n och an = 0 för oändligt m˚anga n < 0 s˚a är z0 en<br />

n=−∞<br />

väsentlig singularitet.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 46 / <strong>89</strong>


Nollställen, poler<br />

∞<br />

Om f (z) = an(z − z0) n d˚a 0 < |z − z0| < R där am = 0 och R > 0,<br />

n=m<br />

(och f (z0) = a0 om m ≥ 0) s˚a är<br />

z0 ett nollställe av ordningen m om m > 0 (dvs. med multipliciteten m),<br />

z0 är en pol av ordningen |m| om m < 0 (dvs. med multipliciteten |m|).<br />

Om f (z) =<br />

∞<br />

an(z − z0) n och an = 0 för oändligt m˚anga n < 0 s˚a är z0 en<br />

n=−∞<br />

väsentlig singularitet.<br />

Exempel<br />

0 är ett nollställe av ordningen 1 till sin(z)<br />

0 är ett nollställe av ordningen 2 till 1 − cos(z)<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 46 / <strong>89</strong>


Lemma<br />

z0 är ett nollställe av ordningen m till f<br />

⇔<br />

f (z0) = f ′ (z0) = . . . = f (m−1) (z0) = 0 och f (m) (z0) = 0<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 47 / <strong>89</strong>


Lemma<br />

z0 är ett nollställe av ordningen m till f<br />

⇔<br />

f (z0) = f ′ (z0) = . . . = f (m−1) (z0) = 0 och f (m) (z0) = 0<br />

⇔<br />

lim<br />

z→z0<br />

f (z)<br />

existerar och är = 0<br />

(z − z0) m<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 47 / <strong>89</strong>


Lemma<br />

z0 är ett nollställe av ordningen m till f<br />

⇔<br />

f (z0) = f ′ (z0) = . . . = f (m−1) (z0) = 0 och f (m) (z0) = 0<br />

⇔<br />

lim<br />

z→z0<br />

⇔<br />

f (z)<br />

existerar och är = 0<br />

(z − z0) m<br />

z0 är en pol av ordningen m till 1<br />

f<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 47 / <strong>89</strong>


Lemma<br />

z0 är ett nollställe av ordningen m till f<br />

⇔<br />

f (z0) = f ′ (z0) = . . . = f (m−1) (z0) = 0 och f (m) (z0) = 0<br />

⇔<br />

lim<br />

z→z0<br />

⇔<br />

f (z)<br />

existerar och är = 0<br />

(z − z0) m<br />

z0 är en pol av ordningen m till 1<br />

f<br />

Lemma<br />

z0 är en pol av ordningen m till f<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 47 / <strong>89</strong>


Lemma<br />

z0 är ett nollställe av ordningen m till f<br />

⇔<br />

f (z0) = f ′ (z0) = . . . = f (m−1) (z0) = 0 och f (m) (z0) = 0<br />

⇔<br />

lim<br />

z→z0<br />

⇔<br />

f (z)<br />

existerar och är = 0<br />

(z − z0) m<br />

z0 är en pol av ordningen m till 1<br />

f<br />

Lemma<br />

z0 är en pol av ordningen m till f<br />

⇔<br />

limz→z0 (z − z0) m f (z) existerar och är = 0<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 47 / <strong>89</strong>


Lemma<br />

z0 är ett nollställe av ordningen m till f<br />

⇔<br />

f (z0) = f ′ (z0) = . . . = f (m−1) (z0) = 0 och f (m) (z0) = 0<br />

⇔<br />

lim<br />

z→z0<br />

⇔<br />

f (z)<br />

existerar och är = 0<br />

(z − z0) m<br />

z0 är en pol av ordningen m till 1<br />

f<br />

Lemma<br />

z0 är en pol av ordningen m till f<br />

⇔<br />

limz→z0 (z − z0) m f (z) existerar och är = 0<br />

⇔<br />

z0 är ett nollställe av ordningen m till 1<br />

f<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 47 / <strong>89</strong>


Lemma<br />

Om f är analytisk i z0 och z0 är ett nollställe av ordningen m s˚a är<br />

<br />

f ′<br />

Res , z0 = m<br />

f<br />

Lemma<br />

Om f är analytisk i mängden { z : 0 < |z − z0| < r } där r > 0 och z0 är<br />

en pol av ordningen m s˚a är<br />

<br />

f ′<br />

Res , z0 = −m<br />

f<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 48 / <strong>89</strong>


Antalet lösningar till en ekvation innanför en sluten stig<br />

Om f är analytisk p˚a och innanför γ ∗ (och den slutna stigen γ g˚ar ett varv<br />

i positiv riktning) och ϕ, där ϕ(t) = f (γ(t)), t ∈ [a, b], g˚ar m varv runt w<br />

s˚a finns det m lösningar till ekvationen f (z) = w innanför γ ∗<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 49 / <strong>89</strong>


Antalet lösningar till en ekvation innanför en sluten stig<br />

Om f är analytisk p˚a och innanför γ ∗ (och den slutna stigen γ g˚ar ett varv<br />

i positiv riktning) och ϕ, där ϕ(t) = f (γ(t)), t ∈ [a, b], g˚ar m varv runt w<br />

s˚a finns det m lösningar till ekvationen f (z) = w innanför γ ∗<br />

Mera exakt<br />

Ifall ν(γ, z) = 1 eller 0 d˚a z ∈ C \ γ ∗ och f är analytisk p˚a och innanför<br />

γ ∗ bortsett fr˚an högst ett ändligt antal poler innanför γ ∗ s˚a är<br />

1<br />

2πi<br />

<br />

γ<br />

f ′ (z) 1<br />

dz =<br />

f (z) − w 2πi<br />

<br />

ϕ<br />

1<br />

dz = ν(ϕ, w),<br />

z − w<br />

där ϕ(t) = f (γ(t)), t ∈ [a, b], antalet lösningar till ekvationen f (z) = w<br />

minus antalet poler innanför γ ∗ (räknade med multiplicitet).<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 49 / <strong>89</strong>


Lemma<br />

Om f är analytisk i punkten z0 och f (z0) = 0 s˚a finns det ett tal r > 0 s˚a<br />

att antingen gäller<br />

f (z) = 0 d˚a |z − z0| < r<br />

eller<br />

f (z) = 0 d˚a 0 < |z − z0| < r.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 50 / <strong>89</strong>


Lemma<br />

Om f är analytisk i punkten z0 och f (z0) = 0 s˚a finns det ett tal r > 0 s˚a<br />

att antingen gäller<br />

Ifall<br />

f (z) = 0 d˚a |z − z0| < r<br />

eller<br />

f (z) = 0 d˚a 0 < |z − z0| < r.<br />

Analytisk fortsättning<br />

Ω är en öppen sammanhängande mängd,<br />

f och g är analytiska i Ω,<br />

f (z) = g(z) d˚a z ∈ A ⊂ Ω<br />

◮ där A är s˚adan att det finns zj ∈ A, zj = z∞, j ≥ 1 och<br />

limj→∞ zj = z∞ ∈ Ω,<br />

◮ tex. s˚a att A = ∅ och antingen är öppen eller A = Ω ∩ γ ∗ där γ är en<br />

stig (s˚a att γ ∗ inte best˚ar av bara en punkt),<br />

s˚a är f (z) = g(z) för alla z ∈ Ω.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 50 / <strong>89</strong>


Exempel<br />

Definiera Γ(z) = ∞<br />

0 e−t t z−1 dt. Man kan visa att Γ d˚a kommer att vara<br />

analytisk d˚a Re (z) > 0 (eftersom |t z−1 | = t Re (z)−1 d˚a t > 0). Hur är det<br />

med Γ(z) för andra värden av z?<br />

Lösning: D˚a Re (z) > 0 f˚ar vi med hjälp av partiell integrering<br />

∞<br />

Γ(z + 1) = e −t t z ∞ −t z<br />

dt = −e t ∞<br />

−<br />

Detta betyder att Γ(z) =<br />

0<br />

Γ(z) =<br />

0<br />

Γ(z + 1)<br />

z<br />

0<br />

−e −t zt z−1 dt = zΓ(z).<br />

d˚a Re (z) > 0 och vi kan definiera<br />

Γ(z + 1)<br />

, z = 0, Re (z) > −1.<br />

z<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 51 / <strong>89</strong>


forts.<br />

Av teoremet om analytisk fortsättning följer att ifall f (z) är en analytisk<br />

funktion i mängden { z : z = 0, Re (z) > −1 } s˚a att f (z) = Γ(z) d˚a<br />

Re (z) > 0 s˚a m˚aste den vara Γ(z+1)<br />

z . Vi kan nu fortsätta p˚a samma sätt<br />

och definiera<br />

Γ(z) =<br />

Γ(z + k)<br />

, z = −j, j = 0, . . . , k−1, Re (z) > −k,<br />

z(z + 1) . . . (z + k − 1)<br />

och vi ser vi f˚ar en funktion Γ som är analytisk i alla punkter utom<br />

{ 0, −1, −2, . . . }.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 52 / <strong>89</strong>


Definition<br />

En funktion u(x1, . . . , xn) är harmonisk i (den öppna) mängden Ω ⊂ R n<br />

om den för alla (x1, . . . , xn) ∈ Ω uppfyller Laplace-ekvationen<br />

∆u = ux1x1 (x1, . . . , xn) + . . . + uxnxn(x1, . . . , xn) = 0.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 53 / <strong>89</strong>


Definition<br />

En funktion u(x1, . . . , xn) är harmonisk i (den öppna) mängden Ω ⊂ R n<br />

om den för alla (x1, . . . , xn) ∈ Ω uppfyller Laplace-ekvationen<br />

Teorem<br />

∆u = ux1x1 (x1, . . . , xn) + . . . + uxnxn(x1, . . . , xn) = 0.<br />

Funktionen u(x, y) är harmonisk i Ω<br />

⇔<br />

u(x, y) = Re (f (x + iy)) där f (z) är en analytisk funktion i Ω (d˚a<br />

x + iy ′′ = ′′ (x, y) ∈ Ω).<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 53 / <strong>89</strong>


Definition<br />

En funktion u(x1, . . . , xn) är harmonisk i (den öppna) mängden Ω ⊂ R n<br />

om den för alla (x1, . . . , xn) ∈ Ω uppfyller Laplace-ekvationen<br />

Teorem<br />

∆u = ux1x1 (x1, . . . , xn) + . . . + uxnxn(x1, . . . , xn) = 0.<br />

Funktionen u(x, y) är harmonisk i Ω<br />

⇔<br />

u(x, y) = Re (f (x + iy)) där f (z) är en analytisk funktion i Ω (d˚a<br />

x + iy ′′ = ′′ (x, y) ∈ Ω).<br />

Lemma<br />

Om v(x, y) är en harmonisk funktion i (den öppna) mängden Ω1 och<br />

f : Ω2 → Ω1 är en analytisk funktion i (den öppna) mängden Ω2 s˚a är<br />

u(x, y) = v(Re (f (x + iy)), Im (f (x + iy))) harmonisk i Ω2.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 53 / <strong>89</strong>


Exempel<br />

Bestäm en harmonisk funktion (uxx + uyy = 0) i första kvadranten<br />

{ (x, y) : x > 0, y > 0 } s˚a att u(x, 0) = 0, x > 0 och u(0, y) = 4, y > 0<br />

genom att visa att v(x, y) = 1 1 x<br />

2 − π arctan( y ) är en harmonisk funktion d˚a<br />

y > 0 s˚a att v(x, 0) = 1 d˚a x < 0 och v(x, 0) = 0 d˚a x > 0 och genom<br />

att använda den analytiska funktionen z ↦→ z2 .<br />

Lösning: När vi deriverar f˚ar vi<br />

och<br />

vx(x, y) = − 1<br />

π<br />

vy (x, y) = − 1<br />

π<br />

1<br />

y<br />

1 + x2<br />

y 2<br />

− x<br />

y 2<br />

1 + x2<br />

y 2<br />

= − 1<br />

π<br />

= 1<br />

π<br />

y<br />

x 2 + y 2 och vxx(x, y) = 1<br />

π<br />

x<br />

x 2 + y 2 och vyy(x, y) = − 1<br />

π<br />

Av detta är det klart att vxx + vyy = 0 s˚a v är harmonisk.<br />

2xy<br />

(x 2 + y 2 ,<br />

) 2<br />

2xy<br />

(x 2 + y 2 ) 2<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 54 / <strong>89</strong>


forts.<br />

Om x > 0 s˚a är limy→0+ x<br />

y = +∞ s˚a att limy→0+ v(x, y) = 1<br />

2<br />

och om x < 0 s˚a är limy→0+ x<br />

y = −∞ s˚a att<br />

1 π − π 2<br />

= 0<br />

limy→0+ v(x, y) = 1 1 π<br />

2 − π (− 2 ) = 1. Eftersom f (z) = z2 är analytisk och<br />

avbildar mängden { z ∈ C : Re (z) > 0, Im (z) > 0 } p˚a mängden<br />

{ z ∈ C : Im (z) > 0 } s˚a är funktionen<br />

u(x, y) = cv(Re (f (x + iy)), Im (f (x + iy))) harmonisk i mängden<br />

{ (x, y) : x > 0, y > 0 }. Dessutom ser vi att mängden { (0, y) : y > 0 }<br />

avbildas mängden { (x, 0) : x < 0 } s˚a att om vi vill att u(0, y) = 4 och s˚a<br />

skall vi välja c = 4. P˚a motsvarande sätt ser vi att mängden<br />

{ (x, 0) : x > 0 } avbildas p˚a mängden { (x, 0) : x > 0 } s˚a att vi ocks˚a<br />

automatiskt f˚ar u(x, 0) = 0 d˚a u definieras som ovan. Eftersom<br />

f (x + iy) = (x + iy) 2 = x 2 + 2ixy − y 2 s˚a blir<br />

u(x, y) = 4v(Re (f (x + iy)), Im (f (x + iy)))<br />

= 4v(x 2 − y 2 , 2xy) = 4 1 1<br />

− 4<br />

2 π arctan<br />

<br />

x 2 − y 2<br />

2xy<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 55 / <strong>89</strong>


Harmoniska funktioner i polära koordinater<br />

Antag att w(r, θ) = v(r cos(θ), r sin(θ)) (dvs. w är v uttryckt med polära<br />

koordinater). D˚a uppfyller v Laplace-ekvationen vxx + vyy = 0, (dvs. den<br />

är harmonisk) om och endast om<br />

wrr + 1<br />

r wr + 1<br />

r 2 wθθ = 0.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 56 / <strong>89</strong>


Harmoniska funktioner i polära koordinater<br />

Antag att w(r, θ) = v(r cos(θ), r sin(θ)) (dvs. w är v uttryckt med polära<br />

koordinater). D˚a uppfyller v Laplace-ekvationen vxx + vyy = 0, (dvs. den<br />

är harmonisk) om och endast om<br />

Medelvärdesegenskapen<br />

Om Ω ⊂ R n är öppen s˚a gäller:<br />

wrr + 1<br />

r wr + 1<br />

r 2 wθθ = 0.<br />

u är tv˚a g˚anger deriverbar och harmonisk (∆u = 0) i Ω<br />

⇔<br />

u är oändligt m˚anga g˚anger deriverbar och harmonisk (∆u = 0) i Ω<br />

⇔<br />

<br />

1<br />

u är kontinuerlig i Ω och u(x0) = R<br />

x=1 1 dS x=1 u(x0 + rx) dS, för<br />

alla x0 och r ≥ 0 s˚a att { x ∈ Rn : x − x0 ≤ r } ⊂ Ω<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 56 / <strong>89</strong>


Maximumprincipen<br />

Om Ω är en öppen begränsad mängd i R n med rand ∂Ω och u är<br />

kontinuerlig i Ω ∪ ∂Ω och harmonisk i Ω s˚a är<br />

max u(x) = max<br />

x∈Ω∪∂Ω x∈∂Ω u(x).<br />

Om Ω dessutom är sammanhängande s˚a gäller att om u(x0) = maxx∈Ω∪∂Ω u(x)<br />

för n˚agot x0 ∈ Ω s˚a är u en konstant funktion.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 57 / <strong>89</strong>


Maximumprincipen<br />

Om Ω är en öppen begränsad mängd i R n med rand ∂Ω och u är<br />

kontinuerlig i Ω ∪ ∂Ω och harmonisk i Ω s˚a är<br />

max u(x) = max<br />

x∈Ω∪∂Ω x∈∂Ω u(x).<br />

Om Ω dessutom är sammanhängande s˚a gäller att om u(x0) = maxx∈Ω∪∂Ω u(x)<br />

för n˚agot x0 ∈ Ω s˚a är u en konstant funktion.<br />

Entydighet<br />

Om Ω är en öppen begränsad mängd i R n med rand ∂Ω och u1 och u2 är<br />

kontinuerliga i Ω ∪ ∂Ω och harmoniska i Ω och u1(x) = u2(x) d˚a x ∈ ∂Ω<br />

s˚a är u1(x) = u2(x) för alla x ∈ Ω.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 57 / <strong>89</strong>


Poissons formel i enhetscirkeln<br />

Om u(x, y) är t.ex. begränsad d˚a x 2 + y 2 < 1 s˚a är u harmonisk i<br />

{ (x, y) : x 2 + y 2 < 1 } ⇔<br />

u(r cos(θ), r sin(θ))<br />

= 1<br />

2π<br />

2π<br />

0<br />

1 − r 2<br />

1 − 2r cos(θ − t) + r 2 u cos(t), sin(t) dt, r < 1.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 58 / <strong>89</strong>


Poissons formel i enhetscirkeln<br />

Om u(x, y) är t.ex. begränsad d˚a x 2 + y 2 < 1 s˚a är u harmonisk i<br />

{ (x, y) : x 2 + y 2 < 1 } ⇔<br />

u(r cos(θ), r sin(θ))<br />

= 1<br />

2π<br />

2π<br />

0<br />

Poissons formel i övre halvplanet<br />

1 − r 2<br />

1 − 2r cos(θ − t) + r 2 u cos(t), sin(t) dt, r < 1.<br />

Om u(x, y) är t.ex. begränsad d˚a x ∈ R och y > 0 s˚a är u harmonisk i<br />

{ (x, y) : x ∈ R, y > 0 } ⇔<br />

u(x, y) = 1<br />

∞ y<br />

π (x − t) 2 u(t, 0) dt, y > 0.<br />

+ y 2<br />

−∞<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 58 / <strong>89</strong>


Strömningsproblem i tv˚a dimensioner<br />

Hastighetsvektorn för en virvelfri strömning av en inkompressibel ”vätska”<br />

med viskositet 0 (dvs. utan ”friktion”) är v(x, y) s˚a att v = ∇Φ och<br />

∇ · v = 0. Dessa villkor innebär att 0 = ∇ · ∇Φ = Φxx + Φyy, dvs. Φ är<br />

harmonisk.<br />

Om nu f (x + iy) = Φ(x, y) + iΨ(x, y) är analytisk s˚a följer av<br />

Cauchy-Riemann ekvationerna att hastighetsvektorn<br />

(Φx, Φy ) = (Re (f ′ ), Im (f ′ )) och (Φx, Φy ) ⊥ (Ψx, Ψy ), dvs.<br />

hastighetsvektorn är parallell med kurvan Ψ = c (där c är en konstant)<br />

som allts˚a är en ”strömningslinje”.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 59 / <strong>89</strong>


Exempel<br />

Vilket strömningsproblem kan beskrivas med funktionen f (z) = z 2 i<br />

mängden { z ∈ C : Im z > 0 }?<br />

Lösning: I dethär fallet är f (x + iy) = x 2 − y 2 + i2xy. Strömningslinjerna<br />

är allts˚a kurvorna xy = c där c är en konstant och hastighetsvektorn är<br />

(2x, −2y) = 2xi − 2yj. Av detta kan vi se att det är fr˚agan om en<br />

strömning mot den reella axeln som delar sig och g˚ar ut mot den positiva<br />

eller negativa sidan.<br />

.<br />

.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 60 / <strong>89</strong>


Exempel<br />

Vilken situation i strömningsmekanik kan beskrivas med funktionen<br />

f (z) = z + 1<br />

z (d˚a |z| ≥ 1)?<br />

Lösning: Vi skriver f = Φ + iΨ och om vi använder polära koordinater,<br />

dvs. väljer z = reiθ s˚a har vi<br />

f (z) = re iθ + 1<br />

r e−iθ <br />

= cos(θ) r + 1<br />

<br />

+ i sin(θ) r −<br />

r<br />

1<br />

<br />

.<br />

r<br />

Eftersom f ′ (z) = 1 − 1<br />

z 2 s˚a ser vi att f ′ (z) → 1 d˚a |z| → ∞, dvs. l˚angt<br />

borta fr˚an origo är hastighetsvektorn i stort sett en enhetsvektor i x-axelns<br />

riktning. Strömningslinjerna bestäms av ekvationerna Ψ = c, dvs.<br />

θ = arcsin( cr<br />

r 2 cr<br />

) och θ = π − arcsin( −1 r 2 ). Om c = 0 gäller antingen<br />

−1<br />

r = 1, θ = 0 eller θ = π och det är fr˚agan om strömning i omr˚adet utanför<br />

enhetscirkeln med mittpunkt i origo.<br />

.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 61 / <strong>89</strong><br />

.


Laplace-transformationen<br />

∞<br />

L(f )(s) = e −st f (t) dt<br />

0<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 62 / <strong>89</strong>


Laplace-transformationen<br />

∞<br />

L(f )(s) = e −st f (t) dt<br />

Exempel<br />

L(1)(s) = 1<br />

s<br />

L(t)(s) = 1<br />

s 2<br />

L(e at )(s) = 1<br />

s−a<br />

L(cos(ωt))(s) = s<br />

s 2 +ω 2<br />

L(sin(ωt))(s) = ω<br />

s 2 +ω 2<br />

0<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 62 / <strong>89</strong>


Laplace-transformationen<br />

∞<br />

L(f )(s) = e −st f (t) dt<br />

Exempel<br />

L(1)(s) = 1<br />

s<br />

L(t)(s) = 1<br />

s 2<br />

L(e at )(s) = 1<br />

s−a<br />

L(cos(ωt))(s) = s<br />

s 2 +ω 2<br />

L(sin(ωt))(s) = ω<br />

s 2 +ω 2<br />

Laplace-transformationen är linjär!<br />

0<br />

L(αf + βg) = αL(f ) + βL(g)<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 62 / <strong>89</strong>


Teorem<br />

Ifall T<br />

0<br />

|f (t)| dt < ∞ för alla T > 0<br />

T<br />

0 e−s0t f (t) dt existerar för n˚agot tal s0 ∈ C<br />

limT →∞<br />

s˚a gäller att<br />

T<br />

F (s) = limT →∞ 0 e−stf (t) dt existerar d˚a Re (s) > Re (s0),<br />

F (s) är analytisk i mängden { s ∈ C : Re (s) > Re (s0) }.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 63 / <strong>89</strong>


Teorem<br />

Ifall T<br />

0<br />

|f (t)| dt < ∞ för alla T > 0<br />

T<br />

0 e−s0t f (t) dt existerar för n˚agot tal s0 ∈ C<br />

limT →∞<br />

s˚a gäller att<br />

T<br />

F (s) = limT →∞ 0 e−stf (t) dt existerar d˚a Re (s) > Re (s0),<br />

F (s) är analytisk i mängden { s ∈ C : Re (s) > Re (s0) }.<br />

Teorem<br />

Laplace-transformationen är entydig, dvs. om L(f )(s) = L(g)(s) d˚a<br />

Re (s) > α s˚a är f(t)=g(t) för nästan alla t ≥ 0.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 63 / <strong>89</strong>


Räkneregler, derivator mm. d˚a F (s) = L(f )(s)<br />

L(f ′ )(s) = sF (s) − f (0)<br />

F ′ (s) = L(−tf (t))(s)<br />

L(f ′′ )(s) = s 2 F (s) − sf (0) − f ′ (0)<br />

L<br />

t<br />

0 f (τ) dτ<br />

<br />

(s) = 1<br />

s F (s)<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 64 / <strong>89</strong>


Räkneregler, derivator mm. d˚a F (s) = L(f )(s)<br />

L(f ′ )(s) = sF (s) − f (0)<br />

F ′ (s) = L(−tf (t))(s)<br />

L(f ′′ )(s) = s 2 F (s) − sf (0) − f ′ (0)<br />

L<br />

t<br />

0 f (τ) dτ<br />

<br />

(s) = 1<br />

s F (s)<br />

Räkneregler, förskjutningsregler mm. d˚a F (s) = L(f )(s)<br />

L(e at f (t))(s) = F (s − a)<br />

L(f (t − a)u(t − a))(s) = e −as F (s)<br />

där u(t) = 1 d˚a t > 0, u(t) = 0 d˚a t < 0 och a ≥ 0.<br />

L f (at) (s) = 1<br />

aF <br />

s<br />

a , a > 0<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 64 / <strong>89</strong>


Exempel<br />

Lös begynnelsevärdesproblemet<br />

y ′′ (t) + 4y ′ (t) + 3y(t) = 4 + 3t, y(0) = 3, y ′ (0) = −6,<br />

med hjälp av Laplace-transformering.<br />

Lösning: L˚at Y (s) = L(y(t)) s˚a att L(y ′ (t)) = sY (s) − y(0) = sY (s) − 3<br />

och<br />

L(y ′′ (t)) = sL(y ′ (t)) − y ′ (0) = s2Y (s) − sy(0) − y ′ (0) = s2Y (s) − 3s + 6.<br />

Eftersom L(4 + 3t) = 4 3<br />

s + s2 = 4s+3<br />

s2 s˚a f˚ar vi genom att ta<br />

Laplace-transformationer av b˚ada sidorna i ekvationen<br />

s 2 Y (s) − 3s + 6 + 4sY (s) − 12 + 3Y (s) =<br />

s˚a att vi f˚ar d˚a vi löser ekvationen<br />

Y (s) =<br />

3s + 6<br />

s2 + 4s + 3 +<br />

4s + 3<br />

s2 (s2 + 4s + 3) .<br />

4s + 3<br />

s 2 ,<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 65 / <strong>89</strong>


forts.<br />

För att kunna göra en partialbr˚aksuppdelning m˚aste vi först faktorisera<br />

s 2 + 4s + 3 och d˚a behöver vi känna till nollställena som är<br />

s = −2 ± √ 4 − 3 = −2 ± 1, dvs. −1 och −3. Partialbr˚aksuppdelningen blir<br />

3s + 6<br />

s 2 + 4s + 3 +<br />

4s + 3<br />

s2 (s2 A B<br />

= +<br />

+ 4s + 3) s + 1 s + 3<br />

C D<br />

+ +<br />

s2 s .<br />

Koefficienterna A och B är residyer till enkla poler s˚a att<br />

<br />

3s + 6<br />

A = lim (s + 1)<br />

s→−1 (s + 1)(s + 3) +<br />

4s + 3<br />

s2 <br />

=<br />

(s + 1)(s + 3)<br />

3 1<br />

− = 1,<br />

2 2<br />

<br />

3s + 6<br />

B = lim (s + 3)<br />

s→−3 (s + 1)(s + 3) +<br />

4s + 3<br />

s2 <br />

=<br />

(s + 1)(s + 3)<br />

3 1<br />

+ = 2, .<br />

2 2<br />

För att f˚a C multiplicerar vi med s2 och l˚ater s → 0:<br />

C = lim s<br />

s→0 2<br />

<br />

3s + 6<br />

(s + 1)(s + 3) +<br />

4s + 3<br />

s2 <br />

(s + 1)(s + 3)<br />

= 1.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 66 / <strong>89</strong>


forts.<br />

Koefficienten D kan t.ex. räknas ut genom att multiplicera med s2 sedan<br />

derivera och sedan l˚ata s → 0:<br />

<br />

d<br />

D = lim s<br />

s→0 ds<br />

2<br />

<br />

3s + 6<br />

(s + 1)(s + 3) +<br />

4s + 3<br />

s2 <br />

(s + 1)(s + 3)<br />

<br />

3s + 6<br />

= lim 2s<br />

s→0 (s + 1)(s + 3) + s2 3(s2 + 4s + 3) − (3s + 6)(2s + 4)<br />

(s2 + 4s + 3) 2<br />

+ 4(s2 + 4s + 3) − (4s + 3)(2s + 4)<br />

(s2 + 4s + 3) 2<br />

<br />

= 0.<br />

Nu ser vi med hjälp av partialbr˚aksuppdelningen att<br />

L −1<br />

<br />

3s + 6<br />

s2 + 4s + 3 +<br />

4s + 3<br />

s2 (s2 <br />

+ 4s + 3)<br />

<br />

1<br />

+ L<br />

s + 1<br />

−1<br />

<br />

2<br />

+ L<br />

s + 3<br />

−1<br />

= L −1<br />

<br />

1<br />

s2 <br />

= e −t + 2e −3t + t.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 67 / <strong>89</strong>


Konvolution (faltning)<br />

(f ∗ g)(t) =<br />

t<br />

0<br />

f (t − τ)g(τ) dτ<br />

L(f ∗ g) = L(f )L(g)<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 68 / <strong>89</strong>


Konvolution (faltning)<br />

Delta-funktionalen<br />

(f ∗ g)(t) =<br />

t<br />

0<br />

f (t − τ)g(τ) dτ<br />

L(f ∗ g) = L(f )L(g)<br />

δT = d<br />

u(t − T )<br />

dt<br />

men u(t − T ) är inte deriverbar s˚a δT är en generaliserad funktion, s˚a att<br />

∞<br />

f (t)δT (dt) = f (T ).<br />

−∞<br />

L(δT )(s) = e −sT , T ≥ 0<br />

(δT ∗ f )(t) = u(t − T )f (t − T ), T ≥ 0<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 68 / <strong>89</strong>


Exempel<br />

Beräkna den inversa transformationen L−1 <br />

1<br />

s2 konvolutionsprodukten.<br />

<br />

1<br />

s+2 med hjälp av<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 69 / <strong>89</strong>


Exempel<br />

Beräkna den inversa transformationen L−1 <br />

1<br />

s2 konvolutionsprodukten.<br />

Lösning: Eftersom L(t) = 1<br />

1<br />

s 2<br />

1<br />

s+2 = L(t ∗ e−2t )(s)<br />

s2 och L(e−2t ) = 1<br />

s+2<br />

<br />

1<br />

s+2 med hjälp av<br />

s˚a är<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 69 / <strong>89</strong>


Exempel<br />

Beräkna den inversa transformationen L−1 <br />

1<br />

s2 konvolutionsprodukten.<br />

Lösning: Eftersom L(t) = 1<br />

1<br />

s 2<br />

1<br />

s+2 = L(t ∗ e−2t )(s) och<br />

(t ∗ e −2t )(t) =<br />

=<br />

= 1<br />

2t −<br />

t<br />

0<br />

t 0<br />

t<br />

1<br />

s2 och L(e−2t ) = 1<br />

s+2<br />

(t − τ)e −2τ dτ<br />

(t − τ) − 1<br />

−2τ<br />

2 e −<br />

0<br />

2e−2τ dτ = 1<br />

2t −<br />

t<br />

0<br />

t 0<br />

<br />

1<br />

s+2 med hjälp av<br />

s˚a är<br />

(−1) − 1<br />

−2τ<br />

2 e dτ<br />

− 1<br />

4<br />

−2τ 1 1<br />

e = 2t + 4e−2t − 1<br />

4 .<br />

(Observera att i uttrycket (t ∗ e −2t )(t) betyder det första t:et funktionen t ↦→ t,<br />

följande i e −2t visar ocks˚a endast att det är fr˚aga om funktionen t ↦→ e −t medan<br />

det sista är ett argument som faktiskt kan ges ett reellt värde.)<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 69 / <strong>89</strong>


Exempel<br />

L˚at<br />

⎧<br />

⎪⎨ t, 0 ≤ t < 3,<br />

f (t) = 6 − t,<br />

⎪⎩<br />

0,<br />

3 ≤ t < 6,<br />

t > 6.<br />

Beräkna Laplace-transformationen av f genom att först skriva ut f med<br />

hjälp av stegfunktionen u och sedan använda förskjutningsregeln.<br />

Lösning: Ifall 0 ≤ a < b s˚a ser vi att<br />

<br />

0, d˚a t < a eller t > b,<br />

u(t − a) − u(t − b) =<br />

1, d˚a a < t < b.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 70 / <strong>89</strong>


forts.<br />

Vi kan nu skriva<br />

f (t) = t u(t) − u(t − 3) + (6 − t) u(t − 3) − u(t − 6) <br />

= tu(t) − 2(t − 3)u(t − 3) + (t − 6)u(t − 6).<br />

Med hjälp av förskjutningsregeln f˚ar vi nu (eftersom L(t)(s) = 1<br />

s 2 ) att<br />

L(f )(s) = 1 1<br />

− 2<br />

s2 s2 e−3s + 1<br />

s2 e−6s .<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 71 / <strong>89</strong>


Exempel<br />

Bestäm en lösning till diffusionsekvationen<br />

ut(t, x) = uxx(t, x), t > 0, x > 0,<br />

d˚a u(0, x) = 0, x > 0, u(t, 0) = g(t) och det finns en konstant a > 0 s˚a<br />

att ∞<br />

0 e−at |g(t)| dt < ∞.<br />

Lösning: Antag att u är en lösning till problemet som är s˚adan att<br />

räkneoperationerna som följer kan genomföras. L˚at<br />

U(s, x) = ∞<br />

0 e−stu(t, x) dx där s > a (dvs. vi använder ett reellt s). Om<br />

vi nu tar Laplace-transformationen av b˚ada sidorna i ekvationen ut = uxx<br />

s˚a f˚ar vi för all s > a, eftersom u(0, x) = 0 och vi antar att vi kan byta ordning<br />

p˚a integreringen och x-deriveringen,<br />

sU(s, x) = Uxx(s, x), x > 0.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 72 / <strong>89</strong>


forts.<br />

Denhär differentialekvationens karakteristiska ekvation är λ 2 − s = 0 s˚a att<br />

λ = ± √ s och den allmänna lösningen är U(s, x) = c1e √ sx + c2e −√ sx . Om<br />

nu U(s, x) för varje x är Laplace-transformationen av en funktion som inte<br />

växer för snabbt s˚a gäller lims→∞ U(s, x) = 0 vilket innebär att c1 = 0. Av<br />

antagandet följer att U(s, 0) = L(g)(s) vilket innebär att c2 = L(g)(s)<br />

och vi har<br />

U(s, x) = L(g)(s)e −x√ s .<br />

Genom att konsultera tabeller eller räkna komplexa integraler kan man<br />

konstatera, att<br />

<br />

1<br />

L<br />

2 √ <br />

1<br />

e− 4t (s) = e<br />

πt3 −√s ,<br />

vilket innebär att e −x√ s = e − √ x 2 s är Laplace-transformationen av<br />

funktionen<br />

1<br />

x 2<br />

1<br />

<br />

2 π( t<br />

x2 ) 3<br />

e<br />

− 1<br />

4 t<br />

x2 = x<br />

2 √ x2<br />

e−<br />

πt3 G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 73 / <strong>89</strong><br />

4t .


forts.<br />

Eftersom Laplace-transformationen av en konvolution är produkten av<br />

Laplace-transformationerna s˚a ser vi att lösningen är (och detta kan<br />

kontrolleras genom en direkt räkning)<br />

u(t, x) =<br />

t<br />

0<br />

x<br />

g(t − τ)<br />

2 √ πτ<br />

3 e− x2<br />

4τ dτ.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 74 / <strong>89</strong>


Fourier-integraler<br />

Om ∞<br />

−∞ |f (t)| dt < ∞ s˚a är Fourier-transformationen av f<br />

∞<br />

ˆf (ω) = e −i2πωt f (t) dt, ω ∈ R<br />

−∞<br />

Även definitionerna ∞<br />

−∞ e−iωtf (t) dt och 1<br />

∞<br />

√<br />

2π −∞ e−iωtf (t) dt används!<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 75 / <strong>89</strong>


Fourier-integraler<br />

Om ∞<br />

−∞ |f (t)| dt < ∞ s˚a är Fourier-transformationen av f<br />

∞<br />

ˆf (ω) = e −i2πωt f (t) dt, ω ∈ R<br />

−∞<br />

Även definitionerna ∞<br />

−∞ e−iωtf (t) dt och 1<br />

∞<br />

√<br />

2π −∞ e−iωtf (t) dt används!<br />

Lemma<br />

Ifall ∞<br />

−∞ |f (t)| dt < ∞ s˚a är ˆf kontinuerlig och lim |ω|→∞ ˆf (ω) = 0<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 75 / <strong>89</strong>


Räkneregler<br />

αf = αˆf<br />

f + g = ˆf + ˆg<br />

ˆg(ω) = e−i2πωt0ˆf (ω) ifall g(t) = f (t − t0)<br />

<br />

f (αt)(ω) = 1<br />

|α| ˆ f ω<br />

α<br />

e −πt2 (ω) = e−πω2 G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 76 / <strong>89</strong>


Teorem<br />

Ifall ∞<br />

−∞ |ˆf (ω)| dω < ∞ s˚a är<br />

∞<br />

f (t) = e i2πωtˆ f (ω) dω.<br />

−∞<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 77 / <strong>89</strong>


Teorem<br />

Ifall ∞<br />

−∞ |ˆf (ω)| dω < ∞ s˚a är<br />

Obs!<br />

∞<br />

f (t) = e i2πωtˆ f (ω) dω.<br />

−∞<br />

Ifall f (t) = 0 d˚a t < 0 s˚a är ˆf (ω) = L(f )(2πiω)<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 77 / <strong>89</strong>


Exempel<br />

Ge en formel med vars hjälp man kan invertera Laplace-transformationen.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 78 / <strong>89</strong>


Exempel<br />

Ge en formel med vars hjälp man kan invertera Laplace-transformationen.<br />

Lösning: L˚at F (s) = L(f )(s) och definiera f (t) = 0 d˚a t < 0. Det innebär<br />

att<br />

∞<br />

F (σ + i2πω) = e −(σ+i2πω)t ∞<br />

f (t) dt = e −i2πωt e −σt f (t) dt = ˆg(ω),<br />

om g(t) = e −σt f (t).<br />

0<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 78 / <strong>89</strong><br />

−∞


Exempel<br />

Ge en formel med vars hjälp man kan invertera Laplace-transformationen.<br />

Lösning: L˚at F (s) = L(f )(s) och definiera f (t) = 0 d˚a t < 0. Det innebär<br />

att<br />

∞<br />

F (σ + i2πω) = e −(σ+i2πω)t ∞<br />

f (t) dt = e −i2πωt e −σt f (t) dt = ˆg(ω),<br />

0<br />

|F (σ + i2πω)| dω < ∞<br />

s˚a f˚ar vi med hjälp av den inversa Fourier-transformationen<br />

g(t) = e −σt ∞<br />

f (t) = e i2πωt ∞<br />

ˆg(ω) dω = e i2πωt F (σ + i2πω) dω.<br />

−∞<br />

om g(t) = e−σtf (t). Om nu ∞<br />

−∞ |ˆg(ω)| dω = ∞<br />

−∞<br />

−∞<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 78 / <strong>89</strong><br />

−∞


Exempel<br />

Ge en formel med vars hjälp man kan invertera Laplace-transformationen.<br />

Lösning: L˚at F (s) = L(f )(s) och definiera f (t) = 0 d˚a t < 0. Det innebär<br />

att<br />

∞<br />

F (σ + i2πω) = e −(σ+i2πω)t ∞<br />

f (t) dt = e −i2πωt e −σt f (t) dt = ˆg(ω),<br />

0<br />

|F (σ + i2πω)| dω < ∞<br />

s˚a f˚ar vi med hjälp av den inversa Fourier-transformationen<br />

g(t) = e −σt ∞<br />

f (t) = e i2πωt ∞<br />

ˆg(ω) dω = e i2πωt F (σ + i2πω) dω.<br />

−∞<br />

om g(t) = e−σtf (t). Om nu ∞<br />

−∞ |ˆg(ω)| dω = ∞<br />

−∞<br />

−∞<br />

Genom att göra ett variabelbyte 2πω = τ kan vi skriva resultatet i formen<br />

f (t) = eσt<br />

∞<br />

e<br />

2π<br />

iτt F (σ + iτ) dτ.<br />

−∞<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 78 / <strong>89</strong><br />

−∞


Konvolution<br />

f ∗ g = ˆf ˆg när (f ∗ g)(t) = ∞<br />

−∞ f (t − τ)g(τ) dτ.<br />

Teorem<br />

Ifall ∞<br />

−∞ |f (t)|2 dt < ∞ s˚a finns det en funktion ˆf s˚a att<br />

∞<br />

|ˆf (ω)| 2 ∞<br />

dω = |f (t)| 2 dt<br />

−∞<br />

lim<br />

T →∞<br />

lim<br />

T →∞<br />

∞ <br />

<br />

<br />

−∞<br />

ˆ f (ω) −<br />

∞ <br />

<br />

f (t) −<br />

−∞<br />

−∞<br />

T<br />

−T<br />

T<br />

−T<br />

e −i2πωt 2<br />

<br />

f (t) dt<br />

dω = 0<br />

e i2πωtˆf<br />

2<br />

<br />

(ω) dω<br />

dt = 0<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 79 / <strong>89</strong>


Fourier-koefficienter<br />

Om 1<br />

0 |f (t)| dt < ∞ och f (t + 1) = f (t) s˚a är<br />

ˆf (n) =<br />

1<br />

0<br />

e −i2πnt f (t) dt<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 80 / <strong>89</strong>


Fourier-koefficienter<br />

Om 1<br />

0 |f (t)| dt < ∞ och f (t + 1) = f (t) s˚a är<br />

ˆf (n) =<br />

1<br />

0<br />

e −i2πnt f (t) dt<br />

Integralen kan ocks˚a räknas över andra intervall, t.ex.<br />

ˆf (n) = 1<br />

2<br />

− 1 e<br />

2<br />

−i2πntf (t) dt.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 80 / <strong>89</strong>


Fourier-koefficienter<br />

Om 1<br />

0 |f (t)| dt < ∞ och f (t + 1) = f (t) s˚a är<br />

ˆf (n) =<br />

1<br />

0<br />

e −i2πnt f (t) dt<br />

Integralen kan ocks˚a räknas över andra intervall, t.ex.<br />

ˆf (n) = 1<br />

2<br />

− 1 e<br />

2<br />

−i2πntf (t) dt.<br />

Om f istället har perioden T s˚a blir<br />

ˆf (n) = 1<br />

T<br />

i2πnt<br />

−<br />

e T f (t) dt =<br />

T 0<br />

1<br />

T<br />

T<br />

2<br />

− T<br />

2<br />

i2πnt<br />

−<br />

e T f (t) dt.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 80 / <strong>89</strong>


Fourier-koefficienter<br />

Om 1<br />

0 |f (t)| dt < ∞ och f (t + 1) = f (t) s˚a är<br />

ˆf (n) =<br />

1<br />

0<br />

e −i2πnt f (t) dt<br />

Integralen kan ocks˚a räknas över andra intervall, t.ex.<br />

ˆf (n) = 1<br />

2<br />

− 1 e<br />

2<br />

−i2πntf (t) dt.<br />

Om f istället har perioden T s˚a blir<br />

Lemma<br />

ˆf (n) = 1<br />

T<br />

i2πnt<br />

−<br />

e T f (t) dt =<br />

T 0<br />

1<br />

T<br />

1<br />

0<br />

T<br />

2<br />

− T<br />

2<br />

i2πnt<br />

−<br />

e T f (t) dt.<br />

|f (t)| dt < ∞ ⇒ lim ˆf (n) = 0<br />

|n|→∞<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 80 / <strong>89</strong>


Obs!<br />

1<br />

0<br />

e i2πnt dt =<br />

<br />

1, d˚a n = 0,<br />

0, d˚a n = 0.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 81 / <strong>89</strong>


Obs!<br />

Obs!<br />

Om<br />

s˚a är<br />

1<br />

0<br />

e i2πnt dt =<br />

f (t) = <br />

<br />

1, d˚a n = 0,<br />

n<br />

0, d˚a n = 0.<br />

ane i2πnt<br />

ˆf (n) = an, n ∈ Z.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 81 / <strong>89</strong>


Obs!<br />

Obs!<br />

Om<br />

s˚a är<br />

Teorem<br />

1<br />

0<br />

e i2πnt dt =<br />

f (t) = <br />

<br />

1, d˚a n = 0,<br />

n<br />

0, d˚a n = 0.<br />

ane i2πnt<br />

ˆf (n) = an, n ∈ Z.<br />

Om f (t + 1) = f (t), 1<br />

0 |f (t)| dt < ∞ och f är deriverbar i punkten t0 s˚a är<br />

lim<br />

N→∞<br />

M→−∞ n=M<br />

N<br />

e i2πnt0ˆf (n) = f (t0)<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 81 / <strong>89</strong>


Varför?<br />

Definiera<br />

g(t) =<br />

Av antagandet följer att 1<br />

t → 0.<br />

f (t+t0)−f (t0)<br />

e −i2πt −1 , t /∈ Z,<br />

i<br />

2π f ′ (t0), t ∈ Z.<br />

0 |g(t)| dt < ∞ eftersom e−i2πt −1<br />

t<br />

→ −i2π d˚a<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 82 / <strong>89</strong>


Varför?<br />

Definiera<br />

g(t) =<br />

f (t+t0)−f (t0)<br />

e −i2πt −1 , t /∈ Z,<br />

i<br />

2π f ′ (t0), t ∈ Z.<br />

Av antagandet följer att 1<br />

0 |g(t)| dt < ∞ eftersom e−i2πt −1<br />

t → −i2π d˚a<br />

t → 0. Eftersom f (t + t0) = g(t)(e−i2πt − 1) + f (t0) s˚a f˚ar vi<br />

ˆf (n)e i2πnt0 =<br />

1<br />

där δ0,0 = 1 och δn,0 = 0, n = 0.<br />

0<br />

e −i2πnt f (t + t0) dt = ˆg(n + 1) − ˆg(n) + f (t0)δn,0,<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 82 / <strong>89</strong>


Varför?<br />

Definiera<br />

g(t) =<br />

f (t+t0)−f (t0)<br />

e −i2πt −1 , t /∈ Z,<br />

i<br />

2π f ′ (t0), t ∈ Z.<br />

Av antagandet följer att 1<br />

0 |g(t)| dt < ∞ eftersom e−i2πt −1<br />

t → −i2π d˚a<br />

t → 0. Eftersom f (t + t0) = g(t)(e−i2πt − 1) + f (t0) s˚a f˚ar vi<br />

ˆf (n)e i2πnt0 =<br />

1<br />

0<br />

e −i2πnt f (t + t0) dt = ˆg(n + 1) − ˆg(n) + f (t0)δn,0,<br />

där δ0,0 = 1 och δn,0 = 0, n = 0. Av detta följer att (d˚a M < 0, N > 0)<br />

N<br />

ˆf (n)e i2πnt0 = ˆg(N + 1) − ˆg(M) + f (t0) → f (t0)<br />

n=M<br />

d˚a M → −∞ och N → ∞ eftersom lim |n|→∞ ˆg(n) = 0, vilket i sin tur är<br />

en följd av att 1<br />

0 |g(t)| dt < ∞.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 82 / <strong>89</strong>


Teorem<br />

Ifall 1<br />

0 |f (t)|2 dt < ∞ s˚a är<br />

1<br />

0<br />

|f (t)| 2 dt =<br />

1<br />

lim<br />

N→∞ 0<br />

M→−∞<br />

∞<br />

n=−∞<br />

<br />

<br />

<br />

f<br />

(t) −<br />

<br />

| ˆ f (n)| 2<br />

N<br />

n=−M<br />

e i2πntˆf<br />

<br />

<br />

<br />

(n) <br />

<br />

2<br />

dt = 0<br />

Omvänt gäller att ifall (cn) ∞ n=−∞ är s˚adan att ∞<br />

n=−∞ |cn| 2 < ∞ s˚a finns<br />

det en funktion f s˚a att 1<br />

0 |f (t)|2 dt < ∞ och ˆf (n) = cn.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 83 / <strong>89</strong>


1<br />

√ N<br />

Den diskreta Fourier-transformationen<br />

ˆF(m) =<br />

N−1 <br />

k=0<br />

i2πmk<br />

−<br />

e N F(k)<br />

N−1 i2πmk<br />

k=0<br />

e− N F(k) eller<br />

Ofta används instället defintionen 1<br />

N<br />

N−1 i2πmk<br />

k=0<br />

e− N F(k) men detta betyder bara att talet N dyker upp p˚a<br />

andra ställen.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 84 / <strong>89</strong>


1<br />

√ N<br />

Den diskreta Fourier-transformationen<br />

ˆF(m) =<br />

N−1 <br />

k=0<br />

i2πmk<br />

−<br />

e N F(k)<br />

N−1 i2πmk<br />

k=0<br />

e− N F(k) eller<br />

Ofta används instället defintionen 1<br />

N<br />

N−1 i2πmk<br />

k=0<br />

e− N F(k) men detta betyder bara att talet N dyker upp p˚a<br />

andra ställen.<br />

Obs!<br />

När F är periodisk, dvs. F(k + N) = F(k) gäller ocks˚a<br />

ˆF(m) =<br />

N+M−1 <br />

k=M<br />

i2πmk<br />

−<br />

e N F(k)<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 84 / <strong>89</strong>


FFT<br />

En algoritm som kan räkna den diskreta Fourier-transformationen med<br />

högst cN log(N) räkneoperationer och inte c ′ N 2 som en direkt räkning kräver.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 85 / <strong>89</strong>


FFT<br />

En algoritm som kan räkna den diskreta Fourier-transformationen med<br />

högst cN log(N) räkneoperationer och inte c ′ N 2 som en direkt räkning kräver.<br />

Den inversa transformationen<br />

F(k) = 1<br />

N−1 <br />

N<br />

m=0<br />

e i2πmk<br />

N<br />

ˆF(k)<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 85 / <strong>89</strong>


FFT<br />

En algoritm som kan räkna den diskreta Fourier-transformationen med<br />

högst cN log(N) räkneoperationer och inte c ′ N 2 som en direkt räkning kräver.<br />

Den inversa transformationen<br />

Konvolutioner<br />

F(k) = 1<br />

N−1 <br />

N<br />

m=0<br />

e i2πmk<br />

N<br />

ˆF(k)<br />

Om F(k + N) = F(k), G(k + N) = G(k) och H(k) = N−1 j=0 F(k − j)G(j)<br />

s˚a är<br />

ˆH(m) = ˆF(m)ˆG(m)<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 85 / <strong>89</strong>


Exempel<br />

Ifall talen ak och bk, k = 0, 1, . . . , är givna kan det i m˚anga fall finnas<br />

behov att räkna ut<br />

ck =<br />

k<br />

ak−jbj, k = 0, 1, . . . .<br />

j=0<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 86 / <strong>89</strong>


Exempel<br />

Ifall talen ak och bk, k = 0, 1, . . . , är givna kan det i m˚anga fall finnas<br />

behov att räkna ut<br />

ck =<br />

k<br />

ak−jbj, k = 0, 1, . . . .<br />

j=0<br />

Man kan först˚as räkna detta direkt men om det skall göras för m˚anga tal<br />

k kan det löna sig att använda den diskreta Fourier-transformationen p˚a<br />

följande sätt: Först m˚aste man bestämma sig för ett tal n s˚a att man<br />

räknar ut ck för k = 0, 1, 2, . . . , n − 1 och för detta behövs talen ak och bk<br />

d˚a k = 0, 1, . . . , n − 1.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 86 / <strong>89</strong>


Exempel<br />

Ifall talen ak och bk, k = 0, 1, . . . , är givna kan det i m˚anga fall finnas<br />

behov att räkna ut<br />

ck =<br />

k<br />

ak−jbj, k = 0, 1, . . . .<br />

j=0<br />

Man kan först˚as räkna detta direkt men om det skall göras för m˚anga tal<br />

k kan det löna sig att använda den diskreta Fourier-transformationen p˚a<br />

följande sätt: Först m˚aste man bestämma sig för ett tal n s˚a att man<br />

räknar ut ck för k = 0, 1, 2, . . . , n − 1 och för detta behövs talen ak och bk<br />

d˚a k = 0, 1, . . . , n − 1. Vi l˚ater N = 2n och definierar A(k) = ak och<br />

B(k) = bk d˚a k = 0, 1, . . . , n − 1, A(k) = B(k) = 0 d˚a k = n, . . . , N − 1<br />

och A(k + N) = A(k) och B(k + N) = B(k). Vi definierar<br />

C(k) =<br />

N−1 <br />

j=0<br />

A(k − j)B(j)<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 86 / <strong>89</strong>


forts.<br />

Eftersom B(j) = bj d˚a j = 0, . . . , n − 1 och B(j) = 0 d˚a j = n, . . . , N − 1<br />

s˚a är C(k) = n−1 j=0 A(k − j)bj. Dessutom är A(k − j) = ak−j d˚a<br />

0 ≤ k − j ≤ n − 1 s˚a att<br />

C(k) =<br />

k<br />

j=0<br />

ak−jbj +<br />

n−1<br />

j=k+1<br />

A(k − j)bj, k = 0, . . . , n − 1.<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 87 / <strong>89</strong>


forts.<br />

Eftersom B(j) = bj d˚a j = 0, . . . , n − 1 och B(j) = 0 d˚a j = n, . . . , N − 1<br />

s˚a är C(k) = n−1 j=0 A(k − j)bj. Dessutom är A(k − j) = ak−j d˚a<br />

0 ≤ k − j ≤ n − 1 s˚a att<br />

C(k) =<br />

k<br />

j=0<br />

ak−jbj +<br />

n−1<br />

j=k+1<br />

A(k − j)bj, k = 0, . . . , n − 1.<br />

D˚a j = k + 1, k + 2, . . . , n − 1 och k = 0, 1, . . . , n − 1 är<br />

k − j = −1, −2, . . . , k − n + 1 ≥ −n + 1 och därför är<br />

n + 1 = 2n − n + 1 = N − n + 1 ≤ k − j + N ≤ N − 1 s˚a att<br />

A(k − j) = A(k − j + N) = 0. Detta innebär att C(k) = ck d˚a<br />

k = 0, 1, . . . , n − 1. Eftersom Ĉ(m) = Â(m)ˆ B(m) s˚a följer det av formeln<br />

inverteringen av den diskreta Fourier-transformationen att<br />

ck = C(k) = 1<br />

N−1 <br />

e<br />

N<br />

i2πmk<br />

N Â(m)ˆB(m) = 1 <br />

ˆB(−k), k = 0, 1, . . . , n − 1.<br />

N<br />

m=0<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 87 / <strong>89</strong>


L˚at<br />

Numerisk beräkning av Fourier-integraler<br />

g(t) =<br />

N−1 <br />

k=0<br />

<br />

t − t0 − k∆t<br />

f (k∆t + t0)p<br />

,<br />

∆t<br />

där ∆t > 0 och p : R → R är s˚adan att p(0) = 1 och p(j) = 0 d˚a j = 0.<br />

D˚a är g är en interpolationsfunktion med g(k∆t + t0) = f (k∆t + t0) d˚a<br />

k = 0, . . . , N − 1 och g(k∆t) = 0 d˚a k < 0 eller k > N − 1. Om nu<br />

och<br />

s˚a är<br />

F(k) = f (k∆t + t0), k = 0, . . . , N − 1,<br />

∆t∆ω = 1<br />

N ,<br />

ˆg(m∆ω) = ∆te −i2πm∆ωt0 ˆF(m)ˆp<br />

<br />

m<br />

<br />

, m ∈ Z.<br />

N<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 88 / <strong>89</strong>


Varför?<br />

Vi väljer ∆ω = 1<br />

N∆t<br />

och d˚a blir<br />

∞<br />

ˆg(m∆ω) = e −i2πm∆ωt g(t) dt<br />

=<br />

=<br />

N−1 <br />

k=0<br />

N−1 <br />

k=0<br />

−∞<br />

∞<br />

F(k) e −i2πm∆ωt <br />

t − t0 − k∆t<br />

p<br />

∆t<br />

−∞<br />

∞<br />

F(k)∆t<br />

= ∆te −i2πm∆ωt0<br />

k=0<br />

dt t = t0 + k∆t + τ∆t<br />

=<br />

e<br />

−∞<br />

−i2πm∆ωt0 −i2πm∆ωk∆t −i2πm∆ω∆tτ<br />

e e p(τ) dτ<br />

N−1 <br />

∞<br />

i2πmk<br />

m<br />

− −i2π<br />

e N F (k) e N<br />

−∞<br />

τ p(τ) dτ<br />

= ∆te −i2πm∆ωt0 ˆF(m)ˆp<br />

<br />

m<br />

<br />

.<br />

N<br />

G. <strong>Gripenberg</strong> (<strong>TKK</strong>) <strong>Mat</strong>-<strong>1.1531</strong> <strong>Gk3</strong>-I <strong>20</strong> <strong>augusti</strong> <strong>20</strong>09 <strong>89</strong> / <strong>89</strong>

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!