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

Дифференциальные уравнения. Учебное пособие-1

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
1 Мб
Скачать

б) Если в граничных условиях * * * 0 и внешняя нагрузка отсутствует, т. е. F* (x) 0 (случай свободных колебаний), то для g(x )

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

Гармонические колебания мембраны возможны в случае, когда граничные условия имеют вид

u sin( t ) *(x, y) ,

а внешняя нагрузка задается выражением F sin( t )F*(x, y) . Тогда

прогиб мембраны следует искать в виде

u sin( t ) g(x, y) .

После подстановки u(x, y,t) в уравнение колебаний мембраны получим для g(x, y) неоднородное уравнение Гельмгольца (эллиптического типа)

g

xx

g

yy

gk 2

 

1

F (x, y),

k 2 2

 

.

 

 

 

 

 

 

*

T0

Граничное условие для g(x, y)

T0

будет иметь вид

g *(x, y) .

В случае свободных колебаний F* * 0 . Приведем в этом случае

однородное уравнение Гельмгольца для определения формы колебаний в полярной системе координат:

grr 1r gr r12 g k 2 g 0 .

Для акустических колебаний сжимаемой среды (газа или жидкости) имеем:

g(x, y, z) sin( t ), gxx gyy gzz k 2 g 0, k 2 2 / a2 ,

где a − скорость звука в данной среде. Тогда из граничных условий для

 

sin( t )

(x, y, z) или

sin( t )

*

(x, y, z)

 

*

n

 

 

 

 

 

 

91

получаем граничные условия для g

 

 

 

 

 

g

 

 

(x, y, z) .

g

(x, y, z) или

 

*

 

*

n

 

 

 

 

 

 

6. Применение операционного исчисления к решению уравнений в частных производных

Многие инженерно-технические задачи приводятся к исследованию линейных уравнений в частных производных второго порядка.

Рассмотрим функцию u u(x, t) , зависящую от пространственной

координаты x и времени t ; t изменяется на полупрямой 0 t , а x – на конечном или бесконечном интервале.

Рассмотрим линейное уравнение в частных производных второго порядка с двумя независимыми переменными x и t

A x,t 2u

2B(x,t)

 

2u

C x,t 2u

D x,t u

 

 

 

2

 

x t

2

x

 

x

E x,t u

 

 

t

 

(2.71)

 

G x,t u F(x,t).

 

 

 

t

 

 

F(x,t) 0. Знак дискриминанта

Уравнение называется однородным, если

B2 AC уравнения (2.71) определяет тип уравнения.

Уравнение (2.71) называется в точке М(x,t) уравнением

гиперболического типа, если

в

точке

М дискриминант

уравнения

B2 AC 0 ; эллиптического

типа,

если

в точке М

B2 AC 0 ;

параболического типа, если в точке М B2 AC 0 .

Уравнения в частных производных могут быть решены операционным методом.

Гиперболические и параболические уравнения описывают так называемые нестационарные задачи (переходные процессы физической системы), решение которых существенно зависит от начальных условий. Исследование переходных (нестационарных) процессов физической системы представляет собой технически важную задачу.

Рассмотрим однородное уравнение вида

A x

2u

C x

2u

D x u

E x u

G x u 0,

(2.72)

x2

t2

 

 

x

t

 

 

где A, C, D, E и G – непрерывные функции на интервале 0 x l . При

A 0 уравнение (2.72) будет гиперболического типа, если C 0, и параболического типа, если C 0 .

92

Нестационарная задача, приводящая к уравнению (2.72), в этом случае формулируется следующим образом. Найти решение u(x,t)

дифференциального уравнения

(2.72)

 

для 0 x l

удовлетворяющее начальным условиям

 

 

u(x,0) (x) , u(x,0)

(x) ,

и граничным условиям

t

 

 

 

 

 

u(0,t) f (t) , u(l,t)

u(l,t)

u(l,t) ,

x

t

 

 

и 0 t ,

(2.73)

(2.74)

где , , – действительные числа (если l , то будет одно граничное условие u(0,t) f (t) ).

Найдем решение нестационарной задачи операционным методом.

Пусть

