|
|
Вот SM меня снабдил основными инструментами для разработки фильтров с произвольными частотными характеристиками. Что получилось?
f = [0 0.2 0.49 0.51 0.8 1 ];
m = [1 0.8 0.5 0.7 0.9 0 ];
b = fir2(64,f,m)
M = max(abs(b))*1.001
%b = b ./ M
Hq = qfilt('fir',{b});
Hq.CoefficientFormat = quantizer('fixed', 'round', 'saturate', [16 15]);
[h,w] = freqz(Hq,128);
plot(f,m,w/pi,(abs(h)));
%plot(f,m,w/pi,20*log10(abs(h)))
Что же меня смущает? Как человек морально устаревший еще в прошлом тысячилетии я много читал... Ну там Основы Теории Цепей, незвестного ныне широкой публике Гоноровского... Так вот столь отвлеченное от Идей Цифровой Обработки Сигналов чтиво подвигло меня на следующие размышления:
Мне, почему-то, очень не нравятся "облизывание" модулем частотной характеристики фильтра прямых отрезков. Что же именно мне не нравиться? Увы, это только интуитивные ощущения. Например, аппроксимация прямоугольного фильтра по-Баттерворту мне милее элиптических фунций Якоби. А вот Гаусов фильтр совсем здорово.
Другими словами, за этим "вылизыванием" я ощущаю появление лишних переходных искажений, что может быть вредным в моменты звуковой атаки, и уж точно приводит к излишней длине импульсной характеристики.
У кого какие мнения?
Вопрос к SM. Сэр, неисключено, что Вы по-молодости увлекались древними окультными радиотехническими практиками, и (чем черт не шутит?) даже перевернули пару страниц из ОТЦ-буки. Поэтому, не подскажите ли неких функций в матлабе, которые также работали бы с произвольной частотной характерисикой. НО!
В fir2 можно задать характерискику из пяти точек, и из пятисот. При этом она может попытаться выдать результат для фильтра первого порядка. Мне кажется надо исходить в данном случае из чисто импирического положения 1 точка - один порядок фильтра. Или же иначе - вести оптимизацию коэффициентов используя ТОЛЬКО ЗАДАННЫЕ точки, а не точки, которые автоматически беруться из "отрезков". Я эти точки НЕ ЗАДАВАЛ. А может есть способ задать функцию-образец не прямыми отрезками, а Гауссианами? А может есть другие мысли?
Прошу пардону за расплывчатость, пока еще сам не осознал чего мне надо, а руководство торопит. :)))
E-mail: 
info@telesys.ru
