|
|
для фильтра из одних полюсов АЧХ можно посчитать как 1/abs(фурье(a)) где а - вектор коэффициентов знаменателя, добитый нулями до необходимой длины фурье. А деление - поэлементное, то есть для каждой точки к делается 1/abs(X(k)). Так тогда в заданной точке можно быстро посчитать коэффициент передачи непосредственно через DFT, зная что коэффициентов "a" всегда 21 штука (первый всегда равен 1), и от всей суммы останется всего ничего слагаемых (остальные нули, так как вектор "а" ими был дополнен).
Просто я что-то никак не сооброжу - матлаб считает АЧХ/ФЧХ вот так: h = fft(b,2*nfft)./fft(a,2*nfft));
где а и b векторы коэффициентов полиномов передаточной функции знаменателя и числителя соответственно. Это разве то-же самое, что и фурье от ИХ?
E-mail: 
info@telesys.ru
