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

УМК

.PDF
Скачиваний:
91
Добавлен:
17.03.2015
Размер:
3.93 Mб
Скачать

 

 

 

1

2π

 

αn

=

 

f (ψ)cos nψ dψ ,

(1.279)

 

 

 

 

 

π 0

 

 

 

1

 

2π

 

βn

=

 

f (ψ)sin nψ dψ.

 

π

 

 

 

0

 

Сопоставляя (1.277) и (1.278), заключаем, что для выполнения граничного

условия (1.277) нужно положить A 0 ,

a n A n , a n Bn равными коэффициентам

Фурье:

 

α0

 

 

 

 

 

 

A 0

=

, a n A n = αn , a n Bn = βn .

 

 

 

2

 

 

 

 

 

 

Следовательно, для внутренней задачи

 

αn

 

= βn ,

A 0 =

α0

, A n

=

, Bn

 

 

 

 

2

 

 

a n

a n

для внешней задачи

 

α0

 

 

 

 

 

 

A 0

=

, A n = αn a n , Bn = βn a n .

 

2

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

u(ρ, ϕ) =

α

0

ρ n

(α

 

 

 

 

 

+

 

 

n

cos nϕ + β

n

2

a

 

n 1

 

 

 

 

 

 

=

 

 

 

 

 

а решение внешней задачи Дирихле в виде ряда

u(ρ, ϕ) =

α

0

a n

(αn cos nϕ + βn

 

+

 

 

2

 

 

n=1

ρ

 

sin nϕ),

(1.280)

sin nϕ),

(1.281)

где α0 , αn , βn определяются как коэффициенты ряда Фурье для функции f (ϕ) по формулам (1.279).

Установлено, что так построенные функции (1.280) и (1.281) удовлетворяют всем условиям задачи.

Преобразуем формулы (1.280), (1.281) к более простому виду. Рассмотрим подробно внутреннюю задачу, а для внешней задачи получим результат по аналогии. Преобразуем ряд (1.280), подставив выражения для коэффициентов Фурье (1.279) и произведя затем перестановку порядка суммирования и интегрирования

 

 

 

 

 

 

 

1

 

2π

 

 

 

 

 

 

 

 

 

 

 

 

u(ρ, ϕ) =

 

 

f (ψ)dψ +

 

 

 

 

 

 

 

 

2π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ

n

 

 

 

 

2π

 

 

 

 

1

2π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

f (ψ)cos nψ dψ cos nϕ +

 

 

f (ψ)sin nψ dψ sin nϕ

 

=

 

 

 

 

 

 

 

n =1 a

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

π 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

1

 

2π

 

 

1

 

 

 

 

 

ρ

 

 

 

 

 

 

 

 

=

 

 

 

f (ψ)

 

 

+

 

 

(cos nψ cos nϕ + sin nψ sin nϕ) dψ =

 

 

 

 

π

 

 

 

 

 

 

 

 

0

 

 

 

2

 

 

 

n =1 a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

1

 

2π

 

 

1

 

 

 

 

 

ρ

 

 

 

 

 

 

 

 

=

 

 

 

f (ψ)

 

 

+

 

 

cos n(ϕ − ψ) dψ.

 

 

 

 

(1.282)

 

π

 

 

 

 

 

 

 

 

 

0

 

 

 

2

 

 

 

n =1 a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Произведем тождественные преобразования в фигурных скобках, обозначив

для краткости ρ = t ,

ϕ − ψ = ω ; при этом воспользуемся формулой

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

(eiθ + eiθ ),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos θ =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

cos n(ϕ − ψ)

=

 

 

+ t n cos nω =

 

 

 

 

+

 

 

 

 

t n (einω + e

inω )=

2

 

2

 

 

 

 

 

 

 

 

 

 

n=1 a

 

 

 

 

 

 

 

 

 

 

 

 

n =1

 

 

 

 

 

 

 

 

2 2 n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

n

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

1 +

(teiω )

+

(te

iω )

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

