GPS Acquisition IF frequency in graph

I am trying to generate GPS Acquisition signals using frequency code phase search acquisition.

1. I generated C/A gold codes and digitised at 16.368e6 MHz(sampling frequency)
2. I generated carrier signals using exponential function at sampling frequency
3. I have input data with IF of 4.092e6 and sampling frequency as mentioned above. my inout data has 163680 samples
4. I multiplied input data with carrier signals which produces 1023 samples by considering 1023 samples from input data and converted into Fourier transform using FFT
5. I converted gold codes into frequency domain using FFT
6. I multiplied gold codes with step 4 to get final output with 1023 samples.

I have converted to absolute values and ploted. Then I got graphs as below. How to plot so that I can get IF in middle of graphs and how to plot correctly?

final_out is output from step 6:

[amp crw] = max(max(abs(final_out')));

[amp ccn] = max(max(abs(final_out)));

figure;

plot(abs(final_out(ccn,1:1023)));

[Mehmed Saad]

1. Use fftshift to discriminate negative and positive frequencies
2. Use db for better visualization
3. Create x-Axis from -fs/2 to +fs/2 (1024 samples)

figure,plot(linspace(-fs/2,+fs/2,1024),db(abs(fftshift(fft(your_data,1024)))))

[Rohit Mallampati]

Sir, Here fs means my sampling frequency (16.368e6 Hz) or chip rate (1.023e6 Hz)?

[Mehmed Saad]

sampling rate