1329
.pdfРассматриваемая задача может быть решена с использованием локально-одномерного метода. Однако в этом случае для достижения необходимой точности приходится уменьшать величину временного интервала и увеличивать число итераций, так как простейшие схемы метода обладают первым порядком точности.
6. ВАРИАЦИОННЫЕ МЕТОДЫ РЕШЕНИЯ КРАЕВЫХ ЗАДАЧ
Конечно-разностные методы предполагают замену дифференциального уравнения системой алгебраических путём аппроксимации дифференциальных операторов конечно-разностными выражениями. Однако такой подход не является единственно возможным. Альтернативный способ получения системы алгебраических уравнений – представление решения краевой задачи в виде разложения по системе базовых функций, удовлетворяющих заданным граничным условиям.
Аналогичный случай рассматривался ранее при решении уравнения Пуассона методом разделения переменных. Однако в отличие от этого метода собственные значения и собственные функции дифференциального оператора в общем случае неизвестны. Поэтому коэффициенты разложения не могут быть рассчитаны с использованием простейших выражений. Их определение связано с решением системы алгебраических уравнений, получаемой для различных методов из определённых условий.
Показано [32, 33, 34], что решение краевой задачи с заданными граничными условиями эквивалентно поиску минимума функционала определённого вида. Эти вопросы рассматриваются разделом математики, называемым вариационным исчислением. Напомним, что функционалом называется объект математики, когда каждой из множества функций, принадлежащих этому объекту, может быть поставлено в соответствие определённое число.
131
Для решения задач вариационного исчисления, в частности нахождения минимума функционала, наиболее часто используются приближённые методы, к которым относятся методы Ритца и Бубнова– Галёркина. Эти методы так же, как и конечно-разностные, приводят к системе алгебраических уравнений, решаемых известными методами. Однако в отличие от последних методов система уравнений находится из условия:
−минимума функционала, если используется метод Ритца;
−ортогональности невязки уравнения к базовым функциям разложения, если используется метод Галёркина.
Метод Ритца применим для решения краевых задач с самосопряжёнными и положительно определёнными операторами дифференциальных уравнений.
Напомним, что самосопряжённым уравнением называется дифференциальное уравнение, записываемое в виде
d |
dy |
− qy = f (x) . |
|
||
|
p |
|
|
(6.1) |
|
|
|
||||
dx |
dx |
|
|
К самосопряжённому виду путём определённых преобразований может быть приведено любое линейное дифференциальное уравнение второго порядка [34]
p |
d 2y |
+ r |
dy |
− qy = f (x) . |
(6.2) |
||
d x |
2 |
dx |
|||||
|
|
|
|
Дифференциальное уравнение считается положительно определённым, когда все собственные значения дифференциального оператора уравнения положительны.
Идея метода Ритца сводится к следующему. Пусть функционал J (u) в области D с граничными условиями на контуре Г, ограничи-
вающем эту область,
u = φ(x, y) на Г |
(6.3) |
записывается в виде
132
J (u) |
= ∫∫ |
x |
y |
(6.4) |
|
F (x, y,u,u′ |
,u′ )dxdy . |
D
Как указывалось выше, решение краевой задачи соответствует поиску минимума этого функционала.
Будем считать, что решение u представляется в виде разложения по системе пробных функций различного вида, удовлетворяющих заданным граничным условиям. Чаще всего в качестве пробных функций используют степенные или тригонометрические ряды с неопределёнными коэффициентами:
u(x, y) = Ф(x, y, c1 , c2 , ..., cn ) . |
(6.5) |
Подставляя указанное разложение (6.5) в выражение функционала (6.4), выполняя дифференцирование, интегрирование и другие, соответствующие функционалу математические преобразования, получим выражение, представляющее определённую функцию коэффициентов разложения c1 , c2 , ..., cn :
J = J (c1 , c2 , ..., cn ) . |
(6.6) |
Для нахождения минимума функционала необходимо приравнять производные полученного выражения по искомым коэффициентам нулю:
∂
=0 , k =1, 2, ..., n, (6.7)
∂ck
врезультате чего получаем относительно их систему алгебраических уравнений. Решение этой системы позволяет определить коэффициенты разложения и восстановить решение краевой задачи (6.5).
Вряде конкретных задач решение может быть упрощено, если пробные функции являются линейными относительно искомых коэффициентов разложения, вследствие чего уравнения системы (6.7) оказываются первой степени, и для получения необходимой точности число уравнений может быть относительно небольшим.
Вметоде Галёркина решение краевой задачи также ищется в ви-
де разложения (6.5), а система алгебраических уравнений определяется из условия ортогональности невязки уравнения краевой задачи.J
133
Если, например, уравнение краевой задачи с однородными граничными условиями записано в виде
L(u) = f (x, y) , |
(6.8) |
решение краевой задачи представлено в виде разложения
|
|
|
n |
|
|
|
|
|
u(x, y) = ∑ciφi(x, y) , |
(6.9) |
|
|
|
|
i =1 |
|
|
где L − дифференциальный оператор; |
φi(x, y) − выбранная система |
||||
базовых функций, то условие ортогональности имеет вид |
|
||||
|
|
∫∫[Lu(x, y) − f (x, y)]φi(x, y)dxdy = |
|
||
|
|
|
D |
|
|
|
|
n |
|
||
|
|
= ∫∫ L∑c jφj(x, y) − f (x, y) φi(x, y)dxdy = 0, |
|
||
|
|
D |
j =1 |
|
|
|
|
|
i =1, 2, ..., n, |
|
(6.10) |
или |
|
|
|
|
|
|
|
n |
|
|
|
∫∫ L |
∑c jφj |
(x, y) φi(x, y)dxdy = ∫∫ f (x, y)φi(x, y)dxdy . |
(6.11) |
||
D |
|
j =1 |
|
D |
|
Опять получаем систему алгебраических уравнений, позволяющую определить коэффициенты разложения сj и восстановить искомую функцию u(x, y) .
Недостатком рассмотренных «классических» способов решения краевых задач является сложность их реализации для получения необходимой точности. Поэтому в настоящее время эти методы совершенствуются с целью упрощения построения системы алгебраических уравнений. Основная идея совершенствования – использовать в качестве базовых функций функции с финитным (конечным) носителем, т.е. функции, отличные от нуля лишь в окрестности рассматриваемой точки. В этом случае число членов каждого уравнения системы алгебраических уравнений для коэффициентов разложения будет определяться числом точек, принадлежащих финитному
134
носителю базовой функции, а члены уравнения, соответствующие точкам за пределами носителя, будут равны нулю.
Примем для упрощения область определения решения краевой задачи 0 ≤ x≤ L одномерной. Разобьём указанную область 0 ≤ x≤ L на конечное число интервалов N и для каждого k-го интервала вве-
дём функцию |
|
|
|
|
|
|
0, 0 ≤ x≤ xk −1 , |
|
|||||
|
|
|
|
|
x≤ xk , |
|
ω1(x), xk −1 ≤ |
||||||
ωk (x) = |
|
|
|
≤ |
|
(6.12) |
ω2(x), xk |
x≤ xk +1 , |
|||||
0, |
x |
k +1 |
≤ |
x≤ |
x = |
1. |
|
|
|
|
N |
|
Эта функция отлична от нуля на интервале xk −1 ≤ x≤ xk +1 , состоит из двух линейных участков и достигает максимума, равного единице при x = xk (рис. 6.1).
Рис. 6.1. Базовая функция с финитным носителем
Подобно тому, как указывалось в п. 1.2, любую непрерывную кусочно-линейную финитную функцию можно представить в виде разложения по системе единичных базовых функций, причём в качестве единичных можно выбрать
135
k |
|
φ(x) = ∑φk ωk (x) , |
(6.13) |
1 |
|
где коэффициенты φk равны значениям функции φ(x) в точках xk .
Рассмотрим свойства базовой функции (6.12), определив её скалярное произведение:
1
∫ωk
0
0, |
n ≤ |
|
|
k− |
2, |
|
||
|
|
|
|
|
|
|
|
|
|
xk − xk −1 |
|
, n = k −1, |
|
||||
|
|
|||||||
|
6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
xk +1 − xk −1 |
|
|
|
|||||
(x)ωn(x)dx = |
|
|
|
|
|
, |
n = k, |
(6.14) |
|
3 |
|
|
|
||||
|
|
|
|
|
|
|
||
xk +1 − xk |
|
|
, |
|
n = k +1, |
|
||
|
|
|
|
|
|
|||
|
6 |
|
|
|
||||
|
k+ |
2. |
|
|||||
0, |
n ≥ |
|
||||||
|
|
|
|
|
|
|
|
Из выражения (6.14) следует, что функция ωk (x) ортогональна всем функциям ωn (x) , кроме ωk −1 , ωk , ωk +1 . Поэтому алгебраические уравнения, описывающие коэффициенты разложения φk , будут со-
держать всего три члена. Таким образом, использование базовых функций с финитными носителями позволяет значительно упростить структуру получаемой системы уравнений и применить для её решения рассмотренные ранее методы.
6.1. ВАРИАЦИОННО-РАЗНОСТНЫЙ МЕТОД РИТЦА
Методы Ритца, использующие функции с финитными носителями, принято обозначать как вариационно-разностные методы. Рассмотрим применение метода Ритца для решения одномерной краевой задачи, записываемой в виде
− |
d |
p(x) |
du |
+ q(x)u = f (x) |
(6.15) |
dx |
|
||||
|
|
dx |
|
с однородными граничными условиями u(0) = 0, u(1) = 0 .
136
Решение краевой задачи соответствует поиску минимума функционала [33, 34]
1
J (x) = ∫ p(x)
0
du 2 |
|
2 |
|
|
||
|
|
|
+ q(x)u |
|
− 2 f (x)u dx . |
(6.16) |
|
|
|||||
dx |
|
|
|
|
||
|
|
|
|
|
|
|
Исследуемая область разбивается на N интервалов − конечных элементов, и решение краевой задачи ищется в виде разложения по системе базовых функций
|
|
N −1 |
|
|
|
|
|
|
u(x) = ∑u k ωk (x) , |
|
(6.17) |
||||||
|
|
k =1 |
|
|
|
|
|
|
где |
|
|
|
|
|
|
|
|
x − xk −1 |
, xk −1 |
≤ x≤ |
xk , |
|
||||
|
|
|
|
|
||||
|
|
|
|
|||||
xk − xk −1 |
|
|
|
|
|
|||
xk +1 − x |
|
≤ |
x≤ |
|
|
|||
ωk (x) = |
|
|
|
, xk |
xk +1 , |
(6.18) |
||
|
|
|
||||||
xk +1 − xk |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
в остальных точках. |
|
||||||
0, |
|
В этом случае минимизируемый функционал может быть представлен суммой отдельных составляющих, записанных для каждого интервала разбиения пространственной координаты:
N −1
J (u) = ∑
k =1
x k +1
∫ p(x) x k −1
du 2 |
|
2 |
|
|
||
|
|
|
+ q(x)u |
|
− 2 f (x)u dx . |
(6.19) |
|
|
|||||
dx |
|
|
|
|
||
|
|
|
|
|
|
|
Запишем уравнение для одной составляющей функционала в точке xi на интервале xi −1 ≤ x≤ xi+1 . Искомое решение на этом интервале
содержит три члена и записывается в виде |
|
u(i) = ui−1ωi−1 + uiωi + ui+1ωi+1 , |
(6.20) |
137
где
ωi −1 = |
|
xi − x |
|
|||||
xi − xi−1 |
|
|||||||
|
|
|
|
|||||
|
x − xi−1 |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
− xi −1 |
|
||||
ωi = |
xi |
|
||||||
|
|
|
|
|
|
|||
|
xi+1 − x |
|
||||||
|
|
x |
+1 |
− xi |
|
|
||
|
i |
|
|
|
|
|||
ωi +1 |
= |
|
x − xi |
|
||||
xi +1 − xi |
|
|||||||
|
|
|
|
на интервале xi −1 ≤ x≤ xi ,
на интервале xi −1 ≤ x≤ xi+1 ,
на интервале xi ≤ x≤ xi+1 .
(6.21)
(6.22)
(6.23)
Интеграл, входящий в выражение (6.19), может быть представлен в виде суммы двух интегралов для отдельных участков интервала
xk +1 |
|
du 2 |
|
|||
∫ |
p(x) |
|
|
+ q(x)u 2 |
− 2 f (x)u dx = |
|
|
||||||
xk −1 |
|
dx |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
xk |
|
du 2 |
2 |
|
|
|
|
|
||||||||||
|
|
|
|
= ∫ |
p(x) |
|
|
|
|
+ q(x)u |
|
|
− 2 f (x)u dx + |
(6.24) |
|||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
dx |
|
|
|
|
|
|
|
|||||||
|
|
|
|
xk −1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
xk +1 |
du 2 |
|
2 |
|
|
|
|
|
|||||||||||
|
|
|
|
+ ∫ |
p(x) |
|
|
|
|
+ q(x)u |
|
− 2 f (x)u dx. |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
xk |
|
dx |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Производные |
dω |
в этих интегралах: |
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
d ωi−1 |
= |
−1 |
|
|
на интервале x |
≤ |
x≤ |
x , |
(6.25) |
||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
dx |
|
|
xi − xi −1 |
|
|
|
|
|
i −1 |
|
|
i |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
ωi |
|
|
|
xi − xi 1 |
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
= |
|
|
− |
на интервале xi −1 ≤ |
|
x≤ |
xi+1 , |
(6.26) |
|||||||||||
|
|
|
|
|
−1 |
|
|||||||||||||||||
|
|
dx |
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
xi+1 − xi |
|
|
|
|
|
|
|
|
|
|
138
|
|
dω |
i +1 |
= |
|
1 |
|
|
на интервале x ≤ x≤ |
x |
. |
(6.27) |
||||||
|
|
|
|
|
|
|
|
|||||||||||
|
|
dx |
|
|
xi +1 − xi |
|
|
i |
i+1 |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Обозначив первый интеграл (6.24) как J1 (u) |
и подставив в него |
|||||||||||||||||
выражения (6.20), будем иметь: |
|
|
|
|
|
|
|
|
||||||||||
xi |
|
|
|
|
|
|
|
|
|
2 |
|
xi − x |
|
|
|
x − xi−1 |
2 |
|
J1(u) = ∫ |
p(x) − |
ui−1 |
|
+ |
|
ui |
|
+ q(x) |
ui−1 + |
ui + |
||||||||
|
|
|
|
|
xi − xi−1 |
|||||||||||||
|
|
|
|
xi − xi−1 |
xi − xi−1 |
xi − xi−1 |
|
|
|
|||||||||
xi−1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
xi − x |
|
|
x − xi−1 |
|
|
|
+2 f (x) |
ui−1 + |
ui dx. |
|||||
− |
|
||||||
|
|
− |
|
||||
|
xi − xi 1 |
|
xi − xi 1 |
|
Аналогично второй интеграл (6.24)
xi+1 |
|
|
|
|
|
|
2 |
|
− x |
|
|
J2 (u) = ∫ |
p(x) |
− |
ui |
+ |
ui+1 |
|
+ q(x) |
xi+1 |
ui + |
||
xi+1 − xi |
|
|
|
||||||||
xi |
|
|
|
|
xi+1 − xi |
xi+1 − xi |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
xi+1 − x |
|
x − xi |
|
|
+2 f (x) |
ui + |
ui+1 dx. |
|||
|
xi+1 − xi |
||||
xi+1 − xi |
|
|
(6.28)
x − xi |
ui+1 2 |
+ |
|
xi+1 − xi |
|||
|
|
(6.29)
Условием минимума функционала является равенство нулю его производных по искомым величинам, в данном случае
d |
( J 1 + J 2 ) = 0 . |
(6.30) |
|
||
dui |
|
Дифференцируя подынтегральное выражение по ui и группируя подобные члены, запишем уравнение (6.30) в виде
|
|
|
ai,i−1ui−1 − ai,iui + ai,i+1ui+1 = f i |
, |
|
|
(6.31) |
|||||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
xi |
|
− p(x) |
|
|
( x |
− x)( x |
− x |
) |
|
|||
ai,i −1 |
= |
∫ |
|
|
+ q(x) |
i |
|
|
|
i−1 |
dx , |
(6.32) |
||
|
|
|
|
( x |
|
|
)2 |
|||||||
( x |
− x |
)2 |
|
− x |
|
|||||||||
|
|
|
|
|
|
|
|
|||||||
|
|
xi−1 |
|
i |
i−1 |
|
|
|
i |
i−1 |
|
|
|
139
|
|
|
|
xi |
|
|
|
p(x) |
|
|
|
|
|
|
( x − x |
−1 |
)2 |
|
|
||||||||||
ai,i = − |
∫ |
|
|
|
|
|
|
|
|
+ q(x) |
|
|
|
i |
|
|
dx − |
|
|||||||||||
( x |
|
− x |
)2 |
( x |
− x |
|
)2 |
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
xi−1 |
|
|
i |
i−1 |
|
|
|
|
|
|
|
i |
|
i−1 |
|
|
|
|
(6.33) |
|||||||
|
|
|
xi+1 |
|
|
|
|
|
|
|
|
|
|
( x |
|
− x)2 |
|
|
|
||||||||||
|
|
|
|
p(x) |
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
− |
∫ |
|
|
|
|
|
|
|
|
+ q(x) |
|
|
i+1 |
|
|
|
|
dx, |
|
|
||||||||
|
|
|
|
|
|
− x )2 |
( x − x |
|
)2 |
|
|
|
|||||||||||||||||
|
|
|
( x |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
xi |
|
|
i +1 |
|
i |
|
|
|
|
|
|
i |
|
i −1 |
|
|
|
|
|
|
|
|||||
|
|
xi+1 |
|
− p(x) |
|
|
|
( |
xi +1 |
− x)( x − x ) |
|
||||||||||||||||||
ai,i +1 = |
|
∫ |
|
|
|
|
+ q(x) |
|
|
|
|
|
|
|
|
i |
dx , |
(6.34) |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
( x |
|
− x )2 |
|
|
||||||||||||||
|
|
|
|
|
− x )2 |
|
|
|
|
|
|
||||||||||||||||||
|
|
|
( x |
|
|
|
|
|
|
|
|
+1 |
|
|
|
|
|||||||||||||
|
|
|
xi |
|
i+1 |
|
i |
|
|
|
|
|
|
|
i |
|
|
i |
|
|
|
|
|||||||
|
|
|
xi |
|
|
|
|
|
|
|
x1 − x |
|
|
xi+1 |
|
|
|
xi +1 − x |
|
|
|
||||||||
f i |
= ∫ |
f (x) |
|
|
|
dx + |
∫ |
|
f (x) |
dx . |
(6.35) |
||||||||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
xi−1 |
|
|
|
xi − xi −1 |
|
xi |
|
|
|
xi+1 − xi |
|
|
||||||||||||||
Рассчитав |
коэффициенты |
уравнений |
|
ai,i |
|
для i =1, 2, ..., N −1 , |
|||||||||||||||||||||||
ai,i −1 для i = 2, 3, ..., N −1 , |
|
ai,i +1 |
для i =1, 2, ..., N − 2 , правую часть f i |
и используя известные граничные условия, получим систему трёхчленных алгебраических уравнений
Au = f , |
(6.36) |
где A − матрица с элементами ai,j; u − вектор неизвестных; |
f − век- |
тор с рассчитанными компонентами fi.
Решив указанную систему методом прогонки, определим значения коэффициентов ui, после чего рассчитывается искомое решение краевой задачи с использованием выражения (6.17).
140