n =1

 

 

n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ряд

 

 

(teiω ) = q n

=

 

 

 

 

 

является

 

 

 

 

бесконечной

 

 

 

 

геометрической

 

 

1

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n =1

 

 

n=1

 

 

 

ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ

 

прогрессией со знаменателем q =

e

i(ϕ−ψ)

,

модуль которого

 

q

 

=

< 1. В силу

 

 

a

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

этого имеем следующее представление:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

n

 

 

1

 

 

 

 

 

teiω

 

 

 

 

teiω

 

 

 

 

1

 

 

 

 

 

 

1 t 2

 

 

 

 

 

 

+

t

 

cos nω =

 

 

1 +

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

iω

 

 

 

 

 

 

 

 

 

 

 

 

 

iω

 

iω

 

 

2

2

 

 

n=1

 

 

 

2

 

 

1 te

 

 

1 te

iω

 

 

2

 

 

1 te

te

+ t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

ρ

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1 t 2

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

2

1 2t cos ω + t 2

2

 

 

 

 

 

ρ

 

 

 

 

 

 

 

+

ρ

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 2 cos(ϕ − ψ)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

=

 

1

 

 

 

 

a 2 − ρ2

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

a 2 2aρcos(ϕ − ψ) + ρ2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

1

2π

 

a 2 − ρ2

 

 

 

 

 

u(ρ, ϕ) =

 

f (ψ)

 

 

 

 

 

 

 

dψ .

(1.283)

2π

a 2

2a

ρ

cos(

)

+ ρ

2

 

 

0

 

ϕ − ψ

 

 

 

 

Эта формула называется формулой Пуассона, а интеграл справа – интегралом Пуассона. Формула (1.283) позволяет записать искомое решение

внутренней задачи Дирихле в форме интеграла, зависящего от параметров ϕ и ρ. Он существует для всех значений ϕ и ρ, 0 ≤ ρ < a и удовлетворяет уравнению (1.263), а также граничному условию (1.264). Заметим, однако, что интеграл (1.283) теряет смысл при ρ = a . Когда говорят, что функция (1.283)

удовлетворяет

 

граничному

условию,

то под этим подразумевают, что

lim u(ρ, ϕ) = f (ϕ). Поэтому решение внутренней задачи записывается так:

ρ→a0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2π

 

a

2

− ρ

2

 

 

 

 

 

 

f (ψ)

 

dψ

 

ρ < a,

 

u(ρ, ϕ) =

 

 

 

при

 

 

ρ2

2aρcos(ϕ − ψ) + a 2

 

 

2π 0

 

при

ρ = a.

(1.284)

 

f (ϕ)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

По аналогии решение внешней задачи имеет вид

 

 

 

 

 

1

2π

 

ρ

2

a

2

 

 

 

 

u(ρ, ϕ) =

 

f (ψ)

 

 

 

dψ

при

ρ > a,

 

 

ρ2

2aρcos(ϕ − ψ) + a 2

 

 

2π 0

 

при

ρ = a.

(1.285)

 

f (ϕ)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Замечание. Интеграл Пуассона дает решение плоской задачи Дирихле для круга. Но он может служить решением и пространственной задачи Дирихле в случае цилиндрической симметрии области, т.е. когда область есть бесконечный цилиндр, а искомая функция u = u(x, y, z), например, температура или электрический потенциал, не зависит от z . Тогда температура или потенциал в любом сечении, перпендикулярном к оси Oz, зависит только от ρ и ϕ. Следовательно, здесь имеет место плоская задача Дирихле.

1.38.2. Решение задачи Дирихле для прямоугольника

y

ϕ(x)

b

 

 

 

F (y)

 

 

 

ψ (y)

 

 

 

0

 

 

 

 

f (x)

a

 

x

 

Поставим

задачу

Дирихле для

прямоугольника

0 x a , 0 ≤ y ≤ b

(рис.

1.12)

при

произвольных

граничных условиях: найти функцию

