Материал: sb000238

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

ность x2 y2 x2 y2

(рис.

1.1). И в этом примере через каждую точку

 

0

0

 

 

 

 

 

 

плоскости, кроме начала координат, проходит одно-единственное решение.

 

Нетрудно написать соответствующее дифференциальное уравнение:

 

 

 

 

y x , или

xdx ydy 0 .

 

 

 

 

 

 

 

y

 

 

 

 

Пример 3. Данный пример дополняет предыдущий (рис. 1.2). В каждой

 

точке вектор совпадает с вектором, начинающимся в этой точке и заканчи-

 

вающимся в начале координат. Решение здесь совсем очевидно это прямые

 

вида y kx , точнее ax by 0. Как и раньше, через каждую точку, кроме на-

 

чала координат, проходит одно решение. Однако в особой точке (в начале

 

координат) картина совсем иная – через нее проходят все решения. Соответ-

 

ствующее дифференциальное уравнение имеет вид y y , или

ydx xdy 0.

 

 

 

 

 

 

x

 

 

 

2

 

 

 

 

2

 

 

 

Рис. 1.1. Геометрический смысл и решения

-2-2 Рис. 1.2. Геометрический смысл и решения

2

-2-2

 

 

x

2

Пример 4. Бывают и более слож

 

уравнения

y

 

 

 

y

y

 

y

 

уравнения

x

 

 

 

 

 

 

Пример 4. Бывают и более сложные примеры. На рис. 1.3 изображен

 

геометрический смысл и решения уравнения y x y2. Известен факт, что

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

ся кривая x y2 (для этого достаточно просто приравнять производную к нулю). Можно также проверить, что кривая перегибов в данном случае (т. е.

6

2

 

1.5

 

1

 

0.5

 

0

 

-0.5

 

-1

 

-1.5

 

-2

Рис. 1.3. Геометрический смысл и решения уравнения y x y2

-2

2

кривая, каждая точка которой является точкой перегиба какого-либо из ре-

шений) это кривая x y2 21y (сначала нужно выразить вторую произ-

водную через x и y , а затем приравнять ее к нулю).

1.3.Условия единственности решения. Задача Коши

Вбольшинстве ситуаций, возникающих при решении дифференциальных уравнений, важно иметь дополнительные условия, гарантирующие единственность решения. Приведенные в 1.2 примеры позволяют предположить, что для этого достаточно знать, что решение проходит через заданную точку области, где определена функция f (x, y) . Эта гипотеза «почти верна», точ-

нее справедлива следующее утверждение [2], называемое теоремой сущест-

вования и единственности решения дифференциального уравнения: если в некоторой области плоскости функции f (x, y) и f y (x, y) непрерывны, то

через каждую точку области проходит одно решение дифференциального уравнения y f (x, y).

Ситуация, когда решение дифференциального уравнения единственно, очень важна, и для нее есть специальное название задача Коши:

y f (x, y), y(x0) y0,

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

7

Теорема единственности и постановка задачи Коши переносятся на дифференциальные уравнения высших порядков:

y(n) f x, y, y ,

, y(n 1) , y(x0) y0,

, y(n 1)(x0) yn 1

и на системы дифференциальных уравнений:

 

y

f (x, y ,

, y

n

), y (x ) y

,

 

1

1

 

1

0

1, 0

 

 

 

 

 

...;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

f (x, y ,

, y

n

), y

n

(x ) y

n, 0

.

 

n

1

 

 

0

 

 

Условия теоремы можно ослабить, но это требует введения дополнительных определений. Доказательство теоремы основано на построении последовательности приближений к решению. Доказательство сходимости такой последовательности трудное и неконструктивное. Устанавливается только факт сходимости, алгоритма поиска решения не предлагается. Однако методы построения приближенного решения служат базой для многочисленных приемов численного решения задачи Коши. Далее рассмотрим наиболее важный алгоритм [3].

 

 

 

1.4. Метод Эйлера (метод касательных)

 

 

Задача Коши

y f (x, y), y(x0) y0 дает возможность описать каса-

 

тельную к решению y y0 f (x0, y0)(x x0). Возьмем число x1, находящее-

 

ся от x0 на малом расстоянии h , и поставим новую задачу Коши:

 

 

 

 

 

 

 

 

 

