LU malkomponaĵo estas metodo de solvado de sistemo de linearaj ekvacioj, kiu havas formon:
Laŭ matrica maniero de skribado
, kaj
- matrico de koeficientoj ,
- vektoro de variabloj,
- vektoro de datumoj.
Laŭ metodo LU matrico de koeficientoj estas malkomponata en multiplikado de du triangulaj matricoj: suba (angle: Lower) kaj supera (angle: Upper).
![{\displaystyle \mathbf {A} =\mathbf {L} \cdot \mathbf {U} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/85b936e264fe11d9ad94edcb090527aab7508d8f)
Sistemo de ekvacioj tiam havas formon:
![{\displaystyle \mathbf {L} \cdot \mathbf {U} \cdot \mathbf {x} =\mathbf {y} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/9879d9f4603c77a25d9a72308f676679e5b78adb)
kaj ĝia solvo povas esti esprimita per du sistemoj de ekvacioj kun triangulaj ekvacioj, kiuj solvas tre facile.
![{\displaystyle \mathbf {L} \cdot \mathbf {z} =\mathbf {y} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/e078627b636ec5eb821b997d07ae0ab0ce11604d)
- La metodo ebligas rapidan kalkuladon de matrico
.
- Kvanto de multiplikoj kiuj estas bezonata por solvi (kalkuli valorojn de vektoro
) estas
kaj adicioj
.
- por komputila kalkulado metodo ŝparas memoron, ĉar ĉiujn valorojn oni povas havi en unu matrico kaj unu vektoro (samaj kiuj enhavas komencajn valorojn).
- Ĝi bezonas malmultajn operacioj ol aliaj metodoj (krom specialaj metodoj)
- Kalkulado de determinanto de matrico A povas kalkuli uzante teoremon de Cauchy:
- kaj uzante de fakto, ke determino de triangula matrico estas multipliko de diagonalaj elementoj de matrico. Alinome:
Ĉefa problemo de ĉi tiu metodo estas malkomponaĵo. Por ke la malkomponaĵo estus unusignifa, decidas ke iu el du matricoj havas diagonalaj elementoj egalaj al unu.
Ekzistas du ĉefaj metodoj por fari tiun:
- Metodo elimina de Gauss
- Metodo de Doolittle (priskribo sube)
Laŭ ĉi tiu metodo egaleco
estas sistemo de
ekvacioj kun
variabloj. La variabloj estas elementoj
por
(elementoj sube diagonalo) kaj
por
(elementoj de diagonalo kaj supere). Kun lemo ke elementoj de diagonalo de matrico L ekvacias 1.
Kalkulado sekvaj elementoj de matricoj
kaj
faras alterne. te. post kalkulado de verso de matrico U kalkulas kolumnon de matrico L kaj denove sekvan verson U.
Ĝeneralaj formuloj por apartaj elementoj de matricoj estas:
por ĉiu
:
dla ![{\displaystyle j=i,i+1\ldots ,n}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f0a2fa03e501ac132388fa2f55dbd2e1ab7e10ac)
dla ![{\displaystyle j=i+1,i+2\ldots ,n}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bda9d4f250a47fbde248a6a029f4a603646066dd)
Laŭ lasta ekvacio metodo ne funkcios, se
.
Kvanto de bezonataj operacioj:
- multiplikaj:
,
- adiciaj:
.
![{\displaystyle {\begin{bmatrix}5&3&2\\1&2&0\\3&0&4\end{bmatrix}}={\begin{bmatrix}1&0&0\\l_{21}&1&0\\l_{31}&l_{32}&1\end{bmatrix}}{\begin{bmatrix}u_{11}&u_{12}&u_{13}\\0&u_{22}&u_{23}\\0&0&u_{33}\\\end{bmatrix}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c7da5e47d32cade5b50ec5eec80e0d80c459bd71)
Unua verso de matrico U:
![{\displaystyle 5=1\cdot u_{11}+0\cdot 0+0\cdot 0\rightarrow u_{11}=5}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8e3b02b80314f84cc3a2b42949ce8c8422e0fab5)
![{\displaystyle 3=1\cdot u_{12}+0\cdot u_{22}+0\cdot 0\rightarrow u_{12}=3}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5e039b2ceb4a86cd1f89cfff50d2206fdfeb53f7)
![{\displaystyle 2=1\cdot u_{13}+0\cdot u_{23}+0\cdot u_{33}\rightarrow u_{13}=2}](https://wikimedia.org/api/rest_v1/media/math/render/svg/85e41e40d851ab68ecfefc841221cab25a2903c4)
![{\displaystyle {\begin{bmatrix}5&3&2\\1&2&0\\3&0&4\end{bmatrix}}={\begin{bmatrix}1&0&0\\l_{21}&1&0\\l_{31}&l_{32}&1\end{bmatrix}}{\begin{bmatrix}5&3&2\\0&u_{22}&u_{23}\\0&0&u_{33}\\\end{bmatrix}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d5ec2c88872e93d94b49c88c75eff85a31e53f33)
Unua kolumno de matrico L:
![{\displaystyle 1=l_{21}\cdot 5+1\cdot 0+0\cdot 0\rightarrow l_{21}={\frac {1}{5}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7ef3657d5a56abfb8c1cc41c15ce1c28463e9932)
![{\displaystyle 3=l_{31}\cdot 5+l_{32}\cdot 0+1\cdot 0\rightarrow l_{31}={\frac {3}{5}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bae05491420de4e197fa6858664815b85af6cc73)
![{\displaystyle {\begin{bmatrix}5&3&2\\1&2&0\\3&0&4\end{bmatrix}}={\begin{bmatrix}1&0&0\\{\frac {1}{5}}&1&0\\{\frac {3}{5}}&l_{32}&1\end{bmatrix}}{\begin{bmatrix}5&3&2\\0&u_{22}&u_{23}\\0&0&u_{33}\\\end{bmatrix}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fd3be2631a8b730021306e1752681f3b25ae8139)
Dua verso de matrico U:
![{\displaystyle 2={\frac {1}{5}}\cdot 3+1\cdot u_{22}+0\cdot 0\rightarrow u_{22}={\frac {7}{5}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/de8ad5a86feed5f084238388d5ec0f095761ead0)
![{\displaystyle 0={\frac {1}{5}}\cdot 2+1\cdot u_{23}+0\cdot u_{33}\rightarrow u_{23}=-{\frac {2}{5}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5eccc3aae4e19189b0922927cab25e2aaa563019)
![{\displaystyle {\begin{bmatrix}5&3&2\\1&2&0\\3&0&4\end{bmatrix}}={\begin{bmatrix}1&0&0\\{\frac {1}{5}}&1&0\\{\frac {3}{5}}&l_{32}&1\end{bmatrix}}{\begin{bmatrix}5&3&2\\0&{\frac {7}{5}}&-{\frac {2}{5}}\\0&0&u_{33}\\\end{bmatrix}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7dc6f2bc91cb61797058e205fe0956f5f1b4bb95)
Dua kolumno de matrico L:
![{\displaystyle 0={\frac {3}{5}}\cdot 3+l_{32}\cdot {\frac {7}{5}}+1\cdot 0\rightarrow l_{32}=-{\frac {9}{7}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5d24f41a1eb7a65f20f00c9c0a2822378bd1ae86)
![{\displaystyle {\begin{bmatrix}5&3&2\\1&2&0\\3&0&4\end{bmatrix}}={\begin{bmatrix}1&0&0\\{\frac {1}{5}}&1&0\\{\frac {3}{5}}&-{\frac {9}{7}}&1\end{bmatrix}}{\begin{bmatrix}5&3&2\\0&{\frac {7}{5}}&-{\frac {2}{5}}\\0&0&u_{33}\\\end{bmatrix}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7139832d0f8378437564215df78bb23506189f12)
Tria verso de matrico U:
![{\displaystyle 4={\frac {3}{5}}\cdot 2+{\frac {9}{7}}\cdot {\frac {2}{5}}+1\cdot u_{33}\rightarrow u_{33}={\frac {16}{7}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/55f7171eaae0ef0d202a6515772a3fa3b0645d31)
![{\displaystyle {\begin{bmatrix}5&3&2\\1&2&0\\3&0&4\end{bmatrix}}={\begin{bmatrix}1&0&0\\{\frac {1}{5}}&1&0\\{\frac {3}{5}}&-{\frac {9}{7}}&1\end{bmatrix}}{\begin{bmatrix}5&3&2\\0&{\frac {7}{5}}&-{\frac {2}{5}}\\0&0&{\frac {16}{7}}\\\end{bmatrix}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b8153afc975117cc1c0bf2420e0661ea6cca5b78)