u = u(x, y),

удовлетворяющую

уравнению Лапласа

 

Рис. 1.12

2 u

+

2 u

= 0 ,

(1.286)

x 2

y2

 

 

 

0 < x < a , 0 < y < b

и двум парам граничных условий

u(0, y) = ψ(y),

u(a, y) = F(y),

0 ≤ y ≤ b ,

(1.287)

u(x,0) = f (x),

u(x, b) = ϕ(y),

0 ≤ x ≤ a ,

(1.288)

причем

 

 

 

ψ(b) = f (0), ψ(b) = ϕ(0),

 

F(0) = f (a ), F(b) = ϕ(a).

(1.289)

Последние условия (1.289) обеспечивают непрерывность граничной функции в вершинах прямоугольника.

Не нарушая общности, можно предположить, что все ее значения в (1.289) равны нулю, т.к. к этому случаю можно привести любую задачу. Для этого надо представить решение в виде

 

 

 

u(x, y) = v(x, y) + ω(x, y),

где

 

 

ω(x, y) = A + Bx + Cy + Dxy

при

A = f (0), B =

1

[f (a)f (0)],

C =

1

[ψ(b) − ψ(0)],

 

 

 

 

a

 

b

D = 1 {[ϕ(a) − ϕ(0)][f (a) f (0)]}.

ab

Для функции v(x, y) получим задачу Дирихле при нулевых значениях в вершинах прямоугольника. Итак, будем решать задачу (1.286) – (1.289), предполагая, что все значения в (1.289) есть нули.

Представим решение в виде

u(x, y) = u1 (x, y) + u 2 (x, y),

(1.290)

где u1 (x, y) и u 2 (x, y) удовлетворяют уравнению Лапласа

и следующим

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

 

 

u1 (0, y) = ψ(y),

u1 (a, y) = F(y),

(1.291)

u1 (x,0) = 0 ,

u1 (x, b) = 0 ,

(1.292)

u 2 (0, y) = 0 ,

u 2 (a, y) = 0 ,

(1.293)

u 2 (x,0) = f (x),

u 2 (x, b) = ϕ(x).

(1.294)

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

Функции u1 (x, y) и u 2 (x, y) можно найти методом разделения переменных. Проведем решение для функции u 2 (x, y).

Ищем частные решения уравнения (1.286), удовлетворяющие условиям

(1.293), в виде

 

 

 

 

 

u 2 (x, y) = X(x)Y(y).

(1.295)

Подстановка (1.295) в (1.286) дает

 

 

X′′(x)

 

Y′′(y)

 

 

 

 

= −

 

= −λ ,

λ = const ,

 

X(x )

Y(y)

откуда имеем

 

 

 

 

 

 

 

Y′′ − λY = 0 ,

(1.296)

 

 

 

X′′ + λX = 0 .

(1.297)

Подстановка граничных условий (1.293) по переменной x в (1.295) приводит к соотношениям

X(0) = 0 , X(a ) = 0 .

(1.298)

Для определения функции X(x ) имеем задачу (1.297), (1.298) – задачу Штурма-Лиувилля, которая, как нам известно (см. п. 1.33.1А), имеет решение вида

 

nπ

2

(x) = sin

nπ

 

λn

=

 

 

, X n

 

x ,

n = 1, 2,...

a

 

 

 

 

 

 

a

 

Общее решение уравнения (1.296) при λ = λn можно записать в виде

Yn (y) = A n ch nπ y + Bn sh nπ y . a a

Перемножив теперь X n (x) и Yn (y), находим совокупность функций u 2n (x, y), удовлетворяющих уравнению (1.286) и граничным условиям (1.293):

 

 

nπ

nπ

 

nπx

 

u 2n

(x, y) = A n ch

 

y + Bn sh

 

y sin

 

,

n = 1, 2,...

 

a

 

 

 

a

 

a

 

Общее решение задачи ищем, как обычно, в виде ряда

в ряд Фурье по синусам:

 

 

 

 

