Suite récurrente linéaire

Restez toujours informé : suivez-nous sur Google (☆)

Introduction

En mathématiques, on appelle suite récurrente linéaire d’ordre p, toute suite à valeurs dans un corps K (généralement $C$ ou $R$ ) définie pour tout $n \geq n_{0}$ par la relation de récurrence suivante :

a₀, a₁, …a_{p − 1} étant p scalaires fixés de K (a₀ non nul), pour tout $n \geq n_{0}$ , on a

u_{n + p} = a_{0} u_{n} + a_{1} u_{n + 1} + \dots + a_{p - 1} u_{n + p - 1}

Une telle suite est entièrement déterminée par la donnée des p premiers termes de la suite et par la relation de récurrence

Les suites récurrentes linéaires d’ordre 1 s’appellent plus simplement des suites géométriques de raison a₀. Les suites récurrentes linéaires d’ordre 2 sont entièrement connues et leur terme général est déterminé en fonction des coefficients a₀ et a₁. Une des suites de ce type est la très célèbre suite de Fibonacci. L’étude des suites récurrentes linéaires d’ordre p fait appel à la notion d’espace vectoriel et au calcul matriciel.

Suite récurrente linéaire d’ordre 1

Voir article détaillé : Suite géométrique

Si la relation de récurrence est u_{n + 1} = q**u_n, le terme général est $u_{n} = u_{n_{0}} q^{n - n_{0}}$

Suite récurrente linéaire d’ordre 2

a et b étant deux scalaires fixés de K avec b non nul, la relation de récurrence est

u_{n + 2} = a**u_{n + 1} + b**u_n (R)

On va prouver que le terme général d'une telle suite est

$λ r_{1}^{n} + μ r_{2}^{n}$ si r₁ et r₂ sont deux racines distinctes du polynôme X − a**X − b
$(λ + μ n) r_{0}^{n}$ si r₀ est racine double du polynôme X − a**X − b
$(λ cos (n θ) + μ sin (n θ)) ρ^{n}$ pour une suite réelle quand ρe et ρe sont les deux racines complexes du polynôme X − a**X − b

On ne perd rien à la généralité de la suite en supposant que celle-ci est définie sur tout $N$ et pas seulement à partir de n₀. En effet, si une suite (u) n’est définie qu’à partir de n₀, elle induit la création d’une suite (v) définie sur $N$ en posant $v_{n} = u_{n + n_{0}}$ .

L’idée est alors de rechercher des suites géométriques vérifiant la récurrence (R). C’est-à-dire chercher des scalaires r tels que la suite $(r^{n})_{n \in N}$ vérifie (R). On démontre aisément que ce problème équivaut à résoudre l’équation du second degré $r^{2} - a r - b = 0$ . Le polynôme $r^{2} - a r - b$ est alors appelé le polynôme caractéristique de la suite. Son discriminant est $Δ = a^{2} + 4 b$ . Il faudra alors distinguer plusieurs cas, selon que le nombre de racine du polynôme caractéristique.

Si le polynôme possède deux racines distinctes

Soient r₁ et r₂ les deux racines distinctes . Les suites $(r_{1}^{n})_{n \in N}$ et $(r_{2}^{n})_{n \in N}$ vérifient (R) ainsi que toute suite dont le terme général serait $λ r_{1}^{n} + μ r_{2}^{n}$ (cela tient au caractère linéaire de la récurrence). A-t-on alors trouvé toutes les suites vérifiant (R) ? Une suite vérifiant (R) étant entièrement déterminée par la donnée de u₀ et u₁, il suffit de prouver que l’on peut toujours trouver λ et μ solutions du système

