[PDF] Analyse Numérique 0 0 Les uii étant non nuls

View - Download Analyse Numérique 0 0

Previous PDF

Next PDF

Universit´e de Nice Sophia-AntipolisLicence L3 Math´ematiques Ann´ee 2008/2009

Analyse Num´eriqueCorrig´e du TD 6

EXERCICE 1

Matrices diagonales, triangulaires

1.1 Matrices diagonales

SoitD= (dii)i=1,...,nune matrice diagonale d"ordren >0. Donner une condition n´ecessaire et suffisante pour queDsoit inversible. On peut repr´esenterDsous forme du tableau suivant : (d 11 ...0 d ii0 d nn))))))))

Commedet D=n?

i=1d ii, on a

Dinversible??det D?= 0??dii?= 0,?i= 1,...,n.

1.2 Matrices triangulaires inf´erieures

SoitL= (lij)i,j=1,...,nune matrice triangulaire inf´erieure d"ordren >0. a. Sous quelle condition n´ecessaire et suffisanteLest-elle inversible?

La matriceLpeut se mettre sous la forme suivant :

(l 11 l

21l22......

l

1ilijlii0

.ln-1n-1 l n1lnj···lnn-1lnn)))))))))))))) d"o`u la matrice triangulaire inf´erieureLpeut ˆetre caract´eris´ee par : 1 Universit´e de Nice Sophia-AntipolisLicence L3 Math´ematiques Ann´ee 2008/2009 lij= 0 sii < j,?i,j= 1,...,n.

Puisquedet L=n?

i=1l ii, on a

Linversible??det L?= 0??lii?= 0,?i= 1,...,n.

b. On suppose que la matrice triangulaire inf´erieureLest inversible. Soitb un vecteur colonne ayantncomposantes. Donner un algorithme qui permet de r´esoudre l"´equation d"inconnuey:

Ly=b.(1.1)

Commelii?= 0,?i= 1,...,n, la r´esolution du syst`eme (1.1) s"´ecrit y1=b1l11, y i=1 lii? b i-i-1? j=1l ijyj? ,?i= 2,...,n. (1.2) Quel est le coˆut de cet algorithme en termes d"op´erations ´el´ementaires (addi- tions, multiplications, divisions) ? Le calcul dey1demande 1 division (div) dans l"algorithme (1.2). Pourifix´e dans{1,...,n}, le calcul deyipar l"algorithme (1.2) requiert 1 division (div), i-1 additions (add) eti-1 multiplications (mult).

Au total le coˆutCLde l"algorithme (1.2) est

C

L= 1 div +n?

i=2? (i-1) add + (i-1) mult + 1 div? = 1 div + n? i=21 div +n-1? k=1kadd +n-1? k=1kmult (n-1)n

2add +(n-1)n2mult +ndiv.

Le nombre d"op´erations ´el´ementairesCLest de l"ordre den2,i.e.CL=O(n2). 2 Universit´e de Nice Sophia-AntipolisLicence L3 Math´ematiques Ann´ee 2008/2009

1.3 Matrices triangulaires sup´erieures

On consid`ere une matrice triangulaire sup´erieureUd"ordren >0 . a. Donner une condition n´ecessaire et suffisante pour queUsoit inversible.

La matriceUpeut se mettre sous la forme suivant :

(u

11u12···u1ju1n-1u1n

u

22u2n-1u2n......

u iiuijuin 0 u n-1n-1un-1n u nn)))))))))))))) d"o`u la matrice triangulaire inf´erieureUpeut ˆetre caract´eris´ee par : u ij= 0 sii > j,?i,j= 1,...,n.

Commedet U=n?

i=1u ii, on a

Uinversible??det U?= 0??uii?= 0,?i= 1,...,n.

b. On suppose que la matrice triangulaire sup´erieureUest inversible. Soity un vecteur colonne donn´e ayantncomposantes. Ecrire un algorithme qui permet de r´esoudre l"´equation d"inconnuex:

U x=y .(1.3)

Lesuii´etant non nuls, l"inconnuexsolution du syst`eme lin´eaire (1.3) est donn´ee par xn=ynunn, x i=1 uii? y i-n? j=i+1u ijyj? ,?i= 1,...,n-1. (1.4)

Donner la complexit´e de cet algorithme.

Le calcul dexnrequiert 1 multiplication (mult) dans l"algorithme (1.4). Pourifix´e dans{1,...,n}, le calcul dexipar l"algorithme (1.4) demande 1 division (div), n-iadditions (add) etn-imultiplications. 3 Universit´e de Nice Sophia-AntipolisLicence L3 Math´ematiques Ann´ee 2008/2009

