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

§2. Линейное однородное уравнение первого порядка. Квазилинейные уравнения с частными производными.

Уравнение Гамильтона-Якоби.

1. Линейное однородное уравнение первого порядка.

n

 

 

Xk

 

 

(1) Lu = ut + fk(t; x)uxk = 0

=1

 

 

или

 

 

(1) Lu = ut + (f; ru) = 0;

(t; x) 2 G ½ Rn+1; L =

@

+ (f; r):

@t

Здесь: f = (f1; :::; fn); fk(t; x) - некоторые известные функции,

r = µ

@

; :::;

@

; u = u(t; x) - искомое решение.

 

 

 

@x1

@xn

Сопутствующая система обыкновенных диф. уравнений:

 

 

 

 

 

(¤)

dx

= f(t; x) (c:c:)

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

Пусть

©(i)(t; x);

 

 

 

- какая-либо система функционально

i = 1; n

 

независимых первых интегралов системы ( ).

 

©

 

 

 

 

 

ª

¤

Замечание. 1) Функция ©(t; x) - называется первым интегралом системы (¤), если она тождественно не равна константе, но в то же время эта функция постоянна вдоль каждого решения

x= x(t) системы (¤).

2)Интегральные кривые системы (¤) x = x(t) называются характеристиками уравнения (1).

3)Об одном теоретическом способе нахождения системы функционально независимых первых интегралов. Пусть x = x(t; x0) -

1

Лекция №2, НГУ, ММФ, 2010

 

 

2

решение задачи Коши

 

 

 

<

dx

 

 

 

 

=t0

= x0;

x0 = (x10; :::; xn0):

x t

8 dt

= f(t; x);

: j

 

 

 

По теореме о неявно заданных функциях векторное уравнение x = x(t; x0) может быть однозначно разрешено относительно x0 : x0 = x0(t; x) и система ©xi0(t; x); i = 1; nª может быть взята в качестве системы функционально независимых первых интегралов векторного уравнения (¤).

Общее решение уравнения (1).

u = F ©©(1)(t; x); :::; ©(n)(t; x)ª; F - произвольная функция (достаточно гладкая).

Свойство любого решения уравнения (1): вдоль характеристики решение u постоянно. Далее уравнение (1) можно еще переписать так:

dudt = 0, где dtd = L - полная производная от u вдоль характеристики.

Задача Коши для уравнения (1):

(

(2)

Lu = 0; ujt=t0 = '(x);

где '(x) - некоторая гладкая функция.

Формализм построения решения задачи Коши (2):

8

>>©(1)(t0; x) = ©(1);

<

а) >.

>

:©(n)(t0; x) = ©(n):

Из этой системы находим зависимость x = X(1); :::; ©(n)):

б) Тогда решение задачи Коши (2) записывается так:

u = '(X(1)(t; x); :::; ©(n)(t; x))):

Лекция №2, НГУ, ММФ, 2010

3

Замечание. Уравнение (1) можно трактовать так: dudt = 0 - пол-

ная производная от u в силу системы (¤) равна 0. Это означает, что вдоль характеристики функция u постоянна.

Рассмотрим вместо уравнения (1) более общее уравнение:

(10) f0(t; x)ut + (f; ru) = 0:

Рассмотрим два предельных случая. 1-ый предельный случай.

Если f0 =6 0, то (10) перепишется в виде (1)

(100) ut + µ 1 f; ru= 0; f0

характеристики которого определяются из соп. системы:

(¤¤)

 

dx

 

1

 

 

 

 

 

=

 

f:

 

dt

f0

Удобно ввести параметр s:

 

 

 

 

 

 

 

 

 

 

 

ds

 

1

; sjt=t0 = 0:

 

 

=

 

 

 

dt

f0

Тогда система (¤¤) перепишется так:

 

8

dt

 

= f0(t; x);

 

 

 

 

 

ds

 

 

>

dx

 

 

 

f

t; x :

 

 

 

 

 

 

> ds

 

 

 

 

 

 

 

:

 

 

 

=

(

)

<

 

 

 

Заметим, что вектор f = (f0; f) = (f0; f1; :::; fn) определяет век-

 

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

x = x(s);

уравнения (10).

