逆合成孔径雷达平动补偿

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 < max1）
            max=max1;
            numi=i;
        end
    end
    ffts0=zeros（2561）;
    if（ numi