Bazo

La kerno de matrico estas ne influata per rudimentaj liniaj operacioj (kio estas aladicio de unu linio al la alia kaj multipliko de iu linio je ne-nula nombro). Tiel eblas per gaŭsa eliminado kalkuli bazon de la kerno:

Estu m×n matrico A. Oni uzu la rudimentajn liniajn operacioj kaj faru gaŭsan eliminadon de A. Tiam al ĉiu i-a kolumno de eliminigita matrico, en kiu estas nur unu valoro 1 kaj ĉiuj ceteraj estas 0, respektivas la dependa variablo x_i; al ĉiu cetera kolumno respektivas la libera variablo x_i.

Por ĉiu dependa variablo, el tiu linio de la matrico kie estas la valoro 1, estas formata ekvacio per preno de la eroj de la linio kiel koeficientoj ĉe la liberaj variabloj kun la kontraŭa sugno. La bazaj vektoroj povas esti ricevitaj per preno de unu libera variablo egala al 1 kaj la aliaj egalaj al 0, tiel kvanto de la bazaj vektoroj egalas al kvanto de la liberaj variabloj.

Ekzemple, supozu ke post gaŭsa eliminado rezultas matrcio

${\displaystyle {\begin{bmatrix}1&0&-4&0&2&-6\\0&1&5&0&-9&3\\0&0&0&1&-7&8\\0&0&0&0&0&0\end{bmatrix}}}$

Tiam x₃, x₅, kaj x₆ estas liberaj, kaj por la dependaj variabloj estas ekvacioj

${\displaystyle x_{1}=4x_{3}-2x_{5}+6x_{6}}$

${\displaystyle x_{2}=-5x_{3}+9x_{5}-3x_{6}}$

${\displaystyle x_{4}=7x_{5}-8x_{6}}$

kaj tri vektoroj ${\displaystyle {\begin{bmatrix}4\\-5\\1\\0\\0\\0\end{bmatrix}}}$, ${\displaystyle {\begin{bmatrix}-2\\9\\0\\7\\1\\0\end{bmatrix}}}$, ${\displaystyle {\begin{bmatrix}6\\-3\\0\\-8\\0\\1\end{bmatrix}}}$ estas bazo de la kerno de A.

Rilato al la linia spaco

Estu m×n matrico A. La produto de A kaj la n-dimensia vektoro x povas esti skribita per la skalara produto de vektoroj:

${\displaystyle A\mathbf {x} ={\begin{bmatrix}\mathbf {a} _{1}\cdot \mathbf {x} \\\mathbf {a} _{2}\cdot \mathbf {x} \\\vdots \\\mathbf {a} _{m}\cdot \mathbf {x} \end{bmatrix}}}$

Ĉi tie a₁, ..., a_m estas la linioj de la matrico A. Tiel x estas en la kerno de A se kaj nur se x estas perpendikulara al ĉiu el la liniaj vektoroj de A (ĉar tio ke la skalara produto de du vektoroj estas egala al nulo estas difino de tio ke ili estas perpendikularaj).

La linia spaco de matrico A estas la lineara generaĵo de la liniaj vektoroj de A. Tiel la kerno de A estas la perpendikulara komplemento al la linia spaco. Tio estas ke vektoro x kuŝas en la kerno de A se kaj nur se ĝi estas perpendikulara al ĉiu vektoro en la linia spaco de A.

La dimensio de la linia spaco de A estas la rango de A.

La rango kaj la dimensio de kerno de A estas ligitaj kiel

rank(A) + dim(ker(A)) = n

La ekvacio pli supre estas sciata kiel la rango-kerna teoremo.

Per SVD

SVD de matrico A komputas unitajn matricojn U kaj V kaj ortangulan diagonalan matricon S de la sama amplekso kiel A kun nenegativaj diagonalaj elementoj, tiaj ke

A=USV

Signifu la kolumnojn de V per

v₁, ..., v_n

kaj la diagonalajn elementojn de S per

s₁, ..., s_min(m,n)

kaj estu

s_{min(m, n)+1} = ... = s_{max(m, n)} = 0

La nombroj s_i estas la singularaj valoroj de A. Tiam la kolumnoj v₁ de V tiaj ke la respektiva s_i=0 formas ortonormalam bazon de la kerno de A.

En cifereca kalkulado, ĉiu el la singularaj valoroj s_i estas konsiderata kiel nulo se ĝi estas pli malgranda ol iu certa malgranda valoro de tolero. Ekzemple en MATLAB la valoro de tolero estas prenita al esti

max (m, n) max {s_i}ε

kie ε estas la maŝina epsilono de la komputilo, kio estas, la plej malgranda nombro tia ke en la flosanta punkto, 1+ε>1. Por la IEEE 64-bita flosanta punkta aranĝo ε≈2,2·10⁻¹⁶.

Per QR-faktorigo

Eblas komputi la kernon per la QR-faktorigo jene.

Estu A m×n matrico kun m≤n. Per la QR-faktorigo de A^* kun elekto de kondukaj kolumnoj oni povas trovi matricojn P, Q, R tiajn ke

A^*P=QR

kie P estas permuta matrico, Q estas n×n unita matrico kaj R estas n×m supra triangula matrico.

Eblas rearanĝi la formulon, konsiderante tion ke P⁻¹=P^* kaj Q⁻¹=Q^*:

Q⁻¹A^*P = R

Q^*A^*P = R

(AQ)^*P = R

(AQ)^* = RP⁻¹

(AQ)^* = RP^*

AQ = (RP^*)^*

AQ = PR^*

Se rango de A estas malpli granda ol n, la lastaj linioj de R estas tute nulaj. Tiam la lastaj kolumnoj de R^* estas tute nulaj, kaj do iuj kolumnoj de PR^* estas tute nulaj.

Se ${\displaystyle Q={\begin{bmatrix}\mathbf {q} _{1}&\ldots &\mathbf {q} _{n}\end{bmatrix}}}$ kaj ${\displaystyle PR^{*}={\begin{bmatrix}\mathbf {b} _{1}&\ldots &\mathbf {b} _{n}\end{bmatrix}}}$ do por ĉiu i=1...n, Aq_i=b_i

Tiel tiuj kolumnoj q_i de Q kiuj respektivas al nulaj kolumnoj b_i de PR^* naskas la kernon de A.

Same kiel en okazo de SVD, en cifereca kalkulado, ĉiu el la b_i estas konsiderata kiel nula se iu ĝia normo estas pli malgranda ol iu elektita malgranda valoro de tolero.