ECE2-B 2019-2020

TP4 : Écriture de sommes finies en Scilab

Pré-requis : l’objectif des premières séances de TP est de faire le point sur des fonctionnalités impor-
tantes de Scilab qui ont été vues en première année. Je vous invite à consulter les chapitres de cours
correspondants sur ma page : support informatique.
On pourra en particulier se reporter au « CH 7 : Les structures itératives ».
I Dans votre dossier Info_2a, créer le dossier TP_4.

I. Avant propos
On considère dans ce TP les suites (un )n2N et (vn ) suivantes :
8
n2 < v0 = 1
8n 2 N, un = n et 2 vn
3 : 8n 2 N, vn+1 =
evn +e vn

Objectif du TP : il s’agit d’explorer les différentes méthodes permettant le calcul des

P
n P
n
sommes partielles d’ordre n : Sn = uk et Tn = vk .
k=0 k=0

II. Calcul des sommes partielles d’ordre n

II.1. Calcul de Sn
Il s’agit ici d’illustrer le cas où la suite (un ) est donnée sous forme explicite.

II.1.a) Méthode itérative

I Écrire une fonction calculSn qui :
⇥ prend en paramètre un entier n,
⇥ renvoie une variable S,
⇥ à l’aide d’une structure itérative, calcule la valeur de Tn et stocke le résultat dans S.

1 function S = calculSn(n)
2 S = 0
3 for i = 0:n
4 S = S + i ^ 2 / (3 ^ i)
5 end
6 endfunction

I Que vaut S0 ? S5 ? S10 ? S100 ? S1000 ?

On trouve : S0 = 0, S5 ' 1.4115, S10 ' 1.4989, S100 ' 1.5 et S1000 ' 1.5.
(noter que l’affichage du résultat se fait avec un nombre fixé de chiffres après la virgule)

1
ECE2-B 2019-2020

II.1.b) Utilisation des fonctionnalités Scilab

I Écrire une fonction premSuiteU qui :
⇥ prend en paramètre un entier n,
⇥ renvoie une variable U, vecteur contenant initialement n + 1 zéros,
⇥ à l’aide d’une structure itérative, calcule les n + 1 premières valeurs de la suite (un ) et stocke le
résultat dans U.

1 function U = premSuiteU(n)
2 U = zeros(1,n+1)
3 for i = 1:n+1
4 U(i) = i ^ 2 / (3 ^ i)
5 end
6 endfunction

I Que réalise l’appel sum(1:5) ? Détailler le rôle de la fonction sum.

La fonction sum prend en paramètre une matrice et renvoie la somme de tous ses
cœfficients. Ici, on considère le vecteur [1, 2, 3, 4, 5] dont la somme vaut 15.

I En déduire un appel permettant de calculer S10 .

sum(premSuiteU(10))

I Que réalise l’appel 1 ./ (1:5) ? Et (4:8) ./ (1:5) ? Détailler le rôle de l’opérateur ./

L’opérateur ./ est l’opérateur de division terme à terme.

• Dans le premier cas, on obtient une valeur approchée de [ 11 , 12 , 13 , 14 , 15 ].
• Dans le deuxième cas, les deux matrices étant de même taille, on obtient une matrice de
même taille où les cœfficients sont obtenus par division des cœfficients en même position.
Plus précisément, on obtient une valeur approchée de [ 41 , 42 , 43 , 44 , 45 ].

I De même, que réalise l’appel [1,2; 1,1] ^ 2 ? Et l’appel [1,2; 1,1]. ^ 2 ? Et (1:5) ^ 2 ?

• L’opérateur . ^ est l’opérateur de puissance terme à terme. L’opérateur ^

est
l’opérateur de puissance mathématique classique.
• L’appel (1:5) ^ 2 est donc impropre : on ne peut mutiplier la matrice [1, 2, 3, 4, 5] par
elle-même du fait de ses dimensions. Cependant, Scilab récupère cette erreur en élevant
au carré tous les cœfficients de [1, 2, 3, 4, 5] ce qui correspond en fait à l’appel (1:5). ^ 2.

