Introduzione Ai Metodi Matematici Della Fisica

Introduzione ai
METODI MATEMATICI DELLA FISICA

M. B. Barbaro, M. Frau, P. Gambino e S. Sciuto
Dipartimento di Fisica, Università di Torino
Dispense del corso 2013
Indice
I Funzioni Analitiche e Equazioni Differenziali in Campo
Complesso
4
1 Analisi Complessa
1.1 Il campo complesso . . . . . . . . . . . . . . . . . . . . . . . . . .
1.1.1 Richiami sui numeri complessi . . . . . . . . . . . . . . . .
1.1.2 Funzioni reali di variabile complessa . . . . . . . . . . . . .
1.2 Funzioni complesse di variabile complessa . . . . . . . . . . . . . .
1.2.1 Derivata di una funzione complessa di variabile complessa
1.2.2 Condizioni di Cauchy-Riemann . . . . . . . . . . . . . . .
1.2.3 Funzioni analitiche . . . . . . . . . . . . . . . . . . . . . .
1.2.4 Esempi . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.3 Integrazione in campo complesso . . . . . . . . . . . . . . . . . .
1.3.1 Curve (richiami) . . . . . . . . . . . . . . . . . . . . . . .
1.3.2 Integrali in campo complesso . . . . . . . . . . . . . . . . .
1.3.3 Teorema di Cauchy . . . . . . . . . . . . . . . . . . . . . .
1.3.4 Rappresentazione integrale di Cauchy . . . . . . . . . . . .
1.4 Serie in campo complesso . . . . . . . . . . . . . . . . . . . . . . .
1.4.1 Serie di potenze . . . . . . . . . . . . . . . . . . . . . . . .
1.4.2 Serie di Taylor . . . . . . . . . . . . . . . . . . . . . . . .
1.4.3 Zeri . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.4.4 Serie di Laurent . . . . . . . . . . . . . . . . . . . . . . . .
1.4.5 Singolarità isolate: poli e singolarità essenziali . . . . . . .
1.5 Residui . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.5.1 Teorema dei residui . . . . . . . . . . . . . . . . . . . . . .
1.5.2 Calcolo dei residui . . . . . . . . . . . . . . . . . . . . . .
1.6 Calcolo di integrali definiti mediante il teorema dei residui . . . .
1.6.1 Integrali trigonometrici . . . . . . . . . . . . . . . . . . . .
1.6.2 Esempi . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.6.3 Lemma di Jordan . . . . . . . . . . . . . . . . . . . . . . .
1.6.4 Esempi . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.6.5 Singolarità sul cammino di integrazione . . . . . . . . . . .
1.7 Studio del punto allinfinito . . . . . . . . . . . . . . . . . . . . .
1
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
5
5
5
9
10
12
12
15
15
19
19
20
24
32
35
35
36
38
40
41
43
44
45
47
47
47
50
52
59
61
1.8
1.7.1 Esempi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.7.2 Calcolo del residuo nel punto allinfinito . . . . . . . . . . . . .
Le funzioni ln z e z nel piano complesso . . . . . . . . . . . . . . . . .
2 Equazioni differenziali in C
2.1 Equazioni differenziali ordinarie II ordine . . . . . . . . . . .
2.1.1 Soluzione nellintorno di un punto regolare . . . . . . .
2.1.2 Soluzione nellintorno di un punto singolare fuchsiano
2.1.3 Esempio: lequazione di Bessel . . . . . . . . . . . . . .
2.1.4 Studio del comportamento allinfinito . . . . . . . . .
2.1.5 Esempi . . . . . . . . . . . . . . . . . . . . . . . . . . .
II
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
Introduzione allanalisi armonica
62
64
68
69
69
72
77
79
81
83
86
3 Serie di Fourier
3.1 Introduzione . . . . . . . . . . . . . . . . . . . . . . . . .
3.2 Funzioni periodiche e serie di Fourier . . . . . . . . . . .
3.2.1 Convergenza puntuale delle serie trigonometriche
di Fourier . . . . . . . . . . . . . . . . . . . . . .
3.2.2 Importanti commenti . . . . . . . . . . . . . . . .
3.2.3 Altri esempi . . . . . . . . . . . . . . . . . . . . .
. . . . . . . .
. . . . . . . .
87
87
92
. . . . . . . . 93
. . . . . . . . 96
. . . . . . . . 101
4 Trasformate Integrali
4.1 Trasformata di Fourier . . . . . . . . . . . . . . . . . .
4.1.1 Esempi . . . . . . . . . . . . . . . . . . . . . . .
4.1.2 Proprietà della trasformata di Fourier . . . . . .
4.1.3 Soluzione di equazioni differenziali mediante la
Fourier . . . . . . . . . . . . . . . . . . . . . . .
4.2 Trasformata di Laplace . . . . . . . . . . . . . . . . . .
4.2.1 Esempi . . . . . . . . . . . . . . . . . . . . . . .
4.2.2 Proprietà della trasformata di Laplace . . . .
4.2.3 Trasformate di Laplace ed equazioni differenziali
cienti costanti . . . . . . . . . . . . . . . . . . .
103
. . . . . . . . . 103
. . . . . . . . . 104
. . . . . . . . . 107
trasformata di
. . . . . . . . . 110
. . . . . . . . . 114
. . . . . . . . . 115
. . . . . . . . . 116
lineari a coeffi. . . . . . . . . 119
5 Spazi L2 e distribuzioni
5.1 Spazi L2 e serie di Fourier . . . . . . . . . . . .
5.1.1 Trasformata di Fourier in L2 . . . . . . . .
5.2 Alcuni sistemi ONC . . . . . . . . . . . . . . . . .
5.3 Sistemi ONC in L2 e distribuzioni. . . . . . . . .
5.3.1 Operatori autoaggiunti e sistemi ortogonali
5.3.2 La di Dirac e cenni sulle distribuzioni . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
124
124
131
134
137
137
139
A Funzioni armoniche
146
B Corollari della rappresentazione integrale di Cauchy
151
Equazioni differenziali del secondordine

156
C.1 Metodo di separazione delle variabili . . . . . . . . . . . . . . . . . . . 158
D Propriet`
a dellintegrale di Lebesgue
162
Parte I
Funzioni Analitiche e Equazioni
Differenziali in Campo Complesso
Capitolo 1
Analisi Complessa
1.1
Il campo complesso
1.1.1
Richiami sui numeri complessi
Un numero complesso z è una coppia ordinata di numeri reali

z = a + ib = (a, b)
a, b R
dove i è lunità immaginaria

i2 = 1 .
a = Re z ,
b = Im z .
Linsieme dei numeri complessi è un campo, indicato

N con C, dato dal prodotto cartesiano del campo reale R con se stesso (C = R R), dotato di due leggi di
composizione interna, laddizione e la moltiplicazione, che godono delle seguenti proprietà:
1) Addizione (+)
Definizione:
z1 = a1 + ib1 ,
z2 = a2 + ib2 z1 + z2 = (a1 + a2 ) + i(b1 + b2 ) .
Proprietà:
Associativa:
z1 + (z2 + z3 ) = (z1 + z2 ) + z3
Commutativa:
z1 + z2 = z2 + z1 .
5
z1 , z2 , z3 C .
Esiste lelemento neutro 0 C, tale che

z+0=0+z =z
z C .
Per ogni z C esiste lelemento inverso z C, tale che

z + (z) = (z) + z = 0 .
Quindi C è un gruppo abeliano rispetto alladdizione, con elemento neutro 0.
2) Moltiplicazione ()
Definizione:
z1 z2 = (a1 + ib1 ) (a2 + ib2 )
= (a1 a2 b1 b2 ) + i(a1 b2 + b1 a2 ) .
Proprietà:
Associativa:
z1 (z2 z3 ) = (z1 z2 ) z3 .
Commutativa:
z1 z2 = z2 z1 .
Esiste lelemento neutro 1 C, tale che
z1=1z =z
z C .
Per ogni z C, z 6= 0 esiste lelemento inverso z 1 C, tale che

z z 1 = z 1 z = 1
z = a + ib z 1 =
1
a
b
= 2
i 2
.
2
a + ib
a +b
a + b2
Quindi C {0} è un gruppo abeliano rispetto alla moltiplicazione, con elemento

neutro 1.
Vale inoltre la proprietà distributiva della moltiplicazione rispetto alladdizione
z1 (z2 + z3 ) = z1 z2 + z1 z3 .
Figura 1.1: Rappresentazione cartesiana del numero complesso z
Rappresentazione geometrica dei numeri complessi

Un numero complesso z si può rappresentare graficamente come un punto nel piano
complesso di z (o piano di Argand), sulle cui ascisse e ordinate si pongono rispettivamente la parte reale e immaginaria di z. La rappresentazione cartesiana di z è
(Fig 1.1)
z = x + iy .
Una rappresentazione equivalente è quella polare (Fig 1.2)
z = rei = r cos + ir sin
con
r = |z| =
p
x2 + y 2
modulo di z
e
= arg z = tan1
y
x
argomento o fase di z .
Graficamente r è il modulo del vettore ~r congiungente lorigine con il punto z, e

è langolo che questo vettore forma con lasse delle ascisse.
Le relazioni fra componenti cartesiane e polari di z sono:
x = r cos
y = r sin
7
Figura 1.2: Rappresentazione polare del numero complesso z
Il complesso coniugato di un numero complesso z è un numero complesso z C

cos` definito:
z = x + iy = rei z = x iy = rei .
(si veda Fig 1.3) Si noti che z z = |z|2 = r2 .
Figura 1.3: Rappresentazione polare del numero complesso z e del suo complesso
coniugato z
1.1.2
Funzioni reali di variabile complessa
Una funzione reale di variabile complessa

f : E R ,
EC
è unapplicazione che associa un numero reale f (z) ad ogni z E, con E sottoinsieme

del campo C.
Definizione di funzione continua.
La funzione f (z) si dice continua nel punto z = z0 se
lim f (z) = f (z0 )
zz0
ossia, ricordando la definizione di limite, se

> 0 > 0 / |f (z) f (z0 )| < ,
z I (z0 ) ,
dove I (z0 ) è un intorno di raggio del punto z0 :

I (z0 ) = {z C/|z z0 | < } .
Esempi
1) La funzione modulo
f (z) = |z|
è una funzione reale di variabile complessa, continua in tutto il piano complesso.
2) Le funzioni
f (z) = Rez
e g(z) = Imz
sono funzioni reali di variabile complessa, continue in tutto il piano complesso.

3) La funzione argomento
(z) = arg z
è una funzione reale di variabile complessa:
C {0} I2 R ,
9
Figura 1.4: Discontinuità dellargomento di z.
dove I2 è un intervallo semiaperto di lunghezza 2. Tale intervallo non è univocamente definito e può essere scelto in infiniti modi diversi, ma in ogni caso
la funzione (z) è discontinua su una semiretta uscente dallorigine del piano
complesso. Si considerino per esempio i due casi rappresentati in Fig. 1.4:
a) I2 = (, ]
b) I2 = [0, 2)
Nel caso a) (z) è discontinua sul semiasse reale negativo. Infatti
(x) =
x R+
ma il limite limzx (z) non è definito, perche i limiti destro e sinistro sono
diversi:
lim (x + i) =
0+
lim (x i) = .
0+
Nel caso b) invece (z) è discontinua sul semiasse reale positivo.
1.2
Funzioni complesse di variabile complessa
Una funzione complessa di variabile complessa

f : E C ,
10
EC
è unapplicazione che associa un numero complesso f (z) ad ogni z E, con E

sottoinsieme del campo C. Useremo la seguente notazione:
f:
z 7 w = f (z) z E , E C , w C ,
z = x + iy
w(z) = u(x, y) + iv(x, y) .
Quindi dare la f è equivalente a specificare due funzioni reali di due variabili reali:
u = u(x, y) e v = v(x, y) .
Analogamente a quanto accade nel caso di funzioni reali, una funzione complessa
di variabile complessa f (z) è detta continua se
lim f (z) = f (z0 )
zz0
ovvero se
> 0 > 0 / |f (z) f (z0 )| < ,
z I (z0 ) ,
dove I (z0 ) è un intorno di raggio del punto z0 e |f (z) f (z0 )| è il modulo del numero
complesso f (z) f (z0 ).
Esempi
1)
f (z) = z 2 = (x + iy)2 = x2 y 2 + 2ixy

u(x, y) = x2 y 2
v(x, y) = 2xy .
In coordinate polari
z = rei
w = r2 e2i .
2)
f (z) = z = x iy
u(x, y) = x
v(x, y) = y .
In coordinate polari
z = rei
w = rei .
11
1.2.1
Derivata di una funzione complessa di variabile complessa
Definizione: una funzione f (z) si dice derivabile nel punto z se esiste il limite per
h 0 (h C) del rapporto incrementale
f (z + h) f (z)
h
considerato come funzione della variabile complessa h. Tale limite deve essere quindi
indipendente dal modo in cui h 0. Esistono infatti infinite direzioni lungo le quali h
può tendere a 0 (si veda Fig 1.5)
Figura 1.5: Direzioni dellincremento h nel rapporto incrementale
La funzione f è derivabile se tutte queste direzioni danno lo stesso risultato per

il limite del rapporto incrementale. In questo caso il limite si chiama derivata di f
rispetto a z:
df (z)
f (z + h) f (z)
= f 0 (z) = lim
.
h0
dz
h
1.2.2
Condizioni di Cauchy-Riemann
Consideriamo una funzione

f (z) = u(x, y) + iv(x, y)
12
tale che nel punto z = x + iy sia la sua parte reale u(x, y) che la sua parte immaginaria
v(x, y) siano di classe C 1 , cioè continue con le loro derivate prime:
u(x, y)
v(x, y)
v(x, y)
u(x, y)
= u0x ,
= u0y ,
= vx0 ,
= vy0 .
x
y
x
y
Teorema 1: le condizioni di Cauchy e Riemann (CR)
u(x, y)
v(x, y)
=
x
y
u(x, y)
v(x, y)
=
y
x
(1.1)
sono condizioni necessarie e sufficienti affinche la funzione f (z) sia derivabile nel punto
z.
Dimostrazione.
Dimostriamo dapprima che le condizioni di (CR) sono necessarie, ossia che
f (z) derivabile (CR) .
Per ipotesi la derivata di f (z)
f (z + h) f (z)
h0
h
f 0 (z) = lim
esiste ed è indipendente dalla direzione di h = hx + ihy . In particolare si potrà scegliere

h puramente reale (h = hx ) o puramente immaginario (h = ihy ).
Se h = hx
f (z + hx ) f (z)
hx 0
hx
u(x + hx , y) + iv(x + hx , y) u(x, y) iv(x, y)
= lim
hx 0
hx
u(x + hx , y) u(x, y)
v(x + hx , y) v(x, y)
= lim
+ i lim
hx 0
hx 0
hx
hx
u(x, y)
v(x, y)
=
+i
.
x
x
f 0 (z) =
lim
13
(1.2)
Se h = ihy
f (z + ihy ) f (z)
hy 0
ihy
u(x, y + hy ) + iv(x, y + hy ) u(x, y) iv(x, y)
= lim
hy 0
ihy
u(x, y + hy ) u(x, y)
v(x, y + hy ) v(x, y)
= lim
+ i lim
hy 0
hy 0
ihy
ihy
u(x, y) v(x, y)
+
.
(1.3)
= i
y
y
Uguagliando ora le parti reali e immaginarie delle espressioni (1.2) e (1.3) per la derivata
f 0 (z) otteniamo le condizioni di Cauchy-Riemann:
f 0 (z) =
lim
u0x = vy0
u0y = vx0 .
Dimostriamo ora che le condizioni di CR sono sufficienti per la derivabilità di f (z),
ossia
(CR) f (z) derivabile .
Consideriamo a questo scopo il rapporto incrementale
f (z + h) f (z)
u(x + hx , y + hy ) + iv(x + hx , y + hy ) u(x, y) iv(x, y)
=
.
h
hx + ihy
(1.4)
Poiche le funzioni u e v sono per ipotesi continue con le loro derivate prime in z, esse
sono differenziabili e si può quindi scrivere nellintorno del punto (x, y):
u(x + hx , y + hy ) = u(x, y) + hx u0x (x, y) + hy u0y (x, y) + o(|h|)
v(x + hx , y + hy ) = v(x, y) + hx vx0 (x, y) + hy vy0 (x, y) + o(|h|) .
Sostituendo questi sviluppi nel rapporto incrementale (1.4) si ottiene
hx u0x (x, y) + ihx vx0 (x, y) + hy u0y (x, y) + ihy vy0 (x, y) + o(|h|)
f (z + h) f (z)
=
.
h
hx + ihy
Utilizzando ora le condizioni di Cauchy-Riemann (1.1) e prendendo il limite per hx , hy
0 si ha
f (z + h) f (z)
lim
=
hx ,hy 0
h
hx u0x (x, y) + ihx vx0 (x, y) hy vx0 (x, y) + ihy u0x (x, y) + o(|h|)
lim
=
hx ,hy 0
hx + ihy
(hx + ihy )[u0x (x, y) + ivx0 (x, y)]
lim
= u0x (x, y) + ivx0 (x, y) = f 0 (z) .
hx ,hy 0
hx + ihy
(1.5)
14
La derivata di f (z) è quindi definita univocamente indipendentemente dalla direzione

di h: la funzione è pertanto derivabile e la sua derivata è
f 0 (z) = u0x (x, y) + ivx0 (x, y) .
[q.e.d.]
Utilizzando le condizioni di Cauchy-Riemann è possibile dare quattro espressioni
equivalenti della derivata di una funzione in termini delle sue parti reale e immaginaria:
f 0 (z) =
=
=
=
u0x (x, y) + ivx0 (x, y)

vy0 (x, y) iu0y (x, y)
u0x (x, y) iu0y (x, y)
vy0 (x, y) + ivx0 (x, y) .
(1.6)
N.B. Dalle ultime due espressioni si deduce che per calcolare la derivata di f (z) è
sufficiente conoscerne o la parte reale u o la parte immaginaria v. Ovvero: nota una
delle due, si può ricavare laltra a meno di una costante.
1.2.3
Funzioni analitiche
Definizione: Una funzione f (z)

f : F C
F C
si dice analitica (o regolare, o olomorfa) in un punto z0 se esiste un intorno I(z0 )

in cui f (z) è derivabile, con derivata continua, in ogni punto z I(z0 ). (N.B. non è
sufficiente che le condizioni di CR siano soddisfatte solo nel punto z = z0 .) La funzione
è analitica in una regione aperta E F se essa è derivabile, con derivata continua, in
ogni punto z E.
` quindi necessario e sufficiente affinche f (z) sia analitica in E che siano soddisfatte
E
le seguenti condizioni in tutti i punti di E:
1) parte reale u(x, y) e parte immaginaria v(x, y) siano di classe C 1 ;
2) siano soddisfatte le condizioni di Cauchy-Riemann (1.1).
I punti in cui f (z) è analitica si dicono punti di analiticità o punti regolari della
funzione. I punti in cui f (z) non è analitica si dicono punti singolari o singolarit`
a
della funzione.
1.2.4
Esempi
Esempio 1
f (z) = costante = c = cx + icy
15
c C, cx , cy R
u(x, y) = cx ,
v(x, y) = cy
u0x = u0y = vx0 = vy0 = 0

La funzione f (z) è continua in tutto il piano complesso C, le funzioni u e v sono
continue e derivabili e le condizioni di Cauchy-Riemann sono verificate z C. Quindi
una funzione costante è analitica in tutto il piano complesso e la sua derivata è zero:
f 0 (z) = u0x + ivx0 = 0
Esempio 2
f (z) = z = x iy
u(x, y) = x ,
u0x = 1 ,
v(x, y) = y
vy0 = 1
La funzione f (z) = z non è analitica. Si può infatti mostrare che il rapporto incrementale dipende dalla direzione dellincremento h. Sia h = ei in rappresentazione
polare. Allora il rapporto incrementale
(z + h) z
h
ei
f (z + h) f (z)
=
=
=
= e2i
h
h
h
ei
dipende dallangolo . La derivata di z non esiste in alcun punto di C.
Esempio 3
f (z) =
u(x, y) =
1
1
x iy
=
= 2
z
x + iy
x + y2
x
,
x2 + y 2
v(x, y) =
y
x2 + y 2
Le funzioni u e v sono continue e derivabili in C {0}. Le condizioni di CR

u0x =
u0y =
y 2 x2
,
(x2 + y 2 )2
2xy
,
+ y 2 )2
(x2
vy0 =
vx0 =
16
y 2 x2
= u0x
(x2 + y 2 )2
(x2
2xy
= u0y
2
2
+y )
sono soddisfatte. La funzione f (z) è quindi analitica in C {0}. La sua derivata è

f 0 (z) = u0x + ivx0 =
y 2 x2 + 2ixy
(y + ix)2
1
=
= 2 .
2
2
2
2
2
(x + y )
(y + ix) (y ix)
z
Vale pertanto la stessa regola di derivazione valida in campo reale

0
1
1
= 2 .
z
z
Si dimostra, esattamente come nel campo reale, che vale
dz n
= nz n1
dz
per n intero qualsiasi 1 .
Le funzioni analitiche in tutto il piano complesso, come, per esempio, i polinomi, la
funzione esponenziale e le funzioni seno e coseno, si chiamano funzioni intere.
Teorema 2: se f1 (z) e f2 (z) sono due funzioni analitiche nel punto z, allora le
funzioni
1) f1 (z) + f2 (z)
2) f1 (z) f2 (z)
3) f1 (z)/f2 (z) se f2 (z) 6= 0
4) f1 (f2 (z))
sono analitiche in z e valgono le seguenti regole di derivazione
a) [f1 (z) + f2 (z)]0 = f10 (z) + f20 (z)
b) [f1 (z) f2 (z)]0 = f10 (z)f2 (z) + f1 (z)f20 (z)
c) [f1 (z)/f2 (z)]0 = [f10 (z)f2 (z) f1 (z)f20 (z)]/[f2 (z)]2
d)
d
[f (f (z))]
dz 1 2
df1 df2
df2 dz
La dimostrazione segue banalmente dalla definizione di derivata o dalla (1.6).

1
La stessa formula vale per ogni esponente, reale o complesso, ma non ne parliamo qui perche non
abbiamo ancora definito z per non intero.
17
Corollario: le funzioni razionali di z sono analitiche in tutto il piano complesso

esclusi gli zeri del denominatore.
Dimostrazione: Poiche le funzioni f1 (z) = 1 e f2 (z) = z sono analitiche in C,
segue dalla proprietà 2) che tutte le potenze di z sono analitiche in C, e quindi per la
proprietà 1) i polinomi di z
Pn (z) =
n
X
ck z k = c0 + c1 z + c2 z 2 + ...cn z n
k=0
sono funzioni ovunque analitiche. Per la proprietà 3) le funzioni razionali (rapporto di

due polinomi Pn e Qm )
R(z) =
Pn (z)
Qm (z)
sono funzioni analitiche in tutto il piano complesso esclusi i punti zi tali che Qm (zi ) = 0.
[q.e.d.]
Teorema 3 : se la parte reale (immaginaria) di una funzione analitica è costante,
anche la sua parte immaginaria (reale) è necessariamente costante. Infatti da u(x, y) =
costante segue u0x = u0y = 0 e quindi dalle condizioni di CR segue che:
vy0 = vx0 = 0 v(x, y) = K 0 = costante .
Pertanto
f (z) = costante .
Nel caso particolare in cui v(x, y) = 0, oppure u(x, y) = 0, ne segue banalmente che
una funzione analitica a valori reali (o immaginari puri) è necessariamente costante.
[q.e.d.]
Teorema 4: una funzione analitica di modulo costante è costante (cioè le sue parti
reale e immaginaria sono separatamente costanti):
f (z) analitica e |f (z)| = cost. f (z) = cost.
Dimostrazione
Per ipotesi
u2 (x, y) + v 2 (x, y) = K .
Se K = 0 la dimostrazione è banale perche ciò implica f (z) = 0; assumiamo quindi nel
seguito K 6= 0.
18
Derivando rispetto a x e a y si ottiene

2ux u + 2vx v = 0
2uy u + 2vy v = 0 .
Moltiplicando la prima equazione per u, la seconda per v, sommando membro a membro
e utilizzando le condizioni di CR, si ottiene
(u2 + v 2 )x u = 0
da cui segue che, poiche u2 + v 2 è per ipotesi costante e diverso da zero, x u = 0.
Analogamente, moltiplicando la prima equazione per v e la seconda per u, si ottiene
(u2 + v 2 )y u = 0.
Questa implica che anche y u = 0 e quindi u(x, y) = costante. Dalle CR segue immediatamente che se le derivate parziali di u sono nulle, anche le derivate parziali di v
sono nulle, e pertanto f (z) = costante, come nel teorema precedente.
[q.e.d.]
1.3
1.3.1
Integrazione in campo complesso

Curve (richiami)
Una curva nel piano complesso è una applicazione continua

: J C
J = [a, b] R
dove J è un intervallo reale limitato e chiuso:

: t z(t) = x(t) + iy(t)
atb.
Lapplicazione associa ad ogni valore del parametro t due funzioni reali x(t) e y(t).
Spesso si considera la curva non solo come lapplicazione appena definita, ma come
limmagine (o sostegno) di tale applicazione, cioè come linsieme di punti
= {z C/z = z(t), t [a, b]} .
Una curva si dice regolare nellintervallo [a, b] se le funzioni x(t) e y(t) hanno
derivate prime continue e non entrambe nulle t [a, b].
Una curva si dice regolare a tratti nellintervallo [a, b] se lintervallo può essere
suddiviso in un numero finito di sottointervalli chiusi in cui la curva sia regolare.
Una curva si dice chiusa se z(a) = z(b). Un caso particolare di curva chiusa è un
punto, cioè una curva di equazione z(t) = costante t [a, b].
19
Una curva si dice semplice se z(t1 ) 6= z(t2 ) t1 6= t2 , con t1 , t2 [a, b). (N.B.
Lintervallo [a, b) è semi-aperto per includere le curve chiuse nella definizione di curve
semplici.) In pratica una curva semplice è una curva che non si interseca con se stessa.
Una curva chiusa e semplice si dice curva di Jordan.
Enunciamo, senza dimostrarlo, il seguente teorema:
Teorema 5: ogni curva di Jordan divide il piano in due regioni, una interna e
una esterna a .
Ad ogni curva chiusa si assegna un verso di percorrenza. Convenzionalmente si
considera come positivo il verso antiorario. Si definisce convenzionalmente interna ad
una curva chiusa la zona lasciata a sinistra se si percorre la curva nel suo verso di
percorrenza (ed esterna quella lasciata a destra).
Due curve di Jordan 1 e 2 si dicono omotopicamente equivalenti (O.E.) in
una regione D se possono essere deformate con continuità luna nellaltra senza uscire
` essenziale specificare la regione D in cui le due curve sono O.E.
da D. N.B. E
Esempio: se D = C ogni curva di Jordan è O.E. a un punto, ma questo non è pi`
u
vero se da C si sottraggono uno o pi`
u punti.
Una regione D C si dice connessa per archi se, z1 , z2 D, esiste una curva
tutta interna a D, che congiunge z1 e z2 .
Una regione S C si dice semplicemente connessa (s.c.) se ogni curva chiusa
contenuta is S è O.E. a un punto. (Definizione alternativa: una regione S è s.c. se
per ogni curva di Jordan contenuta in S la regione interna a è sottoinsieme di S).
Intuitivamente una regione s.c. è una regione senza buchi.
Lemma di Gauss (o teorema di Green): siano P (x, y), Q(x, y) C 1 due funzioni
reali e continue con derivate prime continue in un dominio E semplicemente connesso.
Allora per ogni curva chiusa regolare a tratti contenuta in E
ZZ
I
[P (x, y)dx + Q(x, y)dy] =

Q(x, y) P (x, y)
dxdy ,
x
y
(1.7)
dove S è la regione interna a .
1.3.2
Integrali in campo complesso
Integrale di una funzione di variabile reale a valori complessi Consideriamo

una funzione complessa di una variabile reale w(t) = u(t) + iv(t):
w : [a, b] C
[a, b] R
t [a, b] , w(t) C .
20
Definiamo lintegrale di w(t) in t

Z b
Z b
Z b
w(t)dt =
u(t)dt + i
v(t)dt .
a
Lintegrale esiste se la funzione w(t) è continua o se ha un numero finito di discontinuità

di prima specie.
Lintegrale (alla Riemann) si può interpretare come limite di somme integrali:
Z b
w(t)dt = lim In
n
dove
In =
n
X
w(l )(tl tl1 )
l=1
n
X
u(l )(tl tl1 ) + i
l=1
n
X
v(l )(tl tl1 ) .
l=1
Si divide cioè lintervallo [a, b] in n sottointervalli a = t0 < t1 < t2 < ... < tn = b e si
valuta la funzione w(t) nei punti l interni a ciascun sottointervallo (tl1 < l < tl ).
Dalla disuguaglianza triangolare (|a1 + a2 + ... + an | |a1 | + |a2 | + ... + |an |):
|In |
n
X
|w(l )|(tl tl1 )
l=1
segue, prendendo il limite per n , la relazione

Z b
Z b

w(t)dt
|w(t)|dt .

a
(1.8)
Integrale di una funzione complessa di variabile complessa Sia f (z) una

funzione:
f : D C C
f : z = (x + iy) D 7 w = f (z) = u(x, y) + iv(x, y) C .
Si consideri una curva regolare a tratti nellintervallo [a, b]
: t 7 z(t)
atb
e siano A = z(a) e B = z(b) gli estremi di tale curva (Fig. 1.6)

21
Figura 1.6: Curva aperta che unisce i punti A e B nel piano complesso
Se f (z) è continua z , si definisce lintegrale curvilineo di f (z) tra A e B lungo
A()
dz
dt .
dt
f (z(t))
f (z)dz =
a
(1.9)
N.B. Il secondo membro della (1.9) esiste perche è regolare a tratti (quindi dz/dt
ha un numero finito di discontinuità).
Teorema 6: lintegrale di una funzione continua f (z) = u(x, y) + iv(x, y) lungo
una curva regolare a tratti è dato da
Z
[u(x, y)dx v(x, y)dy] + i
f (z)dz =
A()
A()
[v(x, y)dx + u(x, y)dy] .

A()
(1.10)
Dimostrazione
Dalla definizione (1.9) segue che
22
f (z)dz =
A()
=
+
=

Z b
dz(t)
dx(t)
dy(t)
f (z(t))
dt =
[u(x, y) + iv(x, y)]
+i
dt
dt
dt
dt
a
a

Z b
dy(t)
dx(t)
v(x, y)
dt
u(x, y)
dt
dt
a

Z b
dx(t)
dy(t)
i
v(x, y)
+ u(x, y)
dt
dt
dt
a
Z B
Z B
[v(x, y)dx + u(x, y)dy] .
Z
A()
A()
[q.e.d]
N.B. In generale lintegrale dipende dalla curva e non solo dagli estremi di integrazione.
Valgono per lintegrale (1.10) le proprietà degli integrali curvilinei. In particolare,
se C è un punto sulla curva ,
Z B
Z C
Z B
f (z)dz
f (z)dz +
f (z)dz =
C()
A()
A()
e
Z
f (z)dz =
f (z)dz .
B()
A()
Teorema 7: Disuguaglianza di Darboux Sia M il valore massimo assunto dal

modulo della funzione f (z) lungo la curva :
M = max |f (z)|
z
e l la lunghezza di tra A e B:
Z
l=
a
ds
dt
dt
Z
a
s
dx
dt
2

+
Allora vale la disuguaglianza di Darboux:

Z B

f (z)dz M l .

dy
dt
2
dt .
(1.11)
A()
Dimostrazione: Applicando la disuguaglianza (1.8) alla definizione (1.9) si ottiene

Z B

Z b

dz

dt
f
(z)dz
|f
(z(t))|
dt

A()
a
Z b
dz
dt .
M
dt
a
23
Ora,
s 2 2
dz
dx
dy
ds
=
;
+
dt
dt
dt
dt
pertanto
Z b
dz
dt = l
dt
a
e quindi:
Z

A()

f (z)dz M l .
[q.e.d.]
Consideriamo qualche esempio di integrale in campo complesso. Sia C la circonferenza di raggio R attorno allorigine. Vogliamo calcolare lintegrale della funzione
analitica f (z) = z
I
I=
z dz .
(1.12)
Sostituiamo quindi z = Rei e facciamo riferimento alla (1.9):

Z 2
2
I = iR
e2i d = 0.
(1.13)
Nel caso della funzione non analitica f (z) = z troviamo invece

I
Z 2
2
z dz = iR
d = 2iR2 .
C
1.3.3
(1.14)
Teorema di Cauchy
Teorema 8: sia f (z) una funzione regolare allinterno di un dominio aperto E

semplicemente connesso. Il teorema di Cauchy asserisce che, per ogni curva chiusa,
regolare a tratti, tutta contenuta in E,
I
f (z)dz = 0 .
(1.15)
Lintegrale di una funzione analitica è nullo lungo una qualsiasi curva chiusa omotopicamente equivalente a un punto nel dominio di analiticità della funzione.
24
Figura 1.7: Curva chiusa che non contiene lorigine
Dimostrazione
Dal teorema (1.10) si ha che
I
f (z)dz =
[v(x, y)dx + u(x, y)dy]
e, per il lemma di Gauss (1.7) (si ponga P = u, Q = v nel primo integrale e Q = u,

P = v nel secondo)

ZZ
I
v(x, y) u(x, y)
f (z)dz =
dxdy
x
y
S

ZZ
u(x, y) v(x, y)
dxdy ,
+ i
x
y
S
dove S E è la regione interna a . Poiche f (z) è analitica valgono le condizioni di
Cauchy-Riemann u0x = vy0 e u0y = vx0 . Pertanto
I
f (z)dz = 0 .
[q.e.d.]
In altre parole, f (z)dz = 0 se è contenuta nel dominio E di analiticità di f (z)
ed è deformabile con continuità in un punto senza uscire da E.
In modo pi`
u conciso si può anche dire che la forma differenziale f (z)dz = u(x, y)dx
v(x, y)dy + i[v(x, y)dx + u(x, y)dy] è chiusa in un aperto E: d(f (z)dz) = 0, se valgono
le condizioni di CR; essa diventa esatta se il dominio è semplicemente connesso.
H
25
Figura 1.8: Curve aperte che uniscono i punti A e B in un dominio semplicemente

