Exercices sur les suites

Retour à la page personnelle de Bernard Parisse.

# Chapitre 26 Exercices sur les suites

La commande limite de Xcas donne la limite d’une suite losqu’elle existe !
La commande seqsolve de Xcas donne la valeur d’une suite récurrente de la forme un+1=f(un) ou un+2=f(,un,un+1,...) ou d’un système de suites récurrentes.
Par exemple :
Pour trouver la limite quand n tend vers +∞ de la suite un définie par :
un=∑p=1n1/n+p, on tape :
limite(sum(1/((n+p)),p=1..n),n,inf)
On obtient :
ln(2)
Pour trouver la valeur de la suite un définie par :
un+2=un+2un+1+n+1 u0=0, u1=1, on tape :
seqsolve(x+2*y+n+1,[x,y,n],[0,1])
On obtient :
(-4*n-(sqrt(2)+1)`^`n*(-sqrt(2)*3-2)-(-(sqrt(2))+1)`^`n*(sqrt(2)*3-2)-4)/8
Pour trouver la valeur des suite un et vn définies par :
un+1=un+2vn, vn+1=un+n+1 u0=0, v0=1, on tape :
seqsolve([x+2*y,n+1+x],[x,y,n],[0,1])
On obtient :
[(-2*n-(-1)`^`n+2`^`n*4-3)/2,((-1)`^`n+2*2`^`n-1)/2]
Exercices

1. Soient an, et bn les suites définies par :
   a0, b0 ((a0,b0)≠ (0,0)),
   an+1=an/an2+bn2,
   bn+1=bn/an2+bn2,
   Calculer les 4 premiers termes de ces suites pour a0=a, b0=b,
   On utilise le tableur pour avoir les valeurs des termes de ces suites.
   On tape pour :
   A0 1 (ou a)
   B0 3 (ou b)
   A1 =simplfy(A0/(A0`^`2+B0`^`2))
   B1 =simplfy(B0/(A0`^`2+B0`^`2))
   Puis on remplit vers le bas. On obtient commes valeurs dans la colonne A :
   a, a/(a`^`2+b`^`2), a, a/(a`^`2+b`^`2) ...
   On obtient commes valeurs dans la colonne B :
   b, b/(a`^`2+b`^`2), b, b/(a`^`2+b`^`2) ...
   Soit la fonction de ℝ2 dans ℝ2 :
   g(x,y)=(x/(x2+y2),y/(x2+y2)) Montrons avec Xcas que g∘ g=id
   On tape :
   g(x,y):=(x/(x`^`2+y`^`2),y/(x`^`2+y`^`2))
   simplify((g@g)(x,y))
   On obtient :
   x,y
   La fonction g transforme un vecteur V de ℝ2 en W=V/norm(V)2. On a norm(W)=1/norm(V) donc g transforme W en W/norm(W)2=V.
