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
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>