connesso
Corollario al teorema di Cauchy

Teorema 9: siano 1 e 2 due curve semplici e regolari a tratti che congiungono i
punti A e B e = 1 (2 ) sia tutta contenuta nel dominio semplicemente connesso
di analiticità di f (z) (Fig. 1.8). Allora
Z
f (z)dz
f (z)dz =
A(2 )
A(1 )
ovvero: lintegrale di una funzione analitica non dipende dal cammino di integrazione
purche i cammini siano deformabili con continuità luno nellaltro senza incontrare
singolarità.
Dimostrazione
Z
f (z)dz
A(1 )
Z
Z
+
f (z)dz =
A(2 )
A(1 )
I
f (z)dz =
B(2 )
f (z)dz = 0
per il teorema di Cauchy (1.15).

[q.e.d.]
Esempio: consideriamo lintegrale
I
I=
dz
.
z
La funzione 1/z è analitica in C {0}. Se la regione interna alla curva non contiene
lorigine (Fig. 1.7) lintegrale è nullo per il teorema di Cauchy.
Se invece lorigine è interna a , per esempio è una circonferenza C di raggio R centrata in O (Fig. 1.9) lintegrale è diverso da zero. Calcoliamone il valore. Lequazione
26
Figura 1.9: Curva chiusa che contiene lorigine
della curva C in coordinate polari è

z = z() = Rei = R(cos + i sin )
dz
= R( sin + i cos ) = iz
d
I
I=
C
dz
=
z
Z
0
1 0
z ()d =
z()
0 2
Z
0
iRei
d = 2i
Rei
N.B. Lintegrale non dipende da R.

Ne segue che, se 1 e 2 sono due semicirconferenze centrate nellorigine (Fig. 1.10)
gli integrali
Z B
dz
I1 =
A(1 ) z
e
Z
I2 =
A(2 )
dz
z
non devono necessariamente essere uguali poiche non si può applicare il Corollario del
teorema di Cauchy. Infatti essi valgono
Z
Z
I1 = i
d = i , I2 = i
d = i .
0
27
Figura 1.10: Semicirconferenze centrate nellorigine
N.B. I1 6= I2 perche deformando 1 in 2 si attraversa una singolarità (z = 0).

In modo analogo si può calcolare il seguente integrale
I
I=
C
dz
za
aC
(1.16)
dove la curva C è la circonferenza di raggio R centrata in a (Fig. 1.11).

Infatti, ponendo
z = z() = a + Rei
si ha
z 0 () = iRei = i(z a)
e quindi
Z
I=
0
1
z 0 ()d =
z() a
Z
0
iRei
d = 2i .
Rei
Teorema 10: Teorema di Cauchy generalizzato Sia f (z) una funzione analitica in un dominio D qualsiasi e siano 1 e 2 due curve chiuse omotopicamente
equivalenti in D. In queste ipotesi:
I
I
f (z)dz =
f (z)dz .
1
28
Figura 1.11: Circonferenza centrata nel punto a
Dimostrazione: Per dimostrare il teorema consideriamo 3 casi:

a) D semplicemente connesso
b) D generico, 1 e 2 non si intersechino
c) D generico, 1 e 2 si intersechino
Caso a) In questo caso la dimostrazione è banale perche, per il teorema di Cauchy,
I
I
f (z)dz =
f (z)dz = 0 .
1
Caso b) Effettuiamo due tagli AB e CD (Fig. 1.12). Poiche 1 e 2 sono O.E. in D, la

regione compresa tra le due curve appartiene tutta a D. Si ha allora (per il corollario
al teorema di Cauchy):
Z
=
Z
A(E)
A
A
C
B(F )
B
D(G)
Z
Z
C
A
+
C(H)
.
B
Sommando membro a membro si ottiene

I
I
f (z)dz =
f (z)dz
1
poiche
Z
=
A
e
B
29
.
C
Figura 1.12: Curve 1 e 2 che non si intersecano in un dominio D
Caso c) In questo caso si ha (vedi Fig. 1.13):
Figura 1.13: Curve 1 e 2 che si intersecano in un dominio D
=
Z
A(E)
A
A(F )
A
=
B(G)
.
B(H)
Sommando membro a membro si ottiene

I
I
f (z)dz =
f (z)dz .
1
30
N.B. Non è detto che la regione (AGBF A) appartenga a D (1 e 2 sono O.E. in

D).
[q.e.d.]
Corollario: lintegrale (1.16) vale 2i per ogni curva chiusa che circondi il
punto a:
I
dz
= 2i,
aC
(1.17)
I=
z a
Inoltre, sempre per ogni curva chiusa che circondi il punto a,
I
dz
= 2i n,1 ,
(z a)n
dove nl è la delta di Kronecker definita da

1
se n = l
nl =
0
se n 6= l
nZ
(1.18)
n, l Z .
La (1.18) si dimostra osservando che il cammino può essere deformato in una circonferenza C di raggio R e centro a e ponendo
z a Rei dz = iRei d ;
(1.19)
ne segue immediatamente, per n 6= 1:

I
dz
=
(z a)n
Z
0
i
iRei
d = n1
n
in
R e
R
i(1n)
d =
i
Rn1
2
ei(1n)
= 0,
i(1 n) 0
mentre per n = 1 vale la eq.(1.17)

Abbiamo visto che lintegrale di f (z) tra i punti A e B non dipende dalla curva
dintegrazione se essa rimane allinterno del dominio di analticità D di f e se D è
semplicemente connesso. In questi casi possiamo allora definire una primitiva di f (z)
attraverso lintegrale
Z
z
f (z 0 ) dz 0 .
F (z) =
(1.20)
z0
Si dimostra che F (z) è unica e analitica in D e che F 0 (z) = f (z). Di conseguenza

valgono anche
Z b
f (z) dz = F (b) F (a)
(1.21)
a
31
e le usuali regole di calcolo integrale. Per esempio lintegrale indefinito di z è z 2 /2,

quello di ez è ez , quello di 1/z è ln z. Naturalmente in questultimo caso il dominio
di analiticità C {0} non è semplicemente connesso, e lintegrale su ogni curva chiusa
che contiene lorigine vale 2i. Tuttavia la (1.21) vale ancora, in virt`
u della polidromia
del logaritmo (vedi sezione 1.8): dopo un giro attorno allorigine il logaritmo risulta incrementato di 2i e si dice che siamo passati su un altro ramo della funzione
logaritmo.
La (1.18) ci dice che talvolta lintegrale su un cammino chiuso si annullaHanche se
questo racchiude una singolarità; questo accade per n > 1. Lannullarsi di f (z)dz
infatti non garantisce che la funzione sia analitica allinterno
di . Il teorema di
H
Morera ci assicura invece che se f (z) è continua e f (z) dz = 0, chiusa in D
connesso, allora f (z) è analitica in D.
1.3.4
Rappresentazione integrale di Cauchy
Teorema 11: sia f (z) una funzione analitica in un dominio E aperto semplicemente connesso e una curva di Jordan contenuta in E. Sia S la regione interna
a . Sotto queste ipotesi vale la rappresentazione integrale di Cauchy per la
funzione f (z):
1
f (z) =
2i
f (z 0 ) 0
dz
z0 z
z S , = S .
(1.22)
N.B. Perche valga la (1.22) è essenziale che z appartenga a S, cioè sia interna a .
Infatti,
1) se z la funzione integranda ha una singolarità sul cammino di integrazione e
quindi lintegrale non esiste;
2) se z è esterno a , la funzione integranda f (z 0 )/(z 0 z) è analitica in S e lintegrale
è nullo per il teorema di Cauchy (1.15) ;
3) se z è interno a , la funzione integranda f (z 0 )/(z 0 z) in generale non è analitica
in S, ma può avere una singolarità in z 0 = z.
Dimostrazione
Consideriamo il seguente integrale
I
f (z 0 ) f (z) 0
I(z) =
dz .
z0 z
32
Per il teorema generalizzato di Cauchy può essere deformata in una circonferenza C

centrata in z di raggio arbitrario r, interna a S:
I
f (z 0 ) f (z) 0
I(z) =
dz .
z0 z
C
Per la disuguaglianza di Darboux (1.11),

f (z 0 ) f (z)

2r = 2 max |f (z 0 ) f (z)| ,
|I(z)| max
z 0 C
z 0 C
z0 z
dove si è usato |z 0 z| = r, z 0 C. Poiche la funzione f (z) è continua in S, il secondo
membro può essere reso arbitrariamente piccolo. Infatti, dalla definizione di continuità
di una funzione, limz0 z f (z 0 ) = f (z), segue
> 0 / |f (z 0 ) f (z)| < z 0 I (z) .
Basta allora scegliere r < per dimostrare che > 0 si ha
|I(z)| 2.
Ma lunico numero non negativo minore o uguale a un numero positivo arbitrario è lo
zero, quindi
I(z) = 0,
ossia, utilizzando lintegrale (1.17),
I
I
f (z) 0
f (z 0 ) 0
dz
=
dz = f (z) 2i ,
0
0
z z
z z
da cui la (1.22).
[q.e.d.]
La rappresentazione integrale di Cauchy permette quindi di conoscere i valori di
una funzione analitica in tutta la regione interna ad una curva chiusa una volta noti
i suoi valori nei punti appartenenti alla curva .
Rappresentazione integrale di Cauchy per le derivate
Teorema 12: sia f (z) una funzione analitica in un dominio E aperto semplicemente connesso e una curva di Jordan contenuta in E. Sia S la regione interna
a . In queste ipotesi vale la rappresentazione integrale di Cauchy (1.22); derivando
ripetutamente sotto il segno di integrale e utilizzando lidentità:
dn 1
n!
,
= 0
n
0
dz z z
(z z)n+1
33
si ottiene (in modo non rigoroso) la rappresentazione integrale di Cauchy per le derivate
della funzione f (z):
dn f (z)
n!
=
n
dz
2i
f (z 0 )
dz 0
(z 0 z)n+1
z S .
(1.23)
Unimportantissima conseguenza della (1.23) è che una funzione analitica `

e infinitamente derivabile, ovviamente con derivate continue. Vale la pena sottolineare
che nulla di simile accade per le funzioni reali di variabili reali: una funzione f (x) differenziabile non è necessariamente infinitamente differenziabile (si pensi a f (x) = x|x|,
con derivata 2|x| che non è differenziabile nellorigine).
34
1.4
Serie in campo complesso
1.4.1
Serie di potenze
Una serie di potenze è una serie del tipo
ak (z z0 )k .
k=0
Per le serie di potenze in campo complesso valgono teoremi analoghi a quelli validi in
campo reale:
La regione di convergenza di una serie di potenze in C è un cerchio (centrato
in z0 ), il cui raggio si dice raggio di convergenza della serie. Allinterno di tale
cerchio la serie è uniformemente e assolutamente convergente. Allesterno non
converge mai. Sulla circonferenza può convergere o no, a seconda dei casi, e cè
sempre almeno un punto di non convergenza.
Teorema di Weierstrass: una serie di potenze è, per ogni z interno al cerchio
di convergenza, derivabile termine a termine n volte (con n arbitrario):
dn f (z) X dn (z z0 )k
ak
=
,
ak (z z0 )
f (z) =
n
n
dz
dz
k=0
k=0
k
quindi essa è analitica allinterno del cerchio di convergenza.

Teorema
il raggio di convergenza della serie di poP di Cauchy-Hadamard:
k
tenze k=0 ak (z z0 ) coincide con linverso del massimo fra i punti di accumulazione della successione {|ak |1/k }, ovvero, se il limite esiste, con:
=
lim |ak |1/k
o1
Si può mostrare che la (1.24) è equivalente alla pi`

u comoda:

an
,
= lim
n an+1
sempre che il limite, finito o infinito, esista.
35
(1.24)
(1.25)
1.4.2
Serie di Taylor
Lo sviluppo in serie di Taylor è uno sviluppo in serie di potenze di una funzione nellintorno di un suo punto di analiticità. Sia f (z) una funzione analitica in un dominio
D e sia C un intorno circolare, tutto contenuto in D, di un punto regolare z0 (si veda
` facile dimostrare che la funzione f (z), infinitamente derivabile 2 , può
Fig. 1.14). E
essere rappresentata, per ogni z C, dalla serie di Taylor:
f (z) =
ak (z z0 )k
(1.26)
k=0
con
ak

1 dk f (z)
=
k!
dz k z=z0
I
f (z)
1
dz ,
=
2i c (z z0 )k+1
(1.27)
dove per lultimo passaggio si è usata la rappresentazione di Cauchy (1.23) delle derivate
di una funzione analitica.
Dimostrazione: partendo dalla rappresentazione integrale di Cauchy
I
1
f (z 0 )
dz 0 ,
f (z) =
2i c z 0 z
usando (vedi Fig. 1.14)
n

1
1
1 X z z0
= 0
= 0
z0 z
z z0 (z z0 )
z z0 n=0 z 0 z0
(1.28)
e integrando termine
a termine la serie (uniformemente convergente, poiche |z z0 | <

zz0
0
|z z0 |, quindi z0 z0 < 1) ) si ottiene
I
X
f (z 0 )
n 1
f (z) =
(z z0 )
dz 0 ,
0 z )n+1
2i
(z
0
c
n=0
2
(1.29)
Nel campo reale linfinita derivabilit`

a non è sufficiente per garantire che una funzione sia
sviluppabile in serie di Taylor; nel campo complesso basta invece la analiticità.
36
Figura 1.14: Intorno circolare del punto z0 nel dominio D
che coincide con la (1.26).

[q.e.d]
Poiche il cerchio C deve essere tutto interno al dominio di analiticità di f (z), il suo
raggio r deve essere minore della distanza di z0 dal pi`
u vicino punto singolare di f (z):
il cerchio C è certamente contenuto nel cerchio di convergenza della serie di Taylor
(1.26).
` molto importante sottolineare che esiste una completa equivalenza tra analiticità
E
di una funzione in un punto e sua sviluppabilità in serie di Taylor in un suo intorno:
f (z) analitica in z0
I(z0 ) / z I(z0 ) , f (z) =
ak (z z0 )k .
k=0
Esempio: la serie geometrica
k=0
z k converge uniformemente a
f (z) =
1
1z
nella regione |z| < 1. Il punto z = 1 è infatti un punto singolare di f (z).

Talvolta conosciamo una funzione solo attraverso una serie di potenze, che converge
nel cerchio S1 , ma potrebbe avere un dominio di analiticità pi`
u ampio. Lestensione
tramite serie di Taylor al di fuori delloriginale dominio S1 viene detta continuazione
analitica (alla Weierstrass) della funzione di partenza. Nel caso in Fig.1.15 abbiamo
una serie di potenze attorno allorigine con raggio di convergenza determinato dalla
distanza dalla singolarità z1 . Lespansione di Taylor attorno a z2 S1 converge nel
cerchio S2 , che si estende al di là del dominio S1 (abbiamo qui assunto che z1 sia la
37
Figura 1.15: Continuazione analitica
singolarità pi`
u vicina a z2 , ma potrebbe non essere cos`). In S1 S2 le due serie sono
rappresentazioni diverse della stessa funzione e coincidono. Si dimostra poi che due
funzioni analitiche che coincidono su un insieme continuo di punti, coincidono ovunque
siano entrambe ben definite e rappresentano la stessa funzione. Il procedimento può
essere ripetuto un numero arbitrario di volte, anche in presenza di altre singolarità, fino
a raggiungere qualsiasi punto del dominio di analiticit`
a connesso a quello di partenza.
P
Nellesempio in figura la prima serie è f1 (z) = n=0 (z)n attorno allorigine con
raggio di convergenza 1. Questa serie può essere derivata in z2 interno a S1 ottenendo
un nuovo sviluppo di Taylor attorno a z2 che converge in S2 . Trattandosi di serie
geometriche, sappiamo che in entrambi i casi la somma è f (z) = 1/(1 + z), analitica
in C {0}. Si può pensare a questultima rappresentazione come una continuazione
analitica delle prime due serie.
Riassumendo, la conoscenza di una funzione analitica e delle sue derivate in un
unico punto di analiticità permette, in linea di principio, di ricostruire la funzione in
tutto il suo dominio di analiticità.
1.4.3
Zeri
Un punto regolare z = z0 è uno zero di ordine n della funzione f (z) se:
38
1) La funzione si annulla in z0 :
f (z0 ) = 0
2) Le prime n 1 derivate si annullano in z0 :

dk f (z)
=0,
k = 1, 2, ... , n 1
dz k z=z0
3) La derivata n-esima è diversa da zero in z0 :

dn f (z)
6= 0 .
dz n z=z0
Per esempio la funzione f (z) = z 2 ha uno zero di ordine 2 in z = 0. Infatti:
f 0 (0) = 0 ,
f (0) = 0 ,
f 00 (0) = 2 6= 0 .
Uno zero è un punto regolare di f (z), che sarà quindi rappresentabile tramite uno
sviluppo in serie di Taylor intorno a quel punto:
f (z) =
ak (z z0 )k .
(1.30)
k=0
Se z0 è uno zero di ordine n di f (z), si ha, dalla (1.27), che

a1 = a2 = ... = an1 = 0
an 6= 0 .
Quindi la (1.30) diventa

f (z) =
ak (z z0 )k .
k=n
Cambiando lindice nella sommatoria, k k 0 = k n, si ottiene

f (z) =
ak0 +n (z z0 )k +n
k0 =0
n
= (z z0 )
X
k0 =0
La funzione
g(z) =
ak+n (z z0 )k
k=0
39
ak0 +n (z z0 )k .
è una funzione analitica in z0 , in quanto sviluppabile in serie di Taylor intorno al punto

z = z0 . Inoltre g(z) è non nulla in z0 :
g(z0 ) = an 6= 0 .
Pertanto una funzione f (z) che abbia in z0 uno zero di ordine n può sempre essere
scritta nella forma
f (z) = (z z0 )n g(z) ,
(1.31)
con g(z) analitica e non nulla in z = z0 .
1.4.4
Serie di Laurent
Fra le varie possibili singolarità di una funzione analitica giocano un ruolo particolarmente importante le singolarità isolate: un punto singolare z0 si dice singolarit`
a
isolata della funzione f (z) se esiste un suo intorno privato di z0 in cui f (z) è analitica.
Nellintorno di una singolarità isolata è necessario considerare anche serie a potenze
negative; per esempio, per ogni z 6= 0 vale:
e
1/z
X
z k
k!
k=0
(1.32)
come si vede subito ponendo w = 1/z; evidentemente lorigine è una singolarità isolata
per la funzione e1/z .
Lo sviluppo in serie di Laurent è uno sviluppo in serie di potenze positive e negative
di una funzione f (z) in un intorno bucato I(z0 ) di un suo punto singolare isolato z0
(si veda Fig. 1.16); Il teorema di Laurent dice che se esiste un I(z0 ) in cui f (z) è
analitica, per ogni z I(z0 ), si può scrivere3 :
f (z) =
dk (z z0 )k .
(1.33)
k=
Ricordando la (1.18) è immediato calcolare i coefficienti dk ; basta dividere la (1.33)

per (z z0 )n+1 e integrare su una curva chiusa interna a I(z0 ) e che circondi z0 per
ottenere:
3
La dimostrazione è analoga a quella dello sviluppo in serie di Taylor, con la differenza che la curva
da scegliere per la rappresentazione integrale di Cauchy di f (z) è composta dalle due circonferenze
(percorse in verso opposto) che delimitano una corona circolare centrata in z0 e contenuta in I(z0 ).
40
Figura 1.16: Curva in un intorno bucato del punto z0
1
dk =
2i
f (z)
dz ,
(z z0 )k+1
(1.34)
N.B. Anche per k > 0, il coefficiente dk in generale non è la derivata dk f (z)/dz k

perche f (z) non è analitica nel punto z0 .
Notare inoltre che lo sviluppo in serie di Taylor (1.26) si ottiene come caso particolare dello sviluppo in serie di Laurent (1.33) se invece z0 è un punto regolare di f (z).
Infatti, se z0 è regolare, per il teorema di Cauchy dk = 0 per tutti i k 1 e la (1.33)
si riduce alla (1.26).
1.4.5
Singolarit`
a isolate: poli e singolarit`
a essenziali
Il punto singolare isolato z0 si definisce polo della funzione f (z) se lo sviluppo in serie
di Laurent intorno a z0 possiede un numero finito n di potenze negative; un polo si
dice di ordine n se:
1) I coefficienti ..., dn2 , dn1 si annullano:
dk = 0 ,
k > n
2) Il coefficiente dn è diverso da zero:

dn 6= 0 .
41
Un polo di ordine 1 si dice polo semplice, di ordine 2 polo doppio e cos` via.
Nellintorno di un polo di ordine n lo sviluppo (1.33) si riduce quindi a
f (z) =
dk (z z0 )k
(1.35)
k=n
e, mediante il cambiamento di indice k k 0 = k + n, a

f (z) =
dk0 n (z z0 )k n
k0 =0
= (z z0 )
dkn (z z0 )k .
k=0
Definiamo ora
g(z) =
dkn (z z0 )k .
k=0
La funzione g(z) è data da uno sviluppo in serie di Taylor intorno a z0 : essa è pertanto
analitica in z0 . Inoltre g(z) è diversa da zero in z0 :
g(z0 ) = dn 6= 0 .
Pertanto se z0 è un polo di ordine n della funzione f (z), questa può essere espressa
come
f (z) =
g(z)
,
(z z0 )n
(1.36)
con g(z) analitica e non nulla in z = z0 . Dalle (1.36) e (1.31) segue che il punto z = z0
è uno zero di ordine n per la funzione 1/f (z):
1
= h(z)(z z0 )n ,
f (z)
dove h(z) = 1/g(z) è di nuovo una funzione analitica e non nulla in z0 .
Singolarit`
a essenziali
Il punto z = z0 si definisce invece singolarit`
a essenziale isolata della funzione f (z)
se è un punto singolare isolato e lo sviluppo in serie di Laurent intorno a z0 possiede
un numero infinito di potenze negative.
Come si vede dallo sviluppo in serie di Laurent (1.32), un esempio di singolarità
essenziale isolata è lorigine per la funzione e1/z , e analogamente per le funzioni sin(1/z),
cos(1/z) e simili.
42
Dalle definizioni di polo e singolarità essenziale isolata segue che un polo di ordine
n può essere rimosso moltiplicando la f (z) per (z z0 )n , mentre questo non è possibile
per una singolarità essenziale.
Unaltra importante differenza fra poli e singolarità essenziali è la seguente: è evidente dalla (1.36) che limzz0 f (z) = se il punto z0 è un polo di f (z); se invece z0
è una singolarità essenziale il limite non esiste, perche nellintorno di z0 la funzione
oscilla forsennatamente: per farsene unidea, basta pensare allandamento nellintorno
dellorigine della funzione sin z1 con z reale o immaginario puro.
Pi`
u in generale, si può dimostrare il Teorema di Weierstrass per le singolarit`
a
essenziali isolate: se z = z0 è una singolarità essenziale isolata della funzione f (z),
allora per ogni e piccoli a piacere e per ogni numero complesso c C, esiste un
valore di z I (z0 ) tale che
|f (z) c| < .
In altre parole, il teorema di Weierstrass afferma che in qualunque intorno di una
singolarità essenziale isolata, la funzione f (z) approssima indefinitamente qualunque
valore prefissato c, senza necessariamente raggiungerlo.
1.5
Residui
Sia f (z) una funzione analitica in un dominio D, z0 un punto singolare isolato, una
curva di Jordan, tutta contenuta in D e contenente al suo interno il punto z0 , ma non
altre singolarità (questo è possibile, perche z0 è isolato). Si definisce residuo della
funzione f (z) nel punto z = z0 la quantità
{Resf (z)}z=z0
2i
I
f (z)dz .
(1.37)
Dalla eq.(1.34) che definisce i coefficienti di Laurent, calcolata per k = 1, si vede

subito che vale:
{Resf (z)}z=z0 = d1 .
(1.38)
Quindi il residuo di una funzione in un punto singolare isolato z0 è il coefficiente della

potenza (1) esima del suo sviluppo in serie di Laurent intorno a z0 .
43
Esempio
1
1
{Resf (z)}z=0 =
f (z) =
z
2i
dz
=1,
z
dove è una curva che circonda lorigine.

Ovviamente, se z0 è un punto regolare di f (z) e cerchiamo ugualmente di calcolare
il residuo, troviamo zero per il teorema di Cauchy. Non vale però il viceversa: una
funzione può avere residuo nullo in un punto ed ivi essere singolare, se d1 = 0 ma esiste
qualche dn 6= 0 per n > 1; un esempio caratteristico è f (z) = 1/z 2 che nellorigine ha
un polo doppio, con residuo nullo.
1.5.1
Teorema dei residui
Teorema 13: sia f (z) una funzione analitica in un dominio D, eccetto che in un
numero finito di singolarità isolate. Sia una curva di Jordan contenuta in D, non
passante per alcun punto singolare di f (z). In queste ipotesi vale il teorema dei residui:
I
f (z)dz = 2i
n
X
{Resf (z)}z=zk ,
(1.39)
k=1
dove z1 , z2 , ..., zn sono le singolarità di f (z) interne a .

Dimostrazione
Il teorema dei residui si dimostra facilmente per induzione completa. Infatti la
(1.39) è vera per n=1 per la definizione di residuo. Se le singolarità sono n + 1 isoliamo
la (n + 1)-esima come in Fig. 1.17.
Usando la tesi (1.39) per n singolarità si ottiene
I
I
I
f (z)dz
f (z)dz =
f (z)dz +
= 2i {Resf (z)}z=zn+1 + 2i
n
X
{Resf (z)}z=zk
k=1
= 2i
n+1
X
{Resf (z)}z=zk .
(1.40)
k=1
Quindi la (1.39) è vera per n+1 singolarità.

` utile osservare che il numero di singolarità interne alla curva deve essere finito;
E
se fosse infinito, allinterno di ci sarebbe un punto di accumulazione di singolarità e
quindi una singolarità non isolata, per cui non ha senso definire il residuo.
44
Figura 1.17: Curva che contiene n + 1 singolarità della funzione
1.5.2
Calcolo dei residui
Vediamo ora come si calcola esplicitamente il residuo di una funzione in un suo punto
singolare isolato. Se z0 è una singolarit`
a essenziale non cè altro modo4 che usare le
equazioni (1.37) e (1.38).
Se invece z = z0 è un polo di ordine n di f (z) cè un modo alternativo che richiede
solo di calcolare derivate. Infatti in un intorno di z0 si può scrivere:
f (z) =
g(z)
(z z0 )n
con g(z) analitica e non nulla in z0 . Il residuo è, dalla (1.37),

I
g(z)
1
dz .
{Resf (z)}z=z0 =
2i (z z0 )n
Ora, dalla rappresentazione di Cauchy per le derivate di f (z) (1.23)
k

I
d g(z)
k!
g(z)
=
dz ,
k
dz
2i (z z0 )k+1
z=z0
si ottiene, ponendo k = n 1,
1
2i
g(z)
1
dz =
n
(z z0 )
(n 1)!
dn1 g(z)
dz n1

z=z0
locale, perche vedremo pi`

u avanti che il discorso può essere diverso se si conosce il comportamento
globale della funzione in tutto il piano complesso, punto allinfinito compreso.
45
e quindi, poiche
g(z) = (z z0 )n f (z) ,
{Resf (z)}z=z0
1
=
lim
(n 1)! zz0

dn1
n
[(z z0 ) f (z)] .
dz n1
(1.41)
Nel caso particolare in cui z0 sia un polo semplice (n = 1), si ha

{Resf (z)}z=z0 = lim [(z z0 )f (z)] .
zz0
Ricordando lo sviluppo in serie di Laurent nel caso di un polo semplice (1.35)
f (z) =
dk (z z0 )k
k=1
d1
+ d0 + d1 (z z0 ) + ... ,
z z0
si ottiene
(
{Resf (z)}z=z0 = lim
zz0
d1 +
)
dk (z z0 )k+1
k=0
e quindi
{Resf (z)}z=z0 = d1 ,
(1.42)
a conferma di quanto visto in generale in eq.(1.38).

Esempio: la funzione f (z) = z sin z/(z )3 ha un polo doppio in z = , infatti la
sua serie di Laurent nellintorno di z = è
P
2k+1

k (z)
z
1
1
z sin z
k=0 (1)
(2k+1)!
=
= ( + z )
+ ...
f (z) =
(z )3
(z )3
(z )2 6
2
1
+ ...
(z )2 z
6
da cui segue d1 = 1, che coincide col residuo sul polo secondo la (1.38). Alternativamente definiamo h(z) = (z )3 f (z) e applichiamo la (1.41). Il residuo è
1 00
1
h (z)|z= = (2 cos z z sin z)|z= = 1 .
2
2
46
Si sarà notato che abbiamo applicato la (1.41) come se la funzione avesse un polo triplo.
Infatti la derivazione della (1.41) è valida anche se n è maggiore dellordine del polo.
Nel caso in questione, il calcolo del residuo si può fare usando n = 3 o n = 2, ma è pi`
u
diretto con n = 3.
1.6
Calcolo di integrali definiti mediante il teorema

dei residui
Il teorema dei residui (1.39) è di grande utilità perche permette non solo di calcolare
integrali naturalmente definiti su curve chiuse nel piano complesso, ma anche ampie
classi di integrali definiti sullasse reale, trasformandoli in integrali in campo complesso.
1.6.1
Integrali trigonometrici
Una prima classe da considerare è la seguente5 :

Z 2
f (cos , sin )d .
(1.43)
Per calcolare integrali di questo tipo

a) si usano le formule di Eulero per esprimere sin e cos in funzione di ei (senza
usare la complessa coniugazione!) e si effettua la sostituzione
z = ei d = i
dz
;
z
(1.44)
b) ci si riconduce ad un integrale nel piano z lungo una circonferenza di raggio

unitario;
c) si calcola lintegrale con il teorema dei residui.
1.6.2
Esempi
Esempio 1
Z
I=
0
d
.
5 + 3 cos
Con la sostituzione (1.44) si ha:

ei + ei
1
=
cos =
2
2
5

1
z+
.
z
Naturalmente lintervallo dintegrazione potrebbe essere anche (, ) o qualsiasi altro intervallo

di ampiezza 2.
47
Sostituendo in I:
I
1
dz
C z 5 + 3/2 (z + 1/z)
I
2i
dz
=
2
3 C z + 10
z+1
3
I = i
dove C è una circonferenza di raggio unitario centrata nellorigine. Studiamo ora le

singolarità della funzione integranda:
f (z) =
z2 +
10
z+1=0
3
1
z2 +
10
z
3
+1
1
, z2 = 3
3

10
1
2
z + z+1= z+
(z + 3) .
3
3
z1 =
La funzione f (z) ha due poli semplici in z = 1/3 (interno alla curva C e z = 3

(esterno alla curva C). Pertanto

2i
1
I = 2i Res 2 10
3
z + 3 z + 1 z=1/3
4
z + 1/3

lim
3 z1/3 z + 13 (z + 3)
=
.
2
=
Esempio 2
Z
I=
0
d
,
1 2p cos + p2
Poniamo
z = ei d = i
dz
z
Allora (vedi esempio precedente)

1
cos =
2

1
z+
.
z
48
pC
Sostituendo in I:
I
dz
1
2
C z 1 p (z + 1/z) + p
dz
2
2
C pz (1 + p )z + p
I = i
I
= i
dove C è una circonferenza di raggio unitario centrata nellorigine. La funzione integranda

f (z) =
pz 2
1
(1 + p2 )z + p
ha due poli semplici:

pz 2 (1 + p2 )z + p = 0
z1 = 1/p , z2 = p
pz 2 (1 + p2 )z + p = p (z 1/p) (z p)
e quindi
i
I=
p
I
C
dz
.
(z 1/p) (z p)
Dove sono situati i poli di f (z)?

se
se
|p| < 1
|p| > 1
z = p interno a C , z = 1/p esterno a C

z = p esterno a C , z = 1/p interno a C .
Ne segue che, se |p| < 1,

i
1
I =
2i Res
p
(z p)(z 1/p) z=p

2
zp
= lim
p zp (z p)(z 1/p)
2
=
1 p2
e, se |p| > 1,

i
1
I =
2i Res
p
(z p)(z 1/p) z=1/p

2
z 1/p
=
lim
p z1/p (z p)(z 1/p)
2
.
= 2
p 1
49
Se |p| = 1, lintegrando ha una singolarità sul cammino di integrazione e I non è

definito. Si noti che lesempio 1 è un caso particolare dellesempio 2.
1.6.3
Lemma di Jordan
Molto spesso possono essere calcolati con il metodo dei residui anche integrali estesi a
tutto lasse reale:
Z
g(x) dx ,
(1.45)
I=
dove supporremo che g(x) non abbia singolarità su R. In tal caso la strategia da seguire
è di considerare accanto allintegrale I lintegrale:
Z
Z R
g(x)dx +
g(z)dz,
(1.46)
J(R) =
R
dove R è una semicirconferenza, centrata nellorigine e di raggio R, situata nel semipiano Im z > 0 o Im z < 0 a seconda dei casi (vedi Fig. 1.18).
Figura 1.18: Semicirconferenze di raggio r nel semipiano inferiore (a) e superiore (b)
Lintegrale J(R) è esteso a una curva chiusa e si può quindi calcolare con il metodo
dei residui; dalla conoscenza di J(R) è poi immediato calcolare lintegrale I se la
funzione g(z) è tale che:
Z
lim
g(z)dz = 0 .
(1.47)
R
50
Ciò succede nei casi seguenti6 :

1) La funzione g(z) tende a zero pi`
u velocemente di 1/|z| per z :

