[ht,h]=abcd2hcx(a,b,c,d);

function [ht,h]=abcd2hcx(a,b,c,d);
d=a;
r=75;
sz=length(a);
h=2*ones(sz,1)./(c*r+d+a+b/r);
h(sz+1:2*sz-1)=real(h(sz:-1:2))-sqrt(-1)*imag(h(sz:-1:2));
ht=real(ifft(h));