{λ + μ = u_{0} λ r_{1} + μ r_{2} = u_{1}

Or ce système a pour déterminant r₂ − r₁ non nul. Il est donc toujours possible d’exprimer une suite vérifiant (R) comme combinaison linéaire des suites $(r_{1}^{n})_{n \in N}$ et $(r_{2}^{n})_{n \in N}$

Cette situation se produit pour toute suite à valeurs réelles pour laquelle le discriminant $Δ = a^{2} + 4 b$ est strictement positif, ou pour toute suite à valeurs complexes pour laquelle le discriminant est non nul.

Si le polynôme possède une racine double

Si le discriminant est nul, le problème est tout autre car on ne trouve qu’une seule valeur r₀, donc une seule famille de suites géométriques $(λ r_{0}^{n})_{n \in N}$ vérifiant (R) . L’idée consiste alors à rechercher les suites $(λ_{n})_{n \in N}$ telles que, pour tout entier n, $u_{n} = λ_{n} r_{0}^{n}$ avec $(u_{n})_{n \in N}$ vérifiant (R). Cette méthode s’appelle la méthode de variation de la constante. On s’assure d’abord de l’existence de la suite $(λ_{n})_{n \in N}$ en vérifiant que r₀ n’est jamais nul . La relation de récurrence sur $(u_{n})_{n \in N}$ se traduit par une relation de récurrence sur $(λ_{n})_{n \in N}$ :

r_{0}^{2} λ_{n + 2} = a r_{0} λ_{n + 1} + b λ_{n}

En utilisant ensuite le fait que a + 4b = 0 et que $r_{0} = \frac{a}{2}$ , on obtient la relation caractéristique de toute suite arithmétique :

λ_{n + 2} − λ_{n + 1} = λ_{n + 1} − λ_n

La suite $(λ_{n})_{n \in N}$ est donc une suite arithmétique de terme général

λ_n = λ + μn.

Les suites $(u_{n})_{n \in N}$ vérifiant (R) ont alors pour terme général :

$u_{n} = (λ + μ n) r_{0}^{n}$ .

Ce résultat s'applique pour des suites à valeurs réelles ou complexes pour lesquelles le discrimant du polynôme caractéristique est nul.

Si le polynôme ne possède pas de racine

C'est le cas pour les suites à valeurs réelles pour lesquelles le discriminant du polynôme caractéristique est strictement négatif. L’équation du second degré possède alors dans $C$ deux racines conjuguées.

$r_{1} = ρ e^{i θ}$ et $r_{2} = ρ e^{- i θ}$ .

Les suites de terme général $A ρ^{n} e^{in θ} + B ρ^{n} e^{- in θ}$ sont des suites complexes vérifiant (R). Parmi celles-ci, celles pour lesquelles A et B sont conjugués, sont des suites réelles . Donc les suites de terme général

u_{n} = (λ cos (n θ) + μ sin (n θ)) ρ^{n}

sont des suites réelles vérifiant (R) (on a pris A = λ / 2 − iμ / 2). A-t-on alors trouvé toutes les suites vérifiant (R) ? Une suite vérifiant (R) étant entièrement déterminée par la donnée de u₀ et u₁, il suffit de prouver que l’on peut toujours trouver λ et μ solutions du système

{λ = u_{0} λ ρ cos (θ) + μ ρ sin (θ) = u_{1}

Or ce système a pour déterminant ρsin(θ) non nul. Il est donc toujours possible d’exprimer une suite vérifiant (R) comme combinaison linéaire des suites $(ρ^{n} cos (n θ))_{n \in N}$ et .

Suite récurrente d’ordre p

Sous-espace vectoriel de dimension p

Si on appelle (R_p) la relation de récurrence :

pour tout entier n, $u_{n + p} = a_{0} u_{n} + a_{1} u_{n + 1} + \dots + a_{p - 1} u_{n + p - 1}$

et si on appelle $E_{R_{p}}$ , l’ensemble des suites à valeurs dans K et vérifiant (R_p), on démontre que $E_{R_{p}}$ est un sous-espace vectoriel de l’ensemble des suites à valeurs dans K. Cela tient à la linéarité de la relation de récurrence.

De plus, ce sous espace vectoriel est de dimension p. En effet, il existe un isomorphisme d’espace vectoriel entre $E_{R_{p}}$ et l’ensemble $K^{p}$ : à chaque suite (u) de $E_{R_{p}}$ , on associe le p_uplet $(u_{0}, u_{1}, \dots, u_{p - 1})$ . Il suffit alors de connaître une famille libre de p suites vérifiant (R_p) , l’ensemble $E_{R_{p}}$ est alors engendré par cette famille libre.

Terme général

La recherche du terme général et des suites particulières s’effectue en travaillant sur K . À chaque suite $(u_{n})_{n \in N}$ on associe la suite $(U_{n})_{n \in N}$ telle que

U_{n} = (u_{n}, u_{n + 1}, \dots, u_{n + p - 1})

La relation de récurrence sur $(u_{n})_{n \in N}$ induit une relation de récurrence sur $(U_{n})_{n \in N}$

U_{n + 1} = A**U_n où

A = (010 \dots 0001 \dots 0 ⋮ ⋱ ⋱ \dots ⋮ 0 \dots \dots 0 1 a_{0} a_{1} \dots \dots a_{p - 1})

Le terme général de la suite U est alors déterminé par

U_n = A**U₀

Le problème semble alors terminé. Mais la réelle difficulté consiste alors à calculer A... On préfère plutôt déterminer une base de $E_{R_{p}}$ .

Recherche d'une base

Le polynôme caractéristique de la matrice A est $P (X) = X^{p} - \sum_{i = 0}^{p - 1} a_{i} X^{i}$ . Ce n'est pas un hasard si on le retrouve pour caractériser les suites $u = (u_{n})_{n \in N}$ vérifiant R_p.

On note f la transformation linéaire qui, à une suite $u = (u_{n})_{n \in N}$ associe la suite $v = (v_{n})_{n \in N}$ définie par v_n = u_{n + 1}. La condition u vérifie R_p se traduit alors par P(f)(u) = 0. L'ensemble $E_{R_{p}}$ est donc le noyau de P(f). Si P est un polynôme scindé dans K (ce qui est toujours vrai si $K = C$ ), il existe k racines $r_{1}, r_{2}, \dots, r_{k}$ et k exposants $α_{1}, α_{2}, \dots, α_{k}$ tel que $P = \prod_{i = 1}^{k} (X - r_{i})^{α_{i}}$ . Le noyau de P(f) est alors la somme directe des noyaux des $(f - r_{i} I d)^{α_{i}}$ . Il suffit donc de trouver une base de chacun de ces noyaux pour déterminer une base de $E_{R_{p}}$ .

On peut montrer que toute suite de terme général $Q (n) r_{i}^{n}$ est élément du noyau de $(f - r_{i} I d)^{α_{i}}$ pour peu que le degré de Q soit inférieur strictement à α_i. Cette démonstration se fait par récurrence sur α_i. Comme les suites $(n^{j} r_{i}^{n})_{n \in N}$ , pour j = 0 à α_i − 1 forment une partie libre de α_i éléments, la famille de toutes les suites $(n^{j} r_{i}^{n})_{n \in N}$ , pour j = 0 à α_i − 1 et pour i = 1 à k forme une famille libre de $α_{1} + α_{2} + \dots + α_{k} = p$ éléments de $E_{R_{p}}$ (de dimension p) donc une base de $E_{R_{p}}$ . Les éléments de $E_{R_{p}}$ sont donc des sommes de suites dont le terme général est $Q (n) r_{i}^{n}$ avec degré de Q strictement inférieur à α_i.

Retour à la récurrence d'ordre 2

Si le polynôme caractéristique se scinde en (X − r₁)(X − r₂) alors les polynômes Q sont de degré 0 et les éléments de $E_{R_{2}}$ sont des suites dont le terme général est $λ_{1} r_{1}^{n} + λ_{2} r_{2}^{n}$ .

Si le polynôme caractéristique se scinde en (X − r₀) alors les polynômes Q sont de degré 1 et les éléments de $E_{R_{2}}$ sont des suites dont le terme général est $(λ_{1} n + λ_{2}) r_{0}^{n}$ .