Метод збурення псевдообернених
матриць [1] на основі принципу розщеплення матриць нижче поширюється на
проекційні матриці з метою подальшого використання при розв'язанні задач ідентифікації, нелінійного
регресійного аналізу, апроксимації
функцій і прогнозу.
Відповідно до постановки задачі про аналітичне представлення збурень
псевдообернених матриць [7, 8], будемо розглядати для деякої довільної матриці її псевдообернену
матрицю
, збурену матрицю
,
,
,
збурену псевдообернену матрицю
,
збурену проекційну матрицю
,
а також наступну проекційну матрицю
.
Функції ,
,
мають різний вигляд в
залежності від того, можна або неможливо представити вектори
й
у формі лінійних
комбінацій векторів-стовпчиків або, відповідно, вектор-рядків матриці
.
Розглянемо чотири можливих випадки
залежності векторів і
від елементів матриці
.
Випадок 1.
Вектори і
лінійно незалежні з
векторами-стовпцями і векторами-рядками матриці
відповідно, тобто
. (2.1)
Тоді залежність визначається
наступною теоремою.
Теорема 1. Якщо для матриці ,
,
, виконуються умови (2.1), то
.
Використовуючи співвідношення (2.1) для функцій ,
, їхній вид визначається наслідками з теореми 1.
Наслідок 1.
Якщо виконуються умови теореми, то
. (2.2)
Справедливість цього твердження прямо випливає з теореми 1. Дійсно
.
Наслідок 2. Якщо виконуються умови теореми 1, то
. (2.3)
Наслідок 3. Якщо виконуються умови теореми 1 і
,
,
тобто, вектор є ортогональним до
усіх векторів-стовпців матриці
, а вектор
– до всіх
вектор-рядкам матриці
, то
. (2.4)
Справедливість твердження наслідку 3 випливає з формули (2.3) і
співвідношень
,
,
.
Випадок 2.
Вектор є лінійно залежним
від вектор-стовпців матриці
, а вектор
– лінійно незалежним
від вектор-рядків матриці
, тобто
,
. (2.5)
Тут має місце наступна теорема [8].
Теорема 2. Якщо для матриці ,
,
виконуються умови
(2.5), то
(2.6)
де
.
Наслідок 4. Якщо виконуються умови теореми 2 і вектор є ортогональним до
вектор-рядків матриці
, тобто
, то
, (2.7)
де
. (2.8)
Наслідок 5. Якщо мають місце умови наслідку 4, то
, (2.9)
де визначається по
формулі (2.8).
Доведення наслідку (5) випливає з наступних співвідношень
.
Наслідок 6. Якщо мають місце умови наслідку 4, то
. (2.10)
За умови вирази
й
у цьому випадку
будуть отримані після аналізу випадку лінійної залежності
і
від відповідно
вектор-стовпців і векторів-рядків матриці
.
Випадок 3.
Вектори і
лінійно залежні
відповідно від вектор-стовпців і векторів-рядків матриці
, тобто
,
(2.11)
і при цьому
. (2.12)
У цьому випадку збурення псевдооберненої матриці визначаються відповідно
наступної теореми.
Теорема 3. Якщо для матриці ,
,
виконуються умови
(2.11), (2.12), то мають місце співвідношення
, (2.13)
. (2.14)
Наслідок 7. Якщо виконані умови теореми 3, то
. (2.15)
Справедливість твердження наслідку 7 перевіряється простою підстановкою
формули (2.15) у вираз для матриці
,
де використані властивості
,
( відповідно до
(2.11) ),
( відповідно до
(2.12) ).
Наслідок 8. Якщо виконуються умови теореми 3, то
, (2.16)
де визначається по
формулі (2.14).
При доведенні теореми 3 використовується
наступна лема.
Лема. Квадратна матриця при наявності
умови
має наступну
псевдообернену матрицю
. (2.17)
Ця лема буде використана нижче для обчислення збурень у псевдооберненій матриці, якщо під збуреннями початкової матриці розуміти видалення з неї деякого її стовпця або рядка.
Випадок 4. Нехай мають
місце умови (2.11)
,
але при цьому не виконується рівність (2.12), тобто
. (2.18)
Тоді, як це випливає з доведення теореми 3, виконується умова
. (2.19)
Теорема 4. Якщо для матриці ,
,
виконуються умови
(2.11), (2.18), то мають місце співвідношення
, (2.20)
Доведення цього твердження здійснюється перевіркою умов, яким повинна
задовольняти псевдообернена матриця.
Наслідок 9. Якщо виконані умови теореми 4, то
. (2.21)
Наслідок 10. Якщо виконуються умови теореми 4, то
.(2.22)
Довести останню формулу можна наступним чином
.
Наслідок 11. Якщо мають місце умови теореми 2 і ,
, то
, (2.23)
, (2.24)
, (2.25)
,
. (2.26)
Доведення
наслідку 11 випливає з розбивки вектора на два ортогональні
складові
і застосування
послідовно до матриці
теореми 4, наслідку
9, а потім до матриці
наслідку 4 і 5.
Розглянемо тепер
збурення -матриці
у формі або
поповнення її новим рядком до встановлення її розмірності
або вилучення з неї
однієї з її рядків із зміною розмірності матриці до
. В якості збурень без обмеження спільності будемо розглядати
поповнення матриці або вилучення останнього рядка. При поповненні матриці новим
рядком для визначення збуреної псевдооберненої матриці використовуються добре
відомі формули Гревіля [5]
Нехай
для матриці має місце
співвідношення
, тоді
, (2.27)
, (2.28)
(2.29)
якщо
, то
, (2.30)
, (2.31)
. (2.32)
При
вилученні з матриці останнього рядка
псевдообернена і проекційні матриці набувають
наступні зміни
a) при
, (2.33)
де
, (2.34)
, (2.35)
, (2.36)
,
(2.37)
b)
при
, (2.38)
, (2.39)
, (2.40)
, (2.41)
Слід
зазначити, що умова (2.33) визначає падіння рангу в матриці при вилученні
вектор-рядки , тобто
, (2.42)
а
умова (2.38) – відсутність зниження рангу
, (2.43)
Наведемо
доведення формули (2.35). Якщо
відома
матриця, де – її останній стовпчик. Тоді, відповідно до (2.27), при
одержимо
,
,
тобто,
. Так як
, то
і, відповідно до
(2.17)
,
,
.
Тут
використані співвідношення
.
Подібним
чином доводиться і формула (2.37).
Приведемо
доведення формули (2.39). При відповідно до (2.30)
і
з рівняння
одержимо
.
Формули
(40) і (41) доводяться простим використанням формули (2.39).