- •Оглавление
- •Введение
- •1. Обыкновенные дифференциальные уравнения и системы
- •1.1. Обыкновенные дифференциальные уравнения. Основные понятия
- •1.2. Системы оду. Основные понятия
- •1.3. Связь оду высших порядков и систем оду
- •2.1.2. Уравнения с разделяющимися переменными
- •2.1.3. Однородные уравнения 1-го порядка
- •2.1.4. Уравнения, приводящиеся к однородным
- •2.1.5. Линейные уравнения первого порядка
- •2.1.6. Уравнения Бернулли
- •2.1.7. Уравнения в полных дифференциалах
- •2.2. Обыкновенные дифференциальные уравнения првого порядка. Поведение решений
- •2.2.1. Теорема существования и единственности решения задачи Коши
- •2.2.2. Уравнения первого порядка. Поле направлений
- •2.2.3. Автономные уравнения первого порядка
- •2.2.4. Устойчивость решений обыкновенных дифференциальных уравнений первого порядка
- •2.2.5. Асимптотическая устойчивость решений обыкновенных дифференциальных уравнений первого порядка
- •2.3. Метод изоклин
- •3. Обыкновенные дифференциальные уравнения высших порядков
- •3.1. Обыкновенные дифференциальные уравнения высших порядков. Понижение порядка
- •3.1.1. Понижение порядка обыкновенного дифференциального уравнения. Введение
- •3.1.2. Уравнения, не содержащие независимой переменной
- •3.1.3. Уравнения, не содержащие искомой функции
- •3.1.4. Уравнения с однородной правой частью
- •3.2. Линейные однородные дифференциальные уравнения n-го порядка
- •3.2.1. Линейные однородные дифференциальные уравнения n-го порядка. Введение
- •3.2.2. Свойства решений линейного уравнения. Принцип суперпозиции
- •3.2.3. Существование и единственность решения задачи Коши
- •3.2.4. Линейные уравнения второго порядка. Гармонические колебания
- •3.2.5. Линейные уравнения втрого порядка. Ангармонические колебания
- •3.2.6. Линейные уравнения второго порядка. Уравнение Ньютона
- •3.3. Линейная зависимость и линейная независимость системы функций
- •3.3.1. Линейная зависимость и линейная независимость системы функций
- •3.3.2. Определитель Вронского
- •3.3.3. Исследование линейной независимости системы функций
- •3.3.4. Линейная независимость решений линейного дифференциального уравнения
- •3.4. Структура решения линейного однородного дифференциального уравнения n-го порядка
- •3.4.1. Фундаментальная система решений однородного линейного дифференциального уравнения
- •3.4.2. Структура общего решения линейного однородного уравнения
- •3.4.3. Структура общего решения линейного неоднородного уравнения
- •3.4.4. Метод вариации произвольных постоянных отыскания частного решения
- •3.5. Линейные однородные дифференциальные уравнения с постоянными коэффициентами
- •3.5.1. Решение однородного уравнения с постоянными коэффициентами
- •3.5.2. Метод подбора построения частного решения неоднородного уравнения
- •3.5.3. Уравнение Эйлера
- •4. Системы дифференциальных уравнений
- •4.1. Системы обыкновенных дифференциальных уравнений. Основные понятия
- •4.2. Фазовое пространство. Фазовые траектории
- •4.3. Существование и единственность решения задачи Коши
- •4.4. Интегрирование систем дифференциальных уравнений методом исключения
- •4.4.1. Линейные системы обыкновенных дифференциальных уравнений. Структура
- •4.4.1.1. Линейные системы обыкновенных дифференциальных уравнений. Основные понятия
- •4.4.1.2. Фундаментальная матрица решений однородной линейной системы дифференциальных уравнений
- •4.4.1.3. Структура общего решения однородной линейной системы дифференциальных уравнений
- •4.4.1.4. Структура общего решения неоднородной линейной системы дифференциальных уравнений
- •4.4.1.5. Построение фундаментальной матрицы решений однородной линейной системы дифференциальных уравнений с постоянными коэффициентами методом Эйлера
- •4.4.2. Системы обыкновенных дифференциальных уравнений. Поведение решений
- •4.4.2.1. Устойчивость решений систем дифференциальных уравнений
- •4.4.2.2. Устойчивость и асимптотическая устойчивость по Ляпунову
- •4.4.2.3. Устойчивость положения равновесия линейных систем оду
- •4.4.2.4. Устойчивость точек покоя нелинейных систем по линейному приближению
- •4.4.2.5. Неустойчивость по линейному приближению точек покоя нелинейных систем
- •4.4.3. Автономные системы дифференциальных уравнений
- •4.4.3.1. Автономные системы. Основные понятия
- •4.4.3.2. Свойства фазовых траекторий
- •4.4.3.3. Фазовая плоскость, фазовые кривые, фазовый портрет автономной системы второго порядка
- •4.4.3.4. Векторное поле автономной системы второго порядка
- •4.4.3.5. Точки покоя линейной автономной системы второго порядка с постоянными коэффициентами
- •4.4.4. Численное интегрирование дифференциальных уравнений первого порядка
- •4.4.4.1. Задача Коши. Общие замечания. Постановка задачи
- •4.4.4.2. Метод Эйлера
- •4.4.4.3. Модифицированный метод Эйлера
- •4.4.4.4. Метод Рунге-Кутта
- •2. Решение уравнения модифицированным методом Эйлера:
- •3. Решение уравнения методом Рунге-Кутта:
- •4. Аналитическое решение заданного уравнения:
- •5. Сравнение точного решения и приближенных решений исходного дифференциального уравнения:
- •Библиографический список
- •Обыкновенные дифференциальные уравнения
2.2.2. Уравнения первого порядка. Поле направлений
Рассмотрим обыкновенное дифференциальное уравнение 1-го порядка с правой частью f(x, y), определённой в области D ⊂ R2:
Здесь x независимая переменная (аргумент), y = y(x) неизвестная функция.
Если y = y(x) решение уравнения, то соответствующая интегральная кривая (график решения y = y(x)) в каждой своей точке (x, y(x)) имеет касательную с угловым коэффициентом
Через каждую точку (x, y) области D ⊂ R2 можно провести небольшой отрезок с угловым коэффициентом k .
Выполнив такое построение для всех узлов некоторой прямоугольной сетки в области D ⊂ R2 , получим изображение поля направлений.
Если узлы сетки расположены «достаточно часто» поле направлений дает полную картину поведения интегральных кривых.
На рисунках приведены изображение поля направлений и изображение поля направлений с несколькими интегральными кривыми.
Рассмотрев внимательно рисунки, можно увидеть, что отрезки, изображающие поле направлений действительно указывают направление касательных к интегральным кривым.
Можно видеть, что “аккуратно“ изображенное поле направлений дает достаточно полное представление о поведении интегральных кривых.
В инженерных задачах для того чтобы сформулировать содержательные утверждения об исследуемом процессе, бывает достаточно внимательно изучить поле направлений.
Дифференциальное уравнение задает поле направлений, которое позволяет судить о наиболее характерных особенностях поведения решений уравнения.
Пример №1
Рассмотрим обыкновенное дифференциальное уравнение 1-го порядка
с правой частью, определённой в области y ≠ 0.
На рисунках приведены изображения поля направлений уравнения и поля направлений с несколькими интегральными кривыми.
На рисунке видно, что в точках пересечения с осью абсцисс (y = 0) касательные к интегральным кривым перпендикулярны оси 0y . Это означает, в частности, что в этих точках скорость изменения решения возрастает до бесконечности.
Видно также, что по мере удаления точки (x, y(x)) (x → ∞,y(x) → ∞) скорость изменения решения стабилизируется − есть основания предполагать, что интегральные кривые имеют наклонную асимптоту.
Пример №2
Рассмотрим обыкновенное дифференциальное уравнение 1-го порядка с правой частью, определённой в области x ≠ 0, x2 + (y + 0.5)2 ≥ 0.5:
На рисунках приведено изображение поля направлений уравнения в области, где определена правая часть уравнения.
2.2.3. Автономные уравнения первого порядка
Автономным уравнением первого порядка называется уравнение вида
правая часть f(y) которого не зависит от x.
Решение y(x) ≡ 0 называется неподвижной точкой уравнения.
Для автономного уравнения решение задачи Коши
определяется равенством:
Например, дифференциальное уравнение
описывает движение материальной точки по оси x под действием внешних сил f(x). Координата точки x(t) в момент времени t – решение дифференциального уравнения. Скорость движения точки задана функцией f(x). В момент времени t координата x(t) точки, которая в начальный момент t0 имела координату x0, определяется равенством
При t → ∞ функция x(t) может быть ограниченной, может стремиться к конечному пределу, может ”уходить на бесконечность” (x(t) → ∞) или может быть неограниченной с каким-то более сложным поведением.
Пример №1
Исследуем поведение решений дифференциального уравнения y ' = cos πy .
В какой момент времени x1 решение с начальным условием y(0) = 0 достигнет значения 0.5? Является ли точка y = 0.5 неподвижной точкой? Решения каких типов имеет уравнение?
Время x1, за которое решение задачи Коши y ' = cosπy , y(0) = 0 достигнет значения 0.5, определяется равенством
Этот интеграл расходится. Значит точка, начавшая из нуля движение вдоль оси x, никогда не придет в точку y = 0.5 .
На рисунке видно, что интегральная кривая, проходящая через начало координат асимптотически приближается снизу к прямой y = 0.5 .
Точка y = 0.5 − неподвижная точка уравнения. Действительно, функция y = 0.5 является решением задачи Коши y ' = cos πy , y(0) = 0.5 . Это означает, что точка, начавшая движение из y = 0.5 остается на месте.
Для того чтобы ответить на вопрос о типах решений уравнения, изобразим его поле направлений.
Видно, что уравнение имеет ограниченные решения в областях
Уравнение имеет неподвижные точки: