delta_x = 0.6;
Tp = 0.4e-6;
Wr = 575.3;
Wa = 442.9;
Fc = 242.4e6;
alpha = 33.375e6;
BW = 133.5e6;
Fs = 2560e6;
Dr = 1;
Da = 1;
%Kr = 0.89;
%Ka = 0.89;
R_scene = 1000;
Phi_ac = pi / 2;
Theta_ac = pi / 2;
L = 760.8;
dTheta = 0.231 * pi;
N_points = 3;
Ns = 3;
point_coord = [0 0 0; 200 0 0; 0 -100 0];

Pc = [1 1 1];
k = -L/2:delta_x:L/2-delta_x;
t = 0:1/Fs:Tp-1/Fs;
f = linspace(Fc-BW,Fc,length(t));
Na = length(k);
Nf = length(t);
V = delta_x/Tp;
c = 3e8;
lambda = c/Fc;
Vsum = zeros(length(k),length(t));
for i = 1:Na
    dist_ant = k(i);
    for n = 1:N_points
        yc = R_scene + point_coord(n,2);
        rc = sqrt((point_coord(n,1) - dist_ant)^2 + yc^2);
        Vc = sqrt(Pc(n))/(rc^2)*exp(-1j*4*pi*rc./lambda).*exp(1j*BW/Tp.*(t-2*rc/c).^2);
        Vsum(i,:) = Vsum(i,:) + Vc;
    end
end

Kr = 4*pi*f/c;
Kr = repmat(Kr,Na,1);
Vsum = Vsum.*exp(1j*R_scene*Kr);
%компенсация фазовых набегов %перемещение отражателей по каналам дальности в течение полета – ЛЧМ дальномер
Vsum_compr = fftshift(fft(Vsum,[],2),2);
imagesc(abs(Vsum_compr));

%Vsum = fftshift(Vsum,1);
Vsum = fft(Vsum,[],1);
Vsum = fftshift(Vsum,1);
%результат азмиутального сжатия сигнала
Vsum_compr = fftshift(ifft(Vsum,[],2),2);
imagesc(abs(Vsum_compr));

Kx = linspace(-pi/delta_x,pi/delta_x,L/delta_x).';
Kx = repmat(Kx,1,Nf);
Ky_sq = Kr.^2-Kx.^2;
Neg = logical(Ky_sq < 0);
Ky_sq(Neg) = NaN;
Vsum(Neg) = NaN;
Kr(Neg) = NaN;
Fmf = -R_scene*Kr+R_scene*sqrt(Ky_sq);
Vsum = Vsum.*exp(1j*Fmf);
Vsum_compr = fftshift(ifft(Vsum,[],2),2);
imagesc(abs(Vsum_compr));