g(z) = o
1
|z|

z.
(1.50)
In questo caso la semicirconferenza R può giacere sia nel semipiano Im z > 0 sia
nel semipiano Im z < 0.
Dimostrazione.
Passando a coordinate polari (e supponendo di considerare la semicirconferenza
nel semipiano superiore)
Z
Z

g(z)dz = iR
g Rei ei d
R
e applicando la disuguaglianza (1.8) si ottiene

Z

Z

g Rei d

g(z)dz
R

0
da cui, sfruttando lipotesi

lim R g Rei = 0 ,
si ottiene infine la (1.47). (Notare che non si è fatto altro che ridimostrare la
disuguaglianza di Darboux (1.11) in questo caso particolare.)
[q.e.d.]
2) La funzione integranda g(z) è della forma eiz f (z), con > 0, dove f (z) è una
funzione che tende uniformemente (rispetto allargomento di z) a zero quando |z|
6
In realt`
a dalla (1.47) segue che
Z
lim J(R) = lim
g(x)dx,
(1.48)
che coincide con

Z
I=
R2
lim
R1 ,R2
g(x)dx,
(1.49)
R1
solo nel caso che questultimo integrale esista; se ciò non succede, ma esiste il limite (1.48), allora
questo si chiama valore principale dellintegrale (1.45), si veda anche Sez. 1.6.5
51
tende a infinito e largomento di z è compreso fra 0 e (cioè nel semipiano Im

z 0), ovvero 7
z ,
f (z) = o(1),
0 arg z .
(1.51)
In tal caso, se si sceglie per R una semicirconferenza nel semipiano superiore

(Im z > 0), centrata nellorigine e di raggio R, è facile dimostrare che vale:
Z
lim
eiz f (z)dz = 0 .
(1.52)
3) Con la sostituzione z z, si vede subito che la (1.52) vale anche per < 0,
purche valga la (1.51) con arg z 2 e la semicirconferenza R stia nel
semipiano inferiore.
4) Con la sostituzione z iz si vede subito che vale anche
Z
ez f (z) dz = 0 ,
>0
lim
R
(1.53)
purche valga la (1.51) con 2 arg z 3

e la semicirconferenza R , centrata
2
in qualsiasi punto x0 dellasse reale, stia a sinistra della parallela dellasse immaginario passante per x0 . Ovviamente se < 0 vale un discorso analogo per
una semicirconferenza R che stia a destra della parallela dellasse immaginario
passante per x0 .
Per capire subito su quale semicirconferenza chiudere il cammino per poter applicare
il lemma di Jordan, basta ricordare che essa va scelta in modo che, lungo la sua freccia,
lesponente del fattore che moltiplica f (z) deve essere reale e tendere a per |z|
.
Si indica con lemma di Jordan il contenuto dei punti 2), 3) e 4), ma per comodità
denoteremo con questo termine tutto quanto detto in questo paragrafo.
1.6.4
Esempi
Esempio 1
Z
I=
0
x2
dx
(x2 + 1)(x2 + 4)
Notare che questa condizione è molto meno restrittiva della (1.50); qui è lesponenziale, rapidamente decrescente per Im z +, che si incarica di far tendere rapidamente a zero
lintegrando.
52
La funzione integranda è simmetrica:

f (x) =
x2
= f (x) .
(x2 + 1)(x2 + 4)
Quindi
1
I=
2
x2
dx .
(x2 + 1)(x2 + 4)
Inoltre
|z|
f (z)
1
;
z2
le ipotesi del lemma di Jordan (caso 1) sono soddisfatte in entrambi i semipiani. Possiamo quindi chiudere il cammino di integrazione nel piano complesso come indicato
in Figura 1.18 (scegliamo di chiuderlo nel semipiano positivo).
Indichiamo con CR il cammino chiuso e con R la semicirconferenza. Il lemma di
Jordan ci assicura che
Z
f (z)dz = 0
lim
e quindi
Z
I
f (z)dz = lim
f (z)dz .
CR
Pertanto
1
I = lim
2 R
I
f (z)dz .
CR
Studiamo la funzione f (z):

(z 2 + 1)(z 2 + 4) = 0 z = i , z = 2i .
f (z) ha 4 poli semplici, due nel semipiano Im z > 0 e due nel semipiano Im z < 0.
Quindi
1
I = 2i [{Resf (z)}z=i + {Resf (z)}z=2i ]
2
53
{Resf (z)}z=i = lim(z i)

zi
i
z2
=
(z + 2i)(z 2i)(z 2 + 1)
3
{Resf (z)}z=2i = lim (z 2i)

z2i

I = i
z2
i
=
2
(z + i)(z i)(z + 4)
6
i
i
6 3

=
.
6
Si noti che in questo esempio, e nei successivi esempi 2 e 3, lintegrando è positivo; se

il risultato trovato fosse un numero negativo (o peggio immaginario) si sarebbe certo
commesso un errore di segno (o dimenticato un fattore i).
Esempio 2
Z
I=
dx
,
1 + x2n
f (z) =
n intero positivo
1
1 + z 2n
|z|
1
z 2n
Vale il caso 1). Chiudiamo il cammino di integrazione in Im z > 0:

I
I = lim
CR
dz
.
1 + z 2n
I poli di f (z) sono dati da

1
1 + z 2n = 0 z 2n = 1 z = (1) 2n .
Quanti poli giacciono nel semipiano Im z > 0? Poiche 1 si può rappresentare come
1 = ei(+2k)
con k intero, i poli saranno
z = zk = ei
2k+1
2n
= eik
con
k =
2k + 1
2n
54
cioè
z0 = ei 2n , z1 = ei 2n , z1 = ei 2n , etc.
I poli zk giacciono nel semipiano Im z > 0 se 0 < k < , ovvero se
0<
2k + 1
<1
2n
che, poiche k è intero, equivale a
1
1
<k <n
0 k n 1 .
2
2
La funzione f (z) ha quindi n poli semplici in Im z > 0:

z = zk = eik , k =
2k + 1
,
2n
k = 0, 1, ...(n 1) .
Il residuo di f (z) nel polo zk vale

{Resf (z)}z=zk
1
= lim (z zk )
zzk
1 + z 2n
zk
zk
.
=
=
2n
2nzk
2n

= lim
zzk
1
2nz 2n1
Pertanto
I = 2i
n1
n1
X
zk
i X i (2k+1)
=
.
e 2n
2n
n
k=0
k=0
Poniamo
i
z0 = e 2n z02k+1 = e 2n (2k+1)
Allora
n1
n1
i X (2k+1)
i X 2 k
i 1 z02n
I=
z0
= z0
z0 = z0
n k=0
n k=0
n 1 z02
Ma z02n = 1 e pertanto
I=
2i z0
=
2
n 1 z0
n
e infine
I=
n sin
2n
55
.
1
z0 z01
2i
Esempio 3
Z
I =
cos x
dx .
1 + x2

Lintegrale esiste poiche la funzione integranda è continua sullasse reale ed è O x12 per
x ; non si pu`
o però applicare il caso 1) del lemma di Jordan perche lintegrando
1
non è affatto O z nel piano complesso; infatti
1
I=
2
eiz + eiz
dz ,
1 + z2
(1.54)
quindi esso diverge esponenzialmente per z = iy con y . Invece il primo addendo

dellintegrale (1.54) soddisfa le ipotesi del caso 2) del lemma di Jordan ( = 1 > 0) e
quindi si calcola chiudendo il cammino nel semipiano superiore:

Z +
eiz
eiz
dz = 2i Res
= ;
2
2
1 + z z=i
e
1 + z
il secondo addendo ricade invece nel caso 3) ( = 1 < 0) e quindi si calcola chiudendo
il cammino nel semipiano inferiore

Z + iz
e
eiz
dz = 2i Res
= .
2
2
1 + z z=i
e
1 + z
Pertanto
.
e
I=
Esempio 4
Z
I=
0
sin x
dx .
x
La funzione integranda è pari:

f (x) =
sin x
= f (x) .
x
Quindi
1
I=
2
sin x
dx .
x
Notare che lintegrando non è singolare nellorigine; infatti lo zero semplice del denominatore è compensato da uno zero semplice del numeratore.
56
Figura 1.19: Cammini C1 e C2 che aggirano lorigine
Il lemma di Jordan non è direttamente applicabile, perche

1 eiz eiz
sin z
=
;
z
2i
z
quindi per il primo addendo bisognerebbe applicare il caso 2) ( = 1 > 0) e chiudere
con una semicirconferenza nel semipiano superiore, mentre il secondo ricade nel caso
3) ( = 1 < 0) e bisognerebbe chiudere nel semipiano inferiore. Non si può nemmeno
spezzare lintegrale in una somma di due integrali, perche
Z + iz
e
dz
z
non esiste (lo zero del denominatore non è pi`

u compensato da uno zero del numeratore).
La difficoltà si aggira nel modo seguente: poiche f (z) è ovunque analitica al finito
(si noti che f (z) 1 per z 0), prima di spezzare lintegrale si può deformare il
cammino di integrazione, grazie al teorema di Cauchy. In particolare, i due cammini
C1 e C2 di Figura 1.19 danno lo stesso risultato per I:
Z
Z
1
sin z
1
sin z
dz =
dz .
I=
2 C1 z
2 C2 z
Dopo aver cos` deformato il cammino è possibile spezzare lintegrale in una somma di
due e procedere con il metodo dei residui. Calcoleremo I in due modi diversi:
i) Integriamo su C1 :
Z

Z
1
eiz
eiz
I=
dz
dz .
4i
C1 z
C1 z
iz
iz
Benche la funzione f (z) sia regolare ovunque, le funzioni ez e e z hanno un polo

semplice in z = 0; inoltre esse soddisfano il lemma di Jordan nei semipiani Im z > 0
57
Figura 1.20: Chiusura del cammino che aggira lorigine nel semipiano superiore (a) ed
inferiore (b)
e Im z < 0, rispettivamente. Pertanto si possono chiudere i cammini di integrazione

nelle curve 1 e 2 (vedi Figura 1.20).
I iz
I iz
e
e
1
dz
dz .
I=
4i
1 z
2 z
Si noti che la curva 2 è percorsa in senso orario.
Ora, il primo integrale è nullo per il teorema di Cauchy e il secondo si calcola con
il teorema dei residui, tenendo conto del cambiamento di segno necessario perchè la
curva 2 è percorsa in senso orario:
1
I=
4i
I
2

eiz
1
eiz
= .
dz = + 2i Res
z
4i
z z=0
2
` facile verificare che lo stesso risultato si ottiene integrando su C2 .

E
ii)
Z

Z
1
eiz
eiz
I=
dz
dz .
4i
C1 z
C1 z
Cambiamo variabile nel secondo integrale:
z z ,
C1 C2
Allora
1
I=
4i
Z
C1
eiz
dz
z
Z
C2
eiz
dz
z
58
1
=
4i
Z
Z
+
C1
C2
eiz
dz .
z
(1.55)
Ma
Z
C1
=
C2
Z
=
C2
I
+
C1
dove è una curva chiusa che circonda lorigine. Pertanto

Z
I iz
eiz
e
1
2
dz +
dz .
I=
4i
C1 z
z
(1.56)
Ora,
Z
C1
eiz
dz = 0
z
per il lemma di Jordan (applicabile in Im z > 0) e per il teorema di Cauchy, e

I iz
e
eiz
dz = 2i Res
= 2i .
z z=0
z
Sostituendo infine nella (1.56) si ottiene
I=
1.6.5
.
2
Singolarit`
a sul cammino di integrazione
A differenza di quello in (1.55) lintegrale

Z
cos x
I=
dx
x
(1.57)
ha una singolarità sul cammino dintegrazione e non è definito. Nel caso in questione, e
pi`
u generalmente per integrali logaritmicamente divergenti con singolarità 1/(x x0 ),
si può introdurre il valore principale P dellintegrale (o PV), definito da
Z

Z
Z
cos x
cos x
cos x
P
dx = lim
dx +
dx
(1.58)
0
x
x
x

che risulta finito, per via della cancellazione esatta della singolarità a destra dellorigine
con quella (negativa) a sinistra dellorigine. Nel caso specifico il valore principale è
nullo, perchè i due integrali in (1.58) si cancellano:
Z
Z
Z
cos x
cos(x)
cos x
dx =
d(x) =
dx .
(1.59)
x
x
(x)

59
Si può anche attribuire allintegrale (1.57) un significato deformandone il cammino

dintegrazione sul piano complesso in modo da evitare la singolarità aggirandola in uno
dei due modi in Fig.1.19, per es.:

Z
Z iz
cos z
1
eiz
e
I=
dz =
+
dz .
z
2 C1 z
z
C1
Chiudendo poi in cammino dintegrazione sul semipiano Imz > 0(< 0) per eiz (eiz ) si
ottiene
I iz
I iz
e
e
1
dz +
dz = i
I=
2
1 z
2 z
dove abbiamo usato (1.55). Se si sceglie invece C2 si ottiene I = +i. La parte
reale è nulla e coincide con il valor principale. E importante notare che lintegrale sul
cammino C1 è identico a quello ottenuto su una retta parallela allasse reale e posta i
sopra ad esso, perchè in tutti i punti, tranne nellintorno dellorigine, si può mandare
a zero. Si noti anche che la deformazione del cammino in C1 o C2 è equivalente allo
spostamento della singolarità sotto o sopra lasse reale di una quantità infinitesima, e
cioè
Z
Z
cos x
cos z
dz =
dx
z
x i
C1,2
con > 0.
In generale, lintegrale
Z
f (x)
dx
x x0
(1.60)
dove f (x) è non nulla in x0 , regolare in R, e tale da assicurare la convergenza dellintegrale per x , può essere trattato nello stesso modo:
Z
Z
f (x)
f (x)
dx = P
dx if (x0 )
(1.61)
x x0 i
x x0
dove il valore principale è definito in maniera analoga alla (1.58) e il valore f (x0 )
emerge dal calcolo del residuo in z = x0 (nellesempio (1.57) cos z|z=0 = 1). In termini
di distribuzioni (che studieremo pi`
u avanti) si scrive anche
1
1
=P
i (x x0 )
x x0 i
x x0
60
(1.62)
Figura 1.21: Proiezione stereografica della sfera sul piano complesso.
1.7
Studio del punto allinfinito
Una sfera di Riemann (vedi figura) è topologicamente equivalente al piano complesso

esteso al punto allinfinito: C = C {}; a ogni punto sulla sfera corrisponde infatti
un unico punto sul piano complesso. Al polo Nord della sfera corrisponde il punto
allinfinito, , che ha argomento indefinito e modulo +. Lestensione di C a C è un
esempio di compattificazione.
Il comportamento della funzione f (z) per z si studia effettuando il cambiamento di variabile
1
t=
(1.63)
za
per un opportuno a C, che generalmente si prende uguale a zero, e valutando il comportamento della funzione (t) = f (a + 1/t) per t 0. La sostituzione (1.63) manda
un intorno circolare (di raggio ) dellorigine nel piano complesso di t in un intorno
dellinfinito I (), cioè nellesterno di un cerchio di raggio = 1/ centrato in a.
Per esempio, se t = 0 è un polo di ordine n di (t), z = è un polo di ordine n di
f (z); se t = 0 è uno zero di ordine n di (t), z = è uno zero di ordine n di f (z),
eccetera.
Esempi
La funzione f (z) =
1
z
è regolare allinfinito e ivi ha uno zero semplice.
La funzione f (z) = z 2 ha un polo doppio allinfinito.

La funzione f (z) = ez ha una singolarità essenziale allinfinito, cos` come le
funzioni sin z e analoghe.
La funzione f (z) = sin1 z ha poli semplici nei punti zk = k con k intero qualsiasi; quindi in ogni intorno del punto allinfinito cade almeno un polo (in realtà
ne cadono infiniti): linfinito non è una singolarità isolata, ma un punto di
accumulazione di poli.
61
Naturalmente, se f (z) è regolare allinfinito può essere sviluppata in serie di Taylor

in un intorno dellinfinito I (), cioè allesterno di un cerchio, centrato in un
punto a scelto secondo convenienza (spesso a = 0), e di raggio tale che allesterno
del cerchio la f (z) non abbia singolarità. Tale serie si ottiene sviluppando in serie di
Taylor la funzione (t) = f (a + 1/t) nellintorno del punto t = 0, tornando poi alla
variabile originaria z con la sostituzione t = 1/(z a); lo sviluppo di Taylor nellintorno
del punto allinfinito conterra quindi solo potenze negative di z a, oltre alla potenza
nulla.
Per esempio lo sviluppo (1.32) della funzione e1/z , che abbiamo già visto essere lo
sviluppo di Laurent nellintorno della singolarità essenziale z = 0, può anche essere
letto come lo sviluppo di Taylor nellintorno del punto regolare z = .
Discorso analogo per lo sviluppo di Laurent; solo che stavolta la parte principale
dello sviluppo (cioè quella singolare) conterrà solo potenze positive di (z a), in numero
finito o infinito a seconda se il punto allinfinito è un polo o una singolarità essenziale.
Per ogni funzione intera lo sviluppo di Taylor
f (z) =
an z n
(1.64)
n=0
nellintorno dellorigine può anche essere letto come sviluppo di Laurent intorno allinfinito; quindi:
se ci sono infiniti an 6= 0 linfinito è una singolarità essenziale di f (z);
se an 6= 0 e al = 0, l > n, f (z) è un polinomio di grado n e linfinito è un polo
di ordine n (per n 6= 0) o è regolare (per n = 0).
Ne segue anche che:
una funzione regolare in tutto C e anche allinfinito è necessariamente una costante (in accordo con il teorema di Liouville - vedi Appendice B).
1.7.1
Esempi
Esempio 1: la funzione
f (z) = ez
è regolare in z = 0. Infatti, lo sviluppo in serie che definisce la funzione esponenziale
e =
X
zk
k=0
k!
ha solo potenze positive (sviluppo di Taylor intorno a z = 0). La stessa serie può essere
letta come sviluppo di Laurent attorno alla singolarità essenziale z = .
62
Esempio 2: la funzione
f (z) = e1/z
ha una singolarità essenziale in z = 0. Infatti, lo sviluppo in serie che definisce la

funzione esponenziale
1/z 2
X
(1/z 2 )k
k!
k=0
X
(1)k
k=0
k!
z 2k
ha un numero infinito di potenze negative. La stessa serie può essere letta come sviluppo
di Taylor attorno al punto regolare z = ; essa non contiene infatti potenze positive
di z; pertanto f (z) è analitica in z = .
Esempio 3: consideriamo la funzione
f (z) = ze
1/z
=z
X
(1/z)k
k!
k=0
X
(1)k
k=0
k!
k+1
1
X
(1)k
0
=
zk .
0
(k + 1)!
k0 =
La serie ha un numero infinito di potenze negative e quindi il punto z = 0 è una

singolarità essenziale. Studiamo z = : la serie contiene una potenza positiva (la
prima) di z, quindi la funzione ha un polo di ordine 1 in z = .
Esempio 4: consideriamo la funzione
f (z) = ez/(1z) .
Poniamo z 0 = z 1:
f (z) = e(1+z )/z = e1/z e1 =
1X
(z 1)k
1 X (1/z 0 )k
=
(1)k
e k=0
k!
e k=0
k!
La serie ha un numero infinito di potenze negative e quindi il punto z = 1 (z 0 = 0) è

una singolarità essenziale.
Per z (z 0 = ) la funzione ammette limite (uguale a e1 ) e quindi il punto
allinfinito è regolare; la serie che abbiamo scritto, che è di Laurent attorno al punto
z = 1, può anche essere letta come serie di Taylor nellintorno dellinfinito.
Esempio 5: sia
f (z) = ez1/z = ez e1/z .
I punti z = 0 e z = sono singolarità essenziali.
63
Figura 1.22: Curva che contiene tutte le singolarità al finito della funzione f (z).
1.7.2
Calcolo del residuo nel punto allinfinito
Per valutare il residuo di una funzione f (z) in z = supponiamo che esista una curva
di Jordan che contenga al suo interno tutte le singolarità al finito di f (z) (per esempio
una circonferenza c di raggio sufficientemente grande - Figura 1.22- )
Allora nel punto allinfinito la funzione f (z) o è regolare, o ha una singolarità isolata.
In entrambi i casi definiamo il residuo allinfinito come:
I
1
f (z)dz ,
(1.65)
{Resf (z)}z= =
2i
dove lintegrale è calcolato percorrendo come al solito la curva in senso antiorario; il
segno ricorda che per avere z = al suo interno la curva dovrebbe essere percorsa
in senso orario8 .
Sostituendo nella definizione (1.65) lo sviluppo in serie di Laurent (o di Taylor)
della f (z) attorno al punto allinfinito, e ricordando lintegrale (1.18), si ricava che il
residuo allinfinito è dato dal coefficiente, cambiato di segno, della potenza 1/(z a).
Per esempio, il residuo allinfinito della funzione regolare allinfinito f (z) = 1/z vale
1, mentre quello della funzione f (z) = z (che ha un polo semplice allinfinito) è nullo.
Una conseguenza immediata di quanto abbiamo detto è che una funzione pari ha
sempre residuo nullo allinfinito (sempre che abbia senso definirlo), poiche il suo
sviluppo, di Taylor o di Laurent, in potenze di z non potrà contenere il termine 1/z;
lo stesso succede per una funzione che allinfinito sia O( z12 ).
8
Ricordiamo che un punto z0 si definisce interno ad una curva se, immaginando di percorrere
nel senso di percorrenza indicato, z0 viene lasciato a sinistra.
64
Un altro modo per calcolare il residuo allinfinito si basa sul calcolo diretto dellintegrale (1.65) mediante il cambiamento di variabile:
1
z= ,
t
(1.66)
che manda z in t 0.
Si vede subito che una circonferenza c di raggio R centrata nellorigine del piano z,
di equazione |z| = R, viene trasformata in unanaloga circonferenza c0 di raggio 1/R
nel piano t di equazione |t| = 1/R. Se c è percorsa in senso antiorario, c0 sarà percorsa
in senso orario; infatti quando la fase di z = Rei cresce, quella di t = 1/R ei
diminuisce. La circonferenza cR viene quindi mappata dalla trasformazione (1.66) in
una circonferenza c01/R percorsa in senso opposto. Di conseguenza:
I
1
f (z)dz
{Resf (z)}z= =
2i c
I
1
1 dz
= +
f ( ) dt .
2i c0 t dt
Ora, dz/dt = 1/t2 , e quindi
{Resf (z)}z=
I
1
1 1
=
f
dt
2i c0
t t2

1 1
= Res f
.
t t2 t=0
(1.67)
Si noti che il residuo di f (z) in z = non è uguale al residuo di f (1/t) in t = 0:

{Resf (z)}z= 6= {Resf (1/t)}t=0 .
N.B. Esistono funzioni che, pur essendo regolari in z = , hanno residuo non nullo
allinfinito. Per esempio la funzione
f (z) =
1
z
è regolare in z = (perche la funzione f (1/t) = t è uguale a 0 in t = 0) ma il suo

residuo, calcolato tramite la (1.67), vale

1
Res
z

z=

I
1
1
1
= Res t 2
=
dt = 1 ,
t
2i c0 t
t=0
come abbiamo già visto.

Linteresse principale nel definire il residuo allinfinito sta nel seguente:
65
Teorema 14: se una funzione analitica f (z) possiede solo singolarità isolate in
tutto il piano complesso, punto allinfinito compreso, la somma di tutti i suoi residui,
compreso leventuale residuo allinfinito, è zero.
Dimostrazione. Sia una curva di Jordan che non passa per alcuna singolarità di
f (z). Allora, per il teorema dei residui (1.39), si ha
I
X
f (z)dz = +2i
{Resf (z)}
interni
I
f (z)dz = 2i
{Resf (z)} ,
esterni
da cui si ottiene, sottraendo membro a membro,

X
{Resf (z)} = 0 .
(1.68)
tot
[q.e.d.]
La verifica pi`
u immediata di questo teorema è data dalla solita funzione f (z) = 1/z
che ha residuo +1 nellorigine e 1 allinfinito; conviene richiamare questo esempio elementare ogni volta che non ci si ricordi con quale segno si debba prendere il coefficiente
della potenza 1/(z a) per calcolare il residuo allinfinito.
Esempio 1
A volte può essere conveniente usare leq. (1.68) per semplificare il calcolo di integrali
in campo complesso. Per esempio lintegrale
I
z3
dz con C = {z, |z| = 1}
(1.69)
4
C 2z + 1
richiederebbe di valutare i 4 residui interni alla curva C, nei punti zi soluzioni di
z 4 = 1/2. Utilizzando invece il teorema (1.68) si ha semplicemente

I
1 t13
z3
z3
= +2i lim 2 2
dz = 2i Res 4
4
t0 t 4 + 1
2z + 1 z=
C 2z + 1
t
= i .
(1.70)
Ancor pi`
u semplicemente si trova che il residuo allinfinito dellintegrando è 1/2
guardando allo sviluppo:
z3
1
1
=
+
O(
)
2z 4 + 1
2z
z2
66
Esempio 2
Il teorema 14 permette a volte di calcolare pi`
u facilmente il residuo di una funzione in
una singolarità essenziale. Per esempio il calcolo del residuo della funzione
f (z) =
sin(/z)
z2
(1.71)
nella singolarità essenziale z = 0 è

{Resf (z)}z=0 = {Resf (z)}z=2 == lim sin
z2
= 1 ,
z
dove si è tenuto conto che f (z) = O( z12 ) per z e quindi {Resf (z)}z= = 0. Invece
il calcolo diretto è pi`
u complicato:
#)
2k+1
1 X z l X
Res
(1)k z
2 l=0 2 k=0
(2k + 1)!
(
{Resf (z)}z=0 =
"
z=0
2k+1
(1)
l+1 l2k1,1
(2k
+
1)!
2
l,k=0
X
(1)k 2k+1
= sin = 1 .
(2k + 1)! 2
2
k=0
Corollario del Teorema 14: ogni funzione intera f (z) ha residuo nullo allinfinito.
Dimostrazione. Linfinito può essere punto regolare di f (z) (allora f (z) è costante Teorema di Liouville, Appendice B) o singolarità isolata; in entrambi i casi ha senso
definire il residuo allinfinito. Poichè la somma dei residui fa zero e non ci sono singolarità al finito, si deduce che necessariamente il residuo allinfinito è nullo. In alternativa,
basta vedere che lo sviluppo (1.64) di f (z) in serie di Laurent nellintorno dellinfinito
non contiene la potenza z 1 .
67
1.8
Le funzioni ln z e z nel piano complesso
La funzione logaritmo w = ln z è definita nel campo complesso z 6= 0 dalla equazione

ew = z ,
(1.72)
ovvero
eRe w eiIm
= |z|eiarg
(1.73)
da cui segue, prendendo il modulo di ambo i membri,

eRe
= |z| Re w = ln |z| ,
(1.74)
dove il logaritmo del numero positivo |z| è quello definito in campo reale. Sostituendo
nella (1.73) si ottiene
eiIm w = eiarg z ,
(1.75)
la cui soluzione generale è
Im w = arg z + 2n , n Z .
(1.76)
Si rimuove lambiguità, insita nella stessa definizione di arg z, fissando il valore di n

(per esempio n = 0) e scegliendo un intervallo di ampiezza 2 in cui far variare arg z,
come discusso nel par.1.1.2. Quindi
ln z = ln |z| + i arg z .
(1.77)
La funzione logaritmo è perciò completamente definita solo se si specifica in che intervallo varia arg z ed è discontinua su una semiretta (taglio) uscente dallorigine del
piano complesso, perche tale è la funzione arg z, come si è visto nel paragrafo 1.1.2.
Lorigine è perciò una singolarit`
a non isolata della funzione logaritmo detta punto di diramazione (branching point); lo stesso dicasi per il punto allinfinito. La
funzione z , con C, si definisce come
z = e ln z
(1.78)
` facile
e, salvo che per N, soffre degli stessi problemi della funzione logaritmo9 . E
dimostrare che
d ln z
1
= , z 6= 0
(1.79)
dz
z
e quindi anche
dz
= z 1 , C .
(1.80)
dz
Per N non ci sono ambiguit`

a poiche eN (ln z+2in) = eN ln z .
68
Capitolo 2
Equazioni differenziali in C
2.1
Equazioni differenziali ordinarie II ordine
La forma pi`
u generale di equazione differenziale ordinaria del II ordine omogenea è
A(z)u00 (z) + B(z)u0 (z) + C(z)u(z) = 0 .
(2.1)
Dividendo per A(z) (supposto diverso da zero, altrimenti lequazione sarebbe del I
ordine) si ottiene la cosiddetta forma standard
u00 (z) + P (z)u0 (z) + Q(z)u(z) = 0 .
(2.2)
Note due soluzioni dellequazione omogenea è sempre possibile risolvere, almeno in

linea di principio, lequazione inomogenea
u00 (z) + P (z)u0 (z) + Q(z)u(z) = f (z) .
Condizione necessaria e sufficiente affinchè due soluzioni u1 (z) e u2 (z) della (2.3) siano
linearmente indipendenti è che il wronskiano differisca da zero, ovvero

u1 u2
6= 0
W (z) = det 0
u1 u02
Si noti che il wronskiano è sempre nullo o sempre diverso da zero. Di conseguenza
due soluzioni che si annullano in z0 sono la stessa soluzione. Soluzioni linearmente
indipendenti non hanno zeri in comune.
Punti regolari e singolari di una equazione differenziale
Consideriamo la forma (2.3) di unequazione differenziale del II ordine omogenea. Le
proprietà delle soluzioni dipendono dal comportamento delle funzioni P (z) e Q(z) nel
campo complesso; se esse sono regolari nel punto z = z0 , il punto z0 si dice punto
69
regolare, o ordinario, dellequazione differenziale e qualunque soluzione è regolare in z0 .

Altrimenti il punto z0 si dice punto singolare dellequazione differenziale poichè generalmente le soluzioni saranno ivi singolari. I punti singolari sono a loro volta classificati
in due categorie: singolarit`
a fuchsiane, o regolari, e singolarit`
a essenziali, o
irregolari. Il punto singolare z0 si definisce punto singolare fuchsiano (dal nome del
matematico Fuchs) se in z z0 la funzione P (z) ha al pi`
u un polo semplice e Q(z) al
pi`
u un polo doppio; quindi le funzioni (z z0 )P (z) e (z z0 )2 Q(z) rimangono finite
per z z0 :
lim (z z0 )P (z) = p0
zz0
lim (z z0 )2 Q(z) = q0 ,
zz0
con p0 e q0 finiti; è possibile che uno o anche entrambi siano nulli. Se invece per
esempio P (z) diverge pi`
u velocemente di 1/(z z0 ), in modo tale che (z z0 )P (z)
tenda a infinito per z z0 , oppure se Q(z) diverge pi`
u velocemente di 1/(z z0 )2 ,
2
in modo tale che (z z0 ) Q(z) tenda a infinito per z z0 , il punto z0 è un punto
singolare irregolare, o essenziale. Queste definizioni valgono per tutti i valori finiti di
z0 . Lo studio del punto z verrà trattato separatamente in un prossimo paragrafo.
Esempi
Elenchiamo alcuni esempi di equazioni differenziali ordinarie del II ordine e studiamone
le singolarità al finito.
1) Equazione delloscillatore armonico semplice:
u00 + 2 u = 0
(2.3)
P (z) = 0 , Q(z) = 2
Lequazione è ovunque regolare al finito.
2) Equazione di Legendre:
(1 z 2 )u00 2zu0 + u = 0
(2.4)
2z
,
Q(z)
=
1 z2
1 z2
Lequazione ha due punti singolari fuchsiani in z = 1. Infatti sia P (z) che Q(z)
hanno un polo semplice in z = 1:
P (z) =
lim (z (1))P (z) = 1 = p0
z1
lim (z (1))2 Q(z) = 0 = q0 .
z1
70
3) Equazione di Bessel:
z 2 u00 + zu0 + (z 2 2 )u = 0
(2.5)
1
2
, Q(z) = 1 2
z
z
Lequazione ha una singolarità di tipo fuchsiano in z = 0 con p0 = 1 e q0 = 2 .
P (z) =
4) Equazione di Laguerre
zu00 + (1 z)u0 + au = 0
(2.6)
1
a
1 , Q(z) =
z
z
Lequazione ha una singolarità di tipo fuchsiano in z = 0.
P (z) =
5) Equazione di Hermite:
u00 2zu0 + 2u = 0
(2.7)
P (z) = 2z , Q(z) = 2
Lequazione è regolare al finito.
6) Equazione di Chebyshev:
(1 z 2 )u00 zu0 + n2 u = 0
(2.8)
n2
z
,
Q(z)
=
1 z2
1 z2
Lequazione ha due punti singolari fuchsiani in z = 1.
P (z) =
7) Equazione ipergeometrica:
z(z 1)u00 + [(1 + a + b)z c]u0 + abu = 0
P (z) =
(2.9)
(1 + a + b)z c
ab
, Q(z) =
z(z 1)
z(z 1)
Lequazione ha due punti singolari fuchsiani in z = 0 e z = 1.

8) Equazione ipergeometrica confluente:
zu00 + (c z)u0 au = 0
cz
a
, Q(z) =
z
z
Lequazione ha un punto singolare fuchsiano in z = 0.
P (z) =
71
(2.10)
2.1.1
Soluzione nellintorno di un punto regolare
Se z0 è un punto ordinario di unequazione differenziale nella forma standard, le solu` allora possibile cercare una soluzione che abbia
zioni sono certamente regolari in z0 . E
la forma di uno sviluppo in serie di Taylor intorno al punto z0 ,
u(z) =
ck wk ,
(2.11)
k=0
con w z z0 , e determinare i coefficienti ck mediante la sostituzione della serie

