Поэтому я и написал за желание вкурить это z-преобразование. Острой необходимости пока что не было, но мало ли когда она возникнет?Однако справедливости ради отмечу, что эти алгоритмы фильтрации не являются универсальными.
Каноническая форма БИХ фильтра 10-го порядка требует порядка 16-ти операций, тогда как моя реализация требует 20. Плюс ещё и архитектура у процессоров заточена под умножение с накоплением (количество операций можно уменьшить вдвое), чего в моей реализации не учитывается. Правда, они не будут работать с целочисленной системой, т.к. всё равно делить результат, чтоб влез...Формула требует 5 умножений.
Да как-то я стороною обхожу ресурс DI HALT'а (это ж его ресурс?). Там много интересного можно найти, верно, но для меня открытием явилось то, что там сидят не только просто радиолюбители и просто программисты, которые умеют решать лишь стандартные задачки...Я бы Вам посоветовал на http://we.easyelectronics.ru/ статью написать на эту тему. Успех гарантирован !
Кстати, там же можно обсудить модель динамической головки. Аудиофилов там полно...
На практике RC фильтры часто нужны для сильного подавления шума АЦП не критичных ко времени сигналов в слабых МК. Термодатчик в кристалле AVR, смещение дифференциального усилителя и т.п. Поэтому там берутся вообще доли герца, часто RC второго порядка (он эффективнее остальных при равных затратах).Через чур низкочастотный фильтр получается. При Fsmpl=20 кГц частота среза 17.6 Гц
У меня где-то был кроссовер на фильтрах Бесселя 10-го порядка с фазовращателями, дающий на выходе синфазный сигнал на всех полосах, который я затем смешивал обратно и сравнивал с исходным сигналом. Под трёхполосную АС. Я когда-то его разрабатывал и тестировал в саунфордже, прогоняя через него 192 килосемпловые WAV файлы (частота дискретизации самого усилителя после ресемплинга). Если вам так интересно поковыряться в этом коде, то я предлагаю подождать, пока я откопаю и причешу тот фильтр. Так вот, в нём раздел идёт на границах 200Гц и 2кГц. И в принципе проблем нет с повышением частоты, срезом всего что выше 20кГц и ниже 10Гц и т.п. Арифметика целочисленная. Там я, кстати, заметил что ПК гораздо охотнее жуют плавающую точку (раза в 2 быстрее за фиксированную).
Я тестировал в звуковом редакторе на тестовых сигналах, сгенерированных там же (ну и на реальных музыкальных файлах). На мой взгляд, куда нагляднее и гарантирует отсутствие ошибок.Зная Н(z) легко построить АЧХ и ФЧХ.