14.1.7

[НЕТ ЗАГОЛОВКА]

[24/42%]
Показать
LaTeX
Задача 1585

(Теорема Фредгольма). Неоднородная система линейных уравнений

Ax=A[x1x2⋮xn]=[b1b2⋮bm]=b A x = A\left[\begin{smallmatrix} x_{1} \\ x_{2} \\ \vdots \\ x_{n} \end{smallmatrix}\right] =\left[\begin{smallmatrix} b_{1} \\ b_{2} \\ \vdots \\ b_{m} \end{smallmatrix}\right] = b

совместна тогда и только тогда, когда вектор-столбец b∈Rmb \in \mathbb {R}^{m} ортогонален всем решениям однородной системы A⊤x=0A^{\top } x = 0 с транспонированной матрицей A⊤A^{\top }.

?
Задача 1586

При помощи теоремы Фредгольма исследовать совместность следующих систем линейных неоднородных уравнений:

?
(1)

{2x1+x2+3x3=109,−x1+x2=−50,2x1−x2+x3=103;\begin{cases} 2 x_{1}+x_{2}+3 x_{3} = 109, \\ -x_{1}+x_{2} = -50, \\ 2 x_{1}-x_{2}+x_{3} = 103 ; \end{cases}

(2)

{x1+x2+x3+x4=1317,2x1−x2−x3+2x4=1517,x1+x4=119.\begin{cases} x_{1}+x_{2}+x_{3}+x_{4} = \frac{13}{17}, \\ 2 x_{1}-x_{2}-x_{3}+2 x_{4} = \frac{15}{17}, \\ x_{1}+x_{4} = \frac{1}{19}. \end{cases}

Задача 1587

(Альтернатива Фредгольма). Доказать, что либо неоднородная система линейных уравнений Ax=bA x = b совместна при любой правой части bb, либо однородная система A⊤x=0A^{\top } x = 0 имеет ненулевые решения.

?
Задача 1588

При помощи теоремы Фредгольма исследовать совместность следующих систем линейных неоднородных уравнений при различных значениях λ\lambda :

?
(1)

{x1+x2−x3−x4=λ,x2+x3+2x4=1−λ,x1+x3−x4=1+3λ,3x1+3x2−x4=6;\begin{cases} x_{1}+x_{2}-x_{3}-x_{4} = \lambda , \\ x_{2}+x_{3}+2 x_{4} = 1-\lambda , \\ x_{1}+x_{3}-x_{4} = 1+3 \lambda , \\ 3 x_{1}+3 x_{2}-x_{4} = 6 ; \end{cases}

(2)

{λx1+x2+x3=7,x1+λx2+x3=7,x1+x2+λx3=7.\begin{cases} \lambda x_{1}+x_{2}+x_{3} = 7, \\ x_{1}+\lambda x_{2}+x_{3} = 7, \\ x_{1}+x_{2}+\lambda x_{3} = 7. \end{cases}

Обсуждая в гл. 13 метод наименьших квадратов, мы накладывали некоторое ограничение на ранг матрицы AA несовместной системы линейных уравнений Ax=bA x = b, а именно: система столбцов матрицы AA должна была быть линейно независимой. В общем же случае предлагалось выделить базис системы столбцов и применить к нему метод наименьших квадратов. При этом у исходной системы Ax=bA x = b возникало целое семейство псевдорешений. Уточним теперь некоторые детали.

Задача 1589

Доказать, что множество псевдорешений системы Ax=bA x = b образует аффинное подпространство (линейное многообразие)

P=x′+Ker⁡A P = x^{\prime }+\operatorname {Ker} A

в Rn\mathbb {R}^{n} (x′x^{\prime } — некоторое псевдорешение).

?
Задача 1590

Пусть P=x′+Ker⁡AP = x^{\prime }+\operatorname {Ker} A — аффинное подпространство псевдорешений системы Ax=bA x = b, а x′′∈Px^{\prime \prime } \in P — некоторое псевдорешение. Доказать, что следующие утверждения эквивалентны:

?
(1)

x′′⊥Ker⁡Ax^{\prime \prime } \perp \operatorname {Ker} A;

(2)

∣x′′∣=min⁡x∈P∣x∣\left|x^{\prime \prime }\right| = \min_{x \in P}\left|x\right|;

(3)

x′′∈Im⁡A∗x^{\prime \prime } \in \operatorname {Im} A^{*}.

Доказать, что существует единственный вектор x′′x^{\prime \prime }, определенный условиями 1-3. Он называется нормальным псевдорешением.

Псевдорешения системы Ax=bA x = b можно находить при помощи сингулярных базисов оператора А (см. задачу 1460).

Задача 1591

Рассмотрим систему Ax=bA x = b, где AA-матрица оператора A:Rn→Rm\mathrm{A}: \mathbb {R}^{n} \rightarrow \mathbb {R}^{m}. Пространства Rn\mathbb {R}^{n} и Rm\mathbb {R}^{m} снабжены стандартными скалярными произведениями. Пусть e1,e2,…,ene_{1}, e_{2}, \ldots , e_{n} — канонический базис неотрицательного самосопряженного оператора A∗A\mathbf{A}^{*} \mathbf{A}, а λ1,λ2,…,λr\lambda_{1}, \lambda_{2}, \ldots , \lambda_{r} — его положительные (т. о. ненулевые) собственные значения, отвечающие собственным векторам e1,e2,…,ere_{1}, e_{2}, \ldots , e_{r} соответственно. Доказать, что общее псевдорешение x′x^{\prime } системы Ax=bA x = b имеет вид

x′=(A∗b,e1)λ1e1+…+(A∗b,er)λrer+αr+1er+1+…+αnen x^{\prime } = \frac{\left(\mathbf{A}^{*} b, e_{1}\right)}{\lambda _{1}} e_{1}+\ldots +\frac{\left(\mathbf{A}^{*} b, e_{r}\right)}{\lambda _{r}} e_{r}+\alpha _{r+1} e_{r+1}+\ldots +\alpha _{n} e_{n}

причем нормальное псевдорешение

xnorm ′=(A∗b,e1)λ1e1+…+(A∗b,er)λrer x_{\text{norm }}^{\prime } = \frac{\left(\mathbf{A}^{*} b, e_{1}\right)}{\lambda _{1}} e_{1}+\ldots +\frac{\left(\mathbf{A}^{*} b, e_{r}\right)}{\lambda _{r}} e_{r}

Если мы имеем несовместную систему Ax=bA x = b линейных неоднородных уравнений, то с точки зрения теории линейных операторов речь идет о ситуации, когда оператор A\mathbf{A} необратим, а вектор bb не лежит в его образе Im A. С помощью скалярного произведения удается определить так называемый псевдообратный оператор A+\mathbf{A}^{+}, который будет «наилучшим образом заменять» несуществующий обратный оператор A−1\mathbf{A}^{-1}.

?
Задача 1592

Дан произвольный оператор A:V→W\mathbf{A}: V \rightarrow W, где VV и WW — евклидовы пространства размерностей nn и mm. Выберем в них ортогональные дополнения к ядру Ker A и образу Im A оператора A соответственно:

(Ker⁡A)⊥⊕Ker⁡A=V (\operatorname {Ker} \mathbf{A})^{\perp } \oplus \operatorname {Ker} \mathbf{A} = V

Определим оператор

A~=A∣(Ker⁡A)⊥:(Ker⁡A)⊥→Im⁡A,A~(x)=A(x)∀x∈(Ker⁡A)⊥ \widetilde{\mathbf{A}} = \left.\mathbf{A}\right|_{(\operatorname {Ker} \mathbf{A})^{\perp }}:(\operatorname {Ker} \mathbf{A})^{\perp } \rightarrow \operatorname {Im} \mathbf{A}, \quad \widetilde{\mathbf{A}}(x) = \mathbf{A}(x) \quad \forall x \in (\operatorname {Ker} \mathbf{A})^{\perp }

Доказать, что A~\widetilde{\mathrm{A}} — обратимый оператор. Заметим, что A~\tilde{\mathrm{A}} есть оператор не только с иной областью определения, но и с другой областью значений, чем у исходного оператора А.

?
Задача 1593

Пусть матрица AA системы уравнений определяет некоторое линейное отображение A:Rn→Rm\mathbf{A}: \mathbb {R}^{n} \rightarrow \mathbb {R}^{m}. Пространства Rn\mathbb {R}^{n} и Rm\mathbb {R}^{m} можно считать евклидовыми, снабдив их соответствующими (например стандартными) скалярными произведениями. Определим псевдообратный оператор A+\mathrm{A}^{+}следующим образом:

A+(x)=A~−1(x),x∈Im⁡AA+(x)=0,x∈(Im⁡A)⊥ \begin{array}{ll} \mathbf{A}^{+}(x) = \widetilde{\mathbf{A}}^{-1}(x), & x \in \operatorname {Im} \mathbf{A} \\ \mathbf{A}^{+}(x) = 0, & x \in (\operatorname {Im} \mathbf{A})^{\perp } \end{array}

где A~=A∣(Кке А) \widetilde{A} = \left.\mathrm{A}\right|_{(\text{Кке А) }} — оператор, построенный в предыдущей задаче. Доказать, что:

?
(1)

если оператор А является обратимым, то A+=A−1\mathbf{A}^{+} = \mathrm{A}^{-1};

(2)

вектор x′=A+bx^{\prime } = \mathrm{A}^{+} b представляет собой нормальное псевдорешение системы Ax=bA x = b.

Задача 1594

Доказать, что ядро и образ псевдообратного оператора A+\mathbf{A}^{+} совпадают соответственно с ядром и образом сопряженного A∗\mathrm{A}^{*} оператора.

?
Задача 1595

Пусть e1,e2,…,en∈Ve_{1}, e_{2}, \ldots , e_{n} \in V и e~1,e~2,…,e~m∈W\widetilde{e}_{1}, \widetilde{e}_{2}, \ldots , \widetilde{e}_{m} \in W — сингулярные базисы (см. задачу 1460) оператора A:V→W\mathrm{A}: V \rightarrow W. Найти матрицу псевдообратного оператора A+\mathbf{A}^{+}относительно этих двух базисов.

?
Задача 1596

Пусть A:Rn→Rm\mathrm{A}: \mathbb {R}^{n} \rightarrow \mathbb {R}^{m} — линейное отображение евклидовых пространств, заданное матрицей AA. Определим отображение B:Rm→Rn\mathbf{B}: \mathbb {R}^{m} \rightarrow \mathbb {R}^{n} следующим образом: пусть y∈Rmy \in \mathbb {R}^{m}, тогда значение В (y)∈Rn(y) \in \mathbb {R}^{n} определим как (единственное) нормальное псевдорешение системы Ax=yA x = y. Доказать, что:

?
(1)

В-линейный оператор;

(2)

В совпадает с псевдообратным оператором A+\mathbf{A}^{+}.

Задача 1597

Найти псевдообратный оператор A+\mathrm{A}^{+}для оператора A:R4→R3\mathrm{A}: \mathbb {R}^{4} \rightarrow \mathbb {R}^{3}, заданного матрицей

[1−120−12−3101−11] \left[\begin{smallmatrix} 1 & -1 & 2 & 0 \\ -1 & 2 & -3 & 1 \\ 0 & 1 & -1 & 1 \end{smallmatrix}\right]
?
Задача 1598

Доказать следующие свойства псевдообратных операторов:

?
(1)

(A+)+=A\left(\mathbf{A}^{+}\right)^{+} = \mathbf{A};

(2)

(A∗)+=(A+)∗\left(\mathrm{A}^{*}\right)^{+} = \left(\mathrm{A}^{+}\right)^{*};

(3)

(AA+)∗=AA+\left(\mathrm{AA}^{+}\right)^{*} = \mathrm{AA}^{+}, (AA+)2=AA+\left(\mathrm{AA}^{+}\right)^{2} = \mathrm{AA}^{+}

(4)

(A+A)∗=A+A\left(A^{+} A\right)^{*} = A^{+} A, (A+A)2=A+A\left(\mathbf{A}^{+} \mathbf{A}\right)^{2} = \mathbf{A}^{+} \mathbf{A}.

Задача 1599

Пусть столбцы a1,a2,…,ana_{1}, a_{2}, \ldots , a_{n} матрицы AA оператора A:Rn→Rm\mathbf{A}: \mathbb {R}^{n} \rightarrow \mathbb {R}^{m} линейно независимы. Доказать, что матрица A+A^{+}псевдообратного оператора A+\mathbf{A}^{+}может быть найдена следующим образом: ii-й столбец ai+,i=1,…,ma_{i}^{+}, i = 1, \ldots , m, матрицы A+A^{+}есть псевдорешение системы линейных уравнений

A[x1⋮xn]=[0⋮1⋮0]=ei A\left[\begin{smallmatrix} x_{1} \\ \vdots \\ x_{n} \end{smallmatrix}\right] =\left[\begin{smallmatrix} 0 \\ \vdots \\ 1 \\ \vdots \\ 0 \end{smallmatrix}\right] = e_{i}

(здесь ei−ie_{i}-i-й вектор стандартного базиса Rm\mathbb {R}^{m}), найденное методом наименьших квадратов.

?
Задача 1600

Найти псевдообратный оператор A+\mathrm{A}^{+}для оператора A:R3→R4\mathrm{A}: \mathbb {R}^{3} \rightarrow \mathbb {R}^{4}, заданного матрицей

[1−10−12120−1011] \left[\begin{smallmatrix} 1 & -1 & 0 \\ -1 & 2 & 1 \\ 2 & 0 & -1 \\ 0 & 1 & 1 \end{smallmatrix}\right]
?
Задача 1601

Доказать, что матрица A+A^{+}псевдообратного оператора A+\mathbf{A}^{+} (определенного относительно стандартного скалярного произведения) равна нормальному псевдорешению матричного уравнения AX=EA X = E (рассмотренного как неоднородная система линейных уравнений с nmn m неизвестными).

?
Задача 1602

Пусть ранг матрицы AA оператора A:Rn→Rm\mathbf{A}: \mathbb {R}^{n} \rightarrow \mathbb {R}^{m} равен nn. Доказать, что матрица A+A^{+}псевдообратного оператора A+\mathbf{A}^{+}(определенного относительно стандартного скалярного произведения) может быть вычислена по следующей формуле:

A+=(A⊤A)−1A⊤ A^{+} = \left(A^{\top } A\right)^{-1} A^{\top }

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

?
Задача 1603

Пусть e1,e2,…,ene_{1}, e_{2}, \ldots , e_{n} — ортонормированный базис евклидова пространства Rn\mathbb {R}^{n}. Назовем элементарным вращением Aij(φ)\mathbf{A}_{i j}(\varphi ) поворот на угол φ\varphi в двумерном подпространстве с базисом ei,eje_{i}, e_{j}.

Найти матрицу оператора Aij(φ)\mathbf{A}_{i j}(\varphi ) в базисе e1,e2,…,ene_{1}, e_{2}, \ldots , e_{n}. Доказать, что для любого вектора v∈Rnv \in \mathbb {R}^{n} существует набор из nn - 1 элементарных вращений таких, что их композиция А переводит вектор vv в вектор ∣v∣e1=(∣v∣,0,…,0)\left|v\right| e_{1} = (\left|v\right|, 0, \ldots , 0).

?
Задача 1604

Доказать, что любая ортогональная матрица AA может быть представлена в виде произведения N=n(n−1)2N = \frac{n(n-1)}{2} матриц AiNjN(φN),…A_{i_{N} j_{N}}\left(\varphi_{N}\right), \ldots..., Ai1j1(φ1)A_{i_{1} j_{1}}\left(\varphi_{1}\right) элементарных вращений и диагональной матрицы AdA_{d}, у которой на диагонали стоят {±1}\left\{ \pm 1\right\} :

A=AiNjN(φN)…Ai1j1(φ1)Ad A = A_{i_{N} j_{N}}\left(\varphi _{N}\right) \ldots A_{i_{1} j_{1}}\left(\varphi _{1}\right) A_{d}
?
Задача 1605

Рассмотрим подпространство LL евклидова (эрмитова) пространства Rn(Cn)\mathbb {R}^{n}\left(\mathbb {C}^{n}\right) размерности n−1n-1, заданное уравнением (x,w)=0(x, w) = 0, где ww — некоторый вектор Rn(Cn)\mathbb {R}^{n}\left(\mathbb {C}^{n}\right), причем ∣w∣=1\left|w\right| = 1. Пусть A\mathrm{A} — оператор ортогонального отражения относительно LL. Доказать, что матрица оператора А имеет вид

A=E−2ww∗ A = E-2 w w^{*}

где ww — столбец координат вектора ww, рассмотренный как матрица размера n×1n \times 1. Операторы такого вида будем называть операторами элементарного отражения.

?
Задача 1606

Пусть xx — произвольный вектор единичной длины пространства Rn(Cn)\mathbb {R}^{n}\left(\mathbb {C}^{n}\right). Найти такой оператор отражения относительно гиперплоскости LL (см. предыдущую задачу), т. е. найти такой вектор ww, что

(E−2ww∗)(x)=e1 \left(E-2 w w^{*}\right)(x) = e_{1}

Где геометрически расположено соответствующее подпространство LL?

?
Задача 1607

Доказать, что любая унитарная (ортогональная) матрица может быть представлена в виде произведения n−1n-1 матриц элементарных отражений и диагональной унитарной (ортогональной) матрицы.

?
Задача 1608

Придумать алгоритмы для представления произвольной вещественной квадратной матрицы AA в виде произведения A=URA = U R ортогональной матрицы UU на верхнетреугольную матрицу RR с помощью:

?
(1)

элементарных вращений Aij(φ)A_{i j}(\varphi );

(2)

элементарных отражений E−2ww∗E-2 w w^{*}.

Как с помощью такого представления A=URA = U R решать системы неоднородных линейных уравнений?