Глава 6

Определители

[48/100%]
Показать
LaTeX
§
Задача 6.1.1

Используя определение, вычислите det⁡(A)\operatorname {det}\left(\mathbf{A}\right) для каждой из следующих матриц.

?
(a)

A=(3−21−540216)\mathbf{A} = \begin{pmatrix} 3 & -2 & 1 \\ -5 & 4 & 0 \\ 2 & 1 & 6 \end{pmatrix}.

(b)

A=(211621−221)\mathbf{A} = \begin{pmatrix} 2 & 1 & 1 \\ 6 & 2 & 1 \\ -2 & 2 & 1 \end{pmatrix}.

(c)

A=(00α0β0γ00)\mathbf{A} = \begin{pmatrix} 0 & 0 & \alpha \\ 0 & \beta & 0 \\ \gamma & 0 & 0 \end{pmatrix}.

(d)

A=(a11a12a13a21a22a23a31a32a33)\mathbf{A} = \begin{pmatrix} a_{11} & a_{12} & a_{13} \\ a_{21} & a_{22} & a_{23} \\ a_{31} & a_{32} & a_{33} \end{pmatrix}.

Задача 6.1.2

Чему равен объём параллелепипеда, натянутого на три вектора x→1=(3,0,−4,0)T\overrightarrow {x}_{1} = (3,0,-4,0)^{T}, x→2=(0,2,0,−2)T\overrightarrow {x}_{2} = (0,2,0,-2)^{T} и x→3=(0,1,0,1)T\overrightarrow {x}_{3} = (0,1,0,1)^{T}?

?
Задача 6.1.3

Используя метод Гаусса для приведения A\mathbf{A} к верхнетреугольному виду, вычислите det⁡(A)\operatorname {det}\left(\mathbf{A}\right) для каждой из следующих матриц.

?
(a)

A=(123241144)\mathbf{A} = \begin{pmatrix} 1 & 2 & 3 \\ 2 & 4 & 1 \\ 1 & 4 & 4 \end{pmatrix}.

(b)

A=(135−1423−24)\mathbf{A} = \begin{pmatrix} 1 & 3 & 5 \\ -1 & 4 & 2 \\ 3 & -2 & 4 \end{pmatrix}.

(c)

A=(12−344812−82321−3−11−4)\mathbf{A} = \begin{pmatrix} 1 & 2 & -3 & 4 \\ 4 & 8 & 12 & -8 \\ 2 & 3 & 2 & 1 \\ -3 & -1 & 1 & -4 \end{pmatrix}.

(d)

A=(00−231012−112102−30)\mathbf{A} = \begin{pmatrix} 0 & 0 & -2 & 3 \\ 1 & 0 & 1 & 2 \\ -1 & 1 & 2 & 1 \\ 0 & 2 & -3 & 0 \end{pmatrix}.

(e)

A=(2−1000−12−1000−12−1000−12−1000−11)\mathbf{A} = \begin{pmatrix} 2 & -1 & 0 & 0 & 0 \\ -1 & 2 & -1 & 0 & 0 \\ 0 & -1 & 2 & -1 & 0 \\ 0 & 0 & -1 & 2 & -1 \\ 0 & 0 & 0 & -1 & 1 \end{pmatrix}.

(f)

A=(111⋯1121⋯1113⋯1⋮⋮⋮⋱⋮111⋯n)\mathbf{A} = \begin{pmatrix} 1 & 1 & 1 & \cdots & 1 \\ 1 & 2 & 1 & \cdots & 1 \\ 1 & 1 & 3 & \cdots & 1 \\ \vdots & \vdots & \vdots & \ddots & \vdots \\ 1 & 1 & 1 & \cdots & n \end{pmatrix}.

Задача 6.1.4

Используя определители, вычислите ранг матрицы A=(13−2012−1−1625−6)\mathbf{A} = \begin{pmatrix} 1 & 3 & -2 \\ 0 & 1 & 2 \\ -1 & -1 & 6 \\ 2 & 5 & -6 \end{pmatrix}.

?
Задача 6.1.5

Используя определители, найдите значения α\alpha, при которых следующая система имеет единственное решение.

(1α001−1α01)(x1x2x3)=(−347). \begin{pmatrix} 1 & \alpha & 0 \\ 0 & 1 & -1 \\ \alpha & 0 & 1 \end{pmatrix} \begin{pmatrix} x_{1} \\ x_{2} \\ x_{3} \end{pmatrix} = \begin{pmatrix} -3 \\ 4 \\ 7 \end{pmatrix}.
?
Задача 6.1.6

Если A\mathbf{A} невырождена, объясните, почему det⁡(A−1)=1/det⁡(A)\operatorname {det}\left(\mathbf{A}^{-1}\right) = 1/\operatorname {det}\left(\mathbf{A}\right).

?
Задача 6.1.7

Объясните, почему определители инвариантны относительно преобразований подобия. То есть покажите, что det⁡(P−1AP)=det⁡(A)\operatorname {det}\left(\mathbf{P}^{-1}\mathbf{A}\mathbf{P}\right) = \operatorname {det}\left(\mathbf{A}\right) для всех невырожденных P\mathbf{P}.

