USTHB - Faculté de Mathématiques - 2 Lic TEL - Sections A & B - MNA - 2016/2017

Epreuve …nale

Exercice 1 : On considère le processus itératif

8
< x0 2 [0; 1]
(P) 1
: xn+1 = ' (xn ) = ' (x) =
2 (x2 x + 1)

1. Montrer que le processus (P) converge vers l’unique point …xe de la fonction ' dans l’intervalle [0; 1] :
1
2. Pour x0 = ; calculer x2 et déduire un encadrement de .
2
Exercice 2 : On considère la table de valeurs
1 2
xi 0 3 3
1
9 9
f (xi ) 1 7 7
1

associée à la fonction f dé…nie dans [0; 1] par

1
f (x) =
x2 x+1
On donne
X
3
14 X 3 4
3 X3
98
x2i = ; xi = ; x4i =
i=0
9 i=0 3 i=0
81
1. Orthogonaliser la base f1; x; x2 g de P2 relativement au produit scalaire dé…ni par

X
3
hU; V i = U (xi ) V (xi ) U; V étant deux fonctions.
i=0
_
2. En utilisant cette base orthogonale, trouver P la meilleure approximation au sens des moindres carrés
dans P2 associée à la table précédente.
3-a) Calculer l’erreur quadratique
X3 _ 2
f (xi ) P (xi )
i=0
_
3-b) Que peut-on en déduire ? Donner alors l’expression de l’erreur e(x) = f (x) P(x) pour x 2 [0; 1] en
fonction de la dérivée quatrième de la fonction f:
Exercice 3 : Soit f est une fonction continue dans [0; 1] : On considère la formule de quadrature donnée
par Z 1
f (x) dx ' A0 f ( ) + A1 f (1 )
0
1
où A0 ; A1 sont des réels et est un paramètre réel, 0 < < :
2
1. Trouver les valeurs de A0 et A1 pour que la formule soit exacte dans P1 :
2. Trouver dans ce cas, les valeurs du paramètre pour que la formule soit exacte dans P2 ?
1
3. Quelles formules de quadrature reconnait-on pour = 0 et = ?
3
Barème : 5.5 points + 10 points + 4.5 points
 USTHB - Faculté de Mathématiques - 2 Lic TEL - Sections A & B - MNA - 2016/2017
Correction épreuve …nale

Exercice 1 : 5.5 points

1. On doit véri…er les hypothèses du théorème du point …xe. On a
1 2x 1
'0 (x) = : 2 (0.5)
2 (x x + 1)2
2
00 1 2(x x + 1)2 (2x 1)(2)(2x 1)(x2 x + 1)
' (x) = :
2 (x2 x + 1)4
3x(x 1)
= 2 <0 si x 2 [0; 1] (0.5)+(0.5)
(x x + 1)3

On a donc le tableau de variations généralisé à '0

1
x 0 2
1
'00
1 1 (0.5)
'0 2
& 0 & 2
1 2 1
' 2
% 3
& 2

On obtient
1
K = max j'0 (x)j = <1 (0.5)
x2[0;1] 2
donc ' est contractante dans [0; 1] et

1 2
'([0; 1]) = min '(x); max '(x) = ; [0; 1] (0.5)
x2[0;1] x2[0;1] 2 3

donc [0; 1] est stable par ': Les deux hypothèses du théorème du point …xe sont donc véri…ées et on conclut,
d’une part qu’il existe un unique point …xe de ' dans [0; 1] et d’autre part que le processus itératif
1 1
xn+1 = '(xn ) = : 2
2 xn xn + 1

converge vers pour toute donnée initiale x0 2 [0; 1] : (0.5)

1 2
2. Pour x0 = , on a déja calculé x1 = : De plus
2 3
1 1 1 9 9
x2 = '(x1 ) = : : = = 0; 64286 (0.5)
2 4=9 2=3 + 1 2 7 14

1
Estimation de l’erreur à la deuxième itération K= :
2
1
K 2 9 2 1
jx2 j jx2 x1 j = 1 = ' 0:023810 = (0.5)+(0.5)
1 K 1 2
14 3 42

1 2 1 2 1 2
On a '([0; 1]) = ; donc 2 ; et '0 est négative dans ; donc les termes de la suite (xn )n 1
2 3 2 3 2 3
s’alternent de part et d’autre de ; on a donc d’une part x2 x1 et d’autre part x2 x2 + ; on
obtient l’encadrement (0.5)

0:642860 = x2 x2 + = 0:642860 + 0:023810 = 0:66667

On peut donc écrire que 0:64 0:67 et ' 0:655.

 Exercice 2 : 10 points
1. On note f 0 ; 1 ; 2 g = 0 (x) = 1; 1 (x) = x; 2 (x) = x2 la base canonique de P2 . Utilisons le
procédé d’orthogonalisation de Gramm-Schmidt pour obtenir une base orthogonale (0.5). On a
i = 0 : Initialisation
j=3 j=3
2
X X
2
0 (x) = 0 (x) = 1 et k 0k = 0 (xj ) = 1 = 4 (0.5)
j=0 j=0

i = 1 : Calcul de 1

