TD d’analyse numérique 3 SMA S6

Prof. : Jaouad BENNOUNA

Université Sidi Mohamed Ben Abdelah, Fès

Faculté des Sciences Dhar El Mahraz, Fès

1
 Table des matières

Série N˚1

Exercice 1.
Soit H un espace de Hilbert. K ⊂ H un convexe fermé non vide.
1) montrer que ∀f ∈ H il existe u ∈ K unique tel que :

(1) |f − u| = min|f − v|
v∈K

On notera par |u| la norme associé au produit scalaire (u, u) de l’espace H.

2) Montrer que u est caractérisé par la propriété :

 u ∈ K,
(2)
 (f − u, v − u) ≤ 0 ∀v ∈ K.

On notera que pour la suite u = pK f = projection de f sur K.

3) Montrer que |pK f1 − pK f2 | ≤ |f1 − f2 | ∀f1 , f2 ∈ H
4) Soit M ⊂ H un s.e.v fermé. Soit f ∈ H
Montrer que u = pM f est caractérisé par :

 u∈M
 (f − u, v) = 0 ∀v ∈ M.

Exercice 2.
On considère N points a = x1 < x2 < ... < xN = b de [a, b] et N réels βi , i = 1, ..., N
Déterminer la fonction splines d’interpolation d’ordre 2 tel que :

σ(xi ) = βi 1≤i≤N et σ 00 (a) = σ 00 (b) = 0

Exercice 3.
Soit X = {x1 , ..., xN } ⊂ Ω borné de RN (xi distincts)

1
2 TABLE DES MATIÈRES

Soit m ∈ N∗ . On suppose que X est Pm−1 unisolvant c’est-à-dire ;

∀αj ∈ R, 1 ≤ j ≤ N , il existe un polynôme p de Pm−1 tel que

P (xj ) = αj 1≤j≤N et que m > n/2

Soit β = (β1 , .., βN )

On introduit ρ : H m (Ω) → RN tel que pv = (v(x1 ), v(x2 ), ..., v(xN ))

Kβ = K = {v ∈ H m (Ω) : ρv = β}
K0 = {v ∈ H m (Ω) : ρv = 0}.

On définit sur H m (Ω) le produit scalaire et la norme

m Z b
X 1/2
((u, v))m,Ω = u(i) (t)v (i) (t)dt ||u||m,Ω = ((u, u))m,Ω
i=0 a
le semi produit
Z b scalaire et la semi-norme
1/2
(u, v)m,Ω = u(m) (t)v m (t)dt |u|m,Ω = (u, u)m,Ω
a
N
X
N N
On définit sur R : ∀ξ, η ∈ R < ξ, η >= ξi ηi < ξ >=< ξ, ξ >1/2
i=1

On définit sur H m (Ω) un autre produit scalaire et une autre norme :

(((u, v))) =< ρu, ρv > +(u, v)m,Ω |||u||| = (((u, u)))1/2 = (< ρu >2 +|u|2 )1/2
H m (Ω) muni de la norme ((( , ))) est Hilbertien.
On pose S = {u ∈ H m (Ω) ∀v ∈ K0 , (u, v)m,Ω = 0}.
On admet que la norme ||.||m,Ω est équivalente à |||.|||.
1) montrer que l’opération ρ est linéaire continue ρ ∈ L(H m (Ω), RN ).
Indicateur : H m (Ω) ,→ C 0 (Ω̄) .

2) Montrer que K est fermé convexe non vide de H m (Ω).

3) Montrer que K0 est un s.e.v fermé de H m (Ω).
4) On considère le problème :

 T rouver σ∈K
(1)
 ∀v ∈ K |σ|m,Ω ≤ |v|mΩ

a/− Montrer que le problème (1) admet une solution unique σ.

b/− Montrer que σ est caractérisé par :

(2) σ ∈ K ∩ S

Cours d’analyse numérique SMA S6 Pr. BENNOUNA Jaouad

3 TABLE DES MATIÈRES

ou encore 
 σ∈K
(3)
 ∀v ∈ K (σ, v) = 0
0

5-a/− Montrer que la restriction ρ⊥ de ρ à S est un isomorphisme de S sue RN .

5-b/−En déduire que dimS = N .
On suppose que n = 1 et que Ω =]a, b[
6) Montrer que l’ensemble S est égal à l’ensemble des fonctions splines d’ordre m cette à dire
que s vérifié :

(i) s ∈ P2m−1 (]xi , xi+1 [)

