Эйлер vs РК4 (Маятник)
Методы численного интегрирования — это важнейшие инструменты для решения дифференциальных уравнений, не имеющих аналитического решения, например, для описания движения нелинейного маятника. Данный симулятор наглядно сравнивает два фундаментальных алгоритма: простой (или явный) метод Эйлера и метод Рунге-Кутты четвёртого порядка (РК4). Оба метода решают одно и то же обыкновенное дифференциальное уравнение (ОДУ), выведенное из второго закона Ньютона для маятника: d²θ/dt² = -(g/L) sin(θ), где θ — угловое смещение, g — ускорение свободного падения, а L — длина маятника. Уравнение упрощено путём пренебрежения сопротивлением воздуха и предположения о невесомости стержня, чтобы сосредоточиться исключительно на нелинейной возвращающей силе. Симулятор преобразует это ОДУ второго порядка в систему двух ОДУ первого порядка: dθ/dt = ω и dω/dt = -(g/L) sin(θ). Начиная с одинаковых начальных условий и используя один и тот же фиксированный шаг по времени (h), оба алгоритма продвигают состояние (θ, ω) вперёд. Метод Эйлера использует простую линейную проекцию: θ_{n+1} = θ_n + h * ω_n и ω_{n+1} = ω_n + h * [-(g/L) sin(θ_n)]. В отличие от него, РК4 вычисляет четыре взвешенных промежуточных наклона, что даёт гораздо более точную оценку среднего наклона на интервале. Взаимодействуя с симуляцией, студенты могут наблюдать, как простота метода Эйлера приводит со временем к значительному росту энергии (увеличению амплитуды) — нефизическому артефакту, — в то время как РК4 сохраняет энергию и период почти постоянными гораздо дольше, близко аппроксимируя истинную консервативную систему. Это демонстрирует ключевые концепции вычислительной физики: локальная ошибка усечения, накопление глобальной ошибки, устойчивость и компромисс между вычислительной стоимостью и точностью.
Для кого: Студенты бакалавриата по физике, инженерии или вычислительным наукам, изучающие численные методы решения дифференциальных уравнений или нелинейную динамику.
Ключевые понятия
- Численное интегрирование
- Методы Рунге-Кутты
- Метод Эйлера
- Обыкновенное дифференциальное уравнение (ОДУ)
- Нелинейный маятник
- Ошибка усечения
- Сохранение энергии
- Фазовое пространство
Как это работает
Один и тот же **нелинейный** маятник **θ'' = −(g/L) sin θ**, один и тот же **начальный** угол, один и тот же **фиксированный шаг** **h**: **явный метод Эйлера** (обновление **θ**, затем **ω**) **приобретает** **энергию**, и **фаза** **уходит** сильно при больших **h**; **RK4** остаётся **гораздо** ближе к **высокоточному** эталону при той же **работе** на шаг. Измените **h**, чтобы увидеть, как **Эйлер** **расходится** первым.
Основные формулы
Часто задаваемые вопросы
- Почему маятник, смоделированный методом Эйлера, набирает энергию и раскачивается выше, что, кажется, нарушает законы физики?
- Рост энергии — это численный артефакт, а не физическое явление. Явный метод Эйлера обладает свойством численной неустойчивости для данного типа ОДУ. Его упрощённая линейная экстраполяция последовательно переоценивает изменение угловой скорости, добавляя небольшое количество энергии на каждом шаге. В реальной консервативной системе полная механическая энергия постоянна. Это накопление ошибки подчёркивает ключевое ограничение методов низкого порядка при долгосрочном моделировании.
- Всегда ли РК4 лучше метода Эйлера? Зачем тогда вообще использовать метод Эйлера?
- Хотя РК4 гораздо точнее при заданном шаге, он требует четырёх вычислений функции на шаг против одного у метода Эйлера, что делает его вычислительно более затратным. Метод Эйлера всё ещё используется в ситуациях, где первостепенны простота и скорость, а ошибка допустима, например, в физике некоторых видеоигр реального времени или для очень простых систем с чрезвычайно малым шагом. Он также служит базовой концепцией для понимания более сложных методов.
- Что произойдёт, если я сделаю шаг по времени (h) очень маленьким для обоих методов?
- При уменьшении шага по времени решения обоих методов будут сходиться к истинному решению, а дрейф энергии в методе Эйлера на фиксированном интервале времени станет менее выраженным. Однако, чтобы точность метода Эйлера соответствовала точности РК4 при умеренном шаге, потребуется значительно меньший (и вычислительно неподъёмный) шаг. Это иллюстрирует концепцию 'порядка': ошибка РК4 масштабируется как h^4, а ошибка Эйлера — только как h.
- Моделирует ли эта симуляция реальный маятник идеально?
- Нет, она моделирует идеализированный нелинейный маятник. Ключевые упрощения включают отсутствие трения или сопротивления воздуха (демпфирования), абсолютно жёсткий и невесомый стержень, а также точечную массу груза. Эти предположения позволяют изолировать и изучить основные ошибки численного интегрирования без влияния побочных физических эффектов. В реальном маятнике энергия постепенно уменьшалась бы из-за демпфирования, а не возрастала, как в методе Эйлера.
Ещё из «Визуализация математики»
Другие симуляторы в этой категории — или все 26.
Лотка–Вольтерра
N′ = αN−βNP, P′ = δNP−γP; фазовая плоскость RK4; точка равновесия.
Логистический рост
dN/dt = rN(1−N/K); точная S-образная кривая и ёмкость среды K.
Матрица 2×2 и собственные векторы
Деформация сетки под действием M; стрелки вдоль собственных направлений для вещественных λ.
Сглаживание Савицкого–Голэя
Зашумлённый косинус vs свёртка SG(7,2) — лучше сохраняет пики по сравнению с широким прямоугольным фильтром.
Цепь Маркова (Погода)
Двухсостоянияя цепь Солнце/Дождь: матрица P, стационарное распределение π, эмпирическое и теоретическое.
Градиентный спуск (2D)
Линии уровня функции f(x,y) и траектория (x,y) ← (x,y) − η∇f; чаша или эллиптическая впадина.