Роз’яснення спектру Фур'є: Чому швидкий перетворення……

Дискретна Фур’єрова Перетворення та чому воно повільне

DFT перетворює N часових зразків на N частотних коефіцієнтів, X[k] = Σ x[n]·e^(-j2πkn/N), де кожен вихідний коефіцієнт X[k] вимірює, наскільки частота k/N присутня в сигналі. Обчислення цього безпосередньо потребує підсумовування N членів для кожного з N вихідних значень – N × N = N² комплексних множень. Для N = 1024 це більше одного мільйона операцій; до існування FFT, 1024-точкове DFT на апаратному обладнанні 1960-х років займало декілька хвилин.

Розділяй та володарюй: парна-непарний поділ

Інсайт Коулі та Тьюкі (1965 року) (передбачений Гауссом близько 1805 року) полягав у розділенні N-точкового DFT на дві N/2-точкові DFT, одну над парними індексами та одну над непарними індексами: X[k] = E[k] + W_N^k·O[k]. Оскільки E[k] і O[k] є періодичними з періодом N/2, кожен результат повторно використовується як для k, так і для k+N/2 через оновлення бабочки X[k+N/2] = E[k] − W_N^k·O[k]. Рекурсивне повторення цього розділення дає рекурентну формулу T(N) = 2T(N/2) + O(N), яка розв’язується до O(N log N) — для N = 1024, приблизно 10 000 операцій замість мільйона, що є прискоренням у 100 разів, яке стає ще більш драматичним при більших N.

жива демонстрація · пов'язана симуляція● LIVE

Показники та обернення біт

Термін обертання W_N^k = e^(-j2πk/N) називається показником «Twiddle». Його ключова симетрія, W_N^(k+N/2) = -W_N^k, означає, що кожен «перовий» елемент повторно використовує одне комплексне множення для виробництва двох виходів – зменшуючи кількість множень ще вдвічі, крім заощад від розділяй та володарюй. Оскільки рекурсивний поділ парних/непарних чисел природним чином перемішує порядок вихідних даних, перш ніж запускати «перових» стадій, впроваджується перестановка біт-обернення – замінюючи кожен індекс його бітовими цифрами у зворотному порядку – щоб «перових» стадії могли працювати на місці та все одно забезпечувати природно відсортований спектр на виході.

N (samples)   Naive DFT (N²)   FFT (N log₂N)   Speedup
64            4,096            384             ~10×
1,024         1,048,576        10,240          ~100×
65,536        4,294,967,296    1,048,576       ~4,000×

Де синхронізація з частотним спектром

Аудіо обробка виконує FFT безперервно: вікно N вибірок, перетворення, коригування бін частоти для еквалайзера або шумозаглушення, потім інверсія та перекриття-додавання — аудіо в реальному часі на частоті 44.1 кГц зазвичай використовує N = 512-4096 оновлюється десятки разів на секунду. Завдяки теоремі згортки, ланцюг FFT-множення-інверсії-FFT перетворює O(N²) прямий згорт на O(N log N), основу для ревербування згортки та швидкого множення поліномів. Стиснення блоків 8x8 JPEG використовує пов'язаний лише косинусний перетворення, DCT. OFDM бездротовий зв’язок — Wi-Fi, 4G/5G, DAB радіо — використовує обернену FFT для перетворення символів QAM у частотному діапазоні в трансиверовну форму хвилі, з N варіюється від 64 (Wi-Fi) до 4096 (5G NR).

Frequently asked questions

Чому невмілий DFT є O(N²) і як FFT униклює цю вартість?

Обчислення кожного з N частотних вихідних значень безпосередньо потребує підсумовування N вхідних членів, що дає загалом N помножене на N = N квадрат, тобто множення. Алгоритм швидкого перетворення Фур'є (FFT) Коулі-Тукі рекурсивно розділяє перетворення на парні та непарні індексні частини, повторно використовуючи спільні результати за допомогою періодичних твідових факторів, що зменшує загальну роботу до O(N log N).

Що таке твідова фаза?

Твідова фаза W_N^k = e^(-j2πk/N) - це комплексний поворотний член, який множить зразки на кожному етапі FFT. Її періодичність і 180-градусна симетрія, W_N^(k+N/2) = -W_N^k, дозволяють кожній операції «бабочки» повторно використовувати одне множення для виробництва двох виходів, що є ключовою хитрощаю алгоритму щодо швидкості.

Чому FFT потребує перестановки за бітовим зворотним порядком?

Алгоритм Коулі-Тукі, який працює безпосередньо, природно змішує порядок виводу під час рекурсивного розділення парних та непарних зразків. Попередній перестановки вхідних даних шляхом обернення бітних цифр кожного індексу дозволяє «бабочковим» стадіям працювати безпосередньо і все ще виробляти спектр у природному порядку наприкінці.

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

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

▶ Відкрити симуляцію the simulation

Роз’яснення спектру Фур'є: Чому швидкий перетворення Фур'є змінило світ