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

Pyrkova_--_Metody_reshenia_differentsialnykh_urav

.pdf
Скачиваний:
36
Добавлен:
27.03.2016
Размер:
2.03 Mб
Скачать

Частное решение, соответствующее поставленным на-

 

 

 

x

2

 

 

 

 

y = e

 

+1

 

 

 

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

2

. n

 

 

Пример

3.2.

(7-01)

Решить

задачу

Коши

2 y 2 y′′sin x 2 y 2 ycos x +(y)3 2 y(y)2 sin x = 0 ,

π

 

 

π

 

 

 

 

y

 

=1, y

= −1 .

 

 

 

2

 

 

2

 

 

 

 

d Нетрудно заметить, что исходное уравнение является однородным относительно y, y, y′′ , поэтому делаем замену

y′ = yz , при которой y′′ = yz + yz′ = yz z + yz′ = y(z 2 + z).

Исходное дифференциальное уравнение принимает вид: 2 y 2 y(z 2 + z)sin x 2 y 2 yz cos x +(yz)3 2 y(yz)2 sin x = 0 .

1.y = 0 - не удовлетворяет начальным условиям;

2.при y 0 имеем

2(z 2 + z)sin x 2z cos x + z 3 2z 2 sin x = 0 , или

 

2zsin x 2z cos x = −z 3

- это уравнение Бернулли.

 

Стандартная замена p =

 

1

, p′ = −

2

 

z.

 

 

z 2

z 3

 

 

Уравнение

2zsin x

+

2z cos x

 

=1

переходит

в

z 3

z 3

psin x +2 p cos x =1

 

 

 

 

 

- неоднородное

линейное уравнение

первого порядка по p, которое решаем методом вариации постоянной.

1) Сначала решаем однородное уравнение: psin x +2 p cos x = 0 - переменные разделяются.

 

p 0

 

 

dp

= −

2 cos xdx

дает ln

 

p

 

= −2 ln

 

sin x

 

+ln C0 и

 

 

 

 

 

 

 

 

 

 

p

sin x

p =

 

C

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin 2

x

 

 

 

 

 

 

 

 

 

 

 

44

2) Полагая C = C(x), подставляем

p =

C(x)

 

в линейное

 

 

 

 

неоднородное уравнение

 

 

 

 

 

 

 

 

sin 2 x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

 

sin x

2C cos x

sin x +2

 

 

 

C

cos x =1 дает

 

 

sin 2

 

 

 

 

 

 

 

 

sin 2 x

 

 

 

x

 

 

sin 3 x

 

 

 

 

 

 

 

 

 

 

 

 

 

C′ = sin x , т.е. dC = sin xdx и C(x)

= −cos x +C1

 

 

 

3) Используя найденное значение C(x), получаем

 

 

p =

cos x +C

1

, т.е. z 2 =

1

=

 

 

sin 2

x

 

.

 

 

 

 

 

 

 

sin 2 x

 

 

p

 

cos x +C1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

условия

 

 

π

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

=1, y

= −1 . Таким образом,

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

y

 

 

 

1

 

 

 

 

 

 

 

 

2 π

 

 

 

 

1

 

 

 

 

 

 

 

2

 

 

 

 

= −1 дает z

 

 

 

 

и C1 =1

 

 

z

 

 

 

=

 

 

=

 

 

 

 

 

=1 =

 

 

 

 

 

2

 

π

 

 

1

 

C1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

sin x

 

 

 

 

 

 

 

 

y

 

 

 

 

π

 

 

 

 

 

 

y

= ±

1cos x , т.к. y и

 

при

 

x = 2

 

имеют разные

 

 

 

 

 

 

знаки, то выбираем знак минус:

 

y=

sin x

 

 

- уравнение с

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

1cos x

 

 

 

разделяющимися

 

 

переменными,

т.е.

dy

=

 

sin dx

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

1cos x

 

ln y = −2 1cos x +C2 .

 

π

 

=1

, находим

