Transcript Dvt 256-255-238-406-433 - 432 Approximation de Pi
238 255 256 406
Le nombre π . Algorithmes d'approximation du nombre Exemples de comport t π . Vitesse de convergence, accélération de convergence. asympt. de suites, rapidité de cv. Approximation de Pi par la méthode des polygones réguliers d'Archimède. Rombaldi De Biasi.
433
Approximations du nombre π . Archimède a inventé, vers 250 avant J-C, une méthode originale pour le calcul de la longueur d'un cercle. Il encadre en effet cette valeur par le périmètre d'un polygone régulier inscrit dans ce cercle, et par le périmètre d'un polygone régulier exinscrit. 1°) Considérons
a n
=
n
sin π , pour tout n≥1. Mq la cv vers π est lente. (énoncé modifié – question 1 rajoutée – à retenir) Pourquoi cette suite est-elle inutilisable telle quelle pour calculer une approximation du nombre π ? 2°) On introduit :
x n
=
n
2 sin π 2
n
, pour tout n≥1. Mq cette suite converge vers π géométriquement de raison 1 4 .
I.
3°) On l’accélère en posant
y n
= 4
x n
+ 1 −
x n
3 4°) Qu’obtient-on en itérant ce procédé ?
1 , pour tout n≥1. Mq on a une convergence géométrique de raison 16 .
Outils A. Formule d’Al-Kashi (ou « loi des cosinus », ou « théorème de Pythagore généralisé »). (De Biasi p.254)
Soit un triangle ABC, dans lequel on utilise les notations usuelles exposées sur la figure 1 : d'une part α, β et γ pour les angles et, d'autre part, a , b et c pour les longueurs des côtés respectivement opposés à ces angles. Alors
c
2 =
a
2 +
b
2 − 2
ab
. cos γ Démo :
c
2 =
AB
2 =
CB
−
CA
2 =
CB
2 − 2 +
CA
2 =
CB
2 +
CA
2 − 2 . cos(
ACB
) =
a
2 +
b
2 − 2
ab
. cos γ Il existe de nombreuses autre démonstrations, en particulier purement géométriques, dont l’une due à Al Kashi, une autre à Pythagore, etc…
B. Méthode d’Archimède des polygones réguliers. (pas de source)
On se propose de calculer le périmètre d'un cercle de rayon 1. Pour un polygone régulier à n côtés, l'angle au centre vaut 2 π
n
2 π Par la formule d'Al-Kashi, on a :
c
2 =
n
= 4 sin 2 . Pour le polygone inscrit, on a la figure ci-dessus. π , soit
c
= 2 sin π . Pour le polygone exinscrit, on a la figure : On a donc :
c
= 2 tan π . On en déduit que : 2 sin π ≤ 2 π
n
≤ 2 tan , d’où
n
sin π ≤ π ≤
n
tan
II.
238; 255; 256; 406; 433
Développement
Vitesse de cv de suites / nbre π
1°).
a n
=
n
sin π
, convergence lente. (pas de source, à retenir)
Approximation de π
a n
=
n
sin 0 π
n
− π 3 6 1 .
n
3 +
o
1
n
3 = π − π 3 6 1 .
n
2 +
o
1
n
2 π Rombaldi, Bia. De plus,
e n
=
x n
− π = − π 3 6 .
1
n
2 +
o
1
n
2 . D’où
e n e n
+ 1 ∼
n
2
n
2 + 1 → 1 , donc la convergence est lente. Problème : on utilise π pour calculer les éléments de la suite, mais comme c’est π que l’on souhaite approcher, la définition est inutilisable. En s'inspirant de l'idée développée au paragraphe II.B (256), on va considérer les termes d'indices 2 n .
2°).
x n
=
n
2 sin π 2
n
, convergence géométrique de raison 1/4. (Rombaldi + De Biasi p.141)
Pour n∊ℕ*, considérant le polygone régulier à 2n côtés, on pose : (tend vers π comme sous-suite de la précédente)
x n
=
n
2 sin π 2
n
→ π On a : sin 2 π 2
n
+ 1 = 2 π 2
n
2 2 π 2
n
= 1 − 1 − 2
x n
2
n
2 , Donc
x n
+ 1 = 2
n
+ 1 2 1 − 1 −
x n
2
n
2 =
n
2.2 . 1 − 1 −
x n
2
n
2 , Et cette expression permet un calcul itératif des termes de la suite, sans utiliser le nombre π que l’on cherche à approcher. En utilisant le développement limité : sin
x n
= π − π 3 1 3! 2 2
n
+ π 5 1 5! 2 4
n
+
o
1 2 4
n
= π =
x
−
x
3 3!
+ − π 3!
β 3 2 1 λ 2
x
5 5!
+
n
, il vient le développement asymptotique : + π 5!
γ 5 2 1 µ 4
n
+
o
1 2 4
n
= π +
n
+
n
+
o
Où β, γ sont inconnus, car ils contiennent π que l’on cherche à approximer, et λ = 1 4 et µ 1 = 16 . La suite converge bien vers π, géométriquement de raison λ = 1 . 4
3°).
y n
= 4
x n
+ 1 −
x n
, convergence géométrique de raison 1/16 : méthode d’acc° de Richardson. (Rombaldi)
3 On va introduire un barycentre de x n et x n+1 , pondéré astucieusement pour éliminer β (i.e.une combinaison linéaire), et donc éliminer le facteur λ n . Cela nous amènera à une cv. géométrique de raison µ (plus rapide).
238; 255; 256; 406; 433
x n
= π + β .
n
Vitesse de cv de suites / nbre π + γ .
1 16
n
+
o
1 16
n
x n
+ 1 = π + β 1 . .
4
n
+ γ 1 .
16 .
1 16
n
+
o
1 16
n
+ 1 Approximation de π 4
x n
+ 1 = 4 π + β .
4
n
+ γ 1 . .
4 1 16
n
+
o
1 16
n
+ 1 , pour éliminer dans la soustraction Rombaldi, Bia. 4
x n
+ 1 −
x n
= 3 π γ .
− 3 4 .
1 16
n
+
o
1 16
n
, que l'on divise par 3 pour que ça tende encore vers π
y n
= 4
x n
+ 1 −
x n
3 = .
− 1 .
1 16
n
+
o
1 16
n
γ ' La suite obtenue tend vers π, géométriquement de raison 1/16. On a accéléré la suite (x n ) pour obtenir une convergence géométrique de raison µ 1 = 16 en posant :
y n
= 4
x n
+ 1 −
x n
3
4°). Itération de la méthode d’accélération de Richardson. (Rombaldi)
On a
x n
=
n
2 sin π 2
n
. En effectuant un DL en 0 de
f
= sin ( π
x
) (où
x
1
x
= 2
n
), il vient (pour p∊ℕ*): sin ( π
x
) =
x p
+ ∑ 1
j
= 0
x n
= + 1
p
∑
j
= 0
j
π 2
j
+ 1 ( 2
j
+
x
2
j
+ 2
p
+ 2 ) , d'où le développement asymptotique (en posant
j
( 2 π
j
2
j
+ + 1 ) 1
n
2
j
+
o
1 2
n
2
p
+ 2 = π + 1
p
∑
j
= 0
j
π 2
j
( 2
j
+ ) 1
n
2
j
+
o
1 2
n
2
p
+ 2
x
= 1 2
n
) : donc
x n
= π +
p
+ 1 ∑
j
= 1
j
( 2 π
j
2 +
j
+ 1 ) β
j
1 4
j
λ
j
n
+
o
1 4
p
+ 1
n
. En identifiant avec les notations générales:
x n
= α +
p
+ 1 ∑ β λ
j n j j
= 1 +
o
( λ
p
+ 1
n
) , on a α = π , β
j j
π 2
j
+ 1 ( 2
j
+ et λ =
j
1 . Notons que les coefficients β
j
sont inconnus, puisque l'on cherche à approcher π . D'après les formules de récurrence
x
x x n
,0 =
x n
π =
n
2 sin 2
n
= 1
x n
+ 1,
k
− 1 − 4
k
1 1 − 4
k x
− 1 = 4
k x n
+ 1,
k
− 1 −
x
4
k
− 1
x n
,0 =
x n
=
x n
+ 1,
k
− 1 − λ
k x
1 − λ
k
− 1 de la méthode de Richardson, il vient ici les suites accélératrice: − 1 pour
k
∈ 1;
p
et
n
≥ 1 . D'où:
x n x n
,1
x n
, 2 = = =
n
2 sin 4 16
x n
+ 1 3
x n
+ 1,1 15 − π −
x n
x n
,1
238; 255; 256; 406; 433
Vitesse de cv de suites / nbre π Approximation de π Si l'on pose
f
= sin ( π
x
) , et
x
= , alors du DL
f
= + 1
p
∑
j
= 0
j
π 2
j
+ 1 ( 2
j
+
t
=
x
2 : =
p
+ ∑ 1
j
= 0
j
( 2 π 2
j
+ 1
j
+ ) .
t j
+ = π α + + 1
p
∑
j
= 1
j
( 2 π 2
j
+ 1
j
+ ) .
t
β
j j
+ , on retrouve le "cas général" décrit dans le cours, avec
r
= 1 4 . En effet,
x n
=
f
1
n
.
x
2
j
+ 2
p
+ 2 ) Rombaldi, Bia. , on déduit en posant On peut donc appliquer la formule générale de calcul d'équivalent de l'erreur, et il vient: π −
x
π 2
k
+ 3 ∼ ∞ ( 2
k
+ ( 2 1 )(
k
+ 1 ) , d'où les évaluations d'erreurs: π π π − − −
x n x n
,1
x n
,2 ∼ ∞ ∼ ∞ ∼ ∞ π π π 3 1 6 4
n
5 7 4 5! 16
n
+ 1 1 7! 64
n
+ 1
Attention, le Rombaldi contient quelques erreurs d'indices, il faut vérifier les calculs.
On trouve les valeurs de quelques termes dans le De Biasi:
x
5 = 3,136 548 489 (précision 5.10
-3 )
x
4,2 = 3,141 582 936 (précision 10 -5 )
x
3,3 = 3,141 592 618 (précision 4.10
-8 )
x
2,4 = 3,141 592 653 (précision 10 -9 ) Attention aux notations:
x
n° du terme, n° de la suite .