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 ] 

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



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

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


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

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