2. Soient an, bn, et cn les suites définies par :
   a0=1, b0=3 et c0=8,
   an+1=1/2(bn+cn),
   bn+1=1/2(cn+an),
   cn+1=1/2(an+bn),
   Calculer les 4 premiers termes de ces suites,
   Montrer que an+bn+cn est constant.
   Exprimer an+1−bn+1 en fonction de an−bn,
   Montrer que pour n≥ 1, 2an+1=an+an−1.
   En déduire une expression de an, bn, et cn en fonction de n.
   Étudier la convergence de an, bn, et cn.
   On peut faire le calcul des 4 premiers termes avec le tableur, on obtient : A:=(1,11/2,13/4,35/8,61/16) pour les premiers termes de an,
   B:=(3,9/2,15/4,33/8,63/16) pour les premiers termes de bn,
   C:=(8,2,5,7/2,17/4) pour les premiers termes de cn,
   On tape :
   [A]+[B]+[C]
   On obtient :
   [12,12,12,12,12]
   On a :
   an+1+bn+1+cn+1=1/2(bn+cn+cn+an+an+bn)=an+bn+cn Donc on montre par récurrence que an+bn+cn=12 On a :
   an+1−bn+1=1/2(bn−an)=−1/2(an−bn)
   2an+1=bn+cn=an−1+1/2(bn−1+cn−1)=an−1+an
   Donc on peut définir an par :
   a0=1,a1=11/2 et la relation de récurrence 2an+1=an+an−1.
   L’ensemble E des suites réelles u qui vérifient la relation de récurrence :
   2un+1=un+un−1 pour n>0
   forment un espace vectoriel sur ℝ de dimension 2.
   En effet E est un espace vectoriel sur ℝ et l’application :
   f de E dans ℝ2 qui à u∈ E fait correspondre (u0,u1) est un isomorphisme d’espaces vectoriels.
   On cherche les progressions géométriques un=u0rn qui vérifient :
   2un+1=un+un−1 ont pour raison r=1 ou r=−1/2 donc
   an=α+β(−1/2)n avec :
   α+β=1 et α−β/2=11/2
   On tape :
   solve([x+y=1,x-y/2=11/2],[x,y])
   On obtient :
   [[4,-3]]
   Donc an=4−3(−1/2)n De même :
   bn=α+β(−1/2)n avec :
   α+β=3 et α−β/2=19/2
   On tape :
   solve([x+y=3,x-y/2=9/2],[x,y])
   On obtient :
   [[4,-1]]
   Donc bn=4−(−1/2)n et,
   cn=α+β(−1/2)n avec :
   α+β=8 et α−β/2=2
   On tape :
   solve([x+y=8,x-y/2=2],[x,y])
   On obtient :
   [[4,4]]
   Donc cn=4+4(−1/2)n On vérifie, on tape :
   seqsolve([(b+c)/2,(c+a)/2,(a+b)/2],[a,b,c,n],[1,3,8])
   On obtient :
   [-3*(-1/2)`^`n+4,-(-1/2)`^`n+4,4*(-1/2)`^`n+4]
   Quand n tend vers +∞ on a (−1/2)n tend vers 0 donc :
   an tend vers 4 quand n tend vers +∞
   bn tend vers 4 quand n tend vers +∞
   cn tend vers 4 quand n tend vers +∞ Prolongement : cas général Soient an, bn, et cn les suites définies par :
   a0, b0 et c0,
   an+1=1/2(bn+cn),
   bn+1=1/2(cn+an),
   cn+1=1/2(an+bn),
   On a toujours :
   an+1+bn+1+cn+1=1/2(bn+cn+cn+an+an+bn)=an+bn+cn Donc :
   an+bn+cn=a0+b0+c0 On a toujours :
   an+1−bn+1=1/2(bn−an)=−1/2(an−bn)
   2an+1=bn+cn=an−1+1/2(bn−1+cn−1)=an−1+an
   Donc an est définie par :
   2an+1=an−1+an et a0 et a1=1/2(b0+c0)
   Donc :
   an=α+β(−1/2)n avec :
   α+β=a0 et α−β/2=1/2(b0+c0)
   On tape :
   solve([x+y=a0,x-y/2=1/2*(b0+c0)],[x,y])
   On obtient :
   [[a0/3+b0/3+c0/3,2*a0/3-b0/3-c0/3]]
   Donc :
   an=−1/3(a0+b0+c0)+(2a0−b0−c0)(−1/2)n bn=−1/3(a0+b0+c0)+(2b0−c0−ca0)(−1/2)n cn=−1/3(a0+b0+c0)+(2c0−a0−b0)(−1/2)n On vérifie, on tape :
   seqsolve([(b+c)/2,(c+a)/2,(a+b)/2],[a,b,c,n],[a0,b0,c0])
   On obtient :
   [(2*(-1/2)`^`n*a0+a0-(-1/2)`^`n*b0+b0-(-1/2)`^`n*c0+c0)/3, (-(-1/2)`^`n*a0+a0+2*(-1/2)`^`n*b0+b0-(-1/2)`^`n*c0+c0)/3, (-(-1/2)`^`n*a0+a0-(-1/2)`^`n*b0+b0+2*(-1/2)`^`n*c0+c0)/3]
   On applique la commande factor et on obtient :
   [((2*a0-b0-c0)*(-1/2)`^`n+a0+b0+c0)/3, ((-a0+2*b0-c0)*(-1/2)`^`n+a0+b0+c0)/3, ((-a0-b0+2*c0)*(-1/2)`^`n+a0+b0+c0)/3]
