* Данный текст распознан в автоматическом режиме, поэтому может содержать ошибки

1160 Аналитическое и спектральное моделирование На ФЧХ фильтра можно заметить характерный разрыв, связанный с превы шением фазовым углом граничного значения ?. Такой способ представления фа зового сдвига общепринят, поскольку его изменения стремятся вписать в диапа зон от –? до ?. 15.3.4. Проектирование цифрового фильтра Основной недостаток аналоговых активных фильтров, подобных описанному выше, заключается в их малом порядке. Его повышение, за счет применения мно гих звеньев низкого порядка, ведет к значительному повышению габаритов филь тров и их стоимости. От этого недостатка свободны современные цифровые филь тры, число ячеек которых N даже при однокристальном исполнении может достигать десятков и сотен. Это обеспечивает повышенную частотную селекцию. Спроектируем фильтр (N+1) го порядка класса FIR (Finite Impulse Response, или с конечной импульсной характеристикой). Каждая из N ячеек временной за держки фильтра удовлетворяет следующей зависимости выходного сигнала y от входного x вида: Подключим пакет расширения plots, нужный для графической визуализации проектирования: > restart:with(plots): Warning, the name changecoords has been redefined Зададим исходные данные для проектирования полосового цифрового фильт ра, выделяющего пятую гармонику из входного сигнала в виде зашумленного ме андра с частотой 500 Гц: > > > > > N := 64: fs := 10000: fl := 2300: fh := 2700: m := 10: # # # # # ?enei naeoee oeeuo?a (ia 1 iaiuoa ii?yaea oeeuo?a) ?anoioa eaaioiaaiey Ie?iyy a?aie?iay ?anoioa Aa?oiyy a?aie?iay ?anoioa 2^m > N – ?enei oi?ae aey aiaeeca Вычислим: > T := 2^m-1; T := 1023 > F1 := evalf(fl/fs); F1 := .2300000000 > F2 := evalf(fh/fs); F2 := .2700000000 > Dirac(0) := 1: # Ooieoey Ae?aea > fp1:=2*Pi*F1: fp2:=2*Pi*F2: Зададим характеристику полосового фильтра: > g := (sin(t*fp2)-sin(t*fp1))/(t*Pi);