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

книги / Уравнения математической физики методы решения задач

..pdf
Скачиваний:
14
Добавлен:
12.11.2023
Размер:
5.73 Mб
Скачать

(-у)

- 1

I 1

= у,

2

<р(у) + И -

 

 

ку)

 

1 (- у )~2<р(у) + (- у )~2 <р'(у) ■у +- у/'Г1 Л= У-

2

у

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

_1

¥ = у - ( - у ) 2<р(у)>

.У>

' ¥ f-1 = i + i ( - y f 2iK y)-(-y)'V (y)-

У

Отсюда

<р\у) - (-у)2 => (Ку) = - 1 (-у)2 + с

r t y ) = У + j (-У ) 2 ( - У )2 ч - у р С => <Ку) = | - ( - у ) 2 с .

Далее получаем выражения для

<р(ху),ц/ГУ) . Подставляем их в общее

решение.

 

 

 

V *;

 

 

 

 

 

 

 

 

 

 

2 4 х

\

4f

| х

у

\ J L c =

l x i y + i -

и(х,у) =-

х г(~У)г + J

С +У-

3 4~У

 

 

- у

Зх

V - у

3

Зх

Задача 2.8. Решить задачу Коши

 

 

 

 

 

х2ихх2ху -иху- 3 у 2иуу = 0,

х > 0, у > 0,

 

«Iy=\=f(x),

uy \y=x=g(x),

 

 

 

где / (x), g(x) е С2(0,оо) - заданные функции.

Схема решения. Обычным образом находим общее решение уравнения, получим

и(х,у) = х 4у3 4<р X + ¥ (хъу), (2.2.9) \У )

или

Учитывая предстоящую подстановку, удобнее воспользоваться видом (2.2.9) общего решения. Однако далее мы будем умышленно исходить из вида (2.2.10), предоставляя студентам самостоятельно проделать выкладки для случая (2.2.9) и сравнить полученные результаты.

В случае (2.2.10) система для определения <р,у/ примет вид

х*<р' О + у/ (х3 ) = / ( * ) ,

\ X J

 

I

 

4 (\

 

 

 

 

 

 

 

4

 

J^J+x V[^“j+*V(*3) =£(*)•

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

И х 3) = f ( x)-x*<P\ х )

 

 

(2.2.11)

 

 

 

 

(

 

 

3х 2у/'(х3) = f ' { x ) - ^ x

 

- 5

 

 

*<р( - ) + 1

V

 

 

 

 

 

 

 

UJ

 

 

подставляя полученный результат во второе равенство, получим

1 1 ( 1)

-1

/гПп

1

ч

1 4 Г О

1

-1

/ Р

( Л

4

 

-

+ - х / ' ( х ) - - х 4<р -

\ + - х

V

- = g(x),

- х Аср

+ Х

АФ

\ х )

 

 

4

m

 

 

 

 

\XJ

 

 

 

 

 

 

 

 

 

 

 

—v 4<р‘^ x j

 

 

 

 

 

Замечание. В отличие от предыдущих примеров, где для функций f{x)vig(x) было дано конкретное аналитическое выражение, здесь нельзя брать неопределенный интеграл от обеих частей равенства, иначе решение задачи Коши получим с точностью до константы, а нам нужно единственное решение. Поэтому возьмем интеграл с переменным верхним пределом интегрирования.

3 -1

Итак, домножим обе части последнего равенства на - —х 4 и

проинтегрируем его в пределах от а до х , где а >0 - произвольно.

h i

Л"

-1 d t = \

9>'

у

a L

<p

-<p

= - 7 f<f

+ 7 k ' v w j f ■

\ x j \ a )

4 aJ

4 aJ

(2.2.12)

 

 

Первый интеграл в правой части оставим без изменений, а во втором избавимся от производной под значком интеграла, взяв его по частям:

J* " W

)

= f ' №

*

+ \ j< f * /(£)<*£

 

Отсюда и из (2.2.12) получаем

 

 

 

 

ГО

+

 

+

 

(2.2.13)

;

 

16

и

 

 

 

 

y (£ ) = < f 4 l / ( £ ) - 4 g ( £ ) ] .

 

где С = ^ — - л ~ / ( а ) ,

 

w

 

 

 

 

 

 

Подставляя полученное выражение в (2.2.11), имеем:

 

 

)=

- * ‘ С

+ 1

/ ( х ) -

j'r ( 4 ) d 4

(2.2.14)

 

 

 

 

 

16

 

 

 

 

 

 

а

 

Из (2.2.13) и (2.2.14) получаем следующие аналитические выражения:

 

 

 

 

 

 

X

( у )

 

1

- i

