Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi...
Transcript of Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi...
![Page 1: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/1.jpg)
Metodi Numerici per Equazioni Ellittiche
Vediamo ora di descrivere una tecnica per la risoluzione numerica della
piu semplice equazione ellittica lineare, l’Equazione di Laplace:
uxx + uyy = 0, (x, y) ∈ Ω ⊂ R2. (1)
Se una funzione u(x, y) e di classe C2 in un determinato sottoinsieme
Ω di R2 ed e una soluzione di (1) nello stesso Ω allora prende il nome
di funzione armonica.
Le proprieta di queste funzioni sono state ampiamente studiate poiche
trovano applicazione in numerosi campi della fisica applicata.
1
![Page 2: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/2.jpg)
Tra le principali proprieta di cui godono le funzioni armoniche la prin-
cipale e sicuramente la seguente.
Teorema[Principio del massimo-minimo] Sia Ω una regione limita-
ta e semplicemente connessa e Γ la sua frontiera. Sia Ω = Ω∪Γ. Se
u(x, y) e armonica su Ω e continua su Ω, allora u(x, y) assume i suoi
valori massimo e minimo su Γ.
2
![Page 3: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/3.jpg)
Il Problema di Dirichlet
L’equazione di Laplace puo essere associata ad un problema di Dirich-
let quando, assegnata una funzione f(x, y) di classe C2(Γ), si cerca
una funzione u(x, y) tale che:
1. u(x, y) e continua su Ω ∪ Γ;
2. u(x, y) = f(x, y) per ogni (x, y) ∈ Γ;
3. u(x, y) e armonica nell’insieme Ω.
3
![Page 4: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/4.jpg)
y
xO
ΓΩ
4
![Page 5: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/5.jpg)
y
xO
u(x, y) = f(x, y)Ω
5
![Page 6: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/6.jpg)
y
xO
u(x, y) = f(x, y)uxx + uyy = 0
6
![Page 7: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/7.jpg)
Il Problema di Neumann
L’equazione di Laplace puo essere associata ad un problema di Neu-
mann quando, al posto della condizione 2., si impone che sia
∂u
∂n= f(x, y)
cioe sia assegnata la derivata normale di u(x, y) rispetto alla curva Γ.
Ricordiamo che se nT = (nx, ny), e il vettore normale allora
∂u
∂n=
(
∇Tu)
· n = nx∂u
∂x+ ny
∂u
∂y.
7
![Page 8: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/8.jpg)
Consideriamo ora la risoluzione dell’equazione di Laplace prendendo
Ω uguale al rettangolo [a, b] × [c, d]. In questo caso un metodo e
quello di approssimare l’operatore differenziale dopo avere suddiviso
in modo opportuno l’insieme Ω. Infatti si suddivide l’intervallo [a, b]
in N parti uguali sull’asse x e M sull’asse y ottenendo la reticolazione
di Ω mediante i seguenti punti:
xi = xi−1 + h = a + ih i = 0,1, . . . , N
yj = yj−1 + k = c + jk j = 0,1, . . . , M
dove
h =b − a
N, k =
d − c
M.
8
![Page 9: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/9.jpg)
Ω
Γ(a, c)
(a, d)
(b, c)
(b, d)
9
![Page 10: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/10.jpg)
Abbiamo cosı ottenuto un insieme discreto di punti del piano
RN+1,M+1 =
(xi, yj) ∈ R2|xi = a + ih, i = 0, N, yj = c + jk, j = 0, M
.
La risoluzione numerica del problema di Dirichlet associato consiste
nell’approssimare opportunamente la funzione u(x, y) nei punti ap-
partenenti all’insieme RN+1,M+1.
L’idea alla base del metodo e quella di approssimare le derivate parziali
seconde nei punti del reticolo RN+1,M+1 e imporre che tali approssi-
mazioni soddisfino l’equazione di Laplace.
Poniamo innanzitutto
uij ≃ u(xi, yj), i = 0,1, . . . , N, j = 0,1, . . . , M
10
![Page 11: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/11.jpg)
Per approssimare le derivate parziali seconde uxx(x, y), e uyy(x, y), nel
punto (xi, yj) consideriamo il seguente reticolo:
(xi, yj)
(xi, yj+1)
(xi, yj−1)
(xi+1, yj+1)
(xi−1, yj−1) (xi+1, yj−1)
(xi+1, yj)(xi−1, yj)
(xi−1, yj+1)
11
![Page 12: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/12.jpg)
Per approssimare la derivata parziale seconda uxx(x, y) possiamo sup-
porre che questa sia una derivata del secondo ordine ordinaria rispetto
ad x, considerando la y costante.
Consideriamo i punti (xi−1, yj), (xi, yj) e (xi+1, yj) e, applicando la
formula per l’approssimazione di f ′′(ti), risulta:
∂2u
∂x2(xi, yj) ≃
ui+1,j − 2ui,j + ui−1,j
h2.
12
![Page 13: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/13.jpg)
Analogamente per approssimare uyy(xi, yj) si considera come fosse
una derivata ordinaria fatta rispetto a y, tenendo la x costante e
coinvolgendo i punti del reticolo (xi, yj−1), (xi, yj) e (xi, yj+1) :
∂2u
∂y2(xi, yj) ≃
ui,j+1 − 2ui,j + ui,j−1
k2.
13
![Page 14: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/14.jpg)
Tenendo presente che la funzione u(x, y) e nota sul bordo del ret-
tangolo alcune delle approssimazioni non devono essere calcolate,
infatti:
u0,j = u(x0, yj) = u(a, yj) = f(a, yj), j = 0, . . . , M
ui,0 = u(xi, y0) = u(xi, c) = f(xi, c), i = 0, . . . , N
uN,j = u(xN , yj) = u(b, yj) = f(b, yj), j = 0, . . . , M
ui,M = u(xi, yM) = u(xi, d) = f(xi, d), i = 0, . . . , N.
14
![Page 15: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/15.jpg)
Adesso possiamo imporre che queste approssimazioni soddisfano l’equa-
zione di Laplace
ui+1,j − 2ui,j + ui−1,j
h2+
ui,j+1 − 2ui,j + ui,j−1
k2= 0
(ui+1,j − 2ui,j + ui−1,j)k2 + (ui,j+1 − 2ui,j + ui,j−1)h
2 = 0
h2ui,j−1 + k2ui−1,j − 2(h2 + k2)ui,j + k2ui+1,j + h2ui,j+1 = 0.
Tali relazioni danno luogo ad un insieme di equazioni lineari, una per
ogni punto dell’insieme RN+1,M+1 che non appartiene alla frontiera
Γ.
15
![Page 16: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/16.jpg)
Il valore dell’approssimazione ui,j(x, y), dipende solo dal valore i altri
quattro punti:
ui,j
ui,j+1
ui,j−1
ui+1,jui−1,j
16
![Page 17: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/17.jpg)
Schema a 5 punti per l’equazione di Laplace
I valori della funzione nei vertici delle celle in cui e stato suddiviso
il dominio di integrazione (Rettangolo [a, b] × [c, d]) sono legati dalla
seguente relazione lineare:
h2ui,j−1 + k2ui−1,j − 2(h2 + k2)ui,j + k2ui+1,j + h2ui,j+1 = 0.
per i = 1, . . . , N − 1, e j = 1, . . . , M − 1. Si tratta di (N − 1)(M − 1)
uguaglianze che danno luogo ad un sistema lineare.
17
![Page 18: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/18.jpg)
Lo schema numerico e sintetizzato nel seguente stencil:
−2(h2 + k2) k2k2
h2
h2
18
![Page 19: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/19.jpg)
La prima equazione si ottiene per i = j = 1:
h2u1,0 + k2u0,1 − 2(h2 + k2)u1,1 + k2u2,1 + h2u1,2 = 0
equivalente a
−2(h2 + k2)u1,1 + k2u2,1 + h2u1,2 = −h2u1,0 − k2u0,1.
La seconda equazione si ottiene per i = 2 e j = 1:
h2u2,0 + k2u1,1 − 2(h2 + k2)u2,1 + k2u3,1 + h2u2,2 = 0
equivalente a
k2u1,1 − 2(h2 + k2)u2,1 + k2u3,1 + h2u2,2 = −h2u2,0.
Ogni equazione (i, j) ha al piu 5 coefficienti diversi da 0 di cui 3
coinvolgono 3 incognite numerata consecutivamente (i−1, j), (i, j) e
(i +1, j), una precedente (i, j − 1) e una successiva (i, j +1), distanti
N − 1 incognite (prima e dopo quella di riferimento).
19
![Page 20: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/20.jpg)
Metodi per l’Ordinamento delle incognite
La struttura del sistema lineare dipende dal modo con cui vengono
ordinate le incognite ui,j. Ci sono diversi modi:
1. Ordinamento Lessicografico (o Naturale);
2. Ordinamento Red-Black;
3. Ordinamento Multicolore;
4. Ordinamento Cuthill-McKee;
20
![Page 21: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/21.jpg)
Ordinamento Lessicografico
Si ordinano le incognite partendo dal primo nodo in basso a sinistra e
si procede verso destra, terminata la riga si passa alla successiva riga
in alto.
Nel caso in cui il problema sia di Dirichlet devono essere numerati
solo i punti che non appartengono alla frontiera, poiche negli altri si
suppone che la funzione u(x, y) sia nota.
21
![Page 22: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/22.jpg)
1 2 3 4 5 6 7 8
9 10 11 12 13 14 15 16
17 18 19 20 21 22 23 24
25 26 27 28 29 30 31 32
![Page 23: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/23.jpg)
Ordinamento Red Black
Si dividono le incognite in due colori, a scacchiera, e quindi si procede
con l’ordinamento lessicografico dei due colori in sequenza.
22
![Page 24: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/24.jpg)
1 17 2 18 3 19 4 20
21 5 22 6 23 7 24 8
9 25 10 26 11 27 12 28
29 13 30 14 31 15 32 16
![Page 25: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/25.jpg)
Ordinamento Multicolore
Come l’ordinamento Red-Black con la differenza che i colori sono piu
di due (solitamente quattro o sei).
23
![Page 26: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/26.jpg)
1 9 17 25 2 10 18 26
19 27 3 11 20 28 4 12
5 13 21 29 6 14 22 30
23 31 7 15 24 32 8 16
![Page 27: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/27.jpg)
Ordinamento Cuthill-McKee
Si ordinano le incognite partendo da un nodo arbitrario e numerando
quelli adiacenti che si trovano lungo una direzione prefissata (per
esempio lungo la diagonale).
24
![Page 28: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/28.jpg)
1 2 6 7 14 15 22 23
3 5 8 13 16 21 24 29
4 9 12 17 20 25 28 30
10 11 18 19 26 27 31 32
![Page 29: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/29.jpg)
Struttura delle matrici
• Lo schema a cinque punti definisce un sistema lineare Au = b in cui
la matrice A ha ordine (N − 1)(M − 1);
• I diversi tipi di ordinamento danno luogo a matrici dei coefficienti
che hanno una struttura molto sparsa, cioe la maggior parte degli
elementi sono uguali a zero. Infatti su (N − 1)2(M − 1)2 elementi di
A meno di 5(N − 1)(M − 1) sono diversi da zero.
Consideriamo la costruzione di tali matrici nel caso dell’ordinamento
Lessicografico e Cuthill-McKee.
25
![Page 30: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/30.jpg)
i = j = 1
26
![Page 31: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/31.jpg)
i = 2, j = 1
27
![Page 32: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/32.jpg)
i = 3, j = 1
28
![Page 33: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/33.jpg)
i = N − 1, j = 1
29
![Page 34: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/34.jpg)
i = 1, j = 2
30
![Page 35: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/35.jpg)
i = 2, j = 2
31
![Page 36: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/36.jpg)
i = 3, j = 2
32
![Page 37: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/37.jpg)
i = N − 1, j = M − 1
33
![Page 38: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/38.jpg)
i = j = 1
1 2 6 7 14 15 22 23
3 5 8 13 16 21 24 29
4 9 12 17 20 25 28 30
10 11 18 19 26 27 31 32
34
![Page 39: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/39.jpg)
i = 2, j = 1
1 2 6 7 14 15 22 23
3 5 8 13 16 21 24 29
4 9 12 17 20 25 28 30
10 11 18 19 26 27 31 32
35
![Page 40: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/40.jpg)
i = 1, j = 2
1 2 6 7 14 15 22 23
3 5 8 13 16 21 24 29
4 9 12 17 20 25 28 30
10 11 18 19 26 27 31 32
36
![Page 41: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/41.jpg)
i = 1, j = 3
1 2 6 7 14 15 22 23
3 5 8 13 16 21 24 29
4 9 12 17 20 25 28 30
10 11 18 19 26 27 31 32
37
![Page 42: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/42.jpg)
i = 2, j = 2
1 2 6 7 14 15 22 23
3 5 8 13 16 21 24 29
4 9 12 17 20 25 28 30
10 11 18 19 26 27 31 32
38
![Page 43: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/43.jpg)
i = 3, j = 1
1 2 6 7 14 15 22 23
3 5 8 13 16 21 24 29
4 9 12 17 20 25 28 30
10 11 18 19 26 27 31 32
39
![Page 44: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/44.jpg)
i = 4, j = 1
1 2 6 7 14 15 22 23
3 5 8 13 16 21 24 29
4 9 12 17 20 25 28 30
10 11 18 19 26 27 31 32
40
![Page 45: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/45.jpg)
i = 3, j = 2
1 2 6 7 14 15 22 23
3 5 8 13 16 21 24 29
4 9 12 17 20 25 28 30
10 11 18 19 26 27 31 32
41
![Page 46: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/46.jpg)
i = N − 1, j = M − 1
1 2 6 7 14 15 22 23
3 5 8 13 16 21 24 29
4 9 12 17 20 25 28 30
10 11 18 19 26 27 31 32
42
![Page 47: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/47.jpg)
La struttura della matrice che si ottiene applicando il metodo a 5
punti all’equazione di Laplace dipende da:
• La forma del dominio
• Il tipo di ordinamento utilizzato.
43
![Page 48: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/48.jpg)
La struttura della matrice influisce anche sull’efficienza del metodo
che si deve utilizzare per risolvere il sistema lineare.
Infatti la matrice A ha dimensioni molto elevate ma il numero di
elementi diversi da zero e piuttosto basso.
Se si utilizza un algoritmo di fattorizzazione per risolvere il sistema
(Fattorizzazione LU , di Choleski) si sarebbe auspicabile che anche i
fattori abbiano un numero di elementi diversi da zero che sia il piu
basso possibile.
44
![Page 49: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/49.jpg)
0 5 10 15 20 25 30
0
5
10
15
20
25
30
nz = 136
Ordinamento Lessicografico
45
![Page 50: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/50.jpg)
0 5 10 15 20 25 30
0
5
10
15
20
25
30
nz = 231
Ordinamento Lessicografico−Fattore U
46
![Page 51: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/51.jpg)
0 5 10 15 20 25 30
0
5
10
15
20
25
30
nz = 136
Ordinamento Red−Black
47
![Page 52: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/52.jpg)
0 5 10 15 20 25 30
0
5
10
15
20
25
30
nz = 170
Ordinamento Red−Black−Fattore U
48
![Page 53: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/53.jpg)
0 5 10 15 20 25 30
0
5
10
15
20
25
30
nz = 136
Ordinamento Multicolor
49
![Page 54: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/54.jpg)
0 5 10 15 20 25 30
0
5
10
15
20
25
30
nz = 191
Ordinamento Multicolor−Fattore U
50
![Page 55: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/55.jpg)
0 5 10 15 20 25 30
0
5
10
15
20
25
30
nz = 136
Ordinamento Cuthill−McKee
51
![Page 56: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/56.jpg)
0 5 10 15 20 25 30
0
5
10
15
20
25
30
nz = 148
Ordinamento Cuthill−McKee−Fattore U
52
![Page 57: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/57.jpg)
Il metodo a 5 punti per l’equazione di Laplace su un rettangolo
• Il dominio continuo viene discretizzato e le derivate parziali seconde
sono approssimate negli (N − 1)(M − 1) punti del dominio discreto;
• Applicando la tecnica lessicografica di ordinamento delle incognite
si deve risolvere il sistema lineare
Au = b
dove la matrice A ha la seguente struttura tridiagonale a blocchi
A =
T J
J T J. . . . . . . . .
J T J
J T
53
![Page 58: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/58.jpg)
dove J = h2IN−1, ed IN−1 e la matrice identita di ordine N − 1, e T
e la seguente matrice tridiagonale di dimensione N − 1:
T =
−2(h2 + k2) k2
k2 −2(h2 + k2) k2
. . . . . . . . .
k2 −2(h2 + k2) k2
k2 −2(h2 + k2)
.
![Page 59: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/59.jpg)
• Il sistema da risolvere e di grandi dimensioni e ha struttura sparsa,
cioe solo pochi elementi sono diversi da zero.
• Se il dominio e quadrato oppure si puo scegliere h = k in questo
caso l’espressione del metodo a 5 punti si semplifica ulteriormente:
ui,j−1 + ui−1,j − 4ui,j + ui+1,j + ui,j+1 = 0
e la struttura della matrice A e indipendente da h e k:
A =
T IN−1IN−1 T IN−1
. . . . . . . . .
IN−1 T IN−1IN−1 T
54
![Page 60: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/60.jpg)
mentre T e la seguente matrice tridiagonale di dimensione N − 1:
T =
−4 11 −4 1
.. . . . . . . .
1 −4 11 −4
.
• Per risolvere questo tipo di sistemi sono necessari metodi diversi da
quelli di tipo diretto (occorrerebbe gestire una quantita di memoria
molto elevata).
![Page 61: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/61.jpg)
Metodi di ordine elevato per equazioni ellittiche
• Consideriamo per semplicita l’equazione ellittica di Poisson
∂2u
∂x2+
∂2u
∂y2= g(x, y)
definita sul quadrato Ω = [0,1] × [0,1].
• Il metodo a 5 punti ovviamente puo essere applicato anche al-
l’equazione di Poisson. La matrice dei coefficienti e la stessa che
abbiamo visto nel caso dell’equazione di Laplace, cambia il termine
noto del sistema lineare che e uguale ai valori della funzione g(x, y)
calcolata nei punti del dominio discreto.
55
![Page 62: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/62.jpg)
• Avendo definito l’equazione su un quadrato possiamo scegliere un
numero di punti uguale sia per la variabile x che per y (e quindi h = k),
ottenendo una matrice dei coefficienti indipendente dal valore del pas-
so di discretizzazione.
• Il metodo a 5 punti ha ordine 2 perche l’errore di discretizzazione
dipende da h2.
• Si possono ottenere metodi di ordine superiore cercando di approssi-
mare in modo piu accurato le derivate parziali seconde.
56
![Page 63: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/63.jpg)
Un’approssimazione di ordine 4 per le derivate seconde
Se f(x) e una funzione di classe Cvi([a, b]) allora la seguente formula
approssima la derivata seconda nel punto xi:
f ′′(xi) ≃1
h2
[
−1
12fi−2 +
4
3fi−1 −
5
2fi +
4
3fi+1 −
1
12fi+2
]
con un errore
E(f ′′(xi)) = ch4f(vi)(ξi), c ∈ R,
e pertanto ha ordine 4.
57
![Page 64: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/64.jpg)
Utilizzare questa formula per le derivate parziali seconde porta al
seguente metodo
−1
12ui,j−2 +
4
3ui,j−1 −
1
12ui−2,j +
4
3ui−1,j − 5ui,j +
4
3ui+1,j
−1
12ui+2,j +
4
3ui,j+1 −
1
12ui,j+2 = h2g(xi, yj)
58
![Page 65: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/65.jpg)
Il metodo puo essere descritto dal seguente stencil:
−5
43
43
43
43
− 112
− 112
− 112
− 112 uij = h2g(xi, yj)
59
![Page 66: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/66.jpg)
Osservazioni
• Il metodo e piu preciso del metodo a 5 punti ma ci sono problemi
quando viene applicato in prossimita della frontiera come si evince
dalla seguente immagine.
60
![Page 67: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/67.jpg)
i = j = 1
61
![Page 68: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/68.jpg)
• Il sistema lineare
Au = b
che si deve risolvere ha una struttura piu complessa di quello ottenuto
con il metodo a 5 punti.
62
![Page 69: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/69.jpg)
La matrice dei coefficienti e la seguente
A =
T0 J0 O
J1 T J1 J2J2 J1 T J1 J2
. . . . . . . . . . . . . . .
J2 J1 T J1 J2J2 J1 T J1
O J0 T0
dove, utilizzando l’ordinamento lessicografico, la matrice T ha una
struttura a 5 diagonali (pentadiagonale), la matrice T0 e tridiagonale,
mentre le matrici J0, J1 e J2 sono matrici diagonali.
63
![Page 70: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/70.jpg)
Il Metodo a 9 punti
Utilizzando una diversa approssimazione delle derivate parziali seconde
della funzione u(x, y) si puo derivare uno schema numerico che utilizza
nove punti collocati intorno al nodo di riferimento. Senza descrivere
in dettaglio tale tecnica, lo stencil del metodo applicato con lo stesso
passo di discretizzazione rispetto alle variabili x e y, e il seguente:
64
![Page 71: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/71.jpg)
-103
23
23
23
23
16
16
16
16
uij = h2g(xi, yj)
65
![Page 72: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/72.jpg)
Osservazioni
• Il metodo e preciso quanto il metodo a 5 punti (l’errore dipende
da h2) ma l’uso di maggiori informazioni (l’approssimazione coinvolge
appunto 9 punti) consente di ottenere in pratica risultati migliori;
• Il sistema lineare
Au = b
che si deve risolvere ha una struttura piu complessa di quello ottenuto
con il metodo a 5 punti.
66
![Page 73: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/73.jpg)
La matrice dei coefficienti e la seguente
A =
T J
J T J
J T J. . . . . . . . .
J T J
J T J
J T
dove, utilizzando l’ordinamento lessicografico, le matrici T e J sono
entrambe tridiagonali.
67
![Page 74: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/74.jpg)
L’equazione di Laplace su domini non rettangolari
Il metodo a 5 punti funziona bene quando il dominio Ω e un rettan-
golo, oppure un quadrato o un poligono che puo essere scomposto
come un unione di quadrati o rettangoli.
68
![Page 75: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/75.jpg)
69
![Page 76: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/76.jpg)
70
![Page 77: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/77.jpg)
Un altro caso favorevole e quello in cui il dominio Ω coincide con un
cerchio e Γ e una circonferenza, per esempio:
Ω =
(x, y) ∈ R2 : x2 + y2 < 1
e
Γ =
(x, y) ∈ R2 : x2 + y2 = 1
.
In questo caso il dominio puo essere trasformato in un rettangolo
cambiando le coordinate cartesiane in polari:
Ω =
(ρ, θ) ∈ R2 : 0 ≤ ρ < 1, 0 ≤ θ < 2π
.
e
Γ =
(ρ, θ) ∈ R2 : ρ = 1, 0 ≤ θ < 2π
.
71
![Page 78: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/78.jpg)
ρx
θy72
![Page 79: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/79.jpg)
I problemi che sorgono sono di due tipi:
• Scrivere l’equazione di Laplace in coordinate polari:
∂2v
∂ρ2+
1
ρ2
∂2v
∂θ2+
1
ρ
∂v
∂ρ= 0.
• Condizioni iniziali:
v(ρ,0) = u(ρ,2π), v(0, θ) = costante, v(1, θ) = f(cos θ, sin θ).
dove v(ρ, θ) = u(ρ cos θ, ρ sin θ).
73
![Page 80: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/80.jpg)
Risoluzione dell’Equazione di Laplace su domini irregolari
Analizziamo ora il caso in cui la frontiera del dominio Ω sia una
curva chiusa e regolare senza una particolare forma. In questo caso
si considera un rettangolo [a, b] × [c, d] tale da contenere sia Ω che Γ
e si discretizza tale regione come gia visto in precedenza, definendo
l’insieme discreto:
RN+1,M+1 =
(xi, yj) ∈ [a, b] × [c, d]|i = 0, . . . , N, j = 0, . . . , M
.
74
![Page 81: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/81.jpg)
y
x
Ω
Γ
O
75
![Page 82: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/82.jpg)
y
xa b
c
d
Ω
Γ
O
76
![Page 83: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/83.jpg)
y
xa b
c
d
Ω
Γ
O
77
![Page 84: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/84.jpg)
L’insieme dei punti discreti che appartengono sia al rettangolo [a, b]×
[c, d] che al dominio Ω si denota con
Ωhk = RN+1,M+1 ∩ Ω.
Tale insieme viene detto insieme dei punti interni. Ogni punto interno
(xi, yj) ha quattro punti vicini, cioe (xi±1, yj) e (xi, yj±1).
Un punto vicino ad un punto interno che non appartiene a Ωhk si dice
punto di confine. L’insieme dei punti di confine si indica con Γhk.
I punti angolari sono invece i punti dell’insieme RN+1,M+1 che non
sono ne interni ne di confine, ma risultano essere vertici di una cella
che contiene almeno un punto interno.
78
![Page 85: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/85.jpg)
rs
rsbc
bc
bc
rs
bc
bc
ut
ut
ut
rs
ut
ut
rs
Γ
Ω
Punti di confiners
Punti angolariut
Punti internibc
79
![Page 86: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/86.jpg)
I punti angolari non hanno alcun ruolo nella soluzione numerica di
equazioni ellittiche. La griglia deve essere tale che i segmenti che
congiungono i punti interni devono essere interamente contenuti nel
dominio Ω. Si evita cioe il caso evidenziato dal secondo grafico:
80
![Page 87: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/87.jpg)
(xi+1, yj)(xi, yj)(xi, yj) (xi+1, yj)
81
![Page 88: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/88.jpg)
Di solito tale situazione puo essere evitata scegliendo opportunamente
il passo di discretizzazione oppure effettuando un opportuno cambio
di variabile (per esempio una rotazione degli assi).
Il problema dei domini irregolari sorge quando nelle approssimazioni
delle derivate seconde intervengono valori della funzione calcolati nei
punti di confine. Per ovviare a tale inconveniente ci sono diverse
tecniche, due sono le piu usate:
1. come valori nei punti di confine si utilizzano gli stessi valori assunti
dalla condizione al contorno;
2. si utilizza il valore della condizione al contorno nel punto sulla
curva Γ piu vicino al punto interno.
82
![Page 89: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/89.jpg)
bc bc
bcΓ
(xi, yj)
(xi, yj+1)
(xi+1, yj)
si pone:
u(xi, yj+1) = f(xi, yj+1), u(xi+1, yj) = f(xi+1, yj).
83
![Page 90: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/90.jpg)
Tale assegnazione rappresenta, dal punto di vista matematico, una
forzatura, poiche in realta non e noto neanche se la funzione f(x, y)
sia calcolabile in tali punti.
84
![Page 91: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/91.jpg)
bc bcbc
Γ
xi − h xi xi + h
85
![Page 92: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/92.jpg)
bc bcbc
Γ
xi − h xi xi + hxi + sh
86
![Page 93: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/93.jpg)
Approssimazione della derivata seconda su griglie non uniformi
Consideriamo il problema di approssimare la derivata seconda della
funzione f(x) nel punto di ascissa xi ma considerando i valori della
funzione in punti non equidistanti xi − h e xi + sh, con 0 < s < 1.
Sviluppiamo in serie di Taylor la funzione f(xi + sh) prendendo come
punto iniziale xi:
f(xi+sh) = f(xi)+shf ′(xi)+s2h2
2f ′′(xi)+
s3h3
6f ′′′(ξi), ξi ∈ [xi, xi+sh]
e procediamo in modo analogo per f(xi−1):
f(xi−1) = f(xi) − hf ′(xi) +h2
2f ′′(xi) −
h3
6f ′′′(ηi), ηi ∈ [xi−1, xi].
87
![Page 94: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/94.jpg)
Moltiplichiamo per s quest’ultima espressione
sf(xi−1) = sf(xi) − shf ′(xi) +sh2
2f ′′(xi) −
sh3
6f ′′′(ηi)
e sommiamola con quella di f(xi + sh):
f(xi + sh) + sf(xi−1) = f(xi)(1 + s) +h2
2f ′′(xi)s(1 + s)
+sh3
6
[
s2f ′′′(ξi) − f ′′′(ηi)]
88
![Page 95: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/95.jpg)
Si ricava l’espressione della derivata seconda in xi:
f ′′(xi) = 2f(xi + sh) − f(xi)(1 + s) + sf(xi−1)
sh2(1 + s)+
+h
3(1 + s)
[
f ′′′(ηi) − s2f ′′′(ξi)]
e, trascurando l’ultimo termine, l’approssimazione della derivata se-
conda e:
f ′′(xi) ≃ 2f(xi + sh) − f(xi)(1 + s) + sf(xi−1)
sh2(1 + s)
mentre l’errore vale:
E(f ′′(xi)) =h
3(1 + s)
[
f ′′′(ηi) − s2f ′′′(ξi)]
.
89
![Page 96: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/96.jpg)
Applicando tale risultato ad una funzione in due variabili si otterrebbe
l’approssimazione:
uxx(xi + sh, yj) ≃ 2u(xi + sh, yj) − (1 + s)u(xi, yj) + su(xi−1, yj)
sh2(1 + s).
Ovviamente si procederebbe in modo analogo per approssimare la
derivata seconda uyy(xi, yj + tk), 0 < t < 1:
uyy(xi, yj + tk) ≃ 2u(xi, yj + tk) − (1 + t)u(xi, yj) + tu(xi, yj−1)
tk2(1 + t).
90
![Page 97: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/97.jpg)
Osservazioni
• Per s = 1 (o t = 1) la formula e la stessa gia vista in precedenza;
• L’approssimazione vista e del primo ordine perche l’errore dipende
da h (o da k);
• L’approssimazione della soluzione sara meno precisa sulla frontiera
di Ω;
• La struttura della matrice dei coefficienti non cambia perche e
coinvolto lo stesso numero di punti sia nel caso dei punti interni
che di frontiera.
91
![Page 98: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/98.jpg)
Volendo ottenere un’altra approssimazione si puo utilizzare la seguente
formula, della quale omettiamo la dimostrazione, che utilizza l’ulte-
riore valore f(xi−2) ma e di ordine 2:
f ′′(xi) ≃1
h2
[
s − 1
s + 2f(xi−2) +
2(2 − s)
s + 1f(xi−1) −
3 − s
sf(xi)+
+6f(xi + sh)
s(s + 1)(s + 2)
]
.
92
![Page 99: Metodi Numerici per Equazioni Ellittichetiziano19661.interfree.it/pdf1011/Ellittiche.pdfMetodi Numerici per Equazioni Ellittiche Vediamo ora di descrivere una tecnica per la risoluzione](https://reader035.fdocument.pub/reader035/viewer/2022062610/61126c1aea26ae4e233e09a2/html5/thumbnails/99.jpg)
Osservazioni
• L’approssimazione vista e del secondo ordine perche l’errore dipende
da h2 (o da k2);
• L’approssimazione della soluzione sara precisa allo stesso modo
su tutti i punti del dominio;
• La struttura della matrice dei coefficienti cambia in prossimita dei
punti sulla frontiera perche viene utilizzato un punto in piu.
93