(ii) s ∈ Pm−1 ([a, x1 [) et s ∈ Pm−1 (]xN , b])
(iii) s ∈ C 2m−2 ([a, b])
On admet que s ∈ S ⇒ s ∈ C 2m−2 ([a, b]). et on note que x0 = a et xN +1 = b.
7) On suppose que m = 2.
Préciser la nature de l’espace S dans les 4 situations suivantes :

a/ b/ X = {a, x2 , ..., xN }
b/ c/ X = {x1 , ..., xN , b}
c/ d/ X = {a, x1 , ..., xN , b}

Exercice 4. Soit X ∈ {a, x1 ..., xN , b} ⊂ Ω = [a, b]

On introduit ρ : H 2 (Ω) → RN +4 définie par

ρv = (v 0 (a), v(a), v(x1 ), ......, v(xN ), v(b), v 0 (b))

et β = (β00 , β0 , β1 , ......, βN , βN +1 βN
0
+1 )

On introduit comme précédemment les ensembles K, K0 , S et on montre que le problème


 T rouver σ∈K
(1)
 ∀f ∈ K |σ|2,Ω ≤ |f |2,Ω .

admet une solution unique. Et que


 S est l0 espace des f onctions de classe C 2 sur [a, b] tel que
(2)
 ∀i = 0, ..., N S|]xi ,xi+1 [ ∈ P3 (]xi , xi+1 [).

On définit maintenant sur S la base des B-splines cubiques.

Cours d’analyse numérique SMA S6 Pr. BENNOUNA Jaouad

4 TABLE DES MATIÈRES

On suppose les points d’interpolation régulièrement répartie et on introduit des nœuds supplé-
b−a
mentaires. On pose h = et ∀ − 3 ≤ i ≤ N + 4 xi = a + ih
N +1
Pour −1 ≤ i ≤ N + 2 m considère Bi : r → R défini par :



 (x − xi−2 )3 x ∈ [xi−2 , xi−1 ]


 h3 + 3h2 (x − xi−1 ) + 3h(x − xi−1 )2 − 3(x − xi−1 )3 x ∈ [xi−1, xi ]



1 
Bi (x) = 3 h3 + 3h2 (xi+1 − x) + 3h(xi+1 − x)2 − 3(xi+1 − x)3 x ∈ [xi , xi+1 ]
h 
(xi+1 − x)3 x ∈ [xi+1 , xi+2 ]







 0 ailleurs

1) Tracer la B-splines cubique βi .

2) Vérifié que Bi ∈ C 2 (R) et calculer Bi (xj ), B 0 (xj ) pour j = −3, ..., N + 4.
3) Notant encore Bi la restriction de Bi à [a, b], −1 ≤ i ≤ N +2, on introduit l’espace vectoriel
B engendre par {Bi −1 ≤ i ≤ N +2}. Montrer que B = S et que le système {Bi , 1 ≤ i ≤ N +2}
est une base de S.

Cours d’analyse numérique SMA S6 Pr. BENNOUNA Jaouad

5 TABLE DES MATIÈRES

Corrigé de Série N˚3

Exercice 1 :
1) Inf |f − v| dans R+ minoré par 0 donc Inf |f − v| existe.
v∈K v∈K
Soit α = Inf |f − v| ⇒ ∃xn ⊂ K tel que |f − xn | → α.
v∈K
1 1
En effet ∀ ∃xn ∈ K tel que α ≤ |f − xn | ≤ α + (propriété de la borne Inf).
n n
Montrons que (xn )n>0 est une suite de Couchy :
|xn − xm |2 = |(f − xm ) − (f − xn )|2 = 2|f − xn |2 + 2|f − xm |2 − |2f − xm − xn |2
due à l’identité du parallélogramme
 a + b 2  a − b 2 1
   
2 2


 2  +  2  = 2 |a| + |b|
   

xn + xm
xn ∈ K, xm ∈ K ⇒ ∈ K car K est convexe. On a donc
2
xn + xm
|f − | ≥ α ⇒ |2f − xn − xm | ≥ 2α
2

donc |xn − xm |2 ≤ 2|f − xn |2 + 2|f − xm |2 − 4α2 −→ 4α2 − 4α2 = 0

n,m→∞
⇒ lim |xn − xm | = 0, donc (xn ) suite de Cauchy.
n,m→∞
donc xn → u donc u ∈ K car K fermé.
|f − xn | → |f − u| = α car lim |f − xn | = |lim(f − xn )| = |f − u| par continuité de |.|.
Unicité : Supposons du’il existe u0 ∈ K tel que |f − u0 | = α,

