2014 dxdy logo

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

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


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


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



Начать новую тему Ответить на тему
 
 Устойчивость метода Эйлера
Сообщение19.12.2012, 17:30 


19/12/12
2
Требуется рассмотреть устойчивость этого метода.
Что у нас есть - задача Коши, примененный к ней метод Эйлера.

$\frac{du}{dt}=f(u,t)$

$u_{n+1}=u_n+f_n\cdot\Delta t$

Ошибка на n-ном шаге (то есть отклонение от истинного значения) равна $\varepsilon_n$

Далее расписывается истинное значение функции

$u_{n+1}+\varepsilon_{n+1}=u_n+\varepsilon_n+f_n(u_n+\varepsilon_n)\Delta t = u_n+\varepsilon_n + f_n(u_n)\Delta t+\frac{df}{du}\varepsilon_n\Delta t+o(\varepsilon_n^{2})\Delta t$

Здесь как раз и появляется тот самый момент, который я не до конца понимаю. А именно - раскрытие функции:
$f_n(u_n+\varepsilon_n)=f_n(u_n)+\frac{df}{du}\varepsilon_n + o({\varepsilon_{n}}^{2})$

Пробовал раскладывать по Тейлору с центром разложения в точке $u_n$, получается следующее:
$f_n(u_n+\varepsilon_n)=f_n(u_n)+\frac{df}{du}\varepsilon_n + \frac{d^{2}f}{du^{2}}{\varepsilon_{n}}^{2}+o({\varepsilon_{n}}^{2})$

То разложение, что написано выше, взято с материалов лекции, так что ошибки быть не должно. Почему получается член второго порядка или как от него избавиться? Что я делаю не так?

 i  Deggial: Инструкции по набору формул здесь или здесь (или в этом видеоролике). Исправляйте формулы, иначе тема будет перемещена в Карантин.

 Профиль  
                  
 
 Re: Устойчивость метода Эйлера
Сообщение19.12.2012, 18:07 
Заслуженный участник
Аватара пользователя


23/08/07
5494
Нов-ск
$u_{n+1}=u_n+f_n\cdot\Delta t$

$u_{n+1}+\varepsilon_{n+1}=u_n+\varepsilon_n+f_n(u_n+\varepsilon_n)\Delta t $

Какому из этих уравнений удовлетворяет функция $u_n$?

Исправьте
$f_n(u_n+\varepsilon_n)=f_n(u_n)+\frac{df}{du}\varepsilon_n + o({\varepsilon_{n}}^{2})$
на
$f_n(u_n+\varepsilon_n)=f_n(u_n)+\frac{df}{du}\varepsilon_n + O({\varepsilon_{n}}^{2})$

 Профиль  
                  
 
 Re: Устойчивость метода Эйлера
Сообщение19.12.2012, 19:19 


19/12/12
2
Это одна и та же функция в обоих уравнениях, то есть и первому, и второму, если я правильно понимаю вопрос.
Исправление, конечно, поправит ситуацию именно так, как и надо, но суть вопроса именно в том и состояла - почему в этом месте стоит о малое, или же это утверждение неверно?

 Профиль  
                  
 
 Re: Устойчивость метода Эйлера
Сообщение19.12.2012, 19:26 
Заслуженный участник
Аватара пользователя


23/08/07
5494
Нов-ск
BinGOSU в сообщении #660767 писал(а):
Это одна и та же функция в обоих уравнениях, то есть и первому, и второму, если я правильно понимаю вопрос.
Исправление, конечно, поправит ситуацию именно так, как и надо, но суть вопроса именно в том и состояла - почему в этом месте стоит о малое, или же это утверждение неверно?

Неправильно понимаете, не удовлетворяет она сразу двум уравнеиям. (Для значений точного решения в узлах сетки следовало бы ввести свое обозначение.)
С о-малым равенство неверное.

 Профиль  
                  
 
 Re: Устойчивость метода Эйлера
Сообщение20.12.2012, 03:34 
Заслуженный участник


11/05/08
32166
BinGOSU в сообщении #660694 писал(а):
$u_{n+1}+\varepsilon_{n+1}=u_n+\varepsilon_n+f_n(u_n+\varepsilon_n)\Delta t = u_n+\varepsilon_n + f_n(u_n)\Delta t+\frac{df}{du}\varepsilon_n\Delta t+o(\varepsilon_n^{2})\Delta t$

Тут вообще всё даже не то что неверно, а просто бессмысленно. Бог с ним, с о-маленьким. Бессмысленна уже сама запись $f_n(u_n+\varepsilon_n)$.

 Профиль  
                  
 
 Re: Устойчивость метода Эйлера
Сообщение20.12.2012, 07:04 
Заслуженный участник
Аватара пользователя


23/08/07
5494
Нов-ск
ewert в сообщении #660918 писал(а):
BinGOSU в сообщении #660694 писал(а):
$u_{n+1}+\varepsilon_{n+1}=u_n+\varepsilon_n+f_n(u_n+\varepsilon_n)\Delta t = u_n+\varepsilon_n + f_n(u_n)\Delta t+\frac{df}{du}\varepsilon_n\Delta t+o(\varepsilon_n^{2})\Delta t$

Тут вообще всё даже не то что неверно, а просто бессмысленно. Бог с ним, с о-маленьким. Бессмысленна уже сама запись $f_n(u_n+\varepsilon_n)$.

Смысл очень легко возникает, если не путаться в обозначения. Точное решение (дифференциального уравнения) в узлах обозначить $v_n,$ затем $u_n=v_n+ \varepsilon_n$ подставлять в дискретное уравнение.

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

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



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

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


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

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