I Comment peut-on, par simple manipulation matricielle, créer le vecteur U des 101 premiers éléments
de la suite (un ) ? On commencera par stocker le vecteur 0:100 dans une variable N.

N = 0:100; U = (N . ^ 2) ./ (3 . ^ N)
(l’appel U = (N ^ 2) ./ (3 ^ N) convient aussi)

2
ECE2-B 2019-2020

I Comment obtient-on alors S100 à l’aide de U ? Et comment récupérer S5 ?

• La valeur de S100 est donnée par l’appel sum(U).

• S5 est la somme des 6 premiers éléments de (un ) (ne pas oublier 0). On récupère le
vecteur contenant ces 6 éléments par l’appel U(1:6) et donc S5 par l’appel sum(U(1:6)).

II.2. Calcul de Tn
Il s’agit ici d’illustrer le cas où la suite (un ) est donnée sous forme récurrente.

II.2.a) Méthode itérative

I Écrire une fonction calculTn qui :
⇥ prend en paramètre un entier n,
⇥ renvoie une variable T,

⇥ à l’aide d’une structure itérative, calcule la valeur de Sn et stocke le résultat dans T.

On pourra utiliser une variable auxiliaire v afin de calculer les différents termes de (vn ).

1 function T = calculTn(n)
2 T = 1
3 v = 1
4 for i = 1:n
5 v = 2 ? v / (exp(v) + exp(-v))
6 S = S + v
7 end
8 endfunction

I Que vaut T0 ? T100 ? T10000 ?

On trouve : T0 = 1, T100 ' 18.18, T10000 ' 197.58.

II.2.b) Utilisation des fonctionnalités Scilab

I Écrire une fonction premSuiteV qui :
⇥ prend en paramètre un entier n,
⇥ renvoie une variable V, vecteur contenant initialement n + 1 zéros,
⇥ à l’aide d’une structure itérative, calcule les n + 1 premières valeurs de la suite (vn ) et stocke le
résultat dans V.

1 function V = premSuiteV(n)
2 V = zeros(1,n+1)
3 V(1) = 1
4 for i = 1:n
5 V(i+1) = 2 ? V(i) / (exp(V(i)) + exp(-V(i)))
6 end
7 endfunction

3
ECE2-B 2019-2020

I En déduire un appel permettant de calculer T10 .

sum(premSuiteV(10))

III. Calcul des n premières sommes partielles

P
III.1. Calcul des n premières sommes partielles de un
I Soit n 2 N. Quel lien y a-t-il entre Sn et Sn+1 ?

(n + 1)2
Sn+1 = Sn +
3n+1

I En tirant profit de l’égalité précédente, écrire une fonction premSn qui :

⇥ prend en paramètre une variable n,
⇥ renvoie une variable tabS, vecteur contenant initialement n + 1 zéros,
⇥ à l’aide d’une structure itérative, stocke la valeur de Si dans la ième case de tabS.
On ne devra pas effectuer d’appel à calculSn mais on pourra s’inspirer de son code.

1 function tabS = premSn(n)

2 tabS = zeros(1, n+1)
3 tabS(1) = 0
4 for i = 1:n
5 tabS(i+1) = tabS(i) + i ^ 2 / 3 ^ i
6 end
7 endfunction

I Comme précédemment, on aurait aussi pu tirer parti des fonctionnalités de Scilab.

Que réalise l’appel cumsum(1:5) ? Détailler le rôle de la fonction cumsum.

• La fonction cumsum (Cumulative Sum) prend en paramètre un vecteur u et renvoie un

vecteur v de même taille dont le ième cœfficient est la somme des i premiers éléments de u.
• Ainsi, cumsum(1:5) renvoie le vecteur [1, 3, 6, 10, 15].