u(x,t) ,

u(x,t)

и

2u(x,t)

, рассматриваемые как функции

t ,

x

 

2

 

 

 

 

 

 

 

 

 

x

 

 

 

являются оригиналами. Это означает, что каждому значению x , заданному на интервале 0 x l , при t 0 соответствует изображение, являющееся функцией от p и x , т. е.

U (x, p) e ptu(x,t)dt

0

или

u(x,t) U (x, p) .

Предполагая, что интегрирование и дифференцирование по x при преобразовании Лапласа перестановочны, получаем

 

u(x,t)

 

 

 

 

 

 

 

 

 

 

u(x,t) dt dU (x, p)

 

 

 

 

e ptu(x,t)dt e pt

 

 

 

 

 

 

x

 

 

x

 

 

 

 

0

0

 

 

x

 

 

dx

 

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

2

 

 

 

2

 

 

2

u(x,t)

 

 

 

 

 

 

u(x,t)

 

 

 

 

 

 

e ptu(x,t)dt e pt

 

dt

d U (x, p)

 

 

x

2

 

 

 

x

2

 

 

x

2

dx

2

 

 

 

 

 

 

 

0

0

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(x,t)

 

 

 

dU (x, p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2u(x,t)

 

 

 

d 2U (x, p) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

dx2

 

 

 

 

 

 

 

По свойству дифференцирования оригинала имеем

u(x,t)

 

pU (x, p) u(x,0)

t

 

 

 

 

93

и

 

 

 

 

2u(x,t)

 

p2U (x, p) pu(x,0)

u(x,0)

.

2

 

t

t

 

 

 

 

Операторным уравнением для уравнения (2.72) будет обыкновенное дифференциальное уравнение

A

d 2U (x, p)

 

2

 

 

u(x,0)

 

 

dx

2

C p U (x, p) pu(x,0)

t

 

 

 

 

 

 

 

 

 

 

 

 

 

DdU (x, p) E pU (x, p) u(x,0) GU (x, p) 0

 

или в силу (2.73)

dx

 

 

 

 

 

 

 

 

 

 

 

 

DdU (x, p) U (x, p) Cp2

Ep G

 

Ad 2U (x, p)

(2.75)

 

dx2

 

dx

 

 

 

 

 

 

 

 

Cp (x) C (x) E (x) 0.

 

 

 

 

 

Граничные условия уравнения (2.75) получаем из (2.74). Имеем

 

 

 

 

dU (x, p)

U (x, p)

 

x 0 F( p),

 

 

 

 

 

(2.76)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( p )U (x, p) (x)

 

 

0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x l

 

 

 

 

 

 

 

 

 

 

 

Интегрируя уравнение (2.75) каким-либо методом при граничных условиях (2.76), получаем изображение U (x, p) . Тогда с помощью теоремы

разложения или, пользуясь таблицами формул преобразования Лапласа, находим решение.

Пример 2.5. Стержень длины l находится в состоянии покоя и один его конец закреплен, а к свободному концу приложена сила Asin t ,

направленная по оси стержня.

Решение. Уравнение колебаний стержня имеет вид

2u

a

2

2u

,

(2.77)

t2

 

x2

 

 

 

 

где u u(x,t) – продольное смещение; a2 – коэффициент, зависящий от

материала стержня. Начальные условия

u(x,0) 0 , u(x,0) 0,

t

граничные условия

u(l,t)

 

A

 

 

u(0,t) 0 ,

 

sin t ,

(2.78)

x

 

 

 

E

 

 

 

94

 

 

где Е – модуль упругости. По закону Гука сила, действующая вдоль

стержня, равна Asin t E u

, поэтому u

 

 

 

 

 

 

 

A

sin t .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Имеем

 

 

 

 

 

x

u(x,t)

 

 

 

 

 

 

 

 

 

x

 

x l

 

 

 

 

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U (x, p) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2u(x,t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2u(x,t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p2U (x, p)

pu(x,0)

u(x,0)

 

или

p

2U (x, p) ,

2

 

 

 

 

 

 

t

 

 

 

 

 

 

2

 

 

 

t

 

 

 

 

 

2u(x,t)

 

 

 

 

 

 

d 2U (x, p) .

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

 

 

 

dx2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Операторное уравнение, соответствующее уравнению (2.77), имеет вид

 

 

 

 

d 2U (x, p)

 

p2

U (x, p) 0 .

 

 

 

 

 

 

 

 

 

 

 

(2.79)

 

 

 

 

a2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из (2.78) получаем граничные условия для уравнения (2.79)

 

 

 

 

 

 

 

 

u(x,t)

 

x 0 0

 

 

 

 

 

U (x, p)

 

x 0 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(x,t)

 

 

 

A

sin t

 

 

 

 

 

 

U (x, p)

 

 

 

 

 

 

 

A

 

 

 

 

 

.

(2.80)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

x l

 

E

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

x l

 

 

E p2 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Интегрируем обычным методом уравнение (2.79), получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

p

x

~

 

 

 

p

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

U (x, p) c1e

 

 

 

 

c2e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U (x, p) c ch

x c

2

sh

 

x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

Пользуясь граничными условиями (2.80), определяем произвольные постоянные

 

 

с 0 и

с

 

Aa

 

 

 

 

 

 

1

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

2

p p2 2 ch

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

E

 

 

p

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

Таким образом, операторное уравнение (2.79) имеет решение

 

 

 

 

 

 

 

 

 

Aa

 

 

 

 

 

sh

p

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U (x, p)

 

 

 

 

 

a

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

E

p p2 2 ch

p

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

Обозначим A( p) sh

p

x и

B( p) p p2

2 ch

p

l .

 

 

 

 

 

 

 

 

a

a

 

 

 

 

 

 

 

 

Тогда

U (x, p)

Aa

A( p)

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E B( p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решая уравнение B( p) 0 и принимая во внимание,

что

ch

l

p cosi

l

p ,

a

a

находим нули функции B( p) . Имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p 0,

p i ,

pk i k

, где k

a

 

 

1

 

(k 1,2,...) .

 

 

 

 

l

k

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

95

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Функция U (x, p) имеет простые полюсы

p 0,

p i ,

pk i k , где k

a

1

 

(k 1,2,...) ,

l

k

2

 

 

 

 

 

 

 

определяющиеся нулями функции B( p) (предполагаем, что условие резонанса отсутствует, т. е. k ).

Оригинал u(x,t) для изображения U (x, p) находим по теореме

разложения,

 

учитывая

B ( p) p2 2 ch

l

 

p 2 p2 ch

l

 

 

p

l

 

 

p p2

2 sh

l

p .

 

a

 

a

 

 

a

a

По формуле (1.40) имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A(0)

 

 

 

A(i )

 

 

 

 

 

 

A( i )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A(i k )

 

 

 

 

 

 

 

 

 

 

 

u(x,t) Aa

 

 

 

ei t

 

 

e i t

 

 

 

 

 

ei k t

 

 

 

 

 

 

 

 

 

B (i )

B ( i )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

B (0)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k 1 B (i k )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A( i k )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i sin

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i sin x

 

 

 

 

 

 

 

 

 

 

 

 

 

e i k t

 

Aa

 

 

 

 

 

 

 

 

 

 

a

 

 

 

ei t

 

 

 

 

 

 

 

 

 

a

 

 

 

e i t

 

.

 

 

 

 

 

 

 

2 2 cos l

 

 

2 2 cos l

 

 

 

 

 

 

 

k 1 B ( i k )

 

 

 

 

 

 

 

 

 

 

 

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i sin

 

k

 

x

 

 

 

 

 

 

 

 

 

 

i k t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i sin

 

 

k

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

e

i k t

 

 

 

 

 

i

l

 

k

2

 

2

sin

l

 

 

 

 

 

 

 

i

l

 

 

 

 

2

 

2

sin

l

 

 

 

 

 

 

 

 

 

 

 

 

 

a

a

k

 

 

 

a

k

 

k

 

 

 

 

 

 

 

 

 

 

k 1

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

k 1

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

Aa

sin

x sin t

 

 

 

 

 

 

1 k sin

k

x

 

 

 

i

 

 

t

 

 

 

 

 

1 k sin

k

x

 

 

 

i

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

e

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

e

 

k

 

.

 

E

 

 

2 cos l

 

 

l

 

 

2

2

 

 

 

 

l

 

 

 

2 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

k

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

k 1 a

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

1 a

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

Таким образом, искомое решение имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(x,t)

Aa

 

sin

 

 

a xsin t

 

2a

 

 

1

 

 

sin

 

 

a

 

 

xsin kt

,

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 cos l

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

k

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где k

 

a

 

 

 

 

1

 

 

(k 1,2,...) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

k

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 2.6. Найти распределение температур u(x,t) в линейном

проводнике тепла при граничных и начальных условиях

u(x,0) 0 , u(0,t) u0 , 0 x , t 0 .

Решение. Уравнение теплопроводности имеет вид

u a2 2u .t x2

Имеем

u(x,t) U (x, p) ,

u pU (x, p) ,

t

2u

 

d 2U (x, p) .

 

x2

 

dx2

 

96

Операторное уравнение имеет вид

 

 

 

 

d 2U (x, p)

 

p

U (x, p) 0 .

(2.81)

dx2

a2

 

 

 

Решая уравнение (2.81), получаем

U (x, p) c1e x a p c2e x a p .

По условию задачи, функции u(x,t) и U (x, p) при x являются ограниченными, поэтому с2 0 . Пользуясь граничным условием для

уравнения (2.81) U (x, p) x 0 up0 , находим произвольную постоянную

с

u0

. Тогда U (x, p)

u0

 

e

x p

.

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

p

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пользуясь равенством

1

e

a p

 

a

 

 

 

 

 

 

 

 

 

Erf

 

 

 

(смотри таблицу изображений

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

2

t

 

Раздел I, п. 4.2), находим оригинал для функции U (x, p) .

Итак, решение данного уравнения имеет вид u(x,t) u0 Erf 2ax t .

Здесь и ниже принято обозначение Erf t

2

 

e t2 dt .

 

 

t

Для закрепления материала предлагаются упражнения для самостоятельной работы с ответами.

Упражнения

Задача 1. Найти решения волнового уравнения

2u(x,t)

 

1 2u(x, t)

 

 

 

x2

a2 t2

 

при следующих граничных и начальных условиях:

1) u(0,t) u(l,t) 0 ,

 

2b

x,

0 x

l

,

 

 

l

 

 

2

 

u(x,0)

 

 

l

 

 

 

 

2b

l

x ,

x

l,

 

l

2

 

 

 

 

 

 

 

97

u(x,0)

0

0 x l .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ответ:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2b x

 

ash a

p

 

 

0 x

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

p

2

ch

pl

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U (x, p)

 

 

 

 

 

 

 

 

2a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2b l x

 

a sh a

 

(l x)

 

 

l

x

l,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

p

 

 

 

 

 

 

 

p2 ch

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8b

 

( 1)n 1

 

 

 

 

(2n

1) x

 

 

 

 

 

 

 

(2n 1) at

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin

 

 

 

 

 

 

cos

 

 

 

 

 

,

0 x

 

 

 

,

 

 

 

2

 

(2n

 

1)2

 

 

 

 

 

l

 

 

 

l

 

2

 

 

 

n 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(x,t)

8b

 

( 1)n 1

 

 

 

 

(2n 1) (l

x)

 

 

(2n 1) at

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos

 

 

 

 

 

 

,

 

 

 

x l.

 

 

 

2

 

(2n

 

1)

2

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

l

 

 

2

 

 

 

 

 

n 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2) u(0,t) u(l,t) 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u( x,0) A sin n x ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(x,0)

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0 x l,

 

n натуральное число .

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ответ:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n at

 

n x

 

 

 

 

 

U ( x, p)

 

 

 

 

 

Ap

 

 

 

 

 

sin

n x

; u(x,t) Acos

sin

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

l

 

 

 

l

 

 

 

 

 

p2

n

2

 

2

a

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3) u(0,t) u(l,t) 0 ,

 

 

u(x,0) 0 ,

 

 

 

u(x,0) B sin n x

0 x

t

l

 

 

Ответ:

 

 

 

U ( x, p)

B

 

sin

p 2 n 2 2 a 2

 

 

l 2

l, n натуральное число.

n x ; u(x,t)

Bl

sin n at sin n x .

n a

l

l

l

4) u(0,t) u(l,t) 0 , u(x,0) Ax(l x)

