Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

9060_d71504bbdcb8f95ffb66490c84b51bf8

.pdf
Скачиваний:
56
Добавлен:
04.07.2021
Размер:
8.45 Mб
Скачать

Пример: Разложить в ряд Фурье функцию с периодом 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

Если функция 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).

 

 

 

 

 

 

 

 

Î

x

x + Dx

l

õ

Ðèñ. 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

 

 

 

 

 

 

 

 

 

 

 

 

"

Соседние файлы в предмете Высшая математика