%desired filter length
N = 51;
f = [linspace(-1,-.18,164),linspace(-.1,.3,80),linspace(.38,1,124)];
%frequency grid
d = [zeros(1,164),ones(1,80),zeros(1,124)].*exp(-1i*pi*f*(N-1)/2);
%desired frequency response
w = [10*ones(1,164),ones(1,80),10*ones(1,124)];
A = w(:,ones(1,N)) .* exp(-1i*pi*f*(0:N-1) );
h = A  (w.*d);