[an error occurred while processing this directive]
|
Если у Вас в выражении для H(Z) нет описки в степенях z, т.е. они на самом деле положительные, то нужно
вынести Z^2 за скобки и в числителе и в знаменателе.
Получите в более привычном виде:
H(Z)=(z^(-2)-1.45514z^(-1)+0.57831)/(1-1.45514z^(-1)+0.57831*z^(-2))
Это вполне устойчивый фильтр.
Если для fdatool требуются к-ты числителя и знаменателя H(Z) (я не знаю, что ему требуется), то они:
numerator, то бишь числитель, b=[0.57831,-1.45514,1].
denominator, знаменатель, a=[1,-1.45514,0.57831].
А в разностное уранение 2 последних к-та знаменателя, т.е. рекурсивной части, подставляются с обратным знаком.
Т.е. разностное уравнение Вашего фильтра:
y(n)=0.57831*x[n]-1.45514*x[n-1]+x[n-2]+1.45514*y[n-1]-0.57831*y[n-2].
Повторю - это устойчивый фильтр.
Если же в приведенном Вами выражении степени z на самом деле меньше 0, то фильтр будет неустойчив.
Т.е. если H(Z) на самом деле равна:
H(Z)=(1-1.45514z^(-1)+0.57831z^(-2))/(z^(-2)-1.45514z^(-1)+0.57831)
Такая H(Z) соответствует неустойчивому фильтру.
По той простой причине, что к-т при z^-2 в знаменателе, после всех приведений, больше 1.
А т.к. к-т при z^-2 в знаменателе для БИХ 2-го порядка есть квадрат модуля полюса, то его равенство или превышение 1 означают, что полюс или на единичной окружности или снаружи ее. Что соответствует неустойчивому фильтру.
E-mail: info@telesys.ru