for?j?=?1:256
????ffts（:j）?=?fftshift（fft（s（:j）））;%獲得一次回波中的不同多普勒頻率
end
for?j?=?1:255
????x1?=?abs（ffts（:j））;
????x2?=?abs（ffts（:j+1））;
????n1?=?1:256;
????n2?=?1:256;
????n12?=?1:256;
????numi?=?0;%進行相關運算以對準
????max?=?sum（sigmult（abs（x1）n1abs（x2）n2））/（sum（abs（x1）.^2）*sum（abs（x2）.^2））^（1/2）;
????for?i=-255:1:255
????????[x12n12]=sigshift（x2n2i）;
????????max1=sum（sigmult（abs（x1）n1abs（x2）n2））/（sum（abs（x1）.^2）*sum（abs（x2）.^2））^（1/2）;
????????if?（?max?????????????max=max1;
????????????numi=i;
????????end
????end
????ffts0=zeros（2561）;
????if（?numi?