Vectorizare (matematică)

În matematică, în special în algebră liniară, vectorizarea unei matrice este o transformare liniară care transformă matricea într-un vector. Mai exact, vectorizarea unei matrice A Format:Nowrap, denumită vec(A), este un vector coloană Format:Nowrap obținut prin aranjare coloanelor matricei A una peste alta:

$\mathrm{vec}(A)=[a_{1,1},\dots ,a_{m,1},a_{1,2},\dots ,a_{m,2},\dots ,a_{1,n},\dots ,a_{m,n}{]}^{\mathrm{T}}$

Aici, $a_{i,j}$ reprezintă elementul din linia i și coloana j din A, iar cu indicele superior ${}^{\mathrm{T}}$ este notată transpunerea. Vectorizarea exprimă, prin coordonate, izomorfismul ${𝐑}^{m\times n}:={𝐑}^m\otimes {𝐑}^n\cong {𝐑}^{mn}$ între acestea (adică dintre matrici și vectori) ca spații vectoriale.

De exemplu, pentru matricea 2×2 $A=\left[\begin{array}{cc}a& b\\ c& d\end{array}\right]$, vectorizarea este $\mathrm{vec}(A)=\left[\begin{array}{c}a\\ c\\ b\\ d\end{array}\right]$.

Legătura dintre vectorizarea lui A și vectorizarea transpusei sale este dată de matricea de comutare.

Compatibilitatea cu produsul Kronecker

Vectorizarea este folosită frecvent împreună cu produsul Kronecker pentru a exprima înmulțirea matricilor ca o transformare liniară pe matrici. În special,

$\mathrm{vec}(ABC)=(C^{\mathrm{T}}\otimes A)\mathrm{vec}(B)$

pentru matricele A, B și C cu dimensiunile k×l, l×m și m×n. De exemplu, dacă ${\mathrm{ad}}_A(X)=AX-XA$ (endomorfismul adjunct în Format:Ill-wd Format:Nowrap al tuturor matricilor n×n complexe), atunci $\mathrm{vec}({\mathrm{ad}}_A(X))=(I_n\otimes A-A^{\mathrm{T}}\otimes I_n)\text{vec}(X)$, unde $I_n$ este matricea unitate n×n.

Există alte două formulări utile:

$\begin{array}{cc}\mathrm{vec}(ABC)& =(I_n\otimes AB)\mathrm{vec}(C)=(C^{\mathrm{T}}{B}^{\mathrm{T}}\otimes I_k)\mathrm{vec}(A)\\ \mathrm{vec}(AB)& =(I_m\otimes A)\mathrm{vec}(B)=(B^{\mathrm{T}}\otimes I_k)\mathrm{vec}(A)\end{array}$

Mai general, s-a demonstrat că vectorizarea este autoadjunctă în structura închisă monoidală a oricărei categorii de matrice.

Compatibilitatea cu produsul Hadamard

Vectorizarea este un homomorfism algebric din spațiul matricilor Format:Nowrap cu produsul Hadamard la Cⁿ² cu produsul său Hadamard:

$\mathrm{vec}(A\circ B)=\mathrm{vec}(A)\circ \mathrm{vec}(B).$

Compatibilitatea cu produsul interior

Vectorizarea este o transformare unitară din spațiul matricelor n×n cu produsul interior Frobenius sau Format:Ill-wd la Cⁿ²:

$$\mathrm{tr}(A^{†}B)=\mathrm{vec}(A{)}^{†}\mathrm{vec}(B),$$

Unde cu indicele superior ^† este notată adjuncta.

Vectorizarea ca sumă liniară

Operația de vectorizare a unei matrice poate fi scrisă ca o sumă liniară. Fie X matricea Format:Nowrap care să fie vectorizată și fie e_i vectorul celei de-a i-a bază canonică a spațiului n-dimensional, adică ${𝐞}_i={[0,\dots ,0,1,0,\dots ,0]}^{\mathrm{T}}$. Fie B_i o matrice de blocuri Format:Nowrap definită astfel:

${𝐁}_i=\left[\begin{array}{c}\mathrm{𝟎}\\ ⋮\\ \mathrm{𝟎}\\ {𝐈}_m\\ \mathrm{𝟎}\\ ⋮\\ \mathrm{𝟎}\end{array}\right]={𝐞}_i\otimes {𝐈}_m$

B_i este formată din n matrici de blocuri de dimensiuni Format:Nowrap, plasate una deasupra alteia, și toate aceste matrici sunt nule cu excepția celei de a i-a, care este matricea unitate Format:Nowrap, I_m.

Apoi versiunea vectorizată a lui X poate fi exprimată după cum urmează:

$$\mathrm{vec}(𝐗)=\sum _{i=1}^n{𝐁}_i𝐗{𝐞}_i$$

Înmulțirea lui X cu e_i extrage coloana a i-a, în timp ce înmulțirea cu B_i o pune în poziția dorită în vectorul final.

Alternativ, suma liniară poate fi exprimată folosind produsul Kronecker:

$\mathrm{vec}(𝐗)=\sum _{i=1}^n{𝐞}_i\otimes 𝐗{𝐞}_i$

Semivectorizare

Pentru o matrice simetrică A, vectorul vec(A) conține mai multe informații decât este strict necesar, deoarece matricea este complet determinată de simetrie împreună cu Format:Ill-wd, adică Format:Nowrap de pe și sub diagonala principală. Pentru astfel de matrici, semivectorizarea este uneori mai utilă decât vectorizarea. Semivectorizarea, vech(A), a unei matrice simetrice Format:Nowrap A este vectorul coloană Format:Nowrap obținut prin vectorizarea numai a părții triunghiulare inferioare a lui A:

$\mathrm{vech}(A)=[A_{1,1},\dots ,A_{n,1},A_{2,2},\dots ,A_{n,2},\dots ,A_{n-1,n-1},A_{n,n-1},A_{n,n}{]}^{\mathrm{T}}.$

De exemplu, pentru matricea $A=\left[\begin{array}{cc}a& b\\ b& d\end{array}\right]$ 2×2, semivectorizarea este $\mathrm{vech}(A)=\left[\begin{array}{c}a\\ b\\ d\end{array}\right]$.

Există matrici unice care transformă semivectorizarea unei matrice în vectorizarea acesteia și invers numite, respectiv, matrice de duplicare și matrice de eliminare.

În limbajele de programare

Limbajele de programare care au implementate operații matriciale pot avea mijloace ușoare de vectorizare. În Matlab/GNU Octave o matrice A poate fi vectorizată prin A(:). GNU Octave permite, de asemenea, vectorizarea și semivectorizarea prin vec(A) și, respectiv, vech(A). Julia are și el funcția vec(A).

În Python tablourile NumPy implementează metoda platten, iar în R efectul dorit poate fi obținut cu funcțiile c() sau as.vector(). În R funcția vec() a pachetului „ks” permite vectorizarea, iar funcția vech() implementată în ambele pachete „ks” și „sn” permite semivectorizarea.

Bibliografie

• Format:En icon Jan R. Magnus and Heinz Neudecker (1999), Matrix Differential Calculus with Applications in Statistics and Econometrics, 2nd Ed., Wiley. Format:Isbn.

Format:Portal