ФОРМИРОВАНИЕ СПЕКТРА СИНГУЛЯРНЫХ ЧИСЕЛ КВАДРАТНОЙ МАТРИЦЫ ПРОСТОЙ СТРУКТУРЫ
Т.А. Акунов, Д.С. Бирюков, Н.А. Дударенко, М.В. Полякова, А.В. Ушаков
УДК 112:22.10:241.1
ФОРМИРОВАНИЕ СПЕКТРА СИНГУЛЯРНЫХ ЧИСЕЛ КВАДРАТНОЙ МАТРИЦЫ ПРОСТОЙ СТРУКТУРЫ
Т.А. Акунов, Д.С. Бирюков, Н.А. Дударенко, М.В. Полякова, А.В. Ушаков
Исследуется возможность формирования спектра сингулярных чисел квадратной матрицы простой структуры, для чего используется связь этого спектра со спектром ее собственных значений, которая определяется числом обусловленности модальной матрицы. Задача решается с использованием алгоритмических возможностей обобщенного модального управления. Приводится пример. Ключевые слова: модальное разложение, сингулярное разложение, сингулярные числа, собственные значения, собственные вектора, число обусловленности, обобщенное модальное управление.
Введение. Постановка задачи
Сингулярное разложение, SVD (singular value decomposition) [1–4] матриц прочно вошло в практику специалистов по системному SVD-анализу динамических процессов. Круг задач, решаемых с помощью сингулярного разложения [5] матриц, достаточно широк. Это задачи оценивания параметров, скаляризации векторных процессов в многомерных системах управления путем введения эллипсоидных оценок [6], приближения методом наименьших квадратов. Сингулярное разложение позволяет контролировать плавное изменение ранга матрицы с помощью ее числа обусловленности или функционалов вырождения [7, 8], оно стало конструктивным инструментом при решении задач регрессионного анализа [9], а также задач, основанных на грамианном подходе [10], и задач, связанных с контролем вырождения динамических систем типа многомерный вход – многомерный выход [7]. В теории оптических приборов SVD-анализ матрицы преобразования лучевого вектора [11] позволяет выделять факторы, вариация которых приводит к максимальной разъюстировке оптических систем. Тем не менее, задача формирования матриц с заданным спектром сингулярных чисел пока не решена, в то время как задача формирования матриц с заданным спектром собственных значений, получившая в современной теории динамических систем название «модальное управление», имеет хорошее алгоритмическое обеспечение и достаточно обширную библиографическую поддержку [6, 12, 13]. Решение задачи, поставленной заголовком настоящей работы, опирается на возможности инструментария обобщенного модального управления (ОМУ) [13] как расширенной версии модального управления.
Модальное разложение квадратной матрицы простой структуры
Рассмотрим матрицу N размерности dimN n n , которая характеризуется алгебраическим
спектром N собственных значений i ;i 1, n , вычисляемых в силу соотношения N i : detI N 0; i 1,n , и геометрическим спектром собственных векторов ξi ;i 1,n , удо-
влетворяющих условию ξi argNξi λiξi; i 1,n .
Если матрица N является матрицей простой структуры, то для элементов алгебраического спектра этой матрицы выполняется соотношение
λi λ j при i j;i, j 1, n .
Построим на элементах i ;i 1, n алгебраического спектра матрицы N диагональную матрицу Λ diag i;i 1,n . Нетрудно видеть, что в силу критериев подобия матриц [14] матрицы N и Λ явля-
ются подобными, так что существует неособая n n матрица M , позволяющая записать матричное
соотношение MΛ NM .
(1)
Следует заметить, что столбцы Mi ; i 1,n матрицы
M M1,M2,,Mi ,Mn
являются собственными векторами матрицы N , так что выполняются равенства
(2)
NMi λi Mi ; i 1, n ; Mi ξi .
(3)
Это легко доказывается выделением в левой и правой частях (1) i -ых столбцов Λ i и M i , для ко-
торых получим
MΛi NMi;i 1, n .
(4)
Таким образом, решение матричного уравнения (1) сведено к решению n матрично-векторных
уравнений (4). Для решения этого уравнения относительно столбца M i воспользуемся тем, что Λ i есть i -ый столбец диагональной матрицы Λ , т.е.
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 6 (76)
53
ФОРМИРОВАНИЕ СПЕКТРА СИНГУЛЯРНЫХ ЧИСЕЛ КВАДРАТНОЙ МАТРИЦЫ...
Λi 01i1,λi,01ni T .
(5)
Столбец Λ i вида (5) с учетом представления матрицы M в форме (2) приводит уравнение (4) к
виду (3).
Определение 1. Матрица M приведения подобия в форме (1), составленная из собственных век-
торов Mi ξi матрицы N , называется модальной матрицей (матрицы N ).
□
Определение 2. Модальным разложением матрицы N простой структуры называется представле-
ние ее в виде
N MΛΛ 1 .
□ (6)
Сингулярное разложение квадратной матрицы простой структуры
Определение 3. Сингулярным разложением [5] n n квадратной матрицы N , называется пред-
ставление ее в форме
N UΣ VT ,
(7)
где U – ортогональная n n матрица, столбцы Ui i 1,n которой образуют левый сингулярный базис,
так что выполняются соотношения
UTi U j
ij
1 при 0 при
i i
j j
UUT
UT U I
U
1;
(8)
V – ортогональная n n матрица, столбцы Vi i 1, n которой образуют правый сингулярный базис, так
что выполняются соотношения вида (8)
ViT
Vj
δij
1 при 0 при
i i
j j
VV T
VT V I
V
1;
(9)
Σ – диагональная матрица размерности n n с сингулярными числами i 0 на главной диагонали,
которая принимает вид
Σ diag i :i 1,n .
Оказываются справедливыми следующие положения.
1. Столбцы U i матрицы U являются собственными векторами матрицы NNT , так что выполняется векторно-матричное соотношение
NNT Ui μi Ui;i 1, n;
где μi – собственное значение матрицы NNT , вычисляемое в силу решения характеристического уравнения
det μ I NNT 0;
2. Столбцы V j матрицы V являются собственными векторами матрицы NT N так, что выполняется векторно-матричное соотношение
NT NVj μ j Vj ; j 1, n;
где μ j – собственное значение матрицы NT N , вычисляемое в силу решения характеристического уравнения
det I NT N 0;
3. Сингулярные числа i матрицы N вычисляются в силу соотношения
μ1i/2 ,i 1, n .
Остановимся на свойствах сингулярных чисел матрицы N .
Свойство 1. Сингулярные числа i матрицы N неотрицательны,
спектр N i ;i 1, n сингулярных чисел формируется так,
n
min i
i ;i 1, n
, 1 2
n .
i 0 ;i 1n, . При этом
что
α1
max i
αi ;i 1, n
,
Свойство 2. Характеристические полиномы матриц NT N и NNT в случае квадратной матрицы N совпадают, так что выполняются равенства
54 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2011, № 6 (76)
Т.А. Акунов, Д.С. Бирюков, Н.А. Дударенко, М.В. Полякова, А.В. Ушаков
det μI NNT det μI NT N .
Свойство 3. Матрицы N~ и N~~ , полученные умножением N на ортогональную матрицу R справа или ортогональную матрицу Q слева, обладают тем же спектром сингулярных чисел, что и матрица N , так что выполняются условия
N NR N ; N QN N .
Свойство 4. Если матрица N – диагональная, N diag i; i 1,n , то i i , следовательно
Σ diagi i .
Свойство 5. Если матрица N – симметричная, N NT , и характеризуется спектром собственных
значений N i;i 1,n , то сингулярные числа i этой матрицы удовлетворяют условию i i ,
а, следовательно, Σ diag i i .
Рассмотрим геометрическую интерпретацию сингулярного разложения. Матрица N действует на
элемент Vj j 1, n правого сингулярного базиса так, что он отображается в линейную оболочку, натя нутую на элемент U j j 1, n левого сингулярного базиса матрицы N , при этом выполняется векторно-
матричное равенство
NVj jU j ; j 1,n .
(10)
Линейное отображение
y Nx
(11)
отображает единичную сферу x 1 в эллипсоид, максимальная полуось которого, совпадает с вектором
y α1U1, а минимальная – с вектором y nUn . Действительно, разложим вектор x по элементам правого сингулярного базиса, записав его в форме
x
n
jVj;
j
arg
n
2j 1 .
j1 j1
(12)
x yVn p
1U1
V1
r nUn
xq
x 1
yt
y Nx
Рисунок. Геометрическая интерпретация линейного отображения y Nx сферы x 1
Тогда подстановка (12) в (11) с учетом (10) приводит к результату
y
Nx
n
j 1
j
NV
j
n
j 1
j
α jU j
y y
α1 αn
U1 Un
при при
x V1; x Vn.
Соотношения (13) хорошо иллюстрируются рисунком.
(13)
Связь спектров сингулярных чисел и собственных значений
Покажем, что степень близости спектра сингулярных чисел матрицы N к спектру модулей ее собственных значений в смысле свойств 4 и 1 (при этом носителем первого является диагональная матрица
Σ diag i :i 1,n , а второго – диагональная матрица Λ diag i;i 1,n ) определяется числом обу-
словленности СM модальной матрицы M в силу неравенства
Σ СM Λ ,
(14)
Действительно на основании (6) и (7) становится справедливым матричное соотношение
Σ UT M Λ M1 V , переход в котором к нормам с учетом свойств (8) и (9) приводит к (14).
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 6 (76)
55
ФОРМИРОВАНИЕ СПЕКТРА СИНГУЛЯРНЫХ ЧИСЕЛ КВАДРАТНОЙ МАТРИЦЫ...
Неравенство (14) положим в основу построения алгоритма конструирования матриц простой структуры с желаемым спектром сингулярных чисел путем формирование матрицы N с желаемым спектром собственных значений и спектром собственных векторов, доставляющих модальной матрице M
значение числа обусловленности СM, близкое к единице. Конструирование такой матрицы осуществ-
ляется с использованием возможностей метода обобщенного модального управления [13].
Основной результат. Алгоритм формирования квадратных матриц простой структуры с желаемым спектром сингулярных чисел
К настоящему моменту сформировалась процедура синтеза ОМУ [7, 12], с помощью которой матрице состояния проектируемой системы доставляется желаемая структура как собственных значений, так и собственных векторов.
Суть этой процедуры состоит в следующем. Если записать соотношение подобия (1) в виде однородного матричного уравнения MΛ NM 0 , то его с помощью аддитивно-мультипликативной деком-
позиции матрицы N в форме N A BK , где K HM1 , можно свести к неоднородному матричному
уравнению MΛ AM BH , именуемому матричным уравнением Сильвестра (УС). Решение УС отно-
сительно модальной матрицы M с желаемым числом обусловленности определяется выбором матрицы
H . Таким образом, на основе решения УС и использования соотношения (21) можно предложить сле-
дующий алгоритм синтеза матрицы простой структуры с заданным спектром сингулярных чисел.
Шаг 1. Задать желаемый спектр N αi : i 1,n сингулярных чисел матрицы N , сформиро-
вав матрицу Σ diag i;i 1,n ;
Шаг 2. Построить
Λ diag i : i i;i 1,n ;
на спектре N сингулярных чисел диагональную матрицу
Шаг 3. Построить матрицу M argCM 1, для чего взять произвольную n n матрицу P ,
сформировать ее SVD-разложение
P
U
p
Σ
p
V
T
p
и
положить
M U p или
M Vp , что обеспечит мат-
рице M ортонормированность столбцов, гарантирующую выполнения условия M argCM 1;
Шаг 4. Сформировать произвольную n n матрицу A такую, чтобы алгебраические спектры
A и Λ собственных значений матриц A и Λ не пересекались: A Λ ;
Шаг 5. Сформировать произвольную n r матрицу B такую, чтобы она с матрицей A образо-
вывала управляемую пару A, B; при этом, если rangB r n , то перейти к шагу 6 алгоритма, если
rangB r n , то перейти к шагу 8;
Шаг 6. Оценить принадлежность матрицы AM M образу матрицы B : AMMΛ ImB,
если включение выполняется, то перейти к шагу 7, иначе – к шагу 5;
Шаг 7. Решить матричное уравнение Сильвестра MΛ AM BH относительно матрицы H в
форме H BTB 1BT AM MΛ и перейти к шагу 8;
Шаг 8. Для случая rangB n решить матричное уравнение Сильвестра MΛ AM BH отно-
сительно матрицы H в форме Н argMΛ AM BH B1AM MΛ ;
Шаг 9. Сформировать матрицу N в форме N A BK , где K HM1 ;
Шаг 10. Осуществить сингулярное разложение матрицы N в форме (7) и проверить выполнение
условия
Σ Λ diag i i ;i 1,n .
■
Пример
Постановка задачи. Сконструировать ( 2 2 ) матрицу N со спектром сингулярных чисел
N 1 7;2 1.
Решение. Следуем алгоритму.
1.
Формируем матрицу
Σ
diagα1
7;α2
1
7 0
0 1;
2. Формируем диагональную матрицу Λ diag1 7;2 1;
56 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2011, № 6 (76)
Т.А. Акунов, Д.С. Бирюков, Н.А. Дударенко, М.В. Полякова, А.В. Ушаков
3.
Формируем
произвольную
невырожденную
( 22 )
–
матрицу
P
1 3
2 4
,
строим
ее
сингулярное
разложение
P
UpΣ pVTp
0,4045 0,9145
0,91455,465
0,4045
0
0 0,576 0,366 0,8174
0,8174 0,576 ,
формируем
мо-
дальную
матрицу
M
Up
0,4045 0,9145
0,9145
0,4045
;
4.
Формируем
2 2 матрицу
A
3 0
2 5
:
A
Λ
;
5.
Формируем
произвольную
2 2
матрицу
B
argrangB
n
2
1 1
0 2 ;
6. Выполняем Шаг 7. Решаем матричное уравнение Сильвестра относительно матрицы
Н
α
rgMΛ
AM
BH
B1AM
MΛ
5,874
2,55
2,849
2,638
;
7.
Формируем матрицу N
в
форме
N
A BK
A BHM1
1,982
2,22
2,22 6,02 ;
8. С помощью сингулярного разложения матрицы N строим ее спектр
N 1 7;2 1.
сингулярных
чисел
Заключение
На основе контроля числа обусловленности модальной матрицы с помощью алгоритмических возможностей обобщенного модального управления решена задача формирования матрицы с желаемым спектром сингулярных чисел. Задача решена для матриц простой структуры.
Литература
1. Замарашкин Н.Л., Тыртышников Е.Е. Распределение собственных и сингулярных чисел теплицевых
матриц при ослабленных требованиях к производящей функции // Мат. сборник. – 1997. – Т. 188. – С.
83–92.
2. Wall M.E., Rechtsteiner A., Rocha L.M. Singular value decomposition and principal component analysis // A Practical Approach to Microarray Data Analysis. – Kluwer. – 2003. – P. 91–109.
3. Yeung M.K., Tegner J., Collins J.J. Reverse engineering gene networks using singular value decomposition and robust regression // Proc Natl Acad Sci USA. – 2002. – V. 55. – № 9. – P. 6163–6168.
4. Romo T.D., Clarage J.B., Sorensen D.C., Phillips G.N. Jr. Automatic identification of discrete substates in proteins: singular value decomposition analysis of time-averaged crystallographic refinements // Proteins. – 1995. – V. 22. – P. 311–321.
5. Голуб Дж., Ван Лоун Ч. Матричные вычисления: Пер. с англ. – М.: Мир, 1999. – 548 c. 6. Акунов Т.А., Ушаков А.В. Анализ чувствительности эллипсоидных оценок многомерных процессов
управления // Изв. вузов. Приборостроение. – 1991. – Т. 34. – № 8. – С. 21–27. 7. Дударенко Н.А., Полякова М.В., Ушаков А.В. Вычислительные проблемы формирования функциона-
лов вырождения сложных технических систем с интервальными матричными компонентами // Проблемы управления. – 2011. – № 2. – C. 31–36. 8. Дударенко Н.А., Слита О.В., Ушаков А.В. Математические основы современной теории управления: аппарат метода пространства состояний: Учебное пособие. – СПб: СПбГУ ИТМО, 2005. – 324 с.
9. Jessup E.R., Sorensen D.C. A parallel algorithm for computing the singular-value decomposition of a matrix // Siam Journal on Matrix Analysis and Applications. – 1994. – V. 15. – № 2. – P. 530–548.
10. Бирюков Д.С., Ушаков А.В. Контроль затрат на управление при воспроизведении гармонических экзогенных воздействий: грамианный подход // Научно-технический вестник СПбГУ ИТМО. – 2011. – № 2. – C. 117–122.
11. Джерард А., Берч Дж. Введение в матричную оптику. – М.: Мир, 1978. – 341 с. 12. Григорьев В.В., Дроздов В.Н., Лаврентьев В.В., Ушаков А.В. Синтез дискретных регуляторов при
помощи ЭВМ – Л.: Машиностроение, Ленингр. отд-ние, 1983. – 357 с. 13. Ушаков А.В. Обобщенное модальное управление // Изв. вузов. Приборостроение. – 2000. – Т. 43. –
№ 3. – C. 8–15. 14. Гантмахер Ф.Р. Теория матриц. – М.: Наука, 1973. – 576 с.
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 6 (76)
57
ФОРМИРОВАНИЕ СПЕКТРА СИНГУЛЯРНЫХ ЧИСЕЛ КВАДРАТНОЙ МАТРИЦЫ...
Акунов Таалайбек Абакирович Бирюков Дмитрий Сергеевич Дударенко Наталия Александровна Полякова Майя Вячеславовна Ушаков Анатолий Владимирович
– Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, кандидат технических наук, докторант, takunov@mail.ru
– Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, аспирант,
quaint03@mail.ru – Санкт-Петербургский национальный исследовательский университет
информационных технологий, механики и оптики, кандидат технических наук, доцент, dudarenko@yandex.ru – Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, аспирант,
12noch@mail.ru – Санкт-Петербургский национальный исследовательский университет
информационных технологий, механики и оптики, доктор технических наук, профессор, ushakov-AVG@yandex.ru
58 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2011, № 6 (76)
УДК 112:22.10:241.1
ФОРМИРОВАНИЕ СПЕКТРА СИНГУЛЯРНЫХ ЧИСЕЛ КВАДРАТНОЙ МАТРИЦЫ ПРОСТОЙ СТРУКТУРЫ
Т.А. Акунов, Д.С. Бирюков, Н.А. Дударенко, М.В. Полякова, А.В. Ушаков
Исследуется возможность формирования спектра сингулярных чисел квадратной матрицы простой структуры, для чего используется связь этого спектра со спектром ее собственных значений, которая определяется числом обусловленности модальной матрицы. Задача решается с использованием алгоритмических возможностей обобщенного модального управления. Приводится пример. Ключевые слова: модальное разложение, сингулярное разложение, сингулярные числа, собственные значения, собственные вектора, число обусловленности, обобщенное модальное управление.
Введение. Постановка задачи
Сингулярное разложение, SVD (singular value decomposition) [1–4] матриц прочно вошло в практику специалистов по системному SVD-анализу динамических процессов. Круг задач, решаемых с помощью сингулярного разложения [5] матриц, достаточно широк. Это задачи оценивания параметров, скаляризации векторных процессов в многомерных системах управления путем введения эллипсоидных оценок [6], приближения методом наименьших квадратов. Сингулярное разложение позволяет контролировать плавное изменение ранга матрицы с помощью ее числа обусловленности или функционалов вырождения [7, 8], оно стало конструктивным инструментом при решении задач регрессионного анализа [9], а также задач, основанных на грамианном подходе [10], и задач, связанных с контролем вырождения динамических систем типа многомерный вход – многомерный выход [7]. В теории оптических приборов SVD-анализ матрицы преобразования лучевого вектора [11] позволяет выделять факторы, вариация которых приводит к максимальной разъюстировке оптических систем. Тем не менее, задача формирования матриц с заданным спектром сингулярных чисел пока не решена, в то время как задача формирования матриц с заданным спектром собственных значений, получившая в современной теории динамических систем название «модальное управление», имеет хорошее алгоритмическое обеспечение и достаточно обширную библиографическую поддержку [6, 12, 13]. Решение задачи, поставленной заголовком настоящей работы, опирается на возможности инструментария обобщенного модального управления (ОМУ) [13] как расширенной версии модального управления.
Модальное разложение квадратной матрицы простой структуры
Рассмотрим матрицу N размерности dimN n n , которая характеризуется алгебраическим
спектром N собственных значений i ;i 1, n , вычисляемых в силу соотношения N i : detI N 0; i 1,n , и геометрическим спектром собственных векторов ξi ;i 1,n , удо-
влетворяющих условию ξi argNξi λiξi; i 1,n .
Если матрица N является матрицей простой структуры, то для элементов алгебраического спектра этой матрицы выполняется соотношение
λi λ j при i j;i, j 1, n .
Построим на элементах i ;i 1, n алгебраического спектра матрицы N диагональную матрицу Λ diag i;i 1,n . Нетрудно видеть, что в силу критериев подобия матриц [14] матрицы N и Λ явля-
ются подобными, так что существует неособая n n матрица M , позволяющая записать матричное
соотношение MΛ NM .
(1)
Следует заметить, что столбцы Mi ; i 1,n матрицы
M M1,M2,,Mi ,Mn
являются собственными векторами матрицы N , так что выполняются равенства
(2)
NMi λi Mi ; i 1, n ; Mi ξi .
(3)
Это легко доказывается выделением в левой и правой частях (1) i -ых столбцов Λ i и M i , для ко-
торых получим
MΛi NMi;i 1, n .
(4)
Таким образом, решение матричного уравнения (1) сведено к решению n матрично-векторных
уравнений (4). Для решения этого уравнения относительно столбца M i воспользуемся тем, что Λ i есть i -ый столбец диагональной матрицы Λ , т.е.
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 6 (76)
53
ФОРМИРОВАНИЕ СПЕКТРА СИНГУЛЯРНЫХ ЧИСЕЛ КВАДРАТНОЙ МАТРИЦЫ...
Λi 01i1,λi,01ni T .
(5)
Столбец Λ i вида (5) с учетом представления матрицы M в форме (2) приводит уравнение (4) к
виду (3).
Определение 1. Матрица M приведения подобия в форме (1), составленная из собственных век-
торов Mi ξi матрицы N , называется модальной матрицей (матрицы N ).
□
Определение 2. Модальным разложением матрицы N простой структуры называется представле-
ние ее в виде
N MΛΛ 1 .
□ (6)
Сингулярное разложение квадратной матрицы простой структуры
Определение 3. Сингулярным разложением [5] n n квадратной матрицы N , называется пред-
ставление ее в форме
N UΣ VT ,
(7)
где U – ортогональная n n матрица, столбцы Ui i 1,n которой образуют левый сингулярный базис,
так что выполняются соотношения
UTi U j
ij
1 при 0 при
i i
j j
UUT
UT U I
U
1;
(8)
V – ортогональная n n матрица, столбцы Vi i 1, n которой образуют правый сингулярный базис, так
что выполняются соотношения вида (8)
ViT
Vj
δij
1 при 0 при
i i
j j
VV T
VT V I
V
1;
(9)
Σ – диагональная матрица размерности n n с сингулярными числами i 0 на главной диагонали,
которая принимает вид
Σ diag i :i 1,n .
Оказываются справедливыми следующие положения.
1. Столбцы U i матрицы U являются собственными векторами матрицы NNT , так что выполняется векторно-матричное соотношение
NNT Ui μi Ui;i 1, n;
где μi – собственное значение матрицы NNT , вычисляемое в силу решения характеристического уравнения
det μ I NNT 0;
2. Столбцы V j матрицы V являются собственными векторами матрицы NT N так, что выполняется векторно-матричное соотношение
NT NVj μ j Vj ; j 1, n;
где μ j – собственное значение матрицы NT N , вычисляемое в силу решения характеристического уравнения
det I NT N 0;
3. Сингулярные числа i матрицы N вычисляются в силу соотношения
μ1i/2 ,i 1, n .
Остановимся на свойствах сингулярных чисел матрицы N .
Свойство 1. Сингулярные числа i матрицы N неотрицательны,
спектр N i ;i 1, n сингулярных чисел формируется так,
n
min i
i ;i 1, n
, 1 2
n .
i 0 ;i 1n, . При этом
что
α1
max i
αi ;i 1, n
,
Свойство 2. Характеристические полиномы матриц NT N и NNT в случае квадратной матрицы N совпадают, так что выполняются равенства
54 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2011, № 6 (76)
Т.А. Акунов, Д.С. Бирюков, Н.А. Дударенко, М.В. Полякова, А.В. Ушаков
det μI NNT det μI NT N .
Свойство 3. Матрицы N~ и N~~ , полученные умножением N на ортогональную матрицу R справа или ортогональную матрицу Q слева, обладают тем же спектром сингулярных чисел, что и матрица N , так что выполняются условия
N NR N ; N QN N .
Свойство 4. Если матрица N – диагональная, N diag i; i 1,n , то i i , следовательно
Σ diagi i .
Свойство 5. Если матрица N – симметричная, N NT , и характеризуется спектром собственных
значений N i;i 1,n , то сингулярные числа i этой матрицы удовлетворяют условию i i ,
а, следовательно, Σ diag i i .
Рассмотрим геометрическую интерпретацию сингулярного разложения. Матрица N действует на
элемент Vj j 1, n правого сингулярного базиса так, что он отображается в линейную оболочку, натя нутую на элемент U j j 1, n левого сингулярного базиса матрицы N , при этом выполняется векторно-
матричное равенство
NVj jU j ; j 1,n .
(10)
Линейное отображение
y Nx
(11)
отображает единичную сферу x 1 в эллипсоид, максимальная полуось которого, совпадает с вектором
y α1U1, а минимальная – с вектором y nUn . Действительно, разложим вектор x по элементам правого сингулярного базиса, записав его в форме
x
n
jVj;
j
arg
n
2j 1 .
j1 j1
(12)
x yVn p
1U1
V1
r nUn
xq
x 1
yt
y Nx
Рисунок. Геометрическая интерпретация линейного отображения y Nx сферы x 1
Тогда подстановка (12) в (11) с учетом (10) приводит к результату
y
Nx
n
j 1
j
NV
j
n
j 1
j
α jU j
y y
α1 αn
U1 Un
при при
x V1; x Vn.
Соотношения (13) хорошо иллюстрируются рисунком.
(13)
Связь спектров сингулярных чисел и собственных значений
Покажем, что степень близости спектра сингулярных чисел матрицы N к спектру модулей ее собственных значений в смысле свойств 4 и 1 (при этом носителем первого является диагональная матрица
Σ diag i :i 1,n , а второго – диагональная матрица Λ diag i;i 1,n ) определяется числом обу-
словленности СM модальной матрицы M в силу неравенства
Σ СM Λ ,
(14)
Действительно на основании (6) и (7) становится справедливым матричное соотношение
Σ UT M Λ M1 V , переход в котором к нормам с учетом свойств (8) и (9) приводит к (14).
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 6 (76)
55
ФОРМИРОВАНИЕ СПЕКТРА СИНГУЛЯРНЫХ ЧИСЕЛ КВАДРАТНОЙ МАТРИЦЫ...
Неравенство (14) положим в основу построения алгоритма конструирования матриц простой структуры с желаемым спектром сингулярных чисел путем формирование матрицы N с желаемым спектром собственных значений и спектром собственных векторов, доставляющих модальной матрице M
значение числа обусловленности СM, близкое к единице. Конструирование такой матрицы осуществ-
ляется с использованием возможностей метода обобщенного модального управления [13].
Основной результат. Алгоритм формирования квадратных матриц простой структуры с желаемым спектром сингулярных чисел
К настоящему моменту сформировалась процедура синтеза ОМУ [7, 12], с помощью которой матрице состояния проектируемой системы доставляется желаемая структура как собственных значений, так и собственных векторов.
Суть этой процедуры состоит в следующем. Если записать соотношение подобия (1) в виде однородного матричного уравнения MΛ NM 0 , то его с помощью аддитивно-мультипликативной деком-
позиции матрицы N в форме N A BK , где K HM1 , можно свести к неоднородному матричному
уравнению MΛ AM BH , именуемому матричным уравнением Сильвестра (УС). Решение УС отно-
сительно модальной матрицы M с желаемым числом обусловленности определяется выбором матрицы
H . Таким образом, на основе решения УС и использования соотношения (21) можно предложить сле-
дующий алгоритм синтеза матрицы простой структуры с заданным спектром сингулярных чисел.
Шаг 1. Задать желаемый спектр N αi : i 1,n сингулярных чисел матрицы N , сформиро-
вав матрицу Σ diag i;i 1,n ;
Шаг 2. Построить
Λ diag i : i i;i 1,n ;
на спектре N сингулярных чисел диагональную матрицу
Шаг 3. Построить матрицу M argCM 1, для чего взять произвольную n n матрицу P ,
сформировать ее SVD-разложение
P
U
p
Σ
p
V
T
p
и
положить
M U p или
M Vp , что обеспечит мат-
рице M ортонормированность столбцов, гарантирующую выполнения условия M argCM 1;
Шаг 4. Сформировать произвольную n n матрицу A такую, чтобы алгебраические спектры
A и Λ собственных значений матриц A и Λ не пересекались: A Λ ;
Шаг 5. Сформировать произвольную n r матрицу B такую, чтобы она с матрицей A образо-
вывала управляемую пару A, B; при этом, если rangB r n , то перейти к шагу 6 алгоритма, если
rangB r n , то перейти к шагу 8;
Шаг 6. Оценить принадлежность матрицы AM M образу матрицы B : AMMΛ ImB,
если включение выполняется, то перейти к шагу 7, иначе – к шагу 5;
Шаг 7. Решить матричное уравнение Сильвестра MΛ AM BH относительно матрицы H в
форме H BTB 1BT AM MΛ и перейти к шагу 8;
Шаг 8. Для случая rangB n решить матричное уравнение Сильвестра MΛ AM BH отно-
сительно матрицы H в форме Н argMΛ AM BH B1AM MΛ ;
Шаг 9. Сформировать матрицу N в форме N A BK , где K HM1 ;
Шаг 10. Осуществить сингулярное разложение матрицы N в форме (7) и проверить выполнение
условия
Σ Λ diag i i ;i 1,n .
■
Пример
Постановка задачи. Сконструировать ( 2 2 ) матрицу N со спектром сингулярных чисел
N 1 7;2 1.
Решение. Следуем алгоритму.
1.
Формируем матрицу
Σ
diagα1
7;α2
1
7 0
0 1;
2. Формируем диагональную матрицу Λ diag1 7;2 1;
56 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2011, № 6 (76)
Т.А. Акунов, Д.С. Бирюков, Н.А. Дударенко, М.В. Полякова, А.В. Ушаков
3.
Формируем
произвольную
невырожденную
( 22 )
–
матрицу
P
1 3
2 4
,
строим
ее
сингулярное
разложение
P
UpΣ pVTp
0,4045 0,9145
0,91455,465
0,4045
0
0 0,576 0,366 0,8174
0,8174 0,576 ,
формируем
мо-
дальную
матрицу
M
Up
0,4045 0,9145
0,9145
0,4045
;
4.
Формируем
2 2 матрицу
A
3 0
2 5
:
A
Λ
;
5.
Формируем
произвольную
2 2
матрицу
B
argrangB
n
2
1 1
0 2 ;
6. Выполняем Шаг 7. Решаем матричное уравнение Сильвестра относительно матрицы
Н
α
rgMΛ
AM
BH
B1AM
MΛ
5,874
2,55
2,849
2,638
;
7.
Формируем матрицу N
в
форме
N
A BK
A BHM1
1,982
2,22
2,22 6,02 ;
8. С помощью сингулярного разложения матрицы N строим ее спектр
N 1 7;2 1.
сингулярных
чисел
Заключение
На основе контроля числа обусловленности модальной матрицы с помощью алгоритмических возможностей обобщенного модального управления решена задача формирования матрицы с желаемым спектром сингулярных чисел. Задача решена для матриц простой структуры.
Литература
1. Замарашкин Н.Л., Тыртышников Е.Е. Распределение собственных и сингулярных чисел теплицевых
матриц при ослабленных требованиях к производящей функции // Мат. сборник. – 1997. – Т. 188. – С.
83–92.
2. Wall M.E., Rechtsteiner A., Rocha L.M. Singular value decomposition and principal component analysis // A Practical Approach to Microarray Data Analysis. – Kluwer. – 2003. – P. 91–109.
3. Yeung M.K., Tegner J., Collins J.J. Reverse engineering gene networks using singular value decomposition and robust regression // Proc Natl Acad Sci USA. – 2002. – V. 55. – № 9. – P. 6163–6168.
4. Romo T.D., Clarage J.B., Sorensen D.C., Phillips G.N. Jr. Automatic identification of discrete substates in proteins: singular value decomposition analysis of time-averaged crystallographic refinements // Proteins. – 1995. – V. 22. – P. 311–321.
5. Голуб Дж., Ван Лоун Ч. Матричные вычисления: Пер. с англ. – М.: Мир, 1999. – 548 c. 6. Акунов Т.А., Ушаков А.В. Анализ чувствительности эллипсоидных оценок многомерных процессов
управления // Изв. вузов. Приборостроение. – 1991. – Т. 34. – № 8. – С. 21–27. 7. Дударенко Н.А., Полякова М.В., Ушаков А.В. Вычислительные проблемы формирования функциона-
лов вырождения сложных технических систем с интервальными матричными компонентами // Проблемы управления. – 2011. – № 2. – C. 31–36. 8. Дударенко Н.А., Слита О.В., Ушаков А.В. Математические основы современной теории управления: аппарат метода пространства состояний: Учебное пособие. – СПб: СПбГУ ИТМО, 2005. – 324 с.
9. Jessup E.R., Sorensen D.C. A parallel algorithm for computing the singular-value decomposition of a matrix // Siam Journal on Matrix Analysis and Applications. – 1994. – V. 15. – № 2. – P. 530–548.
10. Бирюков Д.С., Ушаков А.В. Контроль затрат на управление при воспроизведении гармонических экзогенных воздействий: грамианный подход // Научно-технический вестник СПбГУ ИТМО. – 2011. – № 2. – C. 117–122.
11. Джерард А., Берч Дж. Введение в матричную оптику. – М.: Мир, 1978. – 341 с. 12. Григорьев В.В., Дроздов В.Н., Лаврентьев В.В., Ушаков А.В. Синтез дискретных регуляторов при
помощи ЭВМ – Л.: Машиностроение, Ленингр. отд-ние, 1983. – 357 с. 13. Ушаков А.В. Обобщенное модальное управление // Изв. вузов. Приборостроение. – 2000. – Т. 43. –
№ 3. – C. 8–15. 14. Гантмахер Ф.Р. Теория матриц. – М.: Наука, 1973. – 576 с.
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2011, № 6 (76)
57
ФОРМИРОВАНИЕ СПЕКТРА СИНГУЛЯРНЫХ ЧИСЕЛ КВАДРАТНОЙ МАТРИЦЫ...
Акунов Таалайбек Абакирович Бирюков Дмитрий Сергеевич Дударенко Наталия Александровна Полякова Майя Вячеславовна Ушаков Анатолий Владимирович
– Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, кандидат технических наук, докторант, takunov@mail.ru
– Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, аспирант,
quaint03@mail.ru – Санкт-Петербургский национальный исследовательский университет
информационных технологий, механики и оптики, кандидат технических наук, доцент, dudarenko@yandex.ru – Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, аспирант,
12noch@mail.ru – Санкт-Петербургский национальный исследовательский университет
информационных технологий, механики и оптики, доктор технических наук, профессор, ushakov-AVG@yandex.ru
58 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2011, № 6 (76)