3. Soient an, bn, et cn les suites définies par :
   a0=1, b0=0 et c0=0,
   an+1=bn+cn,
   bn+1=cn+an,
   cn+1=an+bn,
   Calculer les 4 premiers termes de ces suites,
   Montrer que pour tout n bn=cn
   Calculer an+bn+cn.
   On a :
   an+1+bn+1+cn+1=(bn+cn+cn+an+an+bn)=2(an+bn+cn)
   donc:
   a1+b1+c1=2(a0+b0+c0)
   a2+b2+c2=2(a1+b1+c1)=22(a0+b0+c0)
   etc... an+bn+cn=2n(a0+b0+c0) On en déduit que :
   an=2n bn+1=bn+2n donc
   bn=2n−1+bn−1=2n−1+2n−2+bn−2=..=∑k=0n−12k= Prolongement : cas général Soient an, bn, et cn les suites définies par :
   a0, b0 et c0,
   an+1=bn+cn,
   bn+1=cn+an,
   cn+1=an+bn,
   On a toujours :
   an+1+bn+1+cn+1=bn+cn+cn+an+an+bn=2(an+bn+cn)
   Donc :
   an+bn+cn=2n(a0+b0+c0) On n’a plus bn=cn mais an+1−bn+1=bn−an=−(an−bn),
   Donc :
   an−bn=(−1)n(a0−b0),
   De même :
   an−cn=(−1)n(a0−c0),
   On a :
   an=1/3((an+bn+cn)+an−bn+an−cn)
   an=1/3(2n(a0+b0+c0)+(−1)n(a0−b0)+(−1)n(a0−c0)
   Donc :
   an=1/3(2n(a0+b0+c0)+(−1)n(2a0−b0−c0)
   De même :
   bn=1/3(2n(a0+b0+c0)+(−1)n(2b0−c0−a0)
   cn=1/3(2n(a0+b0+c0)+(−1)n(2c0−b0−a0)
   On vérifie, on tape :
   seqsolve([b+c,c+a,a+b],[a,b,c,n],[a0,b0,c0])
   On obtient après avoir fait agir la commande factor :
   [((-b0-c0+2*a0)*(-1)`^`n+2`^`n*(b0+c0+a0))/3, ((-a0-c0+2*b0)*(-1)`^`n+2`^`n*(a0+c0+b0))/3,((-a0-b0+2*c0)*(-1)`^`n+2`^`n(a0+b0+c0))/3]
4. Explicitez an, bn, et cn les suites définies par :
   a0, b0 et c0,
   an+1=an+cn,
   bn+1=bn+an,
   cn+1=cn+bn,
   On a toujours :
   an+1+bn+1+cn+1=an+cn+bn+an+cn+bn=2(an+bn+cn)
   Donc :
   an+bn+cn=2n(a0+b0+c0) On a :
   an+1−cn+1+bn+1=2an bn+1−an+1+cn+1=2bn cn+1−bn+1+an+1=2cn Malheureusement, on a moins de chances que dans l’exercice précédent...
   On va donc utiliser une solution matricielle.
   Soit :
   A:=[[1,0,1],[1,1,0],[0,1,1]] On a :
   [an,bn,cn]=A`^`n*[a0,b0,c0]
   Il faut donc calculer An.
   On tape :
   P,B:=jordan(A)
   On obtient :
   [[1,sqrt(3)*(-i)-1,sqrt(3)*(i)-1],[1,sqrt(3)*(i)-1,sqrt(3)*(-i)-1],[1,2,2]],[[2,0,0],[0,(sqrt(3)*(i)+1)/2,0],[0,0,(sqrt(3)*(-i)+1)/2]]
   On tape :
   P1:=simplify(inv(P))
   On obtient :
   [[1/3,1/3,1/3],[((i)*sqrt(3)-1)/12,((-i)*sqrt(3)-1)/12,1/6],[((-i)*sqrt(3)-1)/12,((i)*sqrt(3)-1)/12,1/6]] On tape :
   An:=factor(simplify(exp2trig(pow2exp(P*Bn*P1))))
   On obtient :
   [[(2`^`n+2*cos(n*pi/3))/3, (2`^`n-cos(n*pi/3)-sqrt(3)*sin(n*pi/3))/3, (2`^`(n+1)-2*cos(n*pi/3)+sin(n*pi/3)*2*sqrt(3))/6],
   [(2`^`n-cos(n*pi/3)+sqrt(3)*sin(n*pi/3))/3, (2`^`n+2*cos(n*pi/3))/3, (2`^`(n+1)-2*cos(n*pi/3)-2*sqrt(3)*sin(n*pi/3))/6],
   [(2`^`n-cos(n*pi/3)-sqrt(3)*sin(n*pi/3))/3, (2`^`n-cos(n*pi/3)+sqrt(3)*sin(n*pi/3))/3, (2`^`(n+1)+4*cos(n*pi/3))/6]]
   On tape :
   On obtient :
   On tape en mode complexe :
   simplify(exp2trig(pow2exp(seqsolve([x+z,y+x,z+y],[x,y,z,n],[a0,b0,c0]))))
   On obtient après avoir utilisé factor:
   ((a0+b0+c0)*2`^`n+cos(n*pi/3)*(2*a0-b0-c0)+sqrt(3)*sin(n*pi/3)*(-b0+c0))/3,
   ((a0+b0+c0)*2`^`n-cos(n*pi/3)*(a0-2*b0+c0)+sqrt(3)*(a0-c0)*sin(n*pi/3))/3
   ((a0+b0+c0)*2`^`n-cos(n*pi/3)*(a0+b0-2*c0)+sqrt(3)*(-a0+b0)*sin(n*pi/3))/3
   Si on pose :
   j=(sqrt(3)*(i)−1)/2 on a j2=(−sqrt(3)*(i)−1)/2.
   Donc :
   B=[[2,0,0],[0,−j2,0],[0,0,−j]]
   P=[[1,2j2,2j],[1,2j,2j2],[1,2,2]]
   P1=P−1
   On peux faire les calculs à la main.... :
   B:=[[2,0,0],[0,-j`^`2,0],[0,0,-j]]
   Bn:=matpow(B,n)
   Bn:=[[2`^`n,0,0],[0,(-j`^`2)`^`n,0],[0,0,(-j)`^`n]]
   P:=[[1,2j`^`2,2j],[1,2j,2j`^`2],[1,2,2]]
   P1:=[[1/3,1/3,1/3],[j/6,j`^`2/6,1/6],[j`^`2/6,j/6,1/6]]
   An:=P*Bn*P1
   On obtient pour An après simplification à la main :
   [[u(n),w(n),v(n)],[v(n),u(n),w(n)],[w(n),v(n),u(n)]]
   avec :
   u(n):=1/3*(2`^`n+(-j)`^`n+(-j`^`2)`^`n])
   v(n):=1/3*(2`^`n+j*(-j)`^`n+j`^`2*(-j`^`2)`^`n)
   w(n):=1/3*(2`^`n+j`^`2*(-j)`^`n+j*(-j`^`2)`^`n)
   Donc on a :
   an=a0*u(n)+b0*w(n)+c0*v(n)
   bn=a0*v(n)+b0*u(n)+c0*w(n)
   cn=a0*w(n)+b0*v(n)+c0*u(n)

Prolongements
Soit E un ensemble fini de cardinal n.
Soit an le nombre de parties de E de cardinal congru à 0mod3,
Soit bn le nombre de parties de E de cardinal congru à 1mod3,
Soit cn le nombre de parties de E de cardinal congru à 2mod3,

1. Montrer que :
   an+bn+cn=2n.
2. Montrer que :
   an=1+comb(n,3)+comb(n,6)+... bn=comb(n,1)+comb(n,4)+comb(n,7)+... cn=1+comb(n,2)+comb(n,5)+comb(n,8)+...
3. En développant (1+X)n successivement pour X=1, X=j, X=j2, calculer an, bn, cn.
4. Donner la valeur de a0, b0 et c0. Calculer:
   an+1 en fonction de an et cn,
   bn+1 en fonction de bn et an,
   cn+1 en fonction de cn et bn Et faire le lien avec l’exercice précédent.

1. Soit n=3q+r avec 0≤ r<3.
   Le nombre de parties de E est 2n donc an+bn+cn=2n.
2. Les parties de E de cardinal congru à 0mod3 sont de cardinal p avec p∈ (0,3,6..3q) Les parties de E de cardinal congru à 1mod3 sont de cardinal p avec p∈ (1,4,7..3q−2) si r=0 ou p∈ (1,4,7..3q+1) si 0<r<3 Les parties de E de cardinal congru à 2mod3 sont de cardinal p avec p∈ (2,5,8..3q−1) si r=0 ou r=1 ou p∈ (2,5,8..3q+2) si r=2 Donc an=1+comb(n,3)+comb(n,6)+...comb(n,3q) bn=comb(n,1)+comb(n,4)+comb(n,7)+...+comb(n,3q−2)+comb(n,3q+1) cn=1+comb(n,2)+comb(n,5)+comb(n,8)+...+comb(n,3q−1)+comb(n,3q+2)
3. On a :
   (1+1)n)=1+comb(n,1)+comb(n,2)+comb(n,3)+....+comb(n,n)
   (1+j)n)=1+j*comb(n,1)+j2*comb(n,2)+comb(n,3)+....+jn*comb(n,n)
   (1+j2)n)=1+j2*comb(n,1)+j*comb(n,2)+comb(n,3)+....+j2ncomb(n,n)
   On additionne ces égalités après les avoir mulltiplié successivement par :
   (1,1,1) puis par (1,j2,j) puis par (1,j,j2). On obtient (puisque j3=1 et 1+j+j2=0) :
   3an=2n+(1+j)n+(1+j2)n
   3bn=2n+j2(1+j)n+j(1+j2)n
   3cn=2n+j(1+j)n+j2(1+j2)n
