Ряди Фур’є: Постійні Фазори та Явище Гіббса

Як перетворити хвилю на суму кіл

У 1822 році Джозеф Фуріє зробив заяву, яка шокувала математиків того часу: будь-яка періодична функція, хоч би яку не була, може бути записана як сума звичайних синусів та косинусів. Квадратна хвиля, зубчаста хвиля, крива тиску струменини, що виникла при ударі гітарного струнника – всі вони розкладаються на основну частоту плюс сходинка гармонік у цілих кратностях її. Найчистіший спосіб зрозуміти це — не думати про синуси та косинуси окремо, а замість цього бачити кожен гармонічний як невелику стрілку – фазор, що обертається з постійною швидкістю навколо початку комплексної площини. Додайте кінчики достатньої кількості обертових стрілок, вкладених одна в одну, і сума простежує оригінальну форму хвилі з часом.

Ряди коефіцієнтів: як знайти амплітуду кожного гармонічного

Для функції f(t) з періодом T, амплітуда кожного гармонічного визначається шляхом проектування f на синусоїдний або косинусоїдний гармонічний — внутрішній добуток, який вимірює ступінь перекриття двох форм. Використовуючи тотожність Ейлера e^{iθ} = cos θ + i sin θ, весь ряд стискається в одну компактну комплексну експоненціальну форму, де одне комплексне число cₙ несе як амплітуду, так і фазу гармонічного n:

f(t) = Σ_{n=-∞}^{∞} c_n · e^{i·2πnt/T} c_n = (1/T) ∫₀ᵀ f(t) · e^{-i·2πnt/T} dt // проектування f на гармонічний n |c_n| = амплітуда n-го фазору arg(c_n) = його кут відносно початку (фаза) Кожен член c_n·e^{i2πnt/T} є фазором, що обертається n разів швидше за фундаментальну частоту. З'єднуйте їх кінці до кінців — центр другого фазору знаходиться на кінчику першого, центр третього - на кінчику другого і так далі — і в будь-який момент часу t положення кінчика останньої стрілки точно дорівнює f(t). Це саме те, як працюють класичні «епіцикли Фур'є»: достатньо заглиблених кіл може відтворити хвилю, зірку або портрет.

f(t) = Σ_{n=-∞}^{∞} c_n · e^{i·2πnt/T}

c_n = (1/T) ∫₀ᵀ f(t) · e^{-i·2πnt/T} dt   // project f onto harmonic n
|c_n|   = amplitude of the n-th phasor
arg(c_n) = its starting angle (phase)

Будівництво хвилі-квадрата, по одній гармоніці за раз

Квадрова хвиля містить лише непарні гармоніки, кожна з яких зменшується вчетверо від амплітуди фундаментальної: додавання sin(t), потім sin(3t)/3, потім sin(5t)/5 і так далі, поступово заострюючи округлену синусову хвилю в більш плоскі вершини та крутіші краї. Хвиля типу тетрада використовує всі гармоніки (парні та непарні) з амплітудою 1/n та змінюваним знаком, а трикутна хвиля використовує лише непарні гармоніки, які швидко спадують, як 1/n², тому трикутні хвилі звучать "чистіше", ніж квадратні або тетрадні – їхні вищі гармоніки значно тихіші.

Зміна фази в системі з обмеженою енергією

Спостерігайте за квадратною хвилею уважно біля вертикального стрибка та відбувається щось дивне: часткова сума не просто наближається до стрибка, а перекошується повз верхню точку та недостатньо підходить до нижньої точки, коливається перед тим, як осередитися деінде. Додавання більших гармонік стискає цей ривок ближче і ближче до розриву — але висота першого перекосу завгодно становить приблизно 9% від розміру стрибка, незалежно від того, скільки термінів ви додаєте. Це явище Гіббса, суворо пояснене Джошіа Ховардом Гіббсом у 1899 році: безперервна сума гладких синусоїд може наближатися до стрибка будь-якої міри, але ніколи не може повністю усунути ривок саме на краю. Це математична основа "ривків", які ви бачите навколо гострих країв у зображеннях JPEG та чуєте як попередні ехо в агресивно стиснутому аудіо.

Від серії до трансформації: дозволяємо періоду піти в нескінченність

Серія Фур'є має сенс лише для періодичного сигналу — такого, що повторюється назавжди з періодом T. Більшість цікавих сигналів (один удар, фотографія, ціна акцій) ніколи не повторюються. Щоб їх обробити, уявіть, як розтягуєте T до нескінченності: відстань між сусідніми гармоніями, 1/T, сходиться до нуля, і дискретна драбина коефіцієнтів c_n зливається в безперервну криву. Сума стає інтегралом, а дискретні коефіцієнти c_n стають безперервною функцією F(ω) — перетворення Фур'є:

F(ω) = ∫_{-∞}^{∞} f(t) · e^{-iωt} dt // пряме перетворення (T → ∞) f(t) = (1/2π) ∫_{-∞}^{∞} F(ω) · e^{iωt} dω // обернене перетворення Де серія Фур'є дає вам спектр, зроблений з дискретних шипів у 0, ±1/T, ±2/T…, перетворення Фур'є дає безперервну спектральну щільність — будь-яка частота може з’явитися, з незліченною кількістю енергії в кожній.

F(ω) = ∫_{-∞}^{∞} f(t) · e^{-iωt} dt      // forward transform (T → ∞)
f(t) = (1/2π) ∫_{-∞}^{∞} F(ω) · e^{iωt} dω // inverse transform

Часті запитання

Спробуйте наживо

Усе, що вище, працює прямо у вашому браузері — відкрийте the simulation і змінюйте параметри під час роботи. Нічого не встановлюється, нічого не завантажується на сервер, уся модель живе в одній вкладці.