Так, теперь я окончательно запутался и не понимаю что происходит в мире :) 1. Я утверждал, что мнимая часть ЧХ равна нулю, что вроде как неверно. 2. Я построил АЧХ без мнимой части и получась ерунда, что как бы намекает о существовании мнимой части. 3. Я построил действительную часть, мнимую часть и...

так, теперь еще дополнение. так хитро доказано, что нет мнимой части у ЧХ, но на самом деле конечно же есть и она вносит значительный вклад в АЧХ. Если построить сумму действительной и мнимой частей, то получается такая картинка: http://s017.radikal.ru/i435/1111/5d/f8fb3e216820.jpg теперь можно вспо...

удалить

так, надо точно мне понять. ачх я могу построить используя это выражение: http://dspguru.com/sites/dspguru/img/FirFreq.gif сейчас понял что в моём построении был просчет, но он не сильно меняет ситуацию. и еще - реально(насколько реально) ли построить цифровые фильтры с подавлением в 500 и более дец...

ну дело в том, что ачх сравни фильтру с равновеликими коэффициентами. если взять действительную частотную характеристику, то можно наблюдать главный лепесток в 0, pi, 2*pi, но лепесток в точке pi будет перевернут вниз. то есть конечно да, период ачх 2*pi, но если взять модуль ДЧХ, то картинка предст...

Бахурин Сергей писал(а):Я думаю что в этом кроется ваша ошибка. Я уверен, что
\sum_{n=0}^{N-1} {b(n) \cdot sin(n \cdot w)} = 0
Это не верно

В вашем случае это неверно, так как период частотной характеристики у вас , в моем случае период чх равен , поэтому мнимая частотная характеристика равна нулю.

c *лапшой* разобрался. считаю что такой вид из-за ошибок округления и прочего подобного. если увеличить точность, то отчасти этот эффект пропадет, но я еще не проверял, пока не разобрался как увеличить точность.

1. Где расчет самой |H(w)|? 2 что такое f и q? 3 почему q инкрементируется и откуда вы решили взять коэффициент 2 в строке f+=2*sn[q]*np.cos(xa*(n+1-q*2))? Извините, не подумал написать что к чему. 1.f+=2*sn[q]*np.cos(xa*(n+1-q*2)) - это и есть abs(H(w)). Это тоже самое что и H(j\omega)=\sum_{n=0}^...

Бахурин Сергей писал(а):покажите код как считаете эту картинку из импульсной характеристики

Да реально фигня. Ищите ошибку в расчете. Если что можете выложить свою программу можем вместе поискать. Я так понял она на питоне писана (я правда с ним никак не работал но думаю логику можно проследить потому что там ошибка явно системная) Эта лапша возникает при большом коэффициенте подавления в...