Спотлайт #28 – Акустика та музика: стоячі хвилі, фігури Хладні, кімнатні моди та фізика звуку

Акустика — найдавніша галузь фізики. Піфагор вивчав відношення довжин струн. Галілей скидав маятники, щоб рахувати частоти. Хладні розсипав пісок на скрипкових пластинах і розкрив форми звукових мод. Сабін вимірював реверберацію в лекційних залах Гарварду й заснував архітектурну акустику. Сьогодні ця галузь охоплює проєктування концертних залів, підводний сонар, медичний ультразвук, навушники з шумозаглушенням та кохлеарні імпланти. Все це об'єднує хвильове рівняння — просте диференціальне рівняння в частинних похідних, що керує тим, як поширюються й інтерферують збурення тиску.

Симуляції в цій хвилі досліджують акустику на шести масштабах: одновимірні стоячі хвилі в трубі, двовимірну вузлову геометрію на вібруючій пластині, тривимірні резонанси приміщення, час реверберації та відображення частоти на положення всередині равлика. Кожен масштаб розкриває щось інше про поведінку хвиль при зіткненні з межами.

1. Стоячі хвилі в трубах — гармоніки й граничні умови

Коли синусоїдальна хвиля рухається трубою і відбивається від її кінця, пряма й відбита хвилі накладаються, утворюючи стоячу хвилю. Просторовий візерунок застигає на місці: вузли (нульове зміщення) та пучності (максимальне зміщення) виникають у фіксованих позиціях. Дозволені частоти цілком визначаються довжиною труби та її граничними умовами на кожному кінці.

Хвильове рівняння (1D, тиск):
  ∂²p/∂t² = c² ∂²p/∂x²   де c = швидкість звуку (343 м/с у повітрі при 20°C)

Граничні умови (для тиску p та зміщення u ≡ −∂p/∂x):
  Відкритий кінець:   вузол тиску      p = 0 при x = 0 або L
  Закритий кінець: вузол зміщення  ∂p/∂x = 0 при x = 0 або L

Гармонічні частоти:
  Відкрита-відкрита (обидва кінці відкриті):
    f_n = n · c / (2L),  n = 1, 2, 3, …  (присутні всі гармоніки)
    Пучність на кожному кінці; основний тон = c/2L

  Відкрита-закрита (один кінець відкритий, один закритий):
    f_n = (2n−1) · c / (4L),  n = 1, 2, 3, …  (лише непарні гармоніки)
    Основний тон = c/4L, на октаву нижче за відкриту-відкриту труби тієї ж довжини
    Пояснює, чому закрита органна труба звучить на октаву нижче

  Закрита-закрита (обидва кінці закриті):
    f_n = n · c / (2L),  n = 1, 2, 3, …  (ті самі частоти, що й відкрита-відкрита,
    але вузол на кожному кінці замість пучності)

Стояча хвиля зміщення (відкрита-відкрита, n-та мода):
  u(x, t) = A · cos(nπx/L) · cos(2πf_n t)
  Вузли при x = L/(2n), 3L/(2n), …
  Пучності при x = 0, L/n, 2L/n, …

Реальні інструменти (кінцева поправка):
  Відкритий кінець — не ідеальний вузол тиску; ефективна довжина дорівнює L + 0,6r
  де r = радіус труби (кінцева поправка Релея).
  Враховує часткове випромінювання на відкритому кінці.

Різниця між трубами відкрита-відкрита та відкрита-закрита пояснює тембральну різницю між флейтою й кларнетом. Флейта відкрита з обох кінців і виробляє всі гармоніки; кларнет фактично закритий на тростині й відкритий на розтрубі, тому переходить на третю гармоніку (дуодециму над основним тоном), а не на октаву, а його нижній регістр домінується основним тоном та непарними гармоніками. Парні гармоніки, відсутні в спектрі одинарної тростини, надають кларнету його характерного «порожнистого» тембру.

2. Фігури Хладні — форма звуку

У 1787 році Ернст Хладні розсипав пісок на квадратній металевій пластині, провів смичком по її краю та виявив, що пісок мігрував із вібруючих ділянок і накопичувався вздовж ліній нульового руху — вузлових ліній. Кожна резонансна частота виробляла інший, часто дивовижно красивий геометричний візерунок. Ці фігури, тепер названі фігурами Хладні, — це прямий візуальний відбиток власних мод вібруючої пластини.

Керівне рівняння (тонка пружна пластина, теорія Кірхгофа-Лява):
  D ∇&sup4;w + ρh ∂²w/∂t² = 0
  w(x,y,t) = прогин пластини
  D = Eh³/[12(1−ν²)] = згинальна жорсткість  (E = модуль Юнга, h = товщина, ν = коефіцієнт Пуассона)
  ρ = густина, h = товщина

