- •2. Интерполирование при помощи алгебраических полиномов.
- •3. Интерполяционный полином Лагранжа.
- •4. Интерполяционный полином Ньютона.
- •5. Оценка остаточного члена интерполяционного полинома.
- •6. Численное интегрирование
- •7. Применение полинома Лагранжа
- •8. Формулы прямоугольников
- •9. Формула трапеций
- •10. Формула Симпсона
- •11. Оценка погрешности квадратурных формул
- •16.Аппроксимация обыкновенных и частных производных с помощью разложения в ряд Тейлора и ее порядок.
- •17, 18. Некорректность численного дифференцирования
- •20. Метод прогонки решения систем линейных алгебраических уравнений с трехдиагональной матрицей
- •21. Постановка задач для уравнений в частных производных
- •22. Классификация уравнений в частных производных второго порядка
- •23, 24. Метод сеток (метод конечных разностей) решения краевых задач для уравнений второго порядка в частных производных
- •XI,yi,,.
- •25. Аппроксимация дифференциальных уравнений разностными
22. Классификация уравнений в частных производных второго порядка
Наиболее общее уравнение второго порядка с частными производными в случае двухнезависимых переменныхx,yимеет следующий вид
F(x, y, ux, uy, uxx, uxy, uyy) = 0
а в случае переменных x,t– следующийF(x, t, ux, ut, uxx, uxt, utt) = 0.
В первом случае неизвестной функцией является функция u = u(x,y), а во втором – функцияu = u(x,t).
Различие этих и других уравнений состоит лишь в обозначении переменных.
Уравнения второго порядка классифицируются по знаку выражения (дискриминанта) B2– AC:
B2– AC < 0 –эллиптическийтип,
B2– AC = 0 –параболическийтип,
B2– AC > 0 –гиперболическийтип.
Физические процессы, описываемые разными типами уравнений, существенно отличаются друг от друга. Соответственно, постановки задач для них имеют свои особенности.
Заметим, что уравнение с переменными коэффициентами может иметь разный тип в разных точках x,t.
Уравнения с постоянными коэффициентами приводятся к следующим каноническим формам
uxx+uyy+au=f(x,y) (эллиптическийтип),
utuxx+au=f(t,x) (параболическийтип)
uttuxx+au=f(t,x) (гиперболическийтип),
Будем рассматривать простейшие формы этих уравнений
uxx+uyy= 0 (уравнениеЛапласа),
utuxx= 0 (уравнениетеплопроводности)
uttuxx= 0 (волновоеуравнение),
Для приведенных простейших уравнений второго порядка имеем соответственно
AC< 0,AC= 0,AC> 0,
поэтому первое из них является эллиптическим, второе –параболическим, третье –гиперболическим.
Решения уравнения Лапласа называются гармоническими функциями. Например,u=ax+by+c,u=a(x2y2) +bx+cy+d.
Для уравнения Лапласа задается область с замкнутой границей и рассматриваются три вида краевых условий:
1-го рода (задача Дирихле) u=(x,y), (x,y),
2-го рода (задача Неймана) u/n=(x,y), (x,y),
и 3-го рода (смешанная краевая задача) (u+u/n)=(x,y), (x,y),,константы.
Для волнового уравнения и уравнения теплопроводности задаются начальные условия при t=t0.
При этом, если <x<, то получается задача Коши. Если жеa<x<b, то на концах отрезка [a,b] задаются граничные условия одного из трех указанных видов и задачаназывается начально-краевой или смешанной задачей Коши.
23, 24. Метод сеток (метод конечных разностей) решения краевых задач для уравнений второго порядка в частных производных
Рассмотрение метода конечных разностей начнем со следующей задачи Дирихле для уравнения Лапласа: найти решение u(x,y) уравнения(1)
в некоторой прямоугольной области G = [a, b] x [c, d] плоскости x, y при
, (2)
где функция φ (x, y) непрерывна на границе Г прямоугольника G.
Область Gс границей Г покроем прямоугольной сеткой. Для этого в плоскостиx,yпостроим два семейства параллельных прямых
; .
Здесь x0=a,y0=c– координаты левого нижнего угла прямоугольникаG,xI=b, yJ=d– координаты его правого верхнего угла,h= (xIx0)/Iиl= (yJy0)/J–шаги сеткипо направлениямx иy соответственно.
Узлами сетки являются точки пересечения указанных прямых, имеющие координаты