Matrice - calcul de puissance de A

Réponses à toutes vos questions après le Bac (Fac, Prépa, etc.)
Azuriel
Membre Relatif
Messages: 236
Enregistré le: 28 Déc 2006, 16:48

Matrice - calcul de puissance de A

par Azuriel » 20 Fév 2007, 16:30

Voila mon probleme,
J'ai une matrice A defini par les membres de sa diagonales kisont égaux à b chacun, et sinon pour tout les autres Ai,j , les Ai,j sont égaux à a.

En gros ma matrice à cette tronche (à part que dans l'exo elle appartient à Mn(k))

b a a a
a b a a
a a b a
a a a b etc

On me demande de calculer A puissance m avec m entier naturel mais je ne vois vraiment pas comment faire..Je sais pas si decomposer cette matrice en une identité et l'autre que avec les a peut me servir...

Merci d'avance !



BQss
Membre Irrationnel
Messages: 1202
Enregistré le: 02 Nov 2006, 03:32

par BQss » 20 Fév 2007, 16:39

Azuriel a écrit:Voila mon probleme,
J'ai une matrice A defini par les membres de sa diagonales kisont égaux à b chacun, et sinon pour tout les autres Ai,j , les Ai,j sont égaux à a.

En gros ma matrice à cette tronche (à part que dans l'exo elle appartient à Mn(k))

b a a a
a b a a
a a b a
a a a b etc

On me demande de calculer A puissance m avec m entier naturel mais je ne vois vraiment pas comment faire..Je sais pas si decomposer cette matrice en une identité et l'autre que avec les a peut me servir...

Merci d'avance !



Toute matrice commute avec l'dentité.
decompose ca en (bI + A)^n grace au binome de Newton...

Ps: et tu peux meme decomposer la matrice de diagonale nulle A en B +tB avec B une matrice triangulaire superieure de "a".

fahr451
Membre Transcendant
Messages: 5142
Enregistré le: 05 Déc 2006, 23:50

par fahr451 » 20 Fév 2007, 20:12

bonsoir

dans le même esprit mais plus immédiat

A = (b-a) I + a J

où J est la matrice pleine de 1 , J^2 = nJ et J^p = n^(p-1) J pour p>0

puis formule du binôme

Azuriel
Membre Relatif
Messages: 236
Enregistré le: 28 Déc 2006, 16:48

par Azuriel » 21 Fév 2007, 14:25

Merci beaucoup, j'étais pas sur qu'une decomposition pouvaient servir mais dans ce cas particulier je me disais que peut etre oui. Merci en tout cas.
D'ailleurs j'ai un autre exo qui me pose probleme sur mon TD de matrice, mais que j'avais delaissé mais j'en profite pour vous le sortir au cas où, je ne crois pas que ça soit difficile mais c'est juste que cet exo me derange..je ne dois pas vraiment bien partir. Le voici :

Soit E = {l'ensemble des matrice A apartenant à Mn(K) tel que tout Ai,j > 0 et pour TOUT i , Sum(j=1 à n)Ai,j = 1 }

Montrer que E est stable par multiplication
Determiné les matrices inversibles de E dont l'inverse est egalement dans E.

Merci d'avance !

fahr451
Membre Transcendant
Messages: 5142
Enregistré le: 05 Déc 2006, 23:50

par fahr451 » 21 Fév 2007, 14:36

pour le début il suffit d écrire le coefficient cij du produit AB

cij = sigma aikbkj avec A et B dans E les aik bkj >0 donc cij >0
et
sigma(sur j) cij = sigma( sur j et k ) aik bkj =
sigma (sur k ) aik sigma (surj) bkj = sigma (sur k) aik = 1

donc C est dans E

Azuriel
Membre Relatif
Messages: 236
Enregistré le: 28 Déc 2006, 16:48

par Azuriel » 21 Fév 2007, 15:02

Merci beaucoup, ce qui me posait probleme était surtout pour le fait que la somme soit bien égal à 1.

Sinon pour l'exercice précedent, tu arrives a quoi comme résultat apres avoir calculer, simplifier le binome puis remis sous une forme..
Moi j'ai ça mais ça me parait bizarre :

A^m = [a((b-a) + an)^m-1 ]xJ + [(b-a)^m]/n x J

:hein:

fahr451
Membre Transcendant
Messages: 5142
Enregistré le: 05 Déc 2006, 23:50

par fahr451 » 21 Fév 2007, 15:36

il y a forcément un terme en I qui vaut (b-a)^m I

Azuriel
Membre Relatif
Messages: 236
Enregistré le: 28 Déc 2006, 16:48

par Azuriel » 21 Fév 2007, 15:44

Ba oui c'est celui que j'ai sorti qui est divisé par n..
Je te detailles le calcul et tu me dis là où jai fai une erreur..

A^m = Sum(k=0..m)(k parmi m) [(b-a)I]^k x [aJ]^m-k
A^m = Sum(k=0..m)(k parmi m) (b-a)^k x [aJ]^m-k car I^q = I car I^k=I quel que soit K..donc cela revient juste a multipler par un scalaire la matrice J
A^m = Sum(k=0..m)(k parmi m) [(b-a)I]^k x a^(m-k) x n^(m-k-1) x J
A^m = Sum(k=0..m)(k parmi m) [(b-a)I]^k x a^(m-k) x n^(m-k-1) x J
Et mince, là je ne sais plus comment faire pour continuer ça car je voulais refaire apparaitre un binome mais vu que j'ai pas la meme puissance pour a et n je peux pas,j'avais fais un changement de variable en sortant le m-ieme terme mais en fait ça ne marche pas a cause du (k parmi M)..
Comment jfais ?

fahr451
Membre Transcendant
Messages: 5142
Enregistré le: 05 Déc 2006, 23:50

par fahr451 » 21 Fév 2007, 15:51

attention
J^(k) = n^(k-1) J uniquement pour k>0

il faut sortir de la somme k = 0


le détail

A^(m) =
sigma (k= 1 , m) (k parmi m) a^k n^(k-1) (b-a)^(m-k) + (b-a)^m I

on remet un terme en k = 0
pour faire apparaitre le binôme avec an et b-a


A^m = (1/n) sigma (k = 0 , m) (k parmi m) (an)^k (b-a)^(m-k) J -
(1/n) (b-a)^m J +(b-a)^m I =

(an +b-a)^m J /n +(b-a)^m J/n +(b-a)^m I

Azuriel
Membre Relatif
Messages: 236
Enregistré le: 28 Déc 2006, 16:48

par Azuriel » 21 Fév 2007, 16:03

Ah daccord merci beaucoup, tu as juste oublier de reporter un signe moins de lavan derniere ligne à la derniere pour le terme en k=0 qu'on avait rajouté.

Pour trouver une CNS sur a et b tel que A soit inversible et que notre calcul permettent de calculer cette inverse, mieux vaut repartir de A=(b-a)I +aJ...mais je trouve quelque truc evident mais je ne sais pas comment trouver TOUTES les matrices inversibles..faut il utiliser alors ce qu'on vient de demontrer ?

fahr451
Membre Transcendant
Messages: 5142
Enregistré le: 05 Déc 2006, 23:50

par fahr451 » 21 Fév 2007, 17:30

pour a = b

A = aJ

A^2 = a^2 n J = anA

et si A était inversible en simplifiant par A (multiplier parA^(-1) ) on aurait

A = na I ce qui n'est pas

pour a différent de b

trouve une équation de la forme :

A^2 + alpha A +beta I = 0 avec beta non nul
donc A ( A+alphaI) /beta = I ce qui prouve A inversible d inverse (A+alphaI )/beta

 

Retourner vers ✯✎ Supérieur

Qui est en ligne

Utilisateurs parcourant ce forum : Aucun utilisateur enregistré et 15 invités

Tu pars déja ?



Fais toi aider gratuitement sur Maths-forum !

Créé un compte en 1 minute et pose ta question dans le forum ;-)
Inscription gratuite

Identification

Pas encore inscrit ?

Ou identifiez-vous :

Inscription gratuite