clear; close all; wc = pi/8; T = 1; N=61; l = (N-1)/2; ks = 1:N; h1 = []; h2 = []; h3 = []; h4 = []; w0 = pi/2; for kprime = ks k = kprime-1; hk1 = wc*T/pi*sin((k-l)*wc*T)/(k-l)/wc/T; h1(end+1) = hk1; hk2 = -wc*T/pi*sin((k-l)*wc*T)/(k-l)/w…