|
|
Подсказываю.
Заграждающий фильтр можно следующим образом построить.
1.Разместить 0 на 50 Гц и на -50 Гц и на единичной окружности в Z плоскости. Т.е. комплексно-сопряженный 0 на 50 Гц и единичной окружности.
Это даст числитель перед. ф-ии звена 2-го порядка, при 500 герцовом квантовании, равным:
(1-exp(-j*2pi*0.1)*Z^(-1))*(1-exp(j*2pi*0.1)*Z^(-1))=
=1-2*cos(2pi*0.1)*Z^(-1)+Z^(-2)=1-1.618*Z^(-1)+Z^(-2).
1,-1.618,1 - Ваши к-ты. А 0.1=50/500.
2. Поскольку загражд. фильтр из одного нуля, даже комплексно-сопряженного довольно хреновый, то надо бы и полюса к нему присобачить. С целью компенсации хреновости 0.
Из этого следует, что должна быть комплексно-сопряженная пара полюсов. Один из полюсов должен располагаться на 50 Гц, но не на единичной окружности, а на окружности меньшего радиуса, R.
Второй на -50 герцах и той же окружности.
Т.е. знаменатель перед. ф-ии должен быть равен:
(1-R*exp(-j*2pi*0.1)*Z^(-1))*(1-R*exp(j*2pi*0.1)*Z^(-1))=
=1-2*R*cos(2pi*0.1)*Z^(-1)+R^2*Z^(-2).
Тут всего один неизвестный параметр. R. Он будет полосу определять.
Найти его можно, подставив вместо Z exp(j*2pi*f/Fd).
И взяв модуль (или его квадрат) перед. ф-ии. И выразив R через требуемую полосу. Ну т.к. Вам требуется загр. фильтр из 2-х одинаковых звеньев 2-го порядка, то соответственно нужно брать квадрат модуля от произведения 2-х одинаковых перед. ф-ий.
Делать этого не собираюсь, но замечу, что 1,5164 это как раз 2*R*cos(2pi*0.1). А 0.8783 как раз R^2. Откуда R=0.93718.
А какая ф-ия в матлабе - не знаю.
Можно и по другому соорудить загр. фильтр. Но и так уже много написал.
E-mail: 
info@telesys.ru
