Здравствуйте В общем, сабж Тут делаю прибор - фонокардиограф. Там надо обработать сигнал фильтрами Мааса-Вебера (ниже). Динамический диапазон звука с человеческого сердца очень велик, для удешевления прибора в 1998 году вместо 6 каналов сделан один, на котором есть физический фильтр ВЧ (No1) (который фактически сжимает диапазон) +
12бит АЦП. Затем на компьютере восстанавливается сигнал "антифильтром" (с АЧХ, обратной АЧХ ФВЧ No1) и обрабатывается выбранным врачом фильтром из набора Мааса-Вебера. Но тогда (1998) я не умел разрабатывать аналоговые фильтры (сейчас немного умею), мне разработал коллега (фиговый правда, самовозбуждался, хорошо хоть на частоте, которая не используется и в сердце человека практически =0). Цифровые разработал я разработал по книге Уинтце-Ш*. Правда только Батерворта смог, вместо Бесселя, который хотелось. Процесс я вроде бы записал, но вот не получается. Затем разработку новой версии отдали другим людям, мне только фоно оставили. ФВЧ No1 я пересчитал, сначала в ручную, затем проверил Ворк бенчом и в железе - в общем хороший получился. Антифильтр я пересчитал, получился замечательно. Но вот с фильтрами Мааса-Вебера - проблема. Где-то глуплю. Я б использоавл старые, но новые разработчики поменяли частоту оцифровки - вместо 4кГц сделали 2кГц. В общем, что надо сделать цифровой фильтр с характеристикой: на частоте
140 Гц иметь подавление 20 Дб и наклон 18 Бд на октаву. Частота оцифровки 2000 Гц. Желательно Бессела, затем Батерворт, а там что получится. Как делал я.
1) Рассчитал частоту среза и подавления 6 дб-240 Гц, 96 Дб-7,5 Гц (96 Дб, я понимаю, при 12 бит нонсенс, но это надо только для рассчёта)
2) определил передаточную характеристику аналогового фильтра и нашел А,B,C,D,E,F функции K(p)
3)перевел K(p) в вид K(z)
4) перевёл - не правильно - получил 110 Гц и 13 Дб/октава
Пробовал разные методики, например, из требуемых характеристик нашел характеристики аналогового прообраза, затем нашёл характеристики унифицированного ФНЧ, нашел его K(p), перевел в K(z) - все равно не то - 100 Гц и 12 Дб. Кто-нибудь може рассчитать характеристики требуемого фильтра и рассказать как это делал? Ибо я зашёл в тупик (блин, ну как я тогда за 2 дня всё рассчита - не пойму, хоть и не совсем точно, но в пределах). У Мааса-Вебера фильтров 6 штук, остальные я сам постараюсь поборать, к их точности требования по ТУ ниже.
С уважением, Шаповалов Алексей