nellequazione differenziale. Poiche P (z) e Q(z) sono analitiche in z0 valgono gli sviluppi
in serie di Taylor
P (z) =
l=0
Q(z) =
pl w l
(2.12)
ql w l .
(2.13)
l=0
Sostituendo le (2.12) e (2.13) e le derivate di u(z)

0
u (z) =
u00 (z) =
X
l=1
lcl w
l1
(n + 1)cn+1 wn
l(l 1)cl wl2 =
(n + 1)(n + 2)cn+2 wn
(2.15)
n=0
l=2
nellequazione differenziale (2.3) si ottiene

#
"
n
n
X
X
X
cl qnl wn = 0 .
(l + 1)cl+1 pnl +
(n + 1)(n + 2)cn+2 +
n=0
(2.14)
n=0
(2.16)
l=0
l=0
Una serie di potenze è nulla se e solo se tutti i suoi coefficienti sono nulli, e pertanto
lespressione in parentesi quadra deve annullarsi, n. Si ottengono cos` delle relazioni
di ricorrenza che permettono di determinare i coefficienti ck una volta noti c0 e c1 .
Infatti per n = 0 si ottiene:
2c2 + c1 p0 + c0 q0 = 0 ;
(2.17)
6c3 + c1 p1 + 2c2 p0 + c0 q1 + c1 q0 = 0
(2.18)
per n = 1:
e cos` via. Le costanti arbitrarie c0 e c1 , fissate dalle condizioni iniziali
c0 = u(z0 )
c1 = u0 (z0 ) ,
72
(2.19)
(2.20)
determinano univocamente la soluzione u(z). Se per esempio chiamiamo u1 la soluzione

corrispondente a c0 = 1 e c1 = 0 e u2 quella corrispondente a c0 = 0 e c1 = 1, la
soluzione generale dellequazione differenziale sarà
u(z) = c0 u1 (z) + c1 u2 (z) ;
(2.21)
infatti u1 e u2 sono linearmente indipendenti, essendo il loro Wronskiano diverso da

zero:

1 0
u1 u2

=1.
= det
W (z0 ) = det 0
0 1
u1 u02 z=z0
In generale si può dimostrare che questo metodo fornisce sempre la soluzione generale
nellintorno di un punto regolare z0 e che, per valori generici di c0 e c1 il raggio di
convergenza della serie è uguale alla distanza fra z0 e la singolarità pi`
u vicina dellequazione differenziale (talvolta, ma solo per particolari valori di c0 e c1 , può anche
essere maggiore).
Esempi
1. L equazione delloscillatore armonico semplice
u00 (z) + 2 u(z) = 0 ,
(2.22)
è, come si è detto, regolare per ogni z finito, in particolare per z = 0. Possiamo quindi
cercare una soluzione del tipo
X
ck z k .
u(z) =
k=0
Sostituendo nella (2.16) z0 = 0, pi = 0 e qi = 2 1,0 si ottiene:

ck+2 (k + 1)(k + 2) + 2 ck = 0
k 0 .
(2.23)
da cui segue la relazione di ricorrenza per i coefficienti ck :

ck+2 =
2
ck .
(k + 1)(k + 2)
(2.24)
Dati i coefficienti c0 e c1 , che saranno determinati dalle condizioni al contorno, la (2.24)

permette di costruire tutti i coefficienti delle potenze pari
2
c0
2
2
( 2 )2
=
c2 =
c0
(3)(4)
4!
2
( 2 )3
=
c4 =
c0
(5)(6)
6!
2 n
n ( )
= (1)
c0
(2n)!
c2 =
c4
c6
c2n
73
e delle potenze dispari

2
c1
(2)(3)
2
( 2 )2
=
c3 =
c1
(4)(5)
5!
( 2 )3
2
c5 =
c1
=
(6)(7)
7!
( 2 )n
= (1)n
c1 .
(2n + 1)!
c3 =
c5
c7
c2n+1
Pertanto la soluzione cercata è
u(z) = c0
X
(z)2n c1 X
(z)2n+1
(1)n
+
(1)n
(2n)!
n=0
(2n + 1)!
n=0
= c0 cos(z) + c01 sin(z) ,
(2.25)
che è proprio, come noto, la soluzione dellequazione (2.22).

2. Lequazione di Legendre
(1 z 2 ) u00 (z) 2z u0 (z) + u(z) = 0 ,
(2.26)
compare nella soluzione dellequazione di Laplace in coordinate sferiche e in molte altre

applicazioni. Poichè il punto z = 0 è un punto regolare dellequazione, cercheremo una
soluzione data dalla serie:
u(z) =
ck z k .
(2.27)
k=0
Per determinare i coefficienti calcoliamo le derivate

0
u (z) =
kck z
00
k1
u (z) =
k=0
k(k 1)ck z k2
k=0
e sostituiamole nella (2.26):
k(k 1)ck z k2
k=0
ck [k(k 1) + 2k ] z k = 0 .
(2.28)
k=0
La prima sommatoria si può riscrivere come segue (si noti che i termini k = 0, 1 sono
nulli):
X
k=2
k(k 1)ck z k2 =
(k 0 + 2)(k 0 + 1)ck0 +2 z k ;
k0 =0
74
lequazione (2.28) diventa quindi:
{(k + 2)(k + 1)ck+2 ck [k(k + 1) ]} z k = 0 .
k=0
Uguagliando a zero ogni coefficiente della serie di potenze si ottiene la relazione di

ricorrenza:
ck+2 = ck
k(k + 1)
,
(k + 2)(k + 1)
da cui si ricava per generico
c0
2
(6 )()
6
c2 =
c0
=
12
4!
20
(20 )(6 )()
=
c2 =
c0
30
6!
c2 =
c4
c6
etc...
e
2
c1
6
12
(12 )(2 )
=
c3 =
c1
20
5!
30
(30 )(12 )(2 )
=
c5 =
c1
42
6!
c3 =
c5
c7
etc...
Quindi se scegliamo c0 = 1 e c1 = 0 tutti i coefficienti delle potenze dispari nella (2.27)

sono nulli e la soluzione è pari (u(z) = u(z)), mentre per c0 = 0 e c1 = 1 la soluzione
è dispari (u(z) = u(z)). Notiamo anche che, poiche il il raggio di convergenza della
serie (2.27) è
ck
lim
=1,
(2.29)
k ck+2
ci aspettiamo almeno una singolarità sulla circonferenza |z| = 1: lequazione ha infatti
due punti singolari fuchsiani in z = 1. Le soluzioni trovate valgono C. Ora, se
e solo se = n(n + 1), con n intero positivo o nullo, il coefficiente cn+2 è nullo e cos`
cn+4 = cn+6 = = 0; la serie (2.27) si riduce cos` a un polinomio di grado n:
Pn (z) =
n
X
k=0
75
ck z k
(2.30)
Figura 2.1: I primi 6 polinomi di Legendre. Non hanno zeri fuori da [1, 1].
se n è pari per c0 = 1 e c1 = 0, oppure se n è dispari per c0 = 1 e c1 = 0. Per c0 e

c1 generici la soluzione generale è la combinazione lineare di un polinomio e una serie.
I polinomi (2.30) (opportunamente normalizzati) si chiamano polinomi di Legendre e
sono regolari in tutto il piano complesso. Si noti che la richiesta che u(z) sia finita in
tutto lintervallo [1, 1] seleziona gli = n(n + 1), in accordo con quanto si vedrà in
generale nel capitolo 5 a proposito degli autovalori di operatori differenziali.
3. Lequazione di Hermite
u00 2z u0 + 2 u = 0
(2.31)
compare nello studio delloscillatore armonico quantistico ed è regolare in z = 0.

Cerchiamo una soluzione
X
u(z) =
ck z k .
(2.32)
k=0
Questa, sostituita con le sue derivate nella (2.31), fornisce

ck k(k 1)z k2 2kz k + +2z k = 0 .
k=0
76
(2.33)
Il primo termine della serie può essere riscritto, cambiando lindice di somma da k in
k 2, come visto in precedenza. Sostituendo nella (2.33) si ottiene
z k [ck+2 (k + 2)(k + 1) + ck (2 2k)] = 0 ,
(2.34)
k=0
da cui si ricava la relazione di ricorrenza

ck+2 = ck
2(k )
.
(k + 2)(k + 1)
(2.35)
Leq. (2.35) mostra che la serie (2.32) ha raggio di convergenza infinito. Scegliendo
c0 = 1, c1 = 0 oppure c0 = 1, c1 = 0 si ottengono soluzioni pari o dispari, come si è
visto per i polinomi di Legendre. Se = n è intero positivo o nullo, una delle due serie
(quella pari se è pari, quella dispari se è dispari) si riduce a un polinomio di grado
n, il polinomio di Hermite Hn .
2.1.2
Soluzione nellintorno di un punto singolare fuchsiano
Esistono talvolta soluzioni particolari di unequazione differenziale che sono regolari in

un punto z0 singolare, come abbiamo già visto nel caso dei polinomi di Legendre che
sono regolari in z = 1. Tuttavia, la soluzione generale non può essere regolare in z0
singolare, e viceversa. Supponiamo infatti che ci siano due soluzioni u1 e u2 regolari
in z e che esse siano linearmente indipendenti. Siccome u1,2 sono entrambe soluzioni,
il sistema
u001 (z) + P (z)u01 (z) + Q(z)u1 (z) = 0
u002 (z) + P (z)u02 (z) + Q(z)u2 (z) = 0
permette di ottenere P e Q algebricamente, come rapporto di determinanti nelle funzioni regolari ui , u0i e u00i (i = 1, 2). A denominatore cè il Wronskiano di u1 e u2 ,
ovunque diverso da zero perche u1 e u2 sono indipendenti. Ne consegue che P (z) e
Q(z) sono regolari dove la soluzione generale è regolare.
Sotto quali condizioni è possibile espandere in serie nellintorno di un punto singolare? e che tipo di serie risulta? Vale il Teorema di Fuchs: se z0 è un punto singolare
fuchsiano, esiste sempre almeno una soluzione del tipo:
u1 (z) = (z z0 )
ck (z z0 )k ,
con c0 6= 0 .
(2.36)
k=0
La seconda soluzione o è ancora della forma (2.36) oppure contiene anche un termine
aggiuntivo du1 ln(z z0 ), come vedremo nelleq. (2.42). Le serie che compaiono nella
soluzione hanno raggio di convergenza almeno uguale alla distanza fra z0 e la pi`
u vicina
singolarità dellequazione differenziale.
77
Come si determina lesponente ? Supponiamo che z0 sia un punto singolare

fuchsiano. In questo caso le funzioni P (z) e Q(z) possono essere scritte come
P
pl (z z0 )l
p(z)
P (z) =
= l=0
(2.37)
z z0
z z0
P
l
q(z)
l=0 ql (z z0 )
Q(z) =
=
,
(2.38)
(z z0 )2
(z z0 )2
dove le funzioni p(z) e q(z) sono regolari in z = z0 , e sono state quindi sviluppate in
serie di Taylor intorno a z0 . Sostituendo ora le (2.36), (2.37) e (2.38) nellequazione
(2.3) e moltiplicando per (z z0 )2 si ottiene:
ck ( + k)( + k 1)(z z0 )
k=0
+k
ck ( + k)
k=0
ck
k=0
pl (z z0 )+k+l
l=0
ql (z z0 )+k+l = 0 .
l=0
Uguagliamo ora a zero il coefficiente della potenza (z z0 ) ; poniamo cioè k = l = 0

nellequazione precedente. Otteniamo cos`:
c0 [( 1) + p0 + q0 ] = 0 ,
ovvero, per c0 6= 0,
2 + (p0 1) + q0 = 0 .
(2.39)
Lequazione (2.39), detta equazione indiciale o caratteristica dellequazione differenziale (2.3), è unequazione di secondo grado in e ha quindi due soluzioni, 1 e 2 . Per
risolvere lequazione differenziale occorre quindi risolvere lequazione indiciale (2.39),
dove
p0 = lim (z z0 )P (z)
zz0
e
q0 = lim (z z0 )2 Q(z) ,
zz0
e ricavare gli indici 1 , 2 . Scelti gli indici in modo che Re1 Re2 il teorema di Fuchs
ci assicura che esiste sempre la soluzione particolare
1
u1 (z) = (z z0 )
ck (z z0 )k ,
c0 6= 0
(2.40)
k=0
i cui coefficienti si possono determinare in modo univoco in funzione di c0 sostituendo la

serie nellequazione differenziale e ricavando delle relazioni di ricorrenza. Per risolvere
completamente lequazione differenziale occorre pero ricavare una seconda soluzione,
linearmente indipendente dalla prima. Distinguiamo due casi
78
1) Se le due radici differiscono per un numero non intero, la seconda soluzione è

simile alla prima:
2
u2 (z) = (z z0 )
dk (z z0 )k ,
d0 6= 0
(2.41)
k=0
2) Se le due radici 1 e 2 differiscono per un numero intero n 0 la seconda

soluzione è
2
u2 (z) = (z z0 )
dk (z z0 )k + du1 (z) ln(z z0 ) ,
(2.42)
k=0
dove d e i dk (per k 6= 1 2 )1 si determinano per sostituzione in funzione di d0 ;

può anche succedere che si ottenga d = 0 ma solo se 1 6= 2 .
2.1.3
Esempio: lequazione di Bessel
Lequazione di Bessel
z 2 u00 (z) + zu0 (z) + (z 2 2 )u(z) = 0
(2.43)
si incontra nella soluzione dellequazione di Laplace in coordinate cilindriche e in molte

altre applicazioni. Ha un punto singolare fuchsiano in z = 0. Lequazione indiciale è
2 2 = 0 ,
con soluzioni 1,2 = . Se 2
/ Z allora 1 2
/ N ed entrambe le soluzioni sono
della forma
u (z) = z
ck z k .
(2.44)
k=0
Sostituendo la (2.44) nella (2.43) si ottiene:
ck (k + )(k + 1)z
k+
k=0
ck (k + )z
k+
k=0
ck z
k=0
k++2
ck z k+ = 0
k=0
da cui
c1 (1 + 2)z 1+ +
ck k(k + 2)z k+ +
k=2
ck z k++2 = 0 .
(2.45)
k=0
dn pu`
o essere scelto arbitrariamente poiche la differenza fra due soluzioni del tipo u2 con diversi
valori di dn è proporzionale alla prima soluzione u1 di (2.40).
79
Se nella prima sommatoria effettuiamo il cambiamento di indice k k 2 otteniamo:
X
1+
c1 (1 + 2)z
+
[ck+2 (k + 2)(k + 2 + 2) + ck ]z k++2 = 0 ,
(2.46)
k=0
da cui (essendo 6= 1/2)

