%數據中含有或不含有期望信號對方向圖有影響，當MWF也采用理想數據時，性能相同且采樣數可以大大減小;
%當變換B_0、h0為B、Wq，則為LC-MWF
clc;
clear?all;
close?all;
N=8;%輔助天線的數目
d=0.5;
theta=-90:0.5:90;
theta_signal=0;????????????%來波方向
theta_jamming1=-45;????????%干擾方向
theta_jamming2=25;?????????%干擾方向
theta_jamming3=55;?????????%干擾方向
L=128;%采樣單元數
steer_signal=exp（1i*2*pi*d*sin（theta_signal*pi/180）*[0:N-1]‘）;?%期望信號導向矢量
%干擾信號導向矢量
steer_jamming1=exp（1i*2*pi*d*sin（theta_jamming1*pi/180）*[0:N-1]‘）;?
steer_jamming2=exp（1i*2*pi*d*sin（theta_jamming2*pi/180）*[0:N-1]‘）;
steer_jamming3=exp（1i*2*pi*d*sin（theta_jamming3*pi/180）*[0:N-1]‘）;
%訓練數據中不含期望信號時和含有期望信號時
for?i=1:L;
????amp_signal=sqrt（2*10）*randn（1）;%信號的幅度隨機產生，保證信號之間是不相關的
????amp_jamming1=sqrt（2*10^4）*randn（1）;?%4
????amp_jamming2=sqrt（2*10^3.5）*randn（1）;?%3.5
????amp_jamming3=sqrt（2*10^5）*randn（1）;?%5
????amp_noise=sqrt（2）;
????jam（:i）=amp_jamming1*exp（1i*2*pi*d*sin（theta_jamming1*pi/180）*[0:N-1]‘）...
?????????????+amp_jamming2*exp（1i*2*pi*d*sin（theta_jamming2*pi/180）*[0:N-1]‘）...
?????????????+amp_jamming3*exp（1i*2*pi*d*sin（theta_jamming3*pi/180）*[0:N-1]‘）...
?????????????+amp_noise*（randn（N1）+1i*randn（N1））;%干擾+噪聲
????s（:i）=amp_signal*exp（1i*2*pi*d*sin（theta_signal*pi/180）*[0:N-1]‘）...
???????????+amp_jamming1*exp（1i*2*pi*d*sin（theta_jamming1*pi/180）*[0:N-1]‘）...
???????????+amp_jamming2*exp（1i*2*pi*d*sin（theta_jamming2*pi/180）*[0:N-1]‘）...
???????????+amp_jamming3*exp（1i*2*pi*d*sin（theta_jamming3*pi/180）*[0:N-1]‘）...
???????????+amp_noise*（randn（N1）+1i*randn（N1））;%接收信號（信號+干擾+噪聲）
end
Rs=10*steer_signal*steer_signal‘;%期望信號自相關矩陣，相當于X（t）
Rj=10^4*steer_jamming1*steer_jamming1‘...
???+10^3.5*steer_jamming2*steer_jamming2‘...
???+10^5*steer_jamming3*steer_jamming3‘;
Rn=eye（N）;
Rx=Rs+Rj+Rn;%接收信號協方差矩陣？
%R=s*s‘/L;
%Sx=1/L*s*s‘;?%接收信號協方差矩陣
Sx=1/L*jam*jam‘;
us=d*sin（theta_signal*pi/180）;?%期望信號方向空間頻率;
d0