%?陣列方向圖（使用MVDR方法計算加權向量）
%?陣元數為16，間距為半波長，信號頻率為1.3G，指向角為0度，干擾三個在－40，20，45
%?考慮波束指向方向無信號情況
%?2007.5.22無信號

clear?all
clc
f0?=?1.3*10^9;????????????????????%?信號頻率
fs?=?3*10^9;???????????????????%?采樣頻率
snr?=?[40;40];??????????????????%?干擾的信噪比
w0?=?20/180*pi;??????????????????%?目標指向角
M?=?9;?????????????????????????%?陣元數為M
K?=?2;??????????????????????????%?干擾信號數目
seta?=[20/180*pi45/180*pi];???%?干擾方向??????????
N?=?100;????????????????????????%?采樣點數
c?=?3*10^8;?????????????????????%?光速
d?=?0.5*c/f0;???????????????????%?陣元間距
x?=?zeros（M1）;?????????????????%?數據矢量
R?=?zeros（MM）;?????????????????%?接收數據協方差矩陣
tic
%?信號源
nn?=?1:N;%下面設計的是干擾信號nn為列N的行向量
for?k=1:K
????for?n=1:N
????????fai（1n）?=?rand;
????end
????s（k:）?=?exp（i*2*pi*（f0*nn/fs+fai*k））;%構造信號
end

%?構造方向矢量
thet