%%%GSC廣義旁瓣相消器，阻塞矩陣B2，衛星信號CA碼和載波，一個點頻干擾%%%
clc;
clear?all;
close?all;
ima=sqrt（-1）;

%%天線參數設定%%
N=16;????????????????????????????????%陣元數N
d_lamda=0.5;?????????????????????????%陣元間距d與入射波波長lamda的比值d/lamda
theta=-90:0.1: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*An*（randn（N-1