Пример: Разложить в ряд Фурье функцию с периодом 2l = 2, заданную формулой f(x) = x, x О [-1, 1].
Функция нечетная, поэтому в (31.6) коэффициенты a0 = an = 0,
b |
= 21 x sin npx dx = - |
2x cosnpx |
|
|
1 |
+ |
2 |
1 cosnpx dx = - |
2cosnp |
= |
2( -1)n+1 |
. |
|
|
|
|
|
|
|
n |
ò |
pn |
|
|
|
pn ò |
|
|
|
|
|
np |
|
np |
|
|
0 |
|
|
|
|
|
|
|
|
0 |
|
|
|
0 |
|
¥ |
|
n+1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
( -1) |
|
|
|
|
|
|
В точках непрерывности |
f (x) = |
å |
sin npx, |
à â òî÷- |
|
|
p |
|
|
|
|
|
|
|
|
|
|
|
n=1 |
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ках разрыва х = ±1, ±2, ... имеем S(x) = 0 ![](/html/67644/48/html_BT7NftrvLk.x1Dx/htmlconvd-EQyMuo391xi3.jpg)
Если функция f (x) задана на [0, l], то для ее разложения в ряд Фурье необходимо продолжить сначала по некоторому закону f(x) на [-l, 0] четным или нечетным образом. Например, функцию f(x) = x, x О [0, p] можно продолжить на [-p, 0] нечетным образом (см. рис. 31.1) и по-
лучить разложение в ряд по синусам (31.4) или четным образом (см.
рис. 31.2) и получить разложение в ряд по косинусам (31.5). Можно
осуществить продолжение и другими способами, т.е. получить бес-
конечное множество рядов Фурье, однако на [0, l] все они имеют сум-
ìó f (x).
На практике часто пользуются так называемыми усеченными рядами Фурье:
a0 + ån ak coskx + bk sinkx,
2k=1
где n — конечное число. Сумма такого усеченного ряда лишь приближенно аппроксимирует исходную функцию f(x). При этом ошибка d(x) = |S(x) — f(x)| может быть оценена методами, изложенными в [2. С. 244].
Ряды Фурье широко используются в приложениях, например, при решении линейных дифференциальных уравнений в частных произ-
водных [см. разд. 33], а также в практическом гармоническом анализе [2. С. 245] для представления функций, заданных таблично или графически, в виде аналитических выражений.
Литература: [2. С. 229–254]; [6. С. 410–413]; [10. С. 178– 191];
[11. Ñ. 328–373].
Глава 11
УРАВНЕНИЯ МАТЕМАТИЧЕСКОЙ ФИЗИКИ
32. ОСНОВНЫЕ ТИПЫ УРАВНЕНИЙ МАТЕМАТИЧЕСКОЙ ФИЗИКИ
Опорный конспект ¹ 32
32.1. Понятие об основных УМФ |
|
1. Волновое уравнение |
|
|
¶2u |
= a2 |
¶2u |
, a = const |
(32.1) |
|
¶t2 |
|
|
|
|
|
¶x2 |
|
2. |
Уравнение теплопроводности |
(уравнение Фурье) |
¶u = a |
2 |
¶2u |
, a = const. |
(32.2) |
2 |
¶t |
|
¶x |
|
3. Уравнение Лапласа |
¶2u |
+ |
¶2u |
= 0. |
(32.3) |
¶x2 |
¶y2 |
Задаются начальные и граничные условия, обеспечивающие корректность задач
32.2.Классификация линейных дифференциальных уравнений
âчастных производных II порядка
a11 |
¶2u |
+ 2a12 |
¶2u |
+ a22 |
¶2u |
+ b1 |
¶u |
+ b2 |
¶u |
+ cu=F (x, y),(*) |
¶x2 |
¶x¶y |
¶y2 |
|
|
|
|
|
|
¶x |
¶y |
a11,a12,a22,b1,b2,c — const
a11 (dy)2 + 2a12dy dx + a22(dx)2 = 0 — характеристическое уравнение Ю
dy |
= |
-a12 ± |
a122 - a11a22 |
= a1,2 |
Þ y = a |
|
x + c |
|
|
|
1,2 |
dx |
a11 |
|
1,2 |
|
|
|
1) D = a212 - a11a22>0 — гиперболический тип, приводится к виду
|
¶2u |
= Ô(x, h, u, |
¶u |
, |
¶u |
) |
èëè |
¶2u |
- |
¶2u |
= Ô; |
|
|
|
|
|
|
|
¶x¶h |
|
¶x |
|
¶h |
|
¶x2 |
|
¶h2 |
|
2) D = 0 — параболический тип |
¶2u = Ô; |
|
|
|
|
|
|
|
|
|
¶x2 |
|
3) D < 0 — эллиптический тип ¶2u + ¶2u = Ô ¶x2 ¶h2
УМФ (32.1) — 1-й тип, (32.2) — 2-й тип, (32.3) — 3-й тип. Три вида задач для уравнений:
1)задача Коши для 1, 2-го типов;
2)краевая задача для 2-ãî òèïà;
3)смешанная задача для 1, 2-го типов (задаются начальные
èграничные условия)
32.1. Понятие об уравнениях математической физики. Граничные и начальные условия
Любой технологический процесс протекает в условиях действия основных законов природы, некоторые из них достаточно хорошо изучены и сформулированы в виде математических зависимостей. Часто эти зависимости представляют в форме дифференциальных уравнений в частных производных.
Мы рассмотрим так называемые основные уравнения математической физики (УМФ) [2. С. 294]. С целью упрощения изложения будем рассматривать случай линейных дифференциальных
уравнений 2-го порядка с частными производными для функций двух переменных.
Основными УМФ называют:
1. Волновое уравнение для неизвестной функции u(x,t)
¶2u2 |
= a2 |
¶2u |
, |
a = const |
(32.1) |
2 |
¶t |
|
¶x |
|
|
Таким уравнением описывают процессы поперечных колебаний струны, продольных колебаний стержня, электрических колебаний, колебаний газа и т.д.
2. Уравнение теплопроводности, или уравнение Фурье
¶u = a2 |
¶2u |
, |
a = const, |
(32.2) |
2 |
¶t |
¶x |
|
|
с помощью которого описывают распространение тепла, фильтрацию жидкости и газа в пористой среде и т.д.
3. Уравнение Лапласа для функции u(x,y)
|
¶2u |
+ |
¶2u |
= 0, |
(32.3) |
|
¶x2 |
¶y2 |
|
|
|
|
к которому приводят стационарные задачи об электрических и магнитных полях, гидромеханики, диффузии и т.д.
Соответствующие уравнения для функций с бóльшим числом переменных, например для функций u(x,y,t):
|
¶2u |
= a |
2 |
æ |
¶2u |
+ |
¶2u |
ö |
|
|
|
|
|
ç |
|
|
|
|
|
|
|
|
|
|
÷; |
|
|
2 |
|
|
|
|
|
2 |
|
|
2 |
|
|
¶t |
|
|
|
|
|
ç |
¶x |
|
|
|
|
|
¶y |
|
|
|
÷ |
|
|
|
|
|
|
è |
|
|
|
|
|
|
|
|
ø |
|
¶u |
= a |
2 |
æ |
|
¶2u |
|
+ |
¶2u ö |
|
|
|
|
|
ç |
|
|
|
|
|
|
|
|
|
|
÷; |
|
|
|
|
|
|
|
2 |
|
|
|
|
2 |
|
|
¶t |
|
|
|
|
ç |
|
¶x |
|
|
|
|
|
¶y |
|
÷ |
|
|
|
|
|
|
è |
|
|
|
|
|
|
|
ø |
|
|
¶2u |
+ |
¶2u |
+ |
¶2u |
= 0. |
|
¶x2 |
¶y2 |
¶z2 |
|
|
|
|
|
|
|
|
|
|
Для описания реального физического процесса необходимо знать начальное состояние процесса (начальные условия), а также значения функции и ее производных на границе области (граничные условия).
Так, для уравнения (32.1) необходимы четыре условия: по два на каждую из переменных t и x. Например, задача об определении колеблющейся струны: найти дважды дифференцируемую в области D = {(x, t): 0 £ x £ l, t ³ 0} функцию u(x,t), которая удовлетворяет уравнению (32.1) при заданных условиях:
— начальные lim u(x,t) = j(x), |
lim |
¶u = y(x), |
0 £ x £ l, |
t ®+0 |
t ®+0 |
¶t |
|
— граничные lim u(x,t) = m1(t), |
lim u(x,t) = m2(t), t ³ 0. |
x®+0 |
x®l -0 |
Функция j(х) дает отклонение струны от положения равновесия в т. х в момент t = 0, функция y(х) — скорость струны в т. х в момент t = 0, функции m1(t), m2(t) — законы движения концов струны.
Для уравнения (32.2) начальные и граничные условия записываются в виде
lim u(x,t) = u0(x), 0 £ x £ l; |
t ®+0 |
|
lim u(x,t) =T1(t), |
lim u(x,t) =T2(t), t ³ 0. |
x®+0 |
x®l-0 |
Для уравнения (32.3) задаются только граничные условия, а задача называется краевой. Если ищется функция u(x,y), дважды дифференцируемая и удовлетворяющая в области D уравнению (32.3), принимающая в каждой т. M О L = ¶D заданное значение u|L = y(M), то такая краевая задача называется задачей Дирихле.
Граничные и начальные условия обеспечивают единственность решения поставленных задач [7. С. 46, 196, 297]. Отметим, что с практической точки зрения задача должна быть поставлена корректно, т.е. ее решение не только должно существовать, но и быть единственным и устойчивым (малым изменениям исходных данных должно соответствовать малое изменение решения) [7. С. 297]; [2. С. 325].
32.2. Классификация линейных дифференциальных уравнений в частных производных 2-го порядка
Было отмечено, что для каждого из приведенных УМФ начальные и граничные условия ставятся по разному. Для того чтобы корректно в указанном выше смысле поставить задачу, прежде
следует классифицировать УМФ так, чтобы внутри каждого класса уравнений постановка этих условий была одинаковой.
Классификация УМФ в общем случае представляет собой достаточно сложную задачу [7. С. 18]. Мы же ограничимся лишь случаем линейного дифференциального уравнения в частных производных 2-го порядка с постоянными коэффициентами, общий вид которого
¶2u |
¶2u |
|
¶2u |
|
¶u |
|
¶u |
|
(32.4) |
a11 ¶x2 + 2a12 |
¶x ¶y + a22 |
¶y2 + b1 ¶x + b2 |
¶y + cu = F (x, y), |
|
ãäå a11, a12, a22,b1, b2, c — постоянные; u(x,y), F(x,y) — неизвестная и заданная функции соответственно.
Путем замены в (32.4) переменных х, y на новые x, h можно добиться упрощения этого уравнения [7. С. 12]. Для этого предварительно составляется уравнение [2. С. 307] a11(dy)2 + 2a12dy dx + + a22(dx)2 = 0, которое называется характеристическим. Из него находят
dy |
= |
-a12 ± |
a122 - a11a22 |
= a1,2 Þ y = a1,2x + c1,2, c1,2 = ñînst. |
|
|
|
dx |
a11 |
Возможны три случая:
1)a212 - a11a22>0, уравнение (32.4) может быть путем замены x = y - a1x, h = y - a2x приведено к виду
|
|
|
|
|
|
|
¶2u |
|
+ b1* ¶u |
+ b2* ¶u |
+ c*u = f , |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
¶x ¶h |
¶x |
¶h |
|
¶2u |
|
¶2u |
|
* ¶u |
* ¶u |
|
|
|
|
èëè ¶x2 |
- |
|
+ b1 ¶x |
+ b2 ¶h |
+ c*u = f |
— гиперболический тип |
|
¶h2 |
|
уравнения; |
|
|
|
|
|
|
|
|
|
|
2) a2 |
- a a |
= 0, уравнение (32.4) приводится к виду |
|
12 |
11 |
22 |
|
|
|
|
|
|
|
|
|
¶2u |
* ¶u |
* ¶u |
|
— параболический тип уравнения; |
|
¶x2 + b1 ¶x + b2 |
¶h = f |
|
|
|
|
|
|
3) a212 - a11a22<0, уравнение (32.4) приводится к виду
|
¶2u |
+ |
¶2u |
+ b1* |
¶u |
+ b2* |
¶u |
+ c*u = f |
— эллиптический тип уравне- |
|
¶x2 |
¶h2 |
¶x |
¶h |
|
|
|
|
|
|
íèÿ.
Таким образом, (32.1), (32.2), (32.3) — уравнения гиперболического, параболического и эллиптического типов соответственно.
Различают три вида задач для этих уравнений:
1)задача Коши для уравнений гиперболического и параболи- ческого типов: задаются начальные условия, граничные отсутствуют, область определения уравнения и его решения не ограни- чена;
2)краевая граничная задача для уравнений эллиптического типа: задаются граничные условия на границе L = ¶D области определения неизвестной функции, начальные условия отсутствуют;
3)смешанная задача для уравнений гиперболического и параболического типов (задаются начальные и граничные условия).
32.3.Построение математической модели задачи
îраспространении тепла
Рассмотрим однородный неравномерно нагретый стержень длиной l. Условимся о предположениях, что его боковая поверхность теплонепроницаема и во всех точках поперечного сечения температура одинакова. Расположим ось Ох так, чтобы один конец стержня находился в точке О, а другой — в точке x = l (рис.32.1). Обозначим через u(x,t) температуру стержня в момент t в точке х. Требуется найти зависимость u(x,t), если при t = 0 функция u(x,0) = j(х) (в начальный момент задана температура в различных сечениях стержня) и при х = 0, х = l соответственно u(0,t) = y1(t), u(l,t) = y2(t) (известна температура на концах стержня в любой момент t).
Ðèñ. 32.1
Известно, что скорость распространения тепла q, т.е. количе- ство тепла, протекающего через сечение с абсциссой х за единицу времени, находится по формуле
q = -g ¶u S,
¶x
где S — площадь сечения, g — коэффициент теплопроводности. Рассмотрим элемент стержня между сечениями с абсциссами х и х + Dх. Количество тепла DQ, сообщенное ему за время Dt, определяется формулой
|
æ ¶u(x + Dx,t) |
|
¶u(x,t) ö |
|
DQ =Q(x + Dx) -Q(x) = gS ç |
|
- |
|
÷ Dt. |
|
¶x |
¶x |
|
è |
|
ø |
Применяя формулу Лагранжа (разд.10.1), получим
DQ = gS |
¶2u(x + Q1Dx, t) |
DtDx, 0 < Q1< l. |
|
|
¶x2 |
Пусть выделенный элемент стержня является столь малым, что в каждый момент температуру всех его точек можно считать одинаковой. Тогда сообщенное этому элементу стержня количество тепла DQ за время Dt можно определить по формуле
DQ = ÑrS(u(x, t + Dt) - u(x,t))Dx,
где С и r — соответственно теплоемкость и плотность вещества, из которого состоит стержень.
Применяя формулу Лагранжа, будем иметь
|
DQ = CrS |
¶u(x,t + Q2Dt) |
DtDx, 0 < Q2< l. |
|
|
|
|
|
|
|
|
¶x |
|
|
|
|
Приравнивая выражения для DQ, получаем |
|
C rS |
¶u(x,t + Q2 Dt) |
Dt Dx = |
gS |
¶2u(x + Q1Dx,t) |
Dt Dx. |
|
¶x2 |
|
¶x |
|
|
|
После сокращения на SDtDх и перехода к пределу при Dх ® 0, Dt ® 0 окончательно получаем уравнение
¶u = a2 ¶2u |
, |
a2 = |
g |
. |
|
¶t |
¶x2 |
|
Cr |
Это и будет уравнение распространения тепла в однородном стержне (уравнение теплопроводности) (32.2).
Итак, математическая модель задачи о распространении тепла в стержне имеет следующий вид (смешанная краевая задача):
Найти функцию u(x,t), удовлетворяющую уравнению теплопроводности,
начальному условию u(x,0) = j(x), 0 £ x £ l,
и граничным условиям u(0,t) = y1(t), u(l,t) = y2(t), 0 £ t £ T. (32.4).
Литература: [2. С. 292–306]; [9. С. 83–114].
33.МЕТОДЫ РЕШЕНИЯ УРАВНЕНИЙ МАТЕМАТИЧЕСКОЙ ФИЗИКИ
Опорный конспект ¹ 33
33.1. Mетод Даламбера
Состоит в упрощении уравнения (*) гиперболического типа (см. ОК 32.2) путем замены
x = y - a1x, h = y - a2x. Пример — решение задачи Коши:
ì |
¶2u |
|
2 ¶2u |
|
|
|
|
|
ü |
ï |
|
|
= a |
|
|
|
, |
t > |
0, x Î R, ï |
¶t |
2 |
|
¶x |
2 |
ï |
|
|
|
|
|
|
|
|
|
ï |
í |
|
|
|
|
|
|
|
|
¶u |
|
|
ý Þ |
ïu(x,0) = j(x), |
|
|
|
= y(x),ï |
|
|
|
|
ï |
|
|
|
|
|
|
|
|
¶t |
|
|
ï |
î |
|
|
|
|
|
|
|
|
|
t =0 |
þ |
|
|
|
|
|
|
|
|
1 x+at
u(x,t) = [j(x - at) + j(x + at)]/2 + 2a x-òat y(x)dx
33.2. Метод Фурье
Он основан на разделении переменных в (*) путем замены u(x,y) = X(x)Y(y).
33.2.1. Решение смешанной задачи:
ì |
|
¶2u |
= a |
2 ¶2u |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ü |
|
|
|
|
|
|
|
|
|
ï |
|
|
|
|
|
|
|
|
|
|
|
|
, |
t > 0, 0< x < l, |
|
|
|
|
|
ï |
|
|
|
|
|
|
|
|
|
¶t |
2 |
|
|
|
|
¶x |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ï |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ï |
|
|
|
|
|
|
|
|
|
ï |
|
|
|
|
|
|
= u(l,t) = 0, t ³ 0; |
|
|
|
|
|
|
ï |
Þ |
|
|
|
|
|
|
|
|
íu(x,0) |
|
|
|
|
|
|
|
ý |
|
|
|
|
|
|
|
|
ï |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
¶u |
|
|
|
|
|
|
|
|
|
|
ï |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ïu(x,0) |
= f (x), |
|
|
|
= j(x), 0 |
£ x £ lï |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ï |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
¶t |
t =0 |
|
|
|
|
|
|
ï |
|
|
|
|
|
|
|
|
|
î |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
þ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
¥ |
æ |
|
|
|
|
|
anp |
|
|
|
|
anp |
ö |
np |
|
|
|
|
|
|
|
|
u(x,t) = å |
çCn cos |
|
|
t |
+ Dn |
sin |
|
|
|
|
t |
÷ sin |
|
|
x, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n=1 |
è |
|
|
|
|
|
|
l |
|
|
|
|
l |
|
|
|
ø |
l |
|
|
|
|
|
Þ |
|
|
|
|
2 l |
|
|
|
|
|
|
|
np |
|
|
|
|
|
|
|
|
l j(x)sin |
np |
|
|
|
|
|
|
|
C |
|
|
= |
f (x)sin |
x dx, |
D |
= |
|
2 |
|
x dx |
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
l ò |
|
|
|
|
|
|
|
|
|
l |
|
|
|
n |
|
anp |
ò |
|
|
|
|
l |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
33.2.2. Решение смешанной задачи: |
|
|
|
|
|
|
|
|
|
|
|
|
ì |
|
¶u |
= a |
2 ¶2u |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
¥ |
|
|
æ anp ö2 |
|
|
ï |
|
|
|
|
|
|
|
|
|
|
, t > 0, 0< x < l; |
|
|
|
|
|
|
|
|
- ç |
|
|
÷ t |
np |
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
l |
ï |
|
¶t |
|
|
|
|
|
|
|
¶x |
|
|
|
|
|
|
|
|
|
|
|
|
|
u(x,t) = åBne è |
|
|
ø sin |
|
x, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ï |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n=1 |
|
|
|
|
|
|
l |
íu(0,t) = u(l,t) = 0, t ³ 0, |
Þ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ï |
|
|
|
|
|
|
= y(x), 0 £ x £ l, |
|
|
|
|
|
|
2 l |
|
|
|
np |
|
|
|
ïu(x,0) |
|
Bn = |
|
|
|
y(x)sin |
|
|
x dx |
|
|
|
|
|
|
|
|
|
|
ï |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
l ò0 |
|
|
|
|
l |
|
|
|
î |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
33.2.3. Решение задачи Дирихле в круге: |
|
|
|
|
|
|
|
|
|
¶2u |
|
|
¶ |
2u |
|
|
|
|
|
2 |
2 |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+ |
|
|
|
|
|
|
|
= |
0; |
D: x |
+ y < R ; |
u|¶ |
|
= f(x,y). |
|
|
|
|
|
|
|
|
|
¶x2 |
|
¶y2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
D |
|
|
|
|
|
|
|
|
|
|
|
|