%?script?file:?acquire_compare.m

%?purpose:?產生1ms的GPS數據，設干擾來向為?11.2?度，而權值卻是干擾來向為?13.7?度時求得的。
%??????????比較未加寬零陷和加寬零陷的捕獲效果。（權值在widenull.m里求得）
%----------by?lichan??-----?08年11月29日

%*****************************************************************
N?=?7?;???%-----陣元數目
D?=?[0:N-1]‘;
%--------產生GPS信號---------
t1?=?518417;
t2?=?t1?+?0.001;???????????????????????????????%-----1ms?GPS數據（5000個點）
[svtstsvt1]?=?bbsignal（t1t2）;
fs?=?5e6?;??%-----采樣率
M?=?round（fs?*?（t2-t1））?;???%-----?M=5000
thelta_s?=?50?;??%-----信號來向
a_s?=?exp（-j*pi*D*sin（thelta_s*pi/180））?;???%-----信號導向矢量
s_d?=?a_s?*?svt（21:M）?;???%-----期望信號

%-----?加干擾?------
thelta_j?=?12.2?;???%------干擾來向
a_j?=?exp（-j*pi*D*sin（thelta_j*pi/180））;???%-----干擾的導向矢量
[yy1y2y3pow]?=?jammer;
jam?=?a_j?*?y（11:M）?;

%-----?加噪聲?----
snr?=?-20?;
P_n?=?10.^（-snr/10）?;??%-----噪聲功率
noise?=?zeros（NM）?;
for?i?=?1:N
????noise（i:）?=?sqrt（P_n/2）?*?（randn（1M）?+?j*randn（1M））;
end

s_r?=?s_d?+?jam?+?noise?;

%****************************************************************
%---------?沒有加寬零陷的捕獲?-------------
datalength?=?M?;
CAcode?=?svt1（2:）?;
y?=?w_opt‘?*?s_r?;
figure?;
set（gcf‘color‘‘white‘）?;
fftcode1?=?fft（y（11:datalength））?;
fftcode2?=?fft（CAcode（11:datalength））?;
fftr?=?fftcode1?.*?conj（fftcode2）?;
r?=?fftshift（abs（ifft（fftr）））?;
r?=?r/max（r）?;
t=[-datalength/2:1:datalength/2-1];
plot（tr‘k-‘）;
%?axis（[-50050001]）;
%?xlabel（‘code?phase‘）;
%?ylabel（‘normalized?coherent?integration?‘）;
%?title（‘conventional?method‘）;
xlabel（‘碼位‘）;ylabel（‘歸一化互相關值‘）;
title（‘傳統最小功率法‘）;

%--------?加寬零陷后的捕獲?------------
y1?=?w_wave‘?*?s_r?;
figure?;
set（gcf‘color‘‘white‘）?;
fftcode1?=?fft（y1（11:datalength））?;
fftcode2?=?fft（CAcode（11:datalength））?;
fftr?=?fftcode1?.*?conj（fftcode2）?;
r?=?fftshift（abs（ifft（fftr）））?;
r?=?r/max（r）?;
t=[-datalength/2:1:datalength/2-1];
plot（tr‘k-‘）;
%?axis（[-50050001]）;
%?xlabel（‘code?phase‘）;
%?ylabel（‘normalized?coherent?integration‘）;
%?title（‘new?method‘）;
xlabel（‘碼位‘）;ylabel（‘歸一化互相關值‘）;
title（‘微分約束最小功率法‘）;


%?datalength?=?M?;
%?CAcode?=?svt1（2:）?;
%?y1?=?w_opt‘?*?s_r?;
%?y2?=?w_wave‘?*?s_r?;
%?
%?fftcode?=?fft（CAcode（11:datalength））;
%?fftcode1?=?fft（y1（11:datalength））;
%?fftcode2?=?fft（y2（11:datalength））;
%?fftr1?=?fftcode.*conj（fftcode1）;
%?fftr2?=?fftcode.*conj（fftcode2）;
%?r1?=?fftshift（abs（ifft（fftr1）））;
%?r1?=?r1/max（r1）;
%?r2?=?fftshift（abs（ifft（fftr2）））;
%?r2?=?r2/max（r2）;
%?figure;
%?set（gcf‘color‘‘white‘）;
%?
%?
%?t?=?-datalength/2?:?1?:?datalength/2-1;
%?subplot（211）?plot（tr1‘k-‘）;
%?ylabel（‘normalized?coherent?integration‘）
%?title（‘conventional?method‘）
%?subplot（212）?plot（tr2‘k-‘）;
%?
%?%?title（‘Acquisition‘）;
%?xlabel（‘code?phase‘）;
%?ylabel（‘normalized?coherent?integration‘）
%?title（‘new?method‘）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????3010??2020-10-02?17:59??╡┌╚2╒┬\acquire_compare.m
?????文件????????9221??2020-10-02?17:59??╡┌╚2╒┬\bbsignal.m
?????文件????????1429??2020-10-02?17:59??╡┌╚2╒┬\bishe_example.m
?????文件????????1351??2020-10-02?17:59??╡┌╚2╒┬\codegen.asv
?????文件????????1254??2020-10-02?17:59??╡┌╚2╒┬\codegen.m
?????文件?????????357??2020-10-02?17:59??╡┌╚2╒┬\detobi.asv
?????文件?????????357??2020-10-02?17:59??╡┌╚2╒┬\detobi.m
?????文件????????5991??2020-10-02?17:59??╡┌╚2╒┬\fNavdata.m
?????文件????????1404??2020-10-02?17:59??╡┌╚2╒┬\gpstime.m
?????文件????????2006??2020-10-02?17:59??╡┌╚2╒┬\hidynamic_definition.m
?????文件????????1381??2020-10-02?17:59??╡┌╚2╒┬\jammer.m
?????文件????????3520??2020-10-02?17:59??╡┌╚2╒┬\path.m
?????文件????????4956??2020-10-02?17:59??╡┌╚2╒┬\Position.m
?????文件????????4715??2020-10-02?17:59??╡┌╚2╒┬\widenull.m
?????文件?????????215??2020-10-02?17:59??╡┌╚2╒┬\│╠╨≥╦╡├≈.txt