Подставляя сюда начальное значение y

 

 

2

 

 

 

ln1 = −2

10 +C2 , откуда C2 = 2 и ln y = 2 2

1cos x ,

или, т.к.

y > 0 , получаем частное решение, соответствующее

45

поставленным начальным условиям: y = exp(2 2 1cos x ). o

3.4. Задачи для самостоятельного решения

Решить задачу Коши:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

57.

(7-01) 2 y 2 y′′sin x 2 y 2 ycos x + (y)3 2 y(y)2 sin x = 0 ,

 

π

 

 

 

 

 

 

 

 

π

 

 

 

= −1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

=1 , y

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1)(y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

58.

(7-02)

4 yy

′′

+ (4 y

2

6

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

)

 

 

 

+ (y )

 

= 0 , y(2)=1 , y (2)= −1 .

59.

(7-03) 2 y 2 y′′cos x + 2 y 2 ysin x + (y)3 2 y(y)2 cos x = 0 ,

 

π

 

 

 

 

 

 

 

 

 

π

 

 

 

= −

 

 

 

 

3

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y =1 , y

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

60.

(7-04)

4 yy

′′

+ (y

2

 

 

 

 

 

 

 

 

 

 

6

 

+

(y

2

=

 

 

 

 

 

1

 

 

 

 

9)(y

)

 

)

0 , y(0)=1, y (0)= − 2 .

61.

(6-11)

yy

+ 2x(y

2

 

= xyy

′′

,

 

 

 

 

 

 

 

 

(1)= −1 .

 

 

 

 

)

 

 

 

 

 

y(1)=1, y

 

 

62.

(6-12)

3yy

′′

3(y

 

2

 

=

 

4 y(y

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)

 

 

 

 

)

,

 

y(0)=1, y (0)= −1 .

63.

(6-13)

yy

xyy

′′

+ x(y

 

2

=

 

 

 

 

 

)

2

ln x ,

 

 

 

 

 

 

 

 

)

 

 

 

2x(y

 

y(1)=1, y (1)= −2 .

64.

(6-14)

′′

 

 

 

 

 

1)

2

+ (y

 

 

2

(y

 

 

 

 

 

 

 

3

 

 

 

 

y (y

 

 

 

)

 

 

 

1)= (y ) , y(0)= 2, y

 

(0)= 2 .

65.

(8-21) 2x(y

2

y

′′

 

+

(y

3

 

 

 

 

 

 

 

 

 

2

+ 2xy(y

2

, y(1)= y

 

 

 

 

)

 

 

)= y y

 

)

 

(1)=1 .

66.

(8-22)

′′

 

 

 

 

 

1)+ y

 

 

 

1)

2

 

 

 

 

 

 

2

,

y(0)=

 

 