u(x,0) 0 0

t

,

x l .

98

Ответ:

 

 

 

 

p

x

e

 

p

(l x)

 

1

 

x(l x)

 

 

2

 

 

 

 

 

U (x, p) 2Aa

e

 

a

 

 

a

 

A

 

 

 

 

 

 

 

 

p

 

1

 

 

 

;

 

 

1

e

 

l

p

3

p

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(x,t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

nl x

 

nl x

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A x(l x) a2t2 ( 1)n t

 

a

 

 

t

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

nl x

 

nl x

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( 1)n 1 t

 

a

t

a

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5)

u(0,t)

 

u(l,t)

0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

x

 

n x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(x,0) Acos

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(x,0) 0

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 x l,

 

n натуральное число .

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ответ:

 

 

 

 

 

Ap

 

 

 

 

 

 

 

 

n x

 

 

 

 

 

 

 

 

n at

 

 

 

n x

 

 

 

 

U (x, p)

 

 

 

 

 

 

 

 

cos

; u(x,t)

Acos

cos

.

 

 

 

p2 n

2

 

2

a

2

 

l

l

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(0,t)

 

 

u(l,t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6)

 

0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(x,0) 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(x,0) B cos n x

 

 

0 x l,

n натуральное число .

 

 

 

 

 

t

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ответ:

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