c1 = 0
(2.47)
ck
.
(2.48)
ck+2 =
(k + 2)(k + 2 + 2)
La serie ha quindi solo potenze pari. Notare che se 2 fosse un intero n, diciamo
positivo, la (2.48) non avrebbe senso per 2 = = n/2 e per k = n 2, impedendo
cos` di trovare la seconda soluzione nella forma (2.44) se n è pari. Dalla relazione di
ricorrenza (2.48) ricaviamo infine che tutti i coefficienti di indice dispari sono nulli,
mentre i coefficienti pari sono:
c0
c2 =
2(2 + 2)
c0
c2
=
c4 =
4(4 + 2)
2 4(2 + 2)(4 + 2)
c0
c4
=
c6 =
6(6 + 2)
2 4 6(2 + 2)(4 + 2)(6 + 2)
c0
.
c2n = (1)n
2 4 6...2n(2 + 2)(4 + 2)(6 + 2)...(2n + 2)
Le funzioni di Bessel J (z) sono definite come segue:
J (z) = N u (z) ,
(2.49)
2
dove N è unopportuna costante di normalizzazione. Nel caso particolare = 1/4 le

due soluzioni dellequazione di Bessel sono

z2 z4
J1/2 (z) = N1/2 z 1
+
+
3!
5!

N1/2
sin z
z3 z5
+
+ = N1/2 .
=
z
(2.50)
3!
5!
z
z
e
cos z
J1/2 (z) = N1/2 .
(2.51)
z
Notare che 1 2 = 1 ma non cè il termine con il logaritmo (d = 0). Si può dimostrare
che per tutti gli semi-interi le J sono esprimibili tramite funzioni trigonometriche.
Per esempio:

N3/2 sin z
J3/2 (z) =
cos z
(2.52)
z
z

N3/2 cos z
J3/2 (z) =
sin z .
(2.53)
z
z
80
2.1.4
Studio del comportamento allinfinito
Il comportamento delle soluzioni dellequazione differenziale

u00 (z) + P (z)u0 (z) + Q(z)u(z) = 0
(2.54)
nellintorno del punto z si studia effettuando il cambiamento di variabile t = z1 e

studiando il comportamento di u( 1t ) per t 0. Che forma assume lequazione (2.54)
in termini di t?
du(1/t)
du(1/t) dt
= t2
dt dz
dt
0
du
(1/t)
du(1/t)
d2 u(1/t)
dt
u00 (z) =
= 2t3
+ t4
.
dt dz
dt
dt2
u0 (z) =
Sostituendo nella (2.54) si ottiene

t4
d2 u
du
+ (2t3 t2 P ) + Qu = 0
2
dt
dt
ovvero, in forma standard,

d2 u
+
dt2
2 P
t t2
du Q
+ 4u = 0 .
dt
t
(2.55)
Il punto t = 0 (z ) è un punto ordinario se le funzioni

1
2 P (1/t)
P
=
t
t
t2

Q(1/t)
1
Q
=
t
t4
(2.56)
(2.57)
sono regolari in t = 0, ovvero

2
P (z) = + O
z
1
z2

,
Q(z) = O
1
z4

per z .
(2.58)
Le condizioni necessarie e sufficienti affinchè il punto z = sia fuchsiano sono che le

funzioni

1
(2.59)
p(t) = tP
t

1
2
q(t) = t Q
(2.60)
t
siano regolari in t = 0, ovvero

1
P (z) = O
,
z

Q(z) = O
81
1
z2

per z .
(2.61)
Se linfinito è un punto ordinario, si può cercare una soluzione come sviluppo in serie
di Taylor intorno a z = :
u(z) = u(1/t) =
ck t =
k=0
ck z k
k=0
e determinare i coefficienti ck tramite le relazioni di ricorrenza che si ricavano dalla

sotituzione della serie nellequazione differenziale. Se invece linfinito è un punto singolare fuchsiano, si cercheranno due soluzioni particolari linearmente indipendenti, del
tipo:
u1 (z) = u1 (1/t) = t1
ck tk = z 1
k=0
ck z k ,
c0 6= 0
k=0
e

u2 (1/t) =
P
k
6 n , d0 6= 0
t2
k=0 dk t
P 1k 2 =
2
1 2 = n , b0 6= 0
au1 (1/t) ln t + t
k=0 bk t
(2.62)
ovvero

u2 (z) =
P
k
1 2 =
6 n , d0 6= 0
z 2
k=0 dk z
P
k
b
z
1 2 = n , b0 6= 0 .
au1 (z) ln z + z 2
k=0 k
(2.63)
Gli esponenti 1 e 2 sono le soluzioni dellequazione indiciale relativa allequazione

differenziale (2.55), cioè
2 + (p0 1) + q0 = 0 ,
dove

2 P (1/t)
P (1/t)
= lim t
= 2 lim
2
t0
t0
t
t
t
Q(1/t)
Q(1/t)
= lim t2
= lim
.
4
t0
t0
t
t2

p0
q0
In termini della variabile z:

p0 = 2 lim zP (z)
q0 =
z
2
lim z Q(z) .
Detti ora
p0 =
q0 =
lim zP (z)
lim z 2 Q(z) .
82
lequazione indiciale diventa

2 + (1 p0 ) + q0 = 0 .
Si noti il cambiamento di segno nel termine lineare rispetto allequazione indiciale per
singolarità al finito. I coefficienti ck , dk , bk e a nelle equazioni (2.62) si ottengono
sostituendo le soluzioni nellequazione differenziale e usando gli sviluppi
P (z) =
1 X pn
z n=0 z n
1 X qn
Q(z) = 2
.
z n=0 z n
(2.64)
(2.65)
Le serie (2.62) convergono certamente allesterno di un cerchio centrato nellorigine e

che comprende al suo interno tutte le singolarità dellequazione differenziale.
2.1.5
Esempi
Consideriamo le equazioni (2.3-2.10) e studiamone il comportamento per z .

1) Equazione delloscillatore armonico semplice:
u00 + 2 u = 0
2
2
, Q(1/t)
= 4
P (1/t) =
t
t
Allinfinito lequazione ha una singolarità irregolare.
2) Equazione di Legendre:
(1 z 2 )u00 2zu0 + u = 0
2
2/t
P (1/t) = 2
, Q(1/t)
= 4
t t 1
t t2
Allinfinito lequazione ha una singolarità fuchsiana.
3) Equazione di Bessel:
z 2 u00 + zu0 + (z 2 2 )u = 0
2 1
1
1 2 t2
P (1/t) = =
, Q(1/t)
=
t
t
t
t4
83
4) Equazione di Laguerre
zu00 + (1 z)u0 + au = 0
2 t1
a
P (1/t) = 2 , Q(1/t)
= 3
t
t
t
5) Equazione di Hermite:
u00 2zu0 + 2u = 0
2
2
2
= 4
P (1/t) = + 3 , Q(1/t)
t t
t
6) Equazione di Chebyshev:
(1 z 2 )u00 zu0 + n2 u = 0
1
n2
2
, Q(1/t)
= 2 2
P (1/t) = + 2
t t(t 1)
t (t 1)
7) Equazione ipergeometrica confluente:
zu00 + (c z)u0 au = 0
2 ct 1
a
P (1/t) =
, Q(1/t)
= 3
2
t
t
t
8) Equazione ipergeometrica:
z(z 1)u00 + [(1 + a + b)z c]u0 + abu = 0
2 1 + a + b ct
ab
P (1/t) =
, Q(1/t)
= 2
t
t(1 t)
t (1 t)
84
Lequazione ipergeometrica è un esempio di equazione totalmente fuchsiana con tre

punti singolari in 0, 1, i cui indici valgono rispettivamente (0, 1 c), (0, c a b),
(a, b); questa informazione si racchiude nel simbolo P di Riemann
0
0
a z
P =
.
1c cab b
Per c 6= 0, 1, 2, lequazione ipergeometrica ammette una soluzione regolare
nellorigine detta funzione ipergeometrica:
F (a, b; c; z) =
X
(a)l (b)l
l=0
l!(c)l
zl
(2.66)
dove (a)l a(a + 1) (a + l 1). Lequazione ipergeometrica è particolarmente

importante perche ogni equazione totalmente fuchsiana con 3 punti singolari può essere
ad essa ricondotta con un cambiamento di variabile indipendente del tipo w = z+
e
z+
a meno di potenze a fattore della soluzione.
85
Parte II
Introduzione allanalisi armonica
86
Capitolo 3
Serie di Fourier
3.1
Introduzione
Per presentare largomento della seconda parte di questo corso iniziamo a discutere
un problema fisico molto semplice ma molto significativo. Consideriamo il circuito
oscillante in figura 3.1.
Ricordando che le tensioni ai capi di un condensatore di capacità C, di una resistenza R e di una bobina di induttanza L sono date rispettivamente da q(t)/C, Ri(t)
e Ldi/dt, dove q(t) è la carica su una faccia del condensatore e i(t) = dq/dt la corrente, si ricava facilmente che la tensione in uscita u(t) è legata a quella in entrata f (t)
dallequazione:
du
d2 u
+ RC
+ u(t) = f (t) ,
2
dt
dt
che può essere utile scrivere nella forma
LC
(3.1)
Lt u(t) = f (t) ,
dove Lt è loperatore differenziale
Lt = LC
d2
d
+ RC + 1 .
2
dt
dt
(3.2)
Se la tensione dingresso è sinusoidale

f (t) = V cos(t)
(3.3)
è facile trovare una soluzione dellequazione (3.1) della forma

u(t) = A cos(t) + B sin(t) .
87
(3.4)
Figura 3.1: Circuito oscillante.
Basta infatti sostituire (3.4) nelleq.(3.1) per ottenere il sistema di due equazioni
2 LCA + RCB + A = V
2 LCB RCA + B = 0
(3.5)
(3.6)
nelle due incognite A e B e ricavare quindi

1 2 LC
A=
V ;
(1 2 LC)2 + 2 R2 C 2
B=
(1
RC
V .
+ 2 R2 C 2
2 LC)2
(3.7)
Ponendo A = cos e B = sin , la (3.4) diventa:

u(t) = cos(t ) = Re ei(t) .
Ne segue che, introducendo lampiezza complessa
U = A iB = ei ,
(3.8)
la (3.4) può anche scriversi nella forma:

u(t) = Re(U eit ) .
(3.9)
La soluzione generale delleq. (3.1) si ottiene aggiungendo alla soluzione particolare

(3.4) la soluzione generale dellequazione omogenea associata alleq. (3.1); noi tuttavia
non ci occupiamo del transitorio, che tende esponenzialmente a zero per t e che è
quindi trascurabile per tempi t molto pi`
u grandi della costante di tempo L/R; in altre
u elegantemente,
parole, ci interessano solo le soluzioni periodiche delleq. (3.1). Pi`
possiamo dire che noi intendiamo risolvere lequazione:
Lu = f ,
88
(3.10)
dove loperatore L è definito non solo dalla sua espressione differenziale Lt di eq.(3.2),
ma anche dal suo dominio, specificato dalle condizioni al contorno periodiche:
u(t + T ) = u(t) ,
dove T = 2/ è il periodo della funzione termine noto f (t).
Leq. (3.1) gode di una proprietà molto importante, la linearit`
a:
Propriet`
a P1: se u1 (t) è soluzione delleq. (3.1) con termine noto f1 (t):
LC
d2 u1
du1
+ u1 (t) = f1 (t)
+ RC
2
dt
dt
(3.11)
e analogamente u2 (t) soddisfa:

LC
d2 u2
du2
+ u2 (t) = f2 (t)
+ RC
2
dt
dt
(3.12)
allora
u(t) = 1 u1 (t) + 2 u2 (t)
(3.13)
sarà soluzione delleq. (3.1) con

f (t) = 1 f1 (t) + 2 f2 (t) ,
(3.14)
dove 1 e 2 sono costanti complesse qualsiasi.

Per dimostrare la proprietà P1 basta moltiplicare per 1 leq. (3.11) e per 2
leq. (3.12) e poi sommarle. La linearità delleq. (3.1) può anche esprimersi dicendo che
loperatore differenziale Lt (3.2) è lineare, cioè che vale:
Lt [1 u1 (t) + 2 u2 (t)] = 1 Lt u1 (t) + 2 Lt u2 (t) .
(3.15)
La P1 suggerisce per esempio di scegliere come termine noto anziche la f (t) delleq. (3.3)
la
e (t) = V eit ,
(3.16)
dal momento che cos t = (eit + eit )/2. Allora la soluzione dellequazione
du
d2 u
LC 2 + RC
+ u(t) = e (t)
dt
dt
(3.17)
è ancora pi`
u immediata; infatti se si pone1
u(t) = U ()eit u (t) ,
1
(3.18)
Sia nelle incognite U e u(t) che nella funzione di entrata e(t) si è esplicitamente ricordata la
dipendenza dal parametro per motivi che diventeranno chiari fra poco (vedi eq. (3.23).
89
sostituendo (3.18) nelleq. (3.17) si ottiene:

LC(i)2 + RCi + 1 U ()eit = V eit ,
(3.19)
da cui
U () =
V
1
2 LC
+ iRC
V
ZC ,
R + ZL + ZC
(3.20)
dove
1
(3.21)
iC
sono le impedenze complesse dellinduttanza L e della capacità C. Usando poi lidentità
di Eulero
ZL () = iL ,
cos(t) =
ZC () =
eit + eit
= Re eit ,
2
(3.22)
si può scrivere la f (t) di (3.3) come

1
[e (t) + e (t)] = Re e (t)
per V reale .
(3.23)
2
Allora la proprietà di linearità P1 ci permette di scrivere subito la soluzione delleq. (3.1)
nella forma
f (t) =
1
[u (t) + u (t)] = Re u (t)
(3.24)
2
con u (t) dato dalleq. (3.18)2 . Osservando che la U () delleq. (3.20) coincide con la
U di (3.8), si verifica subito che la soluzione (3.24) coincide con la soluzione (3.4).
La proprietà P1 ha anche una conseguenza molto pi`
u importante e generale:
u(t) =
Propriet`
a P2: se in qualche modo si riesce a scrivere il termine noto delleq. (3.1)
nella forma
f (t) =
Vn ein t ,
(3.25)
allora sarà possibile risolvere leq. (3.1) e la soluzione sarà (sempre a meno del transitorio) ancora della stessa forma:
u(t) =
Un ein t
(3.26)
Si è usata lidentit`
a U () = U () che segue dalla (3.20); con indichiamo la complessa
coniugazione.
2
90
con i coefficienti Un dati dalla

Un =
Vn
Vn
ZC (n ) =
.
R + ZL (n ) + ZC (n )
1 n2 LC + in RC
(3.27)
Lo scopo della seconda parte di questo corso sarà proprio di trovare il modo per
esprimere ampie classi di funzioni f (t) nella forma (3.25); in linguaggio matematico
possiamo dire che oggetto di questo corso è la analisi armonica.
Lanalisi armonica è uno strumento matematico di enorme importanza in fisica: infatti le equazioni differenziali lineari del secondo ordine (in particolare a coefficienti costanti ) intervengono in numerosissimi problemi in ogni campo della fisica,
non solo per i circuiti RLC, ma ogni volta che si vogliano descrivere piccole oscillazioni attorno a una situazione di equilibrio. Inoltre lanalisi armonica, in particolare la
trasformata di Fourier, gioca un ruolo fondamentale nella Meccanica Quantistica.
` utile osservare che alla radice del ruolo privilegiato che giocano le funzioni eit
E
nella soluzione delleq. (3.1) sta il fatto che esse sono autofunzioni delloperatore
differenziale Lt , cioè che vale lequazione
Lt eit = eit ,
(3.28)
= LC 2 + iRC + 1 .
(3.29)
dove lautovalore vale
Nel capitolo successivo vedremo che ulteriori specificazioni delloperatore differenziale

(le condizioni al contorno) faranno s` che non tutti i valori di siano ammissibili.
91
3.2
Funzioni periodiche e serie di Fourier
Una prima classe di funzioni per cui si può effettuare lanalisi armonica (3.25) contiene
le funzioni periodiche (di periodo T ), tali cioè che
f (t + T ) = f (t),
t R .
(3.30)
In tal caso è lecito aspettarsi che gli n delleq. (3.25) siano tutti multipli interi
dellarmonica fondamentale = 2
ovvero
T
n = n
2
,n Z .
T
(3.31)
Infatti dallidentità
e2in cos(2n) + i sin(2n) = 1 , n Z
(3.32)
segue
2
ein (t+T ) ei T
n(t+T )
= ei T
nt i2n
= ein t ,
n Z , t R .
(3.33)
Notare che la scelta di condizioni al contorno periodiche seleziona fra le possibili

soluzioni dellequazione agli autovalori (3.28) solo quelle con n dato dalla (3.31). Senza
occuparci per il momento di discutere la convergenza della serie
f (t) =
Vn ein t
(3.34)
nZ
detta serie trigonometrica di Fourier, vediamo subito come si possono calcolare

i coefficienti Vn nota la f (t); basta moltiplicare ambo i membri per eil t (l Z) e
integrare su un periodo (supponendo di poter integrare termine a termine) per ottenere:
Z T
X Z T
il t
e
f (t)dt =
Vn
eil t ein t dt.
(3.35)
0
Con il cambio di variabile x =
nZ
2t
,
T
suggerito dalla (3.31), e usando lidentità
ei(nl)x dx = 2nl ,
(3.36)
detta relazione di ortogonalit`

a, si ottiene subito:
1
Vl =
T
92
eil t f (t)dt .
(3.37)
Leq. (3.37) è di grande importanza perche ci fornisce i coefficienti della serie di

Fourier (3.34) e quindi, grazie alla proprietà P1, il modo per risolvere lequazione
differenziale (3.1) con termine noto f (t) periodico. I passi sono i seguenti:
con la (3.37) si calcolano i coefficienti Vl della serie di Fourier (3.34) di f (t);
per ognuna delle armoniche, cioè per ognuno dei termini di tale serie, si applica
la procedura che ci ha portato dal termine noto (3.16) alla soluzione (3.18),
ottenendo cos` i coefficienti Un dati dalla (3.27);
la soluzione dellequazione differenziale sarà data perciò dalla serie di Fourier
(3.26).
3.2.1
Convergenza puntuale delle serie trigonometriche

di Fourier
Vogliamo ora discutere le proprietà di convergenza della serie

X
an einx
(3.38)
nZ
con i coefficienti an dati da

1
an =
2
einx f (x) dx .
(3.39)
Per comodità siamo passati alla variabile x = 2t

, in cui il periodo è 2. Notiamo subito
T
che per f (x) periodica, f (x + 2) = f (x), lintegrale può essere esteso a qualsiasi altro
intervallo di ampiezza 2 3 . Notiamo inoltre che affinche la (3.39) abbia senso bisogna
che lintegrale esista; ciò succede certamente se f (x) è sommabile4 , perche tale rimane
dopo essere stata moltiplicata per il fattore einx , il cui modulo vale 1.
3
Vale infatti lidentit`

a, x0 R
Z x0 +2
g(x)dx =
Z
x0
x0 +2
+
x0
g(x)dx
g(x)dx +
x0
g(x)dx +
x0
g(y + 2)dy
(3.40)
(nellultimo integrale si è effettuato il cambio di variabile x = y + 2). Se lintegrando g(x) è una

funzione periodica di periodo 2, allora i due ultimi addendi della (3.40) si cancellano e vale
Z x0 +2
Z
g(x)dx =
g(x)dx,
x0 .
(3.41)
x0
Le due scelte pi`

u consuete sono x0 = oppure x0 = 0.
4
Definiremo pi`
u avanti che cosa significa funzione sommabile; per ora può essere tranquillamente
letto come sinonimo di funzione assolutamente integrabile.
93
Prima di considerare la convergenza della serie conviene considerare il seguente

Lemma di Riemann: Qualunque sia lintervallo (a, b), finito o infinito, per ogni f (x)
sommabile, vale
Z b
Z b
Z b
ikx
lim
f (x)e
dx = lim
f (x) cos(kx) dx = lim
f (x) sin(kx) dx = 0 .
k
(3.42)
Dimostrazione: ci limitiamo a dimostrare tale lemma nel caso particolare in cui f (x)
sia di classe C 1 , cioè continua con la sua derivata prima. In tal caso è lecito integrare
per parti e si ha
Z
Z b

1 b 0
1
ikx
ikx b
f (x) e
f (x)eikx dx ,
(3.43)
f (x)e
dx =
a
ik
ik a
a
da cui

Z b

Z b

1
0
ikx

f (x)e
dx
|f (x)| dx .
|f (b)| + |f (a)| +

|k|
a
a
(3.44)
La parentesi graffa non contiene pi`

u termini dipendenti da k, quindi il secondo membro
tende a zero per k , e di conseguenza anche il primo. Usando le formule di Eulero
si completa la dimostrazione, nel caso particolare di funzioni di classe C 1 ; nel caso
generale la dimostrazione prosegue usando ilR fatto che per ogni funzione sommabile
b
f (x) e ogni > 0 esiste una g C 1 tale che a |f (x) g(x)|dx < .
[q.e.d.]
TEOREMA: Condizione sufficiente affinche la serie (3.38) converga puntualmente
a f (x0 ) è che la funzione f (x) (sommabile nellintervallo (0, 2)) sia di classe C 1
nellintorno del punto x0 .
Dimostrazione Definiamo la ridotta N-esima della serie (3.38) come
N
X
SN (x0 ) =
al eilx0 .
(3.45)
l=N
Sostituendo in (3.45) la definizione (3.39) dei coefficienti al si ottiene

1
SN (x0 ) =
2
dyf (y)
0
N
X
eil(x0 y) .
(3.46)
l=N
Usando lidentità
N
X
il
iN
= e
2N
X
ei
n
n=0
l=N
= eiN
i2(N +1/2)
1 ei(2N +1)
i(N +1/2) 1 e
=
e
1 ei
ei/2 ei/2
sin(N + 1/2)
sin /2
(3.47)
94
la (3.45) diventa
Z 2
1
sin [(N + 1/2)(x0 y)]
SN (x0 ) =
dyf (y)
2 0
sin(x0 y)/2
Z
1
sin [(N + 1/2)t]
=
dtf (x0 + t)
,
2
sin t/2
(3.48)
dove nellultimo passaggio si è posto t = x0 y e si è usata la periodicità dellintegrando

per fissare lintervallo di integrazione5 . Notare che per calcolare limN SN (x0 ) non
si possono applicare direttamente le formule di Riemann (3.42) alla (3.48), poiche la
(x0 +t)
0)
non è integrabile: essa diverge come 2f (x
per t 06 .
funzione fsin
t/2
t
Si noti che lintegrale si pu`
o spezzare come segue:
Z
Z 1 Z 2 Z
=
+
+
(3.49)
1 , 2 (0, 2) si possono applicare le formule di Riemann al primo e terzo integrale, quindi

Z
Z 2
sin [(N + 1/2)t]
sin [(N + 1/2)t]
f (x0 + t)
lim
f (x0 + t)
dt = lim
dt ;
(3.50)
N
N
sin t/2
sin t/2
1
perci`
o la somma della serie nel punto x0 dipende solo dal comportamento locale della funzione
f (x) (sommabile in (0, 2)) in un intorno (arbitrariamente piccolo) del punto x0 . Usando lidentit`
a
1
2
sin [(N + 1/2)t]

dt = 1 ,
sin t/2
(3.51)
che si può verificare direttamente con il metodo dei residui (o anche considerando il
caso particolare della (3.48) per f (x) = 1), si può scrivere

Z
1
1
f (x0 + t) f (x0 )
dt sin N +
SN (x0 ) f (x0 ) =
t
.
(3.52)
2
2
sin t/2
Adesso la funzione che moltiplica sin(N + 1/2)t è sommabile nellintervallo (, );
infatti in tale intervallo sin t/2 si annulla solo nellorigine e
f (x0 + t) f (x0 )
= 2f 0 (x0 ) .
t0
sin t/2
lim
(3.53)
Si può quindi passare al limite per N e applicare le formule di Riemann per

ottenere
f (x0 ) = lim SN (x0 ) .
N
(3.54)
[q.e.d]
sin(N + 1/2)t e sin t/2 hanno periodo 4, ma il loro rapporto ha periodo 2.

è integrabile se f (x0 ) = 0, allora si applicano le formule di Riemann e si ottiene correttamente
limN SN (x0 ) = 0.
6
95
Notare che il Teorema (3.54) può essere esteso al caso in cui nel punto x0 la funzione
abbia una discontinuità di I specie, ma sia di classe C 1 sia in un intorno sinistro che in
un intorno destro di x0 . Al posto della (3.52) si scrive infatti

Z 0
1
1
f (x0 + t) f (x0 )
f (x0 +) + f (x0 )
=
dt sin N +
t
SN (x0 )
2
2
2
sin t/2

Z
1
1
f (x0 + t) f (x0 +)
+
dt sin N +
t
,
2 0
2
sin t/2
(3.55)
dove f (x0 ) e f (x0 +) sono i limiti destro e sinistro nel punto x0 e si è usata lidentità
Z 0
Z
sin (N + 1/2) t
sin (N + 1/2) t
dt =
dt = .
(3.56)
sin t/2
sin t/2
0
Dalle formule di Riemann segue allora:
lim SN (x0 ) =
f (x0 +) + f (x0 )
,
2
(3.57)
di cui la (3.54) è ovviamente un caso particolare. Se il punto x0 cade in uno degli

estremi dellintervallo di definizione della f (x0 ), continua a valere la (3.57) purche la
funzione sia continuata periodicamente: f (x + 2) = f (x).
3.2.2
Importanti commenti
Il teorema appena visto garantisce la convergenza puntuale, mentre per la convergenza

uniforme in [a, b] è necessario che la funzione f (x) sia ivi continua e la sua derivata
continua a tratti. Questione diversa è lintegrabilità termine a termine della serie di
Fourier, che è invece garantita sotto le ipotesi del teorema precedente. Notiamo infatti
che lintegrale della serie ha coefficienti soppressi da k
Z x
+
X
ak ikx x
e |x0
f (x) dx =
(3.58)
ik
x0
k=
e converge pi`
u rapidamente.
In vista di una successiva applicazione fisica, torniamo alla variabile t =
f (t) =
an ein t
T
x:
2
(3.59)
n=
i cui coefficienti sono, secono la (3.37),

Z
1 T /2
am =
f (t)eim t dt .
T T /2
96
(3.60)
Se la funzione f (t) assume valori reali si vede subito che i coefficienti an soddisfano la
relazione seguente:
an = an .
` utile osservare che, posto an = |an |e
E
in
(3.61)
, si può allora scrivere:

an ein t + an ein t = |an | ei(n tn ) + ei(n tn )
= An cos(n t n ),
(3.62)
con
An = 2|an | .
(3.63)
Quindi per f (t) reale la serie (3.59) diventa
f (t) = a0 +
An cos(n t n ) .
(3.64)
n=1
Un integrale dalle importanti applicazioni fisiche è

!
Z
Z T /2
X
X
1
1 T /2
2
im t
|am |2 .
|f (t)| dt =
a e
f (t) =
T T /2
T T /2 m= m
m=
dove si sono sfruttate le identità:
Z T /2
eim t ein t = T mn ,
(3.65)
(3.66)
T /2
note come relazioni di ortogonalit`

a, su cui torneremo pi`
u avanti. La (3.65), che
prende il nome di equazione di Parseval scritta nella base delle funzioni esponenziali,
illustra come ogni componente di Fourier contribuisca separatamente allintegrale; non
ci sono cioè termini di interferenza del tipo am an . Qualora f (t) rappresenti la corrente
elettrica attraverso una resistenza R, la (3.65) moltiplicata per R mostra che la potenza
media dissipata per effetto Joule è uguale alla somma delle potenze dissipate sulle varie
frequenze.
Per f (t) reale la (3.65) diventa infatti:
2
Z

X
X
An
1 T /2
2
2
2
2
|f (t)| dt = a0 + 2
|an | = a0 +
,
(3.67)
T T /2
2
n=1
n=1
dove nellultimo passaggio si è ricordata
la (3.63). In questo caso a0 è la componente
di corrente continua della f (t) e An / 2 il valore efficace della corrente alternata di
pulsazione n .
97
Serie di Fourier e funzioni trigonometriche.

Molto spesso anziche usare il si
stema trigonometrico in forma esponenziale eilx , l Z è utile usare il sistema
trigonometrico tout court:
{1, sin x, cos x, sin(2x), cos(2x), ...} = {sin(nx), cos(nx)} ,
n = 0, 1, 2, ...
(3.68)
` immediato verificare direttamente che le (3.68) formano un sistema di funzioni

E
ortogonali nellintervallo (, ) (o in qualunque altro intervallo di ampiezza 2).
Infatti:
Z
sin(mx) sin(nx)dx = mn (m 6= 0)
Z
cos(mx) cos(nx)dx = mn
(m 6= 0)
= 2mn
Z
(m = 0)
sin(mx) cos(nx)dx = 0 .
(3.69)
Data una f (x) sommabile nellintervallo (, ), anziche la serie (3.38) proviamo a

scrivere
f (x) =
A0 X
[An cos(nx) + Bn sin(nx)] .
+
2
n=1
(3.70)
Se la (3.70) è vera, i coefficienti An e Bn si ottengono moltiplicando la (3.70) rispettivamente per cos(mx) e sin(mx), integrando su x fra e e sfruttando le relazioni di
ortogonalità (3.69), nellipotesi che la serie converga a f (x) e si possa integrare termine
a termine. Si ricava cos`
An
Bn
Z
1
cos(nx)f (x)dx
=

Z
1
=
sin(nx)f (x)dx

n = 0, 1, . . .
(3.71)
n = 1, 2, . . . .
(3.72)
Il coefficiente A0 si ricava integrando la (3.70) tra e .

Per la convergenza puntuale della serie (3.70) valgono esattamente gli stessi teoremi
dimostrati per la serie (3.38), cioè se in un punto x interno allintervallo (, ) la
funzione f (x) è di classe C 1 , cioè continua assieme alla sua derivata prima, allora la
(3.70) è vera nel senso della convergenza puntuale. Se invece nel punto x0 la f (x)
98
ha una discontinuità di prima specie, ma è di classe C 1 sia in un intorno sinistro che

in un intorno destro di x0 , vale allora:
A0 X
f (x0 +) + f (x0 )
+
[An cos(nx0 ) + Bn sin(nx0 )] =
.
2
2
n=1
(3.73)
dove f (x0 +) e f (x0 ) sono rispettivamente i limiti destro e sinistro di f (x) nel punto
x0 .
Se il punto x0 cade in uno degli estremi dellintervallo di definizione continua a
valere quanto abbiamo detto per i punti interni purche la funzione sia continuata
periodicamente su tutto lasse reale secondo la f (x + 2) = f (x).
Lequazione di Parseval, in termini dei coefficienti An e Bn , assume la forma
seguente:
X
2
|An |2 + |Bn |2 .
|f (x)| dx = |A0 | +
2
n=1
Se la funzione f (x) è pari (f (x) = f (x)), i coefficienti Bn sono nulli e la (3.70) si

riduce a una serie di coseni:
A0 X
f (x) =
+
An cos(nx) .
2
n=1
Se la funzione f (x) è dispari (f (x) = f (x)), i coefficienti An sono nulli e la (3.70)
si riduce a una serie di seni:
f (x) =
Bn sin(nx) .
n=1
Esempio
Sviluppiamo in serie di Fourier la funzione a gradino:

1
<x<0
(x) =
+1
0x .
99
Poiche la funzione è dispari, lo sviluppo in serie di Fourier conterrà solo seni (An = 0).
I coefficienti Bn sono:
Z
1
Bn =
sin(nx)(x)dx

Z 0

Z
1
sin(nx)dx +
sin(nx)dx
=
0
Z
2
2
=
sin(nx)dx =
[ cos(nx)]0
0
n

2
0
n pari
n
=
[1 (1) ] =
4
n dispari.
n
n
Pertanto
4
4
4
sin x +
sin(3x) +
sin(5x) + ...
3
5
4 X sin[(2n + 1)x]
=
.
n=0
2n + 1
(x) =
Graficamente, le successive approssimazioni sono mostrate in Fig. (3.2 Notare che

nellorigine e nei punti la somma della serie vale zero, in accordo con la (3.73).
Per n , si verifica il cosiddetto fenomeno di Gibbs, tipico delle serie di Fourier di funzioni
discontinue: in prossimit`
a di una discontinuità di prima specie, la ridotta ennesima della serie di
Fourier presenta un picco, tanto pi`
u stretto quanto maggiore è n, di altezza pari a circa il 9% del
salto. Questo esempio mostra che la serie di Fourier non converge uniformemente alla funzione f (x)
nei punti in cui essa ha una discontinuit`
a di I specie.
Scriviamo adesso esplicitamente la generalizzazione delle eq.(3.70),(3.71) e (3.72) al

caso di funzioni periodiche con periodo T 6= 2. In questo caso si pone
t=
T
x.
2
Se x varia nellintervallo (, ), la variabile t varierà nellintervallo (T /2, T /2). La

serie di Fourier (3.70) diventa cos`:

A0 X
2nt
2nt
f (t) =
+
An cos
+ Bn sin
,
2
T
T
n=1
con
An
Bn

Z
2 T /2
2nt
cos(nx)f (t(x))dx =
cos
f (t)dt
T T /2
T

Z
Z
1
2 T /2
2nt
=
sin(nx)f (t(x))dx =
sin
f (t)dt .

T T /2
T
1
=
100
1.5
N=0
N=1
N=2
N=100
1
0.5
-0.5
-1
-1.5
-4
-3
-2
-1
Figura 3.2:
3.2.3
Altri esempi
Esempio 1. Lo sviluppo in serie di Fourier della funzione

f (x) = | sin x|
<x<
contiene solo coseni, perche f (x) è pari:
f (x) =
A0 X
+
An cos(nx) .
2
n=1
I coefficienti di Fourier sono:

Z
Z
1
2
| sin x| cos(nx)dx =
sin x cos(nx)dx
An =

0

Z
1
1 1 cos(n + 1) cos(n 1) 1
=
[sin(n + 1)x sin(n 1)x] dx =
+
0
n+1
n1
2 1 + cos n
=
se n 6= 1 .
n2 1
101
Se n = 1
A1
2
=

2 sin2 x
sin x cos xdx =
=0.
2 0
Quindi
2 X 1 + (1)n
2
4 X cos 2kx
2
cos nx =
f (x) =
n=2 n2 1
k=1 (2k)2 1

2
4 cos 2x cos 4x cos 6x
=
+
+
+

3
24
35
Esempio 2.
f (x) = |x|
A0
An
<x<
Z
2
=
xdx =
0

Z
Z
2
2
sin nxdx
x cos nxdx =
x sin nx|0
=
0
n
0
2
=
[1 (1)n ]
n2
da cui

4 X cos(2k + 1)x
4
cos 3x cos 5x
f (x) =
=
cos x +
+
+
2 k=0 (2k + 1)2
2
9
25
Esempio 3.
f (x) = x
A0
An
Bn
0 < x < 2
Z
1 2
=
xdx = 2
0

Z
Z 2
1 2
1
2
=
x cos nx =
x sin nx|0
sin nxdx = 0
0
n
0

Z
Z 2
1 2
1
2
2
x sin nx =
x cos nx|0 +
cos nxdx = .
=
0
n
n
0
Quindi
f (x) = 2
X
sin nx
n=1

sin 2x sin 3x
= 2 sin x +
+
+
2
3
Si consiglia agli studenti di disegnare i grafici delle f (x) dei tre esempi proposti.
102
Capitolo 4
Trasformate Integrali
4.1
Trasformata di Fourier
Nel capitolo 1 abbiamo imparato a risolvere leq. (3.1) nel caso in cui il termine noto
sia periodico e abbiamo dovuto poi allargare il discorso per dare un senso preciso
alla convergenza puntuale della serie (trigonometrica) di Fourier. Ma come possiamo
risolvere leq. (3.1) quando il termine noto non è periodico? Come fare a scriverlo come
somma di termini della forma ein t che sono invece periodici?
Il modo migliore è di considerare una funzione non periodica come caso limite di
una periodica con periodo L che tende a infinito. A questo scopo scriviamo i coefficienti
di Fourier
Z
2n
1 L/2 ikn x
e
f (x)dx , con kn
(4.1)
an =
L L/2
L
nella forma
1
an = F (kn )k ,
2
dove
1
F (kn )
2
(4.2)
L/2
eikn x f (x)dx ,
(4.3)
L/2
e
k = kn kn1 =
2
.
L
(4.4)
Allora la serie trigonometrica di Fourier (3.59) si può riscrivere come
1 X ikn x
f (x) =
e F (kn )k .
2 n=
(4.5)
A questo punto effettuiamo il limite L , nel quale k 0 e la serie a secondo

membro della (4.5) può riguardarsi come una somma integrale alla Riemann per la
103
funzione eikx F (k), estesa allintervallo (, +), diviso in infiniti intervalli parziali
di ampiezza k 0. Pertanto la (4.5) diventa
1
f (x) =
2
1
F (k) =
2
F (k)eikx dk ,
(4.6)
f (x)eikx dx .
(4.7)
dove, per la (4.3),
La funzione F (k) si chiama trasformata di Fourier (TF) della funzione f (x) e la

funzione f (x) antitrasformata di Fourier della funzione F (k).
I passaggi che abbiamo fatto per giungere alle (4.6) e (4.7) sono un po disinvolti.
Per essere precisi dobbiamo dimenticare le operazioni di limite e assumere la (4.7) come
definizione della Trasformata di Fourier, sotto la condizione che la f (x) sia sommabile
sullasse reale; la (4.6) va invece scritta pi`
u correttamente come
Z R
1
F (k)eikx dk
(4.8)
lim
f (x) =
2 R R
e si dimostra che vale sotto la condizione (sufficiente) che nellintorno del punto x
la funzione f (x) sia di classe C 1 ; naturalmente se lintegrale (4.6) esiste la (4.8) è
equivalente alla (4.6).
4.1.1
Esempi
Esempio 1. La trasformata di Fourier della funzione

f (x) =
x2
1
,
+ a2
a R+
è
1
F (k) =
2
eikx
dx .
(x + ia)(x ia)
Se k > 0 chiudiamo il cammino di integrazione nel semipiano Imz < 0 e otteniamo:

r ka

1
e
eikz

F (k) = 2i Res
=
.

(z + ia)(z ia) z=ia
2 a
2
104
Se k < 0 chiudiamo invece nel semipiano Imz > 0 e otteniamo:

r ka

1
e
eikz

F (k) = 2i Res
=
.

(z + ia)(z ia) z=ia
2 a
2
Pertanto
r
F (k) =
e|k|a
.
2 a
La verifica della (4.6) è immediata e coinvolge solo un integrale elementare.

Esempio 2. La trasformata di Fourier della funzione

1
|x| < a
f (x) =
0
|x| > a
(4.9)
è
Z +
Z +a
1
1
ikx
f (x)e
dx =
eikx dx
F (k) =
2
2 a
r
1 eika eika
2 sin ka
=
=
.
ik
k
2
(4.10)
Per verificare che lantitrasformata di (4.10) sia effettivamente la (4.9) dobbiamo calcolare lintegrale
Z +R
Z +R
1
1
sin(ka) ikx
ikx
I(x)
F (k)e dk =
lim
lim
e dk
R R
k
2 R R

Z +R ik(x+a)
eik(xa)
e
1
lim
dk
(4.11)
=
2i R R
k
k
deformando il cammino di integrazione come nellesempio 4 del paragrafo 1.6.4, aggirando per esempio lorigine nel semipiano immaginario positivo. Se x > a entrambi
gli integrali in (4.11) ricadono nel caso > 0 del lemma di Jordan e pertanto si può
chiudere il cammino dintegrazione con una semicirconferenza nel semipiano superiore;
allinterno del cammino dintegrazione lintegrando è regolare e pertanto:
I(x) = 0 se x > a .
Analogamente, se x < a conviene aggirare lorigine nel semipiano immaginario negativo, perche entrambi gli integrali in (4.11) ricadono nel caso < 0 del lemma di
Jordan e pertanto si può chiudere il cammino dintegrazione con una semicirconferenza
nel semipiano inferiore, ottenendo di nuovo zero:
I(x) = 0 se x < a .
105
Se invece |x| < a, e si aggira lorigine nel semipiano immaginario positivo, il primo
integrale, che si può chiudere nel semipiano superiore, dà zero, mentre il secondo, che
va chiuso nel semipiano inferiore, dà:

eik(xa)
1
(2i)Res
=1.
I(x) =
2i
k k=0
Abbiamo cos` dimostrato che I(x) = f (x) per ogni x 6= a. Per x = a possiamo
deformare il cammino sopra la singolarità scrivendo, come in Sez. 1.6.5,

Z +R 2ika
e
1
1
I(a) =
lim
dk .
(4.12)
2i R R k + i k + i
Il primo integrale si chiude nel semipiano immaginario positivo e si annulla per il teorema di Cauchy, mentre per il secondo usiamo la (1.61) sapendo che la parte principale
dellintegrale si annulla. Il risultato è
I(a) =
1
1
f (a) + f (a+)
(i) = =
.
2i
2
2
(4.13)
Questo mostra che nei punti di discontinuità di prima specie la situazione è analoga a
quella vista per le serie di Fourier: lantitrasformata dà il valor medio tra i limiti destro
e sinistro della funzione.
Esempio 3. La trasformata di Fourier della funzione gaussiana
f (x) = ex
2 /a2
è
Z +
Z
2
e(ka) /4 + (x/aika/2)2
1
(x/a)2 ikx
e
dx =
e
dx
F (k) =
2
2
Z +
2
e(ka) /4
a
2
2 2
=
a
et dt = ea k /4 ,
2
2
(4.14)
cioè ancora una gaussiana, di larghezza inversamente proporzionale a quella della funzione trasformanda. 1 La verifica della (4.6), che dà lantitrasformata di F (k), è
immediata: non si tratta che di rifare lo stesso conto con a A = 2/a.
1
Con il cambiamento di variabile x t = x/a ika/2, il cammino di integrazione nel piano

complesso di t non è pi`
u lasse reale ma è diventato una retta ad esso parallela; tuttavia, usando la
teoria dellintegrazione in campo complesso, è facile mostrare che ciò non fa differenza.
106
4.1.2
Propriet`
a della trasformata di Fourier
Elenchiamo alcune importanti proprietà delle trasformate di Fourier. Per comodità

introduciamo il simbolo Fk (f ) per indicare la trasformata di Fourier della funzione
f (x):
Z +
1
f (x)eikx dx .
(4.15)
Fk (f )
2
Linearit`
a:
Fk (a1 f1 + a2 f2 ) = a1 Fk (f1 ) + a2 Fk (f2 ) , a1 , a2 C .
(4.16)
Trasformata di Fourier di funzioni a parit`

a definita
Cos` come la serie trigonometrica di Fourier di una funzione pari (dispari) contiene
solo coseni (seni), le trasformate di Fourier di funzioni a parità definita si possono
semplificare come segue:
Z
1
f (x)[cos(kx) i sin(kx)]dx
F (k) =
2
r R
2
f (x) cos(kx)dx
se f (x) = f (x)
0R
=
se f (x) = f (x) .
i 0 f (x) sin(kx)dx
(4.17)
Una ovvia conseguenza delle (4.17) è che la Trasformata di Fourier di una funzione
pari (dispari) è una funzione pari (dispari).
La trasformata di Fourier della derivata f 0 (x) (ammesso che f 0 (x) esista e
sia sommabile) è legata alla trasformata di f (x) dalla relazione:
Fk (f 0 ) = ik Fk (f )
(4.18)
Infatti integrando per parti si ottiene:

Z
1
0
f 0 (x)eikx dx
Fk (f ) =
2

Z
f (x) ikx
ik
= e
+
f (x)eikx dx

2
2
Z
ik
f (x)eikx dx = ikFk (f ) ,
=
2
dove il termine integrato deve essere nullo affinche la trasformata di Fourier esista.
107
La relazione (4.18) può essere iterata per ottenere le trasformate delle derivate
successive:
Fk (f 00 ) = ikFk (f 0 ) = (ik)2 Fk (f )
Fk [f (n) ] = (ik)n Fk (f ) ,
(4.19)
ovviamente supponendo che la funzione f (x) ammetta derivate fino allennesima

e che f (n) (x) sia sommabile sullasse reale.
Moltiplicando ambo i membri della (4.18) per i e usando la linearità della
Trasformata di Fourier si può simbolicamente stabilire la corrispondenza:
d
k
(4.20)
dx
fra loperatore derivata nello spazio delle funzioni f (x) e la semplice moltiplicazione per k nello spazio delle funzioni F (k); tale corrispondenza è di grande
importanza in Meccanica Quantistica.
i
Dalla disuguaglianza
1
|Fk (f )|
2
|f (x)|dx = cost.
(4.21)
segue che la TF di una funzione sommabile è sempre una funzione limitata. Dalla
(4.19) segue quindi che la TF di una funzione n volte derivabile è almeno O(1/k n )
per k ; in breve, quanto pi`
u una funzione è liscia (ovvero quanto maggiore è
il suo ordine di derivabilità) tanto pi`
u velocemente la sua TF va a zero allinfinito.
Nellesempio 2 della sezione precedente f (x) è discontinua e la sua TF è O(1/k).
Negli esempi 1 e 3 abbiamo invece considerato funzioni infinitamente derivabili,
la cui TF è esponenzialmente soppressa a grandi k.
Se largomento della funzione f (x) viene traslato di una costante reale a, per la F vale la
seguente relazione:
Fk [f (x + a)]
= eika Fk [f (x)].
Infatti
Fk [f (x + a)]
Z
1
f (x + a)eikx dx
=
2
Z
Z
0
0
eika
1
f (x0 )eik(x a) dx0 =
f (x0 )eikx dx0
=
2
2
ika
= e Fk (f ) .
Se si moltiplica la funzione f (x) per un esponenziale, la trasformata di Fourier è:

Fk eix f (x) = Fk+ [f (x)] , R
come si pu`
o facilmente verificare a partire dalle definizioni di F.
108
(4.22)
Se si moltiplica la funzione f (x) per x, le trasformata di Fourier diventa:

d
Fk [f (x)]
(4.23)
dk
come si può facilmente verificare derivando sotto il segno, nellipotesi che la funzione xf (x) sia ancora sommabile sullasse reale. Come la (4.18), anche la (4.23)
si può iterare, ottenendo

n
d
n
Fk [x f (x)] =
i
Fk [f (x)]
(4.24)
dk
Fk [xf (x)] = i
sempre nellipotesi che la funzione xn f (x) sia ancora sommabile sullasse reale.
La (4.23) stabilisce la corrispondenza, duale della (4.20),
xi
d
,
dk
fra moltiplicazione per x nello spazio delle f (x) e la derivata nello spazio delle
F (k).
Lequazione (4.24) mostra che quanto pi`
u rapidamente una funzione decresce
allinfinito, tanto pi`
u la sua TF è liscia (cioè maggiormente derivabile). Nellesempio 1 abbiamo infatti visto che la TF di una funzione O(1/x2 ) allinfinito ha
derivata discontinua nellorigine, mentre la TF di una gaussiana è ancora una
gaussiana (infinitamente derivabile).
Se chiamiamo S lo spazio lineare delle funzioni di prova, rapidamente decrescenti e infinitamente derivabili:
S = {f C ; xn f (x) limitata su R, n N} ,
(4.25)
le (4.19) e (4.24) implicano che la TF manda le funzioni di prova (nella variabile

x) in funzioni di prova (nella variabile k):
Fk (f ) S , f S .
(4.26)
Teorema di convoluzione.
Definiamo la convoluzione g = f1 f2 di due funzioni f1 e f2 :
Z
g(x) =
f1 (x0 )f2 (x x0 )dx0 .
(4.27)
` immediato vedere che il prodotto convolutivo è commutativo e associativo:

E
f1 f2 = f2 f1
f1 (f2 f3 ) = (f1 f2 ) f3 .
109
(4.28)
(4.29)
Il teorema di convoluzione afferma che la trasformata di Fourier della convoluzione di due funzioni è (a parte una costante moltiplicativa) il prodotto delle
trasformate di Fourier delle due funzioni:
Fk (g) = 2Fk (f1 )Fk (f2 ) .

(4.30)
Dimostrazione:
Z
1
dx g(x) eikx
Fk (g) =
2
Z
Z
1
=
dx
dx0 f1 (x0 )f2 (x x0 )eikx
2
Z
Z
1
0
0
dx0 f1 (x0 )eikx f2 (x x0 )eik(xx ) .

dx
=
2
Se ora passiamo dalle variabili (x, x0 ) alle variabili (z = x x0 , x0 ) e scambiamo

lordine di integrazione usando il Teorema di Fubini Tonelli (vedi eq.(D.5) in
Appendice D), otteniamo:
Z
Z
1
0
0 ikx0
Fk (g) =
dx f1 (x )e
dzf2 (z)eikz = 2 Fk (f1 ) Fk (f2 ) .
2
[q.e.d.]
4.1.3
Soluzione di equazioni differenziali mediante la trasformata di Fourier
La proprietà (4.20) trasforma unequazione differenziale a coefficienti costanti in

unelementare equazione algebrica lineare; chiamando U (k) e F (k) le Trasformate di
Fourier dellincognita u(x) e del termine noto f (x), lequazione differenziale
au00 (x) + bu0 (x) + cu(x) = f (x),
(4.31)
ak 2 U (k) + ibk U (k) + c U (k) = F (k),
(4.32)
diventa semplicemente:
da cui è immediato ricavare U (k); infine antitrasformando si ricava la funzione incognita

u(x).
Notare come questo procedimento per risolvere equazioni differenziali a coefficienti
costanti mediante la trasformata di Fourier sia lesatto parallelo di quello illustrato
nel Capitolo 2, dopo leq.(3.37), per luso della Serie trigonometrica di Fourier; allora
richiedevamo che termine noto e soluzione fossero funzioni periodiche, ora abbiamo
lasciato cadere questa richiesta; va tuttavia osservato che la procedura appena descritta
110
per risolvere lequazione (4.31) ha senso solo se il termine noto e la soluzione sono
sommabili e ciò è molto restrittivo.
Torneremo su questo punto quando parleremo della Trasformata di Laplace; per
ora limitiamoci a illustrare un esempio di soluzione di unequazione differenziale (alle
derivate parziali) mediante la trasformata di Fourier.
Esempio: lequazione del calore
Risolviamo lequazione di diffusione del calore
1 T (x, t)
2 T (x, t)
=
x2
t
(4.33)
con la condizione iniziale

T (x, 0) = f (x) .
Fisicamente T (x, t) rappresenta la distribuzione di temperatura al tempo t in una
sbarra di lunghezza infinita, e la distribuzione iniziale è data dalla funzione sommabile
f (x)2 . La costante è la conducibilità termica. Moltiplicando lequazione (4.33) per
eikx e integrando su x da a si ottiene:
Z
Z
2
1 ikx T (x, t)
ikx T (x, t)
e
e
dx =
dx .
x2

t
Chiamando
F (k, t) la trasformata di Fourier rispetto a x della T (x, t), ovvero F (k, t) =
R + ikx
1
e
T (x, t) dx, e ricordando la (4.19) si ottiene
2
k 2 F (k, t) =
1
F (k, t) .
t
Questa è unequazione differenziale del primordine in F (k, t), la cui soluzione è

2
F (k, t) = F (k, 0)ek t .

Dalle condizioni iniziali si ha
Z
Z
1
1
ikx
T (x, 0)e
dx =
f (x)eikx dx F (k)
F (k, 0) =
2
2
da cui
1
2
F (k, t) = ek t
2
f (x)eikx dx .
Essendo la lunghezza della sbarra infinita, lequilibrio termico viene raggiunto alla temperatura
data dal limx f (x). Per avere f (x) sommabile, è necessario che lo zero della scala delle temperature
venga fissato a questo valore.
111
Per ricavare T (x, t) antitrasformiamo secondo Fourier

Z
1
F (k, t)eikx dk
T (x, t) =
2
Z
Z
1
0
2
=
dk
dx0 f (x0 )eikx eikx ek t
2
Z
Z
1
0
k2 t
dx0 f (x0 )eik(x x) ,
dke
=
2
integriamo su k
Z
k2 tik(x0 x)
dke

(x0 x)2
(x0 x)2
k2 t+ik(x0 x) 4t + 4t
dke
(x0 x)2
4t
h
k t+ 2i
= e
dke
r
(x0 x)2
e 4t ,
=
t
0 x
x
t
i2
e giungiamo finalmente al risultato:

r
Z
(x0 x)2
1
T (x, t) =
f (x0 )e 4t dx0 .
4t
2
Si può arrivare pi`

u facilmente allo stesso risultato osservando che ek t è la TF di
2
1 ex /(4t) ; quindi F (k, t) `
e il prodotto di due TF e lantitrasformata è la convolu2t
1
zione di f (x) e 2t
ex
La funzione3
2 /(4t)
.
r
1 (x0 x)2
e 4t (t),
(4.34)
4t
detta nucleo del calore (heat kernel), è la funzione di Green o propagatore delleq. (4.33) ed è tale che la
Z
T (x, t) =
G(x, x0 , t)f (x0 )dx0
(4.35)
G(x, x0 , t) =
Con

(t) =
0
1
t<0
t>0
` necessario introdurla perche tutto il discorso fatto

denotiamo la funzione a gradino di Heaviside. E
k2 t
perde completamente senso per t < 0: e
da gaussiana diventa furiosamente crescente per k
se t < 0.
112
descrive la propagazione del calore dal punto x0 al punto x al tempo t > 0.

Se, per esempio, la sorgente è puntiforme (cioè diversa da zero solo nellorigine):
f (x) = (x)
(per la definizione della delta di Dirac (x), vedi pi`
u avanti, paragrafo 5.3.2) allora il calore si propaga
in modo che la temperatura assume una distribuzione gaussiana di larghezza proporzionale a t:

Z
T (x, t) =
r
0
G(x, x , t)(x )dx = G(x, 0, t) =
1 x2
e 4t .
4t
` anche interessante notare che per t 0+ il nucleo del calore (4.34) tende alla delta di Dirac
E
(vedi eq. (5.63)):
lim G(x, x0 , t) = (x x0 ),
(4.36)
t0+
come deve essere affinche la (4.35) riproduca le condizioni iniziali per t 0+.
113
4.2
Trasformata di Laplace
Come abbiamo già detto nel paragrafo precedente, il metodo della TF permette di risolvere equazioni differenziali a coefficienti costanti soltanto in un campo molto ristretto;
la stessa funzione f (x)=1 non potrebbe essere accettata ne come soluzione ne come
termine noto. Daltra parte, se vogliamo risolvere unequazione con condizioni iniziali
` quindi conveniente
al tempo t = 0, ci interessa sapere solo ciò che succede per t 0. E
considerare una sorta di TF definita da un integrale esteso solo al semiasse delle t > 0.
In tal caso, se f (t) è localmente sommabile, cioè è sommabile su ogni intervallo finito
del semiasse reale t 0, e se esistono 0 R, M > 0 e t0 0 tali che t > t0 valga
0
e t |f (t)| < M
(4.37)
la funzione
g (t) et f (t)(t)
è sommabile sullintero asse reale > 0 . Infatti t > t0

0

(0 )t 0 t
e
f (t) e( )t M
|g (t)| = e
(4.38)
e quindi tende esponenzialmente a zero per t +, rimanendo localmente sommabile

se tale era f (t). Quindi g (t) possiede la trasformata di Fourier:
Z
Z
1
1
t
it
f (t)e (t)e
dt =
f (t)e(+i)t dt .
(4.39)
G () =
2
2 0
Introduciamo ora la variabile complessa s = + i e definiamo F (s) = 2G ().

Allora la (4.39) diventa
Z
f (t)est dt
(4.40)
F (s) =
0
e prende in nome di trasformata di Laplace (TL) della funzione f (t); chiamando

ascissa di convergenza 0 R lestremo inferiore degli 0 per cui vale la (4.37),
la F (s) è definita nel semipiano Res > 0 e ivi è analitica, come si vede facilmente
derivando sotto il segno, poiche, grazie alla (4.38), tn f (t) ha la stessa ascissa di convergenza di f (t), per ogni n naturale. Lantitrasformata di Laplace è definita da un
integrale in campo complesso. Infatti antitrasformando la (4.39) si ricava, nellipotesi
che f (t) sia di classe C 1 nellintorno di t,
Z R
1
t
f (t)e (t) =
lim
G ()eit d ,
(4.41)
2 R R
da cui
1
f (t)(t) =
lim
2 R
R
(+i)t
G ()e
R
114
1
d =
2i
F (s)est ds ,
(4.42)
dove il cammino di integrazione è una retta parallela allasse immaginario del piano di s, di equazione Re s= > 0 . Lantitrasformata (4.42) si indica di solito,
sottintendendo la (t), come
Z +i
1
f (t) =
F (s)est ds ,
(4.43)
2i i
che per essere pi`
u precisi andrebbe scritta come
Z +iR
1
F (s)est ds .
f (t) =
lim
2i R iR
(4.44)
Osserviamo subito: si può dimostrare che F (s) = o(1) per s in ogni direzione
del semipiano di analiticità. Per t < 0 si può quindi applicare il lemma di Jordan
(vedi caso 4) chiudendo il cammino con una semicirconferenza nel semipiano a destra
di Re s = , ottenendo f (t) = 0, visto che F (s) è analitica nel semipiano Re s > 0 .
4.2.1
Esempi
La trasformata di Laplace della funzione f (t) = 1 è

Z
1
est dt = .
F (s) =
s
0
Lintegrale converge per Re(s) > 0 (cioè 0 = 0). Il calcolo dellintegrale (4.43)
mediante il metodo dei residui mostra subito che lantitrasformata di Laplace di
1/s è (t).
La trasformata di Laplace della funzione f (t) = t è

Z
Z
Z
st
d 1
d st
1
st
e dt =
= 2
te dt =
e dt =
F (s) =
s
ds 0
ds s
s
0
0
La trasformata di Laplace della funzione f (t) = tn è
Z
Z n
st
n st
n
t e dt = (1)
F (s) =
e dt
sn
0
0

n Z
n
1
n d
st
n d
= (1)
e dt = (1)
n
n
ds 0
ds
s
n!
= n+1 .
s
La trasformata di Laplace della funzione f (t) = cos t è
Z
Z

1 st it
st
F (s) =
e cos tdt =
e
e + eit dt
2 0
0

(si)t

1
e
e(s+i)t
s

.
=
+
= 2

2 (s i) 0
(s + i) 0
s +1
115
Tutte le funzioni discusse in questi esempi hanno ascissa di convergenza 0 = 0; la loro

Trasformata di Laplace è perciò analitica in tutto il semipiano Res > 0; nellesempio
successivo vedremo che non è sempre cos`.
La trasformata di Laplace della funzione f (t) = eat , con a C, è

Z
e(sa)t
1
st at
e e dt =
,
F (s) =
=

(s a) 0
sa
0
dove è stato necessario supporre Re s > Re a per poter affermare che
limt+ e(sa)t = 0; per Re s < Re a lintegrale che definisce la trasformata di
Laplace diverge, quindi lascissa di convergenza della funzione f (t) = eat è 0 =
Re a; ciò concorda con il fatto che la trasformata di Laplace F (s) è analitica nel
semipiano Re s > Re a.
Per tutti questi esempi lasciamo allo studente la verifica delleq. (4.43). Riflettendo
su questi esempi lo studente si convincerà anche della seguente importante proprietà:
Se la TL Ls (f (t)) ha poli con Re s > 0 la funzione f (t) esplode esponenzialmente
per t +; viceversa se Ls (f (t)) ha singolarità solo a sinistra dellasse immaginario allora f (t) decresce esponenzialmente per t +; se i poli sono sullasse
immaginario f (t) può oscillare o crescere come una potenza di t.
Questa proprietà è di grande importanza per le applicazioni a sistemi fisici e fornisce
un criterio di stabilità nel tempo. Quando, per esempio, un sistema di amplificazione
comincia a produrre un sibilo di ampiezza crescente (fortunatamente limitata dalla non
linearità e quindi saturazione del sistema) possiamo dire che una qualche singolarità
della TL della sua funzione di trasferimento ha acquistato parte reale non negativa.
4.2.2
Propriet`
a della trasformata di Laplace
Studiamo ora alcune proprietà delle trasformate di Laplace, analoghe a quelle viste per
le trasformate di Fourier. Indicheremo la trasformata di Laplace (4.40) con il simbolo
Z
f (t)est dt .
Ls [f (t)] = F (s) =
0
La trasformata di Laplace è lineare (per linearità degli integrali):

Ls [a1 f1 (t) + a2 f2 (t)] = a1 Ls [f1 (t)] + a2 Ls [f2 (t)]
La trasformata di Laplace della derivata f 0 (t), se f 0 (t) esiste e ammette TL, è
legata alla trasformata di f (t) dalla relazione:
Ls [f 0 (t)] = sLs [f (t)] f (0) ,
116
(4.45)
come si ottiene integrando per parti:

Z
0
f 0 (t)est dt
Ls [f (t)] =
0
f (t)est 0
+s
Z
= f (0) + s
f (t)est dt
f (t)est dt = sLs [f (t)] f (0) .
` evidente che qui con f (0) si intende il limite destro di f (x) per x 0.
E
Nellipotesi che anche le derivate successive di f (t) esistano e ammettano TL,
la relazione (4.45) può essere iterata per ottenere le trasformate delle derivate
successive:
Ls [f 00 (t)] =
=
000
Ls [f (t)] =
=
(n)
Ls [f (t)] =
=
sLs [f 0 (t)] f 0 (0)

s2 Ls [f (t)] sf (0) f 0 (0)
(4.46)
00
00
sLs [f (t)] f (0)
s3 Ls [f (t)] s2 f (0) sf 0 (0) f 00 (0)
(4.47)
(n1)
(n1)
sLs [f
(t)] f
(0)
n
n1
s Ls [f (t)] s f (0) sn2 f 0 (0) ... sf (n2) (0) f (n1) (0) .
(4.48)
Dalla (4.48) segue che se la funzione f (t) è n volte derivabile e f (n) (t) ammette
TL vale
f (0) f 0 (0)
f (n1) (0) Ls [f (n) ]
Ls [f (t)] =
+ 2 +
+
,
(4.49)
s
s
sn
sn
che dà utili informazioni sullandamento per s della trasformata di Laplace
(teorema di Tauber) a partire dalla funzione f e dalle sue derivate in t = 0+.
La (4.49) si può verificare anche negli esempi precedenti; la TL di tn è O(sn1 )
` da notare
per grande s perchè vanno a zero tutte le prime n derivate in t = 0. E
che, anche se f (t) C , non è affatto detto che la serie che facilmente si deduce
dalla (4.49) converga. In realtà essa è in generale una serie asintotica, che sarà
definita nel corso di Metodi Matematici della Fisica II.
La trasformata di Laplace dellintegrale di una funzione g(t) è legata alla trasformata di g(t)
dalla relazione:
Z x

1
Ls [g(x)] .
(4.50)
Ls
g(t)dt
=
s
0
Partendo dalla formula per la trasformata di Laplace per le derivate (4.45) e ponendo g(x) =
f 0 (x) si ottiene infatti
Z x
f (x) = f (0) +
g(t)dt .
0
117
La (4.45) diventa cos`:

Z
Ls [g(x)] = sLs f (0) +

Z
g(t)dt f (0) = sLs [f (0)] + sLs

g(t)dt f (0) .
Ricordando ora che Ls [1] = 1/s, otteniamo

1
Ls [g(x)] = sf (0) + sLs
s
Z

g(t)dt f (0) ,
ovvero
x
Z
Ls [g(x)] = sLs

g(t)dt ,
da cui segue immediatamente la (4.50).

Se largomento della funzione f (t) viene traslato di una costante a, per la trasformata di Laplace
vale la seguente relazione:

Z a
as
st
Ls [f (t + a)] = e
Ls [f (t)] (a)
f (t)e dt .
(4.51)
0
Per dimostrare la (4.51) bisogna distinguere i due casi a < 0 e a > 0. Ricordiamo infatti che la
trasformata di Laplace è un integrale tra 0 e e che è sottintesa una (t), che implica f (t) = 0
se t < 0. Quindi
Z
Z
0
(t0 =t+a)
Ls [f (t + a)] =
f (t + a)est dt =
f (t0 )es(t a) dt0
0
a
Z
sa
0 st0 0
= e
f (t )e
dt .
a
Ora, se a < 0,
Ls [f (t + a)] = e
sa
f (t )e
st0
dt = e
sa
f (t0 )est dt0 = esa Ls [f (t)] .
Se invece a > 0,
Ls [f (t + a)]
= e
sa
f (t )e
st0
sa
dt = e
f (t )e
a
0
Z a
0
= esa Ls [f (t)] esa
f (t0 )est dt0 .
st0
dt e
sa
f (t0 )est dt0
Se si moltiplica la funzione f (t) per un esponenziale, la trasformata è:

Ls et f (t) = Ls [f (t)] , C
come si pu`
o facilmente verificare derivando sotto il segno la L. Per esempio ricordando che
Ls [sin t] =
1
s2 + 1
segue che
Ls [e2t sin t] =
senza un calcolo esplicito della TL.
118
1
(s 2)2 + 1
Se si moltiplica la funzione f (t) per t, la trasformata diventa:

Ls [tf (t)] =
d
Ls [f (t)] ,
ds
come si può facilmente verificare a partire dalle definizioni di L. Notare che,

diversamente da quanto avviene per la TF, qui siamo sempre sicuri che se f (t)
ammette TL anche tn f (t) la ammette, n N; equivalentemente, ogni TL è
sempre infinitamente derivabile nella sua regione di convergenza, mentre ciò non
1
è affatto detto per la TF (vedi per esempio la TF di a2 +x
2 , par.4.1.1).
Teorema di convoluzione. Sia g(t) la convoluzione di due funzioni f1 e f2 nulle
per t < 0 (la loro TL non dipende da questo). Allora t > 0 abbiamo
Z
Z t
0
0
0
g(t) =
f1 (t )f2 (t t )dt =
f1 (t0 )f2 (t t0 )dt0 ,
(4.52)
dove abbiamo usato lannullarsi di f1,2 per argomenti negativi. Allora

Ls [g(t)] = Ls [f1 f2 ] = Ls [f1 (t)] Ls [f2 (t)] .
(4.53)
La dimostrazione è del tutto analoga a quella vista per le trasformate di Fourier.
4.2.3
Trasformate di Laplace ed equazioni differenziali lineari

a coefficienti costanti
Come si è detto nel par. 4.1.3 il metodo della trasformata di Fourier per risolvere
equazioni differenziali lineari si può applicare solo in un numero ristretto di casi, cioè
quando il termine noto e la soluzione dellequazione differenziale sono sommabili. La
trasformata di Laplace permette di risolvere equazioni differenziali lineari a coefficienti
costanti
du
d2 u
+ c0 u = f (t)
(4.54)
c2 2 + c1
dt
dt
per una classe pi`
u estesa di funzioni e permette inoltre di tener conto automaticamente
delle condizioni iniziali
u(0) = u0

du
= u1 .
dt
(4.55)
t=0
Infatti trasformando secondo Laplace la (4.54) e ponendo Ls [u(t)] U (s) e Ls [f (t)]

F (s) (ammesso che queste esistano) si ottiene, utilizzando le (4.45) e (4.46),

c2 s2 U (s) su0 u1 + c1 [sU (s) u0 ] + c0 U (s) = F (s) ,
(4.56)
119
che dà immediatamente

U (s) =
F (s)
c2 su0 + c2 u1 + c1 u0
+
.
2
2
c2 s + c1 s + c0
c2 s + c1 s + c0
(4.57)
Antitrasformando si ottiene (t)u(t). Lantitrasformata del primo addendo dà la soluzione generale dellomogenea associata (interpretando u0 e u1 come parametri liberi),
mentre lantitrasformata del secondo dà la soluzione particolare dellinomogenea con
condizioni iniziali u(0) = u0 (0) = 0. Notare che la soluzione cos` ottenuta è particolarmente interessante quando la sollecitazione esterna f (t) sia inserita al tempo
t = 0. In questo caso (t)u(t) ci dice cosa succede dopo aver chiuso o aperto linterruttore. Ovviamente la procedura è estendibile a equazioni differenziali lineari di ordine
qualsiasi.
Esempio
Consideriamo il circuito oscillante di Fig. 3.1 e cerchiamo la soluzione dellequazione
differenziale
d2 u
du
+ LC 2 = f (t)
(4.58)
u + RC
dt
dt
con le condizioni iniziali
u0 =
i0 =
lim u(t)
(4.59)
t0+
lim C
t0+
du
.
dt
(4.60)
Trasformando la (4.58) secondo Laplace e chiamando U (s) e F (s) le trasformate di

u(t) e f (t) otteniamo

i0
2
U (s) + RC [sU (s) u0 ] + LC s U (s) su0
= F (s) .
(4.61)
C
Consideriamo due casi:
a) Chiusura del circuito: se il circuito è inizialmente aperto e il condensatore è scarico e lo si chiude allistante t = 0 su un generatore che fornisca una tensione
alternata V eit (in particolare continua, se = 0), il termine noto è:
f (t) = (t)V eit ,
(4.62)
u0 = i0 = 0
(4.63)
le condizioni iniziali sono

e la trasformata di Laplace di f (t) è
Z
F (s) =
V eit est dt =
0
120
V
.
s i
(4.64)
Quindi la (4.61) fornisce in questo caso:

U (s) =
V
.
(s i)(1 + RCs + LCs2 )
(4.65)
Per ottenere u(t) dobbiamo antitrasformare la (4.65), che possiede tre poli semplici in
s0 = i ,
s1,2 =
R
1
2L 2
con =
R2
4
.
2
L
LC
(4.66)
Ora, se > 0, s1 e s2 giacciono

sullasse reale negativo (perche R > 0), se
R
i
< 0, s1 = s2 = 2L + 2 e se = 0 i due poli sono reali e coincidenti

(polo doppio). Lantitrasformata di U vale, secondo la (4.43),
Z r+i
1
u(t) =
est U (s) , con r > 0
(4.67)
2i ri
ovvero, se s1 6= s2 ,

V X
est
u(t) =
Res
LC i=0,1,2
(s i)(s s1 )(s s2 ) s=si

eit
es1 t
e s2 t
V
+
+
.
=
LC (s1 i)(s2 i) (s1 i)(s1 s2 ) (s2 i)(s2 s1 )
(4.68)
A parte il termine sinusoidale
V eit /LC
V eit
V
eit
=
=
, (4.69)
LC (s1 i)(s2 i)
s1 s2 i(s1 + s2 ) 2
1 + iRC 2 LC
che riproduce esattamente la (3.18), la tensione u(t) è quindi una somma di
funzioni sinusoidali smorzate nel caso < 0 , mentre per > 0 è una somma di
` immediato verificare che a t = 0 le condizioni iniziali
esponenziali decrescenti. E
sono soddisfatte:

1
1
1
V
+
+
=0
u(0) =
LC (s1 i)(s2 i) (s1 i)(s1 s2 ) (s2 i)(s2 s1 )

V
ieit
s1 es1 t
s2 es2 t
i(t) =
+
+
L (s1 i)(s2 i) (s1 i)(s1 s2 ) (s2 i)(s2 s1 )
i(0) = 0 ,
mentre per
una volta che il condensatore si sia caricato, quindi
t
+
(ovvero

1
1
R
per t Res1 = Res2 = 2L ) la tensione ai capi del condensatore si riduce a
121
quella sinusoidale (4.69). Se invece s1 = s2 , U (s) ha un polo doppio in s1 e il

residuo vale

d est
test
est
est
es1 t
=
=
(s1 t 1) , (4.70)
Res
=
s(s s1 )2 s=s1
ds s s=s1
s
s2 s=s1
s21
dove si è scelto per semplicità = 0, da cui la soluzione
u(t) =

V 1
s1 t
1
+
e
(ts
1)
,
1
LC s21
che verifica le condizioni iniziali

V 1
[1 1] = 0
LC s21
V 1 s1 t 2
i(t) =
e (ts1 s1 + s1 ) i(0) = 0
LC s21
u(0) =
e per t + si comporta come nel caso precedente.

b) Apertura del circuito: se allistante t = 0 il circuito viene aperto, dopo essere stato per lungo
tempo a contatto con la batteria a tensione costante V , si avrà
f (t) = 0 per t > 0,
con le condizioni iniziali
u0 = V e i0 = 0 .
Leq. (4.61) diventa in questo caso (F (s) = 0):
U (s)(1 + RCs + LCs2 ) = V (RC + sLC)
da cui
U (s) = V
s+ R
L
s2 + R
Ls +
1
LC
=V
s+ R
L
.
(s s1 )(s s2 )
Si noti che, per qualunque valore di , s1 , s2 6= R/L; quindi non cè cancellazione fra numeratore e denominatore e U (s) ha sempre due poli semplici (o un polo doppio). Per s1 6= s2 si
ha
(
)
X
est s + R
L
u(t) = V
Res
(s s1 )(s s2 )
s=s ,s
1
es1 t (s1 + R/L) es2 t (s2 + R/L)

= V
,
s1 s2
che soddisfa le condizioni iniziali:
u(0)
i(t)
s1 s2
=V
s1 s2

CV
R
R
s1 t
s2 t
e s1 s1 +
e s2 s2 +
i(0) = 0 .
s1 s2
L
L
= V
=
122
R
Se invece s1 = s2 = 2L
(
u(t)
Res

= V
est s + R
L
(s s1 )2
)
s=s
d
=V
ds
R
s+
L
1

R
1 + s1 t + t es1 t .
L
st
s=s1
Si verifica facilmente che le condizioni iniziali sono soddisfatte:

u(0)
i(t)

d
R
1 + s1 t + t es1 t
dt
L

R s1 t
R
e
= CV s21 t + s1 t + 2s1 +
L
L

R
i(0) = 2s1 +
=0.
L
=
CV
In entrambi i casi, per t + sia la tensione u(t) che la corrente i(t) tendono esponenzialmente
a zero, come ci si aspetta. Anche qui t + significa t 2L/R.
123
Capitolo 5
Spazi L2 e distribuzioni
Sin dallinizio del Capitolo 3.2 abbiamo visto limportante ruolo giocato dalle relazioni
di ortogonalità (3.36) nel calcolo dei coefficienti di Fourier. In questo capitolo torneremo
su questo argomento trattandolo in modo un po pi`
u approfondito, e ciò ci permetterà
di dare una definizione di convergenza delle serie (e trasformata) di Fourier molto pi`
u
generale ed appropriata. Allargheremo inoltre il campo delle funzioni introducendo il
concetto di distribuzione; riusciremo cos` a derivare funzioni anche nei loro punti di
discontinuità di I specie e ad ampliare di molto il campo delle funzioni che ammettono
TF
5.1
Spazi L2 e serie di Fourier
Linsieme di tutte le funzioni a valori complessi in un intervallo reale [a, b] costituisce

uno spazio vettoriale, ovvero spazio lineare, sui complessi. Ogni funzione rappresenta
quindi un vettore e lo spazio vettoriale si chiama spazio funzionale; la somma di vettori
e il prodotto per un numero complesso sono definiti dalle corrispondenti operazioni
sulle funzioni. Nello spazio funzionale si può inoltre definire il prodotto scalare di
due vettori, f e g, come
Z
(f , g) =
f (x)g(x)dx .
(5.1)
Per aiutare la memoria, è utile notare lanalogia fra lintegrale (5.1) e la consueta
definizione
(f , g) =
fi gi .
(5.2)
di prodotto scalare negli spazi vettoriali (sui complessi) di dimensione finita; negli spazi
funzionali la variabile x gioca un ruolo analogo a quello dellindice i che individua la
124
componente i-esima; la somma su i è sostituita dallintegrale perche lindice x corre su

un insieme continuo, lintervallo [a, b].
Lintegrale (5.1) soddisfa le proprietà del prodotto scalare negli spazi unitari (cioè
negli spazi vettoriali su C dotati appunto di prodotto scalare) :
1) linearit`
a nel secondo fattore:
Z
f (x) [g(x) + h(x)] dx

a
Z b
Z b
=
f (x)g(x)dx +
f (x)h(x)dx
(f , g + h) =
= (f , g) + (f , h) .
2) Hermiticit`
a:
Z
(f , g) =
f (x)g(x)dx
Z
=
f (x)g (x)dx = (g, f ) ,
che si riduce alla commutatività nel caso di spazio vettoriale sui reali.
3) Positivit`
a1 :
Z
(f , f ) =
|f (x)|2 dx 0 ,
(5.3)
La radice quadrata del numero non negativo (f , f ) si dice norma del vettore f :
p
(5.4)
kf k = (f , f )
Dalle proprietà 1) e 2) segue che il prodotto scalare (5.1) è antilineare nel primo
fattore:
(f + g, h) = (f , h) + (g, h) .
Vale inoltre per il prodotto scalare la seguente disuguaglianza, detta disuguaglianza di
Schwartz:
|(f , g)| kf k kgk .
(5.5)
La disuguaglianza di Schwartz si può considerare come unestensione agli spazi funzionali della ben nota disuguaglianza della geometria euclidea |~a ~b| = |ab cos| ab,
dove a e b sono le lunghezze (norme) dei due vettori e langolo compreso.
1
Questa propriet`
a rende chiaro perche nelle definizioni (5.1) e (5.2) le componenti del primo
vettore siano soggette alla complessa coniugazione.
125
` importante osservare che la definizione (5.1) di prodotto scalare ha senso solo se

E
le funzioni che consideriamo sono quadrato sommabili cioè se esiste lintegrale (5.3)
che definisce la norma di un vettore2 .
Il prodotto scalare deve ancora soddisfare una proprietà: leguaglianza nella (5.3)
deve valere se e solo se f (x) = 0. Ma anche una funzione nulla ovunque tranne che in
in alcuni punti isolati ha norma nulla. Occorre essere precisi al riguardo, e dobbiamo
anche dire che nella definizione di prodotto scalare, e quindi di norma, si deve usare
una generalizzazione dellintegrale di Riemann dovuta a Lebesgue. Per i nostri scopi
non è necessario definire lintegrale di Lebesgue, ma è sufficiente enunciarne alcune
proprietà.
Propriet`
a A: ogni funzione assolutamente integrabile alla Riemann, in modo
proprio o improprio, lo è anche alla Lebesgue e i due integrali coincidono.
Viceversa, se in un punto x0 [a, b] una funzione diverge troppo per essere integrabile alla Riemann (per esempio una f (x) = O(x ) per x 0, con 1, su un
intervallo dintegrazione comprendente lorigine) allora non è nemmeno integrabile alla
Lebesgue. Analogamente, se allinfinito una funzione va a zero troppo lentamente per
essere integrabile alla Riemann su intervallo infinito (per esempio una f (x) = O(x )
per x , con 1), allora non è nemmeno integrabile alla Lebesgue.
Propriet`
a B: come per lintegrale di Riemann, ogni funzione assolutamente integrabile è anche integrabile, ma per lintegrale di Lebesgue vale anche il viceversa3 .
Grazie alla proprietà A, per calcolare un integrale alla Lebesgue di fatto si continua a
calcolare il solito integrale di Riemann. Perche allora complicarci la vita con lintegrale
di Lebesgue? La ragione è che esistono delle funzioni piuttosto bizzarre che sono
sommabili, cioè integrabili alla Lebesgue, senza esserlo alla Riemann 4 ; tali funzioni,
del tutto prive di interesse per la fisica, sono però essenziali per rendere completo lo
spazio funzionale, nel senso che preciseremo fra poco.
Un esempio di tali funzioni è la funzione di Dirichlet, cos` definita su tutto lasse
reale:

