Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Уравнения математической физики..pdf
Скачиваний:
29
Добавлен:
15.11.2022
Размер:
5.66 Mб
Скачать

Глава 1. КЛАССИФИКАЦИЯ УРАВНЕНИЙ ВТОРОГО ПОРЯДКА. ОБЩЕЕ РЕШЕНИЕ УРАВНЕНИЯ

 

Обозначим через D область «-мерного евклидова пространства Еп

точек

JC= ( д г | ) с

декартовыми ортогональными

координатами

хь ...,хп, п ^ 2 .

 

 

 

 

 

Пусть F(x,...,pi

h

- заданная действительная функция точек

хе D

и действительных

переменных Piwj n с неотрицательными

целочисленными индексами

П

___

 

ц,..лп, £ ij = к,

к = 0,т,т>\, по крайней

 

 

 

М

 

 

мере, одна из частных производных которой

3F

 

* 0,

>1

 

Уравнение вида

 

 

dPib-in

 

 

 

 

 

 

F

( x

■9 U , ,..•) = О,

xeD .

(1)

 

 

 

cbcj1...дх1”

 

 

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

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

 

Уравнение (1) называется линейным, если F линейно зависит от

 

д^и

0<к<т.

всех частных производных —:----- —,

 

дх1...дх1"

 

 

Линейное уравнение можно записать в виде Lu = / (х), х е D , где

L

дифференциальный оператор первой степени относительно всех

частных производных. Линейное уравнение будет однородным или неоднородным в зависимости от / (х) = 0 или / (х) * 0 .

Уравнение (1) называется квазилинейным, если F линейно зависит

лишь от

дх‘п

1 (Г т .

 

 

 

 

ах{‘ ..илп

я

 

 

 

 

 

 

1.1. Классификация уравнений в точке

 

 

Рассмотрим квазилинейное (линейное относительно всех старших

производных)

дифференциальное

уравнение

второго

порядка

относительно функции и(х) = и(х},х2,...,хп)

 

 

 

 

 

 

д2и

ди

ди

 

 

У &у (*1»Х2>• • • » )

dxjdxj + Ф

Х|,..., Xfl, дх.

дху

= 0 ,

(1.1.1)

‘J=]

 

ау(х)

 

Ч

“лп У

 

где коэффициенты

непрерывно дифференцируемые

функции

независимых переменных х = (х],х2,...,хп) в некоторой области G. Ф

заданная функция своих аргументов. Уравнение (1.1) на самом деле

содержит при i ф j

не отдельные слагаемые ау

д2и

и а

д2и

, а их

 

lJ

dxjdxj

 

Jl dxjdxj

сумму (ау +aji)

д и—. Выражение ау + ajk

можно

разбить

на два

Ji' dxjdxj

слагаемых каким угодно способом, и будем считать, что ау(х)= а^(х).

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

обладают различными аналитическими свойствами.

 

Зафиксируем точку х° е G

и рассмотрим уравнение (1.1.1) как

