็ทšๅฝขไปฃๆ•ฐ ยท ็ฌฌ4ๅ›ž ยท ๅ‹•็”ป2 ยท ใ‚คใƒณใ‚ฟใƒฉใ‚ฏใƒ†ใ‚ฃใƒ–ๆ•™ๆ

A = LU ๅˆ†่งฃ

I้‡่ฆใชๅ…ฌๅผ

ๅ…ฌๅผ ๅ็งฐ ่ชฌๆ˜Ž
A=LUA = LU LU ๅˆ†่งฃ ไธ‹ไธ‰่ง’่กŒๅˆ— LL ใจไธŠไธ‰่ง’่กŒๅˆ— UU ใฎ็ฉ
โ„“ij=aijajj\ell_{ij} = \dfrac{a_{ij}}{a_{jj}} ไน—ๆ•ฐ ๆถˆใ—ใŸใ„ๆˆๅˆ†ใ‚’ใƒ”ใƒœใƒƒใƒˆใงๅ‰ฒใฃใŸใ‚‚ใฎ
LL๏ผšๅฏพ่ง’ๆˆๅˆ†ใฏ 1ใ€ใใฎไธ‹ใซไน—ๆ•ฐ LL ใฎๆง‹้€  LL ใฏๆถˆๅŽปๆณ•ใ‹ใ‚‰ใŸใ ใงๆ‰‹ใซๅ…ฅใ‚‹
A=LDUโ€ฒA = LDU' LDU ๅˆ†่งฃ DD ใŒใƒ”ใƒœใƒƒใƒˆใ‚’ใ‚‚ใกใ€LL ใจ Uโ€ฒU' ใฏใฉใกใ‚‰ใ‚‚ๅฏพ่ง’ๆˆๅˆ†ใŒ 1

IIๅฏ่ฆ–ๅŒ–1 โ€” ๆถˆๅŽปๆณ•ใŒ U ใ‚’ไฝœใ‚Šใ€ไน—ๆ•ฐใŒ L ใ‚’ไฝœใ‚‹

ๆถˆๅŽปใฎๅ„ใ‚นใƒ†ใƒƒใƒ—ใฏ UU ใฎๅฏพ่ง’็ทšใ‚ˆใ‚Šไธ‹ใ‚’ 0 ใซใ—ใฆใ„ใใ€ใใ“ใงไฝฟใฃใŸไน—ๆ•ฐใฏใใฎใพใพ LL ใซๅ…ฅใ‚Šใพใ™ใ€‚

ใ‚นใƒ†ใƒƒใƒ—0 ยท ้–‹ๅง‹
U=AU = Aใ€L=IL = I ใ‹ใ‚‰ๅง‹ใ‚ใพใ™ใ€‚ๆถˆๅŽปๆณ•ใง UU ใฎๅฏพ่ง’็ทšใ‚ˆใ‚Šไธ‹ใ‚’ 0 ใซใ—ใ€ไฝฟใฃใŸไน—ๆ•ฐใฏใใ‚Œใžใ‚Œ LL ใซๅ…ฅใ‚Šใพใ™ใ€‚
ใ‚นใƒ†ใƒƒใƒ—1 ยท (2,1) ใ‚’ๆถˆใ™
โ„“21=42=2\ell_{21} = \tfrac{4}{2} = 2 ใชใฎใงใ€R2โ†R2โˆ’2R1R_2 \leftarrow R_2 - 2R_1 ใจใ—ใพใ™ใ€‚22 ใ‚’ L21L_{21} ใซ่จ˜้Œฒใ—ใพใ™ใ€‚
ใ‚นใƒ†ใƒƒใƒ—2 ยท (3,1) ใ‚’ๆถˆใ™
โ„“31=82=4\ell_{31} = \tfrac{8}{2} = 4 ใชใฎใงใ€R3โ†R3โˆ’4R1R_3 \leftarrow R_3 - 4R_1 ใจใ—ใพใ™ใ€‚44 ใ‚’ L31L_{31} ใซ่จ˜้Œฒใ—ใพใ™ใ€‚
ใ‚นใƒ†ใƒƒใƒ—3 ยท (3,2) ใ‚’ๆถˆใ™
โ„“32=31=3\ell_{32} = \tfrac{3}{1} = 3 ใชใฎใงใ€R3โ†R3โˆ’3R2R_3 \leftarrow R_3 - 3R_2 ใจใ—ใพใ™ใ€‚33 ใ‚’ L32L_{32} ใซ่จ˜้Œฒใ—ใพใ™ใ€‚UU ใฏไธŠไธ‰่ง’่กŒๅˆ—ใซใชใ‚Šใพใ—ใŸ โœ“