y (y

 

(y

 

 

 

 

= (y )

 

2, y (0)= −2 .

67.

(8-23)

xyy

′′

= yy

+ x(y

 

2

+

 

4x

5

y

2

,

y(1)

=

 

 

 

 

 

 

 

)

 

 

 

 

1, y (1)=1 .

68.

(8-24)

2 yy

′′

 

 

 

 

 

 

 

 

 

 

 

 

2

 

(1 + 2 ln y), y(0)

 

 

 

 

 

 

ln y = (y )

 

 

 

= e, y (0)= e .

69.

(8-31)

yy

′′

 

 

= (y

 

2

(y

2

y

+1),

 

y(0)=

 

 

 

 

 

 

 

 

 

)

 

 

 

 

 

1, y

(0)= −3 .

 

 

70.(8-32) y′′y 2 x2 xy 2 y′− x2 (y)2 y e2x (y)3 = 0 , y(2)=1, y(2)= −2e .

71.(8-33) yy′′−(y)3 ctg y + cos2 y sin 4 y = 0 , y(0)= π4 , y(0)= 12 .

46

72.

(8-34)

yy

′′

 

yy

(y

 

2

(e

x

 

+1)= 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

)

 

 

 

y(1)= sh1, y (1)= −1 .

73.

 

 

2

 

 

 

′′

 

 

 

 

 

 

 

2

(y

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

(5-41)

x

 

yy

2x

 

 

 

 

+ xyy

=

0 ,

 

 

 

= − e .

 

 

 

 

 

)

 

 

 

y(e)=1, y

(e)

74.

(5-42)

xyy

′′

2x(y

 

2

 

+ yy

 

= 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

)

 

 

 

 

 

 

y(1)= y (1)=1 .

 

 

 

75.

(5-43)

y

2

y

′′

 

 

 

 

 

 

 

 

 

 

2

 

y

3

= 0

,

 

 

 

 

 

 

 

 

 

 

 

 

y(y )

 

 

 

 

 

y(0)=1, y

(0)=1 .

 

 

 

76.

(5-44)

xyy

′′

 

 

 

 

 

 

2

 

yy

= 0 ,

y(1)

 

 

 

 

 

 

 

 

 

 

+ x(y )

 

 

 

 

 

= 4, y

(1)=1 .

 

 

 

77.

(5-51)

3y

2

 

 

′′

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

= 4 y

3

,

 

 

 

 

 

 

 

 

 

 

 

 

y y

 

+ 2 y(y

)

 

 

 

y(1)= −1, y (1)=1.

78.

(5-52)

4 yy

′′

+ (y

2

 

+ y(y

 

4

 

=

0 ,

 

 

 

 

 

 

 

 

 

 

 

 

)

 

 

 

)

 

 

 

y(1)=1, y (1)= −1 .

79.

(5-53) 2(y +1)y

′′

+

 

 

 

 

 

 

2

 

= 2(y

+1),

y(2)=

0, y

 

 

 

 

 

 

(y )

 

 

 

(2)= −1 .

80.

 

 

 

′′

 

+3(y

2

 

 

 

 

 

 

 

 

2

(y

3

 

y(4)

 

 

 

 

1

 

 

 

 

(5-54)

yy

 

 

)

 

 

 

 

= y

 

 

)

,

 

= 2, y (4)=

4 .

 

 

 

81.

(5-61)

xyy

′′

+ 2x

3

 

 

 

 

 

2

yy

=

0 ,

y(1)= y

 

 

 

 

 

 

 

 

 

 

 

(y )

 

 

 

 

(1)=1 .

 

 

 

82.

(5-62)

y′′+ (y)2 tg y (y)4 sin 2 y = 0 , y(0)= π6 , y(0)= −

83.

(5-63)

xyy

′′

+ x

4

(y

 

2

 

+3yy

=

0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)

 

 

 

y(1)=1, y (1)= 2 .

84.

(5-64)

y

′′

(y

2

 

th y +

 

 

 

 

 

4

sh 2 y =

 

 

 

 

 

 

 

)

 

 

 

(y )

 

 

0 , y(0)= 0, y

(0)=1 .

2 .

3

3.6. Ответы:

57.y = exp{2 2 1 cos x}

58.y = 12 [(7 3x)23 +1]

59.y = exp{ 2 2 sin x}

60.y = 8 3 x 23 3

2

61. y =

2

 

x2 +1

62.y = 1 4 x 34

3

47

63.y =1 2 ln x

64.y =1 + e2x

65.y = e2 x 2

66.

y =

 

 

2e x

 

 

2e x

1

 

 

 

 

 

 

 

x6 1

 

67.

y = e

 

6

 

 

 

 

 

 

x

 

2

 

 

 

 

 

+1

 

68.

 

2

 

y = e

 

 

 

69.

y = 3 1 9x

 

 

 

 

 

 

 

 

 

x

 

70.

y = exp

(4 2x)e 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

71.y = arctg(e x )

72.y = e x +1 (e 1)2

e x 1 2e

73.y = ln1x

74.y = 1 1ln x

75.

x+

1

x2

 

 

 

y = e

2

 

 

76.

y = 2

x2 +3

77.

x 4 3

y =

 

 

 

 

 

3

 

 

 

2

78.y = 5 3x 3

2

79.y = x 4 2 1

2

48

80.y = 3 3x 4

81.y = 2x2

x2 +1

82.y = arcsin 1 x2

83.y = (2x3 1)23

x2

84.

 

x

2

 

x = sh y

y = ln x +

 

+1

 

 

 

 

 

 

49

§ 4. ЭЛЕМЕНТЫ ВАРИАЦИОННОГО ИСЧИСЛЕНИЯ

4.1 Основные понятия

Пусть M - некоторое множество функций. Функционалом J = J (y) называется переменная величина,

зависящая от функции y(x), если каждой функции y(x) M по некоторому правилу поставлено в соответствие число.

Множество M функций y(x), на котором определен функционал J (y), называется его областью определения.

В приложениях часто встречаются функционалы вида

b

 

(4.1)

J (y)= F(x, y(x), y (x))dx ,

a

где F(x, y, p) - заданная дважды непрерывно дифференци-

руемая функция для x [a, b] и (y, p) R(2y, p) - плоскости с декартовыми прямоугольными координатами y, p .

Будем обозначать через C1 [a, b] пространство всех непрерывно дифференцируемых на [a, b] функций с нормой

y

 

 

 

= max

 

y(x)

 

+ max

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

y (x)

 

 

 

 

 

x [a,b]

 

 

 

x [a,b]

 

 

 

 

 

 

 

 

 

 

 

 

 

Говорят, что функция yˆ(x) M , дает (слабый локальный)

минимум функционала J (y), если число ε > 0 такое, чтоy(x) M , для которого y(x)yˆ(x) < ε , выполнено неравенство J (y)J (yˆ ).

Говорят, что функция yˆ(x) M , дает (слабый локальный) максимум функционала J (y), если число ε > 0 такое, чтоy(x) M , для которого y(x)yˆ(x) < ε , выполнено неравенство J (y)J (yˆ ).

50

4.2 Простейшая вариационная задача

Простейшей вариационной задачей называется задача нахождения слабого локального экстремума функционала (4.1) в классе непрерывно дифференцируемых на [a, b] функций

y(x), удовлетворяющих граничным условиям

y(a)= A , y(b)= B .

(4.2)

Если функция yˆ(x)

является решением простейшей ва-

риационной задачи, то она на [a, b] необходимо удовлетворяет уравнению Эйлера

 

F

 

d

 

 

F

= 0

(4.3)

 

y

dx y

 

 

 

 

 

(здесь

 

d

 

 

- полная производная по x ).

 

 

dx

 

 

 

 

 

 

 

 

 

Экстремалью называется всякое решение уравнения Эй-

лера (4.3).

4.3 Задача со свободным концом (концами)

Пусть функционал (4.1) рассматривается при граничном

условии

 

y(a)= A .

(4.4)

Тогда говорят, что

x = a - закрепленный конец, x =b - сво-

бодный конец.

Задачей со свободным концом называется задача нахождения слабого локального экстремума функционала (4.1) в классе непрерывно дифференцируемых на [a, b] функций

y(x), удовлетворяющих условию (4.4).

Если функция yˆ(x) является решением задачи со свободным концом, то она на [a, b] необходимо удовлетворяет

уравнению Эйлера (4.3) и граничному условию при x =b вида

51

F

(b, yˆ(b), yˆ (b))= 0 .

(4.5)

y

 

 

Если функционал J (y) рассматривается при граничном

условии

 

y(b)= B ,

(4.6)

то x = a - свободный конец. Функция

yˆ(x), доставляющая

J (y) слабый локальный экстремум, должна удовлетворять

уравнению Эйлера (4.3), граничному условию (4.6) и граничному условию при x = a :

F

(a, yˆ(a), yˆ (a))= 0 .

(4.7)

y

 

 

Если граничных условий не ставится, то есть оба конца свободные, то yˆ(x) должна удовлетворять уравнению Эйлера (4.3) и граничным условиям (4.5), (4.7).

4.4. Решение уравнения Эйлера

Линейным дифференциальным уравнением n-го порядка с переменными коэффициентами называется дифференциальное уравнение вида

 

 

a0 (x)y (n) +a1 (x)y (n1)

+...+an1 (x)y′+an (x)y = f (x),

(4.8)

где x [a, b] - независимая переменная;

y(x) - искомая функ-

ция; f (x), a0 (x), a1

(x),K, an (x) - заданные на [a, b] функции,

причем x : x [a, b]функцияa0 (x)0 .

 

 

 

 

 

Уравнением Эйлера называется линейное дифференци-

альное уравнение

 

с переменными коэффициентами

вида

a

 

(x)= b

 

x nk , k =

 

,

где b , b ,K, b

 

- заданные

числа,

k

k

0, n

n

 

 

 

 

 

 

 

0 1

 

 

причем b0 0 :

 

 

 

+...+ bn1 xy′ + bn y = f (x).

 

 

 

b0 x n y (n)

+ b1 x n1 y (n1)

(4.9)

 

 

Заменой

x = et

(t = ln x) (4.9) сводится к линейному диф-

ференциальному уравнению с постоянными коэффициента-

52

ми. Действительно

d 2 y

 

d dy

= = e

t

d

t

 

=

 

 

 

 

 

 

e

 

dx 2

 

 

 

 

 

 

dx dx

 

 

dt

 

dydx = dydt dxdt = 1x dydt = et dydt ,

dy

 

2t

d 2 y

 

dy

 

 

 

= e

 

 

 

 

 

 

. Допус-

 

 

 

2

 

dt

 

 

 

dt

 

 

 

 

 

 

 

 

 

dt

 

тим, что k-ая производная имеет вид

d k y

=

dx k

 

 

ekt d k ky +α1

dt

 

1

d k y

 

=

 

 

 

 

 

+α

 

k

 

k

 

x

 

dt

 

 

 

 

 

 

d k 1 y

 

 

 

 

dy

 

 

 

 

 

 

 

 

 

 

 

 

 

+K+α

k 1

 

 

 

=

 

 

 

 

 

 

 

 

 

k 1

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

d k 1 y

 

 

 

 

dy

 

 

 

 

 

 

 

1

 

 

 

 

 

+K+

α

k 1

 

 

 

, где

α

,α

2

,K,α

k 1

-

 

 

k

1

 

 

dt

 

 

dt

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

постоянные.

 

Тогда

 

(k+1)-ая

производная

будет равна

 

d k +1 y

 

 

d

 

d k

y

 

 

 

t d

 

d k y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

=

e

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

k +1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

dx

 

 

dx

 

dx

k

 

 

 

 

 

dt

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

k +1

y

+(α

 

k )

d

k

y

 

 

 

 

 

 

 

 

dy

 

 

= e(k +1)t

 

 

 

1

 

+K−kα

k

1

 

=

 

 

k|+1

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

dt

k

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

d

k +1

y

+(α

 

k )

d

k

 

y

 

 

 

 

 

 

 

 

dy

 

 

 

=

 

 

 

 

 

 

 

 

1

 

 

+K−kα

k 1

 

.

 

 

 

k +1

 

 

k|+1

 

 

k

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

Так как в преобразованном уравнении, в случае отсутствия кратных корней характеристического уравнения, реше-

ния имеют вид y = eλt , следовательно, в исходном уравнении

они имеют вид y = xλ . Поэтому можно непосредственно подставить его в уравнение Эйлера (4.9). Поскольку

xk d k xλ = λ(λ 1)L(λ k +1) при k λ , то характеристиче- dxk

ское уравнение имеет вид

b0 λ(λ 1)L(λ n +1)+K+bn2 λ(λ 1)+bn1λ +bn = 0 . (4.10)

Каждому простому корню λ уравнения (4.10) соответствует частное решение однородного уравнения Эйлера xλ ;

53

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