y = conv(x1,x2); % time domain convolution of x1 and x2 h1 = freqz(x1, 1, w); % DTFT of sequence x1 h2 = freqz(x2, 1, w); % DTFT of sequence x2 ...