๐Ÿ’ก LL ใซใฏ่กŒๅŸบๆœฌๅค‰ๅฝขใŒใใฎใพใพ่จ˜้Œฒใ•ใ‚Œใฆใ„ใ‚‹ใฎใงใ€LULU ใฏใใ‚Œใ‚’ๅ†็”Ÿใ—ใฆ AA ใ‚’็ต„ใฟ็ซ‹ใฆ็›ดใ—ใพใ™ โ€” ใ“ใฎๅˆ†่งฃใฏใ€ๆถˆๅŽปๆณ•ใ‚’้€†ๅ‘ใใซใŸใฉใฃใŸใ‚‚ใฎใงใ™ใ€‚

IIIๅฏ่ฆ–ๅŒ–2 โ€” ้€†่กŒๅˆ—ใฎ็ฉใŒใใ‚Œใ„ใชใพใพใงใ‚ใ‚‹็†็”ฑ

EE ใ‚’ใใฎใพใพๆŽ›ใ‘ใ‚‹ใจ (3,1)(3,1) ใซไบคๅทฎ้ …ใŒ็”Ÿใพใ‚Œใพใ™ใŒใ€้€†่กŒๅˆ—ใ‚’้€†ใฎ้ †ๅบใงๆŽ›ใ‘ใ‚‹ใจใ€ใใ‚Œใ„ใช LL ใŒๅพ—ใ‚‰ใ‚Œใพใ™ใ€‚

IVๅฏ่ฆ–ๅŒ–3 โ€” LDU ๅˆ†่งฃ๏ผšใƒ”ใƒœใƒƒใƒˆใ‚’ D ใซใใใ‚Šๅ‡บใ™

UU ใ‹ใ‚‰ใƒ”ใƒœใƒƒใƒˆใ‚’ๅ–ใ‚Šๅ‡บใ—ใฆๅฏพ่ง’่กŒๅˆ— DD ใซใพใจใ‚ใ‚‹ใจใ€LL ใจ Uโ€ฒU' ใฏใฉใกใ‚‰ใ‚‚ๅฏพ่ง’ๆˆๅˆ†ใŒ 1 ใซใชใ‚Šใพใ™ใ€‚

๐Ÿ’ก AA ใŒๅฏพ็งฐ่กŒๅˆ—ใฎใจใใฏ Uโ€ฒ=LTU' = L^{\mathsf{T}} ใจใชใ‚Šใ€ๅˆ†่งฃใฏใใ‚Œใ„ใช A=LDLTA = LDL^{\mathsf{T}} ใฎๅฝขใซใพใจใพใ‚Šใพใ™ใ€‚

V็ขบ่ชใ‚ฏใ‚คใ‚บ

ๅ•1

่กŒๅˆ— A=[3165]A = \begin{bmatrix} 3 & 1 \\ 6 & 5 \end{bmatrix} ใซใคใ„ใฆใ€ๆถˆๅŽปๆณ•ใฎไน—ๆ•ฐ โ„“21\ell_{21} ใฏใ„ใใคใงใ™ใ‹ใ€‚

โœ… ๆญฃ่งฃใงใ™๏ผ ไน—ๆ•ฐใฏ a21/a11=6/3=2a_{21}/a_{11} = 6/3 = 2 ใงใ™ใ€‚

โŒ ๆƒœใ—ใ„๏ผ ๆ€ใ„ๅ‡บใ—ใฆใใ ใ•ใ„๏ผšไน—ๆ•ฐใฏใ€ๆถˆใ—ใŸใ„ๆˆๅˆ†ใ‚’ใใฎไธŠใฎใƒ”ใƒœใƒƒใƒˆใงๅ‰ฒใฃใŸใ‚‚ใฎใงใ™๏ผšโ„“21=a21/a11\ell_{21} = a_{21} / a_{11}ใ€‚

