Nfreq=64;
Npt=19;
Hd_LP=ideal_response_qs('diam_lp.txt',Nfreq);
h_LP = fwind1(Hd_LP,hamming(Npt));