%%%GSC廣義旁瓣相消器，阻塞矩陣B1，衛(wèi)星信號(hào)CA碼和載波，一個(gè)點(diǎn)頻干擾%%%
clc;
clear?all;
close?all;
ima=sqrt（-1）;

%%天線參數(shù)設(shè)定%%
N=16;????????????????????????????????%陣元數(shù)N
d_lamda=0.5;?????????????????????????%陣元間距d與入射波波長(zhǎng)lamda的比值d/lamda
theta=-90:0.5:90;????????????????????%方向圖自變量theta
theta0=30;???????????????????????????%期望（衛(wèi)星）信號(hào)的角度30
theta1=65;???????????????????????????%點(diǎn)頻干擾方向65
L=1023;??????????????????????????????%快拍數(shù)
f1=1575.42*10^6;?????????????????????%GPS衛(wèi)星載波信號(hào)L1的頻率

%%期望信號(hào)、干擾、噪聲的建模%%
s1=generateCAcode（3）;
dt=1e-6;???????????????????????????????????%CA碼碼元寬度為0.977微秒，相當(dāng)于1e-6秒。
T=1.023*1e-3;??????????????????????????????%
t=dt:dt:T;
Ac=sqrt（2）*10^（-8）;????????????????????????%已知信號(hào)功率-130dBm，折算到載波幅值，信號(hào)功率10^（-16）W，即-160dB。
C1=Ac*sin（2*pi*f1*t）;??????????????????????%載波
s=s1.*C1;??????????????????????????????????%調(diào)制后的衛(wèi)星信號(hào)

Ajam=sqrt（2）*10^（-5.75）;???????????????????%干擾的幅值?已知信干比SIR=-45dB，信號(hào)功率-160dB，折算得到，干擾功率-115dB即10^（-11.5）W。?????
jam=Ajam*cos（2*pi*f1*t）;???????????????????%與期望衛(wèi)星信號(hào)載波同頻率的點(diǎn)頻干擾

An=2*10^（-20.4）;???????????????????????????%高斯白噪聲強(qiáng)度-174dBm/Hz，帶寬按2M考慮
n=0.5*An*（randn（NL）+ima*randn（NL））;??????%高斯白噪聲注意復(fù)隨機(jī)信號(hào)randn+j*randn產(chǎn)生的是強(qiáng)度為2的高斯白噪聲
n0=0.5*A