I think u can use the Multiplication properity of the discreate time fourier transform
which is " Multiplication on time domain is a convolution in frequency domain"
so if u have function
h=xy ---->DTFT----> H(e^jw)=1/(2pi)int(X(e^jw)Y(e^j(w-m))dm,0,2pi)