nπ

 

 

 

nπ

nπx

 

u

2

(x, y) = A

n

ch

 

y + B

n

sh

 

y sin

 

.

(1.299)

a

a

 

 

n 1

 

 

 

 

a

 

 

 

=

 

 

 

 

 

 

 

 

 

 

Удовлетворим теперь граничным условиям (1.294), которые содержат заданные произвольные функции

 

 

 

u 2 (x,0) =

 

 

 

nπ

 

= f (x),

 

 

 

 

 

 

 

 

A n

sin

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n=1

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

nπ

 

 

 

 

nπ

nπx

 

 

 

 

 

u

2

(x, b) = A

n

ch

 

 

b + B

n

sh

 

b sin

 

= ϕ(x).

 

 

 

 

a

a

 

 

 

 

n 1

 

a

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

nπ

 

nπ

 

Отсюда видно, что

постоянные

множители

A n и

A n ch

b + Bn sh

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

a

должны являться коэффициентами разложения функций соответственно f (x ) и

ϕ(x)

A n

=

 

2

af (x)sin

nπ

xdx = f n ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a 0

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

A n ch

nπ

b + Bn sh

nπ

b =

2

a

ϕ(x)sin

nπ

xdx = ϕn .

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

a

 

 

 

a 0

 

 

 

 

 

a

Откуда определяем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

nπb

 

 

= f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

,

B

n

=

 

 

 

 

 

ϕ

n

f

n

ch

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

n

 

 

 

 

sh

nπb

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

Подстановка этих значений постоянных A n и Bn в (1.299) приводит функцию u 2 (x, y) после несложных преобразований к виду

 

 

 

 

 

 

 

nπ(b y)

 

 

 

nπy sin

nπx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

u

2

(x, y) =

f

n

sh

 

 

 

+ ϕ

n

sh

 

 

 

 

 

 

 

 

.

(1.300)

 

a

a

 

 

nπb

 

n 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Совершенно аналогично найдем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

nπ(a x)

 

 

 

nπx sin

nπy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

u

 

(x, y) = ψ

 

 

sh

 

+ F sh

 

 

 

 

 

 

 

 

 

,

(1.301)

 

 

 

b

b

 

 

 

 

 

nπa

 

 

1

n 1

 

 

n

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

где

ψn

=

 

2

b

ψ(y)sin

nπy

dy ,

 

 

 

 

 

 

 

b 0

 

 

b

F

=

2

b F(y)sin

nπy

dy .

 

 

n

 

 

 

 

 

 

b

 

 

b 0

 

 

Окончательное решение задачи: u(x, y) = u1 (x, y) + u 2 (x, y).

1.39 ЗАКЛЮЧЕНИЕ

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

– дуги окружностей с центрами в полюсе и отрезки лучей, выходящих из полюса, и.т.д.

Эти обстоятельства сильно ограничивают применимость метода Фурье. Есть множество задач, где разделение переменных не удается либо потому что граница тела достаточно сложна, либо потому что коэффициенты в исходном дифференциальном уравнении не являются постоянными, либо потому что уравнение не является линейным. Поэтому разработаны и другие методы решения задач математической физики, такие, например, как метод конечных разностей, метод интегральных уравнений, метод Ритца-Галеркина, метод физического моделирования, рассмотрение которых выходит за рамки настоящего курса.

УЧЕБНО - МЕТОДИЧЕСКИЙ КОМПЛЕКС

РАЗДЕЛ 9 «ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ»

2. Методические указания для студентов

2.1 ОБЫКНОВЕННЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ ПЕРВОГО ПОРЯДКА. ОСНОВНЫЕ ПОНЯТИЯ

Уравнение, связывающее независимую переменную, функцию и ее про-

изводные называется обыкновенным дифференциальным уравнением.

F(x, y, y,Ky(n) ) = 0

Наивысший порядок производной, входящей в уравнение,

называется

порядком дифференциального уравнения.

 

Например:

 

 

1)

