clear?all
close?all
clc
j=sqrt（-1）;
M=16;
m=0:M-1;
theta=[0-40];
theta?=?theta*pi/180;
us=0;
A0=exp（j*pi*m‘*sin（theta（1）））;%信號源的方向矢量
Ai=exp（j*pi*m‘*sin（theta（2）））;%干擾的方向矢量
SNR=20;??%?信噪比（可變）
INR=30;??%干噪比
P0=10.^（SNR/10）;?
Pi=10.^（INR/10）;
o=pi/180;
xxl=-90:1:90;
xx=xxl*o;
a=exp（j*pi*m‘*sin（xx））;??????????????????%?掃描矢量;
%%%%%%%%%%%最優波束形成%%%%%%%%%%%%%%%%
L=5000;
Xs=A0*sqrt（P0*2）*wgn（1L1‘complex‘）;%期望
Xi=Ai*sqrt（Pi*2）*wgn（1L1‘complex‘）;%干擾
Xn=sqrt（2）*wgn（ML1‘complex‘）;%噪聲
Xx=Xs+Xi+Xn;
%%%Xx=Xi+Xw;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure
Rs=P0*A0*A0‘;??????????????????????????????????????%期望信號協方差矩陣
Ri=Pi*Ai*Ai‘;??????????????????????????????????????%干