15/12/2025
Metodi HBUM
- IsEsR
HBVM (KS) e un metodo R-K a 10 stadi
, definito dal tableau di Butcher
dove& è il verore delle ascisse bi
E è il vettore dei pesi
Pa(cil = 0 i = 1
, -c
k .
P()] Spax I
-
11
be
= Po(G)
I
1 =
Ps Es =
Sk
Po([k) Po((r)
:
(0) Is = Est2Ys dove
1
Si =
ixo
,
2 (di 21
.
(11 ... 1)EIRK
Abbiamo visto che l'implementazione classica come
↑
S ij exyo + hIsPRImf1
HBVM classico applicata al pb .
=
f(y) ,
telo , h)
significherebbe risolvere =
y(0) yoIR" I
(I ei usaih)
- = dove = =
con u() =
yo +
hj(x)dxy
[
-
= =
exyo + hi
+
1) Im =
Y =
20yo +
hIsIm
=>
8 = IsImf(exyo + hIsQImG) che ha dim = s
,
j =
Es2 * Im f(x) VkYS .
Chiamiamo G(*) =
j -
&SImf(exyothIsQIm) = e irsm il
L'applicazione di un metodo di Newton semplificato prevede che si risolva un sistema lin . Che
ha come matrice dei coeff .
la J(G) con hp
. che 5(f) sia 5(f)(y0)
↑ punto iniziale .
·
calcoliamol derivata di (1)
* Im-Isf'(yd)
.
hs QIm
=>
IQIm hIs'RIs * f'(yd
↓
(0)
ISM Is * Ist Is + s
Ys
1
ElpSXS
+
>SRIs +1
i = 1
j -
, Sj = 1
,
-
,
S+ 1
(b2
Pj z(()
(ISIstlij (pin(C Pic(r))
-
= ...
=
bk
Pj 2([k)-
= be pie(e)pj = (e) (2)
che è una quadratura di ordine 24 325 .
L'integrando ha ordine (i -1) 1) i+
j
+
(j 22S + S+1 2 25 1
- = - - = -
=> la quadrature è esalta
Per cui (2) =
[Dic (* Pj (x)dx
- =
Sij
Allora &SREst =
[[12] quindi ESTREstaYs = Xs
[3)
Risolvo[[QIm-hXsef'(yd]1j8 = -
G(jl)
aggiorno +1 A , E
e
... scegliendo
la matrice
converge velocemente ma richiede di fattorizzateche (3) è una matrice smysm
= che ha un costo
di (sm)
ci ricondurremo a
Em , indipendente quindi da .
s
Supponiamo che f'lyo) sia diagonalizzabile e sia 1 un
generico autovalore .
Se indichiamo al solito hX allora possiamo decomporre il problema
, , q =
, ,
nel sottospazio di competenza, in uno più semplice della forma :
dipenda dalla sol. Stessa
(4)(z 9(s)x -
=
n
↓
non sing Vs
.
(4) e
equivalente a
=
SX5(I-qXs)x ,
dove go e uno scalare da determine
(5)g(X5 -
g])x =
Ma
( seats
seg = 0
Definiamo la funzione peso O(q) = I .
(2-pq =
Questo ci permette di ottenere un interiore pb .
equivalente :
O(q((4) + (1 -
O(a))(5)
[0(q)(I -qXs) + (I -
O(q))g(Xs" 9])]x -
=
0(q)2 + (I -
a(q)(41
!
M(q) =
(seo e aso
2
M(g) N(q) [(1 gq) 0(q)
-
= -
=
Questo induceto
splitting N(q) X =
(N(q) -
M(q) x +
ocqy + (I -
0(q))12
nz + oc(n -
22)
ym (I y2)]
0()[nz
*
= > -
O(q)M(q))x + +
0(g)(y -
Metodo iterativo che
convergerà quando il raggio spett .
di 11-0(q)(M(g) e [1 .
f si
sceglie in modo da minimizz il .
max valore del
raggio Spett. .
FqE . Si arriva dire che
f = min
Il .
MET(Xs)
dipende dal metodo, non dal problema !
Oss che possiamo fare un'ulteriore
.
semplificazione .
se partiamo da Xo = 0
,
non faccio il primo pezzo X" 0
=
L
x
+
=
( acq)[nz"
-
a(q)(x + + o(q)(r ni)) -
aggiorno ne 12
ad
ogn
s
= xr + 1 =
O(q)[n2(xx) +
O(q)(R(XY) -
R2(Xr))] = 0
,
1 . . .
definisce l'iterazione " blended" per risolvere la (1) .
per tornare al metodo di Newton semplificato ,
l'iterazione si modifica come segue
:
·
calcolo 0 =
(Im-hefcya) , 8
:
=
0
·
per l = 0
,
1
, ... poniamo ne =
-G(jb)(e(1))
n? =
(9Xs[m)pe
faccio prod . tensoriale
uscire di dim.
perché
deve .
= ms
gl G (10)[22 [x0(22 22)] [ 2. JeRX
+1
=
+ +
·
-
I =
Is =
Il costo principale qui è la fattorizzazione di 0 che non dipende da s
Funzione Vec : veccal =
=
es T =
[gjz ...
j5 2)
-
= vec(t) =
·
vec[AXB) =
(BA) ·
Vec [])
=
G(f) = T -
FM&s dove F =
[f(y2) ...
f)[m)] emx
se = [Is . . .
Yr] = vec(]) =
Alla luce di ciò, l'iterazione sarà :
2
IRmXS
/Im-hgf'(yd)
-
·
calcolo 0 e No & E
pongo
= =
·
per l = 0
.
1
, ...
poniamo 48 = -( ↑ e)
pe(a
Ne
+ 2
una e se
·
+b + 1 = +1 + [42 + (42 -
4)8]85