Travaux dirigés d’informatique 4

Serie N ◦ 2

Exercice 1 : (Somme des inverses des carrés d’entiers > 0)

N
X 1
∀n ∈ N∗ , on pose : SN = .
n2
n=1

1) Concevoir et écrire une fonction algorithmique non récursive qui calcule le terme de rang N
arbitraire de la suite SN . N.B. Essayer d’envisager plusieurs approches.

2) On sait que (SN ) est la suite des sommes partielles d’une série numérique convergente. Conce-
voir et écrire alors une fonction algorithmique qui somme les termes de cette série jusqu’à un rang
N tel que la somme SN est une approximation de la somme totale de cette série avec une incertitude
absolue < ε, avec ε réel > 0 donné, et qui renvoie cette approximation comme résultat.

Exercice 2 : (Somme des inverses des factorielles - Calcul du nombre e)

N
X 1
Comme dans l’Exercice 1, mais avec :∀ n ∈ N, SN = .
n!
n=0
N.B. Bien faire attention ici à la problématique posée par la présence de la fonction factorielle.

Exercice 3 : (Calcul approchée d’une intégrale par développement en série)

Z 1 x +∞
e −1 X 1
On pose : I = dx, A = .
0 x n · (n !)
n=1

1) Montrer que I ∈ R et A ∈ R.

2) Démontrer que I = A.

3) Utiliser ce dernier résultat pour construire une méthode numérique permettant de calculer une ap-
proximation de I avec une incertitude absolue < ε, avec ε réel > 0 donné, ainsi qu’une approximation
de l’erreur absolue associée à cette approximation de I.

4) Concevoir et écrire alors une fonction algorithmique mettant cette méthode numérique en œuvre.

Exercice 4 : (Calcul approchée d’une intégrale par développement en série)

Z 1 +∞
1 − cos x X (−1)n−1
On pose : I = dx, S = .
0 x2 (2n) ! · (2n − 1)
n=1

TD Info4 page 1 ENSP

 1) Montrer que I ∈ R et S ∈ R.

2) Démontrer que I = S.

3) Utiliser ce dernier résultat pour construire une méthode numérique permettant de calculer une ap-
proximation de I avec une incertitude absolue < ε, avec ε réel > 0 donné, ainsi qu’une approximation
de l’erreur absolue associée à cette approximation de I.

4) Concevoir et écrire alors une fonction algorithmique mettant cette méthode numérique en œuvre.

Exercice 5 :
π +∞ +∞ Z π
(−1)n (−1)n
Z
2 X X 2
On pose : I = cos(sin x)dx, A = n 2
, B= · Jn , où Jn = (sin x)2n dx.
0 4 (n !) (2n) ! 0
n=0 n=0
L’objectif ici est de calculer, avec une précision fixée d’avance, la valeur de l’intégrale I.

1) Démontrer que : a) A ∈ R; b) I = B; c) ∀n ∈ N∗ , (2n) · Jn = (2n − 1) · Jn−1 .

π
2) Déduire l’égalité : I = 2 · A.

3) Utiliser ce dernier résultat pour construire une méthode numérique permettant de calculer une ap-
proximation de I avec une incertitude absolue < ε, avec ε réel > 0 donné, ainsi qu’une approximation
de l’erreur absolue associée à cette approximation de I.

4) Concevoir et écrire alors une fonction algorithmique mettant cette méthode numérique en œuvre.

Exercice 6 :
Z π +∞  2 Z π
2 X n! 2
On pose : I = sin(cos x)dx, A = (−4)n , Jn = (cos x)2n+1 dx,
0 (2n + 1) ! 0
n=0
+∞
X (−1)n · Jn
B= .
(2n + 1) !
n=0
L’objectif ici est de calculer, avec une précision fixée d’avance, la valeur de l’intégrale I.

1) Démontrer que : a) A ∈ R; b) I = B; c) ∀n ∈ N∗ , (2n + 1) · Jn = (2n) · Jn−1 .