่งฃ็ญ”ใ‚’่ฆ‹ใ‚‹

่งฃ็ญ”๏ผš

ไน—ๆ•ฐใฏใ€ๆถˆใ—ใŸใ„ๆˆๅˆ†ใ‚’ใƒ”ใƒœใƒƒใƒˆใงๅ‰ฒใฃใŸใ‚‚ใฎใงใ™๏ผš

โ„“21=a21a11=63=2\ell_{21} = \frac{a_{21}}{a_{11}} = \frac{6}{3} = 2

ไน—ๆ•ฐใฏ 2 ใงใ™ใ€‚2่กŒ็›ฎใ‹ใ‚‰1่กŒ็›ฎใฎ2ๅ€ใ‚’ๅผ•ใ„ใฆใ€ใƒ”ใƒœใƒƒใƒˆใฎไธ‹ใซ 0 ใ‚’ไฝœใ‚Šใพใ™ใ€‚

ๅ•2

LU ๅˆ†่งฃ A=LUA = LU ใซใŠใ„ใฆใ€LL ใฎๅฏพ่ง’ๆˆๅˆ†ใฏใ„ใคใ‚‚ๆฌกใฎใฉใ‚Œใงใ™ใ‹ใ€‚

โœ… ๆญฃ่งฃใงใ™๏ผ LL ใฎๅฏพ่ง’ๆˆๅˆ†ใฏใ„ใคใ‚‚ 1 ใงใ™ โ€” ไน—ๆ•ฐใฏๅฏพ่ง’็ทšใฎไธ‹ใซไธฆใณใ€ใƒ”ใƒœใƒƒใƒˆใฏ UU ใซใ‚ใ‚Šใพใ™ใ€‚

โŒ ๆƒœใ—ใ„๏ผ LL ใฎๆง‹้€ ใ‚’่€ƒใˆใฆใฟใพใ—ใ‚‡ใ†๏ผšใ“ใ‚ŒใฏๆถˆๅŽป่กŒๅˆ—ใฎ้€†่กŒๅˆ—ใฎ็ฉใงใ€ใใฎใฉใ‚Œใ‚‚ใŒๅฏพ่ง’ๆˆๅˆ† 1 ใ‚’ใ‚‚ใกใพใ™ใ€‚

่งฃ็ญ”ใ‚’่ฆ‹ใ‚‹

่งฃ็ญ”๏ผš

LL ใฏไธ‹ไธ‰่ง’่กŒๅˆ—ใงใ€ๅฏพ่ง’ๆˆๅˆ†ใฏใ™ในใฆ 1ใ€ใใฎไธ‹ใซๆถˆๅŽปๆณ•ใฎไน—ๆ•ฐใŒไธฆใณใพใ™๏ผš

L=[100โ„“2110โ„“31โ„“321]L = \begin{bmatrix} 1 & 0 & 0 \\ \ell_{21} & 1 & 0 \\ \ell_{31} & \ell_{32} & 1 \end{bmatrix}

  • ใƒ”ใƒœใƒƒใƒˆใŒ็พใ‚Œใ‚‹ใฎใฏ LL ใงใฏใชใใ€UU ใฎๅฏพ่ง’็ทšไธŠ๏ผˆLDU ๅˆ†่งฃใชใ‚‰ DD ใฎไธญ๏ผ‰ใงใ™ใ€‚
  • ไน—ๆ•ฐใŒๅ…ฅใ‚‹ใฎใฏ LL ใฎๅฏพ่ง’็ทšใฎไธŠใงใฏใชใใ€ใใฎไธ‹ใงใ™ใ€‚

ๅ•3

