2014 dxdy logo

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

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


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


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

Если Вы хотите задать новый вопрос, то не дописывайте его в существующую тему, а создайте новую в корневом разделе "Помогите решить/разобраться (М)".

Если Вы зададите новый вопрос в существующей теме, то в случае нарушения оформления или других правил форума Ваше сообщение и все ответы на него могут быть удалены без предупреждения.

Не ищите на этом форуме халяву, правила запрещают участникам публиковать готовые решения стандартных учебных задач. Автор вопроса обязан привести свои попытки решения и указать конкретные затруднения.

Обязательно просмотрите тему Правила данного раздела, иначе Ваша тема может быть удалена или перемещена в Карантин, а Вы так и не узнаете, почему.



Начать новую тему Ответить на тему
 
 оценка аппроксимации
Сообщение27.03.2011, 19:32 


27/03/11
3
Аппроксимирует ли разностная схема $ \frac{u_{n+1}-u_{n-1}   }{2h}+u_n = nh+1$, $n=1,..., N-1, h=1/n$
$u_0=0, u_1=0$; дифференциальную задачу $ \frac{du}{dx}+u=x+1$, $x$ из $[0;1]$, $u(0)=0$ со вторым порядком по $h$? Если нет, то видоизменить разностную схему так, чтобы она имела второй порядок аппроксимации.

Само уравнение разностной схемы уже аппроксимировано. У него порядок $о(h^2)$. А как найти оценку граничного условия?

 Профиль  
                  
 
 
Сообщение27.03.2011, 21:30 
Заслуженный участник


11/05/08
32166
YanaKass в сообщении #428140 писал(а):
А как найти оценку граничного условия?

По общему правилу: недостающие начальные значения (в Вашем случае это $u_1$) надо аппроксимировать с тем же порядком точности, что и глобальный порядок точности метода. Т.е. достаточно для их нахождения использовать метод на единицу меньшего порядка, что у основного метода. Соответственно, в Вашем случае достаточно использовать для нахождения $u_1$ метод Эйлера (случайно, кстати, решение получится абсолютно точным). Условие же $u_1=0$ откровенно имеет порядок погрешности $O(h)$ (даже не вдаваясь в детали -- просто потому, что $u_0=0$, а производная в нуле ненулевая) и, следовательно, не подходит.

 Профиль  
                  
 
 Re:
Сообщение27.03.2011, 22:55 


27/03/11
3
ewert в сообщении #428175 писал(а):
достаточно использовать для нахождения $u_1$ метод Эйлера (случайно, кстати, решение получится абсолютно точным).

Если я не ошибаюсь, расчет там по формуле $u_{i+1}=u_i+hf(x_i,u_i)$
А что брать за $f(x_i,u_i)$?

 Профиль  
                  
 
 Re: Re:
Сообщение28.03.2011, 17:38 


13/11/09
166
YanaKass в сообщении #428224 писал(а):
Если я не ошибаюсь, расчет там по формуле $u_{i+1}=u_i+hf(x_i,u_i)$
А что брать за $f(x_i,u_i)$?


$u_{1}=u_0+h u'_0 + O(h^2)$. А далее используют дифференциальное уравнение, верное и на границе.

 Профиль  
                  
 
 
Сообщение03.04.2011, 17:07 


27/03/11
3
Спасибо большое=) Все получилось

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

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



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

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


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

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