- h = sqrt(j * c * Tsamp^2 / (D * lambda^2 * z)) * ...
- exp(-j * pi * c * Tsamp^2 * k .^ 2 / (D * lambda^2 * z));
+ h = exp(-j * pi * c * Tsamp^2 * k .^ 2 / (D * lambda^2 * z));
+
+ len_fft = max(length(y), length(h));
+ H = fft(h, len_fft);
+ Y = fft(y, len_fft);
+
+ yCDComp = ifft(H.' .* Y);
+ l = (length(h) - 1) / 2;
+ if l > 0
+ yCDComp = [yCDComp(l:end); yCDComp(1:l-1)];
+ else
+ yCDComp = [yCDComp(end); yCDComp(1:end-1)];
+ end