тор, касательный к

e

(s = s(t)

 

Говорят, что этот вектор задает характеристическое направление в точке (t; x). Если мы решаем задачу Коши для уравнения (10) ((100)) с данными при t = t0 : ujt=t0 = '(x), то гиперплоскость t = t0 ни в одной точке не имеет хар. направления.

Лекция №2, НГУ, ММФ, 2010

4

2-ой предельный случай.

Если f0(t; x) ´ 0, то (10) перепишется так:

(1000) (f(t; x); ru) = 0;

характеристики которого находятся из системы:

8

> dt = 0; < ds

> dx

: ds = f;

т.е. при t = const характеристики расположены в гиперплоскости t = const. Поскольку вдоль каждой такой характеристики u

постоянно, то следовательно задача Коши

(

(f; ru) = 0; ujt=t0 = '(x)

разрешима не при любой функции '(x). Промежуточный случай будет рассмотрен далее на примере.

До сих пор мы рассматривали данные Коши на гиперплоскости t = t0. Рассмотрим теперь так называемую обобщенную задачу

Лекция №2, НГУ, ММФ, 2010

5

Коши, которая ставится так:

(

(3)f0ut + (f; ru) = 0;

uj° = '(t; x); (t; x) 2 °:

Здесь ° - гладкая гиперповерхность с уравнением

e

¯

 

 

 

e

 

ª(t; x) = 0;

 

 

 

 

 

 

 

 

 

 

 

 

 

; :::; ªxn) = (ªt; rª). Сделаем

причем jrªj¯°

6= 0; rª = (ªt:ªx1

в задаче (3) замену¯

независимых переменных:

 

 

 

 

(+) (

x = x;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

» = ª(t; x); u(t; x) = u(»; x);

при этом:

 

 

@u

 

 

@u

 

 

@u

 

 

@u

 

@u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

ª ;

 

 

=

 

+

e

ª

 

;

 

 

 

 

 

 

 

 

 

@xek

 

xk

 

 

 

@t

 

 

t

@x

 

 

 

 

т.е.

 

 

 

 

 

 

e

 

u =

 

k

 

ª:

e

 

 

 

 

 

 

 

 

 

 

 

u +

@u

 

 

 

 

 

 

 

 

 

 

 

r

 

re

 

 

e

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Лекция №2, НГУ, ММФ, 2010

6

> 0

= 0

<0

Следовательно задача (3) перепишется так:

(30)

[f0ªt + (f; rª)]u» + (f; ru) = 0;

(u

»=0

= '(t; x); t = t(0; x)

e

 

j

 

e

(заметим, что из (+)eследует, что t = t(»; x), если ªtj° 6= 0, напри-

мер).

 

 

 

 

Задача (30) однозначно разрешима, если

 

 

[f0ªt + (f; rª)] j° 6= 0;

 

т.е. (f; rª)j° 6= 0; f = (f0; f):

 

Это e e

e

 

 

 

означает, что вектор fe не лежит в касательной гиперплоскости к гиперповерхности ° (иными словами, ни в одной точке поверхность ° не имеет характеристического направления).

Примеры:

1) xut ¡ (t + 1)ux = 0; x 2 R1; уравнение характеристик: dxdt = ¡t +x 1,

т.е. общее решение: u = F(x2 + (t + 1)2):

Найдем решение задачи Коши при t > 0 с начальным условием:

Лекция №2, НГУ, ММФ, 2010

7

ujt=0 = x:

Однако простые рассуждения показывают, что начальное условие можно задавать либо при x < 0, либо при x > 0 (на всей оси t = 0 начальное условие задавать нельзя!). Если начальное условие

задается при x < 0, то решение имеет вид:

p

u = ¡ x2 + (t + 1)2 ¡ 1; t > 0:

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

t

(f ,f )

0

0

1

x

-1

всех x, заключается в том, что в точке (0; 0) линия t = 0 имеет характеристическое направление.

2)ut + ux = 0; u = F(x ¡ t); ª = x ¡ t; f0ªt + f1ªx = 0:

2.Квазилинейные уравнения с частными производны-

ми.

(4) Luk = gk(t; x; u); k = 1; m;

@

L = @t + (f; r); f = (f1; :::; fn);

= f(t; x; u);
= g(t; x; u); g = (g1; :::; gm)

Лекция №2, НГУ, ММФ, 2010

8

fk = fk(t; x; u); k = 1; n;; u = (u1; :::; um).

Система (4) называется квазилинейной. Если f = f(t; x), то система называется почти линейной.

а) Нахождение общего решения.

8 > dx

(5) < dt

> du

:

dt

(соп. сист. об. диф. уравнений).

Интегральные кривые системы dxdt = f(t; x; u) называются харак-

теристиками системы (4). Но в отличии от лин. уравнения (1), в квазилинейном случае нельзя найти характеристики, не зная

решения u = u(t; x). Каждое уравнение системы dudt = g(t; x; u)

называется соотношением на характеристике. Пусть f©(i)(t; x; u); i = 1; n + mg - какая-либо функционально независимая система первых интегралов системы (5). Тогда общее решение системы (4) дается в следующем виде:

8F1 n©(1)(t; x; u); :::; ©(n+m)(t; x; u)o

= 0;

>

 

 

 

 

>

 

 

 

 

>

: : : : : : : :

 

(t; x; u)o

= 0;

>Fm n© (t; x; u); :::; ©

 

>

 

 

 

 

>

(1)

(n+m)

 

 

:

 

 

<

 

 

т.е. функции uk; k = 1; m определяются неявно.

б) Решение задачи Коши

(

(6)

Lu = g;

ujt=t0 = '(x) = ('1(x); :::; 'm(x));

Лекция №2, НГУ, ММФ, 2010

9

строится так:

8

©:

(1): :(t:0:; x;: :u:) = ©(1);

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

(n+m)

 

 

 

 

 

(n+m)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

<

 

 

 

 

 

 

 

т.е.

(1)

 

>

©

 

(t0

; x; u) = © ;

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

:

 

 

 

 

 

 

 

 

x = X

 

; :::);

 

 

 

 

 

 

 

 

 

u = U(1); :::):

