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

Tom_2

.pdf
Скачиваний:
59
Добавлен:
18.02.2016
Размер:
3.2 Mб
Скачать

при свободном внутреннем теплообмене, если в левом конце его (при x = 0 ) поддерживается постоянная температура u0 , а через

правый конец (при x = l > 0 ) происходит теплообмен с окружающей средой, температура которой задана функцией ϕ(t) .

(Указание. Для получения граничного условия при x = l воспользоваться тем, что количество тепла, проходящего через

конец стержня

за время t , равно q = −α

u

 

S

t , где

S

 

x

 

 

 

x=l

 

 

 

 

 

 

 

площадь поперечного сечения стержня, а α

 

– коэффициент

теплопроводности материала стержня. С другой стороны,

q

пропорционально произведению разности температур

u(l,t) −ϕ(t)

на

площадь

 

 

 

S

и на время Dt ,

т.е. q = h[u(l,t) −ϕ(t)]S t , где

 

h ³ 0

зависит от

качества теплоизоляции конца стержня от окружающей среды. Сравнивая оба выражения для q , получаем граничное условие

u(l,t) + h [u(l,t) −ϕ(t)] = 0 ).

x α

7. Используя телеграфные уравнения, поставить задачу об отыскании

закона изменения

напряжения и силы тока в длинной

линии

(0 ≤ x < ∞)

без потерь (т.е. R = G = 0 ), если известны начальные

напряжение

ϕ(x) ,

сила тока ψ (x) , а напряжение в точке

x = 0

постоянно и равно E0 . (Указание. Найти функции u(x,t) и i(x,t) , определенные при 0 £ x < +¥ и 0 £ t < ¥ , являющиеся решениями

уравнения

2w(x,t)

= LC

2w(x,t)

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

x2

t2

 

 

 

условиям

u(x,0) = ϕ(x) ,

i(x,0) =ψ (x) , и граничным u(0,t) = E0 ,

u(+∞,t) = i(+∞,t) = 0 ).

 

 

8.На границе бесконечного изотропного однородного цилиндра, направляющая которого – кривая L – лежит в плоскости,

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

304

распределении температуры внутри цилиндра (плоская задача Дирихле). (Указание. Систему координат выбираем так, что ось Oz параллельна образующей. Область, ограниченную кривой L в плоскости Oxy , обозначим через D. Нужно найти функцию u(x, y) , которая для всех точек M (x; y) D удовлетворяет

уравнению 2u + 2u = 0 , т.е. является гармонической в области

x2 y2

D, а на границе области, т.е. на кривой L, удовлетворяет условию

u(x, y )

 

(x;y) L = f (x, y) , где

f (x, y)

– заданная на L непрерывная

 

 

функция. Для установившегося теплового процесса

u

= 0 ).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

9. Найти решения уравнения

2u

=

2u

при заданных начальных

t2

x2

условиях

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

u (x,0) = sin x ;

а) u(x,0) = sin x

,

 

u (x,0) = 0 ; б) u(x,0) =

 

,

 

 

+ x2

 

 

 

 

x

 

 

t

 

 

 

1

 

t

в) u(x,0) =

 

1

 

,

u (x,0) = cos x .

 

 

 

 

 

 

 

 

 

+ x2

 

 

 

 

 

 

 

 

1

 

t

 

 

 

 

 

 

 

 

 

 

10.Струна, закрепленная на концах x = 0 и x = 1 в начальный момент имеет форму u = h(x4 - 2x3 + x). Найти форму струны для любого момента времени t, если начальные скорости отсутствуют.

11.Струна закреплена в точках x = 0 и x = l . Начальные отклонения точек струны равны нулю, а начальная скорость выражается формулой

 

 

 

ì

æ

l ö

 

 

 

 

 

 

 

 

 

 

 

 

 

ï

π ç x -

 

÷

 

 

 

l

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

ïcos

è

2 ø

при

x -

 

<

 

,

 

 

 

 

 

 

 

 

 

t

 

t=0

= í

h

 

 

 

2

 

2

 

 

ï

 

 

 

 

 

 

l

 

 

 

h

 

 

 

 

 

ï

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при

x - 2

 

> 2 .

 

 

 

î0

 

 

 

 

Найти форму струны для любого момента времени t.

305

12. Найти решение

уравнения

 

 

 

u

=

2u

,

 

удовлетворяющее

 

 

 

t

x2

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ì

-

 

 

