Физика формулы Эйлера: как комплексные экспоненты решают дифференциальные уравнения

Физика формулы Эйлера: как комплексные экспоненты решают дифференциальные уравнения

Спикер: 3Blue1Brown | Длительность: 27:49

Транскрипт

Введение

Это первое видео в трилогии, цель которой — демистифицировать преобразование Лапласа. Хотя само преобразование рассматривается в следующих главах, здесь закладываются ментальные модели и базовые знания, которые делают его понимание максимально доступным.

Главные герои всей трилогии — экспоненциальные функции вида e^(st), где t — время, а s — некоторое число, определяющее конкретную экспоненту. Главная цель видео — показать на примерах из физики, почему полезно разрешить s принимать не только вещественные, но и комплексные значения.

Формула Эйлера через динамику

Всё начинается с ключевого свойства: e^t — собственная производная. Но вместо того чтобы визуализировать производную как наклон графика, предлагается думать о ней динамически. Если e^t задаёт положение точки на числовой прямой, то производная говорит: вектор скорости всегда совпадает с вектором положения.

Для e^(2t) скорость вдвое больше положения — рост ускоряется ещё быстрее. Для e^(−0.5t) скорость направлена в противоположную сторону и вдвое меньше — экспоненциальное затухание.

Ключевой момент: что если константа мнимая, s = i? Умножение на i — это поворот на 90°. Значит, скорость всегда перпендикулярна положению. Единственное движение с таким свойством — вращение по окружности. Начальное положение — 1, значит это единичная окружность, и точка проходит один радиан за единицу времени.

Через π единиц времени точка оказывается в −1 — вот вам формула Эйлера e^(πi) = −1.

Важное замечание: нотация обманчива. При комплексном входе выражение e^x не имеет ничего общего с многократным умножением. Вычисление на самом деле — это подстановка в бесконечный полином (ряд Тейлора). Но фокусировка на свойстве «собственная производная» гораздо полезнее, чем на вычислительной процедуре.

Когда s имеет и вещественную, и мнимую часть (например, s = −0.5 + i), экспонента описывает спираль: вращение с затуханием. Инженеры называют комплексную плоскость значений s «S-плоскостью» — каждая точка кодирует целое поведение: мнимая часть — частоту колебаний, вещественная — рост или затухание.

Гармонический осциллятор

Масса на пружине — центральный пример физики. Сила пропорциональна смещению: F = −kx. По второму закону Ньютона: mx'' = −kx. Часто добавляют демпфирование — член, пропорциональный скорости: mx'' = −kx − μx'.

Физическая интуиция подсказывает: решение — затухающие колебания. Но как решить уравнение формально?

Применяется «странный трюк» — подставляем x = e^(st). Каждая производная добавляет множитель s: x' = s·e^(st), x'' = s²·e^(st). Экспонента выносится за скобку, и остаётся квадратное уравнение относительно s.

Без демпфирования (μ = 0): s² = −k/m, откуда s = ±i√(k/m) = ±iω. Комплексные числа возникают неизбежно! Чем сильнее пружина (больше k), тем быстрее колебания (больше ω) — физическая интуиция подтверждается.

Два комплексных решения e^(iωt) и e^(−iωt) при сложении дают вещественный результат: 2cos(ωt). Линейность уравнения позволяет масштабировать и складывать решения для подбора под начальные условия.

С демпфированием: квадратная формула даёт корни с отрицательной вещественной частью — затухающие колебания. При достаточно сильном демпфировании мнимая часть исчезает — передемпфированный режим, чистое затухание без колебаний.

Обобщение на линейные уравнения произвольного порядка

Тот же трюк работает для любого линейного уравнения n-го порядка с постоянными коэффициентами. Подстановка e^(st) даёт полином степени n. По фундаментальной теореме алгебры он имеет n корней (возможно, комплексных). Каждый корень даёт решение, и общее решение — линейная комбинация всех n экспонент.

Путь к преобразованию Лапласа

