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

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

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

X 1

Я-1 = 0 1 0 5 в ~ х =1 , в =

0 1J

( 1

-1

П

о

 

О

0

\ )

Здесь матрица В~1 является обратной матрице В, которая определяет в равенстве (1.1.7) преобразование старых переменных через новые. Получим

/ = а\\У\ + g = У\ + У2 +2У2Уз + Уз »

где bn = 1, Ъ22 = 1, b23 = 1, Ь33 = 1. Тогда

 

 

 

 

у

з

Л2

 

 

 

 

 

 

 

 

^22

Z

62кУк

={У2 +Уз?>

 

 

 

 

 

 

\ к = 2

 

 

 

 

 

 

И

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g' = y 2 + 2у 2у 3 + у 3 - ( у 2 +у3)2 = 0 .

 

Совершая

 

линейное

преобразование

zt = у,; z2 = у 2 +у3; z3 = у3 с

матрицей

 

 

г\

0

 

 

г\

 

0 "

 

 

 

 

 

0

 

 

С" 1 = 0 1 1 5 С- 1 = 1, с = 0 1 -1

 

 

 

,0

0

 

 

 

 

ч°

0

и

 

 

2

2

 

 

 

 

 

получим /

-

 

 

 

 

 

 

 

 

 

= Z ]

+ z 2 .

 

 

 

 

 

 

 

 

Линейное преобразование, приводящее квадратичную форму к

каноническому виду, будет иметь своей матрицей

 

 

 

 

 

 

