clear all;
Fs = 400;   %частота дискретизации (Гц)
f0 = 100;    %центральная частота фильтра (Гц)
band = 50;    %полоса фильтра (Гц)

N = 256;    %количество к-тов фильтра (число линий задержки равно N-1)

w = hamming(N)'; %окно хемминга

%моменты дискретизации для линейной ФЧХ
t = (-N/2+0.5):(N/2-0.5);
t = t*band/Fs;

h = w.*sin(pi*t)./(pi*t); %расчет ФНЧ c полосой band

h = 2*h/sum(h);              %нормирую

h=h.*cos(2*pi*t*f0/band); %переношу на частоту f0
%в массиве h лежат к-ты полосового фильтра, рассчитанного на частоту f0 с полосой band
% НЕ меняя структуры фильтра (не меняя N) можно расчитать фильтр на любую частоту с любой полосой
%(разумеется чтобы он укладывался от 0 до Fs/2)

% построим АЧХ фильтра при помощи fft

f = (0:8*N-1)*Fs/(8*N);

H = 20*log10(abs(fft([h, zeros(1,7*N)])));

subplot(211), stem(1:N,h), grid;
subplot(212), plot(f,H), grid, axis([0,Fs/2,-140,10]);