test

pss0 = s2((N_f-72-1023):(N_f-72));
p_f = fft(pss0);
p_f2 = [0 p_f(1:36) p(1024+(-35:0))];
H = p_f2 .* conj([0 pss0_f(37:72) pss0_f(1:36)]);