Повторить

Ортогональное дополнение. Метод наименьших квадратов

Прогресс изучения0 из 5

Отмечай теоремы кнопками ✓ (изучено) и ↻ (повторить)

Все теоремы темы5
  1. Свойства ортогонального дополнения
  2. Ортогональная проекция и расстояние до подпространства
  3. Изоморфизм конечномерных евклидовых пространств
  4. Метод наименьших квадратов. Нормальная система
  5. Нормальное псевдорешение и псевдообратная матрица
теорема

Свойства ортогонального дополнения

Формулировка

Для подпространства LL nn-мерного евклидова пространства VV:

  1. L⊥L^\perp — подпространство и dim⁡L⊥=n−dim⁡L\dim L^\perp = n - \dim L;
  2. V=L⊕L⊥V = L\oplus L^\perp;
  3. (L⊥)⊥=L(L^\perp)^\perp = L;
  4. если L=span⁡(a1,…,ak)L = \operatorname{span}(a_1,\dots,a_k), то x∈L⊥  ⟺  (x,ai)=0x\in L^\perp\iff(x,a_i) = 0 при i=1,…,ki = 1,\dots,k.

Доказательство

Развернуть

Пусть e1,…,eke_1,\dots,e_k — ортонормированный базис LL; дополним его до ортонормированного базиса e1,…,ene_1,\dots,e_n пространства VV. Так как xi=(x,ei)x_i = (x,e_i), вектор xx ортогонален LL тогда и только тогда, когда x1=⋯=xk=0x_1 = \dots = x_k = 0, то есть

L⊥=span⁡(ek+1,…,en).L^\perp = \operatorname{span}(e_{k+1},\dots,e_n).

Отсюда следуют п. 1 и 2. Аналогично (L⊥)⊥=span⁡(e1,…,ek)=L(L^\perp)^\perp = \operatorname{span}(e_1,\dots,e_k) = L.

  1. Если xx ортогонален всем aia_i, то он ортогонален и их линейным комбинациям. ■\blacksquare
Источники:Канатников, Крищенко, Линейная алгебра, гл. «Евклидовы пространства»Ильин, Позняк, Линейная алгебра, гл. «Евклидовы пространства»
теорема

Ортогональная проекция и расстояние до подпространства

Формулировка

  1. Ортогональная проекция вектора на подпространство существует и единственна.
  2. Если e1,…,eke_1,\dots,e_k — ортонормированный базис LL, то pr⁡Lx=∑i=1k(x,ei)ei\operatorname{pr}_Lx = \sum_{i=1}^k(x,e_i)e_i.
  3. Если a1,…,aka_1,\dots,a_k — произвольный базис LL, то pr⁡Lx=∑cjaj\operatorname{pr}_Lx = \sum c_ja_j, где коэффициенты — единственное решение системы
∑j=1k(ai,aj)cj=(x,ai),i=1,…,k.\sum_{j=1}^k(a_i,a_j)c_j = (x,a_i),\qquad i = 1,\dots,k.
  1. Проекция — ближайший к xx вектор подпространства: ∣x−u∣>∣x−pr⁡Lx∣|x - u| > |x - \operatorname{pr}_Lx| при u∈Lu\in L, u≠pr⁡Lxu\ne\operatorname{pr}_Lx. Поэтому ρ(x,L)=∣x−pr⁡Lx∣\rho(x,L) = |x - \operatorname{pr}_Lx|.

Доказательство

Развернуть
  1. Следует из V=L⊕L⊥V = L\oplus L^\perp.

  2. Вектор z=x−∑(x,ei)eiz = x - \sum(x,e_i)e_i удовлетворяет (z,ej)=(x,ej)−(x,ej)=0(z,e_j) = (x,e_j) - (x,e_j) = 0, то есть z∈L⊥z\in L^\perp.

  3. Условие x−∑cjaj⊥aix - \sum c_ja_j\perp a_i при всех ii — это указанная система. Её матрица — матрица Грама базиса, она невырождена, поэтому решение единственно.

  4. Пусть y=pr⁡Lxy = \operatorname{pr}_Lx, z=x−yz = x - y. Для u∈Lu\in L вектор y−uy - u лежит в LL и ортогонален zz, поэтому по теореме Пифагора

∣x−u∣2=∣z∣2+∣y−u∣2>∣z∣2при u≠y.|x - u|^2 = |z|^2 + |y - u|^2 > |z|^2\quad\text{при}\ u\ne y.

■\blacksquare

Источники:Канатников, Крищенко, Линейная алгебра, гл. «Евклидовы пространства»Ильин, Позняк, Линейная алгебра, гл. «Евклидовы пространства»
теорема

Изоморфизм конечномерных евклидовых пространств

Формулировка

Конечномерные евклидовы пространства изоморфны тогда и только тогда, когда их размерности равны. В частности, любое nn-мерное евклидово пространство изоморфно Rn\mathbb R^n со стандартным скалярным произведением.

Доказательство

Развернуть

Изоморфизм евклидовых пространств — в частности, изоморфизм линейных, поэтому размерности равны. Обратно, выберем ортонормированные базисы ee в VV и ff в WW и положим φ(∑xiei)=∑xifi\varphi\big(\sum x_ie_i\big) = \sum x_if_i. Это изоморфизм линейных пространств, а скалярное произведение в обоих пространствах равно ∑xiyi\sum x_iy_i в координатах ортонормированного базиса. ■\blacksquare

Источники:Канатников, Крищенко, Линейная алгебра, гл. «Евклидовы пространства»Ильин, Позняк, Линейная алгебра, гл. «Евклидовы пространства»
теорема