1
xQ
d(x) =
0
x
/Q,
dove Q è linsieme dei numeri razionali. Si può dimostrare che tale funzione, che
evidentemente non è integrabile alla Riemann, è però sommabile e che il suo integrale
di Lebesgue vale zero.
2`
E molto facile mostrare che lesistenza della norma di due vettori implica che anche il loro prodotto
scalare esiste.
3
Lassoluta integrabilit`
a implica lintegrabilità alla Lebesgue solo nellipotesi che la funzione sia
misurabile, ma questo è il caso per tutte le funzioni di interesse in fisica; quindi daremo sempre per
scontata la misurabilit`
a, senza nemmeno preoccuparci di definirla.
4
Inoltre lintegrale di Lebesgue gode di proprietà, discusse in Appendice D, che rendono molto pi`
u
semplice derivare sotto il segno, scambiare il limite con lintegrale e scambiare lordine dintegrazione
in integrali multipli.
126
La funzione di Dirichlet è diversa da zero solo sui razionali, cioè su uninfinità

numerabile di punti; il fatto che il suo integrale si annulli è un caso particolare della
pi`
u generale
Propriet`
a C: Condizione necessaria e sufficiente affinche una funzione a valori
reali non negativi abbia integrale di Lebesgue nullo è che essa sia quasi ovunque nulla
nellintervallo (finito o infinito) di integrazione; dicendo che una proprietà vale quasi
ovunque (che si abbrevia con q.o.) su un intervallo si intende che linsieme dei punti
in cui non vale sia di misura nulla, cioè, in pratica, finito o infinito numerabile5 .
Come conseguenza della Proprietà C si può allora affermare che il vettore nullo
(cioè il vettore di norma nulla, che vogliamo sia unico) è rappresentato dallintera
classe delle funzioni quasi ovunque nulle; analogamente, a ogni vettore dello spazio
astratto corrisponde una classe di funzioni quasi ovunque uguali quadrato sommabili;
tale spazio, dotato del prodotto scalare (5.1), 6 si denota con il simbolo L2 (a, b). Si
denota invece con L(a, b) lo spazio delle funzioni sommabili 7
Se e solo se lintervallo (a, b) è finito vale L2 (a, b) L(a, b), ovvero ogni funzione
quadrato sommabile, su un intervallo finito, è ivi anche sommabile; per convincersene
basta lidentità
Z
b
f (x) dx = (1, f ) ,
a
dove con 1 intendiamo il vettore corrispondente alla funzione f (x) = 1, che è sommabile
su ogni intervallo finito. Non è invece vero che ogni funzione sommabile sia quadrato
sommabile; per esempio
1
/ L2 (0, 1) .
f (x) = L(0, 1) ma
x
Su intervallo infinito non esiste invece alcuna relazione di inclusione fra L e L2 ; per
esempio
1
f (x) =
L2 (0, ) ma
/ L(0, ) .
2
1+x
Due funzioni f (x) e g(x) appartenenti a L2 (a, b) si definiscono ortogonali se il loro
prodotto scalare è nullo:
Z b
(f , g) =
f (x)g(x)dx = 0 .
a
5
Nel piano, un insieme di misura nulla può essere costituito non solo da uninfinità numerabile di
punti, ma anche da un numero finito o infinito numerabile di segmenti (o curve differenziabili) e cos`
via.
6
che evidentemente non dipende da quale funzione si sceglie per rappresentare la classe.
7
In L(a, b) non è definito il prodotto scalare, anche se è ancora uno spazio vettoriale normato.
127
Un sistema di funzioni {i (x)} L2 (a, b) si definisce ortonormale (ON) se

Z b
i (x)k (x)dx = ik , i, k .
(i , k ) =
a
Sia {1 2 ..., n } un sistema ortonormale finito di L2 (a, b) e sia f (x) L2 (a, b)

definita da:
f (x) =
n
X
ck k (x)
(5.6)
k=1
o equivalentemente, in linguaggio vettoriale:

f=
n
X
ck k ,
(5.7)
k=1
dove ck sono arbitrari numeri complessi. Le componenti ci del vettore f lungo le

direzioni i si calcolano subito moltiplicando scalarmente la (5.7) per il vettore i ; si
ottiene:
Z b
n
n
X
X
i (x)f (x)dx =
(i , f ) =
ck (i , k ) =
ck ik = ci .
(5.8)
a
k=1
k=1
Questo risultato era ovviamente atteso, poiche stiamo lavorando in uno spazio a numero finito di dimensioni, la varietà n-dimensionale generata dai vettori 1 , n . In
L2 (a, b) esistono però sistemi ON formati da uninfinità numerabile di vettori; per
esempio la eq.(3.36) mostra che in L2 (0, 2) le funzioni trigonometriche (in forma
esponenziale)
eilx
(5.9)
l (x) =
2
formano un sistema ON infinito numerabile. Per estendere la (5.6) al caso in cui n
sia infinito, bisogna affrontare il problema della convergenza della serie, e questo è
largomento centrale di questo paragrafo.
Se {i , i = 1, 2, ...} è un insieme infinito numerabile di funzioni ortonormali di
2
L (a, b) e f (x) L2 (a, b), consideriamo dapprima la proiezione ortogonale fN del
vettore f sul sottospazio generato dai vettori ON 1 , N . Essa è data da
fN (x) =
N
X
N = 1, 2 ,
ci i (x)
(5.10)
i=1
con
Z
ci = (i , f )
i (x)f (x)dx,
128
(5.11)
e si dice ridotta N -esima dello sviluppo in serie di Fourier del vettore f nel sistema
{i }, mentre le costanti ci sono i coefficienti di Fourier dello sviluppo.
Dimostriamo innanzitutto la seguente disuguaglianza, detta disuguaglianza di
Bessel:
X
|ck |2 (f , f ) .
(5.12)
k=1
Da
kf fN k2 0
(5.13)
segue che
kf fN k2 = (f , f ) +
N
X
ck cl (k , l )
k,l=1
= (f , f ) +
N
X
= (f , f )
[ck (f , k ) + ck (k , f )]
k=1
|ck |2
k=1
N
X
N
X
N
X
(ck ck + ck ck )
k=1
|ck |2 0 ,
k=1
da cui
N
X
|ck |2 (f , f ) .
k=1
Questo risultato vale per qualunque N . Facendo tendere N a infinito si ottiene infine la
disuguaglianza di Bessel (5.12): la serie dei moduli quadrati dei coefficienti di Fourier
di un vettore f è minore o uguale alla norma quadrata del vettore.
Se il sistema di vettori {k } è tale che la disuguaglianza di Bessel vale con il segno
= per ogni f , allora il sistema di vettori {k } si dice completo e costituisce una base
ON in L2 (a, b). La (5.12) diventa allora
|ck |2 = (f , f )
(5.14)
k=1
e prende il nome di equazione di Parseval.

In uno spazio unitario a N dimensioni un sistema ON è completo se e solo se è
costituito da N vettori; un insieme infinito numerabile invece non cambia il numero
dei suoi elementi (la sua potenza) anche se da esso ne tolgo uno o pi`
u; quindi lunico
criterio che mi assicura che non mi sia perso dei vettori base, cioè che il sistema sia
completo, è, come abbiamo appena detto, che sia soddisfatta lequazione di Parseval
(5.14), che generalizza il Teorema di Pitagora agli spazi infinito dimensionali.
129
Uno spazio unitario infinito dimensionale in cui esista un sistema ON completo