(\

- 1

п (\

0

0 " (\

- 1

2 N

 

 

С==ЯС =

0

1

0 0

1

- 1

- 0

1

- 1

 

 

 

 

,0

0

и 1° 0

К

,0

0

l j

Всякую квадратичную форму можно представить в виде матричного

т

равенства f = X А Х , где

f *x'

х = х 2 ,

^ r =(x!

*2

 

xn), A

симметрическая

матрица

\ XnJ

 

 

 

 

 

 

 

 

 

 

 

 

квадратичной формы f причем Ат = А .

 

 

 

 

 

В нашем примере

 

 

 

1

1

- П

 

 

 

 

 

 

 

 

 

 

х\

 

 

 

 

/

= (х,,х2,х3

1

2

0

х2

 

 

 

 

 

 

 

 

 

-1 0

2 , КХ2

 

 

 

Линейное

преобразование

X = GZ,5

X T = Z TGT приводит

форму / к

виду

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/

= Z T[pTAG^fZ,

 

 

 

 

где

 

 

 

 

 

1 - Г | 1 - 1

 

 

 

 

1

0

0" 1

 

2 N ( \

0 °1

(G T A G )= - 1

1 0

1

 

2

0

0

1 - 1 = 0 1 0

 

U

- 1

и 1 -1

0

2J1° 0

1 , 10

0

o j

 

 

Таким образом,

 

 

 

 

 

 

 

 

 

 

 

 

 

"1

0

0

"

 

 

 

м

 

 

/

= (21>22>г з)| 0

1

0

 

z2

= (^b Z 2,23l

= z \

+ z 2 ■

 

 

,0

0

0

,

 

 

 

г

2

 

 

 

 

 

 

 

l

o j

 

 

а матрицей линейного преобразования переменных в дифференциальном

уравнении в соответствии с равенством (1.1.3) является матрица GT или, выражая новые переменные через старые,

[а = х

j Р = -х + у

 

у = 2х - у + z.

Делая замену переменных

в исходном уравнении, получаем

( u(a,/3,y)=u(x,y,z)):

 

их =иа -ах + ир

{Зх +иу -ух = йа - й р +2иу ,

uxx = i ux)a

а х + { их ) р Рх + (их)у -Ух =

~ ^аа ^рр

уу ~ 2wар + 4wау — Аиру,

Uу — Up — Uy ,

Uyy — иар — иау — ирр — 2Uyy + 3иру,

ыуу = и рр ~ 2мру +Uyy,

uz =uy ,

uxz = иау —иру +2Uyy,

uzz ~ Ууу •

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

Масс иРР = ^ •

Это в соответствии с классификацией (1.1.8) означает, что исходное уравнение является уравнением параболического типа.

Крассмотренному алгоритму следует добавить, что если

квадратичная форма /

не

содержит квадратов переменных, т.е.

а\\ = а22 ~ ■■• ~ апп = 0 ,

то

предварительно

нужно

совершить

вспомогательное линейное преобразование:

 

 

 

X, =у, - У 2

 

 

*2 =У\ +3'2

xi=yi> / = з,...,«,

приводящее к появлению в форме / квадратов переменных. Квадратичную форму можно привести к каноническому виду

также, например, методом выделения полных квадратов.

1.2. Классификация квазилинейных уравнений второго порядка с двумя независимыми переменными

Пусть а, Ь, с - заданные вещественные функции двух переменных x , y s R , F заданная функция пяти вещественных аргументов. Рассмотрим следующее уравнение относительно неизвестной функции

и(х, у):

aitxx +2buxy +cuyy +F(x,y,u,ux,uy) = 0,x,yi=R.

(1.2.1)

1. Выделяем множество A c R 2 тех значений аргумента (х, у), при которых функции а, Ь, с определены и непрерывны, а также определена и

непрерывна функция F:Ax.R? —>R . Очевидно, на множестве А можно ставить задачу о классических (регулярных) решениях уравнения.

2. Выделяем множество D a А тех точек (х, у), при которых функции а, Ь, с не обращаются одновременно в нуль, то есть

D := {(х, у) е А \ \ а{х, у) \+1Ь(х, у) j +1с(х, у) |* 0 .

Очевидно, в D уравнение (1.2.1) является линейным относительно старших производных уравнением второго порядка с двумя независимыми переменными. Говорят, что точки множества A \D - точки вырождения уравнения (1.2.1), или что на множестве A \D уравнение вырождается.

3. На множестве D классификация уравнения производится по

знаку функции Д:

 

 

 

 

 

 

Д :=Ь2 -ас.

 

(1.2.2)

 

Уравнение (1.2.1) называют уравнением

 

 

a) гиперболического типа, если Д > О,

 

 

 

B) параболического типа, если Д = О,

 

 

 

c) эллиптического типа, если Д < 0.

D+ := {(х, у) е D | Д(х, .у) > 0},

При

этом

множества

DQ := {(х, у) D | Д(х, у) = 0},

D_ := {(х, .у) е D \Д(х, у) <0}

называют

соответственно областью гиперболичности, параболичности или эллиптичности уравнения.

Очевидно, в частном случае уравнения (1.2.1), когда а, Ь, с

константы, |а | + | 6 | + |с|* 0

и F непрерывна на R5, имеем A = D = R2 и

уравнение на всей плоскости имеет только один из типов.

Задача 1.2. Определить тип уравнения

4Uxx + cos х sin2_y

+у ех+уих + 2cos х sin у иу +и2 = 1995,

где у - const. Изобразить графически области, где уравнение имеет тот или иной тип.

Решение.

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

V(x, у) е R

, так как всюду а Ф0.

Изобразить графически области, где уравнение имеет тот или иной тип. а)(1 - х 2)ихх - 2 иу - ( у 2 +\)Uyy +2хих - 2 ху иху = 0 ;

б) 15Uyy + 26sin(x3.y3) = 64 +х тихх, т = 0,1,2,...

 

 

в) 4sinхиху -

8 - 5(их +иу )~ sign х - и ^ - sign у •

;

 

г) ихх + cosхиу - (cos2 х - sin2 x)uyy + 2sin x • uxy = ( - l)slgn(j°4

 

1.3.

Приведение уравнений к каноническому виду

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

 

 

 

 

аихх +2 buxy + cUyy +F(x,y,u,ux,uy ) = 0>

(1.3.1)

в некоторой точке (х, у) е D (D - область, где уравнение определено и не

вырождается).

Рассмотрим

некоторое

преобразование

со : R2 R2 = {£ = <*(х,у),Т] = rj(x,y)}

(1.3.2) переменных,

для которых

якобиан преобразования

 

 

 

 

/ =

 

 

(1.3.3)

 

Лх

 

 

 

(это условие обеспечивает взаимооднозначность преобразования). Можно перейти к новым функциям v = u°co: v(£, 7) = w[x(£, 7), y(g, 7)] и записать новое уравнение, эквивалентное (1.3.1).

Используя известные в курсе математического анализа формулы

получим

 

«х =£cv£ + 7 * V иу

 

 

+ W

(1-3-4)

д

 

 

д

 

 

 

д

 

 

 

 

 

 

 

 

UXX ~

(4xV^ Лх^Ц ) —

4xXVtj

(v £ )

ЛххУт) "** Лх

