2014 dxdy logo

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

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




 
 Разностная схема. Краевые условия второго рода.
Сообщение19.12.2009, 15:26 
Как составлять разностную схемы для краевой задачи с условими второго рода?
Ну к примеру задача такая:
$
-(k(x)u')' + q(x)u = f(x),
u'(a) = uu_a,
u'(b) = uu_b,
$

 
 
 
 Re: Разностная схема. Краевые условия второго рода.
Сообщение19.12.2009, 15:44 
Хм... или я что-то не понимаю, или в условии что-то лишнее.

 
 
 
 Re: Разностная схема. Краевые условия второго рода.
Сообщение19.12.2009, 15:47 
Понятно, что одна пара условий лишняя (судя по заголовку -- первая). Непонятно другое: что понимается под $uu_a$ и $uu_b$?...

 
 
 
 Re: Разностная схема. Краевые условия второго рода.
Сообщение19.12.2009, 16:10 
ewert в сообщении #273015 писал(а):
Понятно, что одна пара условий лишняя (судя по заголовку -- первая). Непонятно другое: что понимается под $uu_a$ и $uu_b$?...

Просто числа

 
 
 
 Re: Разностная схема. Краевые условия второго рода.
Сообщение19.12.2009, 16:25 
Что -- числа?... Почему буква $u$ повторяется только два раза, а не восемнадцать?... И почему вообще именно эта буква?...

Ладно. Поскольку наверняка имеется в виду стандартная разностная схема, имеющая точность $O(h^2)$ -- чтобы эту точность не потерять, используйте для аппроксимации граничных условий трёхточечные разностные формулы: $u'(x_0)=\dfrac{-3u_0+4u_1-u_2}{2h}+O(h^2)$ и $u'(x_n)=\dfrac{u_{n-2}-4u_{n-1}+3u_0}{2h}+O(h^2)$.

Если покажется противным нарушение трёхдиагональности -- подковыряйте пальчиками два первых и два последних уравнения.

 
 
 
 Re: Разностная схема. Краевые условия второго рода.
Сообщение19.12.2009, 17:03 
ewert
Спасибо.
Только в последней формуле наверное вместо $u_0$ должно быть $u_n$.

 
 
 
 Re: Разностная схема. Краевые условия второго рода.
Сообщение19.12.2009, 18:47 
да уж наверное, наспех перередактировал, пардон

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


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