2014 dxdy logo

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

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




 
 оценка аппроксимации
Сообщение27.03.2011, 19:32 
Аппроксимирует ли разностная схема $ \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 
YanaKass в сообщении #428140 писал(а):
А как найти оценку граничного условия?

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

 
 
 
 Re:
Сообщение27.03.2011, 22:55 
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 
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 
Спасибо большое=) Все получилось

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


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