Сигнальные процессоры. Аксенов В.П. - 93 стр.

необходимости увеличения порядка n. Тогда порядок фильтра подбирается
вручную изменением параметра с n+1 на n+2 и т.д.
     » b = fir1(n, Wn, ftype, kaiser (n+2, beta), 'noscale');
     Проверяем влияние округления коэффициентов b из-за конечной разряд-
ности процессора ADSP-2189M
     » [h,f]=freqz(b);                % расчет АЧХ фильтра
     » subplot(1,2,1)
     » plot(f/pi,20*log10(abs(h)))       % 1-ый график АЧХ c коэфф-ми b
     » grid on
     » bq=round(b*65536)/65536;           % округление коэфф-ов b (формат 1.15)
     » [hq,f]=freqz(bq);
     » subplot(1,2,2)
     » plot(f/pi,20*log10(abs(hq)))       % 2-ой график АЧХ c коэфф-ми bq
     » grid on
     Сравнение графиков на рис. 48 показывает, что амплитудно-частотные
характеристики фильтров практически совпадают, дискретностью младшего
разряда можно пренебречь. Формат чисел процессора 1.15 не оказывает влия-
ния на частотную характеристику фильтра. Структурная схема синтезированно-
го фильтра приведена на рис. 49, его уравнения имеют вид
    y(k) = b0 x(k) + b1 x(k – 1) + b2 x(k – 2) + ··· + b37 x(k – 37) + b38 x(k – 38).
     Выполним синтез режекторного БИХ-фильтра с такими же исходными
данными. Функция ellipord определяет порядок эллиптического фильтра n,
функция ellip возвращает его коэффициенты b и a
     » f1 = 2.0*V; f2 = 3.0*V         % f1,f2- частоты входного сигнала
     » f3 = 1.15*f2; f4=1.55*f2; f5 = 1.90*f2; f6 = 2.4*f2;         % точки АЧХ РФ
     » w3 = 2*f3/Fs; w4 = 2*f4/Fs; w5 = 2*f5/Fs; w6 = 2*f6/Fs;
     % нормированные частоты




                                           93