Тогда решение задачи Коши (6) дается в виде:

U(1)(t; x; u); :::) = '(X(1)(t; x; u); :::));

т.е. определяется в неявном виде.

Пример.

8 ut + uux = 0;

> dx

с.с.об.ур. < dt = u ! x ¡ ut = const; > du

: dt = 0 ! u = const:

(©(1) = x ¡ ut;

©(2) = u;

F(x ¡ ut; u) = 0 - общее решение. Задача Коши:

(

ut + uux = 0; ujt=0 = '(x)

имеет решение:

(8) u = '(x ¡ ut):

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

Так задача Коши: ut + ux = 0; ujt=0 = '(x) имеет гладкое решение при всех t > 0; x 2 R1, если функция '(x) непрерывно

Лекция №2, НГУ, ММФ, 2010

10

дифференцируема. Однако в случае задачи Коши (7) дело обстоит сложнее. Оказывается, далеко не всегда можно построить гладкое решение этой задачи при всех t > 0 (даже, если '(x) - гладкая функция). Итак, гладкое решение перестает существовать, как

 

 

 

1

 

 

-

 

 

t=

 

 

-

 

u=1,

x

 

 

 

 

 

-1

t

 

 

 

1

 

x=

(x0 )t+x0

 

 

x+t=1,

 

 

 

u=

 

 

 

-

 

 

 

1

 

 

 

1

 

x0

x

 

 

 

 

только характеристики пересеклись. Из (8) легко получаем

ux(t; x) = ux(t; x0 + '(x0)t) =

'0(x0)

:

1 + t'0(x0)

Следовательно, гладкое решение задачи Коши (7) существует при всех t > 0, если '0(x0) ¸ 0. Если же в некоторой области '0(x0) < 0, то гладкое решение задачи (7) существует при 0 < t < tk, где:

1

tk = sup j'0(x0)j

x0

(sup берется в той области, где '0(x0) · 0).

x0

При t ¸ tk гладкое решение перестает существовать. Явление неограниченного роста градиентов основных величин (например, ux) получило название градиентной катастрофы.

3.Уравнение Гамильтона-Якоби.

(9)ut + H(t; x; ru) = 0;

Соседние файлы в папке 2010-08-29