x

 

 

 

 

 

£ x £ l,

 

 

 

 

 

 

 

 

 

 

 

 

 

ï1

 

 

 

 

при 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

ï

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

í

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u(x,t )

 

t=0 =

f (x) = ï1

 

+

 

 

при

 

- l £ x £ 0,

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ï

 

 

 

 

 

при x ³ l и x £ -l.

 

 

 

 

 

 

 

 

 

 

î0

 

 

 

 

 

(Указание. Решение выражается формулой

 

 

 

 

 

 

 

 

1

 

0

æ

τ ö

(x−τ )2

 

 

 

 

 

 

1

 

 

l æ

 

 

τ ö

(x−τ )2

 

u(x,t) =

 

 

 

 

ò

ç1+

÷e

 

4t

 

 

dτ +

 

 

 

 

 

òç1

-

÷e

 

4t dτ .

 

 

 

 

 

 

 

 

 

 

 

 

π t

 

 

 

 

 

 

2

l

è

l ø

 

 

 

 

 

 

2

πt 0 è

 

 

l ø

 

 

 

Заменой x2-τt упростить ответ).

13. Дан тонкий однородный стержень длины l, изолированный от внешнего пространства, начальная температура которого равна

f (x) = cx(l - x) . Концы стержня поддерживаются при температуре l2

равной нулю. Определить температуру стержня в момент времени t > 0 . (Указание. Закон распределения температуры стержня

описывается уравнением

u

= a

2

2u

, начальным условием

t

 

x2

 

 

 

cx(l - x)

 

 

 

 

 

 

 

 

 

u

 

t=0 = f (x) =

и краевыми условиями u

 

x=0 = u

 

x=l

= 0 ).

 

 

 

 

l2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

14. Решить задачу Дирихле для уравнения Лапласа u = 0

в круге

0 ≤ ρ ≤ 1, если u ρ=1 запишется в виде:

u(ρ,ϕ) = 43 π 2

= ϕ2 = f%(ϕ) .

æ

1

+ 4åρn

ç

 

 

 

2

n=1

è n

 

Указание. Искомое решение

cos nϕ - 4nπ sin nϕ ö÷ .

ø

15.При условии h = 1 в прямоугольнике Q с вершинами в точках A(0;0), B(0;5), C(4;5), D(4;0) и границей S найти приближенное

решение граничной задачи

306

uxx ut = 0, (x;t) Q, u(x,t) = ϕ(x,t), (x;t) S,

если ϕ(x,0) = x, ϕ(0,t) = 0, ϕ(4,t) = 4 .

307

16. Найти для x = 0,1× m ,

m = 0,1,...,10,

t = 0,02 решение уравнения

u

= 2u

+α (x2 - 2t ), 0 £ x £ 1,

0 £ t £ 0,02,

 

t

x2

 

 

 

u(0,t) = 0 ,

 

u(1,t) = αt ,

0 £ t £ 0,02 ,

удовлетворяющее

условиям

 

u(x,0) = 0 ,

0 £ x £ 1 ,

α = 0,5 × k ,

k = 1,2,3,4,5,6 взяв h = 0,1,

τ = 0,005 .

 

 

 

 

 

 

 

 

 

 

 

17. При условии h = 1 найти

приближенное решение задачи Гурса

uxy = 0 ,

0 < x < ¥ ,

0 < y < ∞ ,

 

u(0, y) = ϕ(y) ,

0 < y < ∞ ;

u(x,0) =ψ (x) , 0 < x < ¥ ,

в узлах

(2;2) ,

(2;3) , (2;4) . Отдельно

рассмотреть случаи, если:

 

 

 

 

 

 

 

 

а) ϕ(y) = 0, ψ (x) = x ; б) ϕ( y) = y, ψ (x) = 0 ;

 

 

в) ϕ( y) = y, ψ (x) = x .

 

 

 

 

 

 

 

 

18. Найти решение уравнения

2u(x, y)

+

2u(x, y)

= 0 ,

0 £ x £ 0,8 ,

x2

 

 

y2

 

 

 

 

 

 

 

 

 

 

0 ≤ y ≤ 0,8 , удовлетворяющее граничным условиям

 

u(x,0) = x2 ,

u(x;0,8) = x2 + 0,64,

0 £ x £ 0,8,

 

u(0, y) = y2 ,

u(0,8; y) = y2 + 0,64,

0 £ y £ 0,8

 

