Синтез режекторного КИХ - фильтра

[parenya]

Привет участникам форума!
Мне нужно синтезировать режекторный фильтр с параметрами(вообще задание синтезировать смесь реж-го и ФНЧ с линейной ФЧХ, но не получается именно этот этап): f1=18,f2=21,f3=24,f4=42, fmax=132,все в кГц, a0=3,amin=-35.
Почему-то в полосе задердживания АЧХ никак не опускается ниже 14 дБ. Буду рад, если кто-то сможет решить эту задачу, код:

Используется синтаксис Matlab M

f1=18e3;
f2=21e3;
f3=24e3;
f4=42e3;
f7=135e3;
fd=2*f7;
amin=35;

n=-200:200;
N=171;


CE=find(n==0);


fd=2*f7;
fc1=(f1+f2)/(2*fd)
fc2=(f3+f4)/(2*fd)

h1(CE)=1-2*(fc1-fc2);
for k=1:length(n)

    if k~=CE
      h1(k)=(sin(2*pi*n(k)*fc1)-sin(2*pi*n(k)*fc2))/(n(k)*pi);
    end
end


figure(1)
stem(n,h1);


alf=0.5842*(amin-21)^0.4+0.07886*(amin-21)

w=kaiser(N,alf);

hr01=h1(CE-fix(N/2):CE+fix(N/2));

figure(2)
stem(n(1:length(hr01)),hr01);

hr1=hr01.*transpose(w);

a=1;

[H1,f]=freqz(hr01,a,1000,fd);

figure(3);
plot(f,20*log10(H1))

[Toucan]

parenya, для оформления программного кода следует использовать тег syntax. Поправил.