function [out]=estimation_rotate_header(datai,dataq,iv,qv);
out=zeros(1,128*12);
tempi(1:128*6)=(iv.*datai(1:128*6)+qv.*dataq(1:128*6));
tempq(1:128*6)=(-qv.*datai(1:128*6)+iv.*dataq(1:128*6));
tempi(128*6+1:128*12)=(iv.*datai(128*6+1:128*12)+qv.*dataq(128*6+1:128*12));
tempq(128*6+1:128*12)=(-qv.*datai(128*6+1:128*12)+iv.*dataq(128*6+1:128*12));

temp=tempi+j*tempq;

for k=1:128*6
    if((iv(k)~=0)|(qv(k)~=0))
        out(k)=temp(k)/sqrt((iv(k)^2+qv(k)^2));
        out(k+128*6)=temp(k+128*6)/sqrt((iv(k)^2+qv(k)^2));
    end
end