Тысяча полезных мелочей    

ПРАКТИЧЕСКОЕ ПРИМЕНЕНИЕ ТЕОРЕМ СУЩЕСТВОВАНИЯ И ЕДИНСТВЕННОСТИ

ПРАКТИЧЕСКОЕ ПРИМЕНЕНИЕ ТЕОРЕМ СУЩЕСТВОВАНИЯ И ЕДИНСТВЕННОСТИ

Источник - Запись к врачу без талона
Ахметова Ю.А. Афонасенков О.В. Статья в формате PDF 654 KB

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

Теорема существования. Если в уравнении

(1) функция f определена и непрерывна в некоторой ограниченной области D плоскости (x, y), то для любой точки (x0, y0) € D существует решение y(x) начальной задачи

, y(x0) = y0, (2)

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

Теорема существования и единственности. Если в уравнении  функция f определена и непрерывна в некоторой ограниченной области D плоскости (x,y), причем она удовлетворяет в области D условию Липшица по переменной y, т.е.

где L - положительная постоянная, то для любой точки (x0, y0)€ D существует единственное решение y(x) начальной задачи (2) определенное на интервале, содержащем точку x0.

Теорема о продолжении. При выполнении условий теоремы существования или теоремы существования и единственности всякое решение уравнения с начальными данными (x0, y0) Î D может быть продолжено до точки, сколь угодно близкой к границе области D. При этом в первом случае продолжение, вообще говоря, будет не обязательно единственным, во втором же случае оно единственно.

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

Пример 1. Требуется, используя численный метод интегрирования Эйлера с итерационной схемой

c шагом h =0,1, решить начальную задачу

 y(-1) = 0,21 (3)

на отрезке [-1, 3].

Решение. (c помощью пакета Mathcad)

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

Таким образом, в соответствии с теоремой существования существует решение y(x) начальной задачи (3), определенное на некотором интервале, содержащем точку x0 = -1, и это решение по теореме о продолжении может быть продолжено до значения, близкого к значению y(x) = 0. В результате численного интегрирования получаем решение начальной задачи (3) на некотором интервале (a, b), где a < -1; 1,3 < b < 1,32. Однако, учитывая, что это уравнение с разделяющимися переменными, можно аналитически найти частное решение, удовлетворяющее начальной задаче (3)

Интегрируя, получаем, что

Отсюда следует, что решение начальной задачи (3) существует только для

Оказывается, обращение к теореме существования (и к теореме о продолжении) позволило отсечь отрезок (приблизительно [1,315; 5]), на котором решение исходной начальной задачи (3) заведомо не существует. Одно же только численное интегрирование приводит к ошибочному результату. Дело здесь в том, что при приближении решения y = y(x) к оси Ox угол наклона кривой приближается к 90°. Поэтому пока аргумент x изменяется на величину 0,1 значение y успевает перескочить ось Ox, и мы попадаем на интегральную кривую, отличную от исходной. А это происходит потому, что метод Эйлера учитывает угол наклона только в текущей точке.

Пример 2. Используя метод Эйлера, а затем метод Эйлера-Коши, с шагом h = 0,1 и итерационной схемой

,

где

,

решить начальную задачу

 y(-1) = -1, (4)

на промежутке [-1, 1].

Решение. На базе Mathcad методом Эйлера, а затем методом Эйлера-Коши будем иметь:

 

 

Рис. 2

Получили чертеж (рис. 3) отличный от чертежа, изображенного на (рис. 2). Чтобы лучше разобраться в причине расхождения в результатах, проинтегрируем исходную начальную задачу. Разделяя переменные, имеем

или, окончательно, .

Становится понятно, что решение по методу Эйлера приближает функцию y1(x) = x3, а по улучшенному методу Эйлера - функцию

При этом как, y1 так и y2 являются решениями начальной задачи (4), а значит, для рассматриваемой на промежутке [-1; 1] начальной задачи имеет место неединственность.