4. On a :
   a0=1, b0=c0=0
   Quand on rajoute un élément α à E :
   parmi les sous-ensembles de cardinal congru à 0mod3, il y en a an qui ne contiennent pas α et cn qui contiennent α donc
   an+1=an+cn,
   parmi les sous-ensembles de cardinal congru à 1mod3 il y en a bn qui ne contiennent pas α et an qui contiennent α donc
   bn+1=bn+an,
   parmi les sous-ensembles de cardinal congru à 2mod3 il y en a cn qui ne contiennent pas α et bn qui contiennent α donc
   cn+1=cn+bn,
   Donc an, bn, cn sont définies par :
   a0=1, b0=c0=0
   an+1=an+cn,
   bn+1=bn+an,
   cn+1=cn+bn
   On tape :
   simplify(exp2trig(pow2exp(seqsolve([x+z,y+x,z+y],[x,y,z,n],[1,0,0]))))
   On obtient :
   [(2`^`n+2*cos(n*pi/3))/3,
   (2`^`n-cos(n*pi/3)+(sqrt(3))*sin(n*pi/3))/3),
   (2`^`n-cos(n*pi/3)+(-(sqrt(3))*sin(n*pi/3))/3]
5. On tape :
   On obtient :
6. On tape :
   On obtient :

