%數(shù)據(jù)中含有或不含有期望信號(hào)對(duì)方向圖有影響，當(dāng)MWF也采用理想數(shù)據(jù)時(shí)，性能相同且采樣數(shù)可以大大減小;
%當(dāng)變換B_0、h0為B、Wq，則為LC-MWF
clc;
clear?all;
close?all;
N=8;%輔助天線的數(shù)目
d=0.5;
theta=-90:0.5:90;
theta_signal=0;????????????%來波方向
theta_jamming1=-45;????????%干擾方向
theta_jamming2=25;?????????%干擾方向
theta_jamming3=55;?????????%干擾方向
L=128;%采樣單元數(shù)
steer_signal=exp（1i*2*pi*d*sin（theta_signal*pi/180）*[0:N-1]‘）;?%期望信號(hào)導(dǎo)向矢量
%干擾信號(hào)導(dǎo)向矢量
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]‘）;
%訓(xùn)練數(shù)據(jù)中不含期望信號(hào)時(shí)和含有期望信號(hào)時(shí)
for?i=1:L;
????amp_signal=sqrt（2*10）*randn（1）;%信號(hào)的幅度隨機(jī)產(chǎn)生，保證信號(hào)之間是不相關(guān)的
????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））;%接收信號(hào)（信號(hào)+干擾+噪聲）
end
Rs=10*steer_signal*steer_signal‘;%期望信號(hào)自相關(guān)矩陣，相當(dāng)于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;%接收信號(hào)協(xié)方差矩陣？
%R=s*s‘/L;
%Sx=1/L*s*s‘;?%接收信號(hào)協(xié)方差矩陣
Sx=1/L*jam*jam‘;
us=d*sin（theta_signal*pi/180）;?%期望信號(hào)方向空間頻率;
d0