function y=farrow(x,mu)
    %插值器采用分段抛物线插值滤波器的Fallow结构(四点三阶的分段抛物线，输入x含4个点，索引为-2 -1 0 1，当前定时点mk对应0点。)
    %x：定时恢复输入数据  mu来自NCO的定时误差
    %y: 插值器输出
    alpha=0.5;
    c0=[0 0 1 0];
    c1=[-alpha 1+alpha alpha-1 -alpha];
    c2=[alpha -alpha -alpha alpha];
    c_array=[c0(1) c0(2) c0(3) c0(4);
            c1(1) c1(2) c1(3) c1(4);
            c2(1) c2(2) c2(3) c2(4)];
    m=length(c_array(1,:));
    v0=0;
    v1=0;
    v2=0;
    
    for i=1:1:m
        v0=v0+x(m+1-i)*c_array(1,i);
        v1=v1+x(m+1-i)*c_array(2,i);
        v2=v2+x(m+1-i)*c_array(3,i);        
    end
    y=(v2*mu+v1)*mu+v0;
end