Обращаясь теперь к теореме существования и единственности, отметим, что, так как функция f, заданная равенством, непрерывна во всей плоскости (x, y), то из теоремы существования следует, что существует решение начальной задачи (4), определенное на некотором промежутке, содержащем точку x0 = -1, и это решение по теореме о продолжении может быть продолжено на любой промежуток. Далее, поскольку

,

то функция

удовлетворяет условию Липшица по переменной y в любой области, не содержащей точки оси Ox. Если же область содержит точку оси Ox, то нетрудно показать, что в ней указанная выше функция условию Липшица не удовлетворяет. Поэтому из теоремы существования и единственности (и теоремы о продолжении) следует, что в данном случае решение начальной задачи может быть продолжено единственным образом, по крайней мере, до оси Ox. Но поскольку прямая y = 0 является особой интегральной прямой для дифференциального уравнения

,

то, как только y станет равным нулю решение начальной задачи (4) не может быть единственным образом продолжено за точку O(0, 0).

Рис. 3

Итак, обращение в данном случае к теореме существования и единственности (и теореме о продолжении) позволило разобраться в результатах численного интегрирования. Если речь идет о единственном на промежутке [-1; 1] решении начальной задачи (4), то оно существует и определено лишь на отрезке [-1; 0]. В общем же случае таких решений несколько.

Список литературы

1. Roberts C.E. Jr. Why teach existence and uniqueness theorems in the first course in ordinary differential equations? // Int. J. Math. Educ. Sci. Technol. - 1976. - Vol. 7, № 1. - P. 41-44.



КОМПОНЕНТНОЕ РАВНОВЕСИЕ И УСТОЙЧИВОСТЬ ДРЕВОСТОЯ

КОМПОНЕНТНОЕ РАВНОВЕСИЕ И УСТОЙЧИВОСТЬ ДРЕВОСТОЯ Статья в формате PDF 382 KB...

31 10 2026 18:38:55

Использование углубленных интеграционных методов исследования в изучении роли энергии геннообусловленных патологий

Использование углубленных интеграционных методов исследования в изучении роли энергии геннообусловленных патологий Статья посвящена использованию углубленных интеграционных методов исследования в изучении роли энергии геннообусловленных патологий, влиянию изменения структуры в цепи ДНК на ее энерговоспринимаемость, энергопроводимость, энергоотдаваемость, энергонакопляемость по цепи ДНК и на развитие геннообусловленных патологий, прежде всего, на развитие злокачественных опухолей. ...

26 10 2026 19:22:59

СОЗДАНИЕ БАЗЫ ДАННЫХ ДЛЯ ГЕОИНФОРМАЦИОННОЙ СИСТЕМЫ

СОЗДАНИЕ БАЗЫ ДАННЫХ ДЛЯ ГЕОИНФОРМАЦИОННОЙ СИСТЕМЫ Статья в формате PDF 253 KB...

24 10 2026 13:41:38

ЧАЙКОВСКИЙ ВИТОЛЬД КАЗИМИРОВИЧ

ЧАЙКОВСКИЙ ВИТОЛЬД КАЗИМИРОВИЧ Статья в формате PDF 327 KB...

18 10 2026 22:54:44

Производство цукатов из мякоти плодов и фруктов

Производство цукатов из мякоти плодов и фруктов Статья в формате PDF 103 KB...

16 10 2026 4:27:48

СТЕРЕОТИПЫ РЕЧЕВОГО ПОВЕДЕНИЯ

СТЕРЕОТИПЫ РЕЧЕВОГО ПОВЕДЕНИЯ Статья в формате PDF 327 KB...

14 10 2026 20:24:46

EPIDEMIOLOGY OF ALLERGIC R HINITS IN PUPILS OF THE REPUBLIC SAKHA(YAKUTIA)

EPIDEMIOLOGY OF ALLERGIC R HINITS IN PUPILS OF THE REPUBLIC SAKHA(YAKUTIA) Статья в формате PDF 99 KB...

11 10 2026 10:26:12