y′ + xy x2

= 0 обыкновенное дифференциальное уравнение пер-

 

вого порядка;

 

 

2)

y′′ − y′ = 1

обыкновенное дифференциальное уравнение 2-го по-

 

рядка;

 

 

3)

y2 y′ ′ +′ x5

= 0 обыкновенное дифференциальное

уравнение

третьего порядка.

Рассмотрим обыкновенные дифференциальные уравнения первого порядка, то есть уравнения F(x, y, y) = 0 или в разрешенном относительно yви-

де, y′ = f(x, y) .

Решением дифференциального уравнения F(x, y, y) = 0 называется та-

кая дифференцируемая функция y = ϕ(x), которая при подстановке в уравне-

ние вместо неизвестной функции обращает его в тождество.

Процесс нахождения решения дифференциального уравнения называется

интегрированием дифференциального уравнения.

Общим решением дифференциального уравнения первого порядка y′ = f(x, y) в области D называется функция y = ϕ(x, C), обладающая сле-

дующими свойствами: 1) она является решением данного уравнения при любых значениях произвольной постоянной, принадлежащих некоторому множеству;

2) для любого начального условия y(x0 ) = y0 , такого, что (x0 ; y0 ) D , суще-

ствует единственное значение C = C0 , при котором решение y = ϕ(x; C0 )

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

С геометрической точки зрения общему решению y = ϕ(x, C) на плос-

кости XOY соответствует семейство интегральных кривых, зависящее от одного параметра - произвольной постоянной C.

Равенство Φ(x, y, C) = 0, неявно задающее общее решение называется

общим интегралом уравнения y′ = f (x, y).

Всякое

решение y = ϕ(x; C0 ) , полученное из общего решения

y = ϕ(x; C)

при конкретном значении C = C0 , называется частным реше-

нием.

Частному

решению

удовлетворяющему

начальному

условию

y |x =x0 = y0 ,

на плоскости X0Y соответствует линия

y0

= y0 (x 0 ), проходя-

щая через точку M 0 (x0 ; y0 )Аналогично определяются частные интегралы.

Задача,

в

которой

требуется найти

частное

решение

уравнения

y′ = f(x, y), удовлетворяющее начальному условию

y(x0 ) = y0

называется

задачей Коши.

 

 

 

 

 

 

 

Теорема Коши.

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

непрерывна в некоторой области

D плоскости Х0Y и имеет в этой области непрерывную частную производ-

ную по y , f y(x, y), то,

какова бы ни была точка (x 0 , y0 ) области D , сущест-

вует, и притом единственное, решение y = ϕ(x) уравнения y′ = f (x, y), опре-

деленное в некотором интервале, содержащем точку x 0 .

Особым решением называется такое решение, во всех точках которого условие единственности не выполняется, т.е. в любой окрестности каждой точки (x; y) особого решения существует по крайней мере две интегральные кри-

вые, проходящие через эту точку.

Особые решения не получаются из общего решения дифференциального

уравнения ни при каких значениях произвольной постоянной C ( в том числе и при C = ±∞ ).

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

Например, общее решение уравнения y′ = ±1 y2 записывается в ви-

де y = sin(x + C), − ∞ < C < ∞ . Это семейство имеет две огибающие: y = 1

и y = −1, которые будут особыми решениями данного уравнения.

2.2 ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ С РАЗДЕЛЕННЫМИ И РАЗДЕЛЯЮЩИМИСЯ ПЕРЕМЕННЫМИ

Дифференциальное уравнение y′ = f(x, y) называется уравнением с разделяющимися переменными, если его правая часть есть произведение двух функций, одна из которых зависит от y , а другая от x . y= f (x) g(y).

Предположим, что функции f (x) и g(y) непрерывны на интервале a < x < b, c < y < d и что g(y)≠ 0 .

Умножая обе части уравнения y= f (x) g(y) на dx и деля на g(y), за-

пишем его в виде: dy( ) = f (x)dx . g y

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