Для вільно вібруючої квадратної пластини (сторона a), розділення змінних w = W(x,y)e^(iωt):
  ∇&sup4;W = λ&sup4;W,   де λ&sup4; = ρhω²/D

Просто опертa (SS) прямокутна пластина — точний розв'язок:
  W_mn(x,y) = sin(mπx/a) · sin(nπy/b)
  f_mn = (π/2) · (D/ρh)^(1/2) · [(m/a)² + (n/b)²]
  Вузлові лінії: вертикальні при x = ka/m (k=1..m−1); горизонтальні при y = lb/n (l=1..n−1)

Емпіричне правило Хладні (для квадратних пластин, вільні краї):
  f_mn ∝ (m + 2n)²   ← наближено для вільних граничних умов
  Реальні моди вільної пластини вимагають чисельного МСЕ; вузлові візерунки складніші,
  ніж у випадку простого опирання.

Кількість вузлових ліній:
  Мода (m,n): m−1 вузлових ліній паралельно осі Y, n−1 паралельно осі X
  Або m−1 + n−1 вигнутих вузлових ліній для непрямокутної симетрії.

Виродженість та змішування мод:
  Коли дві моди мають однакову частоту (f_mn = f_nm для квадратної пластини),
  будь-яка лінійна комбінація також є дійсною модою — фактичний візерунок залежить
  від початкових/граничних умов. Це змішування створює повернуті та вигнуті вузлові
  візерунки, що спостерігаються в реальних експериментах Хладні.

Фігури Хладні — це не лише краса: скрипкові майстри використовують їх для перевірки акустичної якості верхньої та нижньої деки. Верхня дека скрипки Страдіварі має конкретні резонансні частоти, які майстер формує й налаштовує, знімаючи деревину, доки фігури Хладні не збігатимуться з бажаними шаблонами. Сучасні конструктори гучномовців використовують метод скінченних елементів для аналізу вузлових візерунків, щоб усунути небажані резонанси в конусах і мембранах.

3. Кімнатні моди — резонанси в прямокутних просторах

Концертний зал, студія звукозапису чи вітальня — це акустичний резонатор. Звукові хвилі, відбиваючись від шести стін, створюють тривимірні стоячі хвилі — кімнатні моди — на конкретних частотах, визначених розмірами приміщення. На цих частотах деякі позиції в кімнаті мають дуже високий тиск (пучності), а інші — майже нульовий тиск (вузли). Слухацький досвід кардинально змінюється залежно від того, де ви сидите.

Модальні частоти (жорсткостінна прямокутна кімната, L×W×H):
  f_lmn = (c/2) · √[(l/L)² + (m/W)² + (n/H)²]
  l, m, n = невід'ємні цілі числа (не всі нулі)
  c = швидкість звуку

Типи мод:
  Осьова (одна ненульова індексна складова):     хвилі лише вздовж однієї осі
    напр. f_100 = c/(2L)         ← найпростіші, найсильніші моди
  Тангенціальна (дві координати):       хвилі вздовж двох осей, що відбиваються від чотирьох стін
  Коса (три координати):        хвилі вздовж усіх трьох осей, що відбиваються від усіх шести стін

Густина мод (кількість мод нижче частоти f):
  N(f) ≈ (4πV/3c³)f³ + (πS/4c²)f² + (L_загальна/8c)f   (формула Вейля)
  V = об'єм кімнати, S = загальна площа стін, L_загальна = сума всіх довжин ребер
  Моди розріджені на низьких частотах, щільні на високих.

Частота Шредера (перехід від модального до статистичного режиму):
  f_S ≈ 2000 · √(RT60/V)   Гц
  Нижче f_S: домінують кімнатні моди; вище f_S: справедливе припущення дифузного поля.
  Типова мала кімната (V=40 м³, RT60=0,4 с): f_S ≈ 200 Гц
  Концертний зал (V=20 000 м³, RT60=2 с): f_S ≈ 20 Гц

Форми мод тиску (жорсткі стіни):
  p_lmn(x,y,z) = A · cos(lπx/L) · cos(mπy/W) · cos(nπz/H)
  Вузлові площини у позиціях x = L(2k+1)/(2l) тощо.

Проблемні інтервали (моди занадто далеко одна від одної):
  Великий інтервал між частотами послідовних мод → нерівномірна басова відповідь.
  Рекомендовані пропорції кімнати (критерій Бонелло):
    L:W:H ≈ 1,6:1,25:1 (або 1,9:1,4:1) — мінімізує збіжні модальні частоти.

Кімнатні моди — головна причина, чому маленькі студії звукозапису звучать погано в басовому діапазоні. Нижче приблизно 200 Гц моди розташовані широко; якщо ви сидите в пучності тиску для осьової моди 80 Гц, бас звучить гулко; у вузлі він зникає повністю. Дизайнери студій розміщують поглиначі, дифузори та басові пастки в кутах (де тиск найвищий для найнижчих мод) і обирають асиметричні розміри приміщення, щоб розподілити модальні частоти.