уравнение с постоянными коэффициентами ау = а(у (х°):

 

0

д и + ... = 0 ,

(1.1.2)

1 «.'J

дх :дх,

 

Ненаписанные члены уравнения не содержат производных второго порядка. Сделаем замену переменных при помощи линейного преобразования:

Ук = с к\х \ + с к2х 2 + • • • + с кпх п

= 1,2

(1.1.3)

Предполагаем, что преобразование

(1.1.3) неособенное, т.е., что

определитель |сЛ(| не равен нулю. Производные по старым переменным выразятся через производные по новым переменным по следующим формулам (и{у\,у2,...уп) = и(хх,Х2,...,хп))\

 

 

 

ди

 

ди

дук

 

 

ди

 

 

 

 

 

 

 

tei

~ 2~i а ,,

^

~ 2-jcki

$Ук

 

 

 

 

 

 

 

ы \дУк dXi

к=1

 

 

 

 

 

 

д2и

_

д2и

 

ду,

_

^

 

 

д2и

 

 

 

 

дх,-9 *j

к М

& к f y i

Ckl

d x j

 

k j l IСк' C>J

 

дУкдУ/

 

 

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

 

 

 

 

 

 

 

п

0 О U

 

п

 

п

 

 

д2и

 

 

 

 

 

Z

 

О

 

 

 

 

...=

 

 

 

а и ^ Г Т 7

+ - - = L

a v

L

c kiCtj

 

+

 

 

 

дукду,

 

 

 

 

 

i,j=\

д х & 1

 

i,j=\

к м

 

 

 

 

(1.1.4)

 

п

( п

 

д2и

 

 

п

 

 

 

 

 

 

,

 

л ()

и

 

 

 

 

= Z

T al ckicij

 

 

 

V

+... - О,

 

 

 

+ ... =

 

a kl

 

 

 

 

kj=\\i,j=\

 

дукду,

 

kj=1

дукду.

 

 

где

новые

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

Щ выражаются через старые

согласно

формулам:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2$ = Z

a fjc kic ij

 

 

 

(1-1-5)

 

 

 

 

 

 

<J=i

 

 

 

 

 

 

 

 

Рассмотрим соответствующую уравнению (1.1.2) квадратичную

форму от вспомогательных переменных t s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0-16)

 

 

 

 

 

 

'.7=1

 

 

 

 

 

 

 

 

Если в квадратичной форме (1.1.6) перейти к новым переменным

rs

при помощи матрицы, транспонированной к матрице

\\ckj\\

в (1.1.3),

выражая старые переменные ts через новые гЛ

 

 

 

 

 

 

t, =сит] +с2,т2 +... + сп,тп

(/ = 1,2,...,и),

(1.1.7)

получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

'>7=1

= Z 4

Z

cfor/ tZ c/yr/

/,y=l

U=1

/=1

n (

n

 

 

= Z

Z 4 jCkiClj

TkTl = Z 5* / ^ ,

Л,/=1\/,у=1

у

A./=l

т.е. преобразованная квадратичная форма (1.1.6) будет иметь как раз коэффициенты 5$, определяемые формулой (1.1.5).

Итак, чтобы упростить уравнение (1.1.1) в точке х° с помощью замены переменных (1.1.3), достаточно упростить в этой точке квадратичную форму (1.1.6) с помощью неособенного линейного преобразования (1.1.7). Но в курсе линейной алгебры доказывается, что всегда существует неособенное преобразование (1.1.7), при котором квадратичная форма (1.1.6) принимает следующий канонический вид:

5 > /2 -

2 >,2, т<п

(1.1.8)

/=1

/=/-+)

 

Кроме того, в силу закона инерции квадратичных форм целые

числа г и я ? не зависят

от преобразования (1.1.7).

Это позволяет

классифицировать дифференциальные уравнения (1.1.1) в зависимости от значений, принимаемых коэффициентами ау в точке х°

Если в квадратичной форме (1.1.8)

т = п и все слагаемые одного

знака (т.е. либо г = т, либо г = 0 ), то

уравнение (1.1.1)

называется

уравнением эллиптического типа\ если

т = п, но имеются

слагаемые

разных знаков (т.е.

1 < г </7- 1), то уравнение (1.1.1) —гиперболического

типа (при г = 1

или г = п 1

нормально-гиперболического типа);

наконец, если т < п, то уравнение (1.1.1) - параболического типа (при

т —п —\ и г -1 или г = п -1 - нормально-параболического типа).

При произвольной неособенной замене независимых переменных

у = у(х) в уравнении (1.1.1), т.е. при

 

Ук=Ук(*\'Х2,...,х„) (к =1,2....,я),

(1.1.9)

D _ д{У\’Уг,--;Уп) Q

 

Ук

д(х],х2,...,хп)

формулы (1.1.5) примут вид

 

 

 

 

 

2И(у )= 2 > *

dxt dxj

 

(1.1.10)

 

 

 

 

 

 

 

 

 

 

 

При

этом,

так

как

D * 0,

то в

некоторой

окрестности можно

выразить переменные х через переменные

у,

х = х(у).

Обозначим

ы(х(д/)) = «(у).

 

 

 

 

 

 

 

 

 

 

Для того чтобы одним и тем же преобразованием (1.1.9) можно

было привести уравнение (1.1.1) к каноническому виду

 

 

 

 

 

г

д2~

т

и U

 

 

 

 

 

 

 

z/=1

 

v-'

 

 

 

 

 

 

 

dyf

i=r+1 dyf

+ . . . =

о ,

 

 

 

 

 

 

 

 

 

 

необходимо, чтобы число условий в соответствии с (1.1.10)

 

 

 

 

ак/= 0, / * к,

1,к = 1,2,...,и;

 

(1.1.11)

 

 

 

a,i =£iau ,

I = 2,3,...,л;

£ ,,*0,

 

 

 

 

где

£•/ = 0, ± 1,

не

превосходило

числа

неизвестных

функций

у к,

к - 1,2,...,л:

 

 

 

 

 

 

 

 

 

 

 

 

 

я(и - 1) + н -1 < я,

т.е. п <2 .

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

В

противном

случае

 

система

уравнений (1.1.11) -

переопределенная и потому, вообще говоря, не разрешима. Если коэффициенты уравнения (1.1.1) постоянны, то, приводя (1.1.1) к канонической форме в одной точке, мы получим уравнение, приведенное к канонической форме во всей области G определения уравнения.

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

и ^ + 2иху - 2иХ2 + 2Uyy + 2uzz = 0 .

Решение. Составим соответствующую уравнению квадратичную форму (характеристическую форму) от вспомогательных переменных xs :

 

 

 

з

 

/

= х 2 + 2xlx 2- 2xixi + 2*2 + 2*3 = ^О уХ ^у

 

 

 

 

<J=i

 

Причем,

так

как

ау = ajt,

то

<*\\ =* 1, «12 = 1, ^i3 = -1, а22 = 2, «зз = 2. Поскольку, например,

яп * 0,

то, как легко

проверить,

выражение

а ^(а п х^ + ai2x2 +

+а\„хп)2

содержит такие же члены с переменной * |, как и форма / = 'y'j aiixjxj , а ',7=1

поэтому разность

п\ 2

/- « и X «i к*к Ук=1

будет квадратичной формой, содержащей лишь переменные х2,...,хп, но не X]. Отсюда

п"\2

/ = «11 2 > 1**А + £ U =1

Если ввести новые переменные

У\ = «11*1 "*■«12*2 + •■• + «!пхп,У1= *,при / = 2,3, то получим

/ = «nVi2 + S >

где g будет теперь квадратичной формой от переменных У2>Уз>--->Уп >т-е- будет зависеть от меньшего, чем п, числа переменных. Повторяя указанную процедуру нужное число раз, мы приведем квадратичную форму/ к каноническому виду.

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

( з

"N2

 

S

а\кхк - (*i +х2 ~ хз)

и

 

g ~~ х\ 2х,х2 2Х|Х3 + 2х2 + 2х3 —(xj + х2 - х3)

= *2 + *3 + 2х2х3. Заменим переменные по формулам

У] =Х]+Х2 - х 2, У2 = х2’

Уз= хз

с матрицей

ю

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иу ,