function outdata=H_interp(indata)
len=length(indata);tmp=zeros(10,1);
if 1==1, %linear interpolation
    %for i=1:5,
     %   outdata(i)=(6-i)*(indata(2)-indata(1))/14+indata(1);
     %end
    for i=1:13
        outdata(5+i)=(i)*(indata(2)-indata(1))/14+indata(1);
    end
    for i=1:13
        outdata(18+i)=i*(indata(3)-indata(2))/14+indata(2);
    end
    for i=1:13
        outdata(31+i)=i*(indata(4)-indata(3))/14+indata(3);
    end
    %for i=1:5
     %   outdata(44+i)=i*(indata(4)-indata(3))/14+indata(4);
     %end
    for i=1:10
        tmp(i)=i*(indata(1)-indata(4))/11+indata(4);
    end
    outdata(45:49)=tmp(1:5);
    outdata(1:5)=tmp(6:10);
end
tmp=zeros(48,1);
tmp(1:24)=outdata(1:24);
tmp(25:48)=outdata(26:49);
outdata=conj(tmp);