4. Реверберація концертного залу — RT60 та геометрія кімнати

Час реверберації RT60 — це час, за який звук згасає на 60 дБ після зупинки джерела. Його визначив і вперше виміряв Уоллес Клемент Сабін у 1900 році, вивчаючи акустику Художнього музею Фогга в Кембриджі. Його емпірична формула, отримана з ретельних експериментів із секундоміром та органною трубою, залишається однією з найкорисніших інженерних формул в архітектурній акустиці.

Формула Сабіна (1900):
  RT60 = 0,161 · V / A
  V = об'єм кімнати (м³)
  A = сумарне поглинання (м² сабінів) = Σ α_i S_i
  α_i = коефіцієнт поглинання поверхні i (0 = ідеальне відбиття, 1 = ідеальне поглинання)
  Дійсна для низького середнього поглинання (α_сер < 0,2)

Формула Айрінга (1930, точніша для вищого поглинання):
  RT60 = −0,161 · V / [S · ln(1 − α_сер)]
  S = загальна площа поверхні
  Зводиться до формули Сабіна при α_сер → 0

Середня довжина вільного пробігу (середня відстань між відбиттями):
  λ_mfp = 4V/S   (дійсна для дифузного поля)
  Кількість відбиттів за секунду: c / λ_mfp = cS/(4V)
  Швидкість згасання енергії: η = −cSα/(4V) → RT60 = 24 ln(10)/η ≈ 0,161 V/(Sα)

Оптимальний RT60 за призначенням приміщення:
  Безлунна камера:     < 0,1 с
  Студія звукозапису:     0,2 – 0,4 с
  Мовлення (клас):   0,4 – 0,8 с
  Оперний театр:          1,2 – 1,4 с
  Симфонічний зал:        1,8 – 2,2 с
  Собор:            4 – 10 с

EDT (раннє час згасання):
  Згасання оцінюється від 0 до −10 дБ (екстрапольовано до 60 дБ).
  Корелює тісніше з суб'єктивною «живістю», ніж повний RT60.
  Великий EDT відносно RT60 свідчить про домінування пізніх відбиттів.

Перше застосування формули Сабіна відбулося при реновації театру Сандерс у 1900 році: додавши м'які сидіння й драпірування, він підвищив поглинання й знизив RT60 із незрозумілих 5,5 секунди до приблизно 1,2 секунди. Бостонський симфонічний зал, спроєктований із застосуванням формул Сабіна у 1900 році, залишається одним із акустично найкращих концертних залів світу — прямою спадщиною застосування фізики до архітектури.

5. Механіка равлика — внутрішнє вухо як частотний аналізатор

Внутрішнє вухо виконує неперервне перетворення Фур'є. Равлик — заповнена рідиною скручена трубка завдовжки близько 35 мм — діє як механічний частотний аналізатор: кожна позиція вздовж базилярної мембрани резонує на конкретній частоті, з високими частотами біля основи й низькими — біля верхівки. Волоскові клітини в кожній позиції перетворюють механічний рух на нейронні сигнали, створюючи тонотопічну карту.

Тонотопічна частотна карта (функція Грінвуда, 1990):
  f(x) = A · (10^(a·x) − k)
  x = дробова відстань від верхівки (0 = верхівка, 1 = основа)
  Людський равлик: A = 165,4, a = 2,1, k = 0,88 (дає f у Гц)
  Верхівка:  ~20 Гц (низька частота)
  Основа:  ~20 000 Гц (висока частота)
  ~3,5 мм на октаву вздовж базилярної мембрани