2) Déduire l’égalité : I = A.

3) Utiliser ce dernier résultat pour construire une méthode numérique permettant de calculer une ap-
proximation de I avec une incertitude absolue < ε, avec ε réel > 0 donné, ainsi qu’une approximation
de l’erreur absolue associée à cette approximation de I.

4) Concevoir et écrire alors une fonction algorithmique mettant cette méthode numérique en œuvre.

Exercice 7 : (Calcul du nombre π - Procédé d’accélération d’Aitken)

+∞
π X (−1)n
On sait que : = .
4 (2n + 1)
n=0
1) Utiliser ce résultat et le procédé d’accélération d’Aitken(Cf. Chapitre sur la Sommation numérique
des séries) pour construire une méthode numérique permettant de calculer une approximation de

TD Info4 page 2 ENSP

 π avec une incertitude absolue < ε, avec ε réel > 0 donné, ainsi qu’une approximation absolue
associée à cette approximation de π.

2) Concevoir et écrire alors une fonction algorithmique mettant cette méthode numérique en œuvre.

3) Manipulation des grands nombres dans un programme informatique.

Dans certaines situations, on a besoin de manipuler des nombres de taille beaucoup plus grande
que ceux connus en standard par un langage de programmation. Pour y arriver, la solution consiste
alors à simuler un système de représentation des nombres du type et de la taille appropriés
dans un programme écrit dans ce langage.

Exercice 8 : (Grands entiers dans un programme informatique)

Ici, nous nous intéressons au cas des entiers. Écrire un programme en L.E.A qui permet de manipuler
les entiers pouvant atteindre jusqu’à 50 chiffres en base 10. Plus précisément :

1) Construire un type auxiliaire approprié pour la situation, et qu’on appellera big− entier.

2) Écrire des fonctions algorithmiques booléennes pour traduire respectivement, dans ce cadre, les com-
paraisons habituelles entre entiers, i.e. =, 6=, <, >, 6, > .

3) Écrire des fonctions algorithmiques pour effectuer respectivement, dans ce cadre, les opérations
arithmétiques usuelles sur les entiers, i.e. opposé, addition, soustraction, multiplication, reste et
quotient de la division euclidienne.

4) Écrire une procédure algorithmique de lecture d’un tel entier.

5) Écrire une procédure algorithmique d’écriture d’un tel entier.

N.B. Cependant, on évitera de bloquer le programme sur la base 10 et la taille 50. Écrire le
programme de telle sorte que ces 2 valeurs puissent être facilement changées.

Exercice 9 : (Grands nombres réels dans un programme informatique)

Ici, nous nous intéressons au cas des nombres réels. Écrire un programme en L.E.A qui simule la
manipulation d’un système de représentation de réels-machine en virgule flottante normalisée avec une
mantisse de 50 chiffres en base 10, et un exposant à 15 chiffres. Plus précisément :

1) Construire un type auxiliaire approprié pour la situation, et qu’on appellera big− réel.

2) Écrire des fonctions algorithmiques booléennes pour traduire respectivement, dans ce cadre, les com-
paraisons habituelles entre réels, i.e. =, 6=, <, >, 6, > .

3) Écrire des fonctions algorithmiques pour effectuer respectivement, dans ce cadre, les opérations
arithmétiques usuelles sur les réels, i.e. opposé, addition, soustraction, multiplication, division.

4) Écrire une procédure algorithmique de lecture d’un tel réel.

5) Écrire une procédure algorithmique d’écriture d’un tel réel.

TD Info4 page 3 ENSP

 6) (Partie subsidiaire) Utiliser ce programme pour reprendre l’Exercice 7 et calculer π avec la
précision maximale qu’il permet. A-t-on obtenu ainsi les 50 premières décimales de ce nombre ? Si
non, où est le problème et comment le contourner ?

N.B. Cependant, on évitera de bloquer le programme sur la base 10 et la taille de mantisse

