2014 dxdy logo

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

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


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


Посмотреть правила форума



Начать новую тему Ответить на тему
 
 Численное интегрирование системы ОДУ методом Рунге-Кутта
Сообщение26.09.2016, 18:11 


14/04/14
36
Задача. Определить одним из численных методов состояние динамического объекта $ x(t)$, модель которого имеет вид: $\dot{x}=Ax,x_0=x(t_0)$ при $ t = 1 c$. Вектор начальных условий $x_0$ совпадает с одним из собственных векторов матрицы А.
$A=$\begin{bmatrix}
 1&  0& 1& \\
 2&  -1& 0& \\
 1&  -1& 2&
\end{bmatrix}$$
Точность определения $\varepsilon= 0.05$. Осуществить проверку аналитическим методом.

Я определил начальные условия(собственный вектор).$\lambda=\begin{pmatrix}
 2& \\
 -1& \\
 1& 
\end{pmatrix}
X=\begin{pmatrix}
 0& \\
 0.949& \\
 0.316& 
\end{pmatrix}$
Как посчитать данное ОДУ с помощью метода Рунге-Кутта.
Почему собственный вектор - это начальные условия? Куда в методе Рунге-Кутта вставлять матрицу А из уравнения Коши которое дано из условия задачи?

 Профиль  
                  
 
 Re: Численное интегрирование системы ОДУ методом Рунге-Кутта
Сообщение26.09.2016, 20:31 
Заслуженный участник
Аватара пользователя


23/07/08
10910
Crna Gora
Про $\TeX$. Знаете, почему у Вас закрывающие скобки векторов и матриц сдвинуты вправо? Потому что Вы в каждой строке ставите лишний знак & в конце.

Про задачу. Давайте для начала правильно найдём собственные значения и соответствующие им собственные векторы. Что такое $\lambda=\begin{pmatrix}2\\-1\\1\end{pmatrix}$ ?

 Профиль  
                  
 
 Re: Численное интегрирование системы ОДУ методом Рунге-Кутта
Сообщение26.09.2016, 21:42 


14/04/14
36
$\lambda_1=1, \lambda_2=-1, \lambda_3=2 $

 Профиль  
                  
 
 Re: Численное интегрирование системы ОДУ методом Рунге-Кутта
Сообщение26.09.2016, 21:47 
Заслуженный участник


09/05/12
25179
Redmal в сообщении #1154924 писал(а):
$\lambda_1=1, \lambda_2=-1, \lambda_3=2 $
Это неправильно.

 Профиль  
                  
 
 Re: Численное интегрирование системы ОДУ методом Рунге-Кутта
Сообщение26.09.2016, 22:32 
Заслуженный участник
Аватара пользователя


01/03/06
13626
Москва
Redmal в сообщении #1154853 писал(а):
Почему собственный вектор - это начальные условия?

Это требование написано в условии задачи:
Redmal в сообщении #1154853 писал(а):
Вектор начальных условий $x_0$ совпадает с одним из собственных векторов матрицы А.

Redmal в сообщении #1154853 писал(а):
Куда в методе Рунге-Кутта вставлять матрицу А из уравнения Коши которое дано из условия задачи?

Пацаны, кто знаИт, скАрее Сдесь пЕшЫте, чО куда поЦтавлять, ато миня учебник лЕстать лАмАит! :D

 Профиль  
                  
 
 Re: Численное интегрирование системы ОДУ методом Рунге-Кутта
Сообщение27.09.2016, 00:17 


14/04/14
36
Ребят я ошибся в первой строке 3-й элемент матрицы не 1, а 0. Пересчитывал, все верно.
$\begin{bmatrix}
 1&  0& 0\\
 2&  -1& 0\\
 1&  -1& 2
\end{bmatrix}$
$h=0.549$
Дальше выбираю один из $\lambda$ и нахожу вектор X, а дальше нужно выбрать шаг и решить с помощью Рунге-Кутта. Как это сделать?

 Профиль  
                  
 
 Re: Численное интегрирование системы ОДУ методом Рунге-Кутта
Сообщение27.09.2016, 00:42 
Заслуженный участник
Аватара пользователя


23/07/08
10910
Crna Gora
Redmal в сообщении #1155026 писал(а):
Дальше выбираю один из $\lambda$ и нахожу вектор X
Пожалуйста, всё-таки напишите, какое собственное значение Вы выбрали и какой ему соответствует собственный вектор.
Redmal в сообщении #1155026 писал(а):
дальше нужно выбрать шаг и решить с помощью Рунге-Кутта. Как это сделать?
За уравнения какого общего вида «берётся» метод Рунге-Кутта?

 Профиль  
                  
 
 Re: Численное интегрирование системы ОДУ методом Рунге-Кутта
Сообщение27.09.2016, 12:08 


14/04/14
36
Цитата:
Пожалуйста, всё-таки напишите, какое собственное значение Вы выбрали и какой ему соответствует собственный вектор.

Redmal в сообщении #1155026 писал(а):
дальше нужно выбрать шаг и решить с помощью Рунге-Кутта. Как это сделать?

Беру элемент $\lambda=-1$
Получаю вектор X, который писал в первом посте.

 Профиль  
                  
 
 Re: Численное интегрирование системы ОДУ методом Рунге-Кутта
Сообщение27.09.2016, 13:41 
Заслуженный участник
Аватара пользователя


23/07/08
10910
Crna Gora
Ага.
У такого выбора есть достоинства (точное решение ограничено) и недостатки (из-за ошибок округления может возникнуть небольшая, но экспоненциально растущая примесь других собственных векторов, которая вскоре забьёт экспоненциально убывающее точное решение).

Вы, наверное, знаете, что собственные векторы определены с точностью до умножения на скаляр $c\neq 0$. Если $Ax=\lambda x$, то $A(cx)=\lambda cx$. Пользуясь этим, Ваш вектор можно упростить до
$\begin{bmatrix}0\\3 \\1\end{bmatrix}$
(проверьте, что этот вектор тоже подходит в качестве собственного для $\lambda=-1$).

Остается вопрос про общую форму уравнения, к которой применим метод. Когда Вы её запишете, я скажу, где в этом выражении учитывается матрица.

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

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



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

Сейчас этот форум просматривают: mihaild


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

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