u + u0
u, u0 ∈ K ⇒ |f − | ≥ α|2f − u − u0 | ≥ 2α
2

|u − u0 | ≤ 2|f − u|2 + 2|f − u0 |2 − 4α2 = 4α2 − 4α2 = 0

D’où u = u0 .
2/- Montrons l’équivalence entre (1) et (2) :
Soit u ∈ K vérifiant (1) et soit w ∈ K on a

v = (1 + t)u + tw pour t ∈]0, 1].

Donc |f − v| ≤ |f − (1 − t)u − tw| = |(f − u) − t(w − u)|,

donc |f − u|2 ≤ |f − v|2 ≤ |f − u|2 − 2t(f − u, w − u) + t2 |w − u|2 .
C’est-à-dire 2(f − u, w − u) ≤ t|w − u|2
quand t → 0 on obtient
(f − u, w − u) ≤ 0, ∀w ∈ K.

Cours d’analyse numérique SMA S6 Pr. BENNOUNA Jaouad

6 TABLE DES MATIÈRES

Inversement : Soit u vérifient (2). Alors on a :

|u − f |2 − |v − f |2 = 2(f − u, v − u) − |u − v|2 ≤ 0, ∀ ∈ K.

D’où |u − f |2 ≤ |v − f |2 ∀ ∈ K ⇒ (1).
3/- Posons u1 = pK f1 et u2 = pK f2 il vient

(f1 − u1 , v − u1 ) ≤ 0 ∀v ∈ K,
(f2 − u2 , v − u2 ) ≤ 0 ∀v ∈ K.

Posons v = u2 dans la 1ère inégalité et v = u1 dans la 2ème on obtient après addition :

(f1 − f2 + u2 − u1 , u2 − u1 ) ≤ 0 ⇒ |u1 − u2 |2 ≤ (f1 − f2 , u1 − u2 )

Par suite |u1 − u2 | ≤ |f1 − f2 |.

4/- D’après (2) on a
(f − u, v − u) ≤ 0, ∀v ∈ M.

Soit v = w + u avec w ∈ M . on aura

(f − u, w) ≤ 0, ∀w ∈ M.

Soit v = −w + u avec w ∈ M . on aura

(f − u, −w) ≤ 0, ∀w ∈ M ⇒ (f − u, w) ≥ 0 ∀w ∈ M.

 (f − u, w) ≤ 0
⇒ ⇒ (f − u, w) = 0 ∀w ∈ M.
 (f − u, w) ≥ 0

Inversement si 
 u ∈ M tel que
 (f − u, w) = 0 ∀v ∈ M.

On a (f − u, v − u) = 0 ∀v ∈ M d’où (3).
Exercice 2
La fonction splines cherchée est un polynôme de d˚3 sur chaque intervalle [xi−1 , xi ]
i = 2, ..., N on pose
σi00 (x) = ax + b avec σi00 (xi ) = mi .
xi − x x − xi−1
D’où σi00 (x) = mi−1 + mi d’où a, b avec hi = xi − xi−1 , mi = σ 00 (xi ) à
hi hi
déterminer .
En intégrant 2 fois on obtient
(xi − x)3 (x − xi−1 )3
σi (x)mi−1 + mi + Ci (xi−1 − x) + Di (xi − x).
6hi 6hi

Cours d’analyse numérique SMA S6 Pr. BENNOUNA Jaouad

7 TABLE DES MATIÈRES

Les conditions σi (xi ) = βi et σi (xxi−1 ) = βi−1 donne

h2i h2
mi−1 + Di hi = βi−1 et mi i − Ci hi = βi
6 6
(xi − x)3 (x − xi−1 )3 xi−1 − x h2 xi − x h2
d’où σi (x) = mi−1 + mi + (mi i βi ) + (βi−1 − mi−1 i )
6hi 6hi hi 6 hi 6
σ ∈ C 2 [a, b] ⇒ σi0 (xi ) = σi+1
0
(xi ) i = 2, N − 1
0 (xi − x)2 (x − xi−1 )2 1 h2i 1 h2i
σ (x) = −mi−1 + mi − (mi − βi ) − (βi−1 − mi−1 )
2hi 2hi hi 6 hi 6
0 hi 1 hi
σi (xi ) = mi + (βi − βi−1 ) + (mi−1 − mi )
2 hi 6