50, et la taille d’exposant 15. Écrire le programme de telle sorte que ces 3 valeurs puissent
être facilement changées.

Exercice 10 : (Opérations de base sur les polynômes)

Écrire un programme en L.E.A. qui permet de manipuler les polynômes de divers degrés pouvant
atteindre 100. Plus précisément :

1) Construire un type auxiliaire approprié pour la situation, et qu’on appellera polynome.

2) Écrire une fonction algorithmique booléenne pour évaluer l’égalité entre 2 polynômes.

3) Écrire des fonctions algorithmiques pour effectuer respectivement, les opérations arithmétiques usuelles
sur les polynômes, i.e. opposé, dérivation, primitive, intégration, addition, soustraction, multiplica-
tion, reste et quotient de la division euclidienne.

4) Évaluation d’un polynôme en un réel donné.

5) Écrire une procédure algorithmique d’écriture d’un tel polynôme.

N.B. Cependant, on évitera de bloquer le programme sur le degrés maximal 100. Écrire le
programme de telle sorte que cette valeur puisse être facilement changée.

Exercice 11 : (Évaluation d’un polynôme pair)

Dans la suite de l’Exercice précédent, concevoir et écrire une fonction algorithmique efficace en un
réel donné d’un polynôme dont on sait d’avance qu’il est pair.

Exercice 12 : (Évaluation d’un polynôme impair)

Dans la suite de l’Exercice précédent, concevoir et écrire une fonction algorithmique efficace en un
réel donné d’un polynôme dont on sait d’avance qu’il est impair.

Exercice 13 : (Évaluation d’un polynôme de Taylor)

L’objectif ici est l’évaluation algorithmique efficace, en un réel x, du polynôme de Taylor d’ordre n
d’une fonction f , centré en un point x0 :
n
X f (k) (x0 )
y = T (x) = (x − x0 )k ,
k!
k=0

Où n, x et x0 sont donnés, ainsi que les coefficients : b0 = f (x0 ), b1 = f 0 (x0 ), · · · , bn = f (n) (x0 ).

TD Info4 page 4 ENSP

 (x − x0 )k
1) On pose uk = . Trouver une relation de récurrence entre uk et uk−1 .
k!
2) Écrire une fonction algorithmique d’évaluation numérique de y qui s’inspire de celui de l’évaluation
d’un polynôme P (x), écrit sous forme canonique, par calcul des puissances successives de x.

3) (a) Trouver l’analogue efficace du schéma de Hörner pour un polynôme de la forme de T (x).

(b) En déduire une fonction algorithmique correspondante pour l’évaluation d’un tel polynôme.

4) Comparer l’efficacité numérique relative des 2 algorithmes obtenus pour évaluer T (x).

Exercice 14 : (Polynôme d’interpolation de Lagrange : Forme de Newton)

1) Écrire une fonction algorithmique en L.E.A. qui prend en entrée :

1. une suite finie de nombre réels distincts : x0 , · · · , xn ;

2. la suite finie des images correspondantes pour une certaine fonction f (inconnue) : y0 , · · · , yn ,
et renvoie comme résultat la suite des différences divisées appropriées pour évaluer le polynôme
d’interpolation de Lagrange de f relativement aux points d’interpolation x0 , · · · , xn .

2) Écrire une fonction algorithmique en L.E.A. qui prend en entrée :

1. une suite finie de nombre réels distincts : x0 , · · · , xn ;

2. la suite finie des images correspondantes pour une certaine fonction f (inconnue) : y0 , · · · , yn ;

3. la suite des différences divisées issu de la fonction algorithmique précédente ;

4. un réel x arbitraire,

et renvoie comme résultat la valeur en x du polynôme d’interpolation de Lagrange de f relativement

aux points d’interpolation x0 , · · · , xn en s’appuyant sur sa forme de Newton.

N.B. Proposer 2 versions de cette dernière fonction algorithmique.

