3.10

LU-разложение

[10/100%]
Показать
LaTeX
Задача 3.10.1

Пусть A=[1454182631630]\mathbf{A} = \begin{bmatrix} 1 & 4 & 5 \\ 4 & 18 & 26 \\ 3 & 16 & 30 \end{bmatrix}.

?
(a)

Найдите LU-множители матрицы A\mathbf{A}.

(b)

Используя LU-множители, решите Ax→1=b→1\mathbf{A}\overrightarrow {x}_{1} = \overrightarrow {b}_{1}, а также Ax→2=b→2\mathbf{A}\overrightarrow {x}_{2} = \overrightarrow {b}_{2}, где

b→1=(60−6)иb→2=(6612). \overrightarrow {b}_{1} = \begin{pmatrix} 6 \\ 0 \\ -6 \end{pmatrix} \qquad \text{и} \qquad \overrightarrow {b}_{2} = \begin{pmatrix} 6 \\ 6 \\ 12 \end{pmatrix}.
(c)

Используя LU-множители, найдите A−1\mathbf{A}^{-1}.

Задача 3.10.2

Пусть A\mathbf{A} и b→\overrightarrow {b} — матрицы

A=[1241736−12323−3202−26]иb→=(17334). \mathbf{A} = \begin{bmatrix} 1 & 2 & 4 & 17 \\ 3 & 6 & -12 & 3 \\ 2 & 3 & -3 & 2 \\ 0 & 2 & -2 & 6 \end{bmatrix} \qquad \text{и} \qquad \overrightarrow {b} = \begin{pmatrix} 17 \\ 3 \\ 3 \\ 4 \end{pmatrix}.
?
(a)

Объясните, почему матрица A\mathbf{A} не имеет LU-разложения.

(b)

Используя частичный выбор ведущего элемента, найдите матрицу перестановки P\mathbf{P}, а также LU-множители, такие что PA=LU\mathbf{P}\mathbf{A} = \mathbf{L}\mathbf{U}.

(c)

Используя информацию, содержащуюся в P\mathbf{P}, L\mathbf{L} и U\mathbf{U}, решите Ax→=b→\mathbf{A}\overrightarrow {x} = \overrightarrow {b}.

Задача 3.10.3

Найдите все значения ξ\xi, при которых A=[ξ201ξ101ξ]\mathbf{A} = \begin{bmatrix} \xi & 2 & 0 \\ 1 & \xi & 1 \\ 0 & 1 & \xi \end{bmatrix} не имеет LU-разложения.

?
Задача 3.10.4

Если A\mathbf{A} — невырожденная матрица, обладающая LU-разложением, докажите, что ведущий элемент, возникающий после (k+1)(k+1) шагов стандартного метода Гаусса с использованием только операций III типа, задаётся формулой

pk+1=ak+1,k+1−c→TAk−1b→, p_{k+1} = a_{k+1,k+1} - \overrightarrow {c}^{T}\mathbf{A}_{k}^{-1}\overrightarrow {b},

где Ak\mathbf{A}_{k} и Ak+1=(Akb→c→Tak+1,k+1)\mathbf{A}_{k+1} = \begin{pmatrix} \mathbf{A}_{k} & \overrightarrow {b} \\ \overrightarrow {c}^{T} & a_{k+1,k+1} \end{pmatrix} — главные ведущие подматрицы порядков kk и k+1k+1 соответственно. Используя это, выведите, что все ведущие элементы должны быть отличны от нуля, если для A\mathbf{A} существует LU-разложение.

?
Задача 3.10.5

Если A\mathbf{A} — матрица, содержащая только целочисленные элементы, и все её ведущие элементы равны 11, объясните, почему A−1\mathbf{A}^{-1} также должна быть целочисленной матрицей.

?
Примечание.
?

Этот факт можно использовать для построения случайных целочисленных матриц, обладающих целочисленными обратными, случайным образом порождая целочисленные матрицы L\mathbf{L} и U\mathbf{U} с единичными диагоналями и затем составляя произведение A=LU\mathbf{A} = \mathbf{L}\mathbf{U}.

Задача 3.10.6

Рассмотрим трёхдиагональную матрицу T=[β1γ100α1β2γ200α2β3γ300α3β4]\mathbf{T} = \begin{bmatrix} \beta_{1} & \gamma_{1} & 0 & 0 \\ \alpha_{1} & \beta_{2} & \gamma_{2} & 0 \\ 0 & \alpha_{2} & \beta_{3} & \gamma_{3} \\ 0 & 0 & \alpha_{3} & \beta_{4} \end{bmatrix}.

?
(a)

Считая, что T\mathbf{T} обладает LU-разложением, убедитесь, что оно задаётся формулами