## 26.1 Exercices sur les séries

La commande sum de Xcas calcule les sommes de certaines séries.
Voici queques exemples :
On tape :
sum(1/2`^`n,n=0..inf)
On obtient :
2
On tape :
sum(n/2`^`n,n=0..inf)
On obtient :
2
On tape :
sum(1/(n*(n+1)),n=1..inf)
On obtient :
1
On tape :
sum(1/((2n+1)*(2n-1)),n=1..inf)
On obtient :
1/2
On tape :
sum(1/n`^`2,n=1..inf)
On obtient :
pi`^`2/6
On tape :
sum((-1)`^`(n+1)/n,n=1..inf)
On obtient :
ln(2)
On tape :
sum((-1)`^`n/(2n+1),n=0..inf)
On obtient :
pi/4
Exercices

1. Montrer que :
   (1−x)*∑k=0nxk=1−xn+1 En déduire que :

   ∑k=pnxk= |
   1−xn+1
   1−x

   Montrer que :

   ∑k=pinfxk=xp |
   1
   1−x

   Application :
   Calculer : ∑k=0inf2−k,∑k=0inf3−k
2. Soit un=n/2n. Monter que :
   4un+1<3un lorsque n>2
   En déduire que un tend vers 0 quand n tend vers +∞.
3. pour ceux qui connaissent le critére de d’Alembert sur les séries à termes positifs,
   Calculer un+1/un (pour ceux qui connaissent le critére de d’Alembert sur les séries à termes positifs).
   En déduire que la série de terme général un est convergente.
4. pour ceux qui ne connaissent pas le critére de d’Alembert,
   Montrer que limN−>∞∑k=0Nun existe et calculer sa valeur.
5. Calculer la somme ∑k=0Nun.

On tape :
u(n):=n/2`^`n
factor(simplify(4u(n+1)-3u(n)))
On obtient :
(-n+2)/2`^`n
On a donc :
un+1/un<3/4
un+1/un=n+1/2n
Donc un+1/un tend vers 1/2 quand n tend vers +∞.
On tape :
limit(u(n+1)/u(n),n,inf)
On obtient :
1/2
Calcul de SN=∑n=0N un
SN=0+∑n=1N 1/2n+∑n=2N 1/2n+...∑n=kN 1/2n+..∑n=NN 1/2n.
Donc
SN=1+1/2N+1/2+1/2N+..1/2k+1/2N+1/2N
SN=∑k=1N1/2k+N/2N
Comme N/2N tend vers 0 quand N tend vers +∞, on en déduit que :
SN tend vers ∑k=0∞1/2k=2.

Retour à la page personnelle de Bernard Parisse.