2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки


Правила форума


В этом разделе нельзя создавать новые темы.

Если Вы хотите задать новый вопрос, то не дописывайте его в существующую тему, а создайте новую в корневом разделе "Помогите решить/разобраться (М)".

Если Вы зададите новый вопрос в существующей теме, то в случае нарушения оформления или других правил форума Ваше сообщение и все ответы на него могут быть удалены без предупреждения.

Не ищите на этом форуме халяву, правила запрещают участникам публиковать готовые решения стандартных учебных задач. Автор вопроса обязан привести свои попытки решения и указать конкретные затруднения.

Обязательно просмотрите тему Правила данного раздела, иначе Ваша тема может быть удалена или перемещена в Карантин, а Вы так и не узнаете, почему.



Начать новую тему Ответить на тему
 
 Поиск количества решений матричного ур-я Риккати
Сообщение13.02.2009, 15:48 
Аватара пользователя


13/02/09
6
Здача: найти все $(n * n)$ матрицы $P$, удовлетворяющие уравнению: $P^2 + P^TB + B^TP + C = 0$, где $B$ - произвольная $(n * n)$ матрица, $C$ - симметрическая $(n * n)$ матрица.
Что сделано: рассмотрел матрицу $A=$$ 
\left( \begin{array}{cc} -B & -E \\ 
C & B^T \end{array} \right)$ размера $(2n * 2n)$.
Теорема: пусть $P_0$ является решением исходного ур-я, тогда подпространство вида $L_P_0=[z \in \mathbb{R}^{2n}: z=(x,Px), x \in \mathbb{R}^n]$ является инвариантным подпространством оператора $A$.
Доказательство: $ A 
\left( \begin{array}{cc} x \\ 
Px \end{array} \right) = \left( \begin{array}{cc} -Bx-Px \\ 
Cx+B^TPx \end{array} \right) = \left( \begin{array}{cc} -Bx-Px \\ 
P(-Bx-Px) \end{array} \right) = \left( \begin{array}{cc} y \\ 
Py \end{array} \right)$, где $y=-Bx-Px$ в силу равенства $C+B^TP=-PB-P^2$.
Из этой теоремы делаю вывод, что ЕСЛИ ВСЕ СОБСТВЕННЫЕ ЗНАЧЕНИЯ МАТРИЦЫ $A$ РАЗЛИЧНЫ(Т.Е. НЕТ КРАТНЫХ КОРНЕЙ ХАР-ОГО УР-Я), то у матрицы $A$ будет конечное число различных инвариантных подпространств размерности $n$ и все их можно найти, приводя матрицу к базису из собственных векторов и пар векторов, линейная оболочка которых инвариантна оператору $A$. Среди этих подпространств нужно найти те, которые однозначно проектируются на линейную оболочку первых $n$ базисных векторов в $\mathbb{R}^{2n}$.

А теперь, собственно, мой вопрос: а что же делать со случаем, когда кратные корни все-таки есть??? Как тогда искать все решения?? Как оценить их количество? Буду очень благодарен за любые советы, комментарии и мнения!!

 Профиль  
                  
 
 
Сообщение14.02.2009, 10:15 
Заслуженный участник


26/12/08
678
В вашем доказательстве инвариантности неточность: $(C+B^TP)x=-(P^2+P^TB)x\ne-P(P+B)x$. Возможно, требовалось решить уравнение $P^TP+P^TB+B^TP+C=0$? Последнее уравнение заменой $S=P+B$, $F=B^TB-C$ сводится к уравнению $S^TS=F$. Если $F<0$, то решений нет, иначе $S=Q\sqrt{F}$, где $Q$ - произвольная ортогональная матрица.

 Профиль  
                  
 
 
Сообщение14.02.2009, 13:39 


29/09/06
4552
Могу лишь сообщить, что имеется книга А.И.Егоров. Уравнение Риккати. Одна глава, по-моему, посвящена матричным уравнениям. Если она Вам неизвестна и желательна, поищите в Интернете (я находил когда-то).

 Профиль  
                  
 
 
Сообщение14.02.2009, 17:52 
Заслуженный участник


11/05/08
32166
Полосин писал(а):
Если $F<0$, то решений нет

Опечатка: имелось в виду "если $F\not\geqslant0$".
Полосин писал(а):
, иначе $S=Q\sqrt{F}$, где $Q$ - произвольная ортогональная матрица.

Если матрица $F$ строго положительна, то безусловно. А вот если она вырождена, то не уверен; мне кажется, что там произвол больше. Хотя могу и ошибаться.

 Профиль  
                  
 
 
Сообщение25.02.2009, 21:04 
Аватара пользователя


13/02/09
6
Полосин писал(а):
В вашем доказательстве инвариантности неточность: $(C+B^TP)x=-(P^2+P^TB)x\ne-P(P+B)x$. Возможно, требовалось решить уравнение $P^TP+P^TB+B^TP+C=0$? Последнее уравнение заменой $S=P+B$, $F=B^TB-C$ сводится к уравнению $S^TS=F$. Если $F<0$, то решений нет, иначе $S=Q\sqrt{F}$, где $Q$ - произвольная ортогональная матрица.


Нет, здесь я имел в виду именно уравнение $P^2+P^TB+B^TP+C=0$, но забыл уточнить, что таким образом можно искать только симметрические его решения. Что же касается уравнения $P^TP+P^TB+B^TP+C=0$, то его я тоже рассматривал, и, действительно, его общее решение имеет вид $P=-B+U\sqrt{B^TB-C}$, где $U$ - произвольная ортогональная матрица, при условии, что матрица $B^TB-C$ - симметрическая, неотрицательно определенная.

Добавлено спустя 4 минуты 37 секунд:

Алексей К. писал(а):
Могу лишь сообщить, что имеется книга А.И.Егоров. Уравнение Риккати. Одна глава, по-моему, посвящена матричным уравнениям. Если она Вам неизвестна и желательна, поищите в Интернете (я находил когда-то).

Благодарю за ссылку! Книжку скачал, изучаю... Но пока ничего полезного для моей конкретной задачи выудить не получается, там больше расписана дифференциальная составляющая уравнений Риккати.
Но все равно большое спасибо, так как за последнее время я просмотрел много литературы на данную тематику, но о существовании этой книжки до сих пор не знал...

 Профиль  
                  
 
 
Сообщение25.02.2009, 21:09 
Заслуженный участник


26/12/08
678
Если вы ищете только симметрические решения, то зачем тогда писать $P^T$?

 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ Сообщений: 6 ] 

Модераторы: Модераторы Математики, Супермодераторы



Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group