L=[1000α1/π11000α2/π21000α3/π31],U=[π1γ1000π2γ2000π3γ3000π4], \mathbf{L} = \begin{bmatrix} 1 & 0 & 0 & 0 \\ \alpha _{1}/\pi _{1} & 1 & 0 & 0 \\ 0 & \alpha _{2}/\pi _{2} & 1 & 0 \\ 0 & 0 & \alpha _{3}/\pi _{3} & 1 \end{bmatrix}, \qquad \mathbf{U} = \begin{bmatrix} \pi _{1} & \gamma _{1} & 0 & 0 \\ 0 & \pi _{2} & \gamma _{2} & 0 \\ 0 & 0 & \pi _{3} & \gamma _{3} \\ 0 & 0 & 0 & \pi _{4} \end{bmatrix},

где числа πi\pi_{i} порождаются рекуррентной формулой

π1=β1иπi+1=βi+1−αiγiπi. \pi _{1} = \beta _{1} \qquad \text{и} \qquad \pi _{i+1} = \beta _{i+1} - \frac{\alpha _{i}\gamma _{i}}{\pi _{i}}.

(Это верно для трёхдиагональных матриц произвольного размера, что делает вычисление LU-множителей таких матриц очень простым.)

(b)

Применяя приведённую выше рекуррентную формулу, найдите LU-разложение матрицы

T=[2−100−12−100−12−100−11]. \mathbf{T} = \begin{bmatrix} 2 & -1 & 0 & 0 \\ -1 & 2 & -1 & 0 \\ 0 & -1 & 2 & -1 \\ 0 & 0 & -1 & 1 \end{bmatrix}.
Задача 3.10.7

Матрица An×n\mathbf{A}_{n \times n} называется ленточной матрицей (band matrix), если aij=0a_{ij} = 0 при ∣i−j∣>w\left|i-j\right| > w для некоторого положительного целого числа ww, называемого шириной ленты (bandwidth). Иными словами, ненулевые элементы A\mathbf{A} ограничены полосой из ww диагональных линий выше и ниже главной диагонали. Например, трёхдиагональные матрицы имеют ширину ленты один, а диагональные матрицы имеют ширину ленты нуль. Если A\mathbf{A} — невырожденная матрица с шириной ленты ww и если A\mathbf{A} имеет LU-разложение A=LU\mathbf{A} = \mathbf{L}\mathbf{U}, то L\mathbf{L} наследует нижнюю ленточную структуру A\mathbf{A}, а U\mathbf{U} наследует верхнюю ленточную структуру в том смысле, что L\mathbf{L} имеет «нижнюю ширину ленты» ww, а U\mathbf{U} имеет «верхнюю ширину ленты» ww. Проиллюстрируйте, почему это так, используя типичную матрицу 5×55 \times 5 с шириной ленты w=2w = 2.

?
Задача 3.10.8
?
(a)

Постройте пример невырожденной симметричной матрицы, не обладающей LU- (или LDU-) разложением.

(b)

Постройте пример невырожденной симметричной матрицы, обладающей LU-разложением, но не являющейся положительно определённой.

Задача 3.10.9

Напомним про LDU-разложение: вынесение диагональных ведущих элементов из верхнетреугольного множителя U\mathbf{U} LU-разложения A=LU\mathbf{A} = \mathbf{L}\mathbf{U} даёт A=LDU′\mathbf{A} = \mathbf{L}\mathbf{D}\mathbf{U}', где D=\diag(u11,u22,…,unn)\mathbf{D} = \diag (u_{11}, u_{22}, \ldots , u_{nn}) — диагональная матрица ведущих элементов, а U′\mathbf{U}' — верхнетреугольная матрица с единицами на диагонали (переобозначая далее U′\mathbf{U}' как U\mathbf{U}, так что и L\mathbf{L}, и U\mathbf{U} в LDU-разложении имеют единичную диагональ).

?
(a)

Найдите LDU-множители для A=[1454182631630]\mathbf{A} = \begin{bmatrix} 1 & 4 & 5 \\ 4 & 18 & 26 \\ 3 & 16 & 30 \end{bmatrix} (это та же матрица, что использовалась в упражнении 3.10.1).

(b)

Докажите, что если матрица обладает LDU-разложением, то LDU-множители определены однозначно.

(c)

Если A\mathbf{A} симметрична и обладает LDU-разложением, объясните, почему оно обязательно имеет вид A=LDLT\mathbf{A} = \mathbf{L}\mathbf{D}\mathbf{L}^{T}.

Задача 3.10.10

Напомним: симметричная матрица A\mathbf{A}, обладающая LU-разложением, в котором все ведущие элементы положительны, называется положительно определённой, и такая матрица однозначно раскладывается как A=RTR\mathbf{A} = \mathbf{R}^{T}\mathbf{R}, где R\mathbf{R} — верхнетреугольная матрица с положительными диагональными элементами; это разложение Холецкого матрицы A\mathbf{A}, а R\mathbf{R} — множитель Холецкого матрицы A\mathbf{A}.

Объясните, почему A=[123281231227]\mathbf{A} = \begin{bmatrix} 1 & 2 & 3 \\ 2 & 8 & 12 \\ 3 & 12 & 27 \end{bmatrix} положительно определена, а затем найдите множитель Холецкого R\mathbf{R}.

?