Отправлено GroundCtrl 24 ноября 2004 г. 15:19 В ответ на: Ответ: отправлено GroundCtrl 24 ноября 2004 г. 14:07

частота дискретизации Fs = 500 Hz.
частота запирания F0 = 10Hz.

"цифровая частота" w0 = 2*pi*F0/Fs = 0.6283 (в радианах)
Нули помещаем в точки удиничной окружности, соответствующие частоте 10 Гц. Т.е.выражение числителя

A(z) = (1-exp(j*w0)*z^-1)*(1-exp(-j*w0)*z^-1) = 1-2*cos(w0)*z^-1+z^-2 =
=1-1.618056*z^-1+z^-2.

Знаменатель строится с теми же коэффициентами, которые указаны по ссылке. Т.е. полюсы располагаются в точках p1 = r*exp(j*w0) и p2 = r*exp(-j*w0), где r = 0.937. Тогда, знаменатель:

B(z) = (1-0.937*exp(j*w0)*z^-1)*(1-exp(-j*w0)*z^-1) = 1-1.516066*z^-1+0,877969*z^-2

• У вас получились коэффициенты совпадающие а у меня нет. Но на графике рисует правильно. Наверное будет работать. Осталось получить от товарища сигнал в отсчетах и прогнать в MATLABе. — misyachniy   (24.11.2004 15:56, пустое)