ๆญฃใ—ใ„ใ‹่ชคใ‚Šใ‹๏ผšๆถˆๅŽป่กŒๅˆ—ใฎ็ฉ E32โ‹…E21E_{32} \cdot E_{21} ใ‚’่จˆ็ฎ—ใ™ใ‚‹ใจใ€(3,1)(3,1) ใฎไฝ็ฝฎใซไฝ™่จˆใชไบคๅทฎ้ …ใŒ็พใ‚Œใ‚‹ใ“ใจใŒใ‚ใ‚‹ใ€‚ใ—ใ‹ใ—ใ€ใใ‚Œใ‚‰ใฎ้€†่กŒๅˆ—ใ‚’้€†ใฎ้ †ๅบใงๆŽ›ใ‘ใŸ E21โˆ’1โ‹…E32โˆ’1E_{21}^{-1} \cdot E_{32}^{-1} ใงใฏใ“ใฎไบคๅทฎ้ …ใŒ็”Ÿใ˜ใšใ€ใใ‚Œใ„ใช่กŒๅˆ— LL ใŒๅพ—ใ‚‰ใ‚Œใ‚‹ใ€‚

โœ… ๆญฃ่งฃใงใ™๏ผ ้€†่กŒๅˆ—ใฎ็ฉใซใฏไบคๅทฎ้ …ใŒ็”Ÿใ˜ใพใ›ใ‚“ โ€” ใ ใ‹ใ‚‰ใ“ใ LL ใ‚’ๆถˆๅŽป่กŒๅˆ—ใฎ้€†่กŒๅˆ—ใฎ็ฉใจใ—ใฆๅฎš็พฉใ—ใ€ไน—ๆ•ฐใŒใใฎใพใพ LL ใซๅŽใพใ‚‹ใฎใงใ™ใ€‚

โŒ ๆƒœใ—ใ„๏ผ ไธŠใฎๅฏ่ฆ–ๅŒ–2ใงไน—ๆ•ฐใ‚’ 0 ไปฅๅค–ใซใ—ใฆใ€2ใคใฎ็ฉใฎ (3,1) ๆˆๅˆ†ใ‚’ๆฏ”ในใฆใฟใพใ—ใ‚‡ใ†ใ€‚

่งฃ็ญ”ใ‚’่ฆ‹ใ‚‹

่งฃ็ญ”๏ผšๆญฃใ—ใ„

E32โ‹…E21E_{32} \cdot E_{21} ใ‚’่จˆ็ฎ—ใ™ใ‚‹ใจใ€E32E_{32} ใฎๆˆๅˆ† (โˆ’m32)(-m_{32}) ใจ E21E_{21} ใฎๆˆๅˆ† (โˆ’m21)(-m_{21}) ใŒๆŽ›ใ‘ๅˆใ‚ใ•ใ‚Œใ€(3,1)(3,1) ใฎไฝ็ฝฎใซไบคๅทฎ้ … m32โ‹…m21m_{32} \cdot m_{21} ใŒ็”Ÿใพใ‚Œใพใ™๏ผš

E32โ‹…E21=[100โˆ’m2110m32m21โˆ’m321]E_{32} \cdot E_{21} = \begin{bmatrix} 1 & 0 & 0 \\ -m_{21} & 1 & 0 \\ m_{32}m_{21} & -m_{32} & 1 \end{bmatrix}

ใ—ใ‹ใ—ใ€้€†่กŒๅˆ—ใ‚’ๆญฃใ—ใ„้ †ๅบใงๆŽ›ใ‘ใ‚Œใฐใ€ใ“ใ‚ŒใฏใพใฃใŸใ่ตทใ“ใ‚Šใพใ›ใ‚“๏ผš

L=E21โˆ’1โ‹…E32โˆ’1=[100m21100m321]L = E_{21}^{-1} \cdot E_{32}^{-1} = \begin{bmatrix} 1 & 0 & 0 \\ m_{21} & 1 & 0 \\ 0 & m_{32} & 1 \end{bmatrix}

ไบคๅทฎ้ …ใฏใ‚ใ‚Šใพใ›ใ‚“ โ€” ใฉใฎไน—ๆ•ฐใ‚‚ๆœฌๆฅใฎไฝ็ฝฎใซใใ‚Œใ„ใซๅŽใพใ‚Šใพใ™ใ€‚ใ“ใ‚Œใฏใ€ๅ„้€†่กŒๅˆ—ใŒๆ›ธใๆ›ใˆใ‚‹่กŒใซใ€ใ‚ใจใฎ้€†่กŒๅˆ—ใŒๆ‰‹ใ‚’่งฆใ‚Œใชใ„ใ‹ใ‚‰ใงใ™ใ€‚

