[f1,f2] = freqspace(9,'meshgrid');
Hd = ones(9); 
r = sqrt(f1.^2 + f2.^2);
Hd(r>0.3) = 0;
colormap(jet(64))
mesh(f1,f2,Hd)
win = fspecial('gaussian',9,2);
win = win ./ max(win(:)); 
%mesh(win)
h = fwind2(Hd,win);
freqz2(h)