dlay.m

本次的算法是相位周期展开聚类

function dsum=dlay(Ap,d,w);
c=340;
%s=w;
n=d*w/c;
k=ceil(n-0.5);
k1=ceil((d*8000/c)-0.5);
a=[];
m=3.1416;
%for k=1:n-1
        %if(w>M(k)&w<=M(k+1))
             for l=0:k
                    ds=(Ap+2*l*m)/w;
                    %if k~=k1
                    a=[a,ds];
                   
             end
             if k~=k1
                 dsum=[a,zeros(1,(k1-k))];
             else
                 dsum=[a,zeros(1,(k1-k))];
                 dsum=a;
             end
               %for k:k1
               %a=[a,0];
               %end
        %end
% end
%dsum=a;