Factorització de rang

De testwiki
Salta a la navegació Salta a la cerca

Donada una matriu A de dimensió m×n amb rang r, una descomposició de rang o factorització de rang de A és un producte A=CF, on C és una matriu m×r i F és una matriu r×n.

Tota matriu de dimensió finita té una factorització de rang: Sigui A una matriu m×n amb rang per columnes r. Per definició, existeixen r columnes linealment independents d'A; de forma equivalent, la dimensió de l'espai de columnes d'A és r. Sigui c1,c2,,cr una base qualsevol de l'espai de columnes d'A, i col·loquem-la com a vectors columna, per formar la matriu C=[c1:c2::cr], de dimensió m×r. Així, qualsevol vector columna d'A és una combinació lineal de les columnes de C. De forma més precisa, si A=[a1:a2::an] és una matriu de dimensió m×n, on aj representa la columna j-sima, llavors

aj=f1jc1+f2jc2++frjcr,

on fij són els coeficients escalars de aj en termes de la base c1,c2,,cr. Això implica que A=CF, on fij és l'element (i,j)-sim de F.

rang(A) = rang (AT)

Una conseqüència immediata de la factorització de rang és que el rang d'A és igual al rang de la seva transposada AT. Com que les columnes d'A són les files d'AT, llavors el rang per columnes d'A és igual al seu rang per files.

Plantilla:Demostració

Factorització de rang per matrius esglaonades per files

A la pràctica, podem construir una factorització de rang de la següent manera: podem calcular B, la forma esglaonada per files d'A. Llavors obtenim C per l'eliminació de les columnes no-pivot d'A, i obtenim F per l'eliminació de les files a 0 de B.

Exemple

Considerem la matriu

A=[1314273915311208][1020011000010000]=B.

B està en forma esglaonada. Llavors obtenim C tot eliminant la tercera columna d'A, l'única que no és una columna pivot, i obtenim F tot eliminant l'última fila de zeros; és a dir:

C=[134279151128],F=[102001100001].

És senzill comprovar que

A=[1314273915311208]=[134279151128][102001100001]=CF.

Demostració

Sigui P una matriu de permutació de dimensió n×n tal que AP=(C,D) en forma particionada per blocs, on les columnes de C són les r columnes pivot d'A. Tota columna de D és una combinació lineal de les columnes de C, de tal manera que existeix una matriu G tal que D=CG, on les columnes de G contenen els coeficients d'aquestes combinacions lineals. Per tant, AP=(C,CG)=C(Ir,G), on Ir és la matriu identitat r×r. Ara veurem que (Ir,G)=FP.

La transformació de AP en la seva forma esglaonada per files és equivalent a multiplicar per l'esquerra per una matriu E que és producte de matrius elementals, de tal forma que EAP=BP=EC(Ir,G), on EC=(Ir0). Ara podem escriure BP=(IrG00), la qual cosa ens permet identificar que (Ir,G)=FP, és a dir, les r files no-nul·les de la forma esglaonada, amb la mateixa permutació de columnes que havíem obtingut per A. Així tenim que AP=CFP, i com que P és invertible, això implica que A=CF, cosa que completa la demostració.

Bibliografia

Plantilla:Millorar referències