|
|
Т.е. y(n)=x(n)+y(n-1). БИХ фильтр с ЧХ=1/(1-exp(-j*2pi*f/Fs)).
И вообще-то Ваш спектрум надо бы умножать 1/(1-exp(-j*2pi*k/N)).
Но спектрум берется от конечной длины, а интегратор - БИХ. Это одна нестыковка. Эта нестыковка приведет к наложениям во временной области после обратного Фурье. Хотя если Вы обратное Фурье делать не будете, то может и пройдет, с учетом следующей нестыковки.
А хрен ее знает, что делать с нулевой частотой. Так она называется. Потому как означенный БИХ имеет аккурат на 0 частоте свой единственный полюс. И бесконечное значение в этом полюсе.
Еще вариант.
Перед Фурье дополните 0 до двойной длины. Возьмите ПФ, умножьте его на ЧХ прямоугольного импульса, длиной = длине сигнала (до дополнения нулями).
Или проинтегрируйте во временной области.
А потом берите Фурье.
|
|
E-mail: 
info@telesys.ru
