uangle1 = unwrap(angle(port(1+port_pilot)));
uangle2 = unwrap(angle(port(1+end-port_pilot)));
for k = 1:num_port_pilot_1side-1
    pha1 = uangle1(k:k+1);
    pha2 = uangle2(k:k+1);
    idx = port_pilot(k:k+1);
    istart = idx(1)+1;
    iend = idx(2)-1;
    phases1 = interp1(idx,pha1,istart:iend)';
    port(1+(istart:iend)) = port(1+(istart:iend)).*exp(-1i*phases1);
    phases2 = interp1(idx,pha2,istart:iend)';
    port(1+end-(istart:iend)) = port(1+end-(istart:iend)).*exp(-1i*phases2);
end