infinito numerabile si dice separabile; si può dimostrare che L2 (a, b) è separabile,
qualunque sia lintervallo (a, b), finito o infinito; in particolare per (a, b) = (0, 2) un
sistema completo è dato dalle funzioni trigonometriche (5.9). Per aiutare la memoria
si può dire, in modo molto rozzo, che uno spazio separabile ha dimensione infinita
numerabile.
Lequazione di Parseval ci permette anche di dare un senso preciso alla convergenza
della serie di Fourier. Infatti da:
kf fN k2 = (f , f )
N
X
|ck |2
(5.15)
k=1
e dallequazione di Parseval (5.14) segue che:

lim kf fN k2 = 0.
(5.16)
Per definizione di limite nello spazio vettoriale astratto , lequazione (5.16) è equivalente alla scrittura:
f = lim fN =
N
ck k .
(5.17)
k=1
Abbiamo cos` definito la somma della serie di Fourier nello spazio vettoriale astratto;
in termini di funzioni la (5.16) significa:
Z b
lim
|f (x) fN (x)|2 dx = 0,
(5.18)
N
dove la ridotta N -esima della serie di Fourier è data dalla (5.10) con i coefficienti (5.11).
La (5.18) si abbrevia nel modo seguente:
f (x) = l . i. m. fN (x) ,
N
(5.19)
dove l.i.m. si legge per limite in media quadratica. La convergenza della serie
di Fourier si può quindi esprimere nei 5 modi (5.14), (5.16), (5.17), (5.18), (5.19)
perfettamente equivalenti.
Molto spesso si scrive pi`
u semplicemente (e lo faremo anche noi):
f (x) =
ci i (x),
(5.20)
i=1
ma lespressione (5.20) è corretta solo se interpretata come unabbreviazione della

(5.19); in generale, la serie dei numeri complessi a secondo membro della (5.20) può
130
non convergere, o convergere a un valore diverso dal numero complesso f (x); infatti le
eq. (5.18) o (5.19) non ci dicono assolutamente nulla della convergenza puntuale della
serie a secondo membro della eq.(5.20), come era da aspettarsi, visto che la funzione
f (x) è solo un rappresentante allinterno di una classe di funzioni quasi ovunque uguali;
quindi, per esempio, se in un punto x0 si aggiunge a f (x0 ) una costante, non cambia
nulla nei coefficienti di Fourier ci dati dalla (5.11) e quindi nel secondo membro della
eq.(5.20), ma ovviamente cambia il primo membro.
La convergenza (5.19) richiede solo che f L2 (a, b) e che il sistema {i } sia ONC;
la convergenza puntuale della (5.20) in uno specifico punto x0 [a, b] richiede invece
tuttaltre condizioni su f , discusse nel Cap. 3.2.
Si può dimostrare che lo spazio L2 (a, P
b) è completo 8 , cioè che per ogni infinità
numerabile di numeri complessi ci tali che i |ci |2 < , esiste una f L2 (a, b) tale che
P
` per questa dimostrazione che è importante che lintegrale sia di Lebesgue
f = ci i (E
e che vengano quindi incluse nello spazio funzionale anche funzioni non integrabili alla
Riemann).
Uno spazio vettoriale sui complessi, dotato di prodotto scalare e completo, si dice
spazio di Hilbert.
La completezza di uno spazio, in particolare dello spazio L2 (a, b), e la sua separabilit`
a, permettono di stabilire una corrispondenza biunivoca fra i suoi elementi f e gli
insiemi infiniti numerabili dei loro coefficienti di Fourier ci , che possiamo considerare
come le componenti di un vettore rispetto alla base i prefissata.
`
P che tale corrispondenza conserva i prodotti scalari: se f (x) =
P E facile mostrare
i di i , allora
i ci i e g(x) =
!
Z b
X
X
X
(f , g) =
f (x)g(x) dx =
ci i , g =
ci (i , g) =
ci di ,
(5.21)
a
i=1
i=1
i=1
dove si sono usate lantilinearità e la continuità9 del prodotto scalare nel primo fattore.
Quindi tutti gli spazi di Hilbert separabili (infinito dimensionali), in particolare gli
spazi L2 (a, b), sono isomorfi fra loro e con lo spazio di Hilbert delle componenti, i
cui
sono i vettori colonna di componenti ci , con i = 1, 2, .., con il vincolo
Pelementi
2
|c
|
<
; le proprietà di uno spazio si traducono in corrispondenti proprietà
i=1 i
dellaltro.
5.1.1
Trasformata di Fourier in L2
Gli spazi L2 intervengono spesso in fisica; in particolare lambiente in cui vive naturalmente la Meccanica Quantistica è quello degli spazi di Hilbert separabili, che
8
Qui la parola completo viene usata in due contesti diversi: la completezza dello spazio non ha
nulla a che fare con la completezza di un sistema di funzioni o vettori.
9
La continuit`
a del prodotto scalare, cioè luguaglianza (limN fN , g) = limN (fN , g), può
essere facilmente dimostrata a partire dalla disuguaglianza di Schwartz (5.5) e dalla definizione (5.16)
di f = limN fN .
131
` quindi concettualpossono essere realizzati da spazi di funzioni quadrato sommabili. E

mente molto importante estendere la definizione di Trasformata di Fourier, che nella
sua forma (4.7) si applica solo a f (x) L(R), anche a f (x) L2 (R).
Su intervallo infinito una funzione può essere quadrato sommabile senza essere
sommabile (vedi esempio 5.6 in Appendice D), ma su intervallo finito L2 (a, b) L(a, b);
quindi per una generica f (x) L2 (R) lintegrale
Z N
1
eikx f (x)dx = FN (k)

(5.22)
2 N
esiste certamente, ma non è detto che ne esista il limite (puntuale) per N . Si
può invece dimostrare che f (x) L2 (R) esiste sempre una F (k) L2 (R) tale che
F (k) = l . i . m. FN (k) ,
(5.23)
dove l.i.m. è stato definito sopra la (5.19), ovvero

Z +
2
|F (k) FN (k)|2 dk = 0 .
lim kF FN k = lim
N
Si può anche dimostrare che la f (x) è la antitrasformata di Fourier della F (k) nello
stesso senso di media quadratica, ovvero
Z +
|f (x) fN (x)|2 dx = 0 ,
(5.24)
f (x) = l . i . m. fN (x) lim
N
dove
1
fN (x) =
2
In fisica si continua a scrivere
Z +
1
F (k) =
eikx f (x)dx ,
2
eikx F (k) dk .
(5.25)
1
f (x) =
2
eikx F (k)dk
(5.26)
anche per f (x) L2 , ma la scrittura corretta è data dalle equazioni (5.22), (5.23),
(5.24),T (5.25). Naturalmente le (5.26) sono anche formalmente corrette se f (x)
L(R) L2 (R), F (k) L(R) e f (x) è di classe C 1 nel punto x. Come lo sviluppo
in serie di Fourier, anche la Trasformata di Fourier conserva i prodotti scalari (ha
senso parlare di questi perche lavoriamo in L2 ); vale cioè luguaglianza di Parseval
generalizzata
Z
Z
+
f (x)g(x) dx ,
F (k)G(k) dk =
dove f, g L2 (R) e F, G L2 (R) sono le loro trasformate di Fourier.

132
(5.27)
La dimostrazione è immediata10 , usando il teorema di Fubini Tonelli (vedi Appendice D) per scambiare lordine di integrazione:
Z +
Z
Z +
1
F (k)G(k) dk =
dk F (k)
dx eikx g(x)
2
Z +
Z +
Z +
1
=
dx g(x)
dk eikx F (k) =
f (x)g(x) dx .
2
Il contenuto di questo paragrafo va sotto il nome di Teorema di Plancherel

per i matematici. Caso particolare della (5.27) è la conservazione della norma; un
interessante esempio a questo proposito è il seguente esempio.
Esempio
La trasformata di Fourier della funzione
f (t) = et/T sin(0 t)(t) ,
che descrive il moto di un oscillatore armonico smorzato, è
Z
Z
1
1
it
f (t)e
dt =
et/T sin(0 t)eit dt
F () =
2
2 0

Z
Z
1
(1/T ii0 )t
(1/T i+i0 )t
=
e
dt
e
dt
2i 2 0
0

1
1
1
,
2i 2 1/T + i( 0 ) 1/T + i( + 0 )
cioè

1
1
1
F () =
.
( 0 ) i/T
2 2 ( + 0 ) i/T
Per dare uninterpretazione fisica alle funzioni f (t) e F (), supponiamo che f (t)
sia il campo elettrico di unonda irradiata. RAllora la potenza irradiata è W |f (t)|2 e
lenergia totale irradiata è proporzionale a 0 |f (t)|2 dt. Dal teorema di Parseval

Z
Z
2
|f (t)| dt =
|F ()|2 d .
0
Quindi |F ()|2 rappresenta (a meno di costanti) lenergia irradiata per intervallo unitario di frequenza:
|F ()|2 =
1
02
2
2 ( 2 2 )2 + 2 02 +
+
0
T2
10
1
T4
Per semplicit`
a qui si usano le ipotesi (non necessarie) che g(x) sia anche sommabile e valga
puntualmente la (4.6).
133
Se T è molto grande (T 01 ), lenergia irradiata per intervallo unitario di frequenza

è fortemente piccata in = 0 e la larghezza del picco è inversamente proporzionale a
T:
02 T 2
.
|F ()| ' 2
(0 2 )2 T 2 + 2(02 + 2 )
2
5.2
Alcuni sistemi ONC
Passiamo ora a considerare alcuni esempi di spazi L2 e di relativi sistemi ONC al loro
interno.
Esempio: le funzioni trigonometriche
Il sistema delle funzioni esponenziali
ikx
e
,
2
k = 0, 1, 2, ...
(5.28)
che avevamo già visto in (5.9), costituisce una base ON in L2 (0, 2). Usando le formule
di Eulero, esso può essere riscritto come sistema trigonometrico:

1 cos nx sin nx
, , ; n = 1, 2,
(5.29)
2
Il sistema delle potenze e i polinomi ortogonali
Oltre al sistema delle funzioni trigonometriche, di cui ci siamo occupati finora, un altro
esempio di sistema completo in uno spazio L2 è il sistema delle potenze
{xn } = {1, x, x2 , x3 , ...}
(5.30)
Si può dimostrare che il sistema (5.30) è completo in L2 in qualsiasi intervallo finito

(a, b), ma evidentemente esso non è ortogonale:
Z
xn xm dx 6= 0
(5.31)
anche per n 6= m.
` possibile tuttavia ortogonalizzarlo, passando dalle funzioni xn a loro combinazioni
E
lineari, che saranno polinomi Pn (x), costruiti in modo tale che11
Z b
Pn (x)Pm (x)dx = hn nm .
a
11
Per ragioni di comodit`

a e tradizione non sempre le costanti positive hn sono scelte uguali a 1.
134
A questo scopo si sceglie P0 (x) = 1, si pone P1 (x) = x + e si fissa in modo che

P1 (x) sia ortogonale a P0 (x); poi si pone P2 (x) = x2 + x + e si fissano e in modo
` ovvio che i coefficienti , ,
che P2 (x) sia ortogonale a P1 (x) e a P0 (x), e cos` via. E
2
... dipenderanno dallintervallo (a, b) che definisce L (a, b); inoltre i polinomi Pn (x)
cos` ottenuti potranno essere moltiplicati per opportune costanti per normalizzarli a 1
o a unaltra costante hn a scelta.
Il sistema di polinomi ortogonali cos` ottenuto è completo in L2 (a, b). Di conseguenza ogni funzione f (x) L2 (a, b) può essere approssimata in media quadratica,
bene quanto si vuole, con una serie di polinomi ortogonali:
f (x) =
n Pn (x) .
n=0
Per esempio nellintervallo (1, 1) la procedura di ortogonalizzazione conduce ai polinomi di Legendre, definiti dalla formula di Rodrigues:
Pn (x) =
1 dn 2
(x 1)n ,
n
n
2 n! dx
dove si è usata la normalizzazione convenzionale:

Z 1
2
Pm (x)Pn (x)dx =
mn .
2n + 1
1
(5.32)
(5.33)
Per dimostrare che (Pm , Pn ) = 0 per n 6= m basta scrivere, per m < n

n

Z 1
Z 1
dn 2
d
n
(Pm , Pn ) = cost
dxPm (x) n (x 1) = cost
dx
P (x) (x2 1)n , (5.34)
n m
dx
dx
1
1
dove si è integrato n volte per parti tenendo conto che i contributi negli estremi si
annullano poiche
+1

dl 2
n
(x
1)
(5.35)
= 0 per ogni l n 1 ,
dxl
1
e accorgersi che
dn
Pm (x) = 0 n m + 1 .
dxn
Dalla (5.32) segue che i primi polinomi di Legendre sono
(5.36)
1
P0 (x) = 1 , P1 (x) = x , P2 (x) = (3x2 1) , .
2
Inoltre è facile dimostrare che i polinomi di Legendre obbediscono allequazione
differenziale di Legendre
(1 x2 )Pn00 (x) 2xPn0 (x) + n(n + 1)Pn (x) = 0 con n = 0, 1, 2,
135
(5.37)
che si può scrivere nella forma compatta

d
d 2
(x 1)
e n = n(n + 1) .
(5.38)
dx
dx
d 2

d
Dimostrazione: Fn (x) = dx
(x 1) dx
Pn (x) è ancora evidentemente un polinomio
di grado n; può quindi essere scritto nella forma
Lx Pn (x) = n Pn (x) con Lx =
Fn (x) =
n
X
(n)
l Pl (x)
(5.39)
l=0
(n)
con gli l proporzionali a

Z 1
Z 1
d
d
dxPl (x) (x2 1) Pn (x)
dxPl (x)Fn (x) =
dx
dx
1
1

Z 1
Z 1
d
d
2
dx (x 1) Pl (x)
=
dxFl (x)Pn (x) ,
Pn (x) =
dx
dx
1
1
dove si è ripetutamente integrato per parti; usando di nuovo la (5.39) e la (5.33) si
(n)
(n)
ricava subito che l = 0 per l < n; che n valga proprio n(n + 1) segue dal conto
esplicito della potenza pi`
u alta di Fn :
d
d 2
d
d 2
(x 1) xn + =
(x 1)nxn1 + = n xn+1 + = n(n + 1)xn +
dx
dx
dx
dx
Polinomi ortogonali su intervallo infinito
Se lintervallo è infinito non è possibile utilizzare direttamente le potenze perche esse
non sono quadrato sommabili (lintegrale (5.31) non esiste). Tuttavia sullintervallo
2
(, +) si può introdurre un fattore di convergenza ex , definire il sistema di
2
funzioni {ex /2 xn }, quadrato sommabili sullasse reale qualunque sia n = 0, 1, 2, ..., e
ortogonalizzarlo secondo il metodo appena descritto. Questo conduce ai polinomi di
Hermite, dati dalla formula di Rodrigues generalizzata
Hn (x) = (1)n ex
dn x2
e
,
dxn
n = 0, 1, 2
(5.40)
e normalizzati come segue:

Z
2
ex Hm (x)Hn (x)dx = 2n n! mn .
Lortogonalità dei polinomi dati dalla (5.40) si dimostra come per i polinomi di
Legendre. I primi polinomi di Hermite sono
136
H0 (x) = 1 , H1 (x) = 2x , H2 (x) = 4x2 2 .

Lequazione differenziale a cui obbediscono i polinomi (5.40) è lequazione differenziale di Hermite:
Hn00 (x) 2xHn0 (x) + 2nHn (x) = 0 con n = 0, 1, 2, , ,
ovvero

d x2 d
e
Hn (x) = 2nHn (x) ,
e
dx
dx
mentre le funzioni associate di Hermite
x2
(5.41)
n (x) = ex
2 /2
Hn (x) ,
n = 0, 1, 2
(5.42)
(5.43)
che formano una base ortogonale in L (R), (n , m ) = nm hn , sono soluzioni dellequazione delloscillatore armonico quantistico:
d2
+ x2 e n = 2n + 1 ,
(5.44)
dx2
che si dimostra in stretta analogia alla (5.38). Lesistenza del sistema ONC (5.43),
ovviamente numerabile, implica che anche lo spazio delle funzioni quadrato sommabili
sullintero asse reale è separabile.
Lx n (x) = n n (x) con Lx =
5.3
5.3.1
Sistemi ONC in L2 e distribuzioni.

Operatori autoaggiunti e sistemi ortogonali
Abbiamo dato qualche esempio di sistemi ONC, ma ci si può chiedere se esista un

modo pi`
u generale per costruirli. Per rispondere a questa domanda osserviamo che i
polinomi di Legendre {Pn } e le funzioni associate di Hermite {n }, che sono sistemi
ONC in un opportuno spazio L2 (a, b), condividono la proprietà di essere autofunzioni
di un operatore differenziale L, la cui espressione è data rispettivamente da eq. (5.38)
e eq. (5.44). Come vedremo tra poco, loperatore L, munito di opportune condizioni
al contorno, è in entrambi i casi autoaggiunto. Ricordiamo la definizione di aggiunto
o hermitiano coniugato A di un operatore lineare A in uno spazio vettoriale finitodimensionale su C dotato di prodotto interno (x, y) (spazio unitario):
(A y, x) = (y, Ax),
per ogni x, y nel dominio di definizione di A. In termini di matrici che rappresentano
A, (A )ij = Aji . Valgono (A + B) = A + B , (AB) = B A , (A) = A per
C, (A ) = A. Se A = A , loperatore A si dice hermitiano o autoaggiunto e
(Ay, x) = (y, Ax).
137
Un operatore hermitiano ha autovalori reali, i suoi autovettori corrispondenti a autovalori differenti sono ortogonali, e si possono scegliere gli autovettori in modo da formare
una base completa ortonormale dello spazio.
Nel caso dello spazio di Hilbert infinito dimensionale che stiamo studiando un
operatore differenziale L è autoaggiunto se
(u, Lv) = (Lu, v) (v, Lu) ,
(5.45)
u, v nel dominio di definizione di L. Questa proprietà non dipende solo da L ma

anche dal suo dominio. Infatti il dominio di L deve comprendere solo quelle funzioni
che, oltre a essere derivabili quanto basta, sono anche limitate e quadrato sommabili,
cosicche si possa ripetutamente integrare per parti in modo da poter dimostrare la
(5.45). Per esempio su L2 (R) loperatore differenziale di cui le funzioni associate di
Hermite sono autofunzioni (autovettori) è autoaggiunto:

Z +
Z +
d2
2
u (x)Lx v(x) dx =
u (x) 2 + x v(x) dx
(u, Lv) =
dx

Z +
2
d
=
2 + x2 u(x) v(x) dx .
dx
d
Lo stesso vale anche per loperatore in (5.38) (verificare per esercizio) e per Lx = i dx
,
ix
le cui autofunzioni sono esponenziali complessi e , con condizioni al contorno
periodiche u() = u(). Infatti

Z +
d
v(x) = i u (x)v(x)|+
dx u (x) i
(u, Lv) =
dx

Z +
Z +
d
d
+
dx i u (x) v(x) =
dx i u(x) v(x)
dx
dx
= (Lu, v) .
Notare limportanza delle condizioni al contorno periodiche per lannullarsi del contributo agli estremi nellintegrazione per parti. Confrontando (n , Ln ) = n (n , n )
con la sua complessa coniugata e usando la (5.45) si vede subito che
gli autovalori di un operatore autoaggiunto sono reali.
` anche facile dimostrare che
E
linsieme delle autofunzioni degli operatori autoaggiunti `
e ortogonale:
da
Ln = n n , Lm = m m ,
moltiplicando scalarmente la prima equazione per m e la seconda per n , usando la
(5.45) e sottraendo dalla prima equazione il complesso coniugato della seconda segue
0 = (n m ) (n , m )
138
da cui
(n , m ) = 0 se n 6= m .
(5.46)
Notare che queste due proprietà sono identiche a quelle degli operatori (matrici) hermitiani negli spazi unitari finito dimensionali. Il discorso invece sulla completezza dei
sistemi di autofunzioni di un operatore autoaggiunto in uno spazio infinito dimensionale è molto pi`
u delicato e non sarà affrontato qui: ci limitiamo a dire che tale sistema è
generalmente completo, come negli esempi visti qui (polinomi di Legendre in L2 (1, 1),
funzioni associate di Hermite in L2 (R), sistema trigonometrico in L2 (, )), ma ciò
non è sempre vero. Perche il sistema di autofunzioni di un operatore autoaggiunto sia
completo è talvolta necessario includervi autofunzioni generalizzate, rappresentate da
distribuzioni anziche da funzioni quadrato sommabili; daremo un esempio nel prossimo paragrafo. Vogliamo terminare questo paragrafo sottolineando limportanza delle
condizioni al contorno. Per esempio lequazione
i
d
u(x) = u(x)
dx
(5.47)
ammette la soluzione eix per ogni C; è solo la richiesta di periodicità u() = u()
che seleziona i Z, e questa ci dà una base ortogonale in L2 (, ); sullintervallo
(, +) i valori reali di sono selezionati dalla richiesta che u(x) sia limitata su
tutto lasse reale.
5.3.2
La di Dirac e cenni sulle distribuzioni
Le funzioni eix , con R, sono autofunzioni non quadrato sommabili delloperatore autoaggiunto id/dx e costituiscono una base generalizzata in L2 (R) nel senso
che ogni f (x) L2 (R) è sviluppabile, per cos` dire, in una loro combinazione lineare
continua mediante lantitrasformata di Fourier (5.26). Non si può invece dire che le
funzioni (x) = eix costituiscano un sistema ortogonale in senso stretto in L2 (R); infatti, proprio perch
sono quadrato sommabili, se cerchiamo di calcolare il prodotto
eRnon
+
scalare , = dxeix eix troviamo un integrale privo di senso. Tuttavia se
ci ricordiamo la definizione (4.7) di Trasformata di Fourier e quella (4.6) di antitrasformata, per ogni funzione f (x) sommabile e di classe C 1 su tutto R possiamo scrivere
Z +
Z +
Z +
1
1
ikx
ikx
dk F (k) e =
dk e
dy eiky f (y) .
(5.48)
f (x) =
2
2

Per la stessa ragione per cui non esiste il prodotto scalare , non è possibile
scambiare lordine di integrazione in (5.48). Se scriviamo però
1
f (x) = lim
N 2
+N
dk e
N
ikx
139
dy eiky f (y)
(5.49)
60
40
20
-0.2
0.0
-0.1
0.1
0.2
Figura 5.1: La funzione dN (t) per N = 50 e 200.
il teorema di Fubini Tonelli ci permette di scambiare lordine di integrazione ottenendo