0 hi+1 1 hi+1
σi+1 (xi ) = −mi + (βi+1 − βi ) + (mi − mi+1 )
2 hi+1 6
d’où en regroupant on obtient :
hi hi + hi+1 hi+1 βi+1 − βi βi − βi−1
mi−1 + mi + mi+1 = − pour i = 2, N − 1
6 3 6 hi+1 hi
m1 = σ 00 (x1 ) = σ 00 (a) = 0 et mN = σ 00 (xN ) = σ 00 (b) = 0.

Les mi i = 2, ..., N − 1 sont les solution du système linéaire Am = B,

βi+1 − βi βi − βi−1
avec m = (m2 , ..., mN −1 )t et B = (b2 , ..., bN −1 )t où bi = − .
hi+1 hi
et A ∈ Mn−2,n−2 R

 h +h h 
2 3 3
0 0 ··· 0 0
 3 6 
 
 
 
 h3 h3 + h4 h4 
 0 ··· 0 0 

 6 3 6 

 
 
 h4 .. ..


0 . . ... :

6
 
 
A=
 

 
 .. .. hN −2 
 : : . . 0 

 6 

 
 
 ... hN −2 hN −2 + hN −1 hN −1 
 : : 
6 3 6
 
 
 
 
hN −1 hN −1 + hN
 
0 0 ... ... 0
6 3

La matrice A est tridiagonale symétrique à diagonale strictement dominante cette à dire

i) les termes aii sont positifs.
n−1
hi + hi+1 X hi + hi+1
(ii) |aii | = > |aij | =
3 j=2
6
j6=i

Cours d’analyse numérique SMA S6 Pr. BENNOUNA Jaouad

8 TABLE DES MATIÈRES

Elle est donc inversible d’où la détermination des mi par suite de σ.

Remarque : les conditions x1 = a et xN = b signifie que

σ ∈ P1 [a, x2 [, σ ∈ P1 ]xN −1 , b] ⇒ σ (j) (a) = σ (j) (b) = 0, j = 2, ..., 2 × 2 − 1

Exercice 3
1) |v(xi )| ≤ ||v||C 0 [a,b] = sup |v(x)| ≤ C||v||m,Ω
x∈[a,b]
−1
2) • K est fermé car ρ {β} = K (ρ continue ⇔ l’image réciproque de tout fermé est fermé).
• K convexe car αK + µk ⊂ K avec α + µ = 1. c’est-à-dire ∀u, v ∈ K αu + µv ∈ K ? avec
α + µ = 1. En effet ρ(αu + µv) = (α + µ)β = β.
• K non vide car :

∀β ∈ R ∃p ∈ Pm−1 tel que p(xi ) = βi 1≤i≤N

donc ρ(p) = β avec p ∈ H m (Ω)

3) K0 fermé. En effet, K0 est un s.e.v c’est-à-dire αu + µv ∈ k0 , ∀u, v ∈ K0 et α, µ ∈ R ?

ρ(αu + µv) = αρ(u) + µρ(v) = 0

4)-a/ Puisque σ ∈ K ρ(σ) = ρ(v) ∀v ∈ K

σ vérifie (1) ⇔ (σ ∈ K et ∀v ∈ K, |||σ||| ≤ |||v|||)

⇔ |||0 − σ||| = min|||0 − v|||.
v∈K

D’après (Exercice 1 série 3) σ est la projection de l’origine sur K, donc unique car K convexe
fermé. ⇒ existence et unicité de σ car σ = proj (0) (pour la norme|||.|||)
K
b/ σ ∈ K par hypothèse σ est caractérise par :

 (((0 − σ, v − σ))) ≤ 0
(d’après Exercice 1)
 ∀v ∈ K

Posons v = w + σ w ∈ K0 (−σ, w + σ − σ)m,Ω ≤ 0 =⇒ (σ, w)m,Ω ≥ 0

Posons v = −w + σ w ∈ K0 (−σ, −w + σ − σ)m,Ω ≤ 0 =⇒ (σ, w)m,Ω ≥ 0
⇒ (σ, v − σ)m,Ω = 0 et comme v − σ ∈ K0 et v ∈ K. Alors (σ, v) = 0 ∀v ∈ K0 .
5/a/ ρ⊥ est un isomorphisme (Rappel : X et Y sont isomorphisme s’il existe une bijective f
de X sur Y tel que f et f −1 sont continue. Si X et Y sont de Banach f est linéaire continue
surjective ⇒ f −1 continue)
Soit ∀v ∈ S tel que ρ(v) = 0 ⇒ v ∈ K0 ∩ S or ∀w ∈ K0 (v, w)m,Ω = 0
⇒ (v, v)m = 0 or |||v||| =< pv >2 +|u|2m,Ω ⇒ |||v||| = 0 ⇒ v = 0 car |||.||| est une norme.