на сетке (xm ; yn ),

xm = 0,2 × m,

yn = 0, 2 × n .

 

 

 

308

309

ГЛАВА 9

ОПЕРАЦИОННОЕ ИСЧИСЛЕНИЕ

§1. Преобразование Лапласа. Основные понятия

исвойства

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

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

1.От искомых функций переходят к некоторым другим функциям – их изображениям.

2.Над изображениями производят операции, соответствующие заданным операциям над самими функциями.

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

В качестве преобразования, позволяющего перейти от функций

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

10. Оригинал и изображение. Оригиналом будем называть комплекснозначную функцию f (t) действительного переменного t,

удовлетворяющую следующим условиям:

1)f (t) = 0 , если t < 0 ;

2)f (t) – кусочно непрерывная, интегрируемая функция на любом

конечном интервале оси Ot ;

3) с возрастанием t модуль функции f (t) растет не быстрее некоторой показательной функции, т.е. существуют числа M > 0 и s0 ³ 0 , такие, что для всех t имеем

f (t) £ Mes0t .

Нижнюю грань множества чисел s0 , удовлетворяющих указанному неравенству, называют показателем роста функции f (t) .

310

Изображением функции f (t) (по Лапласу) называется функция

F( p) комплексного переменного p = s + iδ , определяемая равенством

¥

 

F( p) = ò f (t)e- pt dt,

(1)

0

где интеграл берется по положительной полуоси.

Преобразование (1), ставящее в соответствие оригиналу f (t) его изображение F( p) , называется преобразованием Лапласа. Соответствие

между оригиналом

f (t) и изображением F( p)

записывается в виде

.

.

 

f (t) = F( p) или F( p)= f (t) .

 

×

×

 

Для всякого оригинала f (t) изображение

F( p) определено в

полуплоскости Re p > s0 , где s0 – показатель роста f (t) , и является в

этой полуплоскости аналитической функцией.

Если функция F( p) является изображением функции f (t) , то

lim F( p) = 0 .

p®¥

В самом деле,

lim F( p) = lim

æ¥

f (t)e

- pt

 

ö

=

lim

æ

¥

e

-( p-s )t

ö

=

ç

ò

 

dt ÷

ç M

ò

0

dt ÷

p®¥

p®¥

ç

 

 

 

 

÷

 

p®¥

ç

 

 

÷

 

 

 

è

0

 

 

 

 

 

ø

 

 

è

0

 

 

ø

 

 

 

 

æ

 

 

 

1

 

e-( p

-s0 )t

 

 

¥ ö

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= lim

ç

-M

 

 

 

 

÷

= 0.

 

 

 

p - s

 

 

 

 

 

p®¥

ç

 

 

 

 

 

 

 

÷

 

 

 

 

 

 

 

 

è

 

 

 

 

0

 

 

 

 

 

0 ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 1. Установить, является ли оригиналом функция

0

при t < 0;

 

 

 

 

 

 

 

0

при t < 0;

 

 

а) f (t) = {2

при t > 0;

 

 

 

б)

f (t) = {tgt

при t > 0.

 

 

Решение. а) Данная функция является оригиналом, так как она непрерывна для всех t > 0 и является функцией ограниченного роста,

т.е. может быть представлена при t > 0 в виде 2 = 2 ×e0×t . Отсюда следует, что s0 = 0 , M – любое действительное число, большее или

равное 2 (см. свойство 3)).

б) Функция в этом случае не является оригиналом, т.к. в точках t = π2 + π k , k = 0,1,2,..., которые принадлежат промежутку (0;∞) , она

имеет разрывы второго рода. Пример 2. Пользуясь определением, найти изображения функций

а) функции Хевисайда

ì1,

t > 0,

f (t) = í

t < 0;

 

î0,

311

б)

f (t) = eαt , α – любое число;

 

 

 

 

 

в)

f (t) = tα , α > -1.

 

Решение. а) По формуле (1)

при Re p > 0 (s0 = 0) находим

 

 

 

 

¥

 

- pt

 

b

 

- pt

 

 

 

æ

 

1

 

- pt ö

 

b

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F( p) =

ò

1×e

 

dt = lim

ò

e

 

dt = lim

ç

-

 

e

÷

 

=

 

,

 

 

 

p

p

 

 

 

 

 

 

b®¥

 

 

 

b®¥

è

 

 

ø

 

 

 

 

F( p) =

1

 

0

 

 

 

