GPS_CA.rar

GPS_CA.rar

close all
clear r r1 rb1 rb2 rb3;
choice = 1;

% 捕获程序
% 捕获程序中暂没有考虑多普勒影响
signal = wave;
deltaf = （-1:0.1:1）*1000;
freq_search = deltaf+100*1e6;

Rs = 1023000*96;
deltat = 1/Rs;
% wave = waveI+j*waveQ;

freq_num = length（freq_search）;
bdm = reshape（repmat（SpRs*10/1000/102301）1[]）.*boc11;
bdm1 =reshape（repmat（SpRs*10/1000/102301）1[]）.*boc61; 
bdm2 =reshape（repmat（SdRs*10/1000/102301）1[]）.*boc11; 

bdm3 = 0.5*bdm1+sqrt（1/11）*bdm2;
%% 
%****************
%相关性分析
Cs = reshape（repmat（SdRs*10/1000/102301）1[]）.*boc11;
Cp1 = reshape（repmat（SpRs*10/1000/102301）1[]）.*boc11;
Cp2 = reshape（repmat（SpRs*10/1000/102301）1[]）.*boc61;

Cs = Cs（1:1023*96）;
Cp1 = Cp1（1:1023*96）;
Cp2 = Cp2（1:1023*96）;

t = deltat*[0:length（Cs）-1];
for n = 1:length（deltaf）
bcI = 0.5*Cs+sqrt（1/11）*Cp2;
bcQ = sqrt（29/44）*Cp1;
% rbcI = bcI.*sin（2*pi*deltaf*t）;
rbcQ = -bcI.*sin（2*pi*deltaf（n）*t）+bcQ.*cos（2*pi*deltaf（n）*t）;
rbcI = bcI.*cos（2*pi*deltaf（n）*t）+bcQ.*sin（2*pi*deltaf（n）*t）;
r（n） = sum（rbcQ.*bdm（1:1023*96））;
rb1（n） = sum（rbcI.*bdm1（1:1023*96））;
rb2（n） = sum（rbcI.*bdm2（1:1023*96））;
rb3（n） = sum（rbcI.*bdm3（1:1023*96））;
end
% disp（r）;



fcs = fft（Cs）;
fcp1 = fft（Cp1）;
fcp2 = fft（Cp2）;

rsp1 = ifft（conj（fcs）.*fcp1）;
rsp2 = ifft（conj（fcs）.*fcp2）;
rp12 = ifft（conj（fcp1）.*fcp2）;
figure
subplot（311）;
plot（abs（rsp1）‘r-‘）;
% hold on
subplot（312）
plot（abs（rsp2）‘b-‘）;
subplot（313）
plot（abs（rp12）‘black-‘）;

rs = ifft（conj（fcs）.*fcs）;
rp1 = ifft（conj（fcp1）.*fcp1）;
rp2 = ifft（conj（fcp2）.*fcp2）;
figure
subplot（311）
plot（abs（rs）‘r-‘）;
% hold on
subplot（312）
plot（abs（rp1）‘b-‘）;
subplot（313）
plot（abs（rp2）‘black-‘）;
%***************




if choice ==1 
    for k = 1:freq_num
        rev_sig = imag（wave.*exp（-1j*2*pi*freq_search（k）*deltat*[0:length（signal）-1]））;  %其实这里可能还有相位的影响。
        rev_sig1 = rev_sig（1023*96*5+1:1023*96*6）;
%         sum（rev_sig1.*bdm（1023*96*7+1:1023*96*8））
        rev_sig2 =[rev_sig1zeros（1length（bdm）-length（rev_sig1））]; 
        
        % n = m.*exp（-1j*2*pi*fc1*[0:1/Rs:T-1/Rs]）;
        % figure
        % plot（rev_sig）;
        % end
%         bdm1 = [bdmzeros（1length（rev_sig）-length（bdm））];
        a = fft（rev_sig2）;
        b = fft（bdm）;
        m = ifft（conj（a）.*b）;
        figure
        plot（abs（m））;
        r1（k） = max（abs（m））
        sum（rev_sig1.*bdm（1023*96*5+1:1023*96*6））
        title（‘时间 1ms‘）;
        % rev_sig = wave.*sin（2*pi*freq_search（k）*deltat*[1:length（signal）]）;
    end
    figure
    stem（r‘*‘）
    figure
    stem（r1‘*‘）
elseif choice == 2
    for k = 1:freq_num
        rev_sig = imag（wave.*exp（-1j*2*pi*freq_search（k）*deltat*[0:length（signal）-1]））;  %其实这里可能还有相位的影响。
        % rev_sig = rev_sig（1:1023）;
        % n = m.*exp（-1j*2*pi*fc1*[0:1/Rs:T-1/Rs]）;
        % figure
        % plot（rev_sig）;
        % end
        bdm1 = [bdmzeros（1length（rev_sig）-length（bdm））];
        a = fft（rev_sig）;
        b = fft（b

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        247  2017-11-28 17:13  Boc61.m

     文件        241  2017-12-06 18:22  Boc11.m

     文件        238  2017-12-06 18:30  invert.m

     文件        215  2017-12-06 18:30  E.mat

     文件        215  2017-12-06 18:31  N.mat

     文件        206  2017-12-07 19:37  ca_code_sample.m

     文件        203  2017-12-06 18:30  B1I_code_sample.m

     文件     110474  2017-12-06 18:30  GPS_receiver_interface_L2C_2.fig

     文件      48091  2017-12-06 18:30  GPS_receiver_interface_L2C_2.m

     文件      24176  2017-12-06 18:30  code_B3I.mat

     文件      23996  2017-12-06 18:23  data.mat

     文件      10520  2017-12-06 18:30  code_B1I.mat

     文件       8387  2017-12-06 18:31  Mypolar.m

     文件       7451  2017-12-07 19:37  GPS_CA_tracking_2.m

     文件       6935  2017-12-06 18:30  B1I_tracking.m

     文件       6319  2017-12-28 19:24  GPS_CA_receiver.m

     文件       6269  2017-12-06 18:30  B1I_receiver.m

     文件       4825  2017-12-06 18:30  GPS_CA_acquisition.m

     文件       4746  2017-12-06 18:30  B1I_acquisition.m

     文件       4614  2017-11-28 17:30  acquisition.m

     文件       4610  2017-12-06 18:31  satpos.m

     文件       4598  2017-12-06 18:30  ephemeris.m

     文件       4231  2018-01-13 17:38  leastSquarePos.m

     文件       3411  2017-12-06 18:31  PVT.m

     文件       2900  2017-12-06 18:31  tropo.m

     文件       2892  2017-12-06 18:31  togeod.m

     文件       2529  2017-12-06 18:30  B1I_GetTOW.m

     文件       2529  2017-12-06 18:30  GPS_GetTOW.m

     文件       2508  2017-12-06 19:00  B1C_code.m

     文件       2473  2017-11-28 17:13  BOC.m

............此处省略33个文件信息