ЗАДАЧИ С ПАРАМЕТРАМИ – УЖАС? НЕТ, ПРОСТО ЗАДАЧИ

ЗАДАЧИ С ПАРАМЕТРАМИ – УЖАС? НЕТ, ПРОСТО ЗАДАЧИ Статья в формате PDF 250 KB...

05 10 2026 16:30:36

ВРОЖДЕННЫЕ ПОРОКИ СЕРДЦА У ДЕТЕЙ В ПЕРМСКОМ КРАЕ

ВРОЖДЕННЫЕ ПОРОКИ СЕРДЦА У ДЕТЕЙ В ПЕРМСКОМ КРАЕ Статья в формате PDF 93 KB...

03 10 2026 3:54:39

ВЛИЯНИЕ ВИРУСНЫХ АНТИГЕНОВ НА МОРФОЛОГИЮ ТИМУСА

ВЛИЯНИЕ ВИРУСНЫХ АНТИГЕНОВ НА МОРФОЛОГИЮ ТИМУСА Статья в формате PDF 258 KB...

30 09 2026 7:44:19

МОНИТОРИНГ ПОЧВ ТЕРРИТОРИИ ИРКУТСКОЙ ОБЛАСТИ

МОНИТОРИНГ ПОЧВ ТЕРРИТОРИИ ИРКУТСКОЙ ОБЛАСТИ Статья в формате PDF 131 KB...

28 09 2026 9:38:59

ОГНЕСТОЙКИЕ НЕНАСЫЩЕННЫЕ ПОЛИЭФИРЫ

ОГНЕСТОЙКИЕ НЕНАСЫЩЕННЫЕ ПОЛИЭФИРЫ Статья в формате PDF 109 KB...

27 09 2026 4:37:29

ПЕРИФЕРИЧЕСКИЙ КРОВОТОК У ДЕВУШЕК 18-22 ЛЕТ В УСЛОВИЯХ ИЗОМЕТРИЧЕСКИХ НАГРУЗОК НАРАСТАЮЩЕЙ ВЕЛИЧИНЫ

ПЕРИФЕРИЧЕСКИЙ КРОВОТОК У ДЕВУШЕК 18-22 ЛЕТ В УСЛОВИЯХ ИЗОМЕТРИЧЕСКИХ НАГРУЗОК НАРАСТАЮЩЕЙ ВЕЛИЧИНЫ В исследовании изучались и оценивались адаптивные реакции периферического звена кровообращения у дeвyшек 18-22 лет при локальных изометрических нагрузках нарастающей величины. Хаpaктерно, что с ростом прессорных воздействий на сосуды работающих мышц объемная скорость кровотока оставалась более высокой, чем в покое. В целом качество приспособительных реакций кровообращения было ниже у дeвyшек 18 лет при низких объемах выполняемой работы. ...

25 09 2026 16:48:58

Еще:
Поддержать себя -1 :: Поддержать себя -2 :: Поддержать себя -3 :: Поддержать себя -4 :: Поддержать себя -5 :: Поддержать себя -6 :: Поддержать себя -7 :: Поддержать себя -8 :: Поддержать себя -9 :: Поддержать себя -10 :: Поддержать себя -11 :: Поддержать себя -12 :: Поддержать себя -13 :: Поддержать себя -14 :: Поддержать себя -15 :: Поддержать себя -16 :: Поддержать себя -17 :: Поддержать себя -18 :: Поддержать себя -19 :: Поддержать себя -20 :: Поддержать себя -21 :: Поддержать себя -22 :: Поддержать себя -23 :: Поддержать себя -24 :: Поддержать себя -25 :: Поддержать себя -26 :: Поддержать себя -27 :: Поддержать себя -28 :: Поддержать себя -29 :: Поддержать себя -30 :: Поддержать себя -31 :: Поддержать себя -32 :: Поддержать себя -33 :: Поддержать себя -34 :: Поддержать себя -35 :: Поддержать себя -36 :: Поддержать себя -37 :: Поддержать себя -38 ::