?
Задача 6.1.8

Объясните, почему det⁡(A∗)=det⁡(A)‾\operatorname {det}\left(\mathbf{A}^{*}\right) = \overline{\operatorname {det}\left(\mathbf{A}\right)}.

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

Объясните, почему ∣det⁡(Q)∣=1\left|\operatorname {det}\left(\mathbf{Q}\right)\right| = 1, когда Q\mathbf{Q} унитарна. В частности, det⁡(Q)=±1\operatorname {det}\left(\mathbf{Q}\right) = \pm 1, если Q\mathbf{Q} — ортогональная матрица.

(b)

Как сингулярные числа матрицы A∈Cn×n\mathbf{A} \in \mathbb {C}^{n \times n} связаны с det⁡(A)\operatorname {det}\left(\mathbf{A}\right)?

Задача 6.1.10

Докажите, что если A\mathbf{A} имеет размер m×nm \times n, то det⁡(A∗A)≥0\operatorname {det}\left(\mathbf{A}^{*}\mathbf{A}\right) \geq 0, и объясните, почему det⁡(A∗A)>0\operatorname {det}\left(\mathbf{A}^{*}\mathbf{A}\right) > 0 тогда и только тогда, когда rk⁡(A)=n\operatorname {rk}\left(\mathbf{A}\right) = n.

?
Задача 6.1.11

Если A\mathbf{A} имеет размер n×nn \times n, объясните, почему det⁡(αA)=αndet⁡(A)\operatorname {det}\left(\alpha \mathbf{A}\right) = \alpha^{n}\operatorname {det}\left(\mathbf{A}\right) для всех скаляров α\alpha.

?
Задача 6.1.12

Если A\mathbf{A} — кососимметричная матрица размера n×nn \times n, докажите, что A\mathbf{A} вырождена всякий раз, когда nn нечётно.

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

Используйте упражнение 6.1.11.

Задача 6.1.13

Как построить случайные целочисленные матрицы с det⁡(A)=1\operatorname {det}\left(\mathbf{A}\right) = 1?

?
Задача 6.1.14

Если kk-я строка матрицы An×n\mathbf{A}_{n \times n} записана как сумма Ak∗=x→T+y→T+⋯+z→T\mathbf{A}_{k*} = \overrightarrow {x}^{T} + \overrightarrow {y}^{T} + \cdots + \overrightarrow {z}^{T}, где x→T,y→T,…,z→T\overrightarrow {x}^{T}, \overrightarrow {y}^{T}, \ldots , \overrightarrow {z}^{T} — строчные векторы, объясните, почему

det⁡(A)=det⁡(A1∗⋮x→T⋮An∗)+det⁡(A1∗⋮y→T⋮An∗)+⋯+det⁡(A1∗⋮z→T⋮An∗). \operatorname {det}\left(\mathbf{A}\right) = \det \begin{pmatrix} \mathbf{A}_{1*} \\ \vdots \\ \overrightarrow {x}^{T} \\ \vdots \\ \mathbf{A}_{n*} \end{pmatrix} + \det \begin{pmatrix} \mathbf{A}_{1*} \\ \vdots \\ \overrightarrow {y}^{T} \\ \vdots \\ \mathbf{A}_{n*} \end{pmatrix} + \cdots + \det \begin{pmatrix} \mathbf{A}_{1*} \\ \vdots \\ \overrightarrow {z}^{T} \\ \vdots \\ \mathbf{A}_{n*} \end{pmatrix}.
?
Задача 6.1.15

Неравенство КБШ (неравенство Коши--Буняковского--Шварца) утверждает, что ∣x→∗y→∣≤∥x→∥2∥y→∥2\left|\overrightarrow {x}^{*}\overrightarrow {y}\right| \leq \left\| \overrightarrow {x}\right\|_{2} \left\| \overrightarrow {y}\right\|_{2} для векторов x→,y→∈Cn×1\overrightarrow {x}, \overrightarrow {y} \in \mathbb {C}^{n \times 1}, причём равенство достигается тогда и только тогда, когда y→=αx→\overrightarrow {y} = \alpha \overrightarrow {x} при α=x→∗y→/x→∗x→\alpha = \overrightarrow {x}^{*}\overrightarrow {y}/\overrightarrow {x}^{*}\overrightarrow {x}. Используя упражнение 6.1.10, дайте альтернативное доказательство неравенства КБШ вместе с альтернативным объяснением того, почему равенство выполняется тогда и только тогда, когда y→\overrightarrow {y} является скалярным кратным x→\overrightarrow {x}.

?
Задача 6.1.16

Пусть Ak\mathbf{A}_{k} — ведущая главная подматрица размера k×kk \times k матрицы An×n\mathbf{A}_{n \times n} --- т.е. подматрица, взятая из левого верхнего угла матрицы A\mathbf{A} и состоящая из её первых kk строк и столбцов. Докажите, что если A\mathbf{A} имеет LULU-разложение A=LU\mathbf{A} = \mathbf{L}\mathbf{U}, то det⁡(Ak)=u11u22⋯ukk\operatorname {det}\left(\mathbf{A}_{k}\right) = u_{11}u_{22} \cdots u_{kk}, и выведите отсюда, что kk-й ведущий элемент равен