I Quel appel, tirant parti des fonctionnalités Scilab, permet d’obtenir les 101 premiers éléments de
la suite (Sn ) ? On utilisera la fonction cumsum.

• Si on a codé la fonction premSuiteU, on peut réaliser l’appel :

cumsum(premSuiteU(100)).
• Sinon on peut, commme précédemment, créer le vecteur qui contient les 101 premiers
éléments de la suite (un ) : N = 0:100; U = (N . ^ 2) ./ (3 . ^ N); cumsum(U).

4
ECE2-B 2019-2020

P
III.2. Calcul des n premières sommes partielles de vn
I Soit n 2 N. Quel lien y a-t-il entre Tn et Tn+1 ?

Tn+1 = Tn + vn+1

I En tirant profit de l’égalité précédente, écrire une fonction premTn qui :

⇥ prend en paramètre une variable n,
⇥ renvoie une variable tabT, vecteur contenant initialement n + 1 zéros,
⇥ à l’aide d’une structure itérative, stocke la valeur de Ti dans la ième case de tabT.
On ne devra pas effectuer d’appel à calculTn mais on pourra s’inspirer de son code. On pourra
notamment utiliser une variable v calculant les valeurs successives des termes de la suite (vn ).

1 function tabT = premTn(n)

2 tabT = zeros(1, n+1)
3 v = 1
4 tabT(1) = 1
5 for i = 1:n
6 v = 2 ? v / (exp(v) + exp(-v))
7 tabT(i+1) = tabT(i) + v
8 end
9 endfunction

 À retenir : on a étudié deux méthodes en Scilab pour obtenir les éléments de (Sn ).
1) La suite des sommes partielles étant une suite (grande nouvelle), on peut se servir des
procédés vus en TP2.
2) On peut aussi préférer créer le vecteur des premiers éléments de la suite (un ) et utiliser
les instructions sum ou cumsum suivant ce que l’on cherche à obtenir. Encore une fois,
il faut se reporter au TP2.

P
IV. Tracé des sommes partielles de un
I Compléter le programme suivant afin qu’il permette d’effectuer le tracé des 51 premiers éléments
de (Sn ). On exécutera ce programme.

1 N = 0:50
2 U = (N . ^ 2) ./ (3 . ^ N)
3 tabS = cumsum(U)
4 plot(N, tabS, 'rx')
P
I Quelle conjecture peut-on émettre sur la nature de la série un ?
P
D’après la représentation graphique, on peut émettre l’hypothèse que un est une série
convergente, de somme 32 .

5
ECE2-B 2019-2020

V. Suite des sommes partielles aux concours

V.1. ECRICOME 2015
En TP2, on a déjà introduit l’épreuve ECRICOME 2015 qui commençait par l’étude d’une suite récurrente
(un )n2N⇤ , définie par une relation de récurrence de type un+1 = F (un ).
La première question, consistait à compléter le programme permettant le calcul des 100 premiers
éléments de (un ). On rappelle ce programme ci-dessous.

1 U = zeros(1,100)
2 U(1) = 1
3 for n = 1 : 99
4 U(n+1) = 1 - exp(-U(n))
5 end
6 plot(U,"+")

I On modifie le programme précédent en remplaçant la dernière ligne par :

1 X = 1 : 100
2 S = cumsum(U)
3 Y = log(X)
4 plot2d(X,S)
5 plot2d(X,Y)
Le programme ci-dessus permet d’obtenir la représentation graphique suivante :

I Que représente le vecteur-ligne S ?

Quelle conjecture pouvez-vous émettre sur la nature de la série de terme général un ?

• Le vecteur U contient les 100 premiers éléments de la suite (un ). L’opérateur cumsum
permet de calculer la somme P cumulée de ce vecteur. Ainsi, S contient les 100 premières
sommes partielles de la série un .
P
• D’après le tracé obtenu, on peut émettre l’hypothèse que la série un est divergente.
Plus précisément, que Sn ! +1.
n!+1

