framesamples=floor(Ts./dt);
pdsamples=floor(Tm./dt); % 在脈沖寬度上的采樣點個數，以脈沖的寬度進行滑動相關
bitsamples=framesamples*Ns;%Ts./dt
N=Tf./Tm;%在一bit信息內有N個不同狀態

zp=zeros(1,N);
for n=1:N;
    mx=RX1(1+(n-1)*pdsamples:bitsamples+(n-1)*pdsamples); %RX1為接受信號
    mxk=mx.*mask;  %mask 是本地模板信號
    zp(n)=sum(dt.*mxk); %dt＝1/fc,
end


if bits(1)
    else zp=-zp;
end           %  如果增加駐留時間和發送的bit數時， 此處需要修改

sd0=No/2;
zp1=abs(zp)./sd0(1);

k=30;

zp2=zeros(1,ceil(N/k)*k+k);
zp2(1:max(size(zp1)))=zp1;
y=zeros(1,ceil(N/k));
for i=1:ceil(N/k)
        y(i)=sum(zp2((i-1)*k+1:(i-1)*k+k));
end
 [Ymax,state]=max(y);
 delay2=state*k*Tm;
 state1=(state-1)*k;
 
 k1=60;
 selvalue=Ymax./k;%門限一
pha=find(zp1(state1:state1+k1)>selvalue);% 數組 
 
    phase=pha+state1-1;   
   z=zp1(phase); 
 q=max(z);

   selvalue1=3/2*selvalue;%門限2
   
[z1,seq]=sort(z); %把z按照從大到小的順序排列

   Lnum=max(size(z1));
    index=zeros(1:Lnum);
   phase1=zeros(1:Lum);
   for p=1:Lnum
       if z1(Lnum-p+1)> selvalue1
          index(p)=seq(Lnum-p+1);p
          phase1(p)=(index(p)+state1-1).*Tm-delay; %phase1為最后輸出的相位，沒有規定要多少個 ，過門限的都留下了 
       end
   end