Как Хэмминг изучал цифровые фильтры
Хэмминг пришёл к цифровым фильтрам как математик, а не как инженер-электротехник. Когда он спросил инженеров, почему они используют синусоиды вместо полиномов или функций Бесселя, никто не дал удовлетворительного ответа. Поэтому он вернулся к основам.
Он выявил три независимые причины, по которым комплексные экспоненты доминируют в цифровой обработке сигналов. Каждая причина отдельно оправдывает этот выбор; вместе они делают его практически обязательным.
Причина 1: Стационарность по времени
Большинство систем обработки сигналов не имеют естественного начала отсчёта времени. Фильтр, применённый в полдень, должен вести себя идентично тому же фильтру, применённому в полночь. Это ограничение стационарности по времени вынуждает собственные функции быть комплексными экспонентами.
В математических терминах: если линейная стационарная система имеет вход x(n) = e^{i2πfn}, выход также должен колебаться на частоте f. Только комплексные экспоненты удовлетворяют этому условию.
Причина 2: Линейность
Линейные системы подчиняются принципу суперпозиции. Собственные функции любого линейного оператора — это функции, которые остаются неизменными (только масштабируются) при действии оператора. Для оператора сдвига S: x(n) → x(n−1) собственные функции — это ровно e^{i2πfn}.
Причина 3: Теорема Найквиста
Если непрерывный сигнал не содержит частотных компонент выше f_max, дискретизация его на частоте ≥ 2f_max сохраняет всю информацию. Теорема Найквиста–Шеннона о дискретизации связывает непрерывную и дискретную обработку сигналов чистым образом только для представлений Фурье.
При равномерной дискретизации одна высокая частота отображается в одну более низкую частоту. При полиномиальных базисах одна высокая степень t отображается во множество более низких степеней: беспорядок, который Фурье полностью избегает.
Передаточная функция как собственное значение
Когда e^{i2πfn} входит в линейный стационарный фильтр, выход равен H(f) · e^{i2πfn} для некоторого комплексного числа H(f). Фильтр масштабирует и сдвигает колебание, но не может изменить его частоту.
H(f) собирает всё поведение фильтра на частоте f в одно комплексное число. Для фильтра с коэффициентами c_k:
H(f) = Σ c_k · e^{−i2πfk}
Эта формула делает H(f) преобразованием Фурье последовательности коэффициентов. Каждый частотный канал работает независимо. Фильтр разлагает вход на частотные компоненты, умножает каждую на H(f) и переассемблирует их.
Теорема о дискретизации
Хэмминг заметил, что теорема Найквиста о дискретизации была известна до Найквиста, но ему досталась слава. Он цитировал Пастера: «Удача улыбается подготовленному уму». Человек, который связывает существующую идею с практической необходимостью, получает известность.
Теорема
Если непрерывный сигнал x(t) не содержит частотных компонент выше f_max, то дискретизация его на частоте f_s ≥ 2·f_max сохраняет всю информацию. Исходный сигнал точно восстанавливается из образцов.
Порог f_s / 2 = f_max носит имя Найквиста. Дискретизация ровно на частоте Найквиста (2·f_max) в принципе достаточна, но опасна на практике: любое малое несовпадение приводит к наложению высшей частоты.
Наложение спектра
Когда сигнал содержит частоты выше f_s/2, те частоты складываются в полосу [0, f_s/2]. Синусоида на частоте f = f_s/2 + δ выглядит неотличимой от одной на f_s/2 − δ. Туки придумал термин наложение спектра для названия этого выдачи себя за другую.
Геометрическая картина: комплексные экспоненты на частотах f и f + f_s дают идентичные образцы в целые моменты времени. Они имеют общий алиас.
Выбор частоты дискретизации
Практическая система цифрового аудио должна выбрать частоту дискретизации перед разработкой фильтров. Люди слышат примерно до 20 кГц. Стандартная частота дискретизации CD 44,1 кГц устанавливает частоту Найквиста на 22,05 кГц.
До дискретизации фильтр подавления наложения должен удалить все частоты выше частоты Найквиста. Если даже небольшая компонента 25 кГц попадёт в дискретизатор, она наложится на 44100 − 25000 = 19,1 кГц — слышимо.
Три ограничения аппаратного обеспечения
Хэмминг преподавал более широкий урок наряду с математикой. Цифровые фильтры существуют потому, что аппаратное обеспечение имеет ограничения — и понимание этих ограничений формирует хороший дизайн.
Он выявил три закона природы, которые ограничивают производительность аппаратного обеспечения:
1. Молекулярный размер: схемы не могут бесконечно сжиматься. Ниже определённого масштаба доминируют квантовые эффекты.
2. Скорость света: сигналы путешествуют максимум со скоростью 3×10⁸ м/с. Тактовые циклы быстрее времени прохождения света по чипу приводят к сбоям.
3. Рассеивание тепла: переключение потребляет мощность, которая становится теплом. Плотные, быстрые чипы перегреваются без охлаждения.
Его философия проектирования следовала напрямую: поймите ограничения, затем проектируйте системы, которые удобно работают в их рамках, с запасом на вариацию и ошибку.
Цифровые фильтры переносят вычисления с аппаратного обеспечения (аналоговые схемы) на программное (арифметика на образцах). Этот переход торгует хрупкостью аппаратного обеспечения на численную точность и программируемость — следствие теоремы о дискретизации, не чудо.
Философия проектирования Хэмминга
Фреймворк Хэмминга: цифровой фильтр реализует в программном обеспечении то, что делает аналоговый фильтр в аппаратном обеспечении. Теорема о дискретизации — мост. Как только вы знаете, что мост крепок, вы можете проектировать фильтры, указав желаемую передаточную функцию H(f), затем найдя последовательность коэффициентов, которая её реализует.
Работа инженера становится спецификацией и арифметикой, а не намоткой индукторов и пайкой конденсаторов.