Доповнення Шура

Матеріал з testwiki
Перейти до навігації Перейти до пошуку

У лінійній алгебрі і теорії матриць доповнення Шура для блоку матриці (тобто, підматриці в більшій матриці) визначено так. Припустимо A, B, C, D є матриці відповідно p×p, p×q, q×p і q×q, і D оборотна. Нехай

M=[ABCD]

так що M — це матриця (p+q)×(p+q).

Тоді доповнення Шура для блоку D матриці M це матриця p×p

ABD1C.

Його назвали на честь Ісаї Шура, який використав його для доведення леми Шура, хоча його використовували і до того.[1]

Підґрунтя

Доповнення Щура виникає як результат застосування методу Гауса щодо блоків через множення на матрицю M на блокову нижньотрикутну матрицю

L=[Ip0D1CIq].

Тут Ip позначає одиничну матрицю p×p. Після множення на матрицю L доповнення Щура з'являється у горішньому p×p блоку. Матрицю добутку така

ML=[ABCD][Ip0D1CIq]=[ABD1CB0D]=[IpBD10Iq][ABD1C00D].

Це аналогічно до LDU-розкладу матриці. Тобто, ми щойно показали, що

[ABCD]=[IpBD10Iq][ABD1C00D][Ip0D1CIq],

отже, обернена до M можна представити за участю D−1 і оберненого доповнення Щура (якщо воно існує) як

[ABCD]1=[Ip0D1CIq][(ABD1C)100D1][IpBD10Iq]=[(ABD1C)1(ABD1C)1BD1D1C(ABD1C)1D1+D1C(ABD1C)1BD1].

Якщо M — симетрична додатноозначена матриця, то й так само буде доповнення Щура для D у M.

Якщо p і q дорівнюють 1 (тоюто A, B, C і D є скалярами), то ми отримуємо формулу для обернення матриці 2-на-2:

M1=1ADBC[DBCA]

за умови, що AD − BC не нуль.

Більше того, також чітко видно, що визначник M задається формулою

det(M)=det(D)det(ABD1C)

яка узагальнює формулу визначника у випадку матриць 2-на-2.

Умови на додатню визначеність і додатню напіввизначеність

Нехай X — це симетрична матриця задана так

X=[ABBTC].

Нехай X/A буде доповненням Щура для A в X, тобто

X/A=CBTA1B,

і X/C буде доповненням Щура для C в X, тобто

X/C=ABC1BT.

Тоді

  • X — додатно визначена тоді і тільки тоді коли A і X/A додатно визначені:
X0A0,X/A=CBTA1B0.
  • X — додатно визначена тоді і тільки тоді коли C і X/C додатно визначені:
X0C0,X/C=ABC1BT0.
  • Якщо A — додатно визначена, тоді X — додатно напіввизначена тоді і тільки тоді коли X/A є додатно напіввизначеною:
If A0, then X0X/A=CBTA1B0.
  • Якщо C є додатно визначеною, тоді X — додатно напіввизначеною тоді і тільки тоді коли X/C є додатно напіввизначеною:
If C0, then X0X/C=ABC1BT0.

Перше і третє твердження можна отримати [2][3] через розгляд мінімізатора величини

uTAu+2vTBTu+vTCv,

як функції від v (для фіксованого u).

Далі, оскільки

[ABBTC]0[CBTBA]0

і подібно для додатно напіввизначених матриць, друге (четверте) твердження негайно випливає з першого (відповідно третього) твердження.

Також існує необхідна і достатня умова на додатню напіввизначенність X в термінах узагальненого доповнення Щура.[1] А саме,

  • X0A0,CBTAgB0,(IAAg)B=0 і
  • X0C0,ABCgBT0,(ICCg)BT=0,

де Ag позначає узагальнену обернену матрицю для A.

Див. також

Примітки

Шаблон:Reflist

Джерела

  1. 1,0 1,1 Шаблон:Cite book
  2. Schur Complement Lemma Шаблон:Webarchive на berkeley.edu
  3. Boyd, S. and Vandenberghe, L. (2004), "Convex Optimization", Cambridge University Press (Appendix A.5.5)