Біжуча хвиля (модель Дьордя фон Бекеші):
  Градієнт жорсткості вздовж базилярної мембрани (жорстка біля основи, м'яка біля верхівки)
  змушує вхідні хвилі тиску рухатися від основи до верхівки.
  Кожна хвиля наростає до максимальної амплітуди в характерному місці для своєї частоти,
  а потім швидко згасає — створюючи різкий пік, що кодує частоту.

Крива механічного налаштування:
  Смуга пропускання резонансу базилярної мембрани: Q_мех ~ 3–10 (пасивна)
  Активне підсилення (зовнішні волоскові клітини, соматична рухливість OHC):
    OHC підсилюють рух, загострюючи налаштування до Q ~ 30–100
    Цей «кохлеарний підсилювач» пояснює 40 дБ підсилення та тонку частотну
    дискримінацію здорового вуха.

Двотонове придушення:
  Тон на f_1 може придушувати відповідь на f_2, якщо частоти достатньо близькі.
  Відображає стискаючу нелінійність підсилювача OHC.
  Джерело маскувальних ефектів, що використовуються у аудіокомпресії (MP3, AAC).

Критична смуга пропускання (шкала Барка):
  Частотна роздільна здатність равлика; тони в межах однієї критичної смуги маскують один одного.
  Δf_КС ≈ 25 + 75[1 + 1,4(f/1000)²]^0,69 Гц (формула Цвікера)
  Аудіокодеки використовують це: відкидають замасковані компоненти нижче порогу маскування.

Дьордь фон Бекеші отримав Нобелівську премію з фізіології 1961 року за демонстрацію біжучої хвилі у видалених людських равликах. Кохлеарний підсилювач — активна роль зовнішніх волоскових клітин — був відкритий у 1980-х роках і пояснює, як здорове вухо досягає одночасно надзвичайної чутливості (виявлення вібрацій, менших за діаметр атома водню) та тонкої частотної вибірковості. Кохлеарні імпланти обходять пошкоджені волоскові клітини, безпосередньо стимулюючи слуховий нерв на декількох позиціях електродів вздовж тонотопічної карти, використовуючи той самий просторовий частотний код, що базилярна мембрана створює природно.

6. Акустичні дифузори — розсіювання звуку в просторі

Плоска стіна відбиває звук як дзеркало: вся енергія повертається в одному дзеркальному напрямку. Дифузор на основі квадратичних лишків (QRD) — періодичний масив колодязів різної глибини — рівномірно розсіює вхідний звук по багатьох кутах. Глибини обираються з послідовностей квадратичних лишків — теоретико- числової конструкції, відкритої Манфредом Шредером у 1975 році. Результат: поверхня, що виглядає акустично шорсткою на всіх частотах у межах своєї проєктної смуги пропускання, усуваючи акустичне «сяйво» та тремтливе відлуння.

Послідовність квадратичних лишків (період N, N — просте):
  s_n = n² mod N   (n = 0, 1, …, N−1)
  Приклад N=7: s = [0, 1, 4, 2, 2, 4, 1]
  Кожна послідовність симетрична: s_n = s_{N-n}

Глибина колодязя (максимальна глибина d_max для найнижчої проєктної частоти f_min):
  d_n = s_n · λ_min / (2N)   або   d_n = s_n · c / (2N · f_min)
  Загальний діапазон глибин: d_max = (N−1) · c / (2N · f_min)
  Для N=7, f_min=500 Гц: d_max = 6 · 343/(14·500) ≈ 294 мм

Розсіяний тиск у дальній зоні (наближення Фраунгофера, M періодів):
  p(θ) = Σ_{n=0}^{N−1} r_n · exp(iknd sinθ)   × sin(Mk d_x sinθ/2)/sin(k d_x sinθ/2)
  k = 2πf/c, d_x = ширина колодязя, θ = кут розсіювання
  Рівномірне |p(θ)|² по θ — мета проєктування (рівна енергія в усіх пелюстках решітки)

Чому це працює (теорія чисел):
  ДПФ коефіцієнтів відбиття QRD r_n = exp(iπs_n/N) має сталу величину:
  |ДПФ{r_n}| = стала для всіх частот → плаский кутовий спектр потужності.
  Це зв'язок між сумами Гаусса (теорія чисел) та акустичною дифузією.

Коефіцієнт дифузії d:
  d = [(NΣE_i² − (ΣE_i)²) / ((N−1)(ΣE_i)²)]^(1/2)   (стандарт AES)
  d = 1 для ідеальної дифузії, d = 0 для дзеркального відбиття.
  QRD зазвичай досягає d > 0,8 у межах своєї проєктної смуги пропускання.

Дифузори Шредера тепер є стандартом у студіях звукозапису, концертних залах та мовних установах по всьому світу. Оскільки їхній дизайн вкорінений у теорії чисел, різні прості числа N виробляють різні ширини періодів і частотні діапазони. Реальні установки поєднують кілька періодів поруч; пелюстки решітки від періодичності додаються конструктивно, даючи сильне дифузне розсіювання. Зв'язок між чистою математикою (квадратичні лишки, суми Гаусса) та акустикою приміщень — один із найелегантніших прикладів застосування абстрактної математики в інженерії.

Ширший світ акустики

Ці шість симуляцій охоплюють акустику від масштабу отвору труби до концертного залу, зачіпаючи одну з найбезпосередніше сприйнятих фізик у світі. Звук всюди, і хвильове рівняння, що ним керує, повторюється в усій фізиці: та сама математика описує електромагнітні хвилі, квантово-механічні амплітуди ймовірностей, сейсмічні хвилі в Землі та гравітаційні хвилі в просторі-часі. Розуміння того, як межі, геометрія та властивості матеріалу формують акустичні хвилі, будує інтуїцію, яка безпосередньо переноситься на всі ці домени.