Exercice 15 : (Polynôme d’interpolation de Lagrange : Algorithme d’Aitken)

Écrire une fonction algorithmique en L.E.A. qui prend en entrée :

1) une suite finie de nombre réels distincts : x0 , · · · , xn ;

2) la suite finie des images correspondantes pour une certaine fonction f (inconnue) : y0 , · · · , yn ;

3) un réel x arbitraire,
et renvoie comme résultat la valeur en x du polynôme d’interpolation de Lagrange de f relativement
aux points d’interpolation x0 , · · · , xn , valeur calculée par l’algorithme d’Aitken.

Exercice 16 : (Interpolation de Lagrange adaptative par la forme de Newton)

Écrire une fonction algorithmique en L.E.A. qui prend en entrée :

1) une suite finie de nombre réels distincts : a0 , · · · , aN ;

TD Info4 page 5 ENSP

 2) la suite finie des images correspondantes pour une certaine fonction f (inconnue) : y0 , · · · , yN ;

3) un réel x arbitraire,
et qui calcule une valeur approchée de y = f (x) :

1. en calculant successivement les valeurs respectives en x des polynômes d’interpolation de La-

grange de f relativement à x0 , (x0 , x1 ), (x0 , x1 , x2 ), etc,

2. où les points x0 , x1 , · · · doivent être pris, dans cet ordre, parmi a0 , · · · , aN ; par ordre de
proximité décroissante avec x

3. et les valeurs des polynômes d’interpolation doivent être calculées en utilisant leur forme de
Newton.
Ce calcul itératif doit être répété jusqu’à convergence des valeurs calculées avec une incertitude
relative < 10−3 si c’est possible, et renvoi d’un message d’erreur sinon.

Exercice 17 : (Interpolation de Lagrange par l’algorithme d’Aitken)

Comme dans l’Exercice 16, mais en utilisant plutôt l’algorithme d’Aitken pour l’évaluation en x des
polynômes d’interpolation de Lagrange successifs.

Exercice 18 : (Polynôme d’interpolation d’Hermite : Forme de Newton)

1) Écrire une fonction algorithme en L.E.A qui prend en entrée :

1. Une suite finie de nombres réels distincts : z0 , . . . , zn ;

2. La suite finie des images correspondantes pour une certaine fonction f (inconnue) : y0 , . . . , ym ;

3. La suite finie des images correspondantes pour la fonction dérivée f 0 : d0 , . . . , dm ;

et renvoie comme résultat la suite des différences divisées approchées pour évaluer le polynôme
d’interpolation d’Hermite de f relativement aux couples (z0 ; 1), . . . , (zm ; 1).

2) Écrire une fonction algorithme en L.E.A. qui prend en entrée :

1. Une suite finie de nombres réels distincts : z0 , . . . , zn ;

2. La suite finie des images correspondantes pour une certaine fonction f (inconnue) : y0 , . . . , ym ;

3. La suite finie des images correspondantes pour la fonction dérivée f 0 : d0 , . . . , dm ;

4. la suite des différences divisées issu de la fonction algorithme précédente ;

TD Info4 page 6 ENSP

 5. un réel x arbitraire,

et renvoie comme résultat la valeur en x du polynôme d’interpolation d’Hermite de f relativement

aux points d’interpolation (z0 ; 1), . . . , (zm ; 1) en s’appuyant sur sa forme de Newton.
N.B. Proposer 2 versions de cette dernière fonction algorithme.

Exercice 19 : (Polynôme d’interpolation d’Hermite : Algorithme d’Aitken)

Écrire une fonction algorithmique en L.E.A. qui prend en entrée :

1. Une suite finie de nombres réels distincts : z0 , . . . , zn ;

2. La suite finie des images correspondantes pour une certaine fonction f (inconnue) : y0 , . . . , ym ;

3. La suite finie des images correspondantes pour la fonction dérivée f 0 : d0 , . . . , dm ;

4. un réel x arbitraire,

