Suite récurrente linéaire
# Suite récurrente linéaire

En mathématiques, on appelle suite récurrente linéaire d’ordre p toute suite à valeurs dans un corps commutatif K (par exemple ℝ ou ℂ ; on ne se placera que dans ce cas dans cet article) définie pour tout $n\geq n_{0}$ par une relation de récurrence linéaire de la forme

$\forall n\geq n_{0}\quad u_{n+p}=a_{0}u_{n}+a_{1}u_{n+1}+\cdots +a_{p-1}u_{n+p-1}$

où $a_{0}$, $a_{1}$, …$a_{p-1}$ sont p scalaires fixés de K ($a_{0}$ non nul).

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

Les suites récurrentes linéaires d’ordre 1 sont les suites géométriques.

L'étude des suites récurrentes linéaires d'ordre supérieur se ramène à un problème d'algèbre linéaire. L'expression du terme général d'une telle suite est possible pour peu qu'on soit capable de factoriser un polynôme qui lui est associé, appelé polynôme caractéristique ; le polynôme caractéristique associé à une suite vérifiant la relation de récurrence ci-dessus est :

    $P(X)=X^{p}-\sum _{i=0}^{p-1}a_{i}X^{i}=X^{p}-a_{p-1}X^{p-1}-a_{p-2}X^{p-2}-\dots -a_{1}X-a_{0}.$

Son degré est ainsi égal à l'ordre de la relation de récurrence. En particulier, dans le cas des suites d'ordre 2, le polynôme est de degré 2 et peut donc être factorisé à l'aide d'un calcul de discriminant. Ainsi, le terme général des suites récurrentes linéaires d'ordre 2 peut être exprimé en utilisant seulement les deux premiers termes, quelques valeurs constantes, quelques opérations élémentaires de l'arithmétique (addition, soustraction, multiplication, exponentielle) et les fonctions sinus et cosinus (si le corps des scalaires est le corps des réels). Une des suites de ce type est la célèbre suite de Fibonacci, qui peut s'exprimer à partir de puissances faisant intervenir le nombre d'or.

## Suite récurrente linéaire d’ordre 1

Les suites récurrentes linéaires d'ordre 1 sont les suites géométriques.

Si la relation de récurrence est $u_{n+1}=qu_{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}=au_{n+1}+bu_{n}\quad (R).$

Les scalaires r tels que la suite $(r^{n})_{n\in \mathbb {N} }$ vérifie (R) sont les solutions de l’équation du second degré $r^{2}-ar-b=0$. Le polynôme $X^{2}-aX-b$ est alors appelé le polynôme caractéristique de la suite. Son discriminant est $\Delta =a^{2}+4b$. Il faudra alors distinguer plusieurs cas, selon le nombre de racines du polynôme caractéristique.

Théorème —  Le terme général d'une suite à valeurs dans K et vérifiant (R) est :

* $\lambda r_{1}^{n}+\mu r_{2}^{n}$ si $r_{1}$ et $r_{2}$ sont deux racines distinctes (dans K) du polynôme $X^{2}-aX-b$,
* $(\lambda +\mu n)r_{0}^{n}$ si $r_{0}$ est racine double du polynôme $X^{2}-aX-b$,

avec $\lambda ,\mu$ paramètres dans K déterminés par les deux premières valeurs de la suite.

Le cas 1 se produit par exemple si $K=\mathbb {R}$ et si le discriminant $\Delta =a^{2}+4b$ est strictement positif, ou si $K=\mathbb {C}$ et $\Delta \neq 0$. De plus, si les deux racines $r_{1},r_{2}$ du polynôme $X^{2}-aX-b$ sont deux complexes conjugués $\rho \mathrm {e} ^{\mathrm {i} \theta }$ et $\rho \mathrm {e} ^{-\mathrm {i} \theta }$, alors le terme général d'une telle suite s'écrit également :

* $\rho ^{n}(A\cos(n\theta )+B\sin(n\theta ))$ avec A et B paramètres dans K ($=\mathbb {R}$ ou $\mathbb {C}$, selon qu'on s'intéresse aux suites réelles ou complexes) déterminés par les deux premières valeurs de la suite.

Le cas 2 se produit lorsque $\Delta =0$ et alors, la racine double est $r_{0}={\frac {a}{2}}$.

On ne perd rien à la généralité de la suite en supposant que celle-ci est définie sur tout ℕ et pas seulement à partir de $n_{0}$. En effet, l'étude d'une suite u qui n’est définie qu’à partir de $n_{0}$ se ramène à celle de la suite v définie sur ℕ par $v_{n}=u_{n+n_{0}}$.

### Identités remarquables