1

+-

з

(р —

= С + - х

V

/ f-1

 

с

-

4 у4 1 —

16 а\г<ЛЩ

 

 

4

 

О ,

 

у) = -X V 4С - 1 Д ху5) - h

X-4у 4 Jr(€)d£

 

 

 

 

 

 

а

Подставляя полученные выражения в (2.2.10), получим:

и(*,У) = 7 /

ху3

1

( V

+ ^ У

| У 4|/( £ ) - 4 Д £ ) ] ^

+ тУ /

у

4

V )

4

16

J

 

 

 

 

 

ху%

зз

2.3. Решение задачи Гурса

Задача для уравнения гиперболического типа

ciUxx + 2Ъиху + сиуу + F(x, у,и,их,иу) = 0 (х, y ) e D a R 2,

u\y«,w - /W >

=

 

где q>и у/

две линейно

независимые характеристики уравнения,

называется задачей Гурса в узком смысле.

 

Задача 2.9. Решить задачу Гурса:

 

 

ихх + Зиху - 4Uyy - Ux + иу

= О,

 

и\ .

=5х + е , и\

=1.

 

\у= 4х

\ у = - х

 

Решение. С помощью замены £ = х +у, т?= у - 4х уравнение приводится к

виду

1

, после этого находим общее решение и получаем

v — v£ = 0

 

 

х + у

 

 

 

и = <р(у - 4х)е 5 + ц/{х + у ) .

 

Далее из дополнительных условий получаем

 

 

 

\<Р{0)ех + у/{5х) = 5х + ех

(2 3 1)

 

 

