广义旁瓣相消器

%%%GSC广义旁瓣相消器，阻塞矩阵B1，卫星信号CA码和载波，一个点频干扰%%%
clc;
clear all;
close all;
ima=sqrt（-1）;

%%天线参数设定%%
N=16;                                %阵元数N
d_lamda=0.5;                         %阵元间距d与入射波波长lamda的比值d/lamda
theta=-90:0.5:90;                    %方向图自变量theta
theta0=30;                           %期望（卫星）信号的角度30
theta1=65;                           %点频干扰方向65
L=1023;                              %快拍数
f1=1575.42*10^6;                     %GPS卫星载波信号L1的频率

%%期望信号、干扰、噪声的建模%%
s1=generateCAcode（3）;
dt=1e-6;                                   %CA码码元宽度为0.977微秒，相当于1e-6秒。
T=1.023*1e-3;                              %
t=dt:dt:T;
Ac=sqrt（2）*10^（-8）;                        %已知信号功率-130dBm，折算到载波幅值，信号功率10^（-16）W，即-160dB。
C1=Ac*sin（2*pi*f1*t）;                      %载波
s=s1.*C1;                                  %调制后的卫星信号

Ajam=sqrt（2）*10^（-5.75）;                   %干扰的幅值 已知信干比SIR=-45dB，信号功率-160dB，折算得到，干扰功率-115dB即10^（-11.5）W。     
jam=Ajam*cos（2*pi*f1*t）;                   %与期望卫星信号载波同频率的点频干扰

An=2*10^（-20.4）;                           %高斯白噪声强度-174dBm/Hz，带宽按2M考虑
n=0.5*An*（randn（NL）+ima*randn（NL））;      %高斯白噪声注意复随机信号randn+j*randn产生的是强度为2的高斯白噪声
n0=0.5*A