Si une suite u vérifie

    $\forall n\in \mathbb {N} \quad u_{n+2}=au_{n+1}+bu_{n}\quad (R)$

alors, elle peut être étendue aux indices négatifs et reliée aux puissances de la matrice (appelée matrice compagnon du polynôme caractéristique)

    $P:={\begin{pmatrix}a&b\\1&0\end{pmatrix}}$

(inversible puisque b ≠ 0) par :

    $\forall n\in \mathbb {Z} \quad {\begin{pmatrix}u_{n+1}\\u_{n}\end{pmatrix}}=P^{n}{\begin{pmatrix}u_{1}\\u_{0}\end{pmatrix}}$.

Ceci permet de montrer que pour v égale à u ou à toute autre suite vérifiant la même relation de récurrence (R) et pour tous entiers i, j, k, l et $r^,$ :

$i+j=k+l\Rightarrow u_{i+r}v_{j+r}-u_{k+r}v_{l+r}=(-b)^{r}(u_{i}v_{j}-u_{k}v_{l})$

En particulier :

${\text{si }}u_{0}=0{\text{ alors}}\quad \forall m,n,r\in \mathbb {Z} \quad u_{n}v_{m+r}-u_{r}v_{m+n}=(-b)^{r}u_{n-r}v_{m}$

## Suite récurrente d’ordre p

### Sous-espace vectoriel de dimension p

Si l'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}+\cdots +a_{p-1}u_{n+p-1}$

et si l'on note $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'espace des suites à valeurs dans K. Cela tient à la linéarité de la relation de récurrence.

De plus, ce sous-espace est de dimension p. En effet, il existe un isomorphisme d'espaces vectoriels entre $E_{R_{p}}$ et $K^{p}$ : à chaque suite u de $E_{R_{p}}$, on associe le p-uplet $(u_{0},u_{1},\cdots ,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^{p}$. À chaque suite $(u_{n})_{n\in \mathbb {N} }$ on associe la suite $(U_{n})_{n\in \mathbb {N} }$ définie par

    $U_{n}={\begin{pmatrix}u_{n}\\u_{n+1}\\\vdots \\u_{n+p-1}\end{pmatrix}}.$

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

    $U_{n+1}=AU_{n}$ où
    $A={\begin{pmatrix}0&1&0&\cdots &0\\0&0&1&\cdots &0\\\vdots &\ddots &\ddots &\ddots &\vdots \\0&\cdots &\cdots &0&1\\a_{0}&a_{1}&\cdots &\cdots &a_{p-1}\end{pmatrix}}$

(A est la matrice compagnon du polynôme caractéristique de la suite).

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

    $U_{n}=A^{n}U_{0}.$

Le problème semble alors terminé. Mais la réelle difficulté consiste alors à calculer $A^{n}$… On préfère 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 \mathbb {N} }$ vérifiant $R_{p}$.

On note f la transformation linéaire qui, à une suite $u=(u_{n})_{n\in \mathbb {N} }$ associe la suite $v=(v_{n})_{n\in \mathbb {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 le polynôme P est scindé sur K (ce qui est toujours vrai si K = ℂ), il s'écrit $P=\prod _{i=1}^{k}(X-r_{i})^{\alpha _{i}}$, où $r_{1},r_{2},\dots ,r_{k}$ sont les racines de P et $\alpha _{1},\alpha _{2},\dots ,\alpha _{k}$ leurs ordres de multiplicité respectifs. Le noyau de P(f) est alors la somme directe des noyaux des $(f-r_{i}\mathrm {Id} )^{\alpha _{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}\mathrm {Id} )^{\alpha _{i}}$ pour peu que le degré de Q soit strictement inférieur à $\alpha _{i}$. Cette démonstration se fait par récurrence sur $\alpha _{i}$. Comme les suites $(n^{j}r_{i}^{n})_{n\in \mathbb {N} }$, pour j = 0 à $\alpha _{i}-1$, forment une partie libre de $\alpha _{i}$ éléments, les suites $(n^{j}r_{i}^{n})_{n\in \mathbb {N} }$, pour j de 0 à $\alpha _{i}-1$ et i de 1 à k, forment une famille libre de $\alpha _{1}+\alpha _{2}+\cdots +\alpha _{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 à $\alpha _{i}$.

### Retour à la récurrence d'ordre 2

Si le polynôme caractéristique se scinde en $(X-r_{1})(X-r_{2})$ 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 $\lambda _{1}r_{1}^{n}+\lambda _{2}r_{2}^{n}$.

Si le polynôme caractéristique se scinde en $(X-r_{0})^{2}$ 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 $(\lambda _{1}n+\lambda _{2})r_{0}^{n}$.

