Отправлено -=ВН=- 25 ноября 2004 г. 12:52 В ответ на: Цифровой фильтр отправлено Vovanchik 25 ноября 2004 г. 12:02

Фильтр БИХ.
2-го порядка.
Разностное уравнение:
y(n)=b0*x(n)+b1*x(n-1)+b2*x(n-2)+a1*y(n-1)+a2*y(n-2).
x(n)-значение вх. сигнала на текущем такте n.
x(n-1) - то же, но такт раньше.
x(n-2) - то же, но еще такт раньше.
y(n),y(n-1),y(n-2) - аналогично для вых. сигнала.
К-ты.
b0,b2=1.
b1=-2*cos(2*pi*F0/Fd).
Fd - частота дискретизации.
F0 - центральная частота фильтра, 50 или 100 Гц.
a1=2*R*cos(2*pi*F0/Fd).
a2=-R^2.
R - определяет ширину подавляемой области.
Можно найти эксперментально, можно аналитически.
Экспериментально - начните с R=0.8-0.9.
Более узкая полоса - нужно увеличить R.
Но при увеличении учтите, что R должен оставаться меньше 1.
Более широкая полоса - нужно уменьшить R.
Аналитически.
Фильтр имеет передаточную ф-ю:
H(Z)=(b0+b1*Z^(-1)+b2*Z^(-2))/(1-a1*Z^(-1)+a2*Z^(-2)).
Частотная из нее получается подстановкой exp(j*2*pi*F/Fd) вместо Z.
F от -Fd/2 до +Fd/2.
Найдите квадрат модуля частотной характеристики.
Задайтесь необходимыми уровнями и выразите R через требуемую полосу.
Да, отклик фильтра можно смасштабировать. Для этого достаточно умножить к-ты b0-b2 на одну и ту же величину.
Подавление обоих частот - каскадное соединение пары таких фильтров.
Один настроен на 50, второй на 100.
Добавочно каждый каскад, в свою очередь, может состоять из нескольких идентичных каскадов, т.е. фильтров 2-го порядка, настроенных на одну и ту же частоту.
Описанное - один из простых способов.

• Очепетка. — -=ВН=-   (25.11.2004 13:01, 75 байт)