Отправлено ВН 20 января 2003 г. 17:47 В ответ на: Уточнение: А может есть для этих целей специальные фильтры или методы расчета. отправлено valera_by 20 января 2003 г. 16:06

Да просто считаются фильтры, хоть bandstop, хоть bandpass, хоть ...
Но, раз уж упомянули про резонатор, то может подойдет уж совсем простое решение. Аналог колебательного контура.
1. bandstop получается из bandpass:
Hbs(Z)=1-Hbp(Z). Естественно при условии, что Hbp нормирована к 1.
В противном случае вместо 1 в выражении будет другое число, равное максимуму модуля Hbp.
2.Резонатор, bandpass:
Hbp(Z)=A/(1-B*Z^(-1)).
Коэффициент B здесь комплексный:
B=b*exp(j*2*pi*W0). W0 - цифровая резонансная частота=F0/Fd, b - действительный коэффициент, определяет добротность резонатора. Меньше 1, при 1 резонатор неустойчив. Чем ближе к 1, тем выше добротность.
Коэффициент A действительный и выбирается из условий нормировки.
Чтобы max(|Hbp|) был равен 1, A=1-b.
Тогда Hbs(Z)=(1-A-B*Z^(-1))/(1-B*Z^(-1))=(b-B*Z^(-1))/(1-B*Z^(-1)).
Это комплексный фильтр. На выходе комлексный сигнал. Действительный из него получается игнорированием на выходе мнимой части (но не при вычислениях).
Или использовать другие выражения:
Hbp(Z)=a/(1-b1*Z^(-1)-b2*Z^(-2)). a,b1,b2 - действительные.
b1=2*b*cos(2*pi*W0), b2=-(b^2), a=sqrt(((1+b2)^2)*(1+(b1^2)/(4*b2)))
a - из нормировки к 1.
И Hbs=1-Hbp. Это уже действительный фильтр. На выходе действительный сигнал (при действительном входном). Требования к b аналогичны первому случаю. Т.е. и в том и другом случае вся характеристика определяется двумя величинами - рез. частотой (W0) и к-том b.