y f (x, y), y(x1) y1,

 

 

 

 

 

 

 

 

 

где

1.2

 

 

 

 

 

 

 

 

y1 y0 f (x0, y0)h ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т. е. заменим решение исходной задачи

 

 

 

 

 

 

 

 

 

Коши касательной. Продолжая этот

 

 

 

 

 

 

 

 

 

процесс, получим ломаную Эйлера

 

(x0, y0)

 

 

 

 

 

 

приближенное решение задачи Коши.

 

 

 

 

 

 

 

Можно доказать, что ломаные сходятся

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

к решению дифференциального урав-

0

 

0.2

0

0.2

0.4

0.6

0.8

1

 

-1

-0.8

-0.6 Рис-0.4. 1.4.- Ломаные Эйлера

 

 

нения (при стремлении длины звеньев

 

для задачи Коши

y x

,

y(x0) y0

 

 

 

ломаной к нулю). Сказанное наглядно

 

 

 

 

y

 

 

 

 

продемонстрировано на рис. 1.4 (разби-

 

 

при различных h

 

 

 

 

 

 

 

 

 

 

 

 

рается уравнение из примера 2).

8

Стоит отметить, что в численных методах метод Эйлера не применяют изза плохой сходимости. Однако небольшая модификация метода, позволяющая использовать для построения аналога касательной несколько точек, исправляет этот недостаток. Такие методы называются методами Рунге Кутты.

Явное (аналитическое) решение дифференциального уравнения – ситуация довольно редкая (основным типам уравнений первого порядка, интегрируемым в квадратурах, посвящены 2.1 2.4). Объемный справочник [4] дает правильное представление об этом. Из многочисленных примеров в нем видно, что стоит немного изменить уравнение и возможность его явного решения исчезает. При построении моделей тех или иных процессов, основанных на дифференциальных уравнениях, очень важно на начальном этапе так упростить модель, чтобы описывающее ее дифференциальное уравнение было достаточно простым. К таковым в первую очередь принадлежат линейные дифференциальные уравнения с постоянными коэффициентами. Понимание того, как можно получить решение такого уравнения, необходимо для правильного использования таких моделей.

С точки зрения математики эти уравнения хороши тем, что сравнительно легко переписываются в эквивалентные линейные алгебраические решения. Решению таких дифференциальных уравнений посвящены 3.1 3.3.

2. ИНТЕГРИРУЕМЫЕ УРАВНЕНИЯ ПЕРВОГО ПОРЯДКА

2.1. Уравнения с разделяющимися переменными

Уравнениями с разделяющимися переменными называются уравнения,

которые при помощи тождественных преобразований могут быть приведены к виду

g(x)dx h(y)dy

(отсюда и происходит название переменные как бы «разделились»). Теорема об интеграле с переменным верхним пределом гарантирует наличие первообразных G(x) и H (y) , для которых выполняется равенство

G(x) H (y) C.

Теперь из теоремы о неявной функции вытекает существование функции y y(x), обращающей это равенство в тождество. (На практике зачастую

данный шаг либо затруднителен, либо не приводит к существенному упрощению записи ответа. В этом случае его опускают.)

9

Продемонстрируем приведенный алгоритм на примерах из предыдущего параграфа:

1. y xy .

Распишем производную через дифференциалы. Получаем соотношение dydx xy .

Далее перенесем все y в левую часть, а все x в правую. Заметим, что при

делении на x решения не теряются (кривая x 0 не являлась решением исходного уравнения, поскольку dx стоял в знаменателе). Следовательно,

ydy xdx .

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

y2

 

x2

C , или

x2 y2 C .

2

 

2

 

 

Произвольная постоянная была переобозначена. Это допустимо, поскольку 2C, как и C, принимает любые значения.

2. y xy .

Распишем производную через дифференциалы. Получаем соотношение dydx xy .

Далее перенесем все y в левую часть, а все x в правую. Заметим, что при делении на y могло потеряться решение y 0. Подставив его в исходное уравнение, получим 0 0, т. е. верное соотношение. Таким образом, y 0 действительно является решением исходного уравнения. Следовательно,

dy dx ,

y xy 0,

где квадратными скобками обозначено объединение множеств решений. Вычисляя интегралы, получаем:

ln | y | ln | x | C,

y 0,

10