fixfix
2014 dxdy logo

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

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


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


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



Начать новую тему Ответить на тему
 
 Проблемы с численным решением дифференциального уравнения
Сообщение24.04.2023, 21:10 


04/03/23
12
Добрый вечер! Имеется следующее дифференциальное уравнение.
$k\dddot{y(x)}=\ddot{y(x)}-F(x)$, $F(x)$ - известная громоздкая функция, $k\sim10^{-20}$
Я знаю решение, удовлетворяющее данному диффуру. Хочу сравнить численное и аналитическое решения.
Однако решая данный диффур Рунге-Куттой 4, столкнулся с проблемой - решения не совпадают.
При небольшом изменении шага интегрирования решение вовсе "улетает".
Что это? Несходимость решения или какая-то неустойчивость?
Подскажите, как лучше всего это разрешить? Может, использовать какой-то другой численный метод?

 Профиль  
                  
 
 Re: Проблемы с численным решением дифференциального уравнения
Сообщение24.04.2023, 21:43 
Заслуженный участник


09/05/12
25179
Первая очевидная идея - использовать очень малый шаг по времени. Вторая - отмасштабировать задачу к характерным величинам для данной системы. Третья - использовать неявные методы. Два последних варианта можно (и рекомендуется) использовать вместе.

 Профиль  
                  
 
 Re: Проблемы с численным решением дифференциального уравнения
Сообщение24.04.2023, 21:44 
Заслуженный участник
Аватара пользователя


23/07/08
10910
Crna Gora
AlexFuser, почему Вы не ответили на мои вопросы в предыдущей теме?

 Профиль  
                  
 
 Re: Проблемы с численным решением дифференциального уравнения
Сообщение24.04.2023, 21:49 
Заслуженный участник


23/05/19
1330
Для начала - избавиться от константы. Не знаю, формат какой точности Вы используете для своих данных, но, скорее всего, эта константа просто все зануляет.

 Профиль  
                  
 
 Re: Проблемы с численным решением дифференциального уравнения
Сообщение24.04.2023, 22:05 
Заслуженный участник
Аватара пользователя


30/01/09
7243
Попробовать решить аналитически.
AlexFuser в сообщении #1590975 писал(а):
При небольшом изменении шага интегрирования решение вовсе "улетает".
Что это? Несходимость решения или какая-то неустойчивость?

Уравнение "жёсткое".

-- Пн апр 24, 2023 23:10:57 --

AlexFuser в сообщении #1590975 писал(а):
Хочу сравнить численное и аналитическое решения.

Пропустил этот момент при чтении условия.
мат-ламер в сообщении #1590983 писал(а):
Попробовать решить аналитически.

Эту рекомендацию снимаю.

 Профиль  
                  
 
 Re: Проблемы с численным решением дифференциального уравнения
Сообщение24.04.2023, 22:26 


04/03/23
12
Pphantom в сообщении #1590978 писал(а):
Первая очевидная идея - использовать очень малый шаг по времени. Вторая - отмасштабировать задачу к характерным величинам для данной системы. Третья - использовать неявные методы. Два последних варианта можно (и рекомендуется) использовать вместе.

К сожалению, уменьшение шага ни к чему хорошему не приводит. А неявные методы посмотрю. Спасибо!

-- 24.04.2023, 22:28 --

Dedekind в сообщении #1590980 писал(а):
Для начала - избавиться от константы. Не знаю, формат какой точности Вы используете для своих данных, но, скорее всего, эта константа просто все зануляет.

При освобождении производной наивысшего порядка от константы путем переноса в правую часть получаются решения 20-го порядка

-- 24.04.2023, 22:31 --

мат-ламер в сообщении #1590983 писал(а):
Попробовать решить аналитически.
AlexFuser в сообщении #1590975 писал(а):
При небольшом изменении шага интегрирования решение вовсе "улетает".
Что это? Несходимость решения или какая-то неустойчивость?

Уравнение "жёсткое".

-- Пн апр 24, 2023 23:10:57 --

AlexFuser в сообщении #1590975 писал(а):
Хочу сравнить численное и аналитическое решения.

Пропустил этот момент при чтении условия.
мат-ламер в сообщении #1590983 писал(а):
Попробовать решить аналитически.

Эту рекомендацию снимаю.


Спасибо, попробую неявные методы

 Профиль  
                  
 
 Re: Проблемы с численным решением дифференциального уравнения
Сообщение24.04.2023, 23:07 
Заслуженный участник


23/05/19
1330
AlexFuser в сообщении #1590986 писал(а):
При освобождении производной наивысшего порядка от константы путем переноса в правую часть получаются решения 20-го порядка

А нужно освобождать не путем переноса, а путем подходящего выбора единиц измерения.

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

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



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

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


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

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