2014 dxdy logo

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

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




 
 Алгоритм Рунге-Кутта 4 порядка
Сообщение20.01.2010, 18:30 
Добрый вечер (или день)!

В ходе решения возникла необходимость в программе вернуться на два шага назад.
Решается система 12 ОДУ метод РК 4-порядка.

Нигде не встречал условия на шаг $h$ по сетке для этого метода. Могу ли я положить $h<0$ тем самым (в цикле) вернуться на $2h$ шагов назад??

Нет ли какой скрытой симметрии в весовых коэффициентах $K_j, \, j=1,2,3,4$?

.Ульрих.

 
 
 
 Re: Алгоритм Рунге-Кутта 4 порядка
Сообщение20.01.2010, 21:03 
Чтоб вернуться назад, можно и не вычислять коэффиценты, а хранить в памяти значение по сетке для двух шагов :)

Можно ли шаг взять отрицательным?
Пусть дано уравнение(порядок не обязательно первый):
$\frac{d Y}{d x}=F(x,Y)\quad\eqno(1)$
Найдем его решение "справа налево". Для этого введем переменную $t=-x$ и напишем уравнение для $Y$ и $t$:
$\frac{d Y}{d t}=G(t,Y)$ , где $G(t,Y)=-F(-t,Y)$.
Его тоже решаем методом Рунге-Кутты.
А теперь вопрос: Эквивалентен ли алгоритм решения последнего уравнения алгоритму решения уравнения (1) с отрицательным шагом?

 
 
 
 Re: Алгоритм Рунге-Кутта 4 порядка
Сообщение21.01.2010, 01:10 
jetyb в сообщении #282055 писал(а):
Чтоб вернуться назад, можно и не вычислять коэффиценты, а хранить в памяти значение по сетке для двух шагов :)


Блин (упс).. вот я тоже думаю что можно (и нужно) хранить... Но мозг настолько настроен на сложное... что совершено не повернуть сознание в это русло ))
:lol:
Переклин))


Да, нет никакой разницы - как выглядит правая часть и какой шаг :roll: ...
Я думал для (для крепости) этого утверждения, есть специальная теорема (лемма)
:idea:

 
 
 [ Сообщений: 3 ] 


Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group