x = complex(ResultI, ResultQ);

D = 64;

 N=512; %FFT size
m1=reshape(x,length(x)/D,D); %generating matrix ---> each row showing range info for one period AND each column showing number of periods
[My,Ny]=size(m1');
win=hamming(Ny);
m2=(win*ones(1,My));applying window for sidelobe reduction (in time domain)
M2=(fft(m2,2*N)); %First FFT for range information
M3=fftshift(fft(M2',2*D)); %Second FFT for doppler information
[My,Ny]=size(M3);