Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Examen SR 2021

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 7

**NB : Vous êtes appelés à utiliser :

1. une flèche −→ pour indiquer une indentation,


2. l’importation des modules “numpy“,“matplotlib.pyplot“ et “sympy“ en utilisant les alias
présentés ci-dessous :

[ ]: import numpy as np
import matplotlib.pyplot as plt
import sympy as sp

Exercice 1 : (10 points)


Z 2
Dans cet exercice, on s’intéresse à approcher la valeur de l’intégrale I = f ( x ) dx en utilisant
1
2x
deux formules de quadrature différentes, avec f ( x ) = .
1 + x2

1. Donner la valeur exacte de I en utilisant une fonction pré-définie dans la bibliothèque


sympy.(1 point)

[ ]: #Code:

La valeur exacte de I vaut.............................................................................................................

1
2. On veut approcher la valeur de I moyennant la formule de quadrature suivante:

3
Q( f ) = C0 f (1) + C1 f ( ) + C2 f (2),
2

avec C0 , C1 et C2 des nombres réels.


Sachant que C0 , C1 et C2 vérifient le système d’équations linéaires (S) suivant:

C0 + C1 + C2
 =1
(S) C0 + 32 C1 + 2C2 = 32

C0 + 49 C1 + 4C2 = 73

(a) Vérifier que le système (S) admet dans R3 une unique solution en utilisant une fonction
pré-définie.(1 point)

[ ]: #Code:

Le système (S) admet dans R3 une unique solution car........................................................


(b) Déterminer les valeurs de C0 , C1 et C2 en utilisant la fonction solve du sous-module
linalg de numpy.(1 point)

2
[ ]: #Code

• La valeur de C0 (4 chiffres après la virgule) est ...............................................................


• La valeur de C1 (4 chiffres après la virgule) est ...............................................................
• La valeur de C2 (4 chiffres après la virgule) est ...............................................................
(c) Pour ces valeurs de C0 , C1 et C2 , calculer l’erreur d’intégration de I par la formule de
quadrature Q.(1,5 point)

[ ]: #Code:

L’erreur d’intégration (4 chiffres après la virgule) par la formule de quadrature Q vaut


........................................................................................................................................................

3. On s’intéresse maintenant à approcher la valeur de I par la méthode composite des rectan-


gles à gauche.
On vous rappelle que la méthode composite des rectangles à gauche est utilisée pour ap-
Z b
procher I ( f ) = f ( x ) dx, où f est une fonction continue sur [ a, b], et dont la formule est
a
donnée par:

n −1
c
IRg (f) = h ∑ f ( x k ),
k =0

avec xk = a + kh (0 ≤ k ≤ n) les points d’intégration issus de la subdivision de [ a, b] en n


b−a
sous intervalles avec un pas fixe h = .
n

(a) Ecrire une fonction gauchecomposite(f,a,b,n) prenant en entrée, une fonction f con-
tinue sur [ a, b], les bornes de l’intervalle a et b et le nombre de sous intervalles n et
retourne la valeur approchée de I par IRg c ( f ).(1,5 points)

3
[ ]: def gauchecomposite(f,a,b,n):

(b) En considérant 5 sous-intervalles, calculer l’erreur d’intégration de I par la méthode


composite des rectangles à gauche.(1 point)

[ ]: #Code:

L’erreur d’intégration (4 chiffres après la virgule) par la méthode des rectangles à


gauche pour n = 5 est ...............................................................................................................
(c) Comparer la méthode des rectangles à gauche pour n = 5 et la formule de quadrature
Q pour l’approximation de I en termes de précision.(1 point)

........................................................................................................................................................

........................................................................................................................................................
c par la formule
4. On rappelle que, d’après l’estimation de l’erreur d’intégration de I, notée ERg
composite des rectangles à gauche,
c
ERg ≤ M,
2−2x2
avec M = 1
2n max | f ′ ( x )| et f ′ ( x ) = (1+ x 2 )2
.
x ∈[1,2]

(a) En considérant la fonction g définie sur [1, 2] par g( x ) = | f ′ ( x )|, écrire les instructions
permettant de tracer la courbe de g sur [1, 2].(1 point)

[ ]: x= # discrétiser [1,2] en 100 nœuds.

plt.grid(True)
plt.xlabel('axe des abscisses ')
plt.ylabel('axe des ordonnées')
plt.title('Courbe de g')

4
(b) Déduire graphiquement une valeur approchée de max g( x ), puis donner la valeur de
x ∈[1,2]
M pour n = 5.(1 point)

max g( x ) = ............................. et M = ................................


x ∈[1,2]

Exercice 2 : (10 points)


 
On considère une fonction f continue sur a, b telle que f ( a) ̸= f (b).
On souhaite approcher numériquement la racine d’une équation ( E) : f ( x ) = 0 en utilisant le
schéma itératif suivant:

b−a

n +1 = x n − f ( x n ),
x
( MI )n : f ( b ) − f ( a)
x0 ∈ [ a, b],

avec ( xn ) est une suite convergente vers la solution exacte x ∗ de ( E) ( lim xn = x ∗ ).


n→+∞

1. Pour atteindre cet objectif, compléter la fonction methodeiterative(a,b,f,x0,epsilon)


qui prend les paramètres suivants: a et b les deux extrémités de [ a, b], la fonction f,
la valeur initiale x0 et la tolérance epsilon. Cette fonction doit retourner deux listes
dont la première contient tous les itérés x = [ x0 , x1 , x2 , · · · ] et la deuxième contient
y = [ f ( x0 ), f ( x1 ), f ( x2 ), · · · ] .(3 points)
Dans cette fonction, on utilise le test d’arrêt f( xn ) < epsilon, ∀ n ≥ 0.

[ ]:
def methodeiterative(a,b,x0,f,epsilon):
x=[x0]
y=[f(x0)]

while (abs(y[-1])>epsilon) :

return x,y

5
On considère la fonction f continue sur [1, 2], par

f ( x ) = ln(1 + x2 ) − sin( x ).

NB: La fonction logarithme s’écrit sous python comme suit : np.log() ou sp.log().
2. (a) Ecrire les instructions qui justifient que ( E) admet au moins une solution dans
l’intervalle [1,2].(1,5 point)

[ ]:

Justification: .................................................................................................................................
.........................................................................................................................................................
.........................................................................................................................................................
(b) Donner l’instruction nécessaire pour déterminer l’expression de la dérivée de f , notée
d f .(1 point)

[ ]:

L’expression de la dérivée de f est

d f (x) = · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · ·

(c) Représenter graphiquement d f , la fonction dérivée de f sur [1,2].(0,75 point)

[ ]: # Donner les instructions nécessaires

plt.grid(True)
plt.xlabel('axe des abscisses ')
plt.ylabel('axe des ordonnées')
plt.title('Courbe de df')

6
 
(d) Déduire que ( E) admet une unique solution sur 1, 2 .(1 point)

........................................................................................................................................................

........................................................................................................................................................

Prenant dans la suite x0 = 2 et epsilon = 10−2 .


3. (a) Donner l’instruction qui permet d’afficher la solution approchée de x ⋆ à epsilon près
en utilisant la fonction "méthodeiterative".(1 point)

[ ]: # code

x ⋆ =................................................................................................................................................

(b) Remplir le tableau ci-dessous qui contient les trois premiers itérés (deux chiffres
seulement après la virgule) du schéma ( MI )n .(1,75 points)

X x0 = 2 x1 = x2 = x3 =
Y f ( x0 ) = f ( x1 ) = f ( x2 ) = f ( x3 ) =

Vous aimerez peut-être aussi