Par suite le coˆutCUde l"algorithme (1.4) est

C

U= 1 div +n-1?

i=1? (n-i) add + (n-i) mult + 1 div? n-1? k=1kadd +n-1? k=1kmult + 1 div +n-1? i=11 div (n-1)n

2add +(n-1)n2mult +ndiv.

Le nombre d"op´erations ´el´ementairesCUest de l"ordre den2,i.e.CU=O(n2)..

Vocabulaire

L"algorithme (1.2) pour inverser les syst`emes triangulaires inf´erieurs est ditdescente ousubstitution directe. L"algorithme (1.4) pour r´esoudre les syst`emes triangulaires sup´erieurs est ditremont´eeousubstitution r´etrograde.

EXERCICE 2

M´ethode d"´elimination de Gauss

2.1 Des exemples

Effectuer une ´elimination de Gauss sur les syst`eme lin´eaires suivants (2 4 41 3 11 5 6)) (x 1 x 2 x 3)) =((21 -6)) (1 0 6 28 0-2-2

2 9 1 3

2 1-3 10))))

(x 1 x 2 x 3 x 4)))) =((((6 -2 -8 -4))))

Premier exemple

Nous ´ecrivons le premier syst`eme sous la forme du tableau (2 4 4 21 3 1 11 5 6-6)) L 1 L 2 L 3 •On effectue l"´elimination de Gauss. On a successivement (2 4 4 20 1-1 0

0 3 4-7))

L2←L2-0.5L1

L

3←L3-0.5L1(2.1)

4 Universit´e de Nice Sophia-AntipolisLicence L3 Math´ematiques Ann´ee 2008/2009 (2 4 4 20 1-1 0

0 0 7-7))

L

3←L3-3L2

On obtient alors le syst`eme triangulaire suivant

?2x1+ 4x2+ 4x3= 2 x

2-x3= 0

7x3=-7

En utilisant l"algorithme de remont´ee (1.2) on a successivement x

3=-1,x2=-1,x1= 5.

L"´elimination de Gauss ci-dessus est ditesans permutation. •On peut par exemple `a l"´etape (2.1) ci-dessus, remplacer le pivot 1 par le coefficient 3 dex2de la derni`ere ligne, parce que 3>1 donne plus de stabilit´e num´erique. Dans ce cas on dit que l"on fait une ´elimination de Gauss avecpivot partiel. Dans ce contexte on obtient (2 4 4 20 3 4-7

0 1-1 0))

L3←→L2

(2 4 4 20 3 4-7 0 0-7

373)))

L

3←L3-13×L2

D"o`u on obtient le syst`eme triangulaire sup´erieur suivant ?2x1+ 4x2+ 4x3= 2

3x2+ 4x3=-7

7

3x3=73

En appliquant l"algorithme de remont´ee `a ce syst`eme on obtient x

3=-1,x2=-1,x1= 5.

•On peut enfin par exemple `a l"´etape (2.1) ci-dessus, remplacer le pivot 1 par le coefficient

le plus grand en module dans la sous-matrice 1-1 3 4 Ceci rend la m´ethode plus stable num´eriquement. Ici on trouve 4 comme nouveau pivot. Dans ce cas on dit que l"on fait une ´elimination de Gauss avecpivot partiel. Dans ce 5 Universit´e de Nice Sophia-AntipolisLicence L3 Math´ematiques Ann´ee 2008/2009 contexte on obtient((2 4 4 20 3 4-7

0 1-1 0))

L2←→L3

(2 4 4 20 4 3-7quotesdbs_dbs48.pdfusesText_48

[PDF] analyse numérique exercices corrigés méthode de newton pdf

[PDF] analyse numérique exercices et problèmes corrigés

[PDF] analyse numérique interpolation polynomiale exercices corrigés

[PDF] analyse numérique matricielle exercices corrigés pdf

[PDF] analyse numérique matricielle pdf

[PDF] analyse numérique pour ingénieurs

[PDF] analyse physico chimique du lait cru pdf

[PDF] analyse physico chimique du lait de vache pdf

[PDF] analyse physico chimique du lait en poudre

[PDF] analyse physico chimique du lait ppt

[PDF] analyse physico chimique du miel

[PDF] analyse physico-chimique des dattes

[PDF] analyse psychologique gratuite

[PDF] analyse séquentielle des politiques publiques

[PDF] analyse sociologique film ressources humaines