0

 

 

 

. 1

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т.е.

 

, или, в символьной записи, 1=

 

.

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

×

 

p

 

 

 

 

 

 

 

 

б) Данная функция является оригиналом. По формуле (1) имеем

¥

 

 

 

 

b

 

 

 

 

 

 

 

1

 

F ( p) = òeαt e- pt dt = lim

òe-( p-α )t dt = - lim

 

e-( p-α )t

 

-α

0

b®¥

0

 

 

 

 

b®¥ p

 

æ

 

 

 

 

 

e-( p-α )b ö

 

 

 

 

 

 

 

 

 

1

 

 

 

 

1

 

 

 

= lim

ç

 

 

 

 

-

 

÷

=

 

 

 

,

 

 

 

 

 

 

 

 

 

 

b®¥

ç

 

 

 

 

 

 

p -α

÷

 

p -α

 

 

 

è p -α

 

 

ø

 

 

 

 

если Re( p −α) > 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, eαt

.

 

1

 

 

 

 

 

 

 

 

 

 

=

 

 

 

.

 

 

 

 

 

 

 

p

-α

 

 

 

 

 

 

 

 

 

×

 

 

 

 

 

 

 

 

в) По определению изображения имеем

F( p) = ò tα e- pt dt .

b

=

0

0

Полагая pt = τ , получим F( p) =

1

òe-ττα dτ , где γ – луч с

α +1

 

p

γ

 

 

направлением arg p < π2 (рис. 1).

v

v

 

γ

γ -

 

 

 

γ R

 

 

t

γ1

t

Рис. 1

 

Рис. 2

Функция τ ® e-ττα аналитическая в полуплоскости T = {τ Î | Reτ > 0} . Рассмотрим замкнутый контур L, составленный

из упорядоченного набора (γ1R - ) ориентированных частей (см. рис. 2). По теореме Коши для аналитических функций имеем

312

òe-ττα dτ = ò e-ττα dτ + ò e-ττα dτ + ò e-ττα dτ = 0,

 

 

L

 

 

 

γ1

γ R

γ

 

откуда

 

òe-ττα dτ = ò e-ttα dt + ò e-ττα dτ .

 

 

 

 

 

 

 

 

γ

 

 

γ1

γ R

 

 

 

Учтем,

что

òe-ττα dτ = - ò e-ττα dτ и

τ = t на γ1 , а также

 

 

 

 

γ

 

 

γ

 

 

 

lim e-ττα = 0 ,

lim

 

ò e-ττα dτ = 0 . Следовательно,

 

τ

®+¥

R®+¥

γ R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

òe-ττα dτ = ò e-ttα dt = Г+1),

α +1 > 0 ,

 

 

 

γ

 

 

 

0

 

 

 

где Г – гамма-функция Эйлера.

 

 

 

 

 

 

. Г+1)

 

 

 

Таким образом, tα =

 

 

 

, α > −1,

Re p > 0 .

 

 

pα +1

 

 

 

 

×

 

 

 

 

 

20. Свойства преобразования Лапласа. Будем обозначать через f (t), g(t), ... оригиналы, а через F( p), G( p), ... – их изображения.

Непосредственно из свойств интегралов

¥

¥

F( p) = ò f (t)e- pt dt,

G( p) = ò g(t)e- pt dt, ...

0

0

получаем следующие свойства преобразования Лапласа.

1) Линейность. Для любых комплексных постоянных α и β

.

α f (t) + β g(t) =α F( p) + βG( p) , т.е. линейной комбинации оригиналов

×

соответствует такая же линейная комбинация изображений.

.

1

æ

p ö

 

2) Подобие. Для любого постоянного α > 0 f t)=

 

F ç

 

÷

,

α

 

×

è

α ø

 

т.е. умножение аргумента оригинала на положительное число α приводит к делению изображения аргумента на это число.

Доказательство. В самом деле, полагая αt = τ имеем

.

¥

f t)e- pt dt =

1

¥

 

-

p

τ

 

1

æ p ö

 

 

 

f t) =

 

ò

f (τ )e

α

dτ =

 

F ç

 

÷.

α

α

 

×

ò

 

 

 

 

 

 

è α ø

 

0

 

 

0

 

 

комплексного

3) Смещение. Для

любого

 

числа α

.

eαt f (t)= F( p -α) , т.е. умножение оригинала на функцию eαt влечет

×

за собой смещение переменной p.

313

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