function g = sfactK(h,K,Z)
%
% sfact with K zeros at z==Z
%

v = 1;
for k = 1:K
   v = conv(v,[1 -Z]);
end

g = sfact(deconv(h,xcorr(v)));
g = conv(g,v);
g = g*sqrt(max(h)/sum(g.*g));