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個文件信息