Достоверность решений линейных математических моделей с приближенно заданными исходными данными
На основании аппарата сингулярного разложения матриц и свойств соответствующих проекционных матриц получены новые оценки погрешности нормального псевдорешения в общем случае несовместных систем линейных алгебраических уравнений с возмущенными коэффициентами. Рассмотрены случаи уменьшения ранга возм...
Gespeichert in:
| Datum: | 2004 |
|---|---|
| Hauptverfasser: | , , |
| Format: | Artikel |
| Sprache: | Russisch |
| Veröffentlicht: |
Інститут проблем математичних машин і систем НАН України
2004
|
| Schlagworte: | |
| Online Zugang: | https://nasplib.isofts.kiev.ua/handle/123456789/83910 |
| Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
| Назва журналу: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| Zitieren: | Достоверность решений линейных математических моделей с приближенно заданными исходными данными / А.Н. Химич, С.А. Войцеховский, В.Н. Брусникин // Мат. машини і системи. — 2004. — № 3. — С. 41-47. — Бібліогр.: 6 назв. — рос. |
Institution
Digital Library of Periodicals of National Academy of Sciences of Ukraine| _version_ | 1860182609695342592 |
|---|---|
| author | Химич, А.Н. Войцеховский, С.А. Брусникин, В.Н. |
| author_facet | Химич, А.Н. Войцеховский, С.А. Брусникин, В.Н. |
| citation_txt | Достоверность решений линейных математических моделей с приближенно заданными исходными данными / А.Н. Химич, С.А. Войцеховский, В.Н. Брусникин // Мат. машини і системи. — 2004. — № 3. — С. 41-47. — Бібліогр.: 6 назв. — рос. |
| collection | DSpace DC |
| description | На основании аппарата сингулярного разложения матриц и свойств соответствующих проекционных матриц получены новые оценки погрешности нормального псевдорешения в общем случае несовместных систем линейных алгебраических уравнений с возмущенными коэффициентами. Рассмотрены случаи уменьшения ранга возмущенной матрицы, его сохранения и увеличения.
На основі методу сингулярного розвинення матриць та властивостей відповідних проекційних матриць отримані нові оцінки похибки нормального псевдорозв'язку в загальному випадку несумісних систем лінійних алгебраїчних рівнянь зі збуреними коефіцієнтами. Досліджено випадки зменшення рангу збуреної матриці, його збереження, та збільшення.
On the basis of a method of singular decomposition of the matrices and the propeties of corresponding projection matrices, there are obtained the new error estimates of the normal pseudosolution in the general case of the non-consistent systems of linear algebraic equations with perturbed coefficients. The cases of decreasing rank of the perturbed matrix, its saving and increasing are discussed.
|
| first_indexed | 2025-12-07T18:02:58Z |
| format | Article |
| fulltext |
ISSN 1028-9763. Математичні машини і системи, 2004, № 3
41
УДК 519.6
А.Н. ХИМИЧ, С.А. ВОЙЦЕХОВСКИЙ, В.Н. БРУСНИКИН
ДОСТОВЕРНОСТЬ РЕШЕНИЙ ЛИНЕЙНЫХ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ С
ПРИБЛИЖЕННО ЗАДАННЫМИ ИСХОДНЫМИ ДАННЫМИ
Abstract: On the basis of a method of singular decomposition of the matrices and the propeties of corresponding
projection matrices, there are obtained the new error estimates of the normal pseudosolution in the general case of
the non-consistent systems of linear algebraic equations with perturbed coefficients. The cases of decreasing rank of
the perturbed matrix, its saving and increasing are discussed.
Key words: matrice, error, linear algebraic equations.
Анотація: На основі методу сингулярного розвинення матриць та властивостей відповідних проекційних
матриць отримані нові оцінки похибки нормального псевдорозв'язку в загальному випадку несумісних
систем лінійних алгебраїчних рівнянь зі збуреними коефіцієнтами. Досліджено випадки зменшення рангу
збуреної матриці, його збереження, та збільшення.
Ключові слова: матриця, похибка, лінійні алгебраїчні рівняння.
Аннотация: На основании аппарата сингулярного разложения матриц и свойств соответствующих
проекционных матриц получены новые оценки погрешности нормального псевдорешения в общем случае
несовместных систем линейных алгебраических уравнений с возмущенными коэффициентами.
Рассмотрены случаи уменьшения ранга возмущенной матрицы, его сохранения и увеличения
Ключевые слова: матрица, погрешность, линейные алгебраические уравнения.
Математические модели с приближенно заданными исходными данными необходимо
рассматривать как задачи с априори неизвестными математическими свойствами: свойства
исходной математической модели могут отличаться от свойств математической модели с
приближенно заданными исходными данными. В этих условиях создание математического
аппарата для контроля за качеством получаемого решения является первостепенной задачей. В
работе на основании теории возмущений исследуется погрешность решения линейных
математических моделей, которые формулируются в виде систем линейных алгебраических
уравнений для матриц произвольного вида и ранга.
1. Постановка задачи
Рассмотрим в общем случае несовместную систему линейных алгебраических уравнений с точно
заданными исходными данными (математическую модель с точно заданными исходными данными)
bAx = (1)
с nm × матрицей и систему линейных алгебраических уравнений с приближенно заданными
исходными данными (математическая модель с приближенно заданными исходными данными)
bxА = , ,ААА −=∆ bbb −=∆ . (2)
Будем предполагать, что для погрешности элементов матрицы и правой части
выполняются следующие соотношения:
., bbAA bA ε∆ε∆ ≤≤ (3)
Получим оценки наследственной погрешности в общем случае нормального
псевдорешения систем линейных алгебраических уравнений. Наибольший практический интерес
представляет случай, когда ранг матрицы не сохраняется при возмущении ее коэффициентов. Для
полноты изложения здесь получены также известные результаты [1].
ISSN 1028-9763. Математичні машини і системи, 2004, № 3
42
2. Матрицы произвольного ранга
Сначала рассмотрим случай, когда ранг матрицы не изменяется при возмущении ее элементов.
Теорема 1. Предположим, что ║ А∆ ║║
+
А ║< 1, krankAАrank == . Тогда имеет
место оценка
−
++
−
≤
−
k
k
АbА
А
b
bb
h
h
h
x
xx
k
εεε
ε
2
1
, (4)
где kb – проекция правой части задачи (1) на главное левое сингулярное
подпространство матрицы А [2], т.е. kb ∈ Im A; ( )Ahh = – число обусловленности матрицы
A; символ , если не оговорено дополнительно, обозначает евклидову векторную норму и
согласованную с ней спектральную матричную норму;
+A – псевдообратная матрица Мура-
Пенроуза к матрице А.
Доказательство. Оценка погрешности следует из соотношения
)()( bbАbААxx −+−=− +++
, (5)
где для разности
++ − АА справедливо представление [1]
++ − АА =
++++++ −−−− АААPIQIAААААА
ТТТ ∆∆∆ )()( , (5а)
где ., ААPААQ ++ ==
Учитывая, что
АxbrrQIbQI −=−=− ,)()(
соотношения (5), (5а), получаем
=−xx ).()()( bbАbАААPIrQIАААbААА
ТТТ −+−−−− +++++++ ∆∆∆
Принимая во внимание, что ,1,1 ≤−≤− PIQI а также неравенство
+
+
+
−
≤
АА
А
А
∆1
, (6)
приходим к оценке:
)2(
1 xА
r
А
А
АА
xА
b
А
А
АА
АА
x
xx ∆∆∆
∆
+
+
+
++
−
≤
−
. (7)
Так как kk bbrbАx −== , , то из (7) следует (4), т.е. утверждение теоремы 1.
Далее рассмотрим случай, когда ранг матрицы неполного ранга может измениться при
возмущении ее элементов.
ISSN 1028-9763. Математичні машини і системи, 2004, № 3
43
Теорема 2. Предположим, что выполняется условие
+AA∆ < 1/2, rank A > rank А = k.
Тогда имеет место оценка
−
++
−
≤
−
k
k
АbА
А
k
b
bb
h
h
h
x
xx
k
εεε
ε
2
21
. (8)
Доказательство. Для получения оценки используем способ из [3], основанный на
сингулярном разложении матриц. Представим А в виде
TVUА Σ= . (9)
Как известно, какова бы ни была прямоугольная матрица Α размера nm × , всегда
существует разложение (9), где VU , – ортогональные матрицы; Σ – прямоугольная матрица
размера nm × с неотрицательными элементами на диагонали, которые являются сингулярными
числами матрицы Α .
Подпространство, совпадающее с линейной оболочкой первых к столбцов матрицы V
(U ), называется главным правым (левым) сингулярным подпространством матрицы A .
Отметим также, что главное правое (левое) сингулярное подпространство матрицы A
совпадает с подпространством образов матрицы
TA (матрицы A ).
Наряду с (9) рассмотрим разложение
,T
kk VUА Σ= (10)
где kΣ – прямоугольная матрица, первые к диагональных элементов которой отличны от нуля и
совпадают с соответствующими элементами матрицы Σ , а все остальные элементы равны нулю.
Используем тот факт из [2], что нормальное псевдорешение kx задачи наименьших
квадратов для системы
bxA kk = (11)
является ортогональной проекцией нормального псевдорешения задачи bxА = на главное
правое сингулярное подпространство матрицы A размерности к. Таким способом построенная
матрица (10) имеет ранг к , такой же, как и матрица невозмущенной задачи. Таким образом,
проблема оценки погрешности псевдорешения для матриц, ранг которых изменился, сведена к
случаю, когда ранги матриц одинаковы.
Для оценки погрешности в этом случае используем представления, аналогичные (5), (5а):
)()( bbАbААxx kkк
−+−=− +++
, (11а)
++ − kАА =
++++++ −−−− АААPIQIAААААА
ТТ
k
TТ
kkk ∆∆∆ )()( .
Представление (11а) получено из соотношения
ISSN 1028-9763. Математичні машини і системи, 2004, № 3
44
( )[ ] ( )[ ]( ) ( )[ ]
( ) ( )[ ] ( )[ ]( ) ( )[ ]
( ) ( ) ( ) ( )
( ) ( ) ( ) ( ) ( ) ( )
( ) ( ) ++++++
++++++
+++++
++++++
++++++
−−−+−
=−−−+−=−−+−
−−−+−−−−+
=−+−−++−+−
×−+=−+−−+=−
APIQIAAAAAAA
APIQIAAPQAQIQAPIQAPI
QIAPIQIAPQAPQIAPQAP
QIQAAPIPAAIAAAA
AAIAAQIQAAPIPAA
kkkk
kkkkk
kkkkkkkk
kkkk
kkkk
с использованием следующих проекционных матриц:
TT
k
T
kk
T
k
TT
kk
TT
k
T
kk
T
k
TT
kk AAAAQAAAAPAAAAQAAAAP ++++ ======== ,,, .
Учитывая этот факт и свойства проекционных матриц [1]
11 ≤−≤− kPI,QI ,
из (11а) получим
++
−
≤
− +
+
+
xА
r
А
А
АА
xА
b
А
А
АА
АА
x
xx
k
k
∆∆∆
∆
2
1
. (12)
Для оценки kk AAA −=∆ воспользуемся [4]: .2 ААk ∆∆ ≤ Кроме того, условие теоремы
2/1<∆ +AA приводит к неравенству 1<∆ +AAk , что является необходимым для
корректности выражения (12). Учитывая это, из (12) получаем оценку (8) для наследственной
погрешности нормального псевдорешения, что и требовадось доказать.
Теорема 3. Предположим, что >rankА lАrank = , 2/1<
∆
l
A
σ
. Тогда имеет место
оценка
l
ll
x
xx
σ∆Α
σσ
/21
/1
−
≤
−
−
++
l
l
l
АbA b
bb
l σ
σεεε 12 , (13)
где lx – проекция нормального псевдорешения задачи (1) на правое главное сингулярное
подпространство матрицы А размерности l ; lb – проекция правой части b на главное левое
сингулярное подпространство размерности l матрицы А; iσ – сингулярные числа матрицы А.
Доказательство. Наряду с задачей (1) рассмотрим задачу определения нормального
псевдорешения системы
bxA ll = (14)
с матрицей Τ∑= VUA ll ранга l .
Записывая (11а) для задач (2), (14), ранги матриц которых совпадают,
++ − ААl =
++++++ −−−− l
Т
l
Т
l
TТ
l ААА)PI()QI(AААААА ∆∆∆ ,
ISSN 1028-9763. Математичні машини і системи, 2004, № 3
45
)-()( bbAbAAxx ll
+++ +−=−
и, учитывая свойства проекционных матриц, придем к оценке:
l
l
x
x−x ≤
||||||||1
||||||||
+
+
− ll
ll
AA
АА
∆
(2
||||
||||
lA
А∆
+
||||||||
||||
ll xA
b∆
+
+
ll AA
||||||||||||
||||||||
lll xAA
rА∆
), (15)
откуда следует (13), т.е. утверждение теоремы 3.
3. Матрицы полного ранга
Для матриц полного ранга существенно то обстоятельство, что их ранг не меняется от возмущения
элементов, если выполнено условие ║ А∆ ║║
+
А ║< 1. Кроме того, в дальнейшем будем
использовать следующее свойство матриц полного ранга [2]:
TT AAAA 1)( −+ = для mn ≤ и
1)( −+ = TT AAAA для mn ≥ . (16)
Аналогичные соотношения выполняются и для матрицы A .
Теорема 4. Предположим, что 1<∆ +AA , knm => . Тогда имеет место оценка
−
++
−
≤
−
k
k
АbА
А
b
bb
h
h
h
x
xx
k
εεε
ε1
. (17)
Доказательство. Для доказательства теоремы 4, как и ранее, будем использовать
соотношение (5). В силу (16) IAAP == +
, так что из (5а) имеем равенство
++ − AA = )( QIАААААА
Т −− ++++ ∆∆ ,
используя которое приходим к (17).
Теорема 5 Предположим, что 1<∆ +AA , kmn => . Тогда
x
xx −
≤ ( )bA
А
h
h εε
ε
+
−
2
1
. (18)
Доказательство. Поскольку в этом случае IAAQ == +
, то выражение для
++ − AA в
силу (5а) приобретает вид
++ − АА =
++++ −− АААPIААА
ТТ∆∆ )( .
Дальнейшие выкладки аналогичны предыдущим. В результате приходим к оценке (18).
Теорема 6. Предположим, что 1<∆ +AA , knm == . Тогда
≤
−
x
xx
)(
1 bА
А
h
h εε
ε
+
−
. (19)
Доказательство. Оценка следует из (5) и равенства
ISSN 1028-9763. Математичні машини і системи, 2004, № 3
46
++ − АА =
++
ААА ∆ , (20)
которое следует из (5а) с учетом того, что для невырожденных матриц всегда имеют место
равенства IPIQ == , .
Замечание 1. Оценки (4), (8), (.13), (17), (.18), (19) имеют мажорантный характер на всем
классе матриц для используемой модели погрешности (3). Тем не менее, по крайней мере, для
невырожденных матриц в [5] показана достижимость оценки (19) для первой матричной нормы и
тем самым доказана неулучшаемость оценки в классе невырожденных матриц.
Замечание 2. Для других ограничений на погрешность исходных данных отличных от (3)
могут быть получены оценки, учитывающие их вид. Например, для AА Aε∆ = , bb bε∆ = из (5),
(20) получим оценку
≤
−
x
xx
А
bА
ε
εε
−
+
1
)(
, (21)
не содержащую множитель ( )Ah .
Замечание 3. Связь между числом обусловленности задачи с точными исходными
данными h(A) и числом обусловленности матрицы системы с приближенно заданными исходными
данными )(Ah устанавливает оценка
hАh
Аh
h A
A
A
A
ε
ε
ε
ε
−
+≤≤
+
−
1
1
)(
)(
1
1
, (22)
которая получена для спектральной матричной нормы на основании теории возмущения
сингулярных чисел [1].
Замечание 4. Наряду с известным представлением (5а) имеет место соотношение
++ − АA = -
++++++ −+−− АААPIQIAААААА ∆∆∆ )()( ,
в справедливости которого можно убедиться непосредственной проверкой. Здесь
., ААPААQ ++ ==
Используя выражение для погрешности:
)()( bbАbААxx −+−=− +++
, (23)
и учитывая, что ,,)()( xАbrrQIbQI −=−=− из (23) получим
=− xx )()()( bbАbАААPIrQIАААbААА −+−+−−− +++++++ ∆∆∆ . (24)
Отсюда следует справедливость следующей теоремы.
Теорема 7. Предположим, что выполняется условие 1<+AA∆ и
krankAАrank == . Тогда для погрешности нормального псевдорешения имеет место оценка
ISSN 1028-9763. Математичні машини і системи, 2004, № 3
47
−
++
−
≤
−
k
k
AbА
А
b
bb
Ah
Ah
Ah
x
xx
k
εεε
ε
)(2
)(1
)(
, (25)
где Almbk ∈ – проекция правой части b на образ матрицы A размерности k ,
.; kbA
bbAA
k
εε ≤∆≤∆
Таким образом, оценки наследственной погрешности, правая часть которых определяется
приближенными данными, могут быть получены и без неравенств (22). Оценки, аналогичные (25),
могут быть получены для всех рассмотренных ранее случаев.
Замечание 5. При выполнении условий теоремы 7, используя неравенство
−
+
−
≤
−
x
xx
x
xx
x
xx
1
и неравенство (25) придем к оценке
−
++
−
=
−
≤
−
k
k
AbА
А
b
bb
Ah
Ah
Ah
x
xx
k
εεε
ε
β
β
β
)(2
)(1
)(
,
1
. (26)
Оценки, аналогичные (26), очевидно могут быть получены для всех ранее рассмотренных
случаев.
4. Заключение
В работе на основании аппарата сингулярного разложения матриц и свойств соответствующих
проекционных матриц (11а) получены новые, отличные от [6], оценки наследственной погрешности
в общем случае нормального псевдорешения систем линейных алгебраических уравнений.
Рассмотрены случаи уменьшения ранга возмущенной матрицы, его сохранения и увеличения.
СПИСОК ЛИТЕРАТУРЫ
1. Лоусон Ч., Хенсон Р. Численное решение задач метода наименьших квадратов. – М.: Наука, 1986. – 230 с.
2. Воеводин В.В., Кузнецов Ю.А. Матрицы и вычисления. – М.: Наука, 1984. – 318 с.
3. Химич А.Н. Оценки возмущений для решения задачи наименьших квадратов // Кибернетика и системный
анализ. – 1996. – № 3. – С. 142 – 145.
4. Химич А.Н. Оценки возмущений для решения вырожденных СЛАУ // Математическое и программное
обеспечение прикладных систем новых классов и поколений. – Киев: Ин-т кибернетики им. В.М. Глушкова АН
УССР, 1996. – С. 63 – 67.
5. Молчанов И.Н, Николенко Л.Д., Кириченко М.П. Об одном пакете программ для решения систем линейных
алгебраических уравнений // Кибернетика. – 1972. – № 1. – С. 127 – 134.
6. Химич А.Н. Оценки полной погрешности решения систем линейных алгебраических уравнений для матриц
произвольного ранга // Компьютерная математика: Сборник научных трудов. – Киев: Институт кибернетики им.
В.М. Глушкова НАН Украины, 2002. – № 2. – С. 41 – 49.
|
| id | nasplib_isofts_kiev_ua-123456789-83910 |
| institution | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| issn | 1028-9763 |
| language | Russian |
| last_indexed | 2025-12-07T18:02:58Z |
| publishDate | 2004 |
| publisher | Інститут проблем математичних машин і систем НАН України |
| record_format | dspace |
| spelling | Химич, А.Н. Войцеховский, С.А. Брусникин, В.Н. 2015-06-28T09:44:05Z 2015-06-28T09:44:05Z 2004 Достоверность решений линейных математических моделей с приближенно заданными исходными данными / А.Н. Химич, С.А. Войцеховский, В.Н. Брусникин // Мат. машини і системи. — 2004. — № 3. — С. 41-47. — Бібліогр.: 6 назв. — рос. 1028-9763 https://nasplib.isofts.kiev.ua/handle/123456789/83910 519.6 На основании аппарата сингулярного разложения матриц и свойств соответствующих проекционных матриц получены новые оценки погрешности нормального псевдорешения в общем случае несовместных систем линейных алгебраических уравнений с возмущенными коэффициентами. Рассмотрены случаи уменьшения ранга возмущенной матрицы, его сохранения и увеличения. На основі методу сингулярного розвинення матриць та властивостей відповідних проекційних матриць отримані нові оцінки похибки нормального псевдорозв'язку в загальному випадку несумісних систем лінійних алгебраїчних рівнянь зі збуреними коефіцієнтами. Досліджено випадки зменшення рангу збуреної матриці, його збереження, та збільшення. On the basis of a method of singular decomposition of the matrices and the propeties of corresponding projection matrices, there are obtained the new error estimates of the normal pseudosolution in the general case of the non-consistent systems of linear algebraic equations with perturbed coefficients. The cases of decreasing rank of the perturbed matrix, its saving and increasing are discussed. ru Інститут проблем математичних машин і систем НАН України Програмно-технічні комплекси Достоверность решений линейных математических моделей с приближенно заданными исходными данными Вірогідність розв'язків лінійних математичних моделей з наближено заданими вихідними даними The reliability of solutions of the linear mahematical models with approximately given initial data Article published earlier |
| spellingShingle | Достоверность решений линейных математических моделей с приближенно заданными исходными данными Химич, А.Н. Войцеховский, С.А. Брусникин, В.Н. Програмно-технічні комплекси |
| title | Достоверность решений линейных математических моделей с приближенно заданными исходными данными |
| title_alt | Вірогідність розв'язків лінійних математичних моделей з наближено заданими вихідними даними The reliability of solutions of the linear mahematical models with approximately given initial data |
| title_full | Достоверность решений линейных математических моделей с приближенно заданными исходными данными |
| title_fullStr | Достоверность решений линейных математических моделей с приближенно заданными исходными данными |
| title_full_unstemmed | Достоверность решений линейных математических моделей с приближенно заданными исходными данными |
| title_short | Достоверность решений линейных математических моделей с приближенно заданными исходными данными |
| title_sort | достоверность решений линейных математических моделей с приближенно заданными исходными данными |
| topic | Програмно-технічні комплекси |
| topic_facet | Програмно-технічні комплекси |
| url | https://nasplib.isofts.kiev.ua/handle/123456789/83910 |
| work_keys_str_mv | AT himičan dostovernostʹrešeniilineinyhmatematičeskihmodeleispribližennozadannymiishodnymidannymi AT voicehovskiisa dostovernostʹrešeniilineinyhmatematičeskihmodeleispribližennozadannymiishodnymidannymi AT brusnikinvn dostovernostʹrešeniilineinyhmatematičeskihmodeleispribližennozadannymiishodnymidannymi AT himičan vírogídnístʹrozvâzkívlíníinihmatematičnihmodeleiznabliženozadanimivihídnimidanimi AT voicehovskiisa vírogídnístʹrozvâzkívlíníinihmatematičnihmodeleiznabliženozadanimivihídnimidanimi AT brusnikinvn vírogídnístʹrozvâzkívlíníinihmatematičnihmodeleiznabliženozadanimivihídnimidanimi AT himičan thereliabilityofsolutionsofthelinearmahematicalmodelswithapproximatelygiveninitialdata AT voicehovskiisa thereliabilityofsolutionsofthelinearmahematicalmodelswithapproximatelygiveninitialdata AT brusnikinvn thereliabilityofsolutionsofthelinearmahematicalmodelswithapproximatelygiveninitialdata |