et renvoie comme résultat la valeur en x du polynôme d’interpolation de Lagrange de f relativement aux

points d’interpolation x0 , . . . , xn , valeurs calculée par l’algorithme d’Aitken.
N.B. Évidemment, il faudra d’abord expliquer l’adaptation de l’algorithme d’Aitken dans ce contexte.

Exercice 20 : (Interpolation d’Hermite adaptative par la forme de Newton)

Écrire une fonction algorithmique en L.E.A. qui prend en entrée

1. Une suite finie de nombres réels distincts : a0 , . . . , aN ;

2. La suite finie des images correspondantes pour une certaine fonction f (inconnue) : y0 , . . . , yN ;

3. La suite finie des images correspondantes pour la fonction dérivée f 0 : d0 , . . . , dN ;

4. un réel x arbitraire,

et qui calcule une valeur approchée de y = f (x) :

1. en calculant successivement les valeurs respectives en x des polynômes d’interpolation d’Hermite de

f relativement à x0 , (x0 , x1 ), (x0 , x1 , x2 ), etc,

2. Où les points x0 , x1 , . . . doivent être pris, dans cet ordre, parmi a0 , . . . , aN (après duplication préalable)
par ordre de proximité décroissante avec x,

TD Info4 page 7 ENSP

 3. et les valeurs dezs polynômes d’interpolation doivent être calculées en utilisant leur forme de Newton.

Ce calcul itératif doit être répété jusqu’a convergence des valeurs calculées avec une incertitude relative
< 10−3 si c’est possible, et renvoi d’un message d’erreur sinon.

Exercice 21 : (Interpolation d’Hermite adaptative par l’algorithme d’Aitken)

Comme dans l’Exercice 20, mais en utilisant plutôt l’algorithme d’Aitken pour l’évaluation en x des
polynômes d’interpolation d’Hermite successifs.

Exercice 22 : (Fonction d’approximation affine par morceaux)

Écrire une fonction algorithmique en L.E.A. qui prend en entrée :

1. une suite finie de nombre réels représentant une subdivision d’un intervalle [a, b] : a0 , . . . , aN ;

2. la suite finie des images correspondantes pour une certaine fonction f (inconnue) : y0 , . . . , yN ;

3. un réel x arbitraire,

et renvoie comme résultat la valeur en x de la fonction d’approximation affine par morceaux de f sur
[a, b] s’appuyant sur la subdivision (a0 , . . . , aN ), si x ∈ [a, b], et renvoie un message d’erreur sinon.

Exercice 23 : (Fonction d’approximation parabolique par morceaux)

Écrire une fonction algorithmique en L.E.A. qui prend en entrée :

1. une suite finie de nombre réels représentant une subdivision d’un intervalle [a, b] : a0 , . . . , aN ;

2. la suite finie des images correspondantes pour une certaine fonction f (inconnue) : y0 , . . . , yN ;

3. un réel x arbitraire,

et renvoie comme résultat la valeur en x de la fonction d’approximation parabolique par morceaux de f

sur [a, b] s’appuyant sur la subdivision (a0 , . . . , aN ), si x ∈ [a, b], et renvoie un message d’erreur sinon.

Exercice 24 : (Fonction d’approximation d’Hermite par morceaux)

Écrire une fonction algorithmique en L.E.A. qui prend en entrée :

1. une suite finie de nombre réels représentant une subdivision d’un intervalle [a, b] : a0 , . . . , aN ;

2. la suite finie des images correspondantes pour une certaine fonction f (inconnue) : y0 , . . . , yN ;

TD Info4 page 8 ENSP

 3. la suite finie des images correspondantes pour la fonction dérivée f 0 (inconnue) : d0 , . . . , dN ;

4. un réel x arbitraire,

et renvoie comme résultat la valeur en x de la fonction d’approximation d’Hermite par morceaux de f sur
[a, b] s’appuyant sur ces données, si x ∈ [a, b], et renvoie un message d’erreur sinon.

