4. Краевые задачи для уравнений второго порядка

При математическом описании реальных физических объектов чаще всего приходится иметь дело с дифференциальными уравнениями в обыкновенных или частных производных второго порядка с начальными, краевыми или граничными условиями.

Преобразование их в конечно-разностную систему алгебраических уравнений осуществляется аналогично: для каждой точки в области (интервале) интегрирования, где не задано краевое или граничное значение искомой функции, записывается исходное уравнение, в котором все производные выражены через заранее определенное число близлежащих ординат искомой функции, принадлежащих области, и вычислены все коэффициенты и функции независимых переменных в этой точке. К полученным таким образом уравнениям добавляются соотношения или значения функции и ее производных в точках границы области. В результате будет сформирована алгебраическая система уравнений с числом уравнений и неизвестных, равном общему числу точек области интегрирования.

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

Если в центральных точках точность аппроксимации производных с n точками удовлетворяет поставленным требованиям и эту точность желательно сохранить и в приграничных точках заданных областей, то для последних выбирают аппроксимирующие формулы, построенные для (n+1)-й точки или более.

Рассмотрим примеры аппроксимации дифференциальных уравнений с краевыми условиями конечно-разностной системой алгебраических уравнений. Эти аппроксимации в литературе получили название "разностные схемы". Ниже в четырех таблицах приведены четыре варианта конечно-разностной аппроксимации одной и той же краевой задачи, для которой известно точное решение. Вид уравнения, условия на границе интервала, решение аналитическое и вычисленное в заданных точках с 12 значащими цифрами приведены в правой крайней колонке первой таблицы. В левых колонках первой и в трех остальных таблицах записаны системы алгебраических уравнений, полученных применением трех-, пяти-, пяти-шести- и семи точечной аппроксимации второй производной в заданном уравнении. Справа от уравнений приведены решения алгебраических уравнений тоже с 12-ю значащими цифрами.

Система уравнений с трехточечным представлением производных

Вектор разностного решения с шагом h=0.1

-199+100+0.1=0

0.0186590989712 0.0186415437361

100-199+100+0.2=0

0.0361316064473 0.0360976603850

100-199+100+0.3=0

0.0512427953890 0.0511947672548

100-199+100+0.4=0

0.0628415300546 0.0627828520998

100-199+100+0.5=0

0.0698118753674 0.0697469636621

100-199+100+0.6=0

0.0710840847137 0.0710183518969

100-199+100+0.7=0

0.0656455142231 0.0655851465687

100-199+100+0.8=0

0.0525504484304 0.0525024675253

100-199+0.9=0

0.0309298757856 0.0309018656257


Система уравнений для пяти-точечного

представления производных

Вектор решения

-9940+3000+2000-500+6=0

0.0186406186406

8000-14940+8000-500+12=0

0.0360968696594

-500+8000-14940+8000-500+18=0

0.0511941848390

-500+8000-14940+8000-500+24=0

0.0627825213460

-500+8000-14940+8000-500+30=0

0.0697468774179

-500+8000-14940+8000-5008+36=0

0.0710184988305

-500+8000-14940+8000-500+42=0

0.0655854996422

-500+8000-14940+8000+48=0

0.0525029672554

-500+2000+3000-9940+54=0

0.0309024932693

 

Система уравнений для пяти- и шести точечного представления производных Вектор решения

-3720-1000+3500-1500+250+3=0

0.0186415486274

8000-14940+8000-500+12=0

0.0360976918947

-500+8000-14940+8000-500+18=0

0.0511948294923

-500+8000-14940+8000-500+24=0

0.0627829167486

-500+8000-14940+8000-500+30=0

0.0697469746974

-500+8000-14940+8000-500+36=0

0.0710183243686

-500+8000-14940+8000-500+42=0

0.0655851063829

-500+8000-14940+8000+48=0

0.0525024168959

250-1500+3500-1000-3720+27=0

0.0309018105849
Система уравнений для семиточечного представления производных Вектор решения

-7260-12750+23500-14250+4650-650+9=0

0.0186415513486

11400-20910+10000+750-600+100+18=0

0.0360976659970

-1350+13500-24410+13500-1350+100+27=0

0.0511947713313

10-135+1350-2441+1350-135+10+3.6=0

0.0627828547351

10-135+1350-2441+1350-135+10+4.5=0

0.0697469648318

10-135+1350-2441+1350-135+10+5.4=0

0.0710183515790

100-1350+13500-24410+13500-1350+63=0

0.0655851447467

100-600+750+10000-20910+11400+72=0

0.0525024640963

-650+4650-14250+23500-12750-7260+81=0

0.0309018602217

В этой задаче весь интервал интегрирования [0,1] был разбит на 10 равных частей с шагом h=0.1. Из одиннадцати точек в двух крайних искомая функция x(t) была задана, поэтому уравнения записывались для девяти внутренних точек, в которых значения функции требовалось найти.


Информация о работе «Краевые задачи и разностные схемы»
Раздел: Математика
Количество знаков с пробелами: 26705
Количество таблиц: 20
Количество изображений: 17

Похожие работы

Скачать
74851
18
9

... схема со сглаживанием, схема прямоугольник со сглаживанием, “шахматная ” схема. Произведены некоторые расчеты для одномерного уравнения переноса с переменными и постоянными коэффициентами на неравномерных сетках, с целью определения наиболее устойчивой разностной схемы. Исследование показало, что наиболее устойчивым методом для одномерного уравнения переноса с переменными коэффициентами является ...

Скачать
20167
4
9

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

Скачать
57698
75
8

... ^у^е^о ^ с^-^. Итак решение по Ритцу: ^-i-^ Сравнительная таблица имеет вид: Л. 0 0,5 1 1,5 2 у^ 0 -0,275 -0,3571 -0,2758 0 ^г) о -0,2126 -0,3520 -0,3258 0 50 3.6. Об одном подходе к решению нелинейных вариационных задач В отличии от метода Ритца, искомую функцию в двуточечной вариа­ционной задаче зададим в виде: r-^^f^-^^ При этом граничные условия и{а ) = ...

Скачать
39796
9
22

... суммы и позволит вычислить приближенное значение приращения Dy:  где Метод четвертого порядка для q = 3, имеет вид  где Особо широко известно другое вычислительное правило Рунге-Кутта четвертого порядка точности:  где Метод Рунге-Кутта имеет погрешность четвертого порядка (~ h4 ). Правило Рунге. Если приближенный метод имеет порядок погрешности m, то погрешность можно приближенно ...

0 комментариев


Наверх