Фур’є Епіцикли: Малювання з Обертовими Циклами

Шлях – це складний знаково-векторний сигнал

Візьміть будь-яку закриту криву, яку ви можете намалювати ручкою — силует кота, вашу підпис, химерну зірку — і відтворіть навколо неї N точок з рівненькими кроками часу. Розглядайте кожен пункт (x, y) як окремий комплексний номер x + iy. Таким чином, у вас виникне періодичний, знаково-векторний сигнал довжиною N, і будь-який періодичний сигнал, незалежно від того, наскільки він нерівномірний, може бути точно записаний як сума комплексних експонент: дискретна ряду Фур’є.

Кожен термін є обертовим колом

Дискретна швидкісна перетворення (ШП) розкладає траєкторію на N комплексних коефіцієнтів ck, по одному для кожного цілого частотного k (включно з від'ємними). Кожен коефіцієнт є комплексним числом, отже має величину та фазу — і термін ck·e^(i·k·2πt/N) відслідковує точно коло радіуса |ck|, що обертається з кутовою швидкістю пропорційною k, починаючи з кутового аргументу arg(ck):

ck = (1/N) · Σ_{n=0}^{N-1} траєкторія[n] · exp(-i · 2π · k · n / N) для k = -N/2 .. N/2 reконструкція(t) = Σ_k ck · exp(i · k · 2π · t / N) ↑ коло k: радіус |ck|, швидкість k, кут відліку arg(ck) З’єднайте кола кінець до кінця — центр кола k+1 знаходиться на кінчику кола k — і, прокотившись один повний період, кінець останнього кола точно відтворює ваш оригінальний малюнок. Відсортуйте кола за |k| так, щоб найбільше, найповільніше коло малювалося першим, і кожне наступне менше, швидше коло уточнює форму, і ви отримуєте класичну анімацію: ланцюг вкладених коліс, зовнішній кінець якого відтворює зображення.

c_k = (1/N) · Σ_{n=0}^{N-1} path[n] · exp(-i · 2π · k · n / N)   for k = -N/2 .. N/2

reconstruction(t) = Σ_k  c_k · exp(i · k · 2π · t / N)
                        ↑ circle k: radius |c_k|, speed k, start angle arg(c_k)

Як багато кіл вам дійсно потрібно?

З усіма N колами реконструкція точна — це просто інвертована трансформація, а не наближення. Цікава поведінка виникає при меншій кількості: для плавної, злегка вигнутої лінії коефіцієнти ck зазвичай швидко спадають у міру збільшення |k|, тому лише кілька низькочастотних кіл вже достатньо для захоплення загальної форми та залишається лише тонке тремтіння. Для малюнка з гострими кутами коефіцієнти повільно спадають, оскільки представлення розриву у нахилі гладкими косинусами в своїй основі потребує високих частот — така ж причина, чому серія Фурієвого розкладу прямокутного сигналу потребує нескінченної кількості членів для ідеального відтворення гострого краю, а кільця з незначним перебільшенням, відомим як явище підвищення Гіббса, біля кута, якщо ви зупиняєтеся занадто рано.

Чому від'ємні частоти мають значення

Для дійсної сигналу комплексні коефіцієнти Фур’є при +k та -k є комплексно-спряженими один до одного, і їх парне поєднання дає звичну синус-косинусну ряду Фур’є з першого курсу. 2D траєкторія відрізняється: x + iy справді є комплексним числом, а не дійсним, тому коефіцієнти +k та -k незалежні і обидва несуть реальну інформацію – кола +k обертаються в одному напрямку, а кола -k – в іншому. Видаліть кола з від’ємною частотою, і зазвичай закрита 2D форма дегенерує у форму, яка лише періодично простежує еліпс або лінію, оскільки ви втратили половину інформації, необхідної для порушення симетрії обертання.

Від Птолемеєвих епіцентрів до обробки сигналів

Візуалізація кіл, що обертаються навколо інших кіл, існувала задовго до аналізу Фурієром — птолемеєві астрономи використовували епіцентри та конічні перетини для наближеного моделювання руху планет навколо нерухомої Землі, додаючи більше кіл, коли потрібно було покращити відповідність. Сучасна версія Фурієра використовує ту ж геометричну ідею, але з математичним обґрунтуванням: замість ручного налаштування кіл для узгодження спостережень, DFT (перетворення Фур’є) обчислює точний радіус, швидкість та фазу кожної спіралі безпосередньо з даних. Одна й та сама техніка — розкладання періодичного сигналу на суму обертових фазорів — є основою для цифрових фільтрів, стиснення MP3 та JPEG, а також практично для всіх аналізаторів спектру.

Frequently asked questions

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

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