Cours d’analyse numérique SMA S6 Pr. BENNOUNA Jaouad

9 TABLE DES MATIÈRES

d’après 4/ ∀β ∈ RN ∃!σ ∈ S ρ(σ) = β

∀β ∈ RN ∃σ ∈ S ρ⊥ (β) = σ donc p⊥ surjective.
β = ρ⊥ (v1 ) = ρ⊥ (v2 ) ⇒ ∃!σ tel que σ = v1 = v2 ⇒ ρ⊥ injective. D’où ρ est bijective.
D’après 1) ρ⊥ linéaire continue ⇒ ρ−1
⊥ continue ⇒ ρ⊥ isomorphisme.
5/-b) dimS = N .
6) ⇒) supposons que∀s ∈ S, ∀v ∈ K0 (s, v) = 0 et montrons (i) , (ii) et (iii). En effet,
 ϕ(x) si x ∈]xi , xi+1 [ avec ϕ ∈ C ∞ (]xi , xi+1 [) 0≤i≤N
soit s ∈ S fixé et v =
 0 ailleurs
⇒ v ∈ K0
Alors Z xi+1
0 = (s, v)m,Ω = s(m) (t)ϕ(m) (t)dt
 xi xi+1 Z x(i+1)
(m−1) (m)
= ϕ s − ϕ(m−1) sm+1
Z xi xi
m−1 (1) (2m−1)
= (−1) ϕ s
 Z 
m−1
 (2m−1) xi+1 2m
= (−1) ϕs xi
ϕs , ∀ϕ.

⇒ s(2m) = 0
⇒ s ∈ P2m−1 (]xi , xi+1 [), d’où (i).
(ii)- Soit s ∈ S fixé et considérons

 ϕ x ∈ [a, x1 [ ϕ ∈ C ∞ ([a, x1 [)
v=
 0 ailleurs

⇒ Zv ∈ K0
x1
(m−1)
0= s(m) ϕ(m) dt = −s(m) (a)ϕ(a) + ... + (−1)m s(2m−1) (a)ϕ(a)
a
car d’après (i), s(2m) = 0 sur chaque ]xi , xi+1 [.
∀j = 0, ..., m − 1, On choisit ϕ ∈ C ∞ [a, x1 [ tel que :

(ϕ(j) (a) 6= 0) et (∀k = 0, ..., m − 1/k 6= j ϕ(k) (a) = 0)

On en déduit que s(m) (a) = ...... = s(2m−1) (a) = 0 ⇒ s ∈ Pm−1 [a, x1 [.

Le raisonnement est le même sur ]xN , b] d’où (ii).
(iii) On admet que s ∈ C 2m−2 (Ω) = C 2m−2 ([a, b])
⇐) Montrons que si s est une fonction vérifiant (i) , (ii) et (iii) alors s ∈ S

d’après (iii) s ∈ C 2m−2 [a, b] ⇒ s ∈ C m−1 [a, b]


⇒ s ∈ H m (Ω)
et d’après (i)

Cours d’analyse numérique SMA S6 Pr. BENNOUNA Jaouad

10 TABLE DES MATIÈRES

il reste à montrer que ∀v ∈ K0 (s, v) = 0. En effet,

d’après la proposition ??
Z b N
X
(2m−1) −
(s, v)m,Ω = (m) (m)
s (t)v (t)dt = (−1)m
(s(2m+1) (x+
i )−s (xi )v(xi )
a i=1

avec v(xi ) = 0 i = 1, ..., N . On déduit que (s, v) = 0.

7/ Soit s ∈ S on montrer comme précédemment que dans tous les cas :

s ∈ P3 (]xi , xi+1 [) ∀i = 0, ..., N

Posons Ω∗ = Ω̄ \ X = [a, b] \ X, il est clair que C ∞ (Ω∗ ) ⊂ K0 .

Notant encore : ∀i = 0, ..., N si = s|]xi ,xi+1 [ , on a ∀v ∈ C ∞ (Ω∗ )
N Z
X xi+1
0 = (s, v) = s00i v 00 dt.
i=0 xi

En intégrant par partie on obtient

N
X
0= {s00i (xi+1 )v 0 (xi+1 ) − s00i (xi )v 0 (xi ) − s000 000
i (xi+1 )v(xi+1 ) + si (xi )v(xi )},
i=0

d’où
0 = s000 (a)v(a) − s00 (a)v 0 (a) − s000 (b)v(b) + s00 (b)v 0 (b),

a/ On sait que l’on peut choisir v ∈ C ∞ (Ω∗ ) tel que l’un seulement des 4 nombres v(a), v 0 (a), v(b), v 0 (b)
est non nul.
⇒ s00 (a) = s000 (a) = s00 (b) = s000 (b) = 0

donc s ∈ P1 ([a, x1 [) s ∈ P1 (]xN , b])

b/ On a une condition d’interpolation en a : s(a) = β0
et on obtient s00 (a) = s00 (b) = s000 (b) = 0
donc s ∈ V = {v ∈ P3 ([a, x1 [)/v 00 (b) = 0} s ∈ P1 (]xN ; b])
c/ On a une condition d’interpolation en b : s(b) = βN +1
et on obtient s00 (a) = s000 (a) = s00 (b) = 0
donc s ∈ P1 ([a, x1 [) et s ∈ W = {v ∈ P3 ]xN , b]/v 00 (b) = 0}.
d/ On a des conditions d’interpolation en a et b tel que :

s(a) = β0 et s(b) = βN +1

et on obtient s00 (a) = s00 (b) = 0 , donc s ∈ V, s ∈ W .

Exercice 4 :
1/ Graphe de Bi -spline cubique.

Cours d’analyse numérique SMA S6 Pr. BENNOUNA Jaouad

11 TABLE DES MATIÈRES

2/ Bi ∈ C 2 (R) 


 4 j=i

Bi0 (xj ) = 1 j =i±1



 0 sinon

3
j =i−1


h



Bi0 (xj ) = 0 j=i
3




 − j =i+1
h
3/-a/ Tout élément de B apparient à C 2 ([a, b]) est un polynôme de degré 3 sur chaque intervalle
partielle [xi , xx+1 ], 0 ≤ i ≤ N , donc d’après (2) B⊂S
b/ Montrons que :
N
X +2
(?) ∀s ∈ S ∃σi ∈ R −1≤i≤N +2: s= αi Bi
i=1

On sait que la restriction ρ⊥ de p à S est un isomorphisme de S sur RN +4 . Donc (?) revient à

montrer
∀β ∈ RN +4 ∃s ∈ S tel que ρ(s) = β
N +2
X 
Cette à dire ∀β ∈ RN +4 ∃αi ∈ R −1≤i≤N +2 ρ αi Bi = β
i=1
Notant α = (α−1 , α0 , α1 , ..., αN +2 ) et utilisons la question 1/ on obtient l’équation matricielle

Aα = β.

Avec A est définie par

Cours d’analyse numérique SMA S6 Pr. BENNOUNA Jaouad

12 TABLE DES MATIÈRES

 
h h
− 0 0 ··· 0 0 

 3 3 

 1 4 1 0 ··· 0 0  
 
 0 1 4 1 ... : : 
A=
 
 : ... ... ... 
 : 0 0  

 : ... 
 : 1 4 1 0  
 h h 
0 0 ... 0 − 0
3 3
On vérifie que detA 6= 0 Pour cela dans detA , on ajoute la 1 colonne à la 3ème colonne et la
ère

(N + 4)ème à la (N + 2)ème en développent detA par rapport à la 1ère ligne puis par rapport à la
9
dernier ligne , on voit que detA = − 2 ∆ , où ∆ est la déterminant d’une matrice à diagonale
h
strictement dominante d’une régulière.
donc A est régulière et S ⊂ B
n
X
|Aii | ≥ |Aij |
j=1,j6=i 
 β 0 = α−1 B−1
0
(a) + α0 B00 (a) + α1 B10 (a)
 0


β0 = α−1 B−1 (a) + α0 B0 (a) + α1 B1 (a)



 β = α B (x ) + α B (x ) + α B (x )
1 0 0 1 1 1 1 2 2 1

3 3

 β00 = − α−1 + 0 + α1



 h h


 0

 β = α −1 + 4α 0 + α 1

⇔ β1 = α0 + 4α1 + α2






 :


 :

Cours d’analyse numérique SMA S6 Pr. BENNOUNA Jaouad