6
ECE2-B 2019-2020

I À l’aide de la question 2.h) (consistait à démontrer : 8n 2 N⇤ , un > n1 ) établir la nature de la série

de terme général un .

On sait :
⇥ 8n 2 N, 0 6 n1 6 un .
P 1
⇥ La série n est une série de Riemann d’exposant 1 (1 6 >1).
n>1
P 1
Ainsi, la série n est divergente.
n>1
P
Donc, par le théorème de comparaison des séries à termes positifs, la série un diverge.

V.2. EML 2015

On considère l’application f : R ! R, x 7! f (x) = x3 ex
et la suite réelle (un )n2N définie par : u0 = 1 et 8n 2 N, un+1 = f (un ).
P 1
Il était demandé de démontrer que la série converge (de somme S) et que :
n>1 f (n)

P
n 1 1
8n 2 N⇤ , S 6
k=1 f (k) (e 1) en
I En déduire une fonction Scilab qui calcule une valeur approchée de S à 10 4 près. (cf TP3)

1
• Si on sait : 6 10 4, on obtient par transitivité : |S Sn | 6 10 4.
(e 1) en
1
L’idée est donc de trouver le premier entier tel que : 6 10 4.
(e 1) en
On peut démontrer que cet entier vaut : d4 ln(10) ln(e 1)e (= 9) et ainsi, utiliser une
boucle for pour calculer S9 qui fournit l’approximation souhaitée.

1 function S = calcApprochS()
2 n = ceil(4 ? log(10) - log(exp(1)-1))
3 S = 0
4 for k = 1:n
5 S = S + 1 / (k ^ 3 ? exp(k))
6 end
7 endfunction

• Le deuxième choix est de calculer les valeurs successives de (Sn ) tant que
1 1
6 10 4 n’est pas vérifiée i.e. tant que > 10 4 .
(e 1) e n (e 1) en

1 function S = calcApprochS()
2 n = 1
3 S = 1 / exp(1)
4 while 1 / ((exp(1)-1) ? exp(n)) > 10 ^
(-4)
5 n = n + 1
6 S = S + 1 / (n ^ 3 ? exp(n))
7 end
8 endfunction

7
ECE2-B 2019-2020

V.3. ECRICOME 2018

Pn 1
Pour tout entier naturel n non nul, on pose : un = ln(n).
k=1 k

I Écrire une fonction d’en-tête : function y = u(n) qui prend en argument un entier naturel n non
nul et qui renvoie la valeur de un .

1 function y = u(n)
2 S = 0
3 for k = 1:n
4 S = S + 1/k
5 end
6 y = S - log(n)
7 endfunction

Détaillons les différents éléments de ce code :

P
n
⇥ en ligne 2, on crée la variable S dont le but est de contenir, en fin de programme k.
1
k=1
Cette variable S est donc initialisée à 0.
⇥ de la ligne 3 à la ligne 5, on met à jour la variable S à l’aide d’une boucle.
Pour ce faire, on ajoute au k ème tour de boucle la quantité k1 .
P
n
Ainsi, S contient bien k en sortie de boucle.
1
k=1
P
n
⇥ en ligne 6, on affecte à la variable y la valeur un = 1
k ln(n).
k=1

Commentaire
P
n
Pour le calcul de la somme k,
1
on peut aussi tirer profit des fonctionnalités Scilab :
k=1

S = sum(1 ./ 1:n)

Pour bien comprendre cette instruction, rappelons que :

⇥ l’instruction 1:n permet de créer la matrice ligne (1 2 ... n).
⇥ l’opérateur ./ permet d’effectuer la division terme à terme.
Ainsi, l’instruction 1 ./ 1:n permet de créer la matrice ligne ( 11 1
2 ... 1
n ).
⇥ la fonction sum permet de sommer tous les coefficients d’une matrice.
On obtient donc bien la somme à calculer par cette méthode.