%期刊文章中的快速FICP算法
function Rn = FICP(X1k,X2k,N2)  
%1 要完成從兩個信號的譜得到補零之后的2*N點的相關序列
%FICP最基本的思想是對兩信號的互譜補零從而提高相關函數的分辨率
%N2表示補零之后的點數，Rn表示2*N的相關函數輸出序列。
%  1--根據相關定理計算兩信號N點的互譜R1k
R1k=X2k.*conj(X1k);
%plot(abs(R1k));
%  2--通過譜的共扼對稱性，可以獲得R3k
N=length(R1k);
Nn=2:N;
R3k(1)=0;
R3k(Nn)=conj(R1k(N+2-Nn));
%  3--以下是計算n=(0-(N-1))的相關函數值
n=0:N-1;
Pn=exp(-j*2*pi*n*N/N2);
An=conj(MCZT(conj(R1k),N,N2))/N2;
Bn=conj(MCZT(conj(R3k),N,N2))/N2;
Bn=Pn.*Bn;
R1n=An+Bn;
%  以下是計算n=（-N：-1）的自相關函數
Hk=R1k.*Pn;
En=conj(MCZT(conj(Hk),N,N2))/N2;
Mk=R3k.*Pn;
Mn=conj(MCZT(conj(Mk),N,N2))/N2;
Qn=exp(j*2*pi*N*(N-n)/N2);
Fn=Qn.*Mn;
R2n=En+Fn;
Rn=real([R2n,R1n]);