ๅ•4

L=[1031]L = \begin{bmatrix} 1 & 0 \\ 3 & 1 \end{bmatrix}ใ€U=[4205]U = \begin{bmatrix} 4 & 2 \\ 0 & 5 \end{bmatrix} ใฎใจใใ€A=LUA = LU ใฏใฉใ‚Œใงใ™ใ‹ใ€‚

โœ… ๆญฃ่งฃใงใ™๏ผ AA ใฎ2่กŒ็›ฎ๏ผš[3โ‹…4+1โ‹…0,โ€…โ€Š3โ‹…2+1โ‹…5]=[12,โ€…โ€Š11][3 \cdot 4 + 1 \cdot 0, \; 3 \cdot 2 + 1 \cdot 5] = [12, \; 11]ใ€‚

โŒ ๆƒœใ—ใ„๏ผ ๆ€ใ„ๅ‡บใ—ใฆใใ ใ•ใ„๏ผšAA ใฎๅ„ๆˆๅˆ†ใฏใ€LL ใฎ่กŒใจ UU ใฎๅˆ—ใฎๅ†…็ฉใงใ™ใ€‚(2,2)(2,2) ใฎไฝ็ฝฎใ‚’็ขบใ‹ใ‚ใพใ—ใ‚‡ใ†๏ผš3โ‹…2+1โ‹…5=?3 \cdot 2 + 1 \cdot 5 = ?

่งฃ็ญ”ใ‚’่ฆ‹ใ‚‹

่งฃ็ญ”๏ผš

Lร—UL \times U ใ‚’1่กŒใšใค่จˆ็ฎ—ใ—ใพใ™๏ผš

1่กŒ็›ฎ๏ผš[1,0]โ‹…[4205]=[1โ‹…4+0โ‹…0,โ€…โ€Š1โ‹…2+0โ‹…5]=[4,โ€…โ€Š2][1, 0] \cdot \begin{bmatrix} 4 & 2 \\ 0 & 5 \end{bmatrix} = [1 \cdot 4 + 0 \cdot 0, \; 1 \cdot 2 + 0 \cdot 5] = [4, \; 2]

2่กŒ็›ฎ๏ผš[3,1]โ‹…[4205]=[3โ‹…4+1โ‹…0,โ€…โ€Š3โ‹…2+1โ‹…5]=[12,โ€…โ€Š11][3, 1] \cdot \begin{bmatrix} 4 & 2 \\ 0 & 5 \end{bmatrix} = [3 \cdot 4 + 1 \cdot 0, \; 3 \cdot 2 + 1 \cdot 5] = [12, \; 11]

A=LU=[421211]A = LU = \begin{bmatrix} 4 & 2 \\ 12 & 11 \end{bmatrix}

ใ‚ˆใใ‚ใ‚‹้–“้•ใ„๏ผš

  • [42126]\begin{bmatrix} 4 & 2 \\ 12 & 6 \end{bmatrix}๏ผš(2,2)(2,2) ใฎไฝ็ฝฎใง 1โ‹…51 \cdot 5 ใฎ้ …ใ‚’ๅฟ˜ใ‚Œใ€3โ‹…2=63 \cdot 2 = 6 ใ—ใ‹่จˆ็ฎ—ใ—ใฆใ„ใพใ›ใ‚“ใ€‚
  • [4236]\begin{bmatrix} 4 & 2 \\ 3 & 6 \end{bmatrix}๏ผšๅ†…็ฉใ‚’่จˆ็ฎ—ใ›ใšใซใ€่กŒใฉใ†ใ—ใ‚’ๆˆๅˆ†ใ”ใจใซใŸใ—ๅˆใ‚ใ›ใฆใ„ใพใ™ใ€‚
  • [4205]\begin{bmatrix} 4 & 2 \\ 0 & 5 \end{bmatrix}๏ผšใ“ใ‚Œใฏ UU ใใฎใ‚‚ใฎใงใ™ โ€” LL ใ‚’ๆŽ›ใ‘ใ‚‹ใฎใ‚’ๅฟ˜ใ‚Œใฆใ„ใพใ™ใ€‚

ๆญฃ่งฃๆ•ฐ๏ผš0 / 4