Z +
f (x) = lim
dN (y x)f (y)dy ,
(5.50)
N
dove
1
dN (t) =
2
+N
dk eikt =
1 sin N t
.
t
(5.51)
Ovviamente non ha senso scambiare il limite con lintegrale poiche il limite puntuale
per N di dN (t) non esiste. Tuttavia si può scrivere
lim dN (t) = (t) ,
(5.52)
dove la funzione di Dirac non è una funzione ma un nuovo oggetto chiamato distribuzione; il limite è inteso in senso debole o nel senso delle distribuzioni
e significa esattamente la eq. (5.50); per dare senso alla (5.52) si deve cioè moltiplicare
per una funzione di prova (che per comodità in matematica si sceglie infinitamente
derivabile e decrescente allinfinito pi`
u rapidamente di ogni potenza; chiamiamo S lo
spazio lineare di tali funzioni di prova), integrare sulla variabile e poi fare il limite. Si
scrive allora per ogni funzione di prova g(x)
Z +
Z +
lim
dy dN (y x) g(y) =
(y x) g(y) dy = g(x) .
(5.53)
N
Lultima uguaglianza, che si può anche scrivere

Z +
(t) g(t) dt = g(0) ,
140
(5.54)
va intesa come definizione della di Dirac, dove il segno di integrale è puramente

simbolico. Simbolicamente si scrive anche
Z +
1 sin N t
1
=
(t) = lim
dk eikt ,
(5.55)
N
t
2
dove il limite (anche degli estremi di integrazione) va inteso in senso debole; si deve
cioè prima moltiplicare per una funzione di t, integrare su t, e poi integrare su k. Pi`
u
in generale la di Dirac non è che un esempio (il pi`
u importante) di distribuzione,
definita nel modo seguente:
una distribuzione (temperata) T è un funzionale lineare continuo sullo spazio S
delle funzioni di prova, cioè unapplicazione lineare continua12 di S in C:
T : SC
(5.56)
che si indica con la notazione seguente:

T : g 7 (T, g) C ,
(5.57)
dove g è una qualsiasi funzione di prova.

Il simbolo (T, g) richiama quello di prodotto scalare, e infatti in fisica si scrive abitualmente
Z +
(T, g)
T (x) g(x)dx ;
(5.58)
tuttavia lintegrale a secondo membro non sempre è un vero integrale (di Lebesgue); lo
è se, per esempio, T (x) e g(x) sono entrambe funzioni quadrato sommabili. Nella teoria
delle distribuzioni g(x) è una funzione non solo quadrato sommabile ma che va a zero
pi`
u rapidamente di 1/xn , n N (S L2 (R)); in compenso non è necessario che T (x)
sia quadrato sommabile; basta che sia localmente sommabile e non troppo crescente
per x (ci pensa g(x) ad assicurare la convergenza dellintegrale) o addirittura,
come nel caso della di Dirac, non è nemmeno necessario che sia una funzione; in
questultimo caso la (5.58) non è che un modo simbolico di scrivere la (5.57), dove la
prescrizione per associare ad ogni g S il numero complesso (T, g) è data in qualche
altro modo, per esempio per la di Dirac
: g 7 g(0) C , g(x) S
(5.59)
che i fisici scrivono nella forma (5.54). Il concetto di distribuzione generalizza quindi
quello di funzione (i russi chiamano le distribuzioni funzioni generalizzate). La di Dirac può essere rappresentata come limite debole di una successione non solo mediante
la (5.55)ma in molti altri modi. Per ogni funzione D(x) L(R) e tale che
Z +
D(x) dx = 1
(5.60)
12
Per dare un senso compiuto allaggettivo continua dovremmo definire il concetto di limite in S,
che è molto pi`
u forte del semplice limite puntuale, ma non abbiamo tempo di farlo qui.
141
si può provare
(x) = lim N D(N x) debole;
(5.61)
basta infatti moltiplicare per una funzione di prova qualsiasi g(x) e integrare per
ottenere
Z +
Z +
Z +
y
dy
(x)g(x)dx = lim
N D(N x)g(x)dx = lim
D(y)g
N
N
N
= g(0) ,
(5.62)
dove:
a) si è prima integrato e poi preso il limite, in accordo con la definizione di limite
debole;
b) si è fatto il limite sotto il segno usando il teorema (D.1) dellAppendice D, tenendo
conto che g S M > 0/|g(x)| < M, x R e quindi che |D(y)g(y/N )| <
M |D(y)|.
Tra le possibili rappresentazioni della delta di Dirac citiamo le seguenti (dove si è anche
posto = N1 ):
1.
2.
1
1
N
2
2
2 2
D(x) = ex (x) = lim e(N x) = lim ex /
0+
N
0
1
D(x) =
0
N
lim
N
0

0
= lim
1
x < 12
|x| < 12 (x) =
x > 21 .
0+
(5.63)
1
x < 2N
1
|x| < 2N
1
x > 2N
.

x
/ 2 , 2
x 2 , 2
(5.64)
3.
D(x) =
1
1
N
1
1

(x) = lim
= lim
2
2
2
0+ x + 2
N (N x) + 1
x +1
(5.65)
` importante osservare che il limite debole nella (5.61) non ha niente a che fare con il
E
limite puntuale, che ha tuttavia importanza storica poiche esso coincide con la definizione intuitiva data da Dirac alla sua funzione : una funzione che è nulla dappertutto
salvo che nellorigine, dove
vale infinito, in modo tale che il suo integrale sullasse reale
R +
valga uno (notare che N D(N x)dx = 1, N ). Ritornando al prodotto scalare

, , con (x) = eix , si può quindi scrivere

, = 2( ) ,
(5.66)
142
che generalizza a intervallo infinito la relazione di ortogonalità in L2 (, )

(l , m ) = 2lm
(5.67)
con l (x) = eilx . La di Dirac può quindi essere vista come generalizzazione della di
Kronecker a indice continuo; si confrontino anche le proprietà caratteristiche
lm = 0 , l 6= m
X
lm = 1
( ) = 0 , 6=
Z
d( ) = 1 .
(5.68)
(5.69)
Enunciamo ora unimportante propriet`

a della delta di Dirac. Sia f (x) una funzione con n zeri semplici:
f (x) = 0 per x = x1 , x2 , ..., xn ; f 0 (xi ) 6= 0 .
Allora:
(f (x)) =
n
X
(x xi )

.
df
i=1 dx
(5.70)
x=xi
Per esempio, se f (x) = x2 a2 , la (5.70) fornisce

(x2 a2 ) =
(x a) (x + a)
1
+
=
[(x a) + (x + a)] .
2|a|
2|a|
2|a|
Caso particolare importante della (5.70) è

(ax) =
1
(x) ,
|a|
a R, a 6= 0 .
(5.71)
Si calcoli come esempio

Z 3
Z 3
(x 21 )
1 1
1
(5x 2)
(5x 2)(1 2x)dx =
dx = (5 2) = .
2
2 2
4
0
0
Unaltra proprietà della di Dirac è
f (x)(x x0 ) = f (x0 )(x x0 ) ,
f C
(5.72)
come si dimostra moltiplicando ambo i membri per una qualsiasi funzione di prova g(x)
e integrando su x.
Ogni distribuzione T `
e infinitamente derivabile secondo la definizione

dg
dT
, g = T,
,
(5.73)
dx
dx
che è la generalizzazione dellidentità
Z +
Z +
dT
dg
g(x)dx =
T (x) dx
dx
dx
143
(5.74)
T
valida per ogni funzione di prova g S se T L C 1 (nellintegrazione per parti
il contributo agli estremi sparisce perche g(x) va rapidamente a zero per x ).
La definizione (5.74) di derivata generalizzata permette di derivare nel senso delle distribuzioni anche funzioni non derivabili in senso ordinario, purche localmente
sommabili. Per esempio è possibile derivare anche funzioni dotate di discontinuità di
prima specie, il cui esempo paradigmatico è la funzione di Heaviside:
d
= (x)
dx
(5.75)
nel senso delle distribuzioni13 . Infatti,

Z
Z
Z
d
0
g(x)dx =
(x)g (x)dx =
g 0 (x)dx = g(x)|+
= g(0) .
0
dx
0
(5.76)
La derivata nel senso delle distribuzioni è continua rispetto al limite debole:
dT
dTn
= lim
debole .
(5.77)
n
dx n dx
La (5.77) fornisce un altro modo per ricavare la (5.75): partendo dalla rappresentazione
T = lim Tn debole
1
1
x
+ lim arctan ;
(5.78)
2 0+

della di Heaviside, derivando ambo i membri e scambiando il limite con la derivata
(il che è sempre lecito solo nel senso delle distribuzioni) si ottiene la rappresentazione
(5.65) della .
(x) =
Ogni distribuzione temperata T ammette Trasformata di Foruier, secondo la

definizione:
(Fk (T ), Fk (g)) (T, g) , g S
(5.79)
che generalizza lanaloga identità (5.27) valida per funzioni quadrato sommabili. La
definizione (5.79) è sensata poiche, come dimostrato nel paragrafo 4.1.2, la trasformata
di Fourier manda lo spazio S in se stesso, quindi Fk (g) S, g S, e quindi Fk (T ) è
a sua volta una distribuzione temperata. Grazie alla definizione (5.79) si può calcolare
la TF anche di funzioni che non siano ne sommabili ne quadrato sommabili su R;
basta che siano localmente sommabili e non troppo crescenti14 . In particolare si può
d
In senso ordinario vale invece dx
=q.o. 0.
Non si pu`
o invece calcolare la TF di funzioni f che allinfinito crescano pi`
u di ogni potenza, per
esempio esponenzialmente; queste infatti non sono distribuzioni temperate e per dare senso a
Z +
(f, g) =
f (x)g(x)dx
(5.80)
13
14
bisognerebbe usare funzioni di prova a supporto compatto che non sono però mandate in se stesse
dalla TF
144
calcolare la TF F (k) della funzione f (x) = eik0 x (k0 R); detta

Z +
1
eikx g(x)dx
G(k) =
2
la TF di g(x), per la (5.79) abbiamo
Z
(F, G) = (f, g) =
eik0 x g(x)dx =
2G(k0 ) ,
(5.81)
(5.82)
da cui si deduce

(5.83)
F (k) Fk eik0 x = 2(k k0 ) .
Per k0 = 0 la TF di f (x) = 1 è 2(k). Viceversa scegliendo f (x) = (x x0 ) la

catena di identità
Z +
Z +
1
(F, G) = (f, g) =
(x x0 )g(x)dx = g(x0 ) =
eikx0 G(k)dk (5.84)
2
implica
1
F (k) Fk ((x x0 )) = eikx0 .
(5.85)
2
Notare che le (5.83) e (5.85) potevano anche essere dedotte (in modo rozzo e scorretto
ma mnemonicamente efficace) scrivendo semplicemente
Z +
1
F (k) =
eikx f (x)dx
(5.86)
2
ed utilizzando rispettivamente la rappresentazione (5.55) della di Dirac e la sua
definzione (5.53). Invitiamo lo studente a verificare che tale modo diventa corretto se
si usano invece i limiti deboli
eik0 x =
(x x0 ) =
lim eik0 x ex
2 /N 2
N
2 2
lim eN x ,
N
la TF della gaussiana (4.14), e la proprietà:

La TF nel senso delle distribuzioni è continua rispetto al limite debole:
T = lim TN debole F(T ) = lim F(TN ) debole .
N
(5.87)
Usando infatti le (5.87) lintegrale (5.86) acquista senso proprio, poiche il suo integrando
diventa sommabile.
Da quanto abbiamo sommariamente detto si deduce che la TF è lo strumento pi`
u
utile per risolvere equazioni differenziali a coefficienti costanti, purche venga
estesa alle distribuzioni temperate. Lunica delle tecniche discusse in questo corso che
non venga del tutto riassorbita dalla TF estesa alle distribuzioni è la Trasformata di
Laplace, che permette anche di studiare soluzioni che crescono esponenzialmente per
t (di solito per evitarle!).
145
Appendice A
Funzioni armoniche
Una funzione di due variabili g(x, y) si dice armonica se soddisfa lequazione di Laplace
42 g(x, y) = 0 ,
(A.1)
dove 42 è loperatore Laplaciano in due dimensioni

42 =
2
2
+
.
x2 y 2
Teorema: Se f (z) = u(x, y) + iv(x, y) è una funzione analitica (e le funzioni u e v

sono di classe C 2 1 ) , le funzioni u(x, y) e v(x, y) sono armoniche.
Dimostrazione: Se f (z) è analitica, u e v soddisfano le condizioni di CauchyRiemann:
v(x, y)
u(x, y)
=
x
y
u(x, y)
v(x, y)
=
y
x
(A.2)
(A.3)
Derivando la (A.2) rispetto a x e la (A.3) rispetto a y si ottiene

2 u(x, y)
2 v(x, y)
=
x2
xy
2
u(x, y)
2 v(x, y)
2 v(x, y)
=
.
y 2
yx
xy
Nellultima equazione è lecito scambiare lordine di derivazione perche v(x, y) è di classe
C 2 . Sottraendo membro a membro le precedenti equazioni si ottiene:
1
In realt`
a questa condizione è sempre soddisfatta perche ogni funzione analitica è infinitamente
derivabile.
146
2 u(x, y) 2 u(x, y)
+
= 4u(x, y) = 0 .
x2
y 2
Analogamente, derivando la (A.2) rispetto a y e la (A.3) rispetto a x e sottraendo
membro a membro si ottiene
2 v(x, y) 2 v(x, y)
+
= 4v(x, y) = 0 .
x2
y 2
[q.e.d.]
Data una funzione u(x, y) armonica in una certa regione del piano (x, y) è possibile
costruire (a meno di una costante) la corrispondente funzione armonica v(x, y) tale che
f (z) = u(x, y) + iv(x, y) sia analitica. Infatti, nota u, le condizioni di Cauchy-Riemann
ci consentono di ricavare le derivate parziali vx0 e vy0 e da queste, integrando, la funzione
v(x, y).
Esempio 1
Sia
u(x, y) = cos xey .
Dimostrare che u(x, y) è armonica e costruire la corrispondente funzione analitica f (z).
Per dimostrare che u è armonica occorre dimostrare che essa soddisfa lequazione
di Laplace (A.1):
u(x, y)
= sin xey ,
x
2 u(x, y)
= cos xey ,
x2
u(x, y)
= cos xey
y
2 u(x, y)
= cos xey
y 2
Pertanto 42 u(x, y) = 0 e u è armonica. Dalle condizioni di CR si ricava che

v(x, y)
u(x, y)
=
= sin xey
y
x
v(x, y)
u(x, y)
=
= cos xey .
x
y
Integrando la (A.4) rispetto a y si ottiene
Z
v(x, y) = sin xey dy + g(x) = sin xey + g(x) .
147
(A.4)
(A.5)
Sostituendo questultima nella (A.5)

v(x, y)
= cos xey + g 0 (x) = cos xey .
x
Pertanto
g 0 (x) = 0 g(x) = costante = K v(x, y) = sin xey + K
e la funzione f (z) cercata è (ponendo la costante K = 0)
f (z) = u(x, y) + iv(x, y) = cos xey + i sin xey = eixy = ei(x+iy) = eiz ,
analitica in tutto C.
Esempio 2
La funzione
u(x, y) = x2 y 2
è armonica. Infatti
u00xx = 2 ,
u00yy = 2 42 u(x, y) = 0
Dalle condizioni di CR segue che:

vy0 = u0x = 2x v(x, y) = 2xy + (x)
e
vx0 = u0y = 2y 2y + 0 (x) = 2y (x) = K .
Pertanto
v(x, y) = 2xy + K f (z) = x2 y 2 + 2ixy + iK = (x + iy)2 + iK = z 2 + iK ,
Esempio 3
La funzione
u(x, y) =
x2
x
+ y2
è armonica in R2 {0}. Infatti:

u0x
y 2 x2
,
= 2
(x + y 2 )2
u00xx
148
2x(x2 3y 2 )
=
(x2 + y 2 )3
KR
Figura A.1: Mappa conforme.
2x(x2 3y 2 )
= u00xx .
(x2 + y 2 )3
u0y =
2xy
,
(x2 + y 2 )2
u0y
2xy
v(x, y) =
= 2
(x + y 2 )2
u00yy =
Dalle CR:
vx0
2xy
dx + g(y)
+ y 2 )2
y
+ g(y) .
= 2
x + y2
(x2
Per fissare g(y) usiamo laltra condizione di CR:

1
2y 2
+
+ g 0 (y)
x2 + y 2 (x2 + y 2 )2
y 2 x2
y 2 x2
0
0
+ g (y) = ux = 2
=
(x2 + y 2 )2
(x + y 2 )2
g(y) = K .
vx0 =
Quindi
f (z) = u(x, y) + iv(x, y) =
=
x2
z
1
+
iK
=
+ iK ,
|z|2
z
x
y
i 2
+ iK
2
+y
x + y2
KR.
Esiste uno stretto legame tra funzioni analitiche e armoniche, per le quali valgono
teoremi simili a quelli che abbiamo visto nella prima parte. Le funzioni armoniche
149
~ E
~ = 0 e il
compaiono in molti problemi fisici. In elettrostatica nel vuoto abbiamo
~ = V
~ , da cui V = 0. Anche un fluido
campo elettrico può essere scritto come E
incompressibile in assenza di vortici è descritto da un potenziale di flusso che soddisfa
lequazione di Laplace, e quindi è una funzione armonica.
Dalle condizioni di Cauchy-Riemann abbiamo anche
~ v
~ = u v + u v = 0
u
x x y y
(A.6)
e si dice che u e v sono funzioni armoniche coniugate. Poichè il gradiente è perpendicolare alle linee su cui la funzione è costante, questo implica che le curve di u = cost
sono localmente perpendicolari a quelle di v = cost, vedi Fig.(A.1), esattamente come
perpendicolari sono le rette x = cost e y = cost nel piano complesso di partenza. Attraverso la funzione analitica f (z) = u(x, y) + iv(x, y) si realizza pertanto una mappa
conforme (che rispetta localmente gli angoli) del piano complesso (x, y) nel piano (u, v).
Per concludere notiamo il seguente teorema: Il campo vettoriale ~a = v(x, y)~i +
~ ~a = 0) e irrotazionale (
~ ~a = 0) se e solo se f (z) =
u(x, y) ~j è solenoidale (
u + iv è analitica. La verifica segue dalle Cauchy-Riemann. Un campo solenoidale e
irrotazionale definisce una funzione analitica e viceversa.
150
Appendice B
Corollari della rappresentazione
integrale di Cauchy
1. Principio della media
Il principio della media si ricava dalla rappresentazione integrale di Cauchy (1.22) nel
caso in cui la curva sia una circonferenza. Se f (z) è una funzione analitica in una
regione S semplicemente connessa e C è una circonferenza di raggio r e centro z0 S
interna a S:
C = {z S / |z z0 | = r} ,
allora
1
f (z0 ) =
2
f [z()]d
z C .
(B.1)
Dimostrazione: Applichiamo la rappresentazione integrale di Cauchy con = C e

passiamo a coordinate polari (z C)
z() = z0 + rei dz = irei d = i(z z0 )d
da cui segue
1
i
f (z0 ) =
2i
Z
0
z z0
1
f [z()]d =
z z0
2
f [z()] d .
0
[q.e.d.]
2. Teorema
Ne la parte reale ne la parte immaginaria di una funzione f (z) = u(x, y) + iv(x, y)
analitica in un dominio D s.c. possono avere estremi allinterno di D.
151
Dimostrazione: dimostriamo il teorema per la parte reale u(x, y). La dimostrazione

è analoga per v(x, y). Sia z0 un punto interno a D, tale cioè che esista un intorno I(z0 )
tutto contenuto in D. Mostriamo che possono esistere solo due alternative:
I(z0 ) D / z I(z0 ) u(z) = u(z0 ) ,
a)
esiste cioè un intorno di z0 in cui la funzione u(x, y) è costante;

b)
I(z0 ) D z1 , z2 I(z0 ) / u(z2 ) < u(z0 ) < u(z1 ) ,
cioè in qualunque intorno di z0 la funzione u assume valori sia maggiori che minori di
u(z0 ).
Neghiamo infatti lalternativa b), ammettiamo cioè che esista un intorno I(z0 ) D
tale che z I(z0 ) sia, per esempio, u(z) u(z0 ). Dalla parte reale del principio della
media (B.1) segue allora che
Z 2
1
u(z0 ) =
u[z()]d .
(B.2)
2 0
Daltra parte è vera la relazione
1
u(z0 ) =
2
u(z0 )d .
(B.3)
Uguagliando le (B.2) e (B.3) si ottiene

Z 2
{u[z()] u(z0 )}d = 0 ,
0
che implica, poiche la funzione integranda è per ipotesi continua e non negativa,
u(z) = u(z0 )
ovvero lalternativa a).
Laltro modo di negare lalternativa b) è di affermare che esiste un intorno I(z0 ) D
tale che z I(z0 ) sia u(z) u(z0 ));allora dalle (B.2) e (B.3) deduciamo di nuovo che
Z
{u(z0 ) u[z()]}d = 0;
0
poiche la funzione integranda è per ipotesi continua e non negativa, ne segue ancora
lalternativa a):
u(z) = u(z0 ) .
[q.e.d.]
152
3. Teorema
Sia f (z) una funzione analitica in un dominio D s.c. . Il suo valore assoluto |f (z)| non
può avere massimi allinterno di D e può avere minimi solo nei punti in cui f (z) = 0.
Dimostrazione: Per quanto riguarda il massimo, la dimostrazione è perfettamente
analoga alla precedente. Supponiamo che:
I(z0 ) D / |f (z)| |f (z0 )|
z I(z0 ) .
Per ogni circonferenza C centrata in z0 e interna a I(z0 ) vale, per il principio della
media, la seguente relazione

Z 2
Z 2

1
1

f [z()] d
|f [z()]| d .
(B.4)
|f (z0 )| =
2 0
2 0
` vero inoltre che
E
1
|f (z0 )| =
2
|f (z0 )| d .
(B.5)
Sottraendo membro a membro le (B.4) e (B.5) si ottiene

Z 2
{|f [z()]| |f (z0 )|} d 0
0
Poiche lintegrando è per ipotesi una funzione continua e non positiva, ne segue che
|f (z)| = |f (z0 )|
per ogni z C. Variando il raggio r della circonferenza C in modo da coprire tutto
lintorno I(z0 ) si dimostra che |f (z)| = |f (z0 )| per ogni z I(z0 ).
Abbiamo cos` dimostrato che |f (z)| non può avere un massimo allinterno di D.
Per dimostrare che |f (z)| non può avere minimi, se non nei punti in cui si annulla,
consideriamo la funzione g(z) = 1/f (z), analitica in D esclusi gli zeri zi di f (z). In
questi punti |f (zi )| è ovviamente minima. Se z 6= zi , i minimi di |f (z)| corrispondono
ai massimi di |g(z)|. Ma, come si è appena dimostrato, la funzione |g(z)|, analitica in
D, non può avere massimi in D, e quindi la |f (z)| non può avere minimi.
[q.e.d.]
4. Teorema
Sia f (z) una funzione analitica in un dominio semplicemente connesso S, S una
curva di Jordan di lunghezza l, f (z) limitata sulla curva e M = maxz |f (z)| il
suo valore massimo su , z0 un punto appartenente alla regione interna a e =
minz |z z0 | la distanza minima della curva dal punto z0 . Sotto queste ipotesi:
153
a)
|f (z0 )|
Ml
2
(B.6)
b)
n

d f (z)

dz n
n!
z=z0
Ml
2 n+1
n = 1, 2, ...
(B.7)
Dimostrazione: La dimostrazione segue immediatamente dalle rappresentazioni di

Cauchy (1.22) e (1.23) e dalla disuguaglianza di Darboux (1.11):

I
1
f (z)

dz
|f (z0 )| =
2i z z0

f (z)
1
l
max
2i z z z0
a)
b)
n

d f (z)

dz n
z=z0
1 Ml
2i

I

n!
f (z)

dz
=
2i (z z0 )n+1
n! M l
2 n+1
[q.e.d.]
5. Teorema di Liouville
Una funzione analitica e limitata in tutto il piano complesso C è necessariamente
costante.
Dimostrazione: Poiche f (z) è limitata in C, esiste un M reale tale che |f (z)| M ,
z C. Allora z0 C possiamo applicare il teorema precedente (B.7) nel caso n = 1:

df (z)
Ml

.
dz
2 2
z=z0
Se scegliamo come una circonferenza centrata in z0 di raggio r (l = 2r, = r)
otteniamo

df (z)
M 2r
M

=
.
dz
2
2r
r
z=z0
154
Ora, poiche f (z) è regolare e limitata in tutto il piano complesso, si può scegliere r
arbitrariamente grande, rendendo il rapporto M/r piccolo quanto si vuole; quindi

df (z)

=0
dz
z=z0
z0 C
da cui
df (z)
=0
dz
z C f (z) = costante
z C .
[q.e.d.]
Se definiamo una funzione intera come una funzione regolare in tutto il piano
complesso (cioè priva di singolarità al finito), il teorema di Liouville afferma che una
funzione intera e limitata è costante.
N.B. Lo stesso teorema non vale nel campo reale. Infatti esistono funzioni f (x) di
variabile reale non costanti che sono infinitamente derivabili e limitate in tutto R. Per
esempio le funzioni
f (x) =
1
,
1 + x2
f (x) = ex
sono limitate (f (x) 1) e infinitamente derivabili. Tuttavia le corrispondenti funzioni

nel campo complesso non sono limitate in C. Infatti
f (z) =
1
1 + z2
non è regolare ne limitata in z = i. Invece la funzione

f (z) =
1
1 + |z|2
è limitata in C (f (z) 1) ma non è analitica (poiche u(x, y) = (1 + x2 + y 2 )1

e v(x, y) = 0, le condizioni di Cauchy-Riemann non sono soddisfatte). Per quanto
riguarda
f (z) = ez
è regolare in tutto C ma per z = iy, con y reale,

f (iy) = ey
non è limitata. Invece la funzione

2
f (z) = e|z| = ez
è limitata ma non analitica (si ricordi che z non è analitica).

155
Appendice C
Equazioni differenziali del
secondordine
La maggior parte dei problemi in fisica è formulata in termini di equazioni differenziali,
che sono molto spesso equazioni differenziali alle derivate parziali, che coinvolgono cioè
derivate rispetto a pi`
u di una variabile. Tra queste, le equazioni che si incontrano pi`
u
frequentemente sono:
1) Lequazione di Laplace:
(~r) = 0 ,
(C.1)
dove loperatore Laplaciano (in tre dimensioni) è

= 2 =
2
2
2
+
+
.
x2 y 2 z 2
Lequazione di Laplace compare spesso nello studio di fenomeni elettromagnetici,

nellidrodinamica, nella propagazione del calore e nello studio della gravitazione.
2) Lequazione di Poisson:
(~r) =
,
0
(C.2)
che è una generalizzazione dellequanzione di Laplace (C.1) in presenza di una

sorgente.
3) Lequazione di Helmholtz, o equazione delle onde
(~r) + k 2 (~r) = 0
(C.3)
e lequazione di diffusione indipendente dal tempo

(~r) k 2 (~r) = 0 .
156
(C.4)
Queste equazioni si incontrano nello studio di diversi fenomeni, come la propagazione di onde elastiche nei solidi, la propagazione del suono e lacustica, le onde
elettromagnetiche e i reattori nucleari.
4) Lequazione delle onde dipendente dal tempo, o equazione di dAlembert,
(~r, t) = 0 ,
(C.5)
dove si è introdotto loperatore dalembertiano:

=
2
2
2
2
2
=
.
t2 x2 y 2 z 2
t2
5) Lequazione del potenziale scalare:

(~r, t) =
0
6) Lequazione di Klein-Gordon
(~r, t) = 2 (~r, t) .
(C.6)
7) Lequazione di Schroedinger
(~r, t)
~2
(~r, t) + V (~r)(~r, t) = i~
,
2m
t
che è alla base della meccanica quantistica non relativistica.

8) Le equazioni di Maxwell, che sono un sistema di equazioni alle derivate parziale
accoppiate per i campi elettrico e magnetico.
9) Lequazione di Dirac, che governa la meccanica quantistica relativistica.
Tutte queste equazioni possono essere scritte nella forma
H = F ,
dove H è un operatore differenziale,

H=H
, , , , x, y, z ,
x y z t
F è una funzione nota e è la funzione incognita (scalare o vettoriale). Inoltre le
equazioni (1,3,4,6,7,9) sono lineari: questo significa che, se 1 e 2 sono soluzioni
dellequazione differenziale, qualsiasi loro combinazione lineare
= a1 1 + a2 2
157
ne è ancora soluzione. Le equazioni (1-8) sono tutte equazioni del secondo ordine: esse
contengono cioè derivate di ordine massimo 2 (in realtà le equazioni di Maxwell sono
del primordine, ma contengono due funzioni incognite e possono essere ricondotte,
eliminando una delle due incognite, a equazioni del secondordine). Esistono diversi
metodi per la soluzione di equazioni differenziali alle derivate parziali del secondordine:
1) metodo di separazione delle variabili;
2) metodo delle funzioni di Green;
3) metodo delle trasformate integrali di Fourier e di Laplace;
4) metodi numerici.
C.1
Metodo di separazione delle variabili
Il metodo di separazione delle variabili consiste nel separare una equazione differenziale alle derivate parziali in n variabili in n equazioni differenziali ordinarie (contenenti
cioè derivate totali), ciascuna corrispondente a una variabile. Ogni separazione introduce una costante arbitraria, detta costante di separazione. Se le variabili sono n, si
ottengono n 1 costanti di separazione, determinate dalle condizioni al contorno del
problema. Illustriamo il funzionamento di questo metodo con un esempio: la soluzione
dellequazione di Helmholtz (C.3). In coordinate cartesiane lequazione (C.3) diventa:
2 (x, y, z) 2 (x, y, z) 2 (x, y, z)
+
+
+ k 2 (x, y, z) = 0 .
x2
y 2
z 2
(C.7)
Il metodo di separazione delle variabili consiste nel cercare una soluzione fattorizzata
del tipo
(x, y, z) = X(x)Y (y)Z(z) ,
(C.8)
dove X(x) dipende solo dalla variabile x, Y (y) solo dalla y e Z(z) solo dalla z. Non
è detto, in generale, che una soluzione di questo tipo esista, ma, se esiste, il procedimento è giustificato. Se non esiste, si dovrà risolvere lequazione con un altro metodo.
Sostituiamo dunque la (C.8) nellequazione (C.7):
d2 Y (y)
d2 Z(z)
d2 X(x)
+ X(x)Z(z)
+ X(x)Y (y)
+ k 2 X(x)Y (y)Z(z) = 0
Y (y)Z(z)
2
2
2
dx
dy
dz
e dividiamo per X(x)Y (y)Z(z), supponendo che 6= 0 (soluzione banale dellequazione
differenziale):
1 d2 X(x)
1 d2 Y (y)
1 d2 Z(z)
+
+
+ k2 = 0 .
2
2
2
X(x) dx
Y (y) dy
Z(z) dz
158
Si noti che le derivate sono derivate ordinarie, non parziali. Separiamo ora i termini
che dipendono da x da quelli che non ne dipendono:
1 d2 X(x)
1 d2 Y (y)
1 d2 Z(z)
=
k2 .
X(x) dx2
Y (y) dy 2
Z(z) dz 2
(C.9)
Ora, il primo membro è una funzione della sola variabile x, mentre il secondo membro
dipende solo da y e z. Poichè x, y e z sono variabili indipendenti, questo è possibile solo
se entrambi i membri della (C.9) sono uguali a una costante, la costante di separazione,
che indicheremo con l2 . Otteniamo cos` due equazioni:
1 d2 X(x)
= l2
X(x) dx2
1 d2 Y (y)
1 d2 Z(z)
+
+ k 2 = l2 .
Y (y) dy 2
Z(z) dz 2
Consideriamo ora la seconda di queste equazioni e ripetiamo il procedimemto; separiamo i termini dipendenti da y:
1 d2 Z(z)
1 d2 Y (y)
=
k 2 + l2 .
Y (y) dy 2
Z(z) dz 2
Questa equazione può essere verificata solo se ambo i suoi membri sono uguali a una
costante di separazione, m2 :
1 d2 Y (y)
= m2
2
Y (y) dy
1 d2 Z(z)
= k 2 + m2 + l2 = n2 ,
2
Z(z) dz
dove abbiamo introdotto, per motivi di simmetria formale, la costante n2 , legata a l,
m e k dalla relazione
l2 + m2 + n2 = k 2 .
Riassumendo, abbiamo trasformato lequazione differenziale alle derivate parziali in
n=3 variabili
2 (x, y, z) 2 (x, y, z) 2 (x, y, z)
+
+
+ k 2 (x, y, z) = 0
x2
y 2
z 2
in n=3 equazioni differenziali ordinarie
1 d2 X(x)
= l2
X(x) dx2
1 d2 Y (y)
= m2
Y (y) dy 2
1 d2 Z(z)
= n2
2
Z(z) dz
159
e abbiamo introdotto n 1=2 costanti di separazione l e m (la terza costante, n, non

è indipendente dalle altre due). La nostra soluzione sarà caratterizzata da tre indici
l, m, n:
lmn (x, y, z) = Xl (x)Ym (y)Zn (z)
(C.10)
e varrà per qualunque scelta delle costanti l, m, n purchè l2 +m2 +n2 = k 2 . La soluzione
generale sarà una combinazione lineare delle (C.10):
X
(x, y, z) =
almn lmn (x, y, z) ,
l,m,n
dove i coefficienti almn saranno determinati dalle condizioni al contorno. Risolviamo ora
lequazione di Helmholtz in coordinate sferiche. Cerchiamo una soluzione fattorizzata
del tipo
(r, , ) = R(r)()() .
Lespressione delloperatore laplaciano in coordinate sferiche è:

1 2
1
2
sin
r
+
sin
+
.
= 2
r sin
r
r
sin 2
(C.11)
(C.12)
Le (C.11) e (C.12), sostituite nellequazione (C.3), forniscono:

d2
1
d
d
1
1 d
2 dR
= k 2 .
r
+
sin
+
Rr2 dr
dr
r2 sin d
d
r2 sin2 d2
Si noti che le derivate sono diventate derivate ordinarie, perchè agiscono su funzioni di
una sola variabile. Moltiplicando per r2 sin2 possiamo isolare i termini dipendenti da
:

1 d2
1 d
1
d
d
2
2
2
2 dR
= r sin k
r
2
sin
. (C.13)
d2
Rr2 dr
dr
r sin d
d
Questa equazione è una uguaglianza fra una funzione che dipende solo da e una che
dipende solo da r e : lunica soluzione possibile è che ambo i membri della (C.13)
siano uguali a una costante, che indicheremo con m2 :
1 d2
= m2
d2

1 d
1
d
d
m2
2 dR
r
+
sin
= k 2 .
Rr2 dr
dr
r2 sin d
d
r2 sin2
(C.14)
Moltiplicando ora la seconda equazione per r2 e riarrangiando i termini si ottiene:

1 d
1
d
d
m2
2 dR
2 2
r
+k r =
sin
+
.
R dr
dr
sin d
d
sin2
160
Le variabili r e sono cos` separate. Se uguagliamo il primo e il secondo membro ad

ununica costante Q otteniamo:

QR
1 d
2 dR
2
r
+
k
R
= 0
(C.15)
r2 dr
dr
r2

d
m2
1 d
sin
+ Q = 0 .
(C.16)
sin d
d
sin2
Abbiamo ottenuto cos` tre equazioni differenziali ordinarie (C.14), (C.15) e (C.16), con
lintroduzione di 2 costanti di separazione, m2 e Q. La soluzione di queste equazioni
è discussa nel capitolo 2.1.1. La soluzione generale dellequazione di Helmholtz in
coordinate sferiche avrà la forma:
X
(r, , ) =
aQm RQ (r)Qm ()m () .
Q,m
Mediante il metodo di separazione delle variabili ci siamo quindi ricondotti a equazioni

differenziali ordinarie del secondordine, delle quali ci occupiamo qui.
161
Appendice D
Propriet`
a dellintegrale di Lebesgue
Elenchiamo qui di seguito (senza dimostrarle) alcune proprietà dellintegrale di Lebesgue che lo rendono spesso molto pi`
u maneggevole dellintegrale di Riemann.
Teorema di Lebesgue sullo scambio di limite con integrale (per successioni).
Se
lim fn (x) q.o.
= f (x) ,
con le fn (x) sommabili, ed esiste una F (x) sommabile tale che n N, |fn (x)| q.o.
F (x), allora anche f (x) è sommabile e si può scambiare il limite con lintegrale:
Z
Z
lim fn (x) dx
fn (x) dx =
lim
f (x) dx .
a n
(D.1)
La seconda proprietà è molto simile:

Scambio di limite con integrale (per integrali dipendenti da un parametro).
Se in y0 la funzione f (x, y), intesa come funzione di y, è continua per quasi ogni
x (a, b), cioè se vale
lim f (x, y) =q.o. f (x, y0 ),
(D.2)
yy0
ed esiste una F (x) sommabile tale che in un opportuno intorno di y0 valga

|f (x, y)| q.o. F (x), allora si può scambiare il limite con lintegrale:
Z
lim
yy0
Z
f (x, y) dx =
Z
lim f (x, y) dx
a yy0
Anche la terza proprietà vale sotto condizioni analoghe:
162
f (x, y0 ) dx .
a
(D.3)
Derivazione sotto il segno.

Posto
Z
f (x, y) dx ,
G(y) =
a
se in un opportuno intorno di y0 e per quasi ogni x (a, b) la derivata parziale

fy0 (x, y) esiste e vale
F (x) con F (x) sommabile ,
|fy0 (x, y)| q.o.
allora si può derivare sotto il segno:

Z b
dG
fy0 (x, y0 ) dx .
=
dy y=y0
a
(D.4)
Enunciamo infine
Teorema di Fubini Tonelli sullo scambio dellordine di integrazione.
Se esiste almeno uno degli integrali
dx|f (x, y)|

a
Z
dy
dy|f (x, y)| ,
dx
a
allora vale
Z
dxf (x, y)
dy
dyf (x, y) =
dx
a
(D.5)
Limportanza di questi quattro teoremi diventa evidente se si osserva che essi valgono sia per intervallo finito che infinito e sia per funzioni limitate che non limitate;
gli analoghi teoremi validi per integrali di Riemann, specie per integrali di Riemann
impropri, richiedono condizioni sufficienti estremamente pi`
u restrittive e difficili da
verificare.
Naturalmente nel caso di intervallo finito anche per lintegrale di Lebesgue si
possono avere ulteriori semplificazioni: se e solo se lintervallo è finito, F (x) = M
costante è sommabile; si può quindi dire, come caso particolare della proprietà (D.1),
che condizione sufficiente affinche valga la (D.1) è che le fn (x) siano q.o. uniformemente
limitate, cioè che esista una costante M , indipendente da n, che le maggiori tutte (in
modulo) quasi ovunque.
Affermazioni analoghe si possono fare per le proprietà (D.3) e (D.4).
163

Introduzione Ai Metodi Matematici Della Fisica

Caricato da

Copyright:

Formati disponibili

Introduzione Ai Metodi Matematici Della Fisica

Caricato da

Informazioni sul documento

Descrizione originale:

Copyright

Formati disponibili

Condividi questo documento

Condividi o incorpora il documento

Opzioni di condivisione

Hai trovato utile questo documento?

Questo contenuto è inappropriato?

Copyright:

Formati disponibili

Introduzione Ai Metodi Matematici Della Fisica

Caricato da

Copyright:

Formati disponibili

Introduzione ai

METODI MATEMATICI DELLA FISICA

Dispense del corso 2013

Introduzione allanalisi armonica

B Corollari della rappresentazione integrale di Cauchy

Equazioni differenziali del secondordine

Richiami sui numeri complessi

Un numero complesso z `e una coppia ordinata di numeri reali

dove i `e lunit`a immaginaria

Linsieme dei numeri complessi `e un campo, indicato

z2 = a2 + ib2 z1 + z2 = (a1 + a2 ) + i(b1 + b2 ) .

Esiste lelemento neutro 0 C, tale che

Per ogni z C esiste lelemento inverso z C, tale che

Per ogni z C, z 6= 0 esiste lelemento inverso z 1 C, tale che

Quindi C {0} `e un gruppo abeliano rispetto alla moltiplicazione, con elemento

Figura 1.1: Rappresentazione cartesiana del numero complesso z

Rappresentazione geometrica dei numeri complessi

Graficamente r `e il modulo del vettore ~r congiungente lorigine con il punto z, e

Figura 1.2: Rappresentazione polare del numero complesso z

Il complesso coniugato di un numero complesso z `e un numero complesso z C

Funzioni reali di variabile complessa

Una funzione reale di variabile complessa

`e unapplicazione che associa un numero reale f (z) ad ogni z E, con E sottoinsieme

ossia, ricordando la definizione di limite, se

dove I (z0 ) `e un intorno di raggio del punto z0 :

sono funzioni reali di variabile complessa, continue in tutto il piano complesso.

Figura 1.4: Discontinuit`a dellargomento di z.

Nel caso b) invece (z) `e discontinua sul semiasse reale positivo.

Funzioni complesse di variabile complessa

Una funzione complessa di variabile complessa

`e unapplicazione che associa un numero complesso f (z) ad ogni z E, con E

f (z) = z 2 = (x + iy)2 = x2 y 2 + 2ixy

Derivata di una funzione complessa di variabile complessa

Figura 1.5: Direzioni dellincremento h nel rapporto incrementale

La funzione f `e derivabile se tutte queste direzioni danno lo stesso risultato per

Consideriamo una funzione

esiste ed `e indipendente dalla direzione di h = hx + ihy . In particolare si potr`a scegliere

La derivata di f (z) `e quindi definita univocamente indipendentemente dalla direzione

u0x (x, y) + ivx0 (x, y)

Definizione: Una funzione f (z)

si dice analitica (o regolare, o olomorfa) in un punto z0 se esiste un intorno I(z0 )

u0x = u0y = vx0 = vy0 = 0

Le funzioni u e v sono continue e derivabili in C {0}. Le condizioni di CR

sono soddisfatte. La funzione f (z) `e quindi analitica in C {0}. La sua derivata `e

Vale pertanto la stessa regola di derivazione valida in campo reale

La dimostrazione segue banalmente dalla definizione di derivata o dalla (1.6).

Corollario: le funzioni razionali di z sono analitiche in tutto il piano complesso

sono funzioni ovunque analitiche. Per la propriet`a 3) le funzioni razionali (rapporto di

Derivando rispetto a x e a y si ottiene

Integrazione in campo complesso

Una curva nel piano complesso `e una applicazione continua

dove J `e un intervallo reale limitato e chiuso:

dove S `e la regione interna a .

Integrali in campo complesso

Integrale di una funzione di variabile reale a valori complessi Consideriamo

Definiamo lintegrale di w(t) in t

Lintegrale esiste se la funzione w(t) `e continua o se ha un numero finito di discontinuit`a

w(l )(tl tl1 )

u(l )(tl tl1 ) + i