(V4 ) -

~ 4 x x v Jj

4 х i 4 x v ^

Лх^^Т) )

Лхх^1]

Л х ( 4 x ^ r j ^

Л х ^ щ

)•

Приводя подобные и используя равенство

 

, получаем

выражение

ихх

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

и 7 . Аналогично можно

получить

выражения для

и ^ и

Иуу

 

(предлагаем

сделать это

самостоятельно). Получаем

PNRPU

^44 ^^х^1х^4ч ~^^Лх) Уг/г) ‘зд;х^4 ~^^XXVTJ

Ихх — ( £ х )

Uxy=Zxtyv4S +(%xrly+rix%y)v4ri+T1xJlyVnT + V ^ W ’ О '3'5)

=у) ^y^lyv4rl + ( ^ ) V77 ^yyv4 ^yyVj]'

Умножая каждое из уравнений (1.3.5) на коэффициенты а, 6, с соответственно, складывая их и подставляя в (1.3.1), используя (1.3.4), приходим к уравнению

av4f + 2bv4n +cvnr]+pv4 +qvn +

(U>6)

+ F(x,y,v,<*xv4 +rixvn,Zyv4 +TiyVTJ) = 0,

 

где

 

 

 

a = a{tx)2 +2btxSy +c(ty )2,

 

Ь ~ ®%xVx

b{£xT)y + Tjx4У) + C^xTJy,

 

C =a(7ixf +2bTjxj]y+c{T)y)2,

(1.3.7)

P =

a Zxx +

2 b £ xy + C % y y ,

 

q -

ат/хх + '^br]Xy + crjyy.

 

Далее осталось в выражении (1.3.6) сделать замену переменных (1.3.2). Получим уравнение относительно переменных (4, tj) :

AV44 +2Bv4,1+Cv,in+F(4,T],v,v4>vt1) = 0.

(1.3.8)

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

 

аим +2Ъиху +сиуу +dux +еиу + fu + g = 0,(x,.y) е D

(1.3.9)

где a,b,c,d,e,f,g

известные функции аргументов (х,у), которые

являются частным случаем уравнения (1.3.1). Его можно привести к виду

av44 +2bv4n +cvw +dv4 +ev,?+ Jv + g = 0,

(1.3.10)

где

a = a ( ^ ) 2 + 2 4 ^ + c ( ^ ) 2,

b = a4xTjx + b(4x7]y + Tjx4y ) + c4xj]y,

c=a(TJx)2 +2bqxqy +c(r)y)2,

(1.3.11)

d ~ a Zxx + 2 b 4 x y + c %yy + d % x + e 4 y,

e =ajixx+2brixy + c?jyy + dr]x + ет]у .

От уравнения (1.3.10), используя (1.3.2), легко перейти к уравнению, содержащему только переменные (£,77). Замена переменных (1.3.2), очевидно, имеет смысл только тогда, когда новое уравнение получается проще исходного. Как известно, если преобразование (1.3.2) выбрать

таким образом, чтобы функции £(х,у), г/(х,у)

удовлетворяли уравнению

характеристик

 

a dy - {Ь ± 4 &)сЬс,

(1.3.12 ±)

то мы придем к уравнению вида (1.3.8), который называется каноническим * 0). Процедуру выбора т] и канонический вид для уравнений каждого типа приводим в таблице.______ ________________

Знак Д

Д > 0 -

 

Д = 0 -

Д <0 -

 

гиперболический

параболический

_эллиптический

 

 

тип

 

тип

тип

 

 

 

%(х,у) =С{,

 

£(х,у) = С

4(x,y) + itj(x,y)

^II

 

TJ(X>у) - С2

 

решение

= с

 

р

решения

 

(1.3.12+), т](х,у)

решение

 

 

 

II

р

уравнений

 

выбирается

любого

из

 

 

(1.3.12+)

и

произвольно, но

уравнений

 

 

 

(1.3.12-)

 

чтобы I * о

(1.3.12)

 

Каноничес

 

 

 

v# + V / +

 

кий вид

V4JJ + ф (-) = 0

 

+ ф (‘) = 0

 

 

+ Ф(-) = 0

 

 

 

 

 

 

 

Примечание: здесь 0(-) =0 (^ , TJ,V,V^ , V}J) .

Таким образом, приводить к каноническому виду уравнение (1.3.10) нужно отдельно в каждой области, где уравнение принимает тот или иной тип, и придерживаясь следующего алгоритма.

