高动态 GPS 自适应抗干扰算法研究

高动态 GPS 自适应抗干扰算法研究

% 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  ╡┌╚²╒┬\acquire_compare.m
     文件        9221  2020-10-02 17:59  ╡┌╚²╒┬\bbsignal.m
     文件        1429  2020-10-02 17:59  ╡┌╚²╒┬\bishe_example.m
     文件        1351  2020-10-02 17:59  ╡┌╚²╒┬\codegen.asv
     文件        1254  2020-10-02 17:59  ╡┌╚²╒┬\codegen.m
     文件         357  2020-10-02 17:59  ╡┌╚²╒┬\detobi.asv
     文件         357  2020-10-02 17:59  ╡┌╚²╒┬\detobi.m
     文件        5991  2020-10-02 17:59  ╡┌╚²╒┬\fNavdata.m
     文件        1404  2020-10-02 17:59  ╡┌╚²╒┬\gpstime.m
     文件        2006  2020-10-02 17:59  ╡┌╚²╒┬\hidynamic_definition.m
     文件        1381  2020-10-02 17:59  ╡┌╚²╒┬\jammer.m
     文件        3520  2020-10-02 17:59  ╡┌╚²╒┬\path.m
     文件        4956  2020-10-02 17:59  ╡┌╚²╒┬\Position.m
     文件        4715  2020-10-02 17:59  ╡┌╚²╒┬\widenull.m
     文件         215  2020-10-02 17:59  ╡┌╚²╒┬\│╠╨≥╦╡├≈.txt