Метод наименьших квадратов. Нормальная система

Формулировка

  1. XX — псевдорешение системы AX=BAX = B тогда и только тогда, когда XX — решение нормальной системы
ATAX=ATB.A^\mathsf{T}AX = A^\mathsf{T}B.
  1. Нормальная система всегда совместна.
  2. Если rank⁡A=n\operatorname{rank}A = n (столбцы AA линейно независимы), псевдорешение единственно: X=(ATA)−1ATBX = (A^\mathsf{T}A)^{-1}A^\mathsf{T}B.

Доказательство

Развернуть

Когда XX пробегает Rn\mathbb R^n, вектор AXAX пробегает подпространство L=Im⁡AL = \operatorname{Im}A, порождённое столбцами A1,…,AnA_1,\dots,A_n.

  1. ∣AX−B∣|AX - B| минимально тогда и только тогда, когда AX=pr⁡LBAX = \operatorname{pr}_LB, то есть B−AX⊥LB - AX\perp L. Это равносильно AjT(B−AX)=0A_j^\mathsf{T}(B - AX) = 0 для всех столбцов, то есть AT(B−AX)=0A^\mathsf{T}(B - AX) = 0.

  2. Проекция pr⁡LB\operatorname{pr}_LB лежит в LL, поэтому равна AX0AX_0 для некоторого X0X_0, и X0X_0 — решение нормальной системы.

  3. Если ATAX=0A^\mathsf{T}AX = 0, то ∣AX∣2=XTATAX=0|AX|^2 = X^\mathsf{T}A^\mathsf{T}AX = 0, то есть AX=0AX = 0, и X=0X = 0 при независимых столбцах. Значит, ATAA^\mathsf{T}A невырождена. ■\blacksquare

Источники:Канатников, Крищенко, Линейная алгебра, гл. «Евклидовы пространства»Ильин, Позняк, Линейная алгебра, гл. «Евклидовы пространства»
теорема

Нормальное псевдорешение и псевдообратная матрица

Формулировка

  1. Нормальное псевдорешение существует и единственно; это единственное псевдорешение, лежащее в Im⁡AT\operatorname{Im}A^\mathsf{T} — линейной оболочке строк AA.
  2. Псевдообратная матрица A+A^+ существует и единственна.
  3. Если rank⁡A=n\operatorname{rank}A = n, то A+=(ATA)−1ATA^+ = (A^\mathsf{T}A)^{-1}A^\mathsf{T}; если rank⁡A=m\operatorname{rank}A = m, то A+=AT(AAT)−1A^+ = A^\mathsf{T}(AA^\mathsf{T})^{-1}; если AA квадратная невырожденная, то A+=A−1A^+ = A^{-1}.

Доказательство

Развернуть

Ядро. AX=0AX = 0 означает, что XX ортогонален всем строкам AA, поэтому Ker⁡A=(Im⁡AT)⊥\operatorname{Ker}A = (\operatorname{Im}A^\mathsf{T})^\perp и Rn=Im⁡AT⊕Ker⁡A\mathbb R^n = \operatorname{Im}A^\mathsf{T}\oplus\operatorname{Ker}A.

  1. Псевдорешения — это все XX с AX=pr⁡Im⁡ABAX = \operatorname{pr}_{\operatorname{Im}A}B, то есть множество X0+Ker⁡AX_0 + \operatorname{Ker}A. Разложим X0=X1+X2X_0 = X_1 + X_2, X1∈Im⁡ATX_1\in\operatorname{Im}A^\mathsf{T}, X2∈Ker⁡AX_2\in\operatorname{Ker}A. Тогда X1X_1 — псевдорешение, и для любого другого псевдорешения X1+KX_1 + K, 0≠K∈Ker⁡A0\ne K\in\operatorname{Ker}A, по теореме Пифагора ∣X1+K∣2=∣X1∣2+∣K∣2>∣X1∣2|X_1 + K|^2 = |X_1|^2 + |K|^2 > |X_1|^2. Два псевдорешения из Im⁡AT\operatorname{Im}A^\mathsf{T} различаются на вектор из Im⁡AT∩Ker⁡A={0}\operatorname{Im}A^\mathsf{T}\cap\operatorname{Ker}A = \{0\}.

  2. Отображение B↦X1B\mapsto X_1 — композиция ортогонального проектирования на Im⁡A\operatorname{Im}A и обратного к ограничению A∣Im⁡ATA|_{\operatorname{Im}A^\mathsf{T}}, которое взаимно однозначно отображает Im⁡AT\operatorname{Im}A^\mathsf{T} на Im⁡A\operatorname{Im}A (ядро нулевое, размерности равны рангу). Значит, оно линейно, и его матрица — A+A^+.

  3. При rank⁡A=n\operatorname{rank}A = n псевдорешение единственно (теорема о методе наименьших квадратов). При rank⁡A=m\operatorname{rank}A = m система совместна при любом BB, и вектор X=AT(AAT)−1BX = A^\mathsf{T}(AA^\mathsf{T})^{-1}B — её решение, лежащее в Im⁡AT\operatorname{Im}A^\mathsf{T}; значит, это нормальное псевдорешение. Матрица AATAA^\mathsf{T} невырождена по тем же соображениям, что и ATAA^\mathsf{T}A. ■\blacksquare

Источники:Канатников, Крищенко, Линейная алгебра, гл. «Евклидовы пространства»Ильин, Позняк, Линейная алгебра, гл. «Евклидовы пространства»