NumBands = 10;
L = 256
% h  - ФНЧ прототип. Каждая из полос получается сдвижкой по частоте АЧХ
% ФНЧ прототипа
h = remez(L-1, [0,1/(2*NumBands), 1.5/(2*NumBands), 1],[1, 1, 0, 0]);

h_all = zeros(1, length(h));
A = ones(1, NumBands); % коэф. усиления каждой из полос эквалайзера

for i =2:5
  A(i) = A(i-1)*power(10, -10/20); % -10dB step down
end

for i =6:10
  A(i) = A(i-1)*power(10, 10/20); % +10dB step up  
end    


for i=1:NumBands-1
    h_all = h_all + A(i)*sin( i*pi/2+(1:L)*i*pi/NumBands);
end

h_all = h_all.*h;
%h_all - синтезированный фильтр
figure(3);
freqz(h_all,1)