n x

 

 

 

 

Bl

 

 

 

n at

 

 

n x

 

 

U (x, p)

 

 

 

 

 

 

 

 

 

cos

; u(x,t)

 

sin

cos

.

 

p2 n

2

 

2

a

2

 

l

an

 

l

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(0,t)

 

 

u(l,t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7)

 

0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(x,0) 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 x

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(x,0)

 

 

B cos

 

 

 

l

 

,

 

 

0 x

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,

 

 

x l.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

99

Ответ:

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

cos

(2n 1) x

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U (x, p)

 

 

 

 

 

 

 

4B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 p2

 

 

n 1 (2n

3)(2n 1)(2n 1)

 

p2

2a2 (2n 1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l 2

 

(2n 1) at

 

 

 

t

 

 

 

4l

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

(2n 1) x

 

u(x,t) B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos

 

 

 

sin

 

 

.

2

 

a

2

 

 

 

 

 

 

 

 

 

 

l

 

 

l

 

 

 

 

 

 

 

 

 

n 1 (2n 3)(2n 1)(2n 1)

 

 

 

 

 

 

 

 

 

 

 

Задача 2. Найти решение уравнения теплопроводности

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2u(x,t)

2 u(x,t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

при следующих граничных и начальных условиях:

 

 

 

 

 

 

1) u(0,t) u(l,t) 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(x,0) Asin

n x

 

0 x l,

n натуральное число .

 

 

 

 

 

 

 

 

 

 

Ответ:

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

n x

 

 

 

 

 

n2 2

t

 

 

n x

 

 

 

 

 

U (x, p)

 

 

 

 

 

 

 

 

 

 

sin

; u(x,t)

Ae

 

2

l

2

sin

.

 

 

 

 

p2 n

2

 

2

l

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2l2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2) u(0,t) A,

 

 

u(l,t) 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(x,0)

A

 

(l x)

 

0 x l .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ответ:

U (x, p) lpA (l x) ; u(x,t) Al (l x) (t).

3) u(0,t) 0 , u(l,t)

x u(x,0) Asin (2n

Ответ:

0

 

 

1) x

0 x l,

n натуральное число .

2l

 

 

 

 

 

 

A

 

 

(2n 1) x ; u(x,t) Ae

2n 1 2 2

t sin

(2n 1) x .

U (x, p)

 

 

 

 

sin

4 2l 2

p

2

 

2n 1 2

2

2l

 

 

2l

 

 

4 2l2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

100

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]