clc;?
close?all
clear?all;?
M=20;%%?the?number?of?antennas

L=200;????????????????????????????????????%%?sample?number?
thetas=45;????????????????????????????????%信號(hào)入射角度
thetai=[-30?30];??thetai1=[40?50?];?thetai2=[42?48];????????????????????????%干擾入射角度
n=[0:M-1]‘;????????????????????????????????%n?

vs=exp（-j*pi*n*sin（thetas/180*pi））;?%?vs1=exp（-j*pi*n1*sin（thetas/180*pi））;??vs2=exp（-j*pi*n2*sin（thetas/180*pi））;???????%信號(hào)方向矢量
vi=exp（-j*pi*n*sin（thetai/180*pi））;?vi1=exp（-j*pi*n*sin（thetai1/180*pi））;vi2=exp（-j*pi*n*sin（thetai2/180*pi））;??????%干擾方向矢量
f=16000;??????????????????????????????????%?carrier?frequency
t=[0:1:L-1]/200;?
snr=20;????????????????????????????????%信噪比
inr=15;????????????????????????????????%干噪比

xs=sqrt（10^（snr/10））*vs*exp（j*2*pi*f*t）;?%xs1=sqrt（10^（snr1/10））*vs*exp（j*2*pi*f*t）;?xs2=sqrt（10^（snr2/10））*vs*exp（j*2*pi*f*t）;?%構(gòu)造有用信號(hào)
xi=sqrt（10^（inr/10）/2）*vi*[randn（length（thetai）L）+j*randn（length（thetai）L）];
xi1=sqrt（10^（inr/10）/2）*vi*[randn（length（thetai1）L）+j*randn（length（thetai1）L