INTERPOLAZIONE Quali condizioni si possono richiedere sulla funzione interpolante ? Quali condizioni si possono richiedere sulla funzione interpolante ? Problema generale di INTERPOLAZIONE Dati n punti distinti (x i ,y i ) i=1,..,n si vuole costruire una funzione f(x) tale che nei nodi (x i ) i=1,..n soddisfi a certe condizioni, dette Condizioni di interpolazione
35
Embed
interpolazione - Consorzio Nettunonettuno.unina.it/files/calcolo/interpolazione_liv0.pdf · Interpolazione di Lagrange: formulazione generale Sia F uno spazio di funzioni di variabile
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
1
INTERPOLAZIONE
Quali condizioni si possono richiedere sulla funzione interpolante ?
Quali condizioni si possono richiedere sulla funzione interpolante ?
Problema generale di INTERPOLAZIONE
Dati n punti distinti (xi,yi)i=1,..,n
si vuole costruire una funzione f(x)
tale che nei nodi (xi)i=1,..n
soddisfi a certe condizioni, dette
Condizioni di interpolazione
2
Esempio 1
condizioni di interpolazione:la funzione deve passare per i punti assegnati
in ciascun nodo xi è assegnata una condizione sulla funzione interpolante
x1 x2 xn…y1 y2 y2
Interpolazione di Lagrange
3
Interpolazione di Lagrange: formulazione generale
Sia F uno spazio di funzioni di variabile reale o complessa.
Assegnati:• n valori distinti reali o complessi {xi} i=1,..n
• n valori distinti reali o complessi {yi} i=1,..n
Si cerca una funzione f(x)∈F tale che:
f(xi)=yi i=1,…,n
Condizioni di interpolazione di Lagrange
Esempio 2
condizioni di interpolazione:la funzione deve passare per i punti assegnati
con una fissata tangente
4
N.B.….se è assegnata in un nodo la condizione sulla derivata di ordine q, deve essere assegnata anche la condizione sulle
derivate dall’ordine 0 fino a q-1.
y1ny1
n-1…y12y1
1
y0ny0
n-1…y02y0
1
y32
xn
y2n-1
xn-1
…y22
x2x1
In ciascun nodo sono assegnate condizionisulla funzione f(x) e sulle sue derivate….
Interpolazione di Hermite
Interpolazione di Hermite: Formulazione generale
Sia F uno spazio di funzioni di variabile reale o complessa.
Si cerca una funzione f ∈F tale che:
f(j)(xi)=yji i=1,…,n j=0,1,…,l i-1
Condizioni di interpolazione di Hermite
n numeri interi {li} i=1,..n tali che Σ ni=1 li = m+1
n valori distinti reali o complessi {xi}i=1,..n
m+1 valori distinti reali o complessi {yji} i=1,…,n
j=0,1,…,li-1
Assegnati
5
In ciascun nodo c’è almeno una condizionesulla funzione o sulla derivata di
qualunque ordine
y1n…y1
1
y0ny0
n-1…y02y0
1
y32
xn
y2n-1
xn-1
…y22
x2x1
Interpolazione di Hermite-Birkoff
….quale forma può avere una
funzione interpolante
?
….quale forma può avere una
funzione interpolante
?
OVVERO:
quale spazio di funzioni Futilizzare
?
OVVERO:
quale spazio di funzioni Futilizzare
?
6
Esempio
-0.97
3
-1.58
π/2
-0.39
π/3
-0.31
-1
1.08
-π/2
y 0.41
-2x
-0.97
3
-1.58
π/2
-0.39
π/3
-0.31
-1
1.08
-π/2
y 0.41
-2x
Interpolazione mediante polinomio di quinto grado
7
-0.97
3
-1.58
π/2
-0.39
π/3
-0.31
-1
1.08
-π/2
y 0.41
-2x
Interpolazione mediante polinomio trigonometrico di quinto grado
un confronto ……
Quale funzione scegliere ?
8
• dal costo computazionale relativo alla costruzione ed alla valutazione della funzione interpolante
la scelta della forma della funzione interpolante ovvero
la scelta dello spazio F dipende :
• dalle caratteristiche del problema reale
In generale,
….fissiamo
lo spazio dei polinomi ….fissiamo
lo spazio dei polinomi
…. Costruiamo il polinomio interpolante di Lagrange
…. Costruiamo il polinomio interpolante di Lagrange
….è unico ?
9
Esempio: Siano P1, P2 due punti fissati.
Esistono infinite parabole per P1, P2
Ma una sola retta per P1, P2
dati i due punti che grado deve avere il polinomio affinché la curva
interpolante sia univocamente determinata?
dati i due punti che grado deve avere il polinomio affinché la curva
interpolante sia univocamente determinata?
10
2) Per m =1 (in generale m=n-1 )
Esiste un unico polinomio di grado 1 interpolante i due punti( esiste un’unica retta passante per due punti )
1) Per m < 1 (in generale m< n-1 )
Non esiste alcun polinomio interpolante ( tranne se i due punti hanno la stessa ordinata )
3) Per m > 1 (in generale m >n-1 )
Esistono infiniti polinomi di grado maggiore di 1 interpolanti i due punti
( esistono infinite parabole, cubiche , etc. per due punti )
Quindi...
Teorema
Dati n nodi distinti {xi} i=1,..nn valori corrispondenti {yi} i=1,..n
il polinomio p(x)∈Πm tale che:
p(xi) = yi i=1,…,n
è unico se:m = n-1
11
Siano p(x), q(x) ∈Πn-1 tali che:
p(xi) = yi i=1,..nq(xi) = yi i=1,..n
Dimostriamo che i due polinomi sono identicamente uguali. Posto:
z(x) = p(x)- q(x) (z(x) ∈Πn-1 )
Si ha:z(xi) = 0 i=1,..n
1/2
Dimostrazione 1
1)Se z(xi) = 0 i=1,..n ⇒ il polinomio ha n zeri distinti⇓
z(x) =a(x- x1) (x- x2)… (x- xn) a∈ℜ
2) Un polinomio non nullo di grado al più n-1 ha al più n-1 zeri distinti
(Teorema fondamentale dell’algebra)
3) L’unico polinomio al più di grado n-1 che ha n zeri distinti è il polinomio identicamente nullo
z(x) è il polinomio identicamente nullocioè p(x)=q(x)
2/2
12
Dimostrazione 2
Se p(x) ∈Πn-1 è il polinomio interpolante gli n punti:(xi, yi ) i=1,..n
dovendo soddisfare le condizioni di interpolazionep(xi) = yi i=1,..n
Dimostriamo che il sistema ammette una ed una sola soluzione
Esistenza ed unicità della soluzione del sistema lineare
3/4
4/4
Poiché i nodi sono a due a due distinti alloradet(V)≠0
Esiste un’unica soluzione
c.v.d
Il sistema è non singolare
14
….quali condizioni garantiscono
l’unicità del polinomio
interpolante di Hermite?
….quali condizioni garantiscono
l’unicità del polinomio
interpolante di Hermite?
ESEMPIO 6 Traiettoria di una pallottola
Una pallottola è sparata vero l’alto dalla sommità di un edificio alto 100 m. Dopo 10 sec. dal lancio la pallottola raggiunge la
massima altezza a 590 . Disegnare la traiettoria della pallottola.
QQ11(0,100)(0,100)
QQ22(10,590)(10,590)
15
Come costruire la parabola che descrive la
traiettoria della pallottola ?
Come costruire la parabola che descrive la
traiettoria della pallottola ?
L’espressione della parabola è p(x) = axL’espressione della parabola è p(x) = ax2 2 + + bxbx + c+ c
Q1 (0,100)
•Q2 (10,590)
Sappiamo che p(0)=100Sappiamo che p(0)=100p(10)=590p(10)=590
Sappiamo inoltre che il punto Q2 è il vertice della parabola
Possono esistere Possono esistere Più parabole…Più parabole…
Esiste una sola Esiste una sola parabola…parabola…
I condizione
II condizione
III condizione
16
p(xp(x11)=y)=y11
p(xp(x22)=y)=y22
p’(xp’(x22)=0)=0
Si ottiene il sistema …Si ottiene il sistema …
2ax2ax22+b=0+b=0
axax222 2 + bx+ bx22 + c=y+ c=y22
axax112 2 + bx+ bx11 + c=y+ c=y11
c=100100 a + 10 b + 100=59020a+b=0
b=88b=88c=100c=100
a=4.4a=4.4
UNICA SOLUZIONE
Interpolazione di Hermite
n=2 nodi ( 2 punti per cui deve passare la parabola)
m=3 nodi condizioni di interpolazione(appartenenza dei due punti alla parabola + vertice nel primo)
# condizioni =
# incognite
17
TEOREMA
Dati n nodi distinti {xi} i=1,..nn numeri interi positivi {li} i=1,..n tali che
Σ ni=1 li = m+1
ed m+1 valori {yji} con i=1,…,n j=0,1,…,li-1
Il polinomio p (x) ∈ Πq tale che:
p(j)(xi)=yji
è unico se:q = m
Siano p(x), q(x) ∈Πm tali che:
p(j) (xi) = yji i=1,..n j=0,1,…,li-1
q(j) (xi) = yji i=1,..n j=0,1,…,li-1
dimostriamo che i due polinomi sono identicamente uguali. Posto:
z(x) = p(x)- q(x) (z(x) ∈Πm )
dimostriamo che z(x) è un polinomio identicamente nullo. Si ha:
z(j) (xi) = 0
1/2
Dimostrazione
18
Un polinomio z∈Πm tali che:
z(j) (xi) = yji i=1,..n j=0,1,…,li-1
È necessariamente del tipo:
z(x) =a(x- x1)l1 (x- x2) l2 … (x- xn) ln a∈ℜ
Cioè ammette n zeri distinti x1 x2 … xn rispettivamente con molteplicità :
l1 , l2 ,… ,lnPoiché Σ n
i=1 li=m+1 e siccome z(x) ha la più grado m, z(x) è necessariamente il polinomio identicamente nullo.
c.v.d2/2
In generale…
Se poniamo F = Πmspazio dei polinomi di
grado al più m
Interpolazione polinomiale
di Lagrange di Hermite
19
Sia per l’interpolazione di Lagrange
che per l’interpolazione di Hermite
# condizioni =
# incognite
l’unicità
della soluzione è garantita
Metodi costruttivi del polinomio interpolante di
di LAGRANGE
20
p(x) ∈Πn-1 è p(x)= a0 + a1x +…+ an-1xn-1
METODO DEI COEFFICIENTI INDETERMINATI
Imponendo le condizioni di interpolazionep(xi) = yi i=1,..n
a0 + a1 x1 +…+ an-1 x1 n-1 =y1
a0 + a1 x2 +…+ an-1 x2 n-1 =y2
a0 + a1 xn +…+ an-1 xnn-1 =yn
……………………………
Sistema di equazioni lineari di ordine n
p(x1) = y1
p(x2) = y2
p(xn) = yn
…..…..
La matrice dei coefficienti di questo sistema è:
=
−
−
−
1nnn
1n22
1n11
x...x1............
x...x1x...x1
V ∏>=
−=n
ji1j
ji )x(xdet(V)
Matrice di Vandermonde( mal condizionata!)
21
Formula di Lagrange
Esempio: Siano Esempio: Siano P1(x1,y1), P2(x2,y2) punti del piano.punti del piano.Costruiamo il polinomio Costruiamo il polinomio p(x) interpolante interpolante P1, P2 ..
Poniamo Poniamo p(x) =y1 l1 (x) + y2 l2(x) con lcon l11(x), l(x), l22(x)(x)ppolinomi olinomi ddi primo grado.i primo grado.
Di che tipo devono essere i polinomi lDi che tipo devono essere i polinomi l11(x) e l(x) e l22(x)?(x)?
ll11(x)= a (x(x)= a (x--xx22))ll22(x)= b (x(x)= b (x--xx11))
ll11(x(x11)= a (x)= a (x11--xx22) =1) =1ll22(x(x22)= b (x)= b (x22--xx11)=1)=1