Дифференциальные уравнения
Волновое уравнение
Звук гитары, землетрясение в Японии, сигнал WiFi и свет от звёзд - все описываются одним уравнением. Разница лишь в скорости c и размерности пространства. Волновое уравнение - самое универсальное в физике.
- **Синтез звука:** тембр инструмента = набор гармоник с коэффициентами Aₙ, Bₙ. Физическое моделирование гитары в Logic Pro/Ableton - численное решение волнового уравнения для струны.
- **Сейсмическое моделирование:** FDTD-симуляция распространения сейсмических волн для поиска нефти и предсказания землетрясений.
- **FDTD в электронике:** моделирование сигналов в PCB-платах, антеннах и оптических волноводах - основной инструмент проектирования RF-компонентов.
Предварительные знания
Вывод волнового уравнения
Рассмотрим колебания струны с линейной плотностью ρ под натяжением T. Сила, действующая на элемент [x, x+dx] в поперечном направлении: T·∂²u/∂x²·dx. Второй закон Ньютона: ρ·dx·∂²u/∂t² = T·∂²u/∂x²·dx. Сокращая dx: **∂²u/∂t² = c²∂²u/∂x²**, где **c = √(T/ρ)** - скорость волны.
**Волновое уравнение** - **гиперболическое** PDE. В отличие от теплопроводности (диссипация), волновое уравнение описывает сохраняющиеся колебания. Скорость c - единственный параметр: √(T/ρ) для струны, √(E/ρ) для упругих волн, 1/√(με) для ЭМ волн.
| Система | c (скорость волны) | Типичное значение |
|---|---|---|
| Звук в воздухе | √(γP/ρ) | 343 м/с |
| Звук в воде | √(K/ρ) | 1480 м/с |
| Свет в вакууме | 1/√(ε₀μ₀) | 3×10⁸ м/с |
| Сейсмические P-волны | √((K+4G/3)/ρ) | ~6000 м/с |
| Струна гитары | √(T/ρ) | 200-500 м/с |
Натяжение струны увеличили в 4 раза, плотность уменьшили в 2 раза. Как изменится скорость волны?
Решение Д'Аламбера
Общее решение волнового уравнения ∂²u/∂t² = c²∂²u/∂x² на всей прямой: **u(x,t) = f(x-ct) + g(x+ct)**. Функция f(x-ct) - волна, бегущая вправо со скоростью c. Функция g(x+ct) - волна, бегущая влево. Это **решение Д'Аламбера** (1747).
При начальных условиях u(x,0) = φ(x) и ∂u/∂t(x,0) = ψ(x) решение: **u(x,t) = [φ(x+ct) + φ(x-ct)]/2 + (1/2c)∫_{x-ct}^{x+ct} ψ(s)ds**.
**Принцип Гюйгенса** (в 3D): волна от точечного источника распространяется по сфере. В 3D информация приходит ровно через время r/c - «острый фронт». В 2D (или 1D) есть «послесвечение» - решение ненулевое за фронтом волны. Это отличие 3D от 2D.
Д'Аламбер, Эйлер и спор о струне (1747-1761)
В 1747 Д'Аламбер опубликовал решение u = f(x-ct) + g(x+ct) для колебаний натянутой струны. Эйлер тут же ответил, что f и g могут быть «любыми» функциями, в том числе с углами. Д'Аламбер настаивал на гладкости. Спор тянулся до 1760-х. Бернулли подлил масла в огонь: предложил разлагать решение в ряд синусов - тот самый ряд, который потом стал фурье-рядом. Эта склока заложила язык, которым физики до сих пор описывают акустику, оптику и квантовую механику.
Современное понятие 'обобщённого решения' и теория распределений Соболева-Шварца выросли именно из этого спора.
Начальный импульс u(x,0) = δ(x), ∂u/∂t(x,0) = 0. Какова форма u(x,t) при t > 0?
Стоячие волны
На конечном промежутке [0, L] с условиями u(0,t) = u(L,t) = 0 разделение переменных даёт: **uₙ(x,t) = sin(nπx/L)·[Aₙcos(ωₙt) + Bₙsin(ωₙt)]**, где **ωₙ = cnπ/L** - собственные частоты. Это стоячие волны: узлы неподвижны, пучности колеблются.
**Музыкальные инструменты:** частоты гитарной струны fₙ = cn/(2L). Первая гармоника (основной тон), вторая, третья - обертоны. Тембр инструмента определяется набором Aₙ, Bₙ - коэффициентами Фурье начального воздействия.
Струна длиной L=1 м с c=300 м/с. Какова частота второй гармоники?
Сохранение энергии
Полная энергия волны: **E(t) = (ρ/2)∫[(∂u/∂t)² + c²(∂u/∂x)²]dx**. Первое слагаемое - кинетическая энергия, второе - потенциальная (упругая). Для волнового уравнения без диссипации **E(t) = const** - энергия сохраняется.
**FDTD (Finite-Difference Time-Domain)** - стандартный метод моделирования ЭМ волн, акустики и сейсмики. Используется в проектировании антенн, ультразвуковой диагностике, предсказании землетрясений. scipy.integrate.odeint + пространственная дискретизация = элементарный FDTD.
В численной симуляции волнового уравнения энергия со временем растёт. Что это означает?
Ключевые идеи
- **∂²u/∂t² = c²∂²u/∂x²** - гиперболическое PDE, описывающее консервативное распространение.
- **Решение Д'Аламбера:** u = f(x-ct) + g(x+ct). Два волновых пакета, бегущих в противоположных направлениях.
- **Стоячие волны:** собственные моды на [0,L] с частотами ωₙ = cnπ/L. Основа физики музыкальных инструментов.
- **Сохранение энергии:** E = const без диссипации. Нарушение = признак численной неустойчивости.
Связанные темы
Волновое уравнение - прототип гиперболических PDE:
- Уравнение теплопроводности — Замена ∂²/∂t² → ∂/∂t переходит от волн к диффузии
- Уравнение Лапласа — Стационарное волновое уравнение (ω→0) - уравнение Лапласа
- Классификация PDE — Волновое уравнение - эталонный пример гиперболического типа
Вопросы для размышления
- Принцип Гюйгенса в 3D: волна от точки доходит с острым фронтом. В 2D - есть послесвечение. Почему телефонная связь работает (3D), а 2D-«телефон» давал бы эхо?
- Тембр инструмента определяется коэффициентами Фурье начального условия. Как изменится тембр гитары, если ущипнуть струну в точке L/2 вместо L/4?
- Волновое уравнение описывает и квантовую механику (уравнение Шрёдингера - тоже волновое). В чём принципиальное отличие квантовой волны от классической?