[р(-5х) + ^/(0) = 1

 

После подстановки в общее решение выражений для (р(у-4х)иу/(х + у ) , получим

х + у

х + у

х + у

и(х,у) = х + у + е 5

+е 5 - е

5 [^(0) + ^(0)]. (2.3.2)

Значение постоянной <р(0) + у/(0)

находим с помощью любого из

уравнений системы (2.3.1). Например, второго <р(0) + у/(о) = 1.

С учетом (2.3.2) имеем

 

 

 

х + у

и(х,у) = х + у + е

5

Задача 2.10. Решить задачу Гурса:

 

 

иху+хих =0,

х>0,у>0,

(2.3.3)

и|*=0 =<Р(У)’и\у=о = У'(Х)>

тде<р,ц/ е С2 (х >0) П С(х >0).

Решение. Поскольку функции (рпу/ непрерывны в точке (0,0), то из (2.3.4) имеем и(0,0) = #>(0) = ^(0).

Таким образом, целесообразно рассматривать два случая:

1)$9(0) * у/(0) . Задача (2.3.3), (2.3.4) не имеет решения;

2)<р(0) = ^(0). Тогда вначале ищем решение (2.3.3) в общем виде. Делаем замену v = ux . Отсюда

dv

— + xv = 0,

ду

dv

= -xdy, v

ln|v| = -ху + 1пС(х),

их(х,у) = е~хуС(х).

Интегрируя последнее равенство в пределах от 0 до х, получим:

u(x,y) = u(0,y)+ je

о

X

и{х,у) = (р{у) + \е~&С($Щ.

Тогда из второго начального условия имеем

л

р(0)+ J C K y f = (K*).

Дифференцируя по х, получим

С(х) = i//'(x),

 

 

X

 

 

 

 

л

 

 

 

 

и = <р(у) + JfT^V'(£)</£ = <Р(У) + Jе~&

 

= <р(у) ■

 

 

 

О

 

 

 

X

О

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ е“*У Су)|* + У\е~®у/(£Щ.

 

 

 

 

 

 

 

 

О

х

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и(х, у) = e ^ y /ix ) + (р{у) - <р{0) + у \е~&yi&dB,

 

Задача 2.11. Решить задачу Гурса:

 

 

 

 

 

 

м „,------— (их - и у )= 1,

у < - х ,х >

2,

 

 

ху

х - у

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и\

 

=0, и\

 

 

.2

 

 

 

 

у=-х

= 2 + 2у + \-у"

 

 

 

 

I

 

>х=2

 

'

2

 

 

Решение.

Сделаем

замену

V

С

помощью

указанной

замены

и = •

 

 

 

 

 

х - у

 

 

 

 

 

выразим

иху,их,иу

через

 

x,y,v и

частные

производные v

Затем

подставим полученные выражения в уравнение задачи Гурса. После приведения подобных

vxy = x - y .

Непосредственно интегрируя полученное выражение, имеем

v = ^

+ <р{х) + у/{у) .

 

Далее, учитывая первоначальную подстановку и = ■

получим

 

 

 

х - у

 

 

+

 

х - у

(2.3.5)

Из краевых условий находим

2

 

 

 

 

 

 

(р(х) + у/(-х) = х ,

 

 

 

<р(2) + у/(у)

 

-

1 2

 

у + г±_/r_±sj__ 2 + 2у + —у

 

2 - у

 

 

2

 

Из второго уравнения имеем

y/(-x) = 2 - X + —X1 (2 + x)-<p(2).

2

Подставив полученный результат в первое уравнение системы, имеем

<р(х) = хъ -

2 - х + —х 2

|(2 + х) + <р(2).

V

2

)

Подставляя выражения для

(р{х)\лц/{у)

в (2.3.5), после приведения

подобных получим

 

 

и(х,у) = {х + уУ

Рассмотрим задачу Коши:

 

 

и„ =а2ихя +

t> 0 ,x e R ,

(3.1)

w|,=o =<Р(Х), «*|,=0

x e R

(3-2)

Классическим (регулярным) решением задачи (3.1)-(3.2) называется функция u(x,t) е С2 (/ > 0) П С2 (t > 0), удовлетворяющая всюду при t > 0 уравнению (3.1) и начальным условиям (3.2).

Теорема 3.1. Если <peC2( R ) ,i/s e C \R ) ,fe C l(t> 0), то классическое решение задачи (З.1.), (3.2) существует, единственно и выражается формулой Даламбера

1

 

1

х+ш

 

 

u{x,t) = -

[<р{х+ at) + <р(х - at)] + —

fy ffld g +

 

 

 

. .

° x~at

(3.3)

J

 

. t x+a{t-x)

 

 

v

+ Г- i J/«,r)d £ /r.

l a 0 x-a(l-r)

Для указанных функций <p,y/,f определена правая часть (3.3), поэтому можно формально найти функцию по формуле Даламбера.

Оказывается, найденная функция u{x,t) при каждом фиксированном t> 0 удовлетворяет уравнению (3.1), за исключением разве что конечного числа точек разрыва, а также удовлетворяет условиям (3.2). Более того, найденная функция оказывается обобщенным решением задачи.

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

3.1. Распространение волн в бесконечной струне. Задача Коши

Заметим, что задача Коши (3.1), (3.2) служит математической моделью различных физических явлений. В качестве одного из примеров рассмотрим задачу о малых поперечных колебаниях упругой однородной струны.

Будем считать, что колебания струны совершаются только в

направлении перпендикулярном оси струны х и искомая функция смещения точек струны u(x,t) зависит только от координаты х и времени t . Считается, что постоянная Т и плотность струны р постоянны. Тогда в указанном случае задача (3.1), (3.2) есть математическая модель задачи

о малых

непрерывных колебаниях

упругой однородной струны, где

Т

Fix t)

F(x,t)

- линейная плотность внешних сил,

а2 = — ,f(x ,t) = ------—,где

<р(х) - начальное смещение струны, а <//(х) - начальная скорость.

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

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

 

х е [0,11

смещения определяется равенством <р(х) = • 2 - х ,

х е (\,2\

0,

х R \ [0,2}

В момент времени I = 0 струну отпускают,

в результате чего она

начинает движение. Изобразить графически профиль струны в моменты времени / = — , / = 0,5. Схематично указать процесс дальнейшего

распространения волны, указать его скорость.

1г Решение. Т = 27, р = Ъ=> а - — = 3. Математическая модель колебаний

струны - задача Коши:

 

 

 

 

 

и„ =9и ^ ,

 

t > 0,

хе R,

 

I

,

ч

I

п

(ЗЗЛ )

% =0=Ф)>

*4=0 =0’

 

где <р(х) = Ж[о,1](*) • х + *(1>2](*) • (2 - х ).

 

 

Решение задачи (3.1.1)

есть

частный

случай задачи

(3.1), (3.2) при

f( x ,t ) = y/{t) = 0 и кусочно-гладкой функции . Решение (3.1.1) согласно

формулам Даламбера (3.3)

примет вид

 

и(х, 0

= 1[<Р(Х+ 3t) + <p(x- 3/)].

(3.1.2)

1) t - 0 и(х,0) = <р(х) - начальное положение струны.

 

2) / =

12

1 Удобно сначала построить вспомогательные графики —р(х) сжатием

графика <р(х) вдоль оси ординат, затем сдвигом влево (вправо) на -

построить графики функций — х + - (соответственно — <р X). 2 V

После этого сложить оба графика.

 

1 Г

Р

х

 

Г ’б Г г *

X + — + -(р

3 ) 4

■) 2

 

Построения проводим аналогично предыдущему пункту, с той лишь