Exercice 25 : (Fonction d’approximation spline cubique)

Écrire une fonction algorithmique en L.E.A. qui prend en entrée :

1. une suite finie de nombre réels représentant une subdivision d’un intervalle [a, b] : a0 , . . . , aN ;

2. la suite finie des images correspondantes pour une certaine fonction f (inconnue) : y0 , . . . , yN ;

et construit la fonction d’approximation spline sur [a, b] s’appuyant sur la subdivision (a0 , . . . , aN ).

Exercice 26 : (Polynômes matrices)

N.B.Ceci est une partie d’un problème d’examen, d’où une certaine redondance avec
d’autres Exercices.
Dans ce problème :

1. n et p sont des constantes entières ¿0 fixées ;

2. a0 , . . . , an sont des nombres réels ;

3. A0 , . . . , An sont les matrices carrées réelles d’ordre p.

L’objectif ici est alors d’examiner le problème de l’évaluation algorithmique efficace d’un
polynôme de degré n d’une variable réelle ou matricielle.
I- Calcul de P (x) = an xn + an−1 xn−1 + . . . + a1 x + a0 en un réel x donné.

1) a) Écrire l’algorithme d’évaluation numérique de P (x) basé sur le calcul de puissance de x.

b) Coût de cet algorithme ?

2) a) Écrire un algorithme de calcul de P (x) basé plutôt sur le fait qu’on ait aussi (schéma de
HÖrner) :
P (x) = (. . . ((an x + an−1 )x + an−2 )x + . . . + a1 )x + a0 .

b) Coût numérique de cet algorithme ?

3) Comparer l’efficacité numérique relative de ces 2 algorithmes d’évaluation de P (x).

TD Info4 page 9 ENSP

II- Quelques fonctions algorithmiques matricielles de base.

1) Écrire les fonctions algorithmiques respectives :

a) P ROD− M AT faisant le produit de 2 matrices carrées réelles d’ordre p ;

b) ADD− M AT faisant la somme de 2 matrices carrées réelles d’ordre p ;

c) SCAL− M AT faisant la multiplication d’une matrice carrée réelle d’ordre p par un réel donné ;

d) M AT− V EC faisant le produit d’une matrice carrée réelle d’ordre p par un vecteur, dans cet
ordre par un vecteur, dans cet ordre ;

e) V EC− M AT faisant le produit d’un vecteur par une matrice carrée réelle d’ordre p, dans cet
ordre ;

2) Evaluer le coût numérique de chacune de ces fonctions algorithmiques.

*** NOTA : Dans toute la suite de ce problème, les fonctions algorithmiques de II seront
supposées écrites, et pourront être appelées partout où le besoin s’en ferait sentir.
III-Calcul de F (M ) = an M n + an−1 M n−1 + . . . + a1 M + a0 Ip pour M ∈ Mp (R).
Adapter I dans ce cas (sachant que Ip est la matrice d’ordre P ).
IV-Calcul de G(M ) = An M n + An−1 M n−1 + . . . + A1 M + A0 pour M ∈ Mp (R).
Adapter I dans ce cas.

8) Tracé de courbes dans un plan

9) Tris de tableaux.

10) Résolution numérique des systèmes d’équations linéaires.

11) Intégration numérique.

12) Résolution numérique d’une équation à une inconnue.

Exercice 27 : (Tri d’une suite finie par sélection-échange)

Trier une suite réelle finie, c’est ré-arranger ses éléments dans l’ordre pré-déterminé. Écrire une
fonction algorithme en L.E.A. qui fait ce travail, dans l’ordre croissant, de la manière suivante :

TD Info4 page 10 ENSP

 1) recherche du plus grand élément dans la suite ;

2) positionnement de cet élément là où il faut ;

3) répéter la manœuvre de manière appropriée jusqu’à ce que le résultat soit OK.

N.B. Proposer une version non récursive et une version récursive de cette fonction algorithmique.

TD Info4 page 11 ENSP