function xp=fftxzhijie(duan,cy)
%得到各頻點的復信號
%得到EH-4的X-FILE，duan為Y文件中得到的4096*4信號，T為采樣時間，XP為得到的X文件。
for i=1:size(duan,2)
    x=duan(:,i);
    [yh(:,i),yw(i,:)]=fftt(x,cy);%得到每一道各頻點的HIBERT 譜；
    [h(:,i),w(i,:)]=findxchazhi(yh(:,i),yw(i,:));%得到每一道與S所對應頻點的HIBERT 譜；
end
[m,n]=size(duan);

%求得各頻率點的自HIBERT 自譜或互譜后，再對M點平均
ZHY=(h(:,1).*conj(h(:,1)))/m;
%ZEX=(h(:,2).*conj(h(:,2)))/m;
ZHX=(h(:,3).*conj(h(:,3)))/m;
%ZEY=(h(:,4).*conj(h(:,4)))/m;

HEXHX=(h(:,2).*conj(h(:,3)))/m;
HEXHY=(h(:,2).*conj(h(:,1)))/m;
HEYHX=(h(:,4).*conj(h(:,3)))/m;
HHYHX=(h(:,1).*conj(h(:,3)))/m;
HHXHY=(h(:,3).*conj(h(:,1)))/m;
HEYHY=(h(:,4).*conj(h(:,1)))/m;

xp=zeros(length(w),9);
xp(:,1)=w(1,:)';
xp(:,2)=HEXHX;
xp(:,3)=HEXHY;
xp(:,4)=HEYHX;
xp(:,5)=HEYHY;
xp(:,6)=ZHX;
xp(:,7)=HHYHX;
xp(:,8)=HHXHY;
xp(:,9)=ZHY;
end