h 1; 0i
1 (x) = 1 (x) 2 : 0 (x)
k 0k
j=3 j=3
X X
h 1; 0i = 1 (xj ): 0 (xj ) = xj = 2
j=0 j=0
2 1
1 (x) =x :1 = x (0.5)
4 2
2
Calcul de k 1k on a
j=3 j=3 j=3
2
X X 1 2 X 2 1 5
2
k 1k = 1 (xj ) = (xj ) = xj xj + = (0.5)
j=0 j=0
2 j=0
4 9

i = 2 : Calcul de 2

h 2; 0i h 2; 1i
2 (x) = 2 (x) 2 : 0 (x) 2 : 1 (x)
k 0k k 1k
j=3 j=3
X X 14
h 2; 0i = 2 (xj ): 0 (xj ) = x2j =
j=0 j=0
9
j=3 j=3
X X 1 4 1 14 5
h 2; 1i = 2 (xj ): 1 (xj ) = x2j (xj )= : = (0.5)
j=0 j=0
2 3 2 9 9
14=9 5=9 1 1
2 (x) = x2 :1 :(x ) = x2 x+ (0.5)
4 5=9 2 9
2
Calcul de k 2k : On a
j=3 j=3
2
X X 1
2
1
2
1 1 1
2
4
2
k 2k = 2 (xj ) = x2j xj + = + 2+ 2+ = (0.5)+(0.5)
j=0 j=0
9 9 9 9 9 81

2. Le meilleur polynôme au sens des moindres carrés est donné par

hf; i i
P 2 (x) = 0 0 (x) + 1 1 (x) + 2 2 (x) où i = (0.5)
k i k2
On a
X
i=3 X
i=3
9 32
hf; 0i = f (xi ) 0 (xi ) = f (xi ) = 2(1 + ) = (0.5)
i=0 i=0
7 7
X
i=3 X
i=3
1 X i=3
1X
i=3
hf; 1i = f (xi ) 1 (xi ) = f (xi )(xi )= f (xi )xi f (xi )
i=0 i=0
2 i=0
2 i=0
3 6 1 32
= ( + + 1) = 0 (0.5)+(0.5)
7 7 2 7
X
i=3 X
i=3
1 X
i=3 X
i=3
1X
i=3
2
hf; 2 i = f (xi ) 2 (xi ) = f (xi )(xi xi + ) = f (xi )x2i f (xi )xi + f (xi )
i=0 i=0
9 i=0 i=0
9 i=0
1 4 16 1 32 4
= ( + + 1) + : = (0.5)+(0.5)
7 7 7 9 7 63
On a donc
32 1 8 4 81 9 81 9 9 2
0 =
: = ; 0 = 0; 0 = : = : et P 2 (x) = 1 + x x (0.5)
7 4 7 63 4 7 4 7 7
3-a) Calcul de l’erreur quadratique : on a la table des écarts
1 2
xi 0 3 3
1
9 9
f (xi ) 1 7 7
1
9 (0.5)
P 2 (xi ) 1 7
97 1
i 0 0 0 0
on a donc l’erreur quadratique au carré
X
3 _ 2
f (xi ) P (xi ) = 0!! (0.5)
i=0

3-b) On déduit que f (xi ) = P 2 (xi ) pour i = 0; 1; 2; 3 et comme on a degré de P 2 inferieur ou égal
à 3, alors P 2 n’est autre que le polynôme d’interpolation qui est de degré 2 dans ce cas là (0.5). D’où
l’expression de l’erreur d’aproximation
1 2 f (4) ( x )
"(x) = x(x )(x )(x 1) x 2 [0; 1] (0.5)
3 3 4!
Exercice 3 : 4.5 points
1. Calcul des coe¢ cients A0 ; A1 : La formule est exacte dans P1 si et seulement si elle est exacte pour
f (x) = 1 et f (x) = x (0.5), A0 ; A1 doivent donc véri…er le système :
8 1
> Z
>
>
>
> 1dx = A0 + A1 ( (
>
< A0 + A 1 = 1 A0 + A1 = 1 (1)
0
Z1 (0.5) , 1 , 1 (0.5)
>
> A0 + (1 )A1 = A0 + (1 )A1 = (2)
>
> 2 2
>
> xdx = A0 + (1 )A1
:
0

1 1=2) 1 (
(2) (1) , = (2 1)A1 , A1 == :(0.5)
2 1=2) 2 2(
1
En injectant la valeur de A1 dans (1) ; on obtient A0 = : (0.5)
2
2. La formule étant exacte dans P1 , elle sera exacte dans P2 si elle l’est pour f (x) = x2 : On a alors
Z1
1 1 1 1
x2 dx = ( 2
+ (1 )2 ) , = 2
+ , 2
+ = 0 (0.5)
2 3 2 6
0

1 1 1
On a donc deux racines : (1 p ) (0.5);on retient celle qui est dans l’intervalle 0; i:e:
2 3 2
1 1
= 1 p (0.5). D’où la formule de quadrature, exacte dans P2 :
2 3
Z1
1 1 1 1 1 1 1
f (x)dx ' f :p +f + :p
2 2 2 3 2 2 3
0

Z1
1
3. Pour = 0, on obtient la formule fermée des trapèzes (h = 1) : f (x)dx ' (f (0) + f (1)) (0.5)
2
0
Z1
1 1 1 1 2
Pour = , on obtient la formule ouverte des trapèzes (h = ) : f (x)dx ' f +f (0.5)
3 3 2 3 3
0
——————————————————————————————–