1.Найти £(х, у), т]{х,у) согласно таблице.

2.Найти коэффициенты a,b,c,d,e согласно (1.3.11); причем можно находить не все коэффициенты, а лишь некоторые, ориентируясь по каноническому виду из таблицы, например, для гиперболического

типа достаточно найти b,d,e (т.к. а = с = 0). 3. Выписать (1.3.10).

уравнения.

Замечания. 1. Очевидно, описанная схема годится и для уравнения (1.1), при этом нужно воспользоваться (1.3.6) и (1.3.7).

2. Можно не использовать формулы (1.3.11), а выражать каждый раз их,иу ,ихх,иху,иуу по формулам (1.3.4), (1.3.5), подставляя затем

выражение в (1.3.1). Очевидно, это более громоздкий способ решения, но он тоже приведет к каноническому виду.

3. Уравнение v^ + Ф(-) = 0 можно с помощью замены переменных

а = £ + 7

*

£ - 7 привести к уравнению м>аа + wpp +Ф(-) = 0, которое

2

, Р

2

 

 

 

 

часто

также

называется

каноническим

видом

уравнения

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

 

 

 

4. Если а = 0 , то возможны две ситуации:

 

 

а)

е = 0

и

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

приводится к каноническому виду простым делением на Ь\

 

б)

с * 0 ,

тогда уравнениями (1.3.12) можно воспользоваться,

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

и(х,у) =и(у,х) в исходном уравнении, в

котором коэффициент при ихх уже отличен от нуля.

5. Очевидно, что канонический вид уравнения в каждой точке определяется неоднозначно, хотя бы потому, что неоднозначна замена переменных (например, с точностью до константы), поэтому в ответе необходимо делать сноску, при какой замене переменных уравнение имеет данный канонический вид.

Задача 1.4. Привести уравнение к каноническому виду.

ихх + 2Uyy + 5Uyy - 32и = 0

Решение. Уравнение имеет эллиптический тип всюду в R2 так как Д = -4<0. a = \,b = \,с = 5,d = е = g = Q,f = -3 2 .

Уравнения характеристик:

dy = (1 ± yf-A )dx, dy = (1 ± 2i)dx, C = y - x { \ + 2i).

Отсюда выберем следующую замену:

£ = у - x,

TJ 2х,

 

6 г = "Ъ

Vx = 2’

 

^

= 0.

 

Вторые производные равны нулю.

 

 

Учитывая эллиптичность уравнения, а=с,Ь= 0 . По

формуле (1.3.11)

с =4,b =d =е =0. Выписываем уравнение (1.3.10) 4

+4м,;7? - 32м = 0.

Ответ:

и ^ + ищ - 8м = 0, (х, у) е R 2, при замене переменных £ = у - х, г) = 2х.

Задача 1.5. Привести уравнение к каноническому виду

 

У2ихх + IxyUjy + x2Uyy = 0 .

9

= ху,с = x

9

,d - е =/= g = 0,А = 0 .

Решение: а = у , Ь

 

1.В точке (0, 0) уравнение вырождается.

2.В области D = R2 \ {(0,0)} уравнение имеет параболический тип:

а) если у =0, то из уравнения следует x2Uyy - 0, так как (х, у) е D, то х * 0, получаем канонический вид и ^ =0 при у - 0;

б) пусть у * 0. Уравнение характеристик: у

dy = xydx. Решая его, получим

 

2

2

 

 

 

 

 

 

 

 

С = ——

— . Возьмем £ = у 2 - х2 Очевидно, можно взять rj = x,

так как

в этом случае

- 2х

= -2у*0

 

 

 

 

 

 

 

 

 

1

0

 

 

 

 

 

 

^>х

2х, т]х 1,

2у, Tjy 0, £хх

2, т]хх

0,

0,

0,

 

%уу= 2; т]уу = 1;=>

Параболический

тип

 

=>а=Ь= 0

Из

(1.3.11)

с - у 2,d = -2 у 2 + 2х2,е = 0. Выписываем вид (1.3.10)

 

 

 

 

 

 

y 2vw + 2(x2 - y 2)v£ = 0 ,

 

 

 

 

 

 

v77 + 2(X

^

'

v<? = 0 -

(i-3-17)

Рекомендуем, как и здесь, делать окончательную замену

х,у на £,77 в

конце

решения задачи

после получения вида

(1.3.17). Это избавит от