codeprallacqure.asv

本人做的GPS并行码捕获的算法.用的是1024点FFT,原始数据可以假如不同功率的噪声.

%--------------------------------------------------------------------------
%测试载波对并行码相关的影响；
%-------------------------------------------------------------------------
% clear;
% clc;
% load AllSate32CA.mat;
% SateNum1CA=[AllSate32CA(1,:)*2-1 0];
% Shift10ChipNum1CA=[SateNum1CA(1001:1024),SateNum1CA(1:1000)];
% FftSateNum1CA=fft(SateNum1CA);
% FftShift10ChipNum1CA=fft(Shift10ChipNum1CA);
% ConjFftSateNum1CA=conj(FftSateNum1CA);
% PrData=ConjFftSateNum1CA.*FftShift10ChipNum1CA;
% IffData=ifft(PrData);
% plot(IffData);
%==========================================================================
% clear;
% clc;
% load AllSate32CA.mat;
% BiPorAllSate32CA=AllSate32CA*2-1;
% AllSateCAAdd=zeros(1,1023);
% for  i=1:32    
% AllSateCAAdd=AllSateCAAdd+BiPorAllSate32CA(i,:);
% end
% clear i;  
% AllSateCAAdd=[AllSateCAAdd 0];
% SateNum1CA=[AllSate32CA(1,:)*2-1 0];
% Shift10ChipNum1CA=[SateNum1CA(1001:1024),SateNum1CA(1:1000)];
% FftAllSateCAAdd=fft(AllSateCAAdd);
% FftShift10ChipNum1CA=fft(Shift10ChipNum1CA);
% ConjFftSateNum1CA=conj(FftAllSateCAAdd);
% PrData=ConjFftSateNum1CA.*FftShift10ChipNum1CA;
% IffData=ifft(PrData);
% plot(IffData);
%=========================================================================
clear;
clc;
load AllSate32CA.mat;
SateNum1CA=[AllSate32CA(1,:)*2-1 0];
Shift10ChipNum1CA=[SateNum1CA(1001:1024),(SateNum1CA(1:1000))]*(-1);
noiz = wgn(1,1023,2);      %2 DBW

FftSateNum1CA=fft(SateNum1CA);
FftShift10ChipNum1CA=fft(Shift10ChipNum1CA);
ConjFftSateNum1CA=conj(FftSateNum1CA);
PrData=ConjFftSateNum1CA.*FftShift10ChipNum1CA;
IffData=ifft(PrData);
plot(abs(IffData));