ukk={det⁡(A1)=a11при k=1,det⁡(Ak)/det⁡(Ak−1)при k=2,3,…,n. u_{kk} = \begin{cases} \operatorname {det}\left(\mathbf{A}_{1}\right) = a_{11} & \text{при } k = 1, \\ \operatorname {det}\left(\mathbf{A}_{k}\right)/\operatorname {det}\left(\mathbf{A}_{k-1}\right) & \text{при } k = 2,3,\ldots ,n. \end{cases}
?
Задача 6.1.17

Докажите, что если rk⁡(Am×n)=n\operatorname {rk}\left(\mathbf{A}_{m \times n}\right) = n, то ATA\mathbf{A}^{T}\mathbf{A} имеет LULU-разложение с положительными ведущими элементами --- т.е. ATA\mathbf{A}^{T}\mathbf{A} положительно определена, что означает, что симметричная матрица ATA\mathbf{A}^{T}\mathbf{A} обладает LULU-разложением, в котором каждый ведущий элемент положителен.

?
Задача 6.1.18

Пусть A(x)=(2−x3404−x−51−13−x)\mathbf{A}(x) = \begin{pmatrix} 2-x & 3 & 4 \\ 0 & 4-x & -5 \\ 1 & -1 & 3-x \end{pmatrix}.

?
(a)

Сначала вычислите det⁡(A)\operatorname {det}\left(\mathbf{A}\right), а затем найдите d(det⁡(A))/dxd\bigl(\operatorname {det}\left(\mathbf{A}\right)\bigr)/dx.

(b)

Используйте формулу (6.1.19) для вычисления d(det⁡(A))/dxd\bigl(\operatorname {det}\left(\mathbf{A}\right)\bigr)/dx.

Задача 6.1.19

Когда элементы матрицы A=[aij(x)]\mathbf{A} = [a_{ij}(x)] являются дифференцируемыми функциями от xx, определим dA/dx=[daij/dx]d\mathbf{A}/dx = [da_{ij}/dx] (матрицу производных). Для квадратных матриц всегда ли верно, что d(det⁡(A))/dx=det⁡(dA/dx)d\bigl(\operatorname {det}\left(\mathbf{A}\right)\bigr)/dx = \operatorname {det}\left(d\mathbf{A}/dx\right)?

?
Задача 6.1.20

Для набора функций S={f1(x),f2(x),…,fn(x)}\mathcal{S} = \left\{ f_{1}(x), f_{2}(x), \ldots , f_{n}(x)\right\}, каждая из которых n−1n-1 раз дифференцируема, определитель

w(x)=∣f1(x)f2(x)⋯fn(x)f1′(x)f2′(x)⋯fn′(x)⋮⋮⋱⋮f1(n−1)(x)f2(n−1)(x)⋯fn(n−1)(x)∣ w(x) = \begin{vmatrix} f_{1}(x) & f_{2}(x) & \cdots & f_{n}(x) \\ f_{1}'(x) & f_{2}'(x) & \cdots & f_{n}'(x) \\ \vdots & \vdots & \ddots & \vdots \\ f_{1}^{(n-1)}(x) & f_{2}^{(n-1)}(x) & \cdots & f_{n}^{(n-1)}(x) \end{vmatrix}

называется вронскианом набора S\mathcal{S}. Если S\mathcal{S} — линейно зависимый набор, объясните, почему w(x)=0w(x) = 0 при всех значениях xx.

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

Вспомните пример 4.3.6 (стр. 189).

Задача 6.1.21

Рассмотрим вычисление определителя размера n×nn \times n по определению (6.1.1).

?
(a)

Сколько умножений требуется?

(b)

Предполагая, что компьютер выполняет 1 000 000 умножений в секунду, и пренебрегая всеми остальными операциями, какого наибольшего порядка определитель можно вычислить за один час?

(c)

При тех же условиях, что и в пункте (б), сколько времени потребуется, чтобы вычислить определитель матрицы размера 100×100100 \times 100? Подсказка: 100!≈9.33×10157100! \approx 9.33 \times 10^{157}.

(d)

Если пренебречь всеми остальными операциями, сколько умножений в секунду должен выполнять компьютер, чтобы задача вычисления определителя матрицы размера 100×100100 \times 100 была завершена за 100 лет?

§
Задача 6.2.1

Используя разложение по алгебраическим дополнениям, вычислите каждый из следующих определителей.

?
(a)

∣211621−221∣\begin{vmatrix} 2 & 1 & 1 \\ 6 & 2 & 1 \\ -2 & 2 & 1 \end{vmatrix}.

(b)

∣00−231012−112102−30∣\begin{vmatrix} 0 & 0 & -2 & 3 \\ 1 & 0 & 1 & 2 \\ -1 & 1 & 2 & 1 \\ 0 & 2 & -3 & 0 \end{vmatrix}.

(c)

∣0111101111011110∣\begin{vmatrix} 0 & 1 & 1 & 1 \\ 1 & 0 & 1 & 1 \\ 1 & 1 & 0 & 1 \\ 1 & 1 & 1 & 0 \end{vmatrix}.

Задача 6.2.2

Используя определители, вычислите следующие обратные матрицы.

?
(a)

(211621−221)−1\begin{pmatrix} 2 & 1 & 1 \\ 6 & 2 & 1 \\ -2 & 2 & 1 \end{pmatrix}^{-1}.

(b)

(00−231012−112102−30)−1\begin{pmatrix} 0 & 0 & -2 & 3 \\ 1 & 0 & 1 & 2 \\ -1 & 1 & 2 & 1 \\ 0 & 2 & -3 & 0 \end{pmatrix}^{-1}.

Задача 6.2.3
?
(a)

Используя правило Крамера, решите систему

x1+x2+x3=1,x1+x2=α,x2+x3=β. \begin{aligned} x_{1} + x_{2} + x_{3} & = 1, \\ x_{1} + x_{2} & = \alpha , \\ x_{2} + x_{3} & = \beta . \end{aligned}
(b)

Вычислите lim⁡t→∞x2(t)\lim_{t \to \infty } x_{2}(t), где x2(t)x_{2}(t) определяется системой

x1+tx2+t2x3=t4,t2x1+x2+tx3=t3,tx1+t2x2+x3=0. \begin{aligned} x_{1} + tx_{2} + t^{2}x_{3} & = t^{4}, \\ t^{2}x_{1} + x_{2} + tx_{3} & = t^{3}, \\ tx_{1} + t^{2}x_{2} + x_{3} & = 0. \end{aligned}
Задача 6.2.4

Является ли следующее равенство корректным выводом правила Крамера для решения невырожденной системы Ax→=b→\mathbf{A}\overrightarrow {x} = \overrightarrow {b}, где Ai\mathbf{A}_{i} обозначает матрицу, совпадающую с A\mathbf{A}, за исключением того, что столбец A∗i\mathbf{A}_{*i} заменён на b→\overrightarrow {b}?

det⁡(Ai)det⁡(A)=det⁡(A−1Ai)=det⁡(e→1⋯e→i−1  x→  e→i+1⋯e→n)=xi. \dfrac {\operatorname {det}\left(\mathbf{A}_{i}\right)}{\operatorname {det}\left(\mathbf{A}\right)} = \operatorname {det}\left(\mathbf{A}^{-1}\mathbf{A}_{i}\right) = \operatorname {det}\left(\overrightarrow {e}_{1} \cdots \overrightarrow {e}_{i-1} \; \overrightarrow {x} \; \overrightarrow {e}_{i+1} \cdots \overrightarrow {e}_{n}\right) = x_{i}.
?
Задача 6.2.5
?
(a)

На примере покажите, что det⁡(A+B)≠det⁡(A)+det⁡(B)\operatorname {det}\left(\mathbf{A}+\mathbf{B}\right) \neq \operatorname {det}\left(\mathbf{A}\right)+\operatorname {det}\left(\mathbf{B}\right).

(b)

Используя квадратные матрицы, постройте пример, показывающий, что

det⁡(ABCD)≠det⁡(A)det⁡(D)−det⁡(B)det⁡(C). \det \begin{pmatrix} \mathbf{A} & \mathbf{B} \\ \mathbf{C} & \mathbf{D} \end{pmatrix} \neq \operatorname {det}\left(\mathbf{A}\right)\operatorname {det}\left(\mathbf{D}\right) - \operatorname {det}\left(\mathbf{B}\right)\operatorname {det}\left(\mathbf{C}\right).
Задача 6.2.6

Пусть rk⁡(Bm×n)=n\operatorname {rk}\left(\mathbf{B}_{m \times n}\right) = n, и пусть Q\mathbf{Q} — ортогональный проектор на N(BT)N(\mathbf{B}^{T}), т.е. Q=I−B(BTB)−1BT\mathbf{Q} = \mathbf{I} - \mathbf{B}(\mathbf{B}^{T}\mathbf{B})^{-1}\mathbf{B}^{T} — проектор, аннулирующий R(B)R(\mathbf{B}) и воспроизводящий векторы, уже принадлежащие N(BT)N(\mathbf{B}^{T}). Для A=[B∣c→n×1]\mathbf{A} = [\mathbf{B} \mid \overrightarrow {c}_{n \times 1}] докажите, что c→TQc→=det⁡(ATA)/det⁡(BTB)\overrightarrow {c}^{T}\mathbf{Q}\overrightarrow {c} = \operatorname {det}\left(\mathbf{A}^{T}\mathbf{A}\right)/\operatorname {det}\left(\mathbf{B}^{T}\mathbf{B}\right).

?
Задача 6.2.7

Если An×n\mathbf{A}_{n \times n} — невырожденная матрица, а D\mathbf{D} и C\mathbf{C} — матрицы размера n×kn \times k, объясните, как использовать (6.2.1) для вывода формулы

det⁡(A+CDT)=det⁡(A)det⁡(Ik+DTA−1C). \operatorname {det}\left(\mathbf{A}+\mathbf{C}\mathbf{D}^{T}\right) = \operatorname {det}\left(\mathbf{A}\right)\operatorname {det}\left(\mathbf{I}_{k}+\mathbf{D}^{T}\mathbf{A}^{-1}\mathbf{C}\right).
?
Примечание.
?

Это обобщение формулы обновления ранга один det⁡(A+c→d→T)=det⁡(A)(1+d→TA−1c→)\operatorname {det}\left(\mathbf{A}+\overrightarrow {c}\overrightarrow {d}^{T}\right) = \operatorname {det}\left(\mathbf{A}\right)\bigl(1+\overrightarrow {d}^{T}\mathbf{A}^{-1}\overrightarrow {c}\bigr) (6.2.3), поскольку если c→i\overrightarrow {c}_{i} и d→i\overrightarrow {d}_{i} — ii-е столбцы матриц C\mathbf{C} и D\mathbf{D} соответственно, то A+CDT=A+c→1d→1T+c→2d→2T+⋯+c→kd→kT\mathbf{A}+\mathbf{C}\mathbf{D}^{T} = \mathbf{A}+\overrightarrow {c}_{1}\overrightarrow {d}_{1}^{T}+\overrightarrow {c}_{2}\overrightarrow {d}_{2}^{T}+\cdots +\overrightarrow {c}_{k}\overrightarrow {d}_{k}^{T}.

Задача 6.2.8

Объясните, почему A\mathbf{A} вырождена тогда и только тогда, когда A[adj⁡(A)]=0\mathbf{A}[\operatorname {adj}(\mathbf{A})] = \mathbf{0}.

?
Задача 6.2.9

Для невырожденной линейной системы Ax→=b→\mathbf{A}\overrightarrow {x} = \overrightarrow {b} объясните, почему каждая компонента решения должна непрерывно зависеть от элементов A\mathbf{A}.

?
Задача 6.2.10

Для скаляров α\alpha объясните, почему adj⁡(αA)=αn−1adj⁡(A)\operatorname {adj}(\alpha \mathbf{A}) = \alpha^{n-1}\operatorname {adj}(\mathbf{A}).

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

Вспомните упражнение 6.1.11.

Задача 6.2.11

Для матрицы A\mathbf{A} размера n×nn \times n докажите, что следующие утверждения верны.

?
(a)

Если rk⁡(A)<n−1\operatorname {rk}\left(\mathbf{A}\right) < n-1, то adj⁡(A)=0\operatorname {adj}(\mathbf{A}) = \mathbf{0}.

(b)

Если rk⁡(A)=n−1\operatorname {rk}\left(\mathbf{A}\right) = n-1, то rk⁡(adj⁡(A))=1\operatorname {rk}\left(\operatorname {adj}(\mathbf{A})\right) = 1.

(c)

Если rk⁡(A)=n\operatorname {rk}\left(\mathbf{A}\right) = n, то rk⁡(adj⁡(A))=n\operatorname {rk}\left(\operatorname {adj}(\mathbf{A})\right) = n.

Задача 6.2.12

В 1812 году Коши открыл формулу, утверждающую, что если A\mathbf{A} имеет размер n×nn \times n, то det⁡(adj⁡(A))=[det⁡(A)]n−1\operatorname {det}\left(\operatorname {adj}(\mathbf{A})\right) = [\operatorname {det}\left(\mathbf{A}\right)]^{n-1}. Установите формулу Коши.

?
Задача 6.2.13

Для следующей трёхдиагональной матрицы An\mathbf{A}_{n} положите Dn=det⁡(An)D_{n} = \operatorname {det}\left(\mathbf{A}_{n}\right) и выведите формулу Dn=2Dn−1−Dn−2D_{n} = 2D_{n-1}-D_{n-2}, чтобы заключить, что Dn=n+1D_{n} = n+1.

An=(2−10⋯0−12−1⋯0⋱⋱⋱0⋯−12−10⋯0−12)n×n. \mathbf{A}_{n} = \begin{pmatrix} 2 & -1 & 0 & \cdots & 0 \\ -1 & 2 & -1 & \cdots & 0 \\ & \ddots & \ddots & \ddots & \\ 0 & \cdots & -1 & 2 & -1 \\ 0 & \cdots & 0 & -1 & 2 \end{pmatrix}_{n \times n}.
?
Задача 6.2.14

Рассматривая матрицы с обновлением ранга один, выведите следующие формулы.

?
(a)

∣1+α1α11⋯111+α2α2⋯1⋮⋮⋱⋮11⋯1+αnαn∣=1+∑αi∏αi\begin{vmatrix} \frac{1+\alpha_{1}}{\alpha_{1}} & 1 & \cdots & 1 \\ 1 & \frac{1+\alpha_{2}}{\alpha_{2}} & \cdots & 1 \\ \vdots & \vdots & \ddots & \vdots \\ 1 & 1 & \cdots & \frac{1+\alpha_{n}}{\alpha_{n}} \end{vmatrix} = \dfrac {1+\sum \alpha_{i}}{\prod \alpha_{i}}.

(b)

∣αββ⋯ββαβ⋯βββα⋯β⋮⋮⋮⋱⋮βββ⋯α∣n×n={(α−β)n(1+nβα−β)если α≠β,0если α=β.\begin{vmatrix} \alpha & \beta & \beta & \cdots & \beta \\ \beta & \alpha & \beta & \cdots & \beta \\ \beta & \beta & \alpha & \cdots & \beta \\ \vdots & \vdots & \vdots & \ddots & \vdots \\ \beta & \beta & \beta & \cdots & \alpha \end{vmatrix}_{n \times n} = \begin{cases} (\alpha -\beta )^{n}\left(1+\dfrac {n\beta }{\alpha -\beta }\right) & \text{если } \alpha \neq \beta , \\ 0 & \text{если } \alpha = \beta . \end{cases}

(c)

∣1+α1α2⋯αnα11+α2⋯αn⋮⋮⋱⋮α1α2⋯1+αn∣=1+α1+α2+⋯+αn\begin{vmatrix} 1+\alpha_{1} & \alpha_{2} & \cdots & \alpha_{n} \\ \alpha_{1} & 1+\alpha_{2} & \cdots & \alpha_{n} \\ \vdots & \vdots & \ddots & \vdots \\ \alpha_{1} & \alpha_{2} & \cdots & 1+\alpha_{n} \end{vmatrix} = 1+\alpha_{1}+\alpha_{2}+\cdots +\alpha_{n}.

Задача 6.2.15

Окаймлённая матрица имеет вид B=(Ax→y→Tα)\mathbf{B} = \begin{pmatrix} \mathbf{A} & \overrightarrow {x} \\ \overrightarrow {y}^{T} & \alpha \end{pmatrix}, где An×n\mathbf{A}_{n \times n} невырождена, x→\overrightarrow {x} — столбец, y→T\overrightarrow {y}^{T} — строка, а α\alpha — скаляр. Объясните, почему следующие утверждения должны быть верны.

?
(a)

∣Ax→y→T−1∣=−det⁡(A+x→y→T)\begin{vmatrix} \mathbf{A} & \overrightarrow {x} \\ \overrightarrow {y}^{T} & -1 \end{vmatrix} = -\operatorname {det}\left(\mathbf{A}+\overrightarrow {x}\overrightarrow {y}^{T}\right).

(b)

∣Ax→y→T0∣=−y→Tadj⁡(A)x→\begin{vmatrix} \mathbf{A} & \overrightarrow {x} \\ \overrightarrow {y}^{T} & 0 \end{vmatrix} = -\overrightarrow {y}^{T}\operatorname {adj}(\mathbf{A})\overrightarrow {x}.

Задача 6.2.16

Если B\mathbf{B} имеет размер m×nm \times n, а C\mathbf{C} — размер n×mn \times m, объясните, почему (6.2.1) гарантирует, что λmdet⁡(λIn−CB)=λndet⁡(λIm−BC)\lambda^{m}\operatorname {det}\left(\lambda \mathbf{I}_{n}-\mathbf{C}\mathbf{B}\right) = \lambda^{n}\operatorname {det}\left(\lambda \mathbf{I}_{m}-\mathbf{B}\mathbf{C}\right) верно для всех скаляров λ\lambda.

?
Задача 6.2.17

Для квадратной матрицы A\mathbf{A} и векторов-столбцов c→\overrightarrow {c} и d→\overrightarrow {d} выведите следующие два обобщения формулы обновления ранга 11 det⁡(A+c→d→T)=det⁡(A)(1+d→TA−1c→)\operatorname {det}\left(\mathbf{A}+\overrightarrow {c}\overrightarrow {d}^{T}\right) = \operatorname {det}\left(\mathbf{A}\right)\bigl(1+\overrightarrow {d}^{T}\mathbf{A}^{-1}\overrightarrow {c}\bigr) (6.2.3).

?
(a)

Если Ax→=c→\mathbf{A}\overrightarrow {x} = \overrightarrow {c}, то det⁡(A+c→d→T)=det⁡(A)(1+d→Tx→)\operatorname {det}\left(\mathbf{A}+\overrightarrow {c}\overrightarrow {d}^{T}\right) = \operatorname {det}\left(\mathbf{A}\right)\bigl(1+\overrightarrow {d}^{T}\overrightarrow {x}\bigr).

(b)

Если y→TA=d→T\overrightarrow {y}^{T}\mathbf{A} = \overrightarrow {d}^{T}, то det⁡(A+c→d→T)=det⁡(A)(1+y→Tc→)\operatorname {det}\left(\mathbf{A}+\overrightarrow {c}\overrightarrow {d}^{T}\right) = \operatorname {det}\left(\mathbf{A}\right)\bigl(1+\overrightarrow {y}^{T}\overrightarrow {c}\bigr).

Задача 6.2.18

Опишите определитель элементарного отражателя R=I−2u→u→∗/u→∗u→\mathbf{R} = \mathbf{I} - 2\overrightarrow {u}\overrightarrow {u}^{*}/\overrightarrow {u}^{*}\overrightarrow {u} (отражения относительно гиперплоскости u→⊥\overrightarrow {u}^{\perp }) и плоского вращения Pij\mathbf{P}_{ij} (вращения в плоскости (i,j)(i,j), использующего cos⁡θ\cos \theta и sin⁡θ\sin \theta в позициях (i,i),(i,j),(j,i),(j,j)(i,i),(i,j),(j,i),(j,j) в остальном единичной матрицы), а затем объясните, как найти det⁡(A)\operatorname {det}\left(\mathbf{A}\right) с помощью приведения Хаусхолдера (приведения A\mathbf{A} к верхней трапециевидной форме последовательностью элементарных отражателей) и приведения Гивенса (того же приведения с помощью последовательности плоских вращений вместо этого).

?
Задача 6.2.19

Предположим, что A\mathbf{A} — невырожденная матрица, все элементы которой являются целыми числами. Докажите, что элементы A−1\mathbf{A}^{-1} являются целыми числами тогда и только тогда, когда det⁡(A)=±1\operatorname {det}\left(\mathbf{A}\right) = \pm 1.

?
Задача 6.2.20

Пусть A=I−2u→v→T\mathbf{A} = \mathbf{I} - 2\overrightarrow {u}\overrightarrow {v}^{T} — матрица, в которой u→\overrightarrow {u} и v→\overrightarrow {v} — векторы-столбцы с целыми элементами.

?
(a)

Докажите, что A−1\mathbf{A}^{-1} имеет целые элементы тогда и только тогда, когда v→Tu→=0\overrightarrow {v}^{T}\overrightarrow {u} = 0 или 11.

(b)

Матрица называется инволютивной, если A−1=A\mathbf{A}^{-1} = \mathbf{A}. Объясните, почему A=I−2u→v→T\mathbf{A} = \mathbf{I}-2\overrightarrow {u}\overrightarrow {v}^{T} является инволютивной, когда v→Tu→=1\overrightarrow {v}^{T}\overrightarrow {u} = 1.

Задача 6.2.21

С помощью индукции покажите, что разложение det⁡(An×n)\operatorname {det}\left(\mathbf{A}_{n \times n}\right) по алгебраическим дополнениям требует

c(n)=n! (1+12!+13!+⋯+1(n−1)!) c(n) = n!\, \left(1+\dfrac {1}{2!}+\dfrac {1}{3!}+\cdots +\dfrac {1}{(n-1)!}\right)

умножений при n≥2n \geq 2. Предположив, что компьютер выполняет 1 000 000 умножений в секунду, и пренебрегая всеми остальными операциями, оцените, сколько времени потребуется для вычисления определителя матрицы 100×100100 \times 100 с помощью разложения по алгебраическим дополнениям.

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

Вспомните разложение в ряд для exe^{x} и используйте 100!≈9,33×10157100! \approx 9{,}33 \times 10^{157}.

Задача 6.2.22

Найдите все значения λ\lambda, при которых матрица A−λI\mathbf{A}-\lambda \mathbf{I} вырождена, где

A=(0−3−2252−2−30). \mathbf{A} = \begin{pmatrix} 0 & -3 & -2 \\ 2 & 5 & 2 \\ -2 & -3 & 0 \end{pmatrix}.
?
Примечание.
?

Если p(λ)=λn+αn−1λn−1+⋯+α1λ+α0p(\lambda ) = \lambda^{n}+\alpha_{n-1}\lambda^{n-1}+\cdots +\alpha_{1}\lambda +\alpha_{0} — приведённый многочлен с целыми коэффициентами, то целые корни p(λ)p(\lambda ) образуют подмножество делителей α0\alpha_{0}.

Задача 6.2.23

Предположим, что f1(t),f2(t),…,fn(t)f_{1}(t), f_{2}(t), \ldots , f_{n}(t) — решения линейного дифференциального уравнения nn-го порядка y(n)+p1(t)y(n−1)+⋯+pn−1(t)y′+pn(t)y=0y^{(n)}+p_{1}(t)y^{(n-1)}+\cdots +p_{n-1}(t)y'+p_{n}(t)y = 0, и пусть w(t)w(t) — вронскиан

w(t)=∣f1(t)f2(t)⋯fn(t)f1′(t)f2′(t)⋯fn′(t)⋮⋮⋱⋮f1(n−1)(t)f2(n−1)(t)⋯fn(n−1)(t)∣. w(t) = \begin{vmatrix} f_{1}(t) & f_{2}(t) & \cdots & f_{n}(t) \\ f_{1}'(t) & f_{2}'(t) & \cdots & f_{n}'(t) \\ \vdots & \vdots & \ddots & \vdots \\ f_{1}^{(n-1)}(t) & f_{2}^{(n-1)}(t) & \cdots & f_{n}^{(n-1)}(t) \end{vmatrix}.

Преобразовав уравнение nn-го порядка в систему nn уравнений первого порядка с помощью подстановок x1=y,x2=y′,…,xn=y(n−1)x_{1} = y, x_{2} = y', \ldots , x_{n} = y^{(n-1)}, покажите, что w(t)=w(ξ0) e−∫ξ0tp1(ξ) dξw(t) = w(\xi_{0})\, e^{-\int_{\xi_{0}}^{t} p_{1}(\xi )\, d\xi } при произвольной постоянной ξ0\xi_{0}.

?
Задача 6.2.24

Вычислите определитель Вандермонда, показав, что

∣1x1x12⋯x1n−11x2x22⋯x2n−1⋮⋮⋮⋯⋮1xnxn2⋯xnn−1∣=∏j>i(xj−xi). \begin{vmatrix} 1 & x_{1} & x_{1}^{2} & \cdots & x_{1}^{n-1} \\ 1 & x_{2} & x_{2}^{2} & \cdots & x_{2}^{n-1} \\ \vdots & \vdots & \vdots & \cdots & \vdots \\ 1 & x_{n} & x_{n}^{2} & \cdots & x_{n}^{n-1} \end{vmatrix} = \prod _{j>i} (x_{j}-x_{i}).

Когда это выражение отлично от нуля (сравните с примером 4.3.4)?

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

Для многочлена p(λ)=∣1λλ2⋯λk−11x2x22⋯x2k−1⋮⋮⋮⋯⋮1xkxk2⋯xkk−1∣k×kp(\lambda ) = \begin{vmatrix} 1 & \lambda & \lambda^{2} & \cdots & \lambda^{k-1} \\ 1 & x_{2} & x_{2}^{2} & \cdots & x_{2}^{k-1} \\ \vdots & \vdots & \vdots & \cdots & \vdots \\ 1 & x_{k} & x_{k}^{2} & \cdots & x_{k}^{k-1} \end{vmatrix}_{k \times k} с помощью индукции найдите степень p(λ)p(\lambda ), корни p(λ)p(\lambda ) и коэффициент при λk−1\lambda^{k-1} в p(λ)p(\lambda ).

Задача 6.2.25

Предположим, что каждый элемент An×n=[aij(x)]\mathbf{A}_{n \times n} = [a_{ij}(x)] является дифференцируемой функцией вещественной переменной xx. Используя формулу (6.1.19), выведите формулу

d(det⁡(A))dx=∑j=1n∑i=1ndaijdx A˚ij. \dfrac {d\bigl(\operatorname {det}\left(\mathbf{A}\right)\bigr)}{dx} = \sum _{j=1}^{n} \sum _{i=1}^{n} \dfrac {da_{ij}}{dx}\, \mathring {A}_{ij}.
?
Задача 6.2.26

Рассматривая элементы A\mathbf{A} как независимые переменные и используя формулу (6.1.19), выведите формулу

∂det⁡(A)∂aij=A˚ij. \dfrac {\partial \operatorname {det}\left(\mathbf{A}\right)}{\partial a_{ij}} = \mathring {A}_{ij}.
?
Задача 6.2.27

В 1772 году французский математик Пьер-Симон Лаплас (1749--1827) представил следующее обобщённое разложение по алгебраическим дополнениям. Для матрицы A\mathbf{A} размера n×nn \times n пусть A(i1i2⋯ik∣j1j2⋯jk)\mathbf{A}(i_{1}i_{2} \cdots i_{k} \mid j_{1}j_{2} \cdots j_{k}) обозначает k×kk \times k подматрицу A\mathbf{A}, лежащую на пересечении строк i1,i2,…,iki_{1},i_{2},\ldots ,i_{k} со столбцами j1,j2,…,jkj_{1},j_{2},\ldots ,j_{k}, а M(i1i2⋯ik∣j1j2⋯jk)M(i_{1}i_{2} \cdots i_{k} \mid j_{1}j_{2} \cdots j_{k}) обозначает минор порядка n−k×n−kn-k \times n-k, получаемый удалением строк i1,i2,…,iki_{1},i_{2},\ldots ,i_{k} и столбцов j1,j2,…,jkj_{1},j_{2},\ldots ,j_{k} из A\mathbf{A}. Обобщённым алгебраическим дополнением A(i1⋯ik∣j1⋯jk)\mathbf{A}(i_{1} \cdots i_{k} \mid j_{1} \cdots j_{k}) называется знаковый минор A˚(i1⋯ik∣j1⋯jk)=(−1)i1+⋯+ik+j1+⋯+jkM(i1⋯ik∣j1⋯jk)\mathring {A}(i_{1} \cdots i_{k} \mid j_{1} \cdots j_{k}) = (-1)^{i_{1}+\cdots +i_{k}+j_{1}+\cdots +j_{k}}M(i_{1} \cdots i_{k} \mid j_{1} \cdots j_{k}), и для каждого фиксированного набора индексов строк 1≤i1<⋯<ik≤n1 \leq i_{1} < \cdots < i_{k} \leq n

det⁡(A)=∑1≤j1<⋯<jk≤ndet⁡A(i1⋯ik∣j1⋯jk) A˚(i1⋯ik∣j1⋯jk). \operatorname {det}\left(\mathbf{A}\right) = \sum _{1 \leq j_{1} < \cdots < j_{k} \leq n} \det \mathbf{A}(i_{1} \cdots i_{k} \mid j_{1} \cdots j_{k})\, \mathring {A}(i_{1} \cdots i_{k} \mid j_{1} \cdots j_{k}).

Каждая такая сумма содержит (nk)\binom {n}{k} слагаемых. Используя разложение Лапласа, вычислите определитель матрицы

A=(00−231012−112102−30) \mathbf{A} = \begin{pmatrix} 0 & 0 & -2 & 3 \\ 1 & 0 & 1 & 2 \\ -1 & 1 & 2 & 1 \\ 0 & 2 & -3 & 0 \end{pmatrix}

по первой и третьей строкам.

?