Для нелинейных и неоднородных уравнений (например, вынужденный осциллятор с внешней силой) простой трюк подстановки не работает напрямую. Но решения всё равно выглядят как комбинации экспонент — только с конкретными, а не свободными коэффициентами.

Преобразование Лапласа — систематический инструмент для нахождения этих экспоненциальных компонент. Оно переводит функции на «язык», где e^(st) — базовые единицы, и дифференцирование выглядит как умножение на s. Дифференциальные уравнения превращаются в алгебру.

Практические задания

Задание 1: Динамическое чтение экспоненты

Нарисуйте числовую прямую и отметьте точку в позиции 1. Для каждого из значений s = 2, s = −0.5, s = i, s = −0.5 + i нарисуйте вектор скорости как модифицированную копию вектора положения. Для вещественных s — на прямой, для комплексных — в плоскости. Проследите траекторию точки за несколько шагов. Убедитесь, что вы видите экспоненциальный рост, затухание, вращение и спираль.

Задание 2: Проверка формулы Эйлера через ряд Тейлора

Выпишите первые 8 членов ряда Тейлора для e^x и подставьте x = πi. Вычислите каждый член: помните, что i² = −1, i³ = −i, i⁴ = 1. Сложите вещественные и мнимые части отдельно. Убедитесь, что вещественная часть стремится к −1, а мнимая — к 0. Нарисуйте спиральную сумму частичных сумм в комплексной плоскости.

Задание 3: Решение уравнения пружины подстановкой

Возьмите уравнение x'' + 4x = 0 (здесь k/m = 4). Подставьте x = e^(st) и получите характеристическое уравнение s² + 4 = 0. Найдите корни s = ±2i. Запишите общее решение как c₁·e^(2it) + c₂·e^(−2it). Покажите, что при c₁ = c₂ = 1 получается 2cos(2t). Подберите коэффициенты для начальных условий x(0) = 3, x'(0) = 0.

Задание 4: Исследование затухающего осциллятора

Возьмите уравнение x'' + 2x' + 5x = 0. Подставьте e^(st) и получите s² + 2s + 5 = 0. Примените формулу дискриминанта: D = 4 − 20 = −16. Найдите корни s = −1 ± 2i. Отметьте их на S-плоскости. Объясните, почему вещественная часть −1 означает затухание, а мнимая часть ±2 — частоту колебаний.

Задание 5: Визуализация S-плоскости

Нарисуйте комплексную плоскость и разделите её на области. Отметьте: правая полуплоскость — экспоненциальный рост, левая — затухание, мнимая ось — чистые колебания, начало координат — постоянная функция. Для каждой из точек s = 3, s = −2, s = 4i, s = −1 + 3i, s = −1 − 3i нарисуйте примерный вид графика Re(e^(st)) от времени.

Задание 6: Переход к передемпфированному режиму

Возьмите уравнение x'' + μx' + x = 0 (k/m = 1). Вычислите дискриминант D = μ² − 4. Найдите критическое значение μ = 2, при котором D = 0. Для μ = 1 (недодемпфирование), μ = 2 (критическое) и μ = 4 (передемпфирование) найдите корни и опишите качественное поведение системы.

Лучшие цитаты

«Вы должны думать об этом как о том, что определяет число e.» — 3Blue1Brown

«Возьмите привычку мыслить более гибко.» — 3Blue1Brown

«Единственное движение, удовлетворяющее этому критерию — вращение по окружности.» — 3Blue1Brown

«Когда вы подставляете комплексное значение, выражение на самом деле имеет очень мало общего с многократным умножением и, честно говоря, не так уж много — с числом e.» — 3Blue1Brown

«Хотите вы того или нет, мнимая единица i вошла в игру.» — 3Blue1Brown

«Угадывание и проверка выглядят так, будто от студента требуется знать ответ заранее.» — 3Blue1Brown

«Можно думать об этих функциях e^(st) как об атомах анализа.» — 3Blue1Brown

«Дифференцирование по времени выглядит как умножение на s — и дифференциальные уравнения начинают выглядеть как алгебра.» — 3Blue1Brown