- •Оглавление
- •Введение
- •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. Сравнение точного решения и приближенных решений исходного дифференциального уравнения:
- •Библиографический список
- •Обыкновенные дифференциальные уравнения
3.2.5. Линейные уравнения втрого порядка. Ангармонические колебания
Рассмотрим линейное однородное дифференциальное уравнение 2-го порядка с постоянными коэффициентами a и b
y '' + 2ay ' + by = 0.
Если a2 < b и a > 0 , то общее решение уравнения имеет вид
Его можно записать в виде
Функция y(x) = A exp(− ax) cos(ωx − φ) непериодическая, но ее нули, максимумы и минимумы, повторяются с периодом T = 2π/ω . Этот период равен периоду колебаний гармонического осциллятора с частотой ω.
Колебания y(x) = Aexp (− ax)cos(ωx − φ) называются затухающими колебаниями, y(x) → 0 при x → + ∞ .
Величина Aexp (− ax) называется амплитудой колебаний; величина δ = a коэффициентом затухания.
Коэффициент затухания δ записывают в виде δ =1/τ , τ – время, за которое амплитуда колебаний уменьшится в e раз.
Величина d = δ · T = 2πδ/ω называется логарифмическим декрементом затухания. Она показывает, насколько уменьшится амплитуда за один период. Логарифмический декремент – естественная мера скорости затухания.
Колебания, которые описывает дифференциальное уравнение y '' + 2ay ' + by = 0 , называют ангармоническими колебаниями.
На рисунке приведены графики решений уравнения y '' + y ' + 4y =0 , проходящих через точку (0,0) с различными начальными скоростями y '(0)
Пример №1
Найдем коэффициент затухания и логарифмический декремент для уравнения y '' + y ' + 4y = 0 .
В уравнении y '' + y ' + 4y = 0 параметры a , b и ω имеют значения
Вычислим коэффициент затухания δ и логарифмический декремент d:
Итак, коэффициент затухания δ = 0.5, логарифмический декремент d ≈ 1.62 .
3.2.6. Линейные уравнения второго порядка. Уравнение Ньютона
Уравнением Ньютона называют линейное неоднородное дифференциальное уравнение 2-го порядка с постоянными коэффициентами вида
y '' + ω02y = Fcosωx.
Известно, что однородное уравнение y '' + ω0y = 0 описывает свободные колебания материальной точки с частотой колебаний ω0
y(x) = C1cos(ω0x) + C2sin(ω0x) .
На рисунке приведен график решения однородного уравнения свободных колебаний с частотой ω0 = 0.1
Неоднородное уравнения с правой частью Fcos(ωx) описывает колебания материальной точки под действием внешней периодической силы, частота которой ω.
Если ω0 ≠ ω , то общее решение уравнения имеет вид
Если ω0 = ω , то
На рисунках приведены решения уравнений Ньютона для различных частот свободных колебаний и частот внешней силы. Рисунки иллюстрируют принцип суперпозиции и модулирование амплитуды колебаний.
На приведенном ниже рисунке показан резонансный случай – частота свободных колебаний ω0 =2 совпадает с частотой внешней силы ω =2 .
Пример №1
Исследуем поведение решения задачи Коши для уравнения Ньютона с частотой свободных колебаний ω0 = 5 и частотой внешней силы ω = 5.1:
y '' + 52y = cos5.1x, y(0)=0, y'(0)=1.
Поскольку ω0 ≠ ω (5 ≠ 5.1) , то общее решение уравнения Ньютона имеет вид
Вычислив решение задачи Коши y '' + 52y = cos5.1x, y(0)=0, y'(0)=1 по этой формуле, получим
Для исследования решение удобнее записать в виде y(x) = A(x) cos(5x− φ(x)) − в виде произведения амплитуды A(x) и периодически изменяющейся компоненты cos(5x− φ(x)).
Выполнив несложные, но громоздкие тригонометрические преобразования имеем:
На рисунках приведены: график медленно меняющейся амплитуды A(x), график быстро меняющихся колебаний cos(5x− φ(x)) и график их